Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional generalized linear model”

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 217 records · Page 12

LASSO for CALPHAD Model Selection Enables Data-Efficient Thermodynamic Modeling: An Application in Thermochemical Hydrogen Production Materials

Phenomenological CALPHAD (CALculation of PHAse Diagrams) models, widely used for multicomponent materials, often contain a considerable number of parameters and require fitting using data from a relatively small number of experimental measurements or theoretical calculations. Sometimes these parameters are introduced for the purpose of improving model fits but without clear physical justification, which leads to overparametrized models with poor generalization performance. Automated approaches for optimal model selection based on the available data therefore become critical. Here, in this work, a least absolute shrinkage and selection operator (LASSO)-based approach is developed for model selection by leveraging the linearity of the CALPHAD model with respect to its parameters to convert the model selection and fitting to a LASSO minimization problem. We demonstrate its utility for thermodynamic modeling of thermochemical hydrogen (TCH) production materials using lanthanum strontium manganite (LSM) as an example. Various TCH-relevant properties, including oxygen stoichiometry as a function of oxygen partial pressure, enthalpy of reduction, and entropy of reduction, are successfully predicted with reasonable accuracy using a minimal set of model parameters. Importantly, the model selection and fitting involve minimal human decision; it can therefore be applied to high-throughput DFT defect calculations and yield efficient workflows for TCH material modeling and optimization.

CALPHAD↗

A data-driven peridynamic continuum model for upscaling molecular dynamics

Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials’ behavior, we propose a learning framework to extract, from MD data, an optimal Linear Peridynamic Solid (LPS) model as a surrogate for MD displacements. To maximize the accuracy of the learnt model we allow the peridynamic influence function to be partially negative, while preserving the well-posedness of the resulting model. To achieve this, we provide sufficient well-posedness conditions for discretized LPS models with sign-changing influence functions and develop a constrained optimization algorithm that minimizes the equation residual while enforcing such solvability conditions. This framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our two-dimensional tests show the robustness of the proposed algorithm on validation data sets that include thermal noise, different domain shapes and external loadings, and discretizations substantially different from the ones used for training.

homogenization↗

Static shape determination and control for large space structures. I - The flexible beam. II - A large space antenna

A method for determining and controlling the shape of large, continuous space structures by means of discrete or pointwise observations and control devices is presented. The general linear boundary value problem satisfied by a one-dimensional shape function is defined, and the existence of solutions is studied. The static shape control problems for one-dimensional systems with and without rigid body modes and the static shape estimation problem are presented and solved. Eigenfunction expansions are presented which provide approximations to the algorithm terms when the associated Green's function is not known. An integral operator approach is applied to the multidimensional static problem, and the results are illustrated with a finite element model of the disk of a large space antenna which assumes no rigid body modes. It is shown that the shape control algorithm must be modified for systems with rigid body modes.

Weeks, C. J.↗

Accretion of matter and spectra of binary X-ray sources in massive gravity

We study low-mass binary X-ray sources involving stellar mass black holes within massive gravity. Regarding the accretion disk, we adopt the standard model for an optically thick, cool, and geometrically thin disk by Shakura–Sunyaev. For the gravitational field generated by the black hole, we consider the analogue of the Schwarzschild–de Sitter space–time of Einstein’s theory in massive gravity, for which we found an additional term linear in the radial coordinate. Then, we compute the radial velocity, the energy density and the pressure as a function of the radial coordinate, and the X-ray emission’s soft spectral component expected from the disk. We also investigated in detail the impact of this new geometry. Our result indicates that by using observed spectra from confirmed X-ray binaries involving astrophysical black holes, we can put strong constraints on alternative theories of gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Low-Energy Ion Scattering Intensities from Supported Nanoparticles: The Spherical Cap Model

Supported nanoparticles are of great importance to many technologies like fuel processing and chemical synthesis using catalysts and electrocatalysts, energy storage and generation using fuel cells and batteries, electrochemistry, magnetic information storage, and more. Low-energy ion scattering spectroscopy (LEIS) with noble gas ions like He+ is a powerful tool for the characterization of nanoparticles dispersed across flat support surfaces due to its ability to probe the elemental composition in the topmost atomic layer of a surface, providing quantitative information regarding the size and number density of nanoparticles. In this work, we present a derivation of the LEIS intensities expected from nanoparticles and the support material as a function of the average particle size, their number per unit area, and their contact angle with the support when modeled as spherical caps of the nanoparticle material dispersed over the surface of a flat support. The model assumes that the ion intensities are determined only by the physical blocking of linear ion trajectories and independent of the tilt angle of the local surface relative to the incident and scattered ion directions, an assumption we support by quantitative modeling of published data which tested tilt-angle effects. The model is a generalization to arbitrary contact angles of the hemispherical cap model, which assumes 90° contact angle and has been widely used to model spectroscopic signals in LEIS (and also in Auger and photoelectron spectroscopies) during nanoparticle growth. Here, this new model quantitatively reveals how LEIS signals are sensitive not only to the diameter and number density of the nanoparticle but also to their contact angle (or height/diameter ratio). With the use of additional data (e.g., from microscopy or adsorption microcalorimetry), the model presented here will enable more accurate determination of the average size, shape, and number density of supported nanoparticles based on LEIS intensity measurements.

99 GENERAL AND MISCELLANEOUS↗

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗

A model of the human observer in failure detection tasks

A model for the human observer in failure detection tasks is proposed which consists of two stages: a linear estimator and a decision mechanism. The estimator is a Kalman filter, and the decision mechanism, which is based on Wald's sequential analysis, leads to a decision function which is the integration of the filter residuals. The final result is a simple detection system which depends on only three parameters, and the sensitivity of the model to these parameters is analyzed. The results of an experiment designed to test the validity of the model are reported. The question of open and closed decision intervals as well as the generalization of the model to more complicated cases is discussed.

Gai, E. G.↗

Jet Surface Interaction-Scrubbing Noise

Generation of sound due to scrubbing of a jet flow past a nearby solid surface is investigated within the framework of the generalized acoustic analogy theory. The analysis applies to the boundary layer noise generated at and near a wall, and excludes the scattered noise component that is produced at the leading or the trailing edge. While compressibility effects are relatively unimportant at very low Mach numbers, frictional heat generation and thermal gradient normal to the surface could play important roles in generation and propagation of sound in high speed jets of practical interest. A general expression is given for the spectral density of the far-field sound as governed by the variable density Pridmore- Brown equation. The propagation Green's function should be solved numerically starting with the boundary conditions on the surface and subject to specified mean velocity and temperature profiles between the surface and the observer. The equivalent sources of aerodynamic sound are associated with non-linear momentum flux and enthalpy flux terms that appear in the linearized Navier-Stokes equations. These multi-pole sources should be modeled and evaluated with input from a Reynolds-Averaged Navier-Stokes (RANS) solver with an appropriate turbulence model.

Khavaran, Abbas↗

A general unified approach to modelling switching-converter power stages

A method for modeling switching-converter power stages is developed, whose starting point is the unified state-space representation of the switched networks and whose end result is either a complete state-space description or its equivalent small-signal low-frequency linear circuit model. A new canonical circuit model is proposed, whose fixed topology contains all the essential input-output and control properties of any dc-to-dc switching converter, regardless of its detailed configuration, and by which different converters can be characterized in the form of a table conveniently stored in a computer data bank to provide a useful tool for computer aided design and optimization. The new canonical circuit model predicts that, in general, switching action introduces both zeros and poles into the duty ratio to output transfer function in addition to those from the effective filter network.

Middlebrook, R. D.↗

Dynamic Reconstruction and Multivariable Control for Force-Actuated, Thin Facesheet Adaptive Optics

The Multiple Mirror Telescope (MMT) under development at the University of Arizona takes a new approach in adaptive optics placing a large (0.65 m) force-actuated, thin facesheet deformable mirror at the secondary of an astronomical telescope, thus reducing the effects of emissivity which are important in IR astronomy. However, The large size of the mirror and low stiffness actuators used drive the natural frequencies of the mirror down into the bandwidth of the atmospheric distortion. Conventional adaptive optics takes a quasi-static approach to controlling the, deformable mirror. However, flexibility within the control bandwidth calls for a new approach to adaptive optics. Dynamic influence functions are used to characterize the influence of each actuator on the surface of the deformable mirror. A linearized model of atmospheric distortion is combined with dynamic influence functions to produce a dynamic reconstructor. This dynamic reconstructor is recognized as an optimal control problem. Solving the optimal control problem for a system with hundreds of actuators and sensors is formidable. Exploiting the circularly symmetric geometry of the mirror, and a suitable model of atmospheric distortion, the control problem is divided into a number of smaller decoupled control problems using circulant matrix theory. A hierarchic control scheme which seeks to emulate the quasi-static control approach that is generally used in adaptive optics is compared to the proposed dynamic reconstruction technique. Although dynamic reconstruction requires somewhat more computational power to implement, it achieves better performance with less power usage, and is less sensitive than the hierarchic technique.

Grocott, Simon C. O.↗

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗

Density functional theory of material design: fundamentals and applications—II

Abstract This is the second and the final part of the review on density functional theory (DFT), referred to as DFT-II. In the first review, DFT-I, we have discussed wavefunction-based methods, their complexity, and basics of density functional theory. In DFT-II, we focus on fundamentals of DFT and their implications for the betterment of the theory. We start our presentation with the exact DFT results followed by the concept of exchange-correlation (xc) or Fermi-Coulomb hole and its relationship with xc energy functional. We also provide the exact conditions for the xc-hole, xc-energy and xc-potential along with their physical interpretation. Next, we describe the extension of DFT for non-integer number of electrons, the piecewise linearity of total energy and discontinuity of chemical potential at integer particle numbers, and derivative discontinuity of the xc potential, which has consequences on fundamental gap of solids. After that, we present how one obtains more accurate xc energy functionals by going beyond the LDA. We discuss the gradient expansion approximation (GEA), generalized gradient approximation (GGA), and hybrid functional approaches to designing better xc energy functionals that give accurate total energies. However, these functionals fail to predict properties like the ionization potential and the band gap. Thus, we next describe different methods of modelling these potentials and results of their application for calculation of the band gaps of different solids to highlight accuracy of different xc potentials. Finally, we conclude with a glimpse on orbital-free density functional theory and the machine learning approach.

36 MATERIALS SCIENCE↗

A Hydrodynamic Theory for Spatially Inhomogeneous Semiconductor Lasers: Numerical Results - 2

We present numerical results of the diffusion coefficients (DCs) in the coupled diffusion model derived in the preceding paper for a semiconductor quantum well. These include self and mutual DCs in the general two-component case, as well as density- and temperature-related DCs under the single-component approximation. The results are analyzed from the viewpoint of free Fermi gas theory with many-body effects incorporated. We discuss in detail the dependence of these DCs on densities and temperatures in order to identify different roles played by the free carrier contributions including carrier statistics and carrier-LO phonon scattering, and many-body corrections including bandgap renormalization and electron-hole (e-h) scattering. In the general two-component case, it is found that the self- and mutual- diffusion coefficients are determined mainly by the free carrier contributions, but with significant many-body corrections near the critical density. Carrier-LO phonon scattering is dominant at low density, but e-h scattering becomes important in determining their density dependence above the critical electron density. In the single-component case, it is found that many-body effects suppress the density coefficients but enhance the temperature coefficients. The modification is of the order of 10% and reaches a maximum of over 20% for the density coefficients. Overall, temperature elevation enhances the diffusive capability or DCs of carriers linearly, and such an enhancement grows with density. Finally, the complete dataset of various DCs as functions of carrier densities and temperatures provides necessary ingredients for future applications of the model to various spatially inhomogeneous optoelectronic devices.

Li, Jianzhong↗

Unsteady aerodynamic modeling and active aeroelastic control

Unsteady aerodynamic modeling techniques are developed and applied to the study of active control of elastic vehicles. The problem of active control of a supercritical flutter mode poses a definite design goal stability, and is treated in detail. The transfer functions relating the arbitrary airfoil motions to the airloads are derived from the Laplace transforms of the linearized airload expressions for incompressible two dimensional flow. The transfer function relating the motions to the circulatory part of these loads is recognized as the Theodorsen function extended to complex values of reduced frequency, and is termed the generalized Theodorsen function. Inversion of the Laplace transforms yields exact transient airloads and airfoil motions. Exact root loci of aeroelastic modes are calculated, providing quantitative information regarding subcritical and supercritical flutter conditions.

Edwards, J. W.↗

A Self-Stabilizing Synchronization Protocol for Arbitrary Digraphs

This paper presents a self-stabilizing distributed clock synchronization protocol in the absence of faults in the system. It is focused on the distributed clock synchronization of an arbitrary, non-partitioned digraph ranging from fully connected to 1-connected networks of nodes while allowing for differences in the network elements. This protocol does not rely on assumptions about the initial state of the system, other than the presence of at least one node, and no central clock or a centrally generated signal, pulse, or message is used. Nodes are anonymous, i.e., they do not have unique identities. There is no theoretical limit on the maximum number of participating nodes. The only constraint on the behavior of the node is that the interactions with other nodes are restricted to defined links and interfaces. This protocol deterministically converges within a time bound that is a linear function of the self-stabilization period. We present an outline of a deductive proof of the correctness of the protocol. A bounded model of the protocol was mechanically verified for a variety of topologies. Results of the mechanical proof of the correctness of the protocol are provided. The model checking results have verified the correctness of the protocol as they apply to the networks with unidirectional and bidirectional links. In addition, the results confirm the claims of determinism and linear convergence. As a result, we conjecture that the protocol solves the general case of this problem. We also present several variations of the protocol and discuss that this synchronization protocol is indeed an emergent system.

Malekpour, Mahyar R.↗

Tracing the Mott-Hubbard transition in one-dimensional Hubbard models without Umklapp scattering

Here we apply the density-matrix renormalization group (DMRG) method to a one-dimensional Hubbard model that lacks Umklapp scattering and thus provides an ideal case to study the Mott-Hubbard transition analytically and numerically. The model has a linear dispersion and displays a metal-to-insulator transition when the Hubbard interaction U equals the bandwidth U c = W, where the single-particle gap opens linearly, Δ(U ≥ W) = U - W. The simple nature of the elementary excitations permits us to determine numerically with high accuracy the critical interaction strength and the gap function in the thermodynamic limit. The jump discontinuity of the momentum distribution $n_k$ at the Fermi wave number $k_F$ = 0 cannot be used to locate accurately U c from finite-size systems. However, the slope of $n_k$ at the band edges $k_B$ = ±π, reveals the formation of a single-particle bound state which can be used to determine U c reliably from $n_k$ using accurate finite-size data.

1-dimensional systems↗

Data-driven reduced-order models for port-Hamiltonian systems with operator inference

Hamiltonian operator inference has been developed in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. The method constructs a low-dimensional model using only data and knowledge of the functional form of the Hamiltonian. The resulting ROMs preserve the intrinsic structure of the system, ensuring that the mechanical and physical properties of the system are maintained. In this work, we extend this approach to port-Hamiltonian systems, which generalize Hamiltonian systems by including energy dissipation, external input, and output. Based on snapshots of the system’s state and output, together with the information about the functional form of the Hamiltonian, reduced operators are inferred through optimization and are then used to construct data-driven ROMs. To further alleviate the complexity of evaluating nonlinear terms in the ROMs, a hyper-reduction method via discrete empirical interpolation is applied. Accordingly, we derive error estimates for the ROM approximations of the state and output. Lastly, we demonstrate the structure preservation, as well as the accuracy of the proposed port-Hamiltonian operator inference framework, through numerical experiments on a linear mass–spring-damper problem and a nonlinear Toda lattice problem.

97 MATHEMATICS AND COMPUTING↗

Highly accurate and constrained density functional obtained with differentiable programming

Using an end-to-end differentiable implementation of the Kohn-Sham self-consistent field equations, we obtain a highly accurate neural network–based exchange and correlation (XC) functional of the electronic density. The functional is optimized using information on both energy and density while exact constraints are enforced through an appropriate neural network architecture. Here we evaluate our model against different families of XC approximations and show that at the meta-GGA level our functional exhibits unprecedented accuracy for both energy and density predictions. For nonempirical functionals, there is a strong linear correlation between energy and density errors. We use this correlation to define an XC functional quality metric that includes both energy and density errors, leading to an improved way to rank different approximations.

36 MATERIALS SCIENCE↗