Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability density function control”

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 55 records · Page 3

Stochastic Framework for Optimal Control of Planetary Reentry Trajectories Under Multilevel Uncertainties

We present a novel stochastic optimal control framework that accounts for various types of uncertainties, with application to reentry trajectory planning. The formulation of the optimal trajectory control problem is presented in the context of an indirect method where a functional objective associated with the terminal vehicle speed is to be minimized. Uncertain input parameters in the optimal trajectory control model, including aerodynamic parameters and initial and terminal conditions, are modeled as aleatory random variables, while the statistical parameters of these aleatory distributions are themselves random variables. The parametric and model uncertainties are simultaneously propagated through an extended polynomial chaos expansion (EPCE) formalism. Several metrics are described to evaluate response statistics and presented as insightful tools for robust decision making. Specifically, the response probability density function (PDF) reflecting influence of both epistemic and aleatory uncertainties is obtained. By sampling over the random variables representing model error, an ensemble of response PDFs is generated and the associated failure probability is estimated as a random variable with its own polynomial chaos expansion. Besides, the sensitivity index functions of response PDF with respect to the statistical parameters are evaluated. Coupling parametric and model uncertainties within the EPCE framework leads to a robust and efficient paradigm for multilevel uncertainty propagation and PDF characterization in general optimal control problems.

Engineering↗

Analysis of homogeneous turbulent reacting flows

Full turbulence simulations at low Reynolds numbers were made for the single-step, irreversible, bimolecular reaction between non-premixed reactants in isochoric, decaying homogeneous turbulence. Various initial conditions for the scalar field were used in the simulations to control the initial scalar dissipation length scale, and simulations were also made for temperature-dependent reaction rates and for non-stoichiometric and unequal diffusivity conditions. Joint probability density functions (pdf's), conditional pdf's, and various statistical quantities appearing in the moment equations were computed. Preliminary analysis of the results indicates that compressive strain-rate correlates better than other dynamical quantities with local reaction rate, and the locations of peak reaction rates seem to be insensitive to the scalar field initial conditions.

Leonard, A. D.↗

The non-Gaussian joint probability density function of slope and elevation for a nonlinear gravity wave field

On the basis of the mapping method developed by Huang et al. (1983), an analytic expression for the non-Gaussian joint probability density function of slope and elevation for nonlinear gravity waves is derived. Various conditional and marginal density functions are also obtained through the joint density function. The analytic results are compared with a series of carefully controlled laboratory observations, and good agreement is noted. Furthermore, the laboratory wind wave field observations indicate that the capillary or capillary-gravity waves may not be the dominant components in determining the total roughness of the wave field. Thus, the analytic results, though derived specifically for the gravity waves, may have more general applications.

Huang, N. E.↗

Turbulence control on airborne laser platform

An avctive flow control device to generate large-scale, periodic structures in a turbulent shear flow is developed. Together with adaptive optics, the device may be used on airborne laser platforms to reduce or eliminate optical distortion caused by the turbulence in the aircraft's boundary layer. A flat plate towed in a water channel is used as a test bed. A cyclic jet issuing from a spanwise slot is used to collect the turbulent boundary layer for a finite time during its 'on' period. When the jet is turned 'off', all of the turbulent fluid is released instantaneously in one large eddy that convects downstream. Flow visualization and hot-film probe measurements are used together with pattern recognition algorithms to demonstrate the viability of the flow control method. The instantaneous velocity signal is used to compute important statistical quantities of the random velocity field, such as the mean, the root-mean-square, the spectral distribution, and the probability density function. When optimized for a given boundary layer, the cyclic jet produces periodic structures that are qualitatively similar to the random, naturally occurring ones. These structures seem to trigger the onset of bursting events near the wall. Thus, the present device generates periodic structures in both the outer and inner regions of a turbulent boundary layer.

Gad-El-hak, Mohamed↗

Establishing Models for Digital Twin of Hydropower Systems Using Probability Density Function Shaping

This paper introduces a digital twin modeling method for hydropower systems with Kaplan turbines using probability density function (PDF) shaping. We first use multilayer perceptron (MLP) model to build the discretized openloop Kaplan unit, where the MLP is trained by historical data. Then we use a proportional integral double derivative (PIDD) controller and a lead-lag exciter to test the obtained digital twin model in a closed-loop fashion. Simulation results show that the proposed digital twin modeling method can accurately capture the dynamics of the Kaplan hydropower unit. Finally, we show that the obtained digital twin can help to optimize the PIDD parameters. Compared with the original PIDD controller, the optimized one can achieve an over 90% improvement on the mean square tracking error.

Yin, Zhun [New York University]↗

A reaction mechanism for plasma electrolysis of AgNO 3 forming silver nanoclusters and nanoparticles

In plasma-driven solution electrolysis (PDSE), gas-phase plasma-produced species interact with an electrolytic solution to produce, for example, nanoparticles. An atmospheric pressure plasma jet (APPJ) directed onto a liquid solution containing a metallic salt will promote reduction of metallic ions in solution, generating metallic clusters that nucleate to form nanoparticles. In this article, results from a computational investigation are discussed of a PDSE process in which a radio-frequency APPJ sustained in helium impinges on a silver nitrate solution, resulting in growth of silver nanoparticles. A reaction mechanism was developed and implemented in a global plasma chemistry model to predict nanoparticle growth. To develop the reaction mechanism, density functional theory was used to generate probable silver growth pathways up to Ag 9 . Neutral clusters larger than Ag 9 were classified as nanoparticles. Kinetic reaction rate coefficients for thermodynamically favorable growth pathways were estimated based on an existing, empirically determined base reaction mechanism for smaller Ag particle interactions. These rates were used in conjunction with diffusion-controlled reaction rate coefficients that were calculated for other Ag species. The role of anions in reduction of Agn ions in forming nanoparticles is also discussed. Oxygen containing impurities or admixtures to the helium, air entrainment into the APPJ, and dissociation of saturated water vapor above the solution can produce additional reactive oxygen species in solution, resulting in the production of anions and O 2 – in particular. For a given molarity, delivering a sufficient fluence of reducing species will produce similar nanoparticle densities and sizes for all applied power levels. In conclusion, comparisons are made to alternate models for nanoparticle formation, including charged nanoparticles and use of direct current plasmas.

36 MATERIALS SCIENCE↗

Variation of fan tone steadiness for several inflow conditions

An amplitude probability density function analysis technique for quantifying the degree of fan noise tone steadiness has been applied to data from a fan tested under a variety of inflow conditions. The test conditions included typical static operation, inflow control by a honeycomb/screen device and forward velocity in a wind tunnel simulating flight. The ratio of mean square sinusoidal-to-random signal content in the fundamental and second harmonic tones was found to vary by more than an order-of-magnitude. Some implications of these results concerning the nature of fan noise generation mechanisms are discussed.

Balombin, J. R.↗

Variation of fan tone steadiness for several inflow conditions

An amplitude probability density function analysis technique for quantifying the degree of fan noise tone steadiness has been applied to data from a fan tested under a variety of inflow conditions. The test conditions included typical static operation, inflow control by a honeycomb/screen device and forward velocity in a wind tunnel simulating flight. The ratio of mean square sinusoidal-to-random signal content in the fundamental and second harmonic tones was found to vary by more than an order-of-magnitude. Some implications of these results concerning the nature of fan noise generation mechanisms are discussed.

Balombin, J. R.↗

CryoFold: Determining protein structures and data-guided ensembles from cryo-EM density maps

Cryoelectron microscopy requires molecular modeling for refinement of structures. Ensemble models arrive at low free-energy molecular structures, but are computationally expensive and limited to resolving only small proteins. Here, we introduce CryoFold, a pipeline of molecular dynamics simulations that determines ensembles of protein structures by integrating density data of varying sparsity at 3–5 Å resolution with sequence information and coarse-grained topological knowledge of the protein folds. We present six examples, folding proteins between 72 and 2,000 residues, including large membrane and multi-domain systems, and results from two Electron Microscopy Data Bank (EMDB) competitions. Driven by data from a single state, CryoFold discovers ensembles of common low-energy models together with rare low-probability structures that capture the equilibrium distribution of proteins constrained by the density maps. Many of these conformations are experimentally validated and functionally relevant. We arrive at a set of best practices for data-guided protein folding that are controlled using a Python graphical user interface (GUI).

59 BASIC BIOLOGICAL SCIENCES↗

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↗

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗

Understanding the Microphysical Control and Spatial‐Temporal Variability of Warm Rain Probability Using CloudSat and MODIS Observations

Abstract By combining measurements from MODIS and the CloudSat radar, we develop a parameterization scheme to quantify the combined microphysical controls by liquid water path (LWP) and cloud droplet number concentration (CDNC) of the probability of precipitation (PoP) in marine low cloud over tropical oceans. We demonstrate that the spatial‐temporal variation of grid‐mean in‐cloud can be largely explained by the variation of the joint probability density function of LWP and CDNC in the phase space specified by the bivariate PoP (LWP and CDNC) function. Through a series of sensitivity tests guided by this understanding, we find that in the Southeastern Pacific and Atlantic the stratocumulus to cumulus transition of the is mainly due to the variation of CDNC while the annual cycle is mainly due to the variation of LWP. The results of this study provide a viable way to diagnose the root cause of warm rain problems in global climate models.

54 ENVIRONMENTAL SCIENCES↗

Turbulence control on an airborne laser platform

An active flow control device to generate large-scale, periodic structures in a turbulent shear flow is developed. Together with adaptive optics, the device may be used on airborne laser platforms to reduce or eliminate optical distortion caused by the turbulence in the aircraft's boundary layer. A cyclic jet issuing from a spanwise slot is used to collect the turbulent boundary layer for a finite time and then release all of the flow instantaneously in one large eddy that convects downstream. Flow visualization and hot-film probe measurements are used together with pattern recognition algorithms to demonstrate the viability of the flow control method. A flat plate towed in a water channel is used as a test bed. The instantaneous velocity signal is used to compute important statistical quantities of the random velocity field, such as the mean, the root-mean-square, the spectral distribution, and the probability density function. When optimized for a given boundary layer, it is shown that the cyclic jet will produce periodic structures that are similar to the random, naturally occurring ones. These structures seem to trigger the onset of bursting events near the wall of the plate. Thus, the present device generates periodic structures in both the outer and inner regions of a turbulent boundary layer.

Gad-El-hak, Mohamed↗

Remote-sensing based technique to account for sub-grid scale variability of land surface properties

A method has been presented for the representation of sub-grid scale variability of surface properties within a land surface processes model. The method uses remotely-sensed data to directly or indirectly estimate probability density functions (PDF's) or key surface variables. Application of this technique in a coupled land surface-atmosphere model requires only grid-scale values of the variables of interest, obtained from low-resolution satellite imagery or surface/remote sensing data assimilation. The PDF's of each controlling surface property are superimposed on the respective grid-scale values to simulate sub-grid scale heterogeneity. Sensitivity studies will be carried out to ascertain the relative importance of the heterogeneity of several variables, and the degree to which non-linear property-process interactions impact large-scale fluxes.

Crosson, William L.↗

Understanding the Microphysical Control and Spatial-Temporal Variability of Warm Rain Probability Using CloudSat and MODIS Observations

By combining measurements from MODIS and the CloudSat radar, we develop a parameterization scheme to quantify the combined microphysical controls by liquid water path (LWP) and cloud droplet number concentration (CDNC) of the probability of precipitation (PoP) in marine low cloud over tropical oceans. We demonstrate that the spatial-temporal variation of grid-mean in-cloud can be largely explained by the variation of the joint probability density function of LWP and CDNC in the phase space specified by the bivariate PoP (LWP and CDNC) function. Through a series of sensitivity tests guided by this understanding, we find that in the Southeastern Pacific and Atlantic the stratocumulus to cumulus transition of the is mainly due to the variation of CDNC while the annual cycle is mainly due to the variation of LWP. The results of this study provide a viable way to diagnose the root cause of warm rain problems in global climate models.

Zhibo Zhang↗

Theoretic aspects of the identification of the parameters in the optimal control model

The identification of the parameters of the optimal control model from input-output data of the human operator is considered. Accepting the basic structure of the model as a cascade of a full-order observer and a feedback law, and suppressing the inherent optimality of the human controller, the parameters to be identified are the feedback matrix, the observer gain matrix, and the intensity matrices of the observation noise and the motor noise. The identification of the parameters is a statistical problem, because the system and output are corrupted by noise, and therefore the solution must be based on the statistics (probability density function) of the input and output data of the human operator. However, based on the statistics of the input-output data of the human operator, no distinction can be made between the observation and the motor noise, which shows that the model suffers from overparameterization.

Vanwijk, R. A.↗

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↗

Preliminary results in using Deep Learning to emulate BLOB, a nuclear interaction model

Purpose: A reliable model to simulate nuclear interactions is fundamental for Ion-therapy. We already showed how BLOB (“Boltzmann-Langevin One Body”), a model developed to simulate heavy ion interactions up to few hundreds of MeV/u, could simulate also 12 C reactions in the same energy domain. However, its computation time is too long for any medical application. For this reason we present the possibility of emulating it with a Deep Learning algorithm. Methods: The BLOB final state is a Probability Density Function (PDF) of finding a nucleon in a position of the phase space. We discretised this PDF and trained a Variational Auto-Encoder (VAE) to reproduce such a discrete PDF. As a proof of concept, we developed and trained a VAE to emulate BLOB in simulating the interactions of 12 C with 12 C at 62 MeV/u. To have more control on the generation, we forced the VAE latent space to be organised with respect to the impact parameter (b) training a classifier of b jointly with the VAE. Results: In this work, the distributions obtained from the VAE are similar to the input ones and the computation time needed to use the VAE as a generator is negligible. Conclusions: We show that it is possible to use a Deep Learning approach to emulate a model developed to simulate nuclear reactions in the energy range of interest for Ion-therapy. We foresee the implementation of the generation part in C++ and to interface it with the most used Monte Carlo toolkit: Geant4.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗