SEARCH · Engineering Papers

Results for “RANS”

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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Numerical challenges for energy conservation in N -body simulations of collapsing self-interacting dark matter halos

Dark matter (DM) halos can be subject to gravothermal collapse if the DM is not collisionless, but engaged in strong self-interactions instead. When the scattering is able to efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. However, a comparison with observations requires accurate predictions. One method to obtain such predictions is to use N-body simulations. Collapsed halos are extreme systems that pose severe challenges when applying state-of-the-art codes to model self-interacting dark matter (SIDM). In this work, we investigate the root of such problems, with a focus on energy non-conservation. Moreover, we discuss possible strategies to avoid them. We ran N-body simulations, both with and without SIDM, of an isolated DM-only halo and we adjusted the numerical parameters to check the accuracy of the simulation. We find that not only the numerical scheme for SIDM can lead to energy non-conservation, but also the modelling of gravitational interaction and the time integration are problematic. The main issues we find are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; and (c) SIDM velocity kicks breaking the time symmetry. Tuning the parameters of the simulation to achieve a high level of accuracy allows us to conserve energy not only at early stages of the evolution, but also later on. However, the cost of the simulations becomes prohibitively large as a result. Some of the problems that make the simulations of the gravothermal collapse phase inaccurate can be overcome by choosing appropriate numerical schemes. However, other issues still pose a challenge. Our findings motivate further works on addressing the challenges in simulating strong DM self-interactions.

dark matter↗

DOE OSTI2025

The aluminium-26 distribution in a cosmological simulation of a Milky Way-type Galaxy

Context . The 1.8 MeV γ -rays corresponding to the decay of the radioactive isotope 26 Al (with a half-life of 0.72 Myr ) have been observed by the SPI detector on the INTEGRAL spacecraft and extensively used as a tracer of star formation and current nucleosynthetic activity in the Milky Way Galaxy. Further information is encoded in the observation related to the higher 26 Al content found in regions of the Galaxy with the highest line-of-sight (LoS) velocity relative to an observer located in the Solar System. However, this feature remains unexplained. Aims . We ran a cosmological “zoom-in” chemodynamical simulation of a Milky Way-type galaxy, including the production and decays of radioactive nuclei in a fully self-consistent way. We then analyzed the results to follow the evolution of 26 Al throughout the lifetime of the simulated galaxy to provide a new method for interpreting the 26 Al observations. Methods . We included the massive star sources of 26 Al in the Galaxy and its radioactive decay into a state-of-the-art galactic chemical evolution model, coupled with cosmological growth and hydrodynamics. This approach allowed us to follow the spatial and temporal evolution of the 26 Al content in the simulated galaxy. Results . Our results are in agreement with the observations with respect to the fact that gas particles in the simulation with relatively higher 26 Al content also have the highest LoS velocities. On the other hand, gas particles with relatively lower 26 Al content (i.e., not bright enough to be observed) generally display the lowest LoS velocities. However, this result is not conclusive because the overall rotational velocity of our simulated galaxy is higher than that observed for cold CO gas in the Milky Way Galaxy. Furthermore, we found no significant correlation between gas temperature, rotational velocity, and 26 Al content at any given radius. We also found the presence of transient 26 Al-rich spots at low LoS velocities and we show that one such spot had been captured by the INTEGRAL/SPI data. Based on our model, we present a prediction for the detection of 1.8 MeV γ -rays by the future COSI mission. We find that according to our model, the new instrument will be able to observe similar 26 Al-emission patterns to those seen by INTEGRAL/SPI.

Wehmeyer, B. (ORCID:0009000254149292)↗

DOE OSTI2025

Impact of T - and ρ -dependent decay rates and new (n, γ ) cross-sections on the s process in low-mass asymptotic giant branch stars

Aims. We study the impact of nuclear input related to weak-decay rates and neutron-capture reactions on predictions for the slow neutron-capture process (s process) in asymptotic giant branch (AGB) stars. We provide the first database of surface abundances and stellar yields of the isotopes heavier than iron from the Monash models. Methods. We ran nucleosynthesis calculations with the Monash post-processing code for seven stellar structure evolution models of low-mass AGB stars with three different sets of nuclear inputs. The reference set has constant decay rates and represents the set used in the previous Monash publications. The second set contains the temperature and density dependence of β decays and electron captures based on the default rates of nuclear NETwork GENerator (NETGEN). In the third set, we further update 92 neutron-capture rates based on re-evaluated experimental cross sections from the ASTrophysical Rate and rAw data Library. We compare and discuss the predictions of the sets relative to each other in terms of isotopic surface abundances and total stellar yields. We also compare the results to isotopic ratios measured in presolar stardust silicon carbide (SiC) grains from AGB stars. Results. The new sets of models result in a ∼66% solar s-process contribution to the p-nucleus 152 Gd, confirming that this isotope is predominantly made by the s process. The nuclear input updates result in predictions for the 80 Kr/ 82 Kr ratio in the He intershell and surface 64 Ni/ 58 Ni, 94 Mo/ 96 Mo, and 137 Ba/ 136 Ba ratios that are more consistent with the corresponding ratios measured in stardust; however, the new predicted 138 Ba/ 136 Ba ratios are higher than the typical values of the SiC grains. The W isotopic anomalies are in agreement with data from the analyses of other meteoritic inclusions. We confirm that the production of 176 Lu and 205 Pb is affected by too large uncertainties in their decay rates from NETGEN.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Chemical templates of the Central Molecular Zone

Context . The Central Molecular Zone (CMZ) of the Milky Way exhibits extreme conditions, including high gas densities, elevated temperatures, enhanced cosmic-ray ionization rates, and large-scale dynamics. This makes it a perfect laboratory for astrochemical studies. With large-scale molecular surveys revealing increasing chemical and physical complexity in the CMZ, it is essential to develop robust methods to decode the chemical information embedded in this extreme region. Aims . A key step to interpreting the molecular richness found in the CMZ is building chemical templates tailored to its diverse conditions. In particular, understanding how CMZ environments affect shock and protostellar chemistry is crucial. The combined impact of high ionization, elevated temperatures, and dense gas remains insufficiently explored for observable tracers. Methods . For this study, we utilized UCLCHEM , a gas-grain time-dependent chemical model, to link physical conditions with their corresponding molecular signatures and identify key tracers of temperature, density, ionization, and shock activity. To achieve this, we ran a grid of models of shocks and protostellar objects representative of typical CMZ conditions, focusing on 24 species, including complex organic molecules. Results . Shocked and protostellar environments show distinct evolutionary timescales (≲10 4 vs. ≳10 4 years); 300 K emerges as a key temperature threshold for chemical differentiation. We find that cosmic-ray ionization and temperature are the main drivers of chemical trends. HCO + , H 2 CO, and CH 3 SH trace ionization, while HCO, HCO + , CH 3 SH, CH 3 NCO, and HCOOCH 3 show consistent abundance contrasts between shocks and protostellar regions over similar temperature ranges. Conclusions . We characterized the behavior of 24 species in protostellar and shock-related environments. While our models underpredict some complex organics in shocks, they reproduce observed trends for most species, supporting scenarios involving a need for recurring shocks in Galactic Center clouds and enhanced ionization toward Sgr B2(N2). Future work should assess the role of shock recurrence and metallicity in shaping chemistry.

Galaxy: center↗

DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗

DOE OSTI2021

Ultrasonic oscillatory two-phase flow in microchannels

Experimental and numerical investigations are performed to provide an assessment of the transport behavior of an ultrasonic oscillatory two-phase flow in a microchannel. The work is inspired by the flow observed in an innovative ultrasonic fabric drying device using a piezoelectric bimorph transducer with microchannels, where a water-air two-phase flow is transported by harmonically oscillating microchannels. The flow exhibits highly unsteady behavior as the water and air interact with each other during the vibration cycles, making it significantly different from the well-studied steady flow in microchannels. Here, the computational fluid dynamics (CFD) modeling is realized by combing the turbulence Reynolds-averaged Navier-Stokes (RANS) k – ω model with the phase-field method to resolve the dynamics of the two-phase flow. The numerical results are qualitatively validated by the experiment. Through parametric studies, we specifically examined the effects of vibration conditions (i.e., frequency and amplitude), microchannel taper angle, and wall surface contact angle (i.e., wettability) on the flow rate through the microchannel. The results will advance the potential applications where oscillatory or general unsteady microchannel two-phase flows may be present.

42 ENGINEERING↗

DOE OSTI2021

Computing the Paschen curve for argon with speed-limited particle-in-cell simulation

Upon inclusion of collisions, the speed-limited particle-in-cell (SLPIC) simulation method successfully computed the Paschen curve for argon. The 1D3V simulations modeled an electron cascade across an argon-filled capacitor, including electron-neutral ionization, electron-neutral elastic collisions, electron-neutral excitation, and ion-induced secondary electron emission. In electrical breakdown, the timescale difference between ion and electron motion makes traditional PIC methods computationally slow. To decrease this timescale difference and speed up computation, here we used SLPIC, a time-domain algorithm that limits the speed of the fastest electrons in the simulation. The SLPIC algorithm facilitates a straightforward, fully kinetic treatment of dynamics and collisions. SLPIC was as accurate as PIC, but ran up to 200 times faster. SLPIC accurately computed the Paschen curve for argon over three orders of magnitude in pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Multiscale porous elastomer substrates for multifunctional on-skin electronics with passive-cooling capabilities

In addition to mechanical compliance, achieving the full potential of on-skin electronics needs the introduction of other features. For example, substantial progress has been achieved in creating biodegradable, self-healing, or breathable, on-skin electronics. However, the research of making on-skin electronics with passive-cooling capabilities, which can reduce energy consumption and improve user comfort, is still rare. Herein, we report the development of multifunctional on-skin electronics, which can passively cool human bodies without needing any energy consumption. This property is inherited from multiscale porous polystyrene-block-poly(ethylene-ran-butylene)-block-polystyrene (SEBS) supporting substrates. The multiscale pores of SEBS substrates, with characteristic sizes ranging from around 0.2 to 7 µm, can effectively backscatter sunlight to minimize heat absorption but are too small to reflect human-body midinfrared radiation to retain heat dissipation, thereby delivering around 6 °C cooling effects under a solar intensity of 840 W∙m -2 . Other desired properties, rooted in multiscale porous SEBS substrates, include high breathability and outstanding waterproofing. The proof-of-concept bioelectronic devices include electrophysiological sensors, temperature sensors, hydration sensors, pressure sensors, and electrical stimulators, which are made via spray printing of silver nanowires on multiscale porous SEBS substrates. The devices show comparable electrical performances with conventional, rigid, nonporous ones. Also, their applications in cuffless blood pressure measurement, interactive virtual reality, and human–machine interface are demonstrated. Notably, the enabled on-skin devices are dissolvable in several organic solvents and can be recycled to reduce electronic waste and manufacturing cost. Such on-skin electronics can serve as the basis for future multifunctional smart textiles with passive-cooling functionalities.

42 ENGINEERING↗

DOE OSTI2023

Mechanism of RanGTP priming H2A-H2B release from Kap114 in an atypical RanGTP•Kap114•H2A-H2B complex

Previously, we showed that the nuclear import receptor Importin-9 wraps around the H2A-H2B core to chaperone and transport it from the cytoplasm to the nucleus. However, unlike most nuclear import systems where RanGTP dissociates cargoes from their importins, RanGTP binds stably to the Importin-9•H2A-H2B complex, and formation of the ternary RanGTP•Importin-9•H2A-H2B complex facilitates H2A-H2B release to the assembling nucleosome. It was unclear how RanGTP and the cargo H2A-H2B can bind simultaneously to an importin, and how interactions of the three components position H2A-H2B for release. Here, we show cryo-EM structures of Importin-9•RanGTP and of its yeast homolog Kap114, including Kap114•RanGTP, Kap114•H2A-H2B, and RanGTP•Kap114•H2A-H2B, to explain how the conserved Kap114 binds H2A-H2B and RanGTP simultaneously and how the GTPase primes histone transfer to the nucleosome. In the ternary complex, RanGTP binds to the N-terminal repeats of Kap114 in the same manner as in the Kap114/Importin-9•RanGTP complex, and H2A-H2B binds via its acidic patch to the Kap114 C-terminal repeats much like in the Kap114/Importin-9•H2A-H2B complex. Ran binds to a different conformation of Kap114 in the ternary RanGTP•Kap114•H2A-H2B complex. Here, Kap114 no longer contacts the H2A-H2B surface proximal to the H2A docking domain that drives nucleosome assembly, positioning it for transfer to the assembling nucleosome or to dedicated H2A-H2B chaperones in the nucleus.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Droplet and Particle Methods to Investigate Turbulent Particle Laden Jets

The SARS-CoV-2 pandemic has heightened the interest in particle-laden turbulent jets generated by breathing, talking, coughing and sneezing, and how these can contribute to disease transmission. We present quantitative measurement methods for such flows, while exploring and offering improvements for common shortcomings. We generate jets consisting of either liquid droplets or solid particles in an isothermal, quiescent and electrically isopotential experimental chamber that was constructed to control the effects of ambient forcing on jet behavior. For liquid droplets, we find promise in surface deposition analysis based on fluorescent tracer use. For particles, we explore the performance of commercially available adhesive sampling strips and develop conductive grounded carbon tape based sampling strips. We explore ways in which the smallest of thermal gradients or electrostatic charge issues can affect particle dispersion, and suggest practical methods to address these issues. The developed methods are applied to study the simultaneous deposition of (Formula presented.) 25, 50 and 200 μm solid particles from a particle laden turbulent jet with a mean velocity of 33.2 m/s. The deposition location as a function of particle size was compared to results from a simple numerical RANS model, and illustrates ways in which imprecise initial or boundary conditions can lead to a notable deviation from experimental results. The differences in deposition pattern seen in experimental and numerical results despite a carefully controlled environment and characterized particle ejection indicate the need for a more stringent numerical model validation, especially when studying fate and transport of mid-range (neither purely aerosol or ballistic) sized particles.

42 ENGINEERING↗

DOE OSTI2022

Mechanical loading prediction through accelerometry data during walking and running

ABSTRACT Currently, there is no way to assess mechanical loading variables such as peak ground reaction forces (pGRF) and peak loading rate (pLR) in clinical settings. The purpose of this study was to develop accelerometry‐based equations to predict both pGRF and pLR during walking and running. One hundred and thirty one subjects (79 females; 76.9 ± 19.6 kg) walked and ran at different speeds (2–14 km·h −1 ) on a force plate–instrumented treadmill while wearing accelerometers at their ankle, lower back and hip. Regression equations were developed to predict pGRF and pLR from accelerometry data. Leave‐one‐out cross‐validation was used to calculate prediction accuracy and Bland–Altman plots. Our pGRF prediction equation was compared with a reference equation previously published. Body mass and peak acceleration were included for pGRF prediction and body mass and peak acceleration rate for pLR prediction. All pGRF equation coefficients of determination were above 0.96, and a good agreement between actual and predicted pGRF was observed, with a mean absolute percent error (MAPE) below 7.3%. Accuracy indices from our equations were better than previously developed equations. All pLR prediction equations presented a lower accuracy compared to those developed to predict pGRF. Walking and running pGRF can be predicted with high accuracy by accelerometry‐based equations, representing an easy way to determine mechanical loading in free‐living conditions. The pLR prediction equations yielded a somewhat lower prediction accuracy compared with the pGRF equations.

Veras, Lucas↗

DOE OSTI2021

Verification of induction zone models for wind farm annual energy production estimation

Fast engineering wake models are the backbone of wind farm annual energy production (AEP) estimators, whereas the addition of induction zone models in existing tools is a more recent response to rising concerns over wind farm blockage associated losses. Here, "blockage" describes the combined induction fields of all wind turbines inside a farm. Unlike the term might suggest, blockage not only reduces flow speeds, but also increases them; for instance along the outer edges of wind farms. Evaluating the overall impact on AEP of these gains and losses necessitates accurate wind farm induction models. Whilst engineering wake models are tuned to predict wind farm performance, existing induction zone models are all derived for accurately predicting the near field induction not the far field value—important for wind farm simulations. This paper presents the induction models implemented in the wind farm simulation tool PyWake, as well as results from novel analytical models and compares their far field predictions to RANS-AD simulations of different turbines. We demonstrate that when including blockage in AEP simulations, the downstream speed-ups need to be included to avoid an unrealistic bias toward AEP loss and that wake expansion significantly impacts induction at rated thrust levels.

17 WIND ENERGY↗

DOE OSTI2026

EPICS for small-scale laboratories with Python soft IOCs

While the Experimental Physics and Industrial Control System (EPICS) is widely used at large laboratories for slow controls and instrumentation, the deployment of a full EPICS installation can be difficult, with a steep learning curve to new users. Taking advantage of the pythonSoftIOC module, we developed an EPICS slow controls implementation for Jefferson Lab's Hall B cryotarget written entirely in Python and based on software IOCs that communicate with instruments over Ethernet. Here, this system ran successfully, interfacing with Jefferson Lab's full EPICS network, and we offer it as an example of the capabilities of pythonSoftIOC to build lightweight, yet robust and flexible instrumentation platforms that would be easily adapted for use at a small-scale laboratory. University groups can use these examples to build complete slow controls systems, from device communication to data archiving and display, using open-source, mature EPICS tools and student-friendly Python as an alternative to expensive and proprietary systems such as LabVIEW.

Computing↗

DOE OSTI2025

Plasma phosphorylated tau217 strongly associates with memory deficits in the Alzheimer’s disease spectrum

Abstract Plasma phosphorylated tau (p-tau) biomarkers open unprecedented opportunities for identifying carriers of Alzheimer’s disease pathophysiology in early disease stages using minimally invasive techniques. Plasma p-tau biomarkers are believed to reflect tau phosphorylation and secretion. However, it remains unclear to what extent the magnitude of plasma p-tau abnormalities reflects neuronal network disturbance in the form of cognitive impairment. To address this question, we included 103 cognitively unimpaired elderly and 40 cognitively impaired, amyloid-β-positive individuals from the TRIAD cohort, in addition to 336 cognitively unimpaired and 216 cognitively impaired, amyloid-β-positive older adults from the BioFINDER-2 cohort. Participants had tau PET scans, amyloid PET scans or amyloid CSF, p-tau217, p-tau181 and p-tau231 blood measures, structural T1-MRI and cognitive assessments. In this cross-sectional study, we used regression models and correlation analyses to assess the relationship between plasma biomarkers and cognitive scores. Furthermore, we applied receiver operating characteristic curves to assess cognitive impairment across plasma biomarkers. Finally, we categorized participants into amyloid (A), p-tau (T1) and tau PET (T2) positive (+) or negative (−) profiles and ran non-parametric comparisons to assess differences across cognitive domains. We found that plasma p-tau217 was more associated with cognitive performance than p-tau181 and p-tau231 and that this relationship was particularly strong for memory scores (TRIAD: βp-tau217 = −0.53, βp-tau181 = −0.35 and βp-tau231 = −0.24; BioFINDER-2: βp-tau217 = −0.52, βp-tau181 = −0.24 and βp-tau231 = −0.29). Associations in amyloid-β-positive participants resembled these results, but other cognitive scores also showed strong associations in cognitively impaired individuals. Moreover, plasma p-tau217 outperformed plasma p-tau181 and plasma p-tau231 in identifying memory impairment (area under the curve values for TRIAD: p-tau217 = 0.86, p-tau181 = 0.77 and p-tau231 = 0.75; and for BioFINDER-2: p-tau217 = 0.86, p-tau181 = 0.76 and p-tau231 = 0.81) and in identifying executive function impairment only in the BioFINDER-2 cohort (p-tau217 = 0.82, p-tau181 = 0.76 and p-tau231 = 0.76). Lastly, we showed that subtle memory deficits were present in A+T1+T2− participants for plasma p-tau217 (P = 0.007) and plasma p-tau181 (P = 0.01) in the TRIAD cohort and for all biomarkers across cognitive domains in A+T1+T2− and A+T1+T2− individuals (P < 0.001 in all) in the BioFINDER-2 cohort. The A+T1+T2− individuals showed cognitive deficits in both cohorts (P < 0.001 in all). Together, our results suggest that plasma p-tau217 stands out as a biomarker capable of identifying memory deficits attributable to Alzheimer’s disease and that memory impairment certainly occurs in amyloid-β- and plasma p-tau-positive individuals who have no significant amounts of tau in the neocortex.

Neurosciences & Neurology↗

DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗

DOE OSTI2023

Measurement of the branching fraction and search for CP violation in D 0 → K $^{0}_{S}$ K $^{0}_{S}$ ⁢π + ⁢π - decays at Belle

We measure the branching fraction for the Cabibbo-suppressed decay D 0 → K $^{0}_{S}$ K $^{0}_{S}$ ⁢π + ⁢π - and search for CP violation via a measurement of the CP asymmetry A CP as well as the T-odd triple-product asymmetry a$^{T}_{CP}$. We use 922 fb -1 of data recorded by the Belle experiment, which ran at the KEKB asymmetric-energy e + ⁢e - collider. The branching fraction is measured relative to the Cabibbo-favored normalization channel D 0 → K $^{0}_{S}$ ⁢π + ⁢π - ; the result is $\mathscr{B}$⁡( D 0 → K $^{0}_{S}$ K $^{0}_{S}$ ⁢π + ⁢π - )=[4.79±0.08⁢(stat)±0.10⁢(syst)±0.31⁢(norm)]×10 -4 , where the first uncertainty is statistical, the second is systematic, and the third is from uncertainty in the normalization channel. We also measure A CP =[-2.51±1.44⁢(stat)$^{+0.11}_{-0.10}$⁢(syst)]%, and a$^{T}_{CP}$=[-1.95±1.42⁢(stat)$^{+0.14}_{-0.12}⁢$(syst)]%. These results show no evidence of CP violation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of three Rayleigh-Taylor mixing configurations with gravity reversal

In this work, high-fidelity large-eddy simulation (LES) is performed of Rayleigh-Taylor (RT) mixing in three different configurations involving gravity reversal. In each configuration, LES results are compared with one-dimensional Reynolds-averaged Navier-Stokes (RANS) results, and a deficiency in a commonly used transport equation for the mass-flux velocity, $a_j$, is identified. In the first configuration, a classical two-component RT mixing layer is allowed to develop before it is subjected to rapid acceleration reversal. In the second configuration, a three-component RT mixing layer with an intermediate density layer is allowed to develop before being subjected to rapid acceleration reversal. Finally, in the third configuration, a light layer is interposed between two heavy layers; in this configuration, only one interface is RT-unstable at a time as it undergoes rapid acceleration reversal. In all cases, a commonly used buoyancy production closure in the $a_j$ transport equation is shown to lead to significant over-prediction of mixing layer growth after gravity reversal. An alternative formulation for this closure is then presented which is shown to more accurately capture the stabilization effect of gravity reversal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Analysis of real-fluid thermodynamic effects on turbulent statistics in transcritical channel flows

Wall-bounded turbulence at high-pressure transcritical conditions with intense density fluctuations are encountered in many technical applications. In this study, we analyze the turbulent energy transport in transcritical channel flows specifically focusing on dissipation rate, turbulent kinetic energy budgets, heat fluxes, and momentum-fluctuation statistics; results from this analysis are used to guide the development of turbulent scaling laws. Furthermore, we find that the dissipation rate of turbulent kinetic energy is dominated by the enstrophy in the logarithmic layer, and the fluctuating viscosity results in the reduced tilting of the vortical structures and the attenuation of streamwise vorticity in the near-wall layer; the fluctuating viscosity attenuates the dissipation rate by reducing the shear strain and the enstrophy production. Local equilibrium of the turbulent kinetic energy exists in the logarithmic layer. We show that the real-fluid thermodynamic effects significantly change the turbulent heat flux correlated with the sweep and the ejection events; the density changes alter the turbulent transport and result in noticeable magnitudes of density-fluctuation-related momentum-fluctuation statistics. From these results, scaling laws for the turbulent length scales and turbulent kinetic energy budgets are proposed, thereby contributing to improvement of the wall models in large-eddy simulations and Reynolds-averaged Navier-Stokes (RANS) simulations.

42 ENGINEERING↗