Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FUNCTIONAL MODELS”

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 379 records · Page 21

Explainable Machine Learning for Functional Data

Black-box machine learning models are recognized as useful tools for prediction applications, but the algorithmic complexity of some models causes interpretation challenges. Explainability methods have been proposed to provide insight into these models, but there is little research focused on supervised modeling with functional data inputs. We argue that, especially in applications of high consequence, it is important to explicitly model the functional dependence in a black-box analysis to not obscure or misrepresent patterns in explanations. As such, we propose the V ariable importance E xplainable E lastic S hape A nalysis (VEESA) pipeline for training supervised machine learning models with functional inputs. The pipeline is an analysis process that includes the data preprocessing, modeling, and post-hoc explanations. The preprocessing is done using elastic functional principal components analysis, which accounts for vertical and horizontal variability in functional data and, ultimately, allows for explanations in the original data space that identify the important functional variability without bias due to correlated variables. Here, we demonstrate the pipeline on two high-consequence applications: explosives classification for national security and inkjet printer identification in forensic science. The applications exhibit the VEESA pipeline’s ability to provide an understanding of the characteristics of the functional data useful for prediction. Code for implementing the pipeline is available in the veesa R package (and supplemental python code).

Elastic Shape Analysis↗

Calculations of multiquark functions in effective models of strong interaction

In this paper we present our results of the investigation of multiquark equations in the Nambu-Jona-Lasinio model with chiral symmetry of SU(2) group in the mean-field expansion. To formulate the mean-field expansion we have used an iteration scheme of solution of the Schwinger-Dyson equations with the fermion bilocal source. We have considered the equations for Green functions of the Nambu-Jona-Lasinio model up to third step for this iteration scheme. To calculate the high-order corrections to the mean-field approximation, we propose the method of the Legendre transformation with respect to the bilocal source, which allows effectively to take into account the symmetry constraints related with the chiral Ward identity. We discuss also the problem of calculating the multiquark functions in the mean-field expansion for Nambu-Jona-Lasinio-type models with other types of the multifermion sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A compendium of human gene functions derived from evolutionary modelling

A comprehensive, computable representation of the functional repertoire of all macromolecules encoded within the human genome is a foundational resource for biology and biomedical research. The Gene Ontology Consortium has been working towards this goal by generating a structured body of information about gene functions, which now includes experimental findings reported in more than 175,000 publications for human genes and genes in experimentally tractable model organisms 1,2 . Here, we describe the results of a large, international effort to integrate all of these findings to create a representation of human gene functions that is as complete and accurate as possible. Specifically, we apply an expert-curated, explicit evolutionary modelling approach to all human protein-coding genes. This approach integrates available experimental information across families of related genes into models that reconstruct the gain and loss of functional characteristics over evolutionary time. The models and the resulting set of 68,667 integrated gene functions cover approximately 82% of human protein-coding genes. The functional repertoire reveals a marked preponderance of molecular regulatory functions, and the models provide insights into the evolutionary origins of human gene functions. We show that our set of descriptions of functions can improve the widely used genomic technique of Gene Ontology enrichment analysis. The experimental evidence for each functional characteristic is recorded, thereby enabling the scientific community to help review and improve the resource, which we have made publicly available.

59 BASIC BIOLOGICAL SCIENCES↗

A photoemission moments model using density functional and transfer matrix methods applied to coating layers on surfaces: Theory

Recent experimental measurements of a bulk material covered with a small number of graphene layers reported by Yamaguchi et al. [NPJ 2D Mater. Appl. 1, 12 (2017)] (on bialkali) and Liu et al.[Appl. Phys. Lett. 110, 041607 (2017)] (on copper) and the needs of emission models in beam optics codes have lead to substantial changes in a Moments model of photoemission. The changes account for (i) a barrier profile and density of states factor based on density functional theory (DFT) evaluations, (ii) a Drude-Lorentz model of the optical constants and laser penetration depth, and (iii) a transmission probability evaluated by an Airy Transfer Matrix Approach. Importantly, the DFT results lead to a surface barrier profile of a shape similar to both resonant barriers and reflectionless wells: the associated quantum mechanical transmission probabilities are shown to be comparable to those recently required to enable the Moments (and Three Step) model to match experimental data but for reasons very different than the assumption by conventional wisdom that a barrier is responsible. The substantial modifications of the Moments model components, motivated by computational materials methods, are developed. The results prepare the Moments model for use in treating heterostructures and discrete energy level systems (e.g., quantum dots) proposed for decoupling the opposing metrics of performance that undermine the performance of advanced light sources like the x-ray Free Electron Laser. The consequences of the modified components on quan-tum yield, emittance, and emission models needed by beam optics codes are discussed. Published by AIP Publishing. https://doi.org/10.1063/1.5008600

36 MATERIALS SCIENCE↗

Revisiting density functionals for the primitive model of electric double layers

Density functional theory (DFT) calculations are typically based on approximate functionals that link the free energy of a multi-body system of interest with the underlying one-body density distributions. Whereas good performance is often proclaimed for new developments, it is difficult to vindicate the theoretical merits relative to alternative versions without extensive comparison with the numerical results from molecular simulations. Besides, approximate functionals may defy statistical-mechanical sum rules and result in thermodynamic inconsistency. Here we compare systematically several versions of density functionals for ionic distributions near a charged surface using the primitive model of electric double layers. We find that the theoretical performance is sensitive not only to the specific forms of the density functional but also to the range of parameter space and the precise properties under consideration. In general, incorporation of the thermodynamic sum rule into the DFT calculations shows significant improvements for both electrochemical properties and ionic distributions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Renormalized site density functional theory for models of ion hydration

Development of accurate statistical mechanics models of molecular liquid systems is a problem of great practical and fundamental importance. Site-density functional theory (SDFT) is one of the promising directions in this area, but its success hinges upon the ability to efficiently reconcile co-existence of two distinct intra- and inter-molecular interaction regimes in a molecular liquid. Renormalized formulation of site-density functional theory (RSDFT), that we have recently developed, resolves this problem by introducing an additional potential field variable that decouples two interaction scales and maps the molecular liquid problem onto the effective simple liquid mixture. This work provides a critical assessment of RSDFT for the hydrated ion system - a problem, which historically has always been one of the most difficult cases for SDFT applications. Using two-site model of water, we perform a comprehensive analysis of hydrated alkali metal and halogen ions, including both structural and free energy based characteristics. The results indicate that RSDFT provides a significant improvement over conventional 3D-RISM implementations, and may prove useful in coarse grained simulations based on two-site solvent models.

Chuev, Gennady N.↗

A Corrected Score Function Framework for Modelling Circadian Gene Expression

Many biological processes display oscillatory behaviour based on an approximately 24 h internal timing system specific to each individual. One process of particular interest is gene expression, for which several circadian transcriptomic studies have identified associations between gene expression during a 24 h period and an individual's health. A challenge with analysing data from these studies is that each individual's internal timing system is offset relative to the 24 h day-night cycle, where day–night cycle time is recorded for each collected sample. Laboratory procedures can accurately determine each individual's offset and determine the internal time of sample collection. However, these laboratory procedures are labour-intensive and expensive. Here, in this paper, we propose a corrected score function framework to obtain a regression model of gene expression given internal time when the offset of each individual is too burdensome to determine. A feature of this framework is that it does not require the probability distribution generating offsets to be symmetric with a mean of zero. Simulation studies validate the use of this corrected score function framework for cosinor regression, which is prevalent in circadian transcriptomic studies. Illustrations with data from three circadian transcriptomic studies further demonstrate that the proposed framework consistently mitigates bias relative to using a score function that does not account for this offset.

59 BASIC BIOLOGICAL SCIENCES↗

The Aemulus Project. II. Emulating the Halo Mass Function

Existing models for the dependence of the halo mass function on cosmological parameters will become a limiting source of systematic uncertainty for cluster cosmology in the near future. We present a halo mass function emulator and demonstrate improved accuracy relative to state-of-the-art analytic models. In this work, mass is defined using an overdensity criteria of 200 relative to the mean background density. Our emulator is constructed from the AEMULUS simulations, a suite of 40 N-body simulations with snapshots from z = 3 to z = 0. These simulations cover the flat wCDM parameter space allowed by recent cosmic microwave background, baryon acoustic oscillation and SNe Ia results, varying the parameters w, Omega(m), Omega(b), sigma(8), N-eff, n(s), and H-0. We validate our emulator using five realizations of seven different cosmologies, for a total of 35 test simulations. These test simulations were not used in constructing the emulator, and were run with fully independent initial conditions. We use our test simulations to characterize the modeling uncertainty of the emulator, and introduce a novel way of marginalizing over the associated systematic uncertainty. We confirm nonuniversality in our halo mass function emulator as a function of both cosmological parameters and redshift. Our emulator achieves better than 1% precision over much of the relevant parameter space, and we demonstrate that the systematic uncertainty in our emulator will remain a negligible source of error for cluster abundance studies through at least the LSST Year 1 data set.

large-scale structure of universe↗

Fully consistent Eulerian Monte Carlo fields method for solving probability density function transport equations in turbulence modeling

The probability density function (PDF) transport equation method is a sophisticated model for the closure of turbulent mixing and turbulent reactive flows. An efficient solution approach for solving the PDF transport equation has been vital for the method to be widely used in applications. The Eulerian Monte Carlo fields (EMCF) method has been developed to solve the PDF transport equation efficiently for decades. A recent work by Wang et al. revealed a serious issue of the EMCF method for not being fully consistent with the PDF transport equation for which the method is designed to solve. This work advances the state of the art by introducing fully consistent EMCF methods for solving the PDF transport equation. The fully consistent EMCF formulations are derived for two different PDF equation forms. The consistency of the EMCF formulations is mathematically confirmed by examining the derived moment transport equations from the EMCF formulations and from the PDF transport equation. The method of manufactured solutions is employed to further verify the consistency and convergence of the different EMCF formulations numerically. The newly introduced EMCF formulations bring the EMCF method to full consistency with the PDF transport equations for the first time.

42 ENGINEERING↗

Advances on statistical/thermodynamical models for unpolarized structure functions

During the eights and nineties many statistical/thermodynamical models were proposed to describe the nucleons' structure functions and distribution of the quarks in the hadrons. Most of these models describe the compound quarks and gluons inside the nucleon as a Fermi / Bose gas respectively, confined in a MIT bag with continuous energy levels. Another models considers discrete spectrum. Some interesting features of the nucleons are obtained by these models, like the sea asymmetries {sup -}d/{sup -}u and {sup -}d-{sup -}u.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Call for Participation: Collaborative Benchmarking of Functional-Structural Root Architecture Models. The Case of Root Water Uptake

Three-dimensional models of root growth, architecture and function are becoming important tools that aid the design of agricultural management schemes and the selection of beneficial root traits. However, while benchmarking is common in many disciplines that use numerical models, such as natural and engineering sciences, functional-structural root architecture models have never been systematically compared. The following reasons might induce disagreement between the simulation results of different models: different representation of root growth, sink term of root water and solute uptake and representation of the rhizosphere. Presently, the extent of discrepancies is unknown, and a framework for quantitatively comparing functional-structural root architecture models is required. We propose, in a first step, to define benchmarking scenarios that test individual components of complex models: root architecture, water flow in soil and water flow in roots. While the latter two will focus mainly on comparing numerical aspects, the root architectural models have to be compared at a conceptual level as they generally differ in process representation. Therefore, defining common inputs that allow recreating reference root systems in all models will be a key challenge. In a second step, benchmarking scenarios for the coupled problems are defined. We expect that the results of step 1 will enable us to better interpret differences found in step 2. This benchmarking will result in a better understanding of the different models and contribute toward improving them. Improved models will allow us to simulate various scenarios with greater confidence and avoid bugs, numerical errors or conceptual misunderstandings. This work will set a standard for future model development.

Plant Sciences↗

Peculiarities of beta functions in sigma models

In this paper we consider perturbation theory in generic two-dimensional sigma models in the so-called first-order formalism, using the coordinate regularization approach. Our goal is to analyze the first-order formalism in application to β functions and compare its results with the standard geometric calculations. Already in the second loop, we observe deviations from the geometric results that cannot be explained by the regularization/renormalization scheme choices. Moreover, in certain cases the first-order calculations produce results that are not symmetric under the classical diffeomorphisms of the target space. Although we could not present the full solution to this remarkable phenomenon, we found some indirect arguments indicating that an anomaly similar to that established in supersymmetric Yang-Mills theory might manifest itself starting from the second loop. We discuss why the difference between two answers might be an infrared effect, similar to that in β functions in supersymmetric Yang-Mills theories.

Sigma Models↗

Water balance model (WBM) v.1.0.0: a scalable gridded global hydrologic model with water-tracking functionality

Abstract. This paper describes the University of New Hampshire Water Balance Model, WBM, a process-based gridded global hydrologic model that simulates the land surface components of the global water cycle and includes water extraction for use in agriculture and domestic sectors. The WBM was first published in 1989; here, we describe the first fully open-source WBM version (v.1.0.0). Earlier descriptions of WBM methods provide the foundation for the most recent model version that is detailed here. We present an overview of the model functionality, utility, and evaluation of simulated global river discharge and irrigation water use. This new version adds a novel suite of water source tracking modules that enable the analysis of flow-path histories on water supply. A key feature of WBM v.1.0.0 is the ability to identify the partitioning of sources for each stock or flux within the model. Three different categories of tracking are available: (1) primary inputs of water to the surface of the terrestrial hydrologic cycle (liquid precipitation, snowmelt, glacier melt, and unsustainable groundwater); (2) water that has been extracted for human use and returned to the terrestrial hydrologic system; and (3) runoff originating from user-defined spatial land units. Such component tracking provides a more fully transparent model in that users can identify the underlying mechanisms generating the simulated behavior. We find that WBM v.1.0.0 simulates global river discharge and irrigation water withdrawals well, even with default parameter settings, and for the first time, we are able to show how the simulation arrives at these fluxes by using the novel tracking functions.

54 ENVIRONMENTAL SCIENCES↗

Large eddy simulation of a reacting spray flame with multiple realizations under compression ignition engine conditions

An n-dodecane spray flame (Spray A from Engine Combustion Network) was simulated using a detailed combustion model along with a dynamic structure LES model to evaluate its performance at engine-relevant conditions and understand the transient behavior of this turbulent flame. The liquid spray was treated with a traditional Lagrangian method and the gas-phase reaction was modeled using a detailed combustion model. A 103-species skeletal mechanism was used for the n-dodecane chemical kinetic model. Significantly different flame structures and ignition processes are observed for the LES compared to those of RANS predictions. The LES data suggests that the first ignition initiates in lean mixture and propagates to rich mixture, and the main ignition happens in rich mixture, preferable less than 0.14 in mixture fraction space. LES was observed to have multiple ignition spots in the mixing layer simultaneously while the main ignition initiates in a clearly asymmetric fashion. The temporal flame development also indicates the flame stabilization mechanism is auto-ignition controlled and modulated by flame propagation. Soot predictions by LES present much better agreement with experiments compared to RANS both qualitatively and quantitatively. Multiple realizations for LES were performed to understand the realization to realization variation and to establish best practices for ensemble-averaging diesel spray flames. The relevance index analysis suggests that an average of 2 and 5 realizations can reach 99\% of similarity to the target average of 16 realizations on the temperature and mixture fraction fields, respectively. However, more realizations are necessary for OH and soot mass fraction due to their high fluctuations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From model to crop: functional characterization of SPL8 in M. truncatula led to genetic improvement of biomass yield and abiotic stress tolerance in alfalfa

Biomass yield, salt tolerance and drought tolerance are important targets for alfalfa ( Medicago sativa L. ) improvement. Medicago truncatula has been developed into a model plant for alfalfa and other legumes. By screening a Tnt1 retrotransposon-tagged M. truncatula mutant population, we identified three mutants with enhanced branching. Branch development determines shoot architecture which affects important plant functions such as light acquisition, resource use and ultimately impacts biomass production. Molecular analyses revealed that the mutations were caused by Tnt1 insertions in the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 8 (SPL8) gene. The M. truncatula spl8 mutants had increased biomass yield, while overexpression of SPL8 in M. truncatula suppressed branching and reduced biomass yield. Scanning electron microscopy (SEM) analysis showed that SPL8 inhibited branching by directly suppressing axillary bud formation. Based on the M. truncatula SPL8 sequence, alfalfa SPL8 (MsSPL8) was cloned and transgenic alfalfa plants were produced. MsSPL8 down-regulated or up-regulated alfalfa plants exhibited similar phenotypes to the M. truncatula mutants or overexpression lines, respectively. Specifically, the MsSPL8 down-regulated alfalfa plants showed up to 43% increase in biomass yield in the first harvest. The impact was even more prominent in the second harvest, with up to 86% increase in biomass production compared to the control. Furthermore, downregulation of MsSPL8 led to enhanced salt and drought tolerance in transgenic alfalfa. Results from this research offer a valuable approach to simultaneously improve biomass production and abiotic stress tolerance in legumes.

09 BIOMASS FUELS↗

Towards a liana plant functional type for vegetation models

Lianas (woody climbers) are crucial components of tropical forests and they have been increasingly recognized to have profound effects on tropical forest carbon dynamics. Despite their importance, lianas' representation in vegetation models remains limited, partly due to the complexity of liana-tree dynamics and the diversity in liana life history strategies. This paper provides a comprehensive review of advances and challenges for mechanistically representing lianas in forest ecosystem models and a proposed path towards effectively representing lianas in these models. Defining a liana plant functional type is a significant challenge because of the high morphological and physiological diversity amongst liana species, and because of their structural association with trees. Here, we identify critical liana traits that likely should contribute to establishing a liana plant functional type, along with key processes to properly represent lianas in ecosystem models. Subsequently, we discuss a variety of possible liana implementation strategies with their associated strengths, limitations, computational costs and data requirements. A fundamental redesign of the tree-centric demographic vegetation models seems appropriate to accommodate the unique growth and competition strategies of lianas. We illustrate the potential of such models with a single-site case study where we disentangle putative mechanisms of liana increasing abundance. Furthermore, we underscore the critical need for comprehensive liana demographic and functional data (including long-term, physiological, and pantropical observations) for the qualitative implementation and evaluation in the proposed modeling efforts. Currently, there is a scarcity of liana data and the data that do exist have a neotropical bias. We finally introduce a new liana functional trait database that can centralize existing liana trait data, incentivize improved data gathering and thus facilitate model development and scientific analyses.

54 ENVIRONMENTAL SCIENCES↗

Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment

We have developed efficient numerical algorithms for the solution of 3D steady-state Poisson-Nernst-Planck equations (PNP) with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are discretized by finite difference scheme and solved iteratively by Gummel method with relaxation. The Nernst-Planck equations are transformed into Laplace equations through the Slotboom transformation. Algebraic multigrid method is then applied to efficiently solve the Poisson equation and the transformed Nernst-Planck equations. A novel strategy for calculating excess chemical potentials through fast Fourier transforms is proposed which reduces computational complexity from O(N2) to O(NlogN) where N is the number of grid points. Integrals involving Dirac delta function are evaluated directly by coordinate transformation which yields more accurate result compared to applying numerical quadrature to an approximated delta function. Numerical results for ion and electron transport in solid electrolyte for Li ion batteries are shown to be in good agreement with the experimental data and the results from previous studies.

Meng, Da↗