A second-order system model for pneumatic instrumentation lines.
Pneumatic instrumentation line transfer function approximated by system model for simplifying reactor control system analysis
Samuelson, R. D.↗
SEARCH · Engineering Papers
Results for “functional 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.
Specialization Restricts the Evolutionary Paths Available to Yeast Sugar Transporters
Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence–function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ~400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.
59 BASIC BIOLOGICAL SCIENCES↗
Kullback-Leibler information function and the sequential selection of experiments to discriminate among several linear models
A sequential adaptive experimental design procedure for a related problem is studied. It is assumed that a finite set of potential linear models relating certain controlled variables to an observed variable is postulated, and that exactly one of these models is correct. The problem is to sequentially design most informative experiments so that the correct model equation can be determined with as little experimentation as possible. Discussion includes: structure of the linear models; prerequisite distribution theory; entropy functions and the Kullback-Leibler information function; the sequential decision procedure; and computer simulation results. An example of application is given.
Sidik, S. M.↗
An Empirical Model For Intrinsic Alignments: Insights From Cosmological Simulations
We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce new, simulation-based predictions for IA and galaxy clustering. Our model is probabilistically formulated, and rests upon the assumption that the orientations of galaxies exhibit a correlation with their host dark matter (sub)halo orientation or with their position within the halo. We examine the necessary components and phenomenology of such a model by considering the alignments between (sub)halos in a cosmological dark matter only simulation. We then validate this model for a realistic galaxy population in a set of simulations in the Illustris-TNG suite. We create an HOD mock with Illustris-like correlations using our method, constraining the associated IA model parameters, with the $\mathcal{X}$$^{2}_{dof}$ between our model’s correlations and those of Illustris matching as closely as 1.4 and 1.1 for orientation–position and orientation–orientation correlation functions, respectively. By modeling the misalignment between galaxies and their host halo, we show that the 3-dimensional two-point position and orientation correlation functions of simulated (sub)halos and galaxies can be accurately reproduced from quasi-linear scales down to 0.1 $\mathcal{h}$ –1 Mpc. We also find evidence for environmental influence on IA within a halo. Our publicly-available software provides a key component enabling efficient determination of Bayesian posteriors on IA model parameters using observational measurements of galaxy-orientation correlation functions in the highly nonlinear regime.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
The total microbiome functions in bacteria and fungi
Unveiling the relationship between phylogeny and function of the microbiome is crucial to determine its contribution to ecosystem functioning. However, while there is a considerable amount of information on microbial phylogenetic diversity, our understanding of its relationship to functional diversity is still scarce. Here we predicted the total microbiome functions of bacteria and fungi on Earth using the total known functions from level 3 of KEGG Orthology by modelling the increase of functions with increasing diversity of bacteria or fungi. For bacteria and fungi, the unsaturated model described the data significantly better (for both P?<2.2e-16), suggesting the presence of two types of functions. Widespread functions ubiquitous in every living organism that make up two thirds of our current knowledge of microbiome functions are separated from rare functions from specialised enzymes present in only a few species. Given previous estimates on species richness, we predicted a global total of 35.5 million functions in bacteria and 3.2 million in fungi; of which only 0.02% and 0.14% are known today. Our approach highlights the necessity of novel and more sophisticated methods to unveil the entirety of rare functions to fully understand the involvement of the microbiome in ecosystem functioning. SIGNIFICANCE: The functionality of and within a microbial community is generally inferred based on the taxonomic annotation of the organism. However, our understanding of functional diversity and how it relates to taxonomy is still limited. Here we predict the total microbiome functionality in bacteria and fungi on Earth using known and annotated protein-coding sequences in species accumulation curves. Our estimates reveal that the majority of functionality (>99%) could be assigned to yet unknown and rare functions, highlighting that our current knowledge is incomplete and functional inference is thus lackluster.
Bacteria↗
Explainable and Differentiable Reinforcement Learning for Multi-objective Optimization in Particle Accelerators
Operating particle accelerators involves optimizing multiple goals simultaneously, which can be challenging due to trade-offs among objectives. While evolutionary algorithms like the genetic algorithm (GA) have been used for various Multi-Objective Optimization (MOO) tasks, they are not inherently suited for complex control problems. This talk highlights two variations of Reinforcement Learning (RL) for concurrently optimizing heat load and trip rates at the Continuous Electron Beam Accelerator Facility (CEBAF). The problem involves strict constraints on individual states, actions, and overall energy requirements of the beam. First, this talk highlights how differentiability can be harnessed through a Deep Differentiable Reinforcement Learning (DDRL) approach to address MOO issues within particle accelerators. We examine the DDRL method alongside Model Free Reinforcement Learning (MFRL), GA, and Bayesian Optimization (BO). The performance of these methods is assessed by generating a Pareto-front for two objectives. Our findings indicate that DDRL excels in handling high-dimensional problems more effectively than MFRL, BO, and GA. Next, we will show integration of explainable physics-based constraints into RL algorithms to enhance trans- parency and trust in decision-making processes by enabling users to verify that agents adhere to established physical principles. This surrogate function can be modeled using neural networks or sparse dictionary mod- els. By examining the mathematical form of the learned constraint function, we are able to confirm the agent has learned to use the established physics of each environment provided but the surrogate model. In addi- tion, we find that the introduction of a mathematical functional dictionary based surrogate model enables our reinforcement learning algorithms to reliably converge for difficult high-dimensional accelerator controls environments.
Rajput, Kishansingh [Thomas Jefferson National Acc↗
Efficient Calibration of Expensive Computational Models
Accounting for uncertainty when calibrating expensive computational models is a common challenge faced by scientists and engineers. Often Bayesian techniques are adopted to estimate a probability density function over the model parameters given noisy empirical data. The methods used to perform this type of probabilistic calibration are computationally prohibitive in that they require a large number of evaluations of the expensive model. In these cases, surrogate modeling -- that is, using a fast-to-evaluate, lower fidelity stand-in for the original computational model -- may be the only option to alleviate this computational burden. However, the upfront cost of generating training data to build a surrogate model can itself be expensive. As such, it is important to be judicious when selecting training points at which the full-fidelity model is evaluated. Here, an active learning approach is proposed that enables efficient selection of training points using approximate samples of the calibrated parameter probability density function. In this way, the training points can be concentrated in regions where the calibration algorithm requires high model accuracy.
active learning↗
PDF approach for compressible turbulent reacting flows
The objective of the present work is to develop a probability density function (pdf) turbulence model for compressible reacting flows for use with a CFD flow solver. The probability density function of the species mass fraction and enthalpy are obtained by solving a pdf evolution equation using a Monte Carlo scheme. The pdf solution procedure is coupled with a compressible CFD flow solver which provides the velocity and pressure fields. A modeled pdf equation for compressible flows, capable of capturing shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed, and an averaging procedure is developed to provide smooth Monte-Carlo solutions to ensure convergence. Two supersonic diffusion flames are studied using the proposed pdf model and the results are compared with experimental data; marked improvements over CFD solutions without pdf are observed. Preliminary applications of pdf to 3D flows are also reported.
Hsu, A. T.↗
PHASM: A Toolkit for Creating AI Surrogate Models within Legacy Codebases
PHASM (“Parallel Hardware viA Surrogate Models”) is a software toolkit for creating AI-based surrogate models of scientific code. AI-based surrogate models are widely used for creating fast and inverse simulations. PHASM anticipates an additional future use case: adapting legacy code to modern hardware. While data centers are investing in heterogeneous hardware such as GPUs and FPGAs, many established scientific codebases remain unable to take advantage of the hardware’s higher parallelism without undergoing a costly rewrite. An alternative is to train a AI-based surrogate model to mimic computationally intensive functions in the code, and run the surrogate instead. PHASM formalizes a development lifecycle for such surrogate models, including discovering functions amenable to replacement with a surrogate model, predicting the resulting performance, identifying the function’s space of inputs and outputs, binding the model to the code, and managing model versions. A suite of software tools for facilitating these steps was written and validated against a set of model problems.
97 MATHEMATICS AND COMPUTING↗
Multifidelity multiobjective optimization for wake-steering strategies
Abstract. Wake steering is an emerging wind power plant control strategy where upstream turbines are intentionally yawed out of perpendicular alignment with the incoming wind, thereby “steering” wakes away from downstream turbines. However, trade-offs between the gains in power production and fatigue loads induced by this control strategy are the subject of continuing investigation. In this study, we present a multifidelity multiobjective optimization approach for exploring the Pareto front of trade-offs between power and loading during wake steering. A large eddy simulation is used as the high-fidelity model, where an actuator line representation is used to model wind turbine blades and a rainflow-counting algorithm is used to compute damage equivalent loads. A coarser simulation with a simpler loads model is employed as a supplementary low-fidelity model. Multifidelity Bayesian optimization is performed to iteratively learn both a surrogate of the low-fidelity model and an additive discrepancy function, which maps the low-fidelity model to the high-fidelity model. Each optimization uses the expected hypervolume improvement acquisition function, weighted by the total cost of a proposed model evaluation in the multifidelity case. The multifidelity approach is able to capture the logit function shape of the Pareto frontier at a computational cost only 30 % that of the single-fidelity approach. Additionally, we provide physical insights into the vortical structures in the wake that contribute to the Pareto front shape.
17 WIND ENERGY↗
Excited states in warm and hot dense matter
Accurate modeling of warm and hot dense matter is challenging in part due to the multitude of excited states that must be considered. Here, we present a variational framework that models these excited states. In this framework an excited state is defined by a set of effective one-electron occupation factors, and the corresponding energy is defined by the effective one-body energy with an exchange and correlation term. Here, the variational framework is applied to an atom-in-plasma model (a generalization of the so-called average atom model). Comparisons with a density functional theory based average atom model generally reveal good agreement in the calculated pressure, but our model also gives access to the excitation energies and charge state distributions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A High-resolution Model of Field-aligned Currents Through Empirical Orthogonal Functions Analysis (MFACE)
Ten years of CHAMP magnetic field measurements are integrated into MFACE, a model of field-aligned currents (FACs) using empirical orthogonal functions (EOFs). EOF1 gives the basic Region-1/Region-2 pattern varying mainly with the interplanetary magnetic field Bz component. EOF2 captures separately the cusp current signature and By-related variability. Compared to existing models, MFACE yields significantly better spatial resolution, reproduces typically observed FAC thickness and intensity, improves on the magnetic local time (MLT) distribution, and gives the seasonal dependence of FAC latitudes and the NBZ current signature. MFACE further reveals systematic dependences on By, including 1) Region-1/Region-2 topology modifications around noon; 2) imbalance between upward and downward maximum current density; 3) MLT location of the Harang discontinuity. Furthermore, our procedure allows quantifying response times of FACs to solar wind driving at the bow shock nose: we obtain 20 minutes and 35-40 minutes lags for the FAC density and latitude, respectively.
CHAMP↗
Predicting fatigue from heart rate signatures using functional logistic regression
Physical fatigue can have adverse effects on humans in extreme environments. Therefore, being able to predict fatigue using easy to measure metrics such as heart rate (HR) signatures has potential to have an impact in real-life scenarios. We apply a functional logistic regression model that uses HR signatures to predict physical fatigue, where physical fatigue is defined in a data-driven manner. Data were collected using commercially available wearable devices on 47 participants hiking the 20.7-mile Grand Canyon rim-to-rim trail in a single day. Fitted model provides good predictions and interpretable parameters for real-life application.
60 APPLIED LIFE SCIENCES↗
Error rate information in attention allocation pilot models
The Northrop urgency decision pilot model was used in a command tracking task to compare the optimized performance of multiaxis attention allocation pilot models whose urgency functions were (1) based on tracking error alone, and (2) based on both tracking error and error rate. A matrix of system dynamics and command inputs was employed, to create both symmetric and asymmetric two axis compensatory tracking tasks. All tasks were single loop on each axis. Analysis showed that a model that allocates control attention through nonlinear urgency functions using only error information could not achieve performance of the full model whose attention shifting algorithm included both error and error rate terms. Subsequent to this analysis, tracking performance predictions for the full model were verified by piloted flight simulation. Complete model and simulation data are presented.
Faulkner, W. H.↗
HFBTHO-AD: Differentiation of a nuclear energy density functional code
The HFBTHO code implements a nuclear energy density functional solver to model the structure of atomic nuclei. HFBTHO has previously been used to calibrate energy functionals and perform sensitivity analysis by using derivative-free methods. To enable derivative-based optimization and uncertainty quantification approaches, we must compute the derivatives of HFBTHO outputs with respect to the parameters of the energy functional, which are a subset of all input parameters of the code. Here, we use the algorithmic/automatic differentiation (AD) tool Tapenade to differentiate HFBTHO. We compare the derivatives obtained using AD against finite-difference approximation and examine the performance of the derivative computation.
Algorithmic differentiation↗
Reynolds-Averaged Turbulence Modeling Using Deep Learning with Local Flow Features: An Empirical Approach
Reynolds-Averaged Navier-Stoke (RANS) models offer an alternative avenue in predicting flow characteristics when the corresponding experiments are difficult to achieve due to geometry complexity, limited budget, or knowledge. RANS models require the knowledge of subgrid scale physics to solve conservation equations for mass, energy, and momentum. Mechanistic turbulence models, such as k-epsilon, are generally evaluated and calibrated for specific flow conditions with various degrees of uncertainty. These models have limited capability to assimilate a substantial amount of data due to model form constraints. Meanwhile, deep learning (DL) has been proven to be universal approximators with the potential to assimilate available, relevant, and adequately evaluated data. Moreover, deep neural networks (DNNs) can create surrogate models without knowing function forms. Such a data-driven approach can be used in updating fluid models based on observations as opposed to hard-wiring models with precalibrated correlations. The paper presents progress in applying DNNs to model Reynolds stress using two machine learning (ML) frameworks. A novel flow feature coverage mapping is proposed to quantify the physics coverage of DL-based closures. It can be used to examine the sufficiency of training data and input flow features for data-driven turbulence models. The case of a backward-facing step is formulated to demonstrate that not only can DNNs discover underlying correlation behind fluid data but also they can be implemented in RANS to predict flow characteristics without numerical stability issues. Finally, the presented research is a crucial stepping-stone toward the data-driven turbulence modeling, which potentially benefits the design of data-driven experiments that can be used to validate fluid models with ML-based fluid closures.
42 ENGINEERING↗
Increased seminal root number associated with domestication improves nitrogen and phosphorus acquisition in maize seedlings
Abstract Background and Aims Domesticated maize (Zea mays ssp. mays) generally forms between two and six seminal roots, while its wild ancestor, Mexican annual teosinte (Zea mays ssp. parviglumis), typically lacks seminal roots. Maize also produces larger seeds than teosinte, and it generally has higher growth rates as a seedling. Maize was originally domesticated in the tropical soils of southern Mexico, but it was later brought to the Mexican highlands before spreading to other parts of the continent, where it experienced different soil resource constraints. The aims of this study were to understand the impacts of increased seminal root number on seedling nitrogen and phosphorus acquisition and to model how differences in maize and teosinte phenotypes might have contributed to increased seminal root number in domesticated maize. Methods Seedling root system architectural models of a teosinte accession and a maize landrace were constructed by parameterizing the functional–structural plant model OpenSimRoot using plants grown in mesocosms. Seedling growth was simulated in a low-phosphorus environment, multiple low-nitrogen environments, and at variable planting densities. Models were also constructed to combine individual components of the maize and teosinte phenotypes. Key Results Seminal roots contributed ~35 % of the nitrogen and phosphorus acquired by maize landrace seedlings in the first 25 d after planting. Increased seminal root number improved plant nitrogen acquisition under low-nitrogen environments with varying precipitation patterns, fertilization rates, soil textures and planting densities. Models suggested that the optimal number of seminal roots for nutrient acquisition in teosinte is constrained by its limited seed carbohydrate reserves. Conclusions Seminal roots can improve the acquisition of both nitrogen and phosphorus in maize seedlings, and the increase in seed size associated with maize domestication may have facilitated increased seminal root number.
59 BASIC BIOLOGICAL SCIENCES↗
Adaptive force-position control for teleoperated manipulators
An adaptive controller with self-tuning can be designed for teleoperated robotic manipulators by determining a time-series model for the function of the teleoperator. Specifically, the position and force exerted by the operator are modelled for determining the derived values for the trajectory of the end-effector of the manipulator. Thus, the adaptive controller can be designed by following the steps which have previously been presented for the controller design of the gross motion.
Koivo, A. J.↗