Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PROBABILITY DENSITY FUNCTIONS”

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 271 records · Page 15

Ground Heat Flux Reconstruction Using Bayesian Uncertainty Quantification Machinery and Surrogate Modeling

Ground heat flux (G 0 ) is a key component of the land-surface energy balance of high-latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming, G 0 is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstruct G 0 across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using available G 0 data (measured or modeled) for snow-free period as a reference. When observed G 0 is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state-of-the-art uncertainty quantification methods, the developed G 0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.

54 ENVIRONMENTAL SCIENCES↗

Gas Transfer Across Air‐Water Interfaces in Inland Waters: From Micro‐Eddies to Super‐Statistics

In inland water covering lakes, reservoirs, and ponds, the gas exchange of slightly soluble gases such as carbon dioxide, dimethyl sulfide, methane, or oxygen across a clean and nearly flat air‐water interface is routinely described using a water‐side mean gas transfer velocity $\overline{k_{L}}$, where overline indicates time or ensemble averaging. The micro‐eddy surface renewal model predicts $\overline{k_{L}}$ = α o Sc -1/2 ($v\bar{ϵ}$) 1/4 , where Sc is the molecular Schmidt number, $v$ is the water kinematic viscosity, and $\bar{ϵ}$ is the waterside mean turbulent kinetic energy dissipation rate at or near the interface. While α o = 0.39 - 0.46 has been reported across a number of data sets, others report large scatter or variability around this value range. It is shown here that this scatter can be partly explained by high temporal variability in instantaneous ϵ around $\bar{ϵ}$, a mechanism that was not previously considered. As the coefficient of variation (CV e ) in ϵ increases, α o must be adjusted by a multiplier (1 = CV e 2 ) -3/32 that was derived from a log‐normal model for the probability density function of ϵ. Reported variations in α o with a macro‐scale Reynolds number can also be partly attributed to intermittency effects in ϵ. Such intermittency is characterized by the long‐range (i.e., power‐law decay) spatial auto‐correlation function of ϵ. That α o varies with a macro‐scale Reynolds number does not necessarily violate the micro‐eddy model. Instead, it points to a coordination between the macro‐ and micro‐scales arising from the transfer of energy across scales in the energy cascade.

Batchelor scale↗

Methods to Evaluate Subcolumn Profiles Based on Two-Point Diagnostics

In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.

54 ENVIRONMENTAL SCIENCES↗

Emergence of unprecedented climate change in projected future precipitation

The future time of emergence when precipitation changes due to anthropogenic influences begins to continuously exceed the previous maximum value is defined as the ‘tipping year’ Historical experiments and future experiments simulated by state-of-the-art climate models were utilized. A total of 510,000 time series from year 1856 to 2095 were generated by sampling the natural internal variability in precipitation. The time evolutions of internal variability in the whole time period were estimated from the combination of past and future experiments with preindustrial control experiments. A large ensemble size enabled an estimation of the probability density function of the tipping year at each grid point, providing precise information on the uncertainty of the projection. The tipping year of average precipitation emerges earlier in high latitudes than in lower latitudes. In some regions in lower latitudes and mid-latitudes, the tipping year of intense precipitation emerges faster than that of average precipitation. The tipping years of average and intense precipitation are earlier for higher anthropogenic forcing scenarios than for lower scenarios. The global average of the tipping year for intense precipitation might be attributed to the enhancement of the thermodynamic effect (moisture) rather than the dynamic effect (vertical motion).

54 ENVIRONMENTAL SCIENCES↗

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Learning the temporal evolution of multivariate densities via normalizing flows

In this work, we propose a method to learn multivariate probability distributions using sample path data from stochastic differential equations. Specifically, we consider temporally evolving probability distributions (e.g., those produced by integrating local or nonlocal Fokker–Planck equations). Here, we analyze this evolution through machine learning assisted construction of a time-dependent mapping that takes a reference distribution (say, a Gaussian) to each and every instance of our evolving distribution. If the reference distribution is the initial condition of a Fokker–Planck equation, what we learn is the time-T map of the corresponding solution. Specifically, the learned map is a multivariate normalizing flow that deforms the support of the reference density to the support of each and every density snapshot in time. We demonstrate that this approach can approximate probability density function evolutions in time from observed sampled data for systems driven by both Brownian and Lévy noise. We present examples with two- and three-dimensional, uni- and multimodal distributions to validate the method.

97 MATHEMATICS AND COMPUTING↗

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗

Lagrangian conditional statistics and flow topology in edge plasma turbulence

Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. Here, the density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantitative x ray phase contrast imaging of oblique shock wave–interface interactions

Oblique shock wave–interface interactions of gases and liquids have been extensively studied in shock tubes using optical methods to measure equation-of-state (EOS) parameters. However, this is difficult with solids due to their opaqueness to visible light. X ray phase contrast imaging (XPCI) has the penetrative strength to probe solids while still being sensitive to mass density and enhancing the visibility of material boundaries. We investigate the accuracy and repeatability of measuring the mean value of the average mass density (areal density divided by thickness) over region S (B S ) and flow deflection angle (θ) from XPCI images of a sample. To that end, a Hough transform-based method for measuring θ is developed. To measure B S , the XPCI image intensity probability density function (PDF) is modeled accounting for the spatial distribution of x ray energy, scintillator response, and pulse-to-pulse variation in the x ray intensity. In addition, a Monte Carlo-based algorithm for computing the B S PDF is developed. Both methods are validated on an impact-generated oblique shock wave interacting at a solid polymer-to-polymer interface. Furthermore this is accomplished through a comparison to hydrodynamic simulations using well-established EOS. Under the modeling framework for the XPCI image intensity, B S is computed with an accuracy of < 0.1% and precision of 3%–5%, while θ has an uncertainty of 0.2°, respectively. This shows that the XPCI-based model that is developed here could be an invaluable tool for high-fidelity testing of hydrodynamic models in shock polar configurations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Examination of probability distribution of mixture fraction in LES/FDF modelling of a turbulent partially premixed jet flame

An accurate prediction of the probability density function (PDF) of the mixture fraction is crucial to the prediction of combustion since mixing plays an important role in turbulent non-premixed and partially premixed flames. This work provides an assessment of the large-eddy simulation (LES)/filtered density function (FDF) method for the prediction of the PDF of the mixture fraction. The advantage of the LES/FDF method is that it provides the full predictions of the statistical distribution of scalars including the mixture fraction. The predictive accuracy of the method for the PDF is yet to be fully validated. Assessing the prediction of the PDF of the mixture fraction, a conserved scalar, is an important starting point. The Sydney/Sandia inhomogeneous inlet jet flame is used as a test case. A quick comparison shows that the LES/FDF predicted PDF shapes of the mixture fraction deviate significantly from the commonly presumed Beta-PDF as well as from the experimental data in the flame. Here, to examine the source of the discrepancy, we clarify the different PDF definitions used in the comparison among the predictions, measurements, and the presumed shape PDFs. The discrepancy observed from the comparison is largely reconciled by clarifying the difference between the PDFs that are examined. The PDF of the resolved mixture fraction is shown to be close to the Beta-PDF in both the measurements and predictions, while the PDF directly deduced from the LES/FDF particles deviates significantly from the Beta-PDF. A multimodal PDF analysis and a pseudo convergence analysis are conducted to provide plausible evidence to support the predicted multimodal PDF shapes. The sub-filter scale FDF is shown to be close to the Beta-PDF too through the construction of a synthesized PDF, which supports the common presumed Beta-PDF assumption used in the presumed PDF methods when combined with LES.

42 ENGINEERING↗

Transport and losses of fusion-born alpha particles in the presence of tearing modes using the new Toroidal Accelerated Particle Simulator (TAPaS)

Here, the transport and losses of fusion-born alpha particles are studied in the presence of a single-helicity tearing mode, characterized by (m = 2, n = 1). The analysis is performed by means of the recently developed Toroidal Accelerated Particle Simulator (TAPaS). Although such modes have usually been believed to result only in a local flattening of the radial profiles, it is shown that the density profile can exhibit a global modification leading to significant losses of alpha particles. This is due to the fact that, although the magnetic field does not exhibit any chaotic behaviour, the trajectories of alpha particles do, as revealed by their Poincaré maps. Such results are in qualitative agreement with past observations and simulations of energetic particles generated by neutral beam injection in TFTR, DIII-D and AUG tokamaks. In-depth analysis is carried out to characterize the impact of the tearing mode on the transport and losses of fusion-born alpha-particles with a realistic density profile. The impact of the amplitude is evidenced. Moreover, the effect of the island rotation frequency is assessed based on a detailed analysis of the linear resonances in phase-space, in agreement with the simulation results. Finally, the probability density function of the exit time has been computed and the transport of alpha particles has been found to be anomalous.

anomalous transport↗

Turbulence statistical analysis of the L-H transition and RMPs in KSTAR

Here, we investigate the turbulence statistics associated with low-to-high confinement (L-H) transitions and externally applied resonant magnetic perturbations (RMPs) in KSTAR. Time-series fluctuations of electron density n e , electron temperature T e , and the time derivative of the poloidal magnetic field dB θ /dt (Mirnov coils) are analysed using information-geometric measures (information rate Γ and information length $\mathcal{L}$ = ∫ Γ dt), together with kurtosis κ and variance σ 2 . In low-density upper single-null plasmas (n e ~ 1.2 x 10 19 m -3 ), a ~80 kHz magnetic mode coupling n e , T e , dB θ /dt and emerges prior to the L-H transition and persists into the edge-localised modes H-mode. Edge-localised RMPs (ERMPs) suppress this coherent mode but enhance intermittency, producing frequent bursts that abruptly reshape the time-dependent probability density functions (PDFs) and generate large spikes in Γ (with smaller changes in κ), signalling ERMP-driven departures from quasi-stationarity. The impact of ERMPs on background fluctuation levels depends on density, radial location, and the fluctuating variable itself ($\tilde{n}$, $\tilde{T}$, $\dot{B}$ θ ), whereas $\mathcal{L}$ provides a robust, regime-agnostic measure of cumulative statistical reorganisation and spatial decorrelation. In particular, at low density we observe weaker coupling between $\tilde{n}$ and $\tilde{T}$, along with a tendency toward decreased radial correlation-most clearly for $\tilde{T}$-under ERMPs. Overall, information geometry cleanly captures intermittent events, quantifies non-equilibrium PDF evolution, and offers a compact, cross-diagnostic metric for assessing resonant magnetic perturbation effects on edge transport and correlation across densities, radial locations, and confinement states.

Kim, Eun-jin [Coventry Univ. (United Kingdom); Seo↗

Disentangling core and edge mechanisms of the density limit in DIII-D negative triangularity plasmas

The density limit is investigated in the DIII-D negative triangularity plasmas which lack a standard H-mode edge. We find the limit may not be a singular disruptive boundary but a multifaceted density saturation phenomenon governed by distinct core and edge transport mechanisms. Sustained, non-disruptive operation is achieved at densities up to 1.8 times the Greenwald limit ($n_\mathrm{G}$) until the termination of auxiliary heating. Systematic power scans show distinct power scalings for the core ($n_e \propto P_\mathrm{SOL}^{0.27\pm0.03}$) and edge ($n_e \propto P_\mathrm{SOL}^{0.42\pm0.04}$) density limits. The edge density saturation is triggered by the onset of a non-disruptive, high-field side radiation front and the associated cooling, which clamps the edge density below . In contrast, the core density continues to rise until it saturates, a state characterized by enhanced core turbulence. Core transport evolves from a diffusive to an intermittent, avalanche-like state, as indicated by heavy-tailed probability density functions (kurtosis ≈ 6), increased Hurst exponents, and a $1/f$-type power spectrum. These findings suggest that the density limit in the low-confinement regime is determined by a combination of edge radiative cooling and core turbulent transport. This distinction provides separate targets for control strategies aimed at extending the operational space of future fusion devices.

density limit↗

Using intrusive approaches as a step towards accounting for stochasticity in wind turbine design

Current wind turbine design methods require tens of thousands of time-domain simulations and use different random seeds to account for the stochasticity of the environmental conditions. The account of stochasticity is nonintrusive because the sampling method calls a deterministic model multiple times without changing its underlying equations. In this work, we investigate and demonstrate using simple proof of concepts how intrusive approaches can be used to directly account for stochasticity in the equations representing a mechanical system. Our long term goal is to apply such methodology to the design of wind turbines without requiring an excessive number of simulations. Intrusive methods manipulate stochastic variables directly to provide the probability density functions (PDFs) of the states and outputs at any time as functions of the PDFs of the inputs. We illustrate how different methods can be used with a reduced-order model of a wind turbine with one degree of freedom and for linear and nonlinear models. We discuss how the methods can be extended and what it will take to apply them to a level of fidelity similar to current state-of-the-art wind turbine design tools.

17 WIND ENERGY↗

Stellar angular momentum distribution linked to galaxy morphology

ABSTRACT We study the spatially resolved stellar specific angular momentum j* in a high-quality sample of 24 Calar Alto Legacy Integral Field Area galaxies covering a broad range of visual morphology, accounting for stellar velocity and velocity dispersion. The shape of the spaxelwise probability density function of normalized s = j*/j*mean, PDF(s), deviates significantly from the near-universal initial distribution expected of baryons in a dark matter halo and can be explained by the expected baryonic effects in galaxy formation that remove and redistribute angular momentum. Further we find that the observed shape of the PDF(s) correlates significantly with photometric morphology, where late-type galaxies have a PDF(s) that is similar to a normal distribution, whereas early types have a strongly skewed PDF(s) resulting from an excess of low-angular momentum material. Galaxies that are known to host pseudo-bulges (bulge Sérsic index nb < 2.2) tend to have less skewed bulge PDF(s), with skewness (b1rb) ≲ 0.8. The PDF(s) encodes both kinematic and photometric information and appears to be a robust tracer of morphology. Its use is motivated by the desire to move away from traditional component-based classifications which are subject to observer bias, to classification on a galaxy’s fundamental (stellar mass and angular momentum) properties. In future, PDF(s) may also be useful as a kinematic decomposition tool.

Sweet, Sarah M.↗

Evaluation of probabilistic photometric redshift estimation approaches for The Rubin Observatory Legacy Survey of Space and Time (LSST)

Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced forLarge Synoptic Survey Telescope The Rubin Observatory Legacy Survey of Space and Time (lsst) Dark Energy Science Collaboration (desc). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. Additionally, we report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.

79 ASTRONOMY AND ASTROPHYSICS↗

The PDF perspective on the tracer-matter connection: Lagrangian bias and non-Poissonian shot noise

ABSTRACT We study the connection of matter density and its tracers from the probability density function (PDF) perspective. One aspect of this connection is the conditional expectation value 〈δtracer|δm〉 when averaging both tracer and matter density over some scale. We present a new way to incorporate a Lagrangian bias expansion of this expectation value into standard frameworks for modelling the PDF of density fluctuations and counts-in-cells statistics. Using N-body simulations and mock galaxy catalogues we confirm the accuracy of this expansion and compare it to the more commonly used Eulerian parametrization. For haloes hosting typical luminous red galaxies, the Lagrangian model provides a significantly better description of 〈δtracer|δm〉 at second order in perturbations. A second aspect of the matter-tracer connection is shot-noise, i.e. the scatter of tracer density around 〈δtracer|δm〉. It is well known that this noise can be significantly non-Poissonian and we validate the performance of a more general, two-parameter shot-noise model for different tracers and simulations. Both parts of our analysis are meant to pave the way for forthcoming applications to survey data.

Friedrich, Oliver↗