Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Black hole horizon edge partition functions

We extend a formula for 1-loop black hole determinants by Denef, Hartnoll, and Sachdev (DHS) to spinning fields on any (d + 1)-dimensional static spherically symmetric black hole. By carefully analyzing the regularity condition imposed on the Euclidean eigenfunctions, we reveal an unambiguous bulk-edge split in the 1-loop Euclidean partition function for tensor fields of arbitrary integer spin: the bulk part captures the “renormalized” thermal canonical partition function recently discussed in [1]; the edge part is related to quasinormal modes (QNMs) that fail to analytically continue to a subset of Euclidean modes with enhanced fall-offs near the origin. Since the edge part takes the form of a path integral on S d–1 , this suggests that these are associated with degrees of freedom living on the bifurcation surface in the Lorentzian two-sided black hole geometry. For massive higher spin on static BTZ and massive vector on Nariai black holes, we find that the edge partition function is related to the QNMs with lowest overtone numbers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

What is the gradient of a scalar function defined on a subspace of square matrices?

We illustrate a technique to calculate the gradient of scalar functions that are defined on any arbitrary matrix subspace. It generalizes our earlier work titled “What is the gradient of a scalar function of a symmetric matrix ?”(Indian Journal of Pure and Applied Mathematics (2022), doi:10.1007/s13226-022-00313-x), in which we considered the special case of the subspace of symmetric matrices. Extant methods to calculate the gradient in such cases have an inherent flaw which leads to spurious results that populate several publications, as well as respected textbooks and handbooks on matrix calculus. We examine these sources and results in a rigorous and concrete mathematical setting of a finite-dimensional inner-product space and discover the inherent flaw and also a remedy. We demonstrate two ways to calculate the derivative/gradient and second derivative for scalar functions of matrices defined over an arbitrary matrix subspace; the first method is by considering any (differentiable) extension to the space of square matrices and projection of its gradient onto the given subspace. The second method utilizes an ordered basis and computes each component of the gradient through evaluation of the directional derivative. All the ideas presented are illustrated by non-trivial examples, namely, considering the subspace of 3 × 3 circulant and Toeplitz matrices and presenting the results of gradient-descent with both the spurious and correct gradients. Moreover, our bibliography makes it clear that a rigorous approach to matrix calculus is not common in practice, and our presentation of matrix calculus in the language of inner-product spaces will be significant and meaningful for applied mathematicians, engineers and researchers working in inter-disciplinary fields to avoid the conceptual pitfalls that exist.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

What is the gradient of a scalar function defined on a subspace of square matrices?

We illustrate a technique to calculate the gradient of scalar functions that are defined on any arbitrary matrix subspace. It generalizes our earlier work titled “What is the gradient of a scalar function of a symmetric matrix ?”(Indian Journal of Pure and Applied Mathematics (2022), https://doi.org/10.1007/s13226-022-00313-x), in which we considered the special case of the subspace of symmetric matrices. Extant methods to calculate the gradient in such cases have an inherent flaw which leads to spurious results that populate several publications, as well as respected textbooks and handbooks on matrix calculus. Here, we examine these sources and results in a rigorous and concrete mathematical setting of a finite-dimensional inner-product space and discover the inherent flaw and also a remedy. We demonstrate two ways to calculate the derivative/gradient and second derivative for scalar functions of matrices defined over an arbitrary matrix subspace; the first method is by considering any (differentiable) extension to the space of square matrices and projection of its gradient onto the given subspace. The second method utilizes an ordered basis and computes each component of the gradient through evaluation of the directional derivative. All the ideas presented are illustrated by non-trivial examples, namely, considering the subspace of 3 x 3 circulant and Toeplitz matrices and presenting the results of gradient-descent with both the spurious and correct gradients. Moreover, our bibliography makes it clear that a rigorous approach to matrix calculus is not common in practice, and our presentation of matrix calculus in the language of inner-product spaces will be significant and meaningful for applied mathematicians, engineers and researchers working in inter-disciplinary fields to avoid the conceptual pitfalls that exist.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Solution of an integral equation encountered in studies on radiative transfer in completely absorbing leaf canopies

The leaf normal distribution function is related to the extinction coefficient of direct solar radiation by a Fredholm integral equation of the first kind. The integral equation is solved by the method of constrained least squares. The solutions are stable and correct only in extreme cases and, in general, it is not possible to recover an arbitrary leaf normal distribution function.

Myneni, Ranga B.↗
NASA NTRS1950

Method for determining the frequency-response characteristics of an element or system from the system transient output response to a known input function

A method is presented for the determination of the frequency-response characteristics of an element or system by utilizing the transient output response to a known but arbitrary input to the system. Since the application of special inputs, such as step functions or sinusoids, is often imperfect or impractical, a method for utilizing arbitrary inputs is desirable. Simple flight-test data may be reduced by this method to give the frequency response of an aircraft. Examples are given as determinations of aircraft frequency responses; however, the method can be applied to any type of dynamic system, such as automatic-control components, vibration-absorption equipment, and many types of instruments. The method requires that the arbitrary input function tend to a finite value after a finite time and that the system or element output be measured as a representative quantity having a static sensitivity. (author)

Curfman, Howard J↗
NASA NTRS1986

Arbitrary motion unsteady aerodynamics and its application to rotary-wing aeroelasticity

This paper presents in a unified manner recent research on arbitrary motion unsteady aerodynamics with an emphasis on applications to a number of rotary wing aeroelastic problems. The term arbitrary motion is used to denote growing or decaying oscillations with a certain frequency. The specific topics treated in this paper are: (1) generalization of Greenberg's theory and its application to hingeless rotor aeroelastic stability; (2) description of a new technique for formulating finite state approximations to steady aerodynamic theories, which are suitable for both fixed-wing and rotary-wing applications, and its application to generalize Loewy's theory; (3) comparison of fixed wing and rotary-wing indicial response functions; (4) influence of arbitrary motion aerodynamics, as represented by dynamic inflow on a helicopter in ground resonance. From the discussion of these topics a number of useful conclusions on the fundamental nature of arbitrary motion aerodynamics and its role in rotary-wing aeroelasticity are obtained.

Friedmann, P. P.↗
DOE OSTI2021

The inflationary wavefunction from analyticity and factorization

Here, we study the analytic properties of tree-level wavefunction coefficients in quasi-de Sitter space. We focus on theories which spontaneously break dS boost symmetries and can produce significant non-Gaussianities. The corresponding inflationary correlators are (approximately) scale invariant, but are not invariant under the full conformal group. We derive cutting rules and dispersion formulas for the late-time wavefunction coefficients by using factorization and analyticity properties of the dS bulk-to-bulk propagator. This gives a unitarity method which is valid at tree-level for general n-point functions and for fields of arbitrary mass. Using the cutting rules and dispersion formulas, we are able to compute n-point functions by gluing together lower-point functions. As an application, we study general four-point, scalar exchange diagrams in the EFT of inflation. We show that exchange diagrams constructed from boost-breaking interactions can be written as a finite sum over residues. Finally, we explain how the dS identities used in this work are related by analytic continuation to analogous identities in Anti-de Sitter space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Functor system: a new on-the-fly take on Material Properties based on C++ functions

In the context of solving multiphysics problems, the discretization of the partial differential equations (PDE) at hand often takes the spotlight. However, for most engineering users and even application developers, the discretization of the equations has already been performed. Instead, they are tasked with implementing specific closure relations and material properties. MOOSE has long enabled this using the Materials system. This system relied on the pre-computation of all properties before they are used in the PDE or in postprocessing. In this talk we will introduce the Functor system, which was deployed in MOOSE in 2021, then present a few applications of functors in flow modeling simulations by the NEAMS program. Functors first offer great flexibility in their evaluation. Rather than storing various arrays for material properties, they are evaluated on the fly at the location and state, e.g. current or old value, requested. Unlike regular material properties, several operations such as the time derivative, the divergence and the curl can be requested from a functor. Similar to material properties, functors can be made to depend on arbitrary combinations of variables, functions, postprocessors and other properties. However, unlike material properties, any of these can be substituted for a functor material property. Thanks to this, objects no longer need to be duplicated based on the types of their parameters.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS1987

The LQR/LTR procedure for multivariable feedback control design

This paper provides a tutorial overview of the linear quadratic Gaussian with loop transfer recovery (LQG/LTR) design procedure for linear multivariable feedback systems. LQR/LTR is interpreted as the solution of a specific weighted H-squared tradeoff between transfer functions in the frequency domain. Properties of this solution are examined for both minimum-phase and nonminimum-phase systems. This leads to a formal weight augmentation procedure for the minimum-phase case which permits essentially arbitrary specification of system sensitivity functions in terms of the weights. While such arbitrary specifications are not possible for nonminimum-phase problems, a direct relationship between weights and sensitivities is developed for nonminimum-phase SISO and certain nonminimum-phase MIMO cases which guides the weight selection process.

Stein, Gunter↗
DOE OSTI2025

General field evaluation in high-order meshes on GPUs

Robust and scalable function evaluation at any arbitrary point in the finite/spectral element mesh is required for querying the partial differential equation solution at points of interest, comparison of solution between different meshes, and Lagrangian particle tracking. This is a challenging problem, particularly for high-order unstructured meshes partitioned in parallel with MPI, as it requires identifying the element that overlaps a given point and computing the corresponding reference space coordinates. Here, we present a robust and efficient technique for general field evaluation in large-scale high-order meshes with quadrilaterals and hexahedra. In the proposed method, a combination of globally partitioned and processor-local maps are used to first determine a list of candidate MPI ranks, and then locally candidate elements that could contain a given point. Next, element-wise bounding boxes further reduce the list of candidate elements. Finally, Newton’s method with trust region is used to determine the overlapping element and corresponding reference space coordinates. Since GPU-based architectures have become popular for accelerating computational analyses using meshes with tensor-product elements, specialized kernels have been developed to utilize the proposed methodology on GPUs. The method is also extended to enable general field evaluation on surface meshes. The paper concludes by demonstrating the use of the proposed method in various applications ranging from mesh-to-mesh transfer during r-adaptivity to Lagrangian particle tracking.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Improved surrogates in inertial confinement fusion with manifold and cycle consistencies

Neural networks have become the method of choice in surrogate modeling because of their ability to characterize arbitrary, high-dimensional functions in a data-driven fashion. This paper advocates for the training of surrogates that are 1) consistent with the physical manifold, resulting in physically meaningful predictions, and 2) cyclically consistent with a jointly trained inverse model; i.e., backmapping predictions through the inverse results in the original input parameters. We find that these two consistencies lead to surrogates that are superior in terms of predictive performance, are more resilient to sampling artifacts, and tend to be more data efficient. Using inertial confinement fusion (ICF) as a test-bed problem, we model a one-dimensional semianalytic numerical simulator and demonstrate the effectiveness of our approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Constraints on OPF Surrogates for Learning Stable Local Volt/Var Controllers

We consider the problem of learning local Volt/Var controllers in distribution grids (DGs). Our approach starts from learning separable surrogates that take both local voltages and reactive powers as arguments and predict the reactive power setpoints that approximate optimal power flow (OPF) solutions. We propose an incremental control algorithm and identify two different sets of slope conditions on the local surrogates such that the network is collectively steered toward desired configurations asymptotically. Our results reveal the trade-offs between each set of conditions, with coupled voltage-power slope constraints allowing an arbitrary shape of surrogate functions but risking limitations on exploiting generation capabilities, and reactive power slope constraints taking full advantage of generation capabilities but constraining the shape of surrogate functions. AC power flow simulations on the IEEE 37-bus feeder illustrate their guaranteed stability properties and respective advantages in two DG scenarios.

asymptotic stability↗
NASA NTRS1974

A method of revitalizing Sampson's theory of the Galilean satellites

A method is developed by which Sampson's theory of motion of Jupiter's Galilean satellites can be revitalized by use of algebraic manipulation software on a digital computer. The technique seeks (1) to remove algebraic errors existing in the current Sampson theory, (2) to introduce some neglected effects due to solar interactions and the 3-7 commensurability between the outer two satellites, (3) to allow for nonzero amplitude and phase of the libration, (4) to allow future revision of the arbitrary constants of integration, (5) to express the final results as analytic functions of variations in the numerous arbitrary constants of integration and arbitrary parameters, and (6) to provide analytic partial derivatives by means of which the numerical values of coefficients in the expressions for the coordinates can be adjusted. The level of precision desired is one arc second (Jovicentric) for the coordinates (2, 3.3, 5.2, and 9.1 km, respectively, for satellites I through IV).

Lieske, J. H.↗
NASA NTRS1980

On the equilibrium structures of self-gravitating masses of gas containing axisymmetric magnetic fields

The general equations describing the equilibrium shapes of self-gravitating gas clouds containing axisymmetric magnetic fields are presented. The general equations admit of a large class of solutions. It is shown that if one additional (ad hoc) asumption is made that the mass be spherically symmetrically distributed, then the gas pressure and the boundary conditions are sufficiently constraining that the general topological structure of the solution is effectively determined. The further assumption of isothermal conditions for this case demands that all solutions possess force-free axisymmetric magnetic fields. It is also shown how the construction of aspherical (but axisymmetric) configurations can be achieved in some special cases, and it is demonstrated that the detailed form of the possible equilibrium shapes depends upon the arbitrary choice of the functional form of the variation of the gas pressure along the field lines.

Lerche, I.↗
NASA NTRS1991

Temperature dependence of optical properties of GaAs

The effect of temperature on the optical properties of GaAs was investigated using spectroscopic ellipsometry measurements, between room temperature and about 610 C in increments of 50 C, of pseudodielectric functions and related optical constants of GaAs. A quantitative analysis of the pseudodielectric function spectrum was carried out using a harmonic-oscillator approximation (HOA) to fit the measured dielectric functions. Good fits were obtained with this model, which provides a convenient means of reproducing the GaAs dielectric function at any temperature, by using the temperature-dependent oscillator parameters. The HOA analysis also provides information about band-gap variation with temperature. Using the measured optical constants at a number of fixed temperatures, an algorithm was developed for computing the dielectric function spectrum at an arbitrary temperature in the range 22-610 C.

Yao, Huade↗
NASA NTRS1995

Potential theoretic methods for far field sound radiation calculations

In the area of computational acoustics, procedures which accurately predict the far-field sound radiation are much sought after. A systematic development of such procedures are found in a sequence of papers by Atassi. The method presented here is an alternate approach to predicting far field sound based on simple layer potential theoretic methods. The main advantages of this method are: it requires only a simple free space Green's function, it can accommodate arbitrary shapes of Kirchoff surfaces, and is readily extendable to three-dimensional problems. Moreover, the procedure presented here, though tested for unsteady lifting airfoil problems, can easily be adapted to other areas of interest, such as jet noise radiation problems. Results are presented for lifting airfoil problems and comparisons are made with the results reported by Atassi. Direct comparisons are also made for the flat plate case.

Hariharan, S. I.↗
NASA NTRS1990

Dielectric function of InGaAs in the visible

Measurements are reported of the dielectric function of thermodynamically stable In(x)Ga(1-x)As in the composition range 0.3 equal to or less than X = to or less than 0.7. The optically thick samples of InGaAs were made by molecular beam epitaxy (MBE) in the range 0.4 = to or less than X = to or less than 0.7 and by metal-organic chemical vapor deposition (MOCVD) for X = 0.3. The MBE made samples, usually 1 micron thick, were grown on semi-insulating InP and included a strain release structure. The MOCVD sample was grown on GaAs and was 2 microns thick. The dielectric functions were measured by variable angle spectroscopic ellipsometry in the range 1.55 to 4.4 eV. The data was analyzed assuming an optically thick InGaAs material with an oxide layer on top. The thickness of this layer was estimated by comparing the results for the InP lattice matched material, i.e., X = 0.53, with results published in the literature. The top oxide layer mathematically for X = 0.3 and X = 0.53 was removed to get the dielectric function of the bare InGaAs. In addition, the dielectric function of GaAs in vacuum, after a protective arsenic layer was removed. The dielectric functions for X = 0, 0.3, and 0.53 together with the X = 1 result from the literature to evaluate an algorithm for calculating the dielectric function of InGaAs for an arbitrary value of X(0 = to or less than X = to or less than 1) were used. Results of the dielectric function calculated using the algorithm were compared with experimental data.

Alterovitz, S. A.↗
NASA NTRS1990

Dielectric function of InGaAs in the visible

Measurements are reported of the dielectric function of thermodynamically stable In(x)Ga(1-x)As in the composition range 0.3 equal to or less than X = to or less than 0.7. The optically thick samples of InGaAs were made by molecular beam epitaxy (MBE) in the range 0.4 = to or less than X = to or less than 0.7 and by metal-organic chemical vapor deposition (MOCVD) for X = 0.3. The MBE made samples, usually 1 micron thick, were grown on semi-insulating InP and included a strain release structure. The MOCVD sample was grown on GaAs and was 2 microns thick. The dielectric functions were measured by variable angle spectroscopic ellipsometry in the range 1.55 to 4.4 eV. The data was analyzed assuming an optically thick InGaAs material with an oxide layer on top. The thickness of this layer was estimated by comparing the results for the InP lattice matched material, i.e., X = 0.53, with results published in the literature. The top oxide layer mathematically for X = 0.3 and X = 0.53 was removed to get the dielectric function of the bare InGaAs. In addition, the dielectric function of GaAs in vacuum, after a protective arsenic layer was removed. The dielectric functions for X = 0, 0.3, and 0.53 together with the X = 1 result from the literature to evaluate an algorithm for calculating the dielectric function of InGaAs for an arbitrary value of X (0 = to or less than X = to or less than 1) were used. Results of the dielectric function calculated using the algorithm were compared with experimental data.

Alterovitz, S. A.↗