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 523 records · Page 29

Screened Range-Separated Hybrid Functional with Polarizable Continuum Model Overcomes Challenges in Describing Triplet Excitations in the Condensed Phase Using TDDFT

Long range-corrected (LRC) or range-separated hybrid (RSH) functionals where the long-range (LR) limit of electronic interactions is set to the exact exchange have been shown to correct the tendency of traditional density functional theory (DFT) to underestimate the frontier orbital gap. Consequently, the use of such functionals in calculating electronic excited states using linear response based time-dependent DFT (TDDFT) has been successful in correcting the tendency for underestimating the energies of charge transfer states by DFT-based calculations. More recently formulations of functionals that attenuate the LR limit to address condensed-phase effects to polarize the electronic density have been reported. In particular screened RSH (SRSH) combined with polarizable continuum model (PCM) was benchmarked successfully in reproducing the fundamental gap and charge transfer state energies of molecular systems in the condensed phase. Here we use SRSH-PCM to address triplet excited states, and show its success in obtaining correspondence of the low-lying triplet states to the singlet–triplet gap in a similar way that the fundamental orbital gap corresponds to electron removal and addition energies. Importantly, the accuracy of the SRSH-PCM in calculating triplet excitations stands on the polarization consistent framework in addressing the scalar dielectric constant and without affecting the optimal tuning by triplet energies. Finally, the prospect of even further improving the SRSH-PCM accuracy in calculating triplet states can be achieved by optimal tuning on the basis of the spin multiplicity gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Linear functional minimization for inverse modeling

In this paper, we present a novel inverse modeling strategy to estimate spatially distributed parameters of nonlinear models. The maximum a posteriori (MAP) estimators of these parameters are based on a likelihood functional, which contains spatially discrete measurements of the system parameters and spatiotemporally discrete measurements of the transient system states. The piecewise continuity prior for the parameters is expressed via Total Variation (TV) regularization. The MAP estimator is computed by minimizing a nonquadratic objective equipped with the TV operator. We apply this inversion algorithm to estimate hydraulic conductivity of a synthetic confined aquifer from measurements of conductivity and hydraulic head. The synthetic conductivity field is composed of a low-conductivity heterogeneous intrusion into a high-conductivity heterogeneous medium. Our algorithm accurately reconstructs the location, orientation, and extent of the intrusion from the steady-state data only. Finally, addition of transient measurements of hydraulic head improves the parameter estimation, accurately reconstructing the conductivity field in the vicinity of observation locations.

47 OTHER INSTRUMENTATION↗

Validation of sea ice models using an uncertainty-based distance metric for multiple model variables: NEW METRIC FOR SEA ICE MODEL VALIDATION

Here, we implement a variance-based distance metric (D n ) to objectively assess skill of sea ice models when multiple output variables or uncertainties in both model predictions and observations need to be considered. The metric compares observations and model data pairs on common spatial and temporal grids improving upon highly aggregated metrics (e.g., total sea ice extent or volume) by capturing the spatial character of model skill. The D n metric is a gamma-distributed statistic that is more general than the χ 2 statistic commonly used to assess model fit, which requires the assumption that the model is unbiased and can only incorporate observational error in the analysis. The D n statistic does not assume that the model is unbiased, and allows the incorporation of multiple observational data sets for the same variable and simultaneously for different variables, along with different types of variances that can characterize uncertainties in both observations and the model. This approach represents a step to establish a systematic framework for probabilistic validation of sea ice models. The methodology is also useful for model tuning by using the D n metric as a cost function and incorporating model parametric uncertainty as part of a scheme to optimize model functionality. We apply this approach to evaluate different configurations of the standalone Los Alamos sea ice model (CICE) encompassing the parametric uncertainty in the model, and to find new sets of model configurations that produce better agreement than previous configurations between model and observational estimates of sea ice concentration and thickness.

58 GEOSCIENCES↗

Implementation and evaluation of the unified stomatal optimization approach in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES)

Abstract. Stomata play a central role in regulating the exchange of carbon dioxide and water vapor between ecosystems and the atmosphere. Their function is represented in land surface models (LSMs) by conductance models. The Functionally Assembled Terrestrial Ecosystem Simulator (FATES) is a dynamic vegetation demography model that can simulate both detailed plant demographic and physiological dynamics. To evaluate the effect of stomatal conductance model formulation on forest water and carbon fluxes in FATES, we implemented an optimality-based stomatal conductance model – the Medlyn (MED) model – that simulates the relationship between photosynthesis (A) and stomatal conductance to water vapor (gsw) as an alternative to the FATES default Ball–Woodrow–Berry (BWB) model. To evaluate how the behavior of FATES is affected by stomatal model choice, we conducted a model sensitivity analysis to explore the response of gsw to climate forcing, including atmospheric CO2 concentration, air temperature, radiation, and vapor pressure deficit in the air (VPDa). We found that modeled gsw values varied greatly between the BWB and MED formulations due to the different default stomatal slope parameters (g1). After harmonizing g1 and holding the stomatal intercept parameter (g0) constant for both model formulations, we found that the divergence in modeled gsw was limited to conditions when the VPDa exceeded 1.5 kPa. We then evaluated model simulation results against measurements from a wet evergreen forest in Panama. Results showed that both the MED and BWB model formulations were able to capture the magnitude and diurnal changes of measured gsw and A but underestimated both by about 30 % when the soil was predicted to be very dry. Comparison of modeled soil water content from FATES to a reanalysis product showed that FATES captured soil drying well, but translation of drying soil to modeled physiology reduced the models' ability to match observations. Our study suggests that the parameterization of stomatal conductance models and current model response to drought are the critical areas for improving model simulation of CO2 and water fluxes in tropical forests.

54 ENVIRONMENTAL SCIENCES↗

Radiative strength functions from the energy-localized Brink-Axel hypothesis

Radiative strength functions (RSFs) model the bulk electromagnetic response of highly excited nuclei and are critical inputs for statistical reaction codes. In this paper, we present a definition of the RSF that is consistent with Hauser-Feshbach reaction codes and that can be efficiently computed with the shell model using the Lanczos strength-function (LSF) method. Here, we introduce a variant of the shell-model LSF method that exploits the energy-localized Brink-Axel hypothesis, which makes it possible to compute both electric and magnetic RSFs across all energies relevant to capture reactions. We verify agreement with the conventional definition of RSFs with benchmark calculations of 24 Mg and then present novel results for 56 Fe. For 56 Fe we find that (i) the 𝑀⁢1 RSF shape evolves smoothly with excitation energy, consistent with the energy-localized Brink-Axel hypothesis; (ii) both 𝑀⁢1 and 𝐸⁢1 transitions contribute significantly to the radiative strength below the photoabsorption threshold; and (iii) within the sdpf model space, the strength below 3 MeV observed in Oslo-type experiments cannot be fully reproduced. These results pave the way for a coherent microscopic description of RSFs and further motivate the use of energy-dependent RSFs in modern reaction codes.

Physics - Nuclear physics and radiation physics↗

Covariant Closed String Coherent States

We give the first construction of covariant coherent closed string states, which may be identified with fundamental cosmic strings. We outline the requirements for a string state to describe a cosmic string, and provide an explicit and simple map that relates three different descriptions: classical strings, light cone gauge quantum states, and covariant vertex operators. The resulting coherent state vertex operators have a classical interpretation and are in one-to-one correspondence with arbitrary classical closed string loops.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predictive Models for Regional Hepatic Function Based on 99mTc-IDA SPECT and Local Radiation Dose for Physiologic Adaptive Radiation Therapy

Purpose: High-dose radiation therapy (RT) for intrahepatic cancer is limited by the development of liver injury. This study investigated whether regional hepatic function assessed before and during the course of RT using 99mTc-labeled iminodiacetic acid (IDA) single photon emission computed tomography (SPECT) could predict regional liver function reserve after RT. Methods and Materials: Fourteen patients treated with RT for intrahepatic cancers underwent dynamic 99mTc-IDA SPECT scans before RT, during, and 1 month after completion of RT. Indocyanine green (ICG) tests, a measure of overall liver function, were performed within 1 day of each scan. Three-dimensional volumetric hepatic extraction fraction (HEF) images of the liver were estimated by deconvolution analysis. After coregistration of the CT/SPECT and the treatment planning CT, HEF dose–response functions during and after RT were generated. The volumetric mean of the HEFs in the whole liver was correlated with ICG clearance time. Three models, dose, priori, and adaptive models, were developed using multivariate linear regression to assess whether the regional HEFs measured before and during RT helped predict regional hepatic function after RT. Results: The mean of the volumetric liver HEFs was significantly correlated with ICG clearance half-life time (r=−0.80, P<.0001), for all time points. Linear correlations between local doses and regional HEFs 1 month after RT were significant in 12 patients. In the priori model, regional HEF after RT was predicted by the planned dose and regional HEF assessed before RT (R=0.71, P<.0001). In the adaptive model, regional HEF after RT was predicted by regional HEF reassessed during RT and the remaining planned local dose (R=0.83, P<.0001). Conclusions: 99mTc-IDA SPECT obtained during RT could be used to assess regional hepatic function and helped predict post-RT regional liver function reserve. This could support individualized adaptive radiation treatment strategies to maximize tumor control and minimize the risk of liver damage.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Non-singlet spin structure function in the valon model and low-x-scaling behavior of g{sub 1}{sup NS} and g{sub 1}{sup p}

A next-to-leading order QCD calculation of non-singlet spin structure function, g{sub 1}{sup NS} is presented within the valon representation of hadrons. In the valon model, it is assumed that a nucleon is composed of three dressed valence quarks: the valons. Each valon has its own internal structure, the valence quark with its associated sea quarks and gluons. The results are in good agreement with all available data from SMC, E143, HERMES, and with the newly released data from COMPASS experiments. It appears that the small-x tail of g{sub 1}{sup NS} can be described by a single Regge-type exchange. The relevant parameter of this exchange is given. Finally we show that the polarized proton structure function has a scaling behavior at small x. The relevant parameters of this behavior are given, too.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Two-Stage Estimation and Variance Modeling for Latency-Constrained Variational Quantum Algorithms

The quantum approximate optimization algorithm (QAOA) has enjoyed increasing attention in noisy, intermediate-scale quantum computing with its application to combinatorial optimization problems. QAOA has the potential to demonstrate a quantum advantage for NP-hard combinatorial optimization problems. As a hybrid quantum-classical algorithm, the classical component of QAOA resembles a simulation optimization problem in which the simulation outcomes are attainable only through a quantum computer. The simulation that derives from QAOA exhibits two unique features that can have a substantial impact on the optimization process: (i) the variance of the stochastic objective values typically decreases in proportion to the optimality gap, and (ii) querying samples from a quantum computer introduces an additional latency overhead. In this paper, we introduce a novel stochastic trust-region method derived from a derivative-free, adaptive sampling trust-region optimization method intended to efficiently solve the classical optimization problem in QAOA by explicitly taking into account the two mentioned characteristics. The key idea behind the proposed algorithm involves constructing two separate local models in each iteration: a model of the objective function and a model of the variance of the objective function. Exploiting the variance model allows us to restrict the number of communications with the quantum computer and also helps navigate the nonconvex objective landscapes typical in QAOA optimization problems. In conclusion, we numerically demonstrate the superiority of our proposed algorithm using the SimOpt library and Qiskit when we consider a metric of computational burden that explicitly accounts for communication costs.

Derivative-free Optimization↗

PRIMA-X - Performance Retargeting of Instrumentation, Measurement, and Analysis Technologies for Exascale Computing

The PRIMA-X (Performance Retargeting of Instrumentation, Measurement, and Analysis Technologies for Exascale Computing) project is the successor of the DOE PRIMA (Performance Refactoring of Instrumentation, Measurement, and Analysis Technologies for Petascale Computing) project, which addressed the challenge of creating a core measurement infrastructure that would serve as a common platform for both integrating leading parallel performance systems (notably TAU and Scalasca) and developing next-generation scalable performance tools. The PRIMA-X project shifts the focus away from refactorization of robust performance tools towards a re-targeting of the parallel performance measurement and analysis architecture for extreme scales. The massive concurrency, asynchronous execution dynamics, hardware heterogeneity, and multi-objective prerequisites (performance, power, resilience) that identify exascale systems introduce fundamental constraints on the ability to carry forward existing performance methodologies. In particular, there must be a deemphasis of per-thread observation techniques to significantly reduce the otherwise unsustainable flood of redundant performance data. Instead, it will be necessary to assimilate multi-level resource observations into macroscopic performance views, from which resilient performance metrics can be attributed to the computational features of the application. This requires a scalable framework for node-level and system-wide monitoring and runtime analyses of dynamic performance information. Also, the interest in optimizing parallelism parameters with respect to performance and energy drives the integration of tool capabilities in the exascale environment further. Initially, PRIMA-X was a collaborative project between the University of Oregon (lead institution) and the German Research School for Simulation Sciences (GRS). Because Prof. Wolf, the PI at GRS, accepted a position as full professor at Technische Universitaet Darmstadt (TU Darmstadt) starting February 1st, 2015, the work of GRS was continued at TU Darmstadt. This report covers the work at TU Darmstadt after the transition. The first main accomplishment of TU Darmstadt is the development of different techniques to aggregate performance data on the level of threads. TU Darmstadt evaluated different schemes designed at GRS during the first phase of this project and integrated them into Score-P, a widely used performance measurement framework for HPC applications. The second main accomplishment is a substantial increase of Score-P’s scalability, achieved by improving its internal representation of the underlying system. Third, we extended our performance-modeling tool Extra-P to automatically create models with more than one input parameter (e.g., process count and input size), which was previously impossible. These models allow the performance behavior of an application to be extrapolated to larger machines. Specific applications of the new multi-parameter capability realized in this project include modeling the isoefficiency function of task-based applications and a novel co-design methodology. Further extensions of Extra-P comprise online modeling functionality, automatic refinement of the model search space, and segmented performance modeling. Further contributions include new locality metrics for both shared- and distributed-memory applications, a new way of measuring the impact of inter-application inference on the execution time of an application, and finally an automatic method that selects appropriate OpenMP constructs for the parallelization of sequential code and classifies the affected variables based on their sharing semantics.

97 MATHEMATICS AND COMPUTING↗

Multi-Omics Driven Metabolic Network Reconstruction and Analysis of Lignocellulosic Carbon Utilization in Rhodosporidium toruloides

An oleaginous yeast Rhodosporidium toruloides is a promising host for converting lignocellulosic biomass to bioproducts and biofuels. In this work, we performed multi-omics analysis of lignocellulosic carbon utilization in R. toruloides and reconstructed the genome-scale metabolic network of R. toruloides . High-quality metabolic network models for model organisms and orthologous protein mapping were used to build a draft metabolic network reconstruction. The reconstruction was manually curated to build a metabolic model using functional annotation and multi-omics data including transcriptomics, proteomics, metabolomics, and RB-TDNA sequencing. The multi-omics data and metabolic model were used to investigate R. toruloides metabolism including lipid accumulation and lignocellulosic carbon utilization. The developed metabolic model was validated against high-throughput growth phenotyping and gene fitness data, and further refined to resolve the inconsistencies between prediction and data. We believe that this is the most complete and accurate metabolic network model available for R. toruloides to date.

09 BIOMASS FUELS↗

Local Structure Analysis and Modelling of Lignin‐Based Carbon Composites through the Hierarchical Decomposition of the Radial Distribution Function

Abstract Carbonized lignin has been proposed as a sustainable and domestic source of activated, amorphous, graphitic, and nanostructured carbon for many industrial applications as the structure can be tuned through processing conditions. However, the inherent variability of lignin and its complex physicochemical structure resulting from feedstock and pulping selection make the Process‐Structure‐Property‐Performance (PSPP) relationships hard to define. In this work, radial distribution functions (RDFs) from synchrotron X‐ray and neutron scattering of lignin‐based carbon composites (LBCCs) are investigated using the Hierarchical Decomposition of the Radial Distribution Function (HDRDF) modelling method to characterize the local atomic environment and develop quantitative PSPP relationships. PSPP relationships for LBCCs defined by this work include crystallite size dependence on lignin feedstock as well as increasing crystalline volume fraction, nanoscale composite density, and crystallite size with increasing reduction temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient Simulation of Cascading Outages Using an Energy Function-Embedded Quasi-Steady-State Model

Here, this paper proposed an energy function-embedded quasi-steady-state model for efficient simulation of cascading outages on a power grid while addressing transient stability concerns. Compared to quasi-steady-state models, the proposed model incorporates short-term dynamic simulation and an energy function method to efficiently evaluate the transient stability of a power grid together with outage propagation without transient stability simulation. Cascading outage simulation using the proposed model conducts three steps for each disturbance such as a line outage. First, it performs time-domain simulation for a short term to obtain a post-disturbance trajectory. Second, along the trajectory, the system state with the local maximum potential energy is found and used as the initial point to search for a relevant unstable equilibrium by Newton's method. Third, the transient energy margin is estimated based on this unstable equilibrium to predict an out-of-step condition with generators. The proposed energy function-embedded quasi-steady-state model is tested in terms of its accuracy and time performance on an NPCC 140-bus power system and compared to a quasi-steady-state model embedding transient stability simulation.

Guo, Zhenping [Univ. of Tennessee, Knoxville, TN (↗

β-function formalism for inflationary models with a non minimal coupling with gravity

We discuss the introduction of a non minimal coupling between the inflaton and gravity in terms of our recently proposed β-function formalism for inflation. Via a field redefinition we reduce to the case of minimally coupled theories. The universal attractor at strong coupling has a simple explanation in terms of the new field. Generalizations are discussed and the possibility of evading the universal attractor is shown.

79 ASTRONOMY AND ASTROPHYSICS↗

Mechanistic Study of Functional Electrolyte Solvents for High-Voltage Lithium Batteries

The pervasive use of Ni-rich cathode active materials, e.g., LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811), for high-energy-density Li-ion batteries (LIBs) has been hindered by rapid battery capacity decay when cycled with high charge cutoff voltages due to electrolyte decomposition in the conventional carbonate solvent-based electrolytes, oxidative parasitic side reactions at the electrolyte/cathode interface, and irreversible phase changes in the cathode active materials leading to dissolution of transition metals into the electrolytes. Various functional electrolyte solvents have been studied to tackle the above technical challenges, yet the roles of individual solvents in the performance of LIBs remain poorly understood. Here, in this study, we systematically investigate electrochemical performance mechanisms of fluorinated and organosilicon single solvents and cosolvents, for the first time, in high-voltage Li/NMC811 batteries, using electrochemical and analytical characterizations and density functional theory modeling. We observe that some unique combinations of the functional solvents can lead to exceptionally stable high-voltage cycle performance in the Ni-rich cathode-based LIBs. Our mechanistic study reveals that the synergistic effect of solvents plays a vital role in enabling electrochemical stability at both the Ni-rich cathode and the Li metal anode. Understanding the electrochemical performance mechanisms of functional solvents can greatly help in designing and formulating advanced electrolytes that enable the development of high-voltage, high-energy-density, long-cycle-life lithium batteries.

density functional theory modeling↗