Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Development and validation of the cavitation-induced erosion risk assessment tool

This work presents the development of a cavitation-induced erosion risk assessment (CIERA) tool that links multiphase flow simulation predictions with the progress towards material erosion. To develop a robust erosion modeling tool, the cavitation and erosion predictions for pressurized diesel fuel flow within a channel geometry were validated over a range of Reynolds and cavitation number conditions in two different aluminum channel geometries, one featuring a rounded inlet corner and the other with a sharp inlet corner. The multiphase flow development within the channel was modeled using a compressible mixture model, where phase change was represented with the homogeneous relaxation model and the turbulent flow evolution was modeled using a dynamic structure approach for Large Eddy Simulations. To improve representation of the incubation period before material rupture over existing approaches, a physics-based metric was derived based on the cumulative energy absorbed by the solid material from repeated hydrodynamic impacts. When the average peak pressure was related to the incubation period, the incubation period and its sensitivity to changes in flow conditions were found to be overpredicted. In contrast, predictions from CIERA provided a more accurate means to qualitatively and quantitatively predict the influence of flow conditions on the incubation period before material erosion. When the predicted stored energy was related to the solid material properties to estimate the incubation period, multiphase flow simulations demonstrated accurate representation of the sensitivity of erosion severity to changes in flow conditions. The use of CIERA led to quantitative agreement of the predicted incubation period within 5% of the experimentally measured incubation period.

42 ENGINEERING↗
DOE OSTI2020

Thermal Overloading Risk Mitigation With a Semi-Analytical Probabilistic Model on Branch Current

A semi-analytical formulation is presented in this paper for the probability computation of branch current in multiphase systems. The developed formula is derived based on the linear power flow model in rectangular coordinates. The system uncertainty injections can be renewable energy resources or loads and are modeled using a Gaussian mixture model (GMM). The developed formula can be used to compute the line current violation probability as well as integrate into optimal power flow problem as chance-constraint relaxation. Here, the proposed formula is first compared with the Matlab embedded numerical integration function to show its performance. Besides, the semi-analytical formula is validated and compared with the Monte Carlo simulation method using the IEEE 123-bus system, EPRI Ckt5, and Ckt7 systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1983

Establishment of a center of excellence for applied mathematical and statistical research

The state of the art was assessed with regards to efforts in support of the crop production estimation problem and alternative generic proportion estimation techniques were investigated. Topics covered include modeling the greeness profile (Badhwarmos model), parameter estimation using mixture models such as CLASSY, and minimum distance estimation as an alternative to maximum likelihood estimation. Approaches to the problem of obtaining proportion estimates when the underlying distributions are asymmetric are examined including the properties of Weibull distribution.

Woodward, W. A.↗
NASA NTRS1985

Geologic mapping using Landsat MSS and TM images - Removing vegetation by modeling spectral mixtures

Semi-arid regions of the island of Hawaii and of the Tucson Mountains, AZ, having less than 50 percent green vegetation cover were mapped into rock/soil types by first identifying spectral mixtures of vegetation and rock/soil, and then removing the vegetation information. Landsat MSS images were used for Hawaii, and a Landsat TM image was used for the Tucson Mountains. Following analysis of the images field transects were done in test areas to determine percent vegetation cover and vegetation community characteristics, and to observe rock and soil types. For both Hawaii and the Tucson Mountains 20 percent or more of green vegetation caused potential misidentification of the rock/soil units unless spectral mixing was taken into account.

Adams, J. B.↗
NASA NTRS2012

Using Statistical Multivariable Models to Understand the Relationship Between Interplanetary Coronal Mass Ejecta and Magnetic Flux Ropes

In-situ measurements of interplanetary coronal mass ejections (ICMEs) display a wide range of properties. A distinct subset, "magnetic clouds" (MCs), are readily identifiable by a smooth rotation in an enhanced magnetic field, together with an unusually low solar wind proton temperature. In this study, we analyze Ulysses spacecraft measurements to systematically investigate five possible explanations for why some ICMEs are observed to be MCs and others are not: i) An observational selection effect; that is, all ICMEs do in fact contain MCs, but the trajectory of the spacecraft through the ICME determines whether the MC is actually encountered; ii) interactions of an erupting flux rope (PR) with itself or between neighboring FRs, which produce complex structures in which the coherent magnetic structure has been destroyed; iii) an evolutionary process, such as relaxation to a low plasma-beta state that leads to the formation of an MC; iv) the existence of two (or more) intrinsic initiation mechanisms, some of which produce MCs and some that do not; or v) MCs are just an easily identifiable limit in an otherwise corntinuous spectrum of structures. We apply quantitative statistical models to assess these ideas. In particular, we use the Akaike information criterion (AIC) to rank the candidate models and a Gaussian mixture model (GMM) to uncover any intrinsic clustering of the data. Using a logistic regression, we find that plasma-beta, CME width, and the ratio O(sup 7) / O(sup 6) are the most significant predictor variables for the presence of an MC. Moreover, the propensity for an event to be identified as an MC decreases with heliocentric distance. These results tend to refute ideas ii) and iii). GMM clustering analysis further identifies three distinct groups of ICMEs; two of which match (at the 86% level) with events independently identified as MCs, and a third that matches with non-MCs (68 % overlap), Thus, idea v) is not supported. Choosing between ideas i) and iv) is more challenging, since they may effectively be indistinguishable from one another by a single in-situ spacecraft. We offer some suggestions on how future studies may address this.

Riley, P.↗
DOE OSTI2022

Unsupervised probabilistic models for sequential Electronic Health Records

We develop an unsupervised probabilistic model for heterogeneous Electronic Health Record (EHR) data. Utilizing a mixture model formulation, our approach directly models sequences of arbitrary length, such as medications and laboratory results. This allows for subgrouping and incorporation of the dynamics underlying heterogeneous data types. The model consists of a layered set of latent variables that encode underlying structure in the data. These variables represent subject subgroups at the top layer, and unobserved states for sequences in the second layer. We train this model on episodic data from subjects receiving medical care in the Kaiser Permanente Northern California integrated healthcare delivery system. The resulting properties of the trained model generate novel insight from these complex and multifaceted data. In addition, we show how the model can be used to analyze sequences that contribute to assessment of mortality likelihood.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Systems Engineering of Rhodococcus opacus to Enable Production of Drop-in Fuels from Lignocellulose

Production of drop-in fuels from lignocellulose using Rhodococcus opacus PD630 (hereafter R. opacus) is a challenging goal. During the grant period we have pushed the field forward significantly in several areas of research. Towards the end goal of accelerating the adoption of R. opacus in biofuel production, during the grant period we have expanded the phenotypic characterization of R. opacus grown in single aromatic (model lignocellulosic) compounds or their mixtures, modeling the growth conditions in lignocellulosic biomass. Harnessing the power of adaptive evolution, we produced evolved R. opacus isolates with superior lignin valorization capabilities and identified differentially expressed genes and pathways after adaptation. We used next generation multi-omic techniques such as genomic, transcriptomic, and metabolomic analyses, to identify the catabolic pathways used by R. opacus to degrade aromatic compounds and funnel these degradation products into central metabolism, as well as the aromatic transport genes required for increased tolerance and utilization. Taking this information one step further, we identified endogenous transcription factors and regulatory mechanisms important for degradation of five model aromatic compounds. To accurately estimate R. opacus growth and consumption on model lignin compounds we pioneered the use of novel extraction procedures prior to GC-MS analysis. Alongside 13 C-metabolic flux analysis, we have elucidated the metabolic routes preferred by Rhodococcus opacus during aromatic compound degradation. Finally, we used in tandem lipidomics and high-resolution mass spectrometry to identify the modulation of mycolic acids and phospholipid membrane composition modification as a strategy for aromatic tolerance in R. opacus. Being a non-model organism, R. opacus lacks the breadth of tools and technical foundation which drive biofuel research in more well-understood microbes such as Escherichia coli. To reduce this burden for use, we designed and produced new tools for genomic manipulation and engineering in R. opacus. These engineering breakthroughs support efficient genomic editing, enabling gene overexpression, repression, and genetic alteration. Using these tools, we have generated synthetically engineered strains with increased lipogenesis and growth, both positive traits required for increased lignin valorization. Optimizing engineered strains for biofuel production from lignocellulose requires extremely sophisticated synthetic rewiring of metabolism. To facilitate systems-level reorganization of metabolism in R. opacus, we created a genome-scale model that accurately predicts metabolic flux and growth rates on the aromatic compound phenol. Lignin requires extensive pre-treatment before biological degradation by R. opacus. Towards an eventual goal of degrading real-world lignin, we developed new depolymerization processes to generate lignin breakdown products (LBP). We optimized LBP storage and composition analysis techniques, enabling accurate prediction of specific LBP compound integration into cell wall components. Overall, through the work funded by this grant we generated 20 manuscripts (17 published, 3 in review/preparation), methods for increased accuracy in metabolomics of aromatic compounds, multiple genetic tools for altering the R. opacus genome, genome scale models for predicting flux through metabolic pathways, as well as multi-omic data for community use. The work funded by this grant has increased the knowledge of aromatic degradation in bacteria and advanced our efforts to optimize R. opacus for lignin valorization.

09 BIOMASS FUELS↗
NASA NTRS1976

Numerical analysis of a two-dimensional nonsteady detonations

In the present work a system of two-dimensional nonsteady hydrodynamic and chemical kinetic equations was numerically integrated for an exothermic system. Assumed two-step reaction model simulates practically an oxyhydrogen mixture. The calculation starts from a plane Chapman-Jouguet detonation as an initial condition. Two-dimensional disturbances are generated by artificially placing nonuniformities ahead of the detonation front. Regardless of the difference of the given initial disturbances, a fixed number of triple shock waves were produced for a fixed combination of mixture model and geometry when the transition period was over. This shows that for a given detonation tube geometry any exothermic system has its own characteristic multidimensional structure. The obtained number of triple shock waves contained in the detonation front was in agreement with existing experimental observations under the same condition.

Taki, S.↗
DOE OSTI2020

Nonlinear sparse Bayesian learning for physics-based models

This paper addresses the issue of overfitting while calibrating unknown parameters of over-parameterized physics-based models with noisy and incomplete observations. Here, a semi-analytical Bayesian framework of nonlinear sparse Bayesian learning (NSBL) is proposed to identify sparsity among model parameters during Bayesian inversion. NSBL offers significant advantages over machine learning algorithm of sparse Bayesian learning (SBL) for physics-based models, such as 1) the likelihood function or the posterior parameter distribution is not required to be Gaussian, and 2) prior parameter knowledge is incorporated into sparse learning (i.e. not all parameters are treated as questionable). NSBL employs the concept of automatic relevance determination (ARD) to facilitate sparsity among questionable parameters through parameterized prior distributions. The analytical tractability of NSBL is enabled by employing Gaussian ARD priors and by building a Gaussian mixture-model approximation of the posterior parameter distribution that excludes the contribution of ARD priors. Subsequently, type-II maximum likelihood is executed using Newton's method whereby the evidence and its gradient and Hessian information are computed in a semi-analytical fashion. We show numerically and analytically that SBL is a special case of NSBL for linear regression models. Subsequently, a linear regression example involving multimodality in both parameter posterior pdf and model evidence is considered to demonstrate the performance of NSBL in cases where SBL is inapplicable. Next, NSBL is applied to identify sparsity among the damping coefficients of a mass-spring-damper model of a shear building frame. These numerical studies demonstrate the robustness and efficiency of NSBL in alleviating overfitting during Bayesian inversion of nonlinear physics-based models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Spatial patterns in occupancy and density of larval lampreys in freshwater habitats restored to a Stage 0 condition

Abstract We examined occupancy and density of larval lampreys ( Entosphenus tridentatus and Lampetra spp.) in two rivers in Oregon (USA) restored to a Stage 0 condition 1–5 years prior, using a multiscale occupancy model and a zero‐inflated Poisson mixture model. We sampled lampreys using backpack electrofishing in randomly distributed, paired, 1‐m 2 quadrats and recorded environmental data. Probabilities of occupancy and density were higher when water velocity was low, the substrate was noncompacted, and sediment was dominated by fines (<4 mm). At mean water depth (0.34 m) and velocity (0.09 m/s), estimated densities in occupied quadrats were 4.8 lampreys/m 2 (95%: 3.4–6.9) when the substrate was compacted, and fines were not dominant, and 21.1 lampreys/m 2 (95%: 17.7–25.3) when the substrate was noncompacted and fines were dominant. Probabilities of detecting occupancy in a 1‐m 2 quadrat sampled by backpack electrofishing were 0.76 (95%: 0.64–0.87) when captured after visual observation and 0.80 (95%: 0.71–0.88) with blind sweeps (i.e., constantly moving the net regardless of observation). The probability of capturing a single lamprey in a quadrat sampled by blind sweeps was 0.32 (95%: 0.27–0.37). Sampling in paired 1‐m 2 quadrats facilitated concurrent examination of patterns in occupancy and density while accounting for capture probability, which could aid temporal monitoring of restored habitats. To the best of our knowledge, this is the first study to document occupancy and estimate densities of larval lampreys in habitats that underwent valley floor restoration to Stage 0. We observed both lamprey genera within 5 years of restoration. Aquatic restoration that increases low‐velocity, noncompacted, fine sediment habitats could benefit lampreys.

Harris, Julianne E.↗
NASA NTRS1974

An approximate inviscid radiating flow-field analysis for sphere-cone Venusian entry vehicles

A computational method based on the earlier work of one of the authors has been improved and modified to provide an accurate and rapid prediction of the radiating flow about blunt conical probes entering the atmosphere of Venus. Special features include a more general approximate equation of state, a more rapidly converging iteration scheme, a 38-step absorption coefficient model for a Venus atmosphere gas mixture model, and a shock-shape formula which yields a computed body shape in excellent agreement with the desired body. The method was used to calculate distributions of radiative heating for a large spherically blunted conical probe vehicle (1.219 m base diameter) and for three small spherically blunted conical probe vehicles (0.724 m base diameter). Solutions obtained with the presently developed approximate method required short computing times of the order of 1 to 3 minutes for each complete solution and gave agreements of 10% to 15% with the more exact predictions.

Falanga, R. A.↗
NASA NTRSDate not supplied

Cryogenic Extension of NASA Species Polynomials Using Hydrogen and Oxygen at Stoichiometry

NASA has been conducting research in Rotating Detonation Rocket Engines (RDREs) for several years. A recent test at the Marshall Spaceflight Center has successfully shown operation of a RDRE that ran for 251 seconds which had combustion chamber cryogenic inlet conditions [1]. However, an overall pressure gain was not observed in the performance of the engine. Such an indication could mean an optimal chamber design is yet to be discovered. Such a design could significantly reduce parasitic pressure loses and improve overall engine efficiency beyond that of conventional rockets. To investigate different chamber configurations, new CFD capability is needed to accurately capture cryogenic thermophysical and transport property data. Such an effort is currently part of NASA’s Early Career Initiative (ECI) program. This paper focuses on the thermophysical property mixture model approach that makes use of NASA’s polynomial fits of species and another cryogenic model chosen from NIST’s Refprop program. Presented in this paper are results of a 1D denotation CFD simulation of stoichiometric Hydrogen and Oxygen mixture at a cryogenic upstream condition.

Combustion↗
NASA NTRSDate not supplied

Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching

Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.

line chilldown↗
NASA NTRS2007

Identification and Control of Aircrafts using Multiple Models and Adaptive Critics

We compared two possible implementations of local linear models for control: one approach is based on a self-organizing map (SOM) to cluster the dynamics followed by a set of linear models operating at each cluster. Therefore the gating function is hard (a single local model will represent the regional dynamics). This simplifies the controller design since there is a one to one mapping between controllers and local models. The second approach uses a soft gate using a probabilistic framework based on a Gaussian Mixture Model (also called a dynamic mixture of experts). In this approach several models may be active at a given time, we can expect a smaller number of models, but the controller design is more involved, with potentially better noise rejection characteristics. Our experiments showed that the SOM provides overall best performance in high SNRs, but the performance degrades faster than with the GMM for the same noise conditions. The SOM approach required about an order of magnitude more models than the GMM, so in terms of implementation cost, the GMM is preferable. The design of the SOM is straight forward, while the design of the GMM controllers, although still reasonable, is more involved and needs more care in the selection of the parameters. Either one of these locally linear approaches outperform global nonlinear controllers based on neural networks, such as the time delay neural network (TDNN). Therefore, in essence the local model approach warrants practical implementations. In order to call the attention of the control community for this design methodology we extended successfully the multiple model approach to PID controllers (still today the most widely used control scheme in the industry), and wrote a paper on this subject. The echo state network (ESN) is a recurrent neural network with the special characteristics that only the output parameters are trained. The recurrent connections are preset according to the problem domain and are fixed. In a nutshell, the states of the reservoir of recurrent processing elements implement a projection space, where the desired response is optimally projected. This architecture trades training efficiency by a large increase in the dimension of the recurrent layer. However, the power of the recurrent neural networks can be brought to bear on practical difficult problems. Our goal was to implement an adaptive critic architecture implementing Bellman s approach to optimal control. However, we could only characterize the ESN performance as a critic in value function evaluation, which is just one of the pieces of the overall adaptive critic controller. The results were very convincing, and the simplicity of the implementation was unparalleled.

Principe, Jose C.↗
NASA NTRS2014

Combined Uncertainty and A-Posteriori Error Bound Estimates for CFD Calculations: Theory and Implementation

Simulation codes often utilize finite-dimensional approximation resulting in numerical error. Some examples include, numerical methods utilizing grids and finite-dimensional basis functions, particle methods using a finite number of particles. These same simulation codes also often contain sources of uncertainty, for example, uncertain parameters and fields associated with the imposition of initial and boundary data,uncertain physical model parameters such as chemical reaction rates, mixture model parameters, material property parameters, etc.

CFD↗
DOE OSTI2024

Model and remote-sensing-guided experimental design and hypothesis generation for monitoring snow-soil–plant interactions

In this study, we develop a machine-learning (ML)-enabled strategy for selecting hillslope-scale ecohydrological monitoring sites within snow-dominated mountainous watersheds, with a particular focus on snow-soil–plant interactions. Data layers rely on spatial data layers from both remote sensing and hydrological model simulations. Specifically, a Landsat-based foresummer drought sensitivity index is used to define the dependency of the annual peak plant productivity on the Palmer drought severity index in the early growing season. Hydrological simulations provide the spatiotemporal dynamics of near-surface soil moisture and snow depth. In this framework, a regression analysis identifies the key hydrological variables relevant to the spatial heterogeneity of drought sensitivity. We then apply unsupervised clustering to these key variables, using the Gaussian mixture model, to group hillslopes into several zones that have divergent relationships regarding soil moisture, snow dynamics, and drought sensitivity. Using the datasets collected in the East River Watershed (Crested Butte, Colorado, United States), results show that drought sensitivity is significantly correlated with model-derived soil moisture and snow-free timing over space and time. The relationship is, however, non-linear, such that the correlation decreases above a threshold elevation and in a heavy snow year due to large snowpacks, lateral flow, and soil storage limitations. Clustering is then able to define the zones that have high or low sensitivity to drought, as well as the mid-elevation regions where sensitivity is associated with the topographic aspect and net potential radiation. In addition, the algorithm identifies the most representative hillslopes with road/trail access within each zone for installing monitoring sites. Our method also aims to significantly increase the use of ML and model-simulation results to guide critical zone and watershed monitoring activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Statistical Behavior of Low-Amplitude Power System Point-on-Wave Measurements

The power grid is undergoing massive changes to ensure resiliency and reliability in a more decentralized world. Distributed energy resources are becoming a prominent source of generation, potentially leading to a lack of centralized generation sources. Due to these new behaviors and system topologies, it is important to install measurement devices that are 1) accurate and 2) self-aware of their measurement quality. In this paper, a residential-scale microgrid is used to generate voltage and current waveforms, captured by Verivolt and National Instruments measurement equipment. A least-squares approach is used to separate the “clean” signals from the noise. Finally, Gaussian mixture modeling is used to approximate noise distributions, and it is shown these higher-order distribution estimates are a better fit to voltage and current noise profiles than single-mode Gaussian estimates.

24 POWER TRANSMISSION AND DISTRIBUTION↗