Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric invariants”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary schemes. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗
NASA NTRS2013

A Numerical Combination of Extended Boundary Condition Method and Invariant Imbedding Method Applied to Light Scattering by Large Spheroids and Cylinders

The extended boundary condition method (EBCM) and invariant imbedding method (IIM) are two fundamentally different T-matrix methods for the solution of light scattering by nonspherical particles. The standard EBCM is very efficient but encounters a loss of precision when the particle size is large, the maximum size being sensitive to the particle aspect ratio. The IIM can be applied to particles in a relatively large size parameter range but requires extensive computational time due to the number of spherical layers in the particle volume discretization. A numerical combination of the EBCM and the IIM (hereafter, the EBCM+IIM) is proposed to overcome the aforementioned disadvantages of each method. Even though the EBCM can fail to obtain the T-matrix of a considered particle, it is valuable for decreasing the computational domain (i.e., the number of spherical layers) of the IIM by providing the initial T-matrix associated with an iterative procedure in the IIM. The EBCM+IIM is demonstrated to be more efficient than the IIM in obtaining the optical properties of large size parameter particles beyond the convergence limit of the EBCM. The numerical performance of the EBCM+IIM is illustrated through representative calculations in spheroidal and cylindrical particle cases.

invariant imbeddings↗
DOE OSTI2014

Three-dimensional tori and Arnold tongues

This study analyzes an Arnold resonance web, which includes complicated quasi-periodic bifurcations, by conducting a Lyapunov analysis for a coupled delayed logistic map. The map can exhibit a two-dimensional invariant torus (IT), which corresponds to a three-dimensional torus in vector fields. Numerous one-dimensional invariant closed curves (ICCs), which correspond to two-dimensional tori in vector fields, exist in a very complicated but reasonable manner inside an IT-generating region. Periodic solutions emerge at the intersections of two different thin ICC-generating regions, which we call ICC-Arnold tongues, because all three independent-frequency components of the IT become rational at the intersections. Additionally, we observe a significant bifurcation structure where conventional Arnold tongues transit to ICC-Arnold tongues through a Neimark-Sacker bifurcation in the neighborhood of a quasi-periodic Hopf bifurcation (or a quasi-periodic Neimark-Sacker bifurcation) boundary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Frame invariant neural network closures for Kraichnan turbulence

Numerical simulations of geophysical and atmospheric flows have to rely on parameterizations of subgrid scale processes due to their limited spatial resolution. Despite substantial progress in developing parameterization (or closure) models for subgrid scale (SGS) processes using physical insights and mathematical approximations, they remain imperfect and can lead to inaccurate predictions. In recent years, machine learning has been successful in extracting complex patterns from high-resolution spatio-temporal data, leading to improved parameterization models, and ultimately better coarse grid prediction. However, the inability to satisfy known physics and poor generalization hinders the application of these models for real-world problems. In this work, we put forth a frame invariant closure approach to improve the accuracy and generalizability of deep learning-based subgrid scale closure models by embedding physical symmetries directly into the structure of the neural network. Specifically, we utilized specialized layers within the convolutional neural network in such a way that desired constraints are theoretically guaranteed without the need for any regularization terms. We demonstrate our framework for a two-dimensional decaying turbulence test case mostly characterized by the forward enstrophy cascade. We show that our frame invariant SGS model (i) accurately predicts the subgrid scale source term, (ii) respects the physical symmetries such as translation, Galilean, and rotation invariance, and (iii) is numerically stable when implemented in coarse-grid simulation with generalization to different initial conditions and Reynolds number. This work opens up a possibility of connecting physics-based theories and data-driven modeling paradigms, and thus represents a promising step towards the development of physically consistent data-driven turbulence closure models.

42 ENGINEERING↗
DOE OSTI2015

Defects in the AdS/CFT correspondence

We present initial results of the studies of localized line-like sources for a scalar field in anti-de Sitter space. Such configurations reveal some subtleties and require novel numerical analysis. Dual field theory configurations involve simple models of lattices and defects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Communication: Studies of the Lennard-Jones fluid in 2, 3, and 4 dimensions highlight the need for a liquid-state 1/d expansion

The recent theoretical prediction by Maimbourg and Kurchan [e-print http://arxiv.org/abs/1603.05023 (2016)] that for regular pair-potential systems the virial potential-energy correlation coefficient increases towards unity as the dimension d goes to infinity is investigated for the standard 12-6 Lennard-Jones fluid. This is done by computer simulations for d = 2, 3, 4 going from the critical point along the critical isotherm/isochore to higher density/temperature. In both cases the virial potential-energy correlation coefficient increases significantly. For a given density and temperature relative to the critical point, with increasing number of dimension the Lennard-Jones system conforms better to the hidden-scale-invariance property characterized by high virial potential-energy correlations (a property that leads to the existence of isomorphs in the thermodynamic phase diagram, implying that it becomes effectively one-dimensional in regard to structure and dynamics). The present paper also gives the first numerical demonstration of isomorph invariance of structure and dynamics in four dimensions. Our findings emphasize the need for a universally applicable 1/d expansion in liquid-state theory; we conjecture that the systems known to obey hidden scale invariance in three dimensions are those for which the yet-to-be-developed 1/d expansion converges rapidly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A gauge-compatible Hamiltonian splitting algorithm for particle-in-cell simulations using finite element exterior calculus

A particle-in-cell algorithm is derived with a canonical Poisson structure in the formalism of finite element exterior calculus. The resulting method belongs to the class of gauge-compatible splitting algorithms, which exactly preserve gauge symmetries and their associated conservation laws via the momentum map. We numerically demonstrate this time invariance of the momentum map and its usefulness in establishing precise initial conditions with a desired initial electric field and fixed background charge. The restriction of this canonical, finite element Poisson structure to the 1X2P $1\frac {1}{2}$ -dimensional phase space is also considered and simulated numerically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Rotation invariants of vector fields from orthogonal moments

Vector field images are a type of new multidimensional data that appear in many engineering areas. Although the vector fields can be visualized as images, they differ from graylevel and color images in several aspects. In order to analyze them, special methods and algorithms must be originally developed or substantially adapted from the traditional image processing area. Here, we propose a method for the description and matching of vector field patterns under an unknown rotation of the field. Rotation of a vector field is so-called total rotation, where the action is applied not only on the spatial coordinates but also on the field values. Invariants of vector fields with respect to total rotation constructed from orthogonal Gaussian–Hermite moments and Zernike moments are introduced. Their numerical stability is shown to be better than that of the invariants published so far. We demonstrate their usefulness in a real world template matching application of rotated vector fields.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Galactic mass loss - A mild evolutionary correction to the angular size test

Galactic mass loss due to gas production by an evolving stellar population and removal due to supernovae or some other energy source can alter the size of a galaxy over a long time through the adiabatic invariants of the orbits. Numerical and analytic results show that if the mass profile changes homologously, the scale size of the galaxy does also. Under conventional assumptions about galaxy age and initial mass function, this bloating is smaller than the differences between open and flat cosmologies at z = 0.5. If, however, some elliptical galaxies maintain star formation until recent epochs and then lose a (correspondingly) larger fraction of their mass, the effect may be important. It may also be important in decreasing the size of cD galaxies as a result of gas accretion.

Richstone, D. O.↗
DOE OSTI2017

Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method

The Molecular Beam Epitaxial model is derived from the variation of a free energy, that consists of either a fourth order Ginzburg–Landau double well potential or a nonlinear logarithmic potential in terms of the gradient of a height function. One challenge in solving the MBE model numerically is how to develop proper temporal discretization for the nonlinear terms in order to preserve energy stability at the time-discrete level. In this paper, we resolve this issue by developing a first and second order time-stepping scheme based on the “Invariant Energy Quadratization” (IEQ) method. The novelty is that all nonlinear terms are treated semi-explicitly, and the resulted semi-discrete equations form a linear system at each time step. Moreover, the linear operator is symmetric positive definite and thus can be solved efficiently. We then prove that all proposed schemes are unconditionally energy stable. The semi-discrete schemes are further discretized in space using finite difference methods and implemented on GPUs for high-performance computing. Various 2D and 3D numerical examples are presented to demonstrate stability and accuracy of the proposed schemes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Invariant Features for Accurate Predictions of Quantum Chemical UV-vis Spectra of Organic Molecules

Including invariance of global properties of a phys-ical system as an intrinsic feature in graph neural networks (GNNs) enhances the model's robustness and generalizability and reduces the amount of training data required to obtain a desired accuracy for predictions of these properties. Existing open source GNN libraries construct invariant features only for specific GNN architectures. This precludes the generalization of invariant features to arbitrary message passing neural network (MPNN) layers which, in turn, precludes the use of these libraries for new, user-specified predictive tasks. To address this limitation, we implement invariant MPNNs into the flexible and scalable HydraGNN architecture. HydraGNN enables a seamless switch between various MPNNs in a unified layer sequence and allows for a fair comparison between the predictive performance of different MPNNs. We trained this enhanced HydraGNN archi-tecture on the ultraviolet-visible (UV-vis) spectrum of GDB-9 molecules, a feature that describes the molecule's electronic exci-tation modes, computed with time-dependent density functional tight binding (TD-DFTB) and available open source through the GDB-9-Ex dataset. We assess the robustness (i.e., accuracy and generalizability) of the predictions obtained using different invariant MPNNs with respect to different values of the full width at half maximum (FWHM) for the Gaussian smearing of the theoretical peaks. Our numerical results show that incorporating invariance in the HydraGNN architecture significantly enhances both accuracy and generalizability in predicting UV-vis spectra of organic molecules.

Baker, Justin↗
DOE OSTI2015

Onsets of hierarchy truncation and self–consistent Born approximation with quantum mechanics prescriptions invariance

The issue of efficient hierarchy truncation is related to many approximate theories. In this paper, we revisit this issue from both the numerical efficiency and quantum mechanics prescription invariance aspects. The latter requires that the truncation approximation made in Schrödinger picture, such as the quantum master equations and their self–consistent–Born–approximation improvements, should be transferable to their Heisenberg–picture correspondences, without further approximations. We address this issue with the dissipaton equation of motion (DEOM), which is a unique theory for the dynamics of not only reduced systems but also hybrid bath environments. We also highlight the DEOM theory is not only about how its dynamical variables evolve in time, but also the underlying dissipaton algebra. We demonstrate this unique feature of DEOM with model systems and report some intriguing nonlinear Fano interferences characteristics that are experimentally measurable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

On estimating the curvature attributes and strain invariants of deformed surface through radial basis functions

Numerical methods in approximation of derivatives on regular or irregular gridded and scattered data points are of significance in the study of numerous problems in engineering and geoscience. A mathematically accurate and stable method is assessed to compute the differential quantities, while being flexible and easy to implement. Strain-invariant parameters in deformation analysis and curvature attributes in geometrical analysis of approximated surfaces are independent quantities which could be derived numerically through radial basis functions (RBFs) for scattered data points. Approximation of a function or its derivatives through RBFs is highly correlated to its shape parameter which depends on the number and distribution of the particles on the support domain. A procedure to find the optimal shape parameter for RBFs in the influence domain of each nodal point together with the presentation of a manner of calculating the curvature attributes and strain parameters through Gaussian RBFs is proposed here. Results indicate a significant improvement in the approximation accuracy of strain parameters and curvature attributes, while this approximation is more stable when RBFs with optimal shape parameter are implemented rather than the traditional moving least squares.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Quantum modularity and complex Chern-Simons theory

The Quantum Modularity Conjecture of Zagier predicts the existence of a formal power series with arithmetically interesting coefficients that appears in the asymptotics of the Kashaev invariant at each root of unity. Our goal is to construct a power series from a Neumann-Zagier datum (i.e., an ideal triangulation of the knot complement and a geometric solution to the gluing equations) and a complex root of unity $\zeta$. We prove that the coefficients of our series lie in the trace field of the knot, adjoined a complex root of unity. We conjecture that our series are those that appear in the Quantum Modularity Conjecture and confirm that they match the numerical asymptotics of the Kashaev invariant (at various roots of unity) computed by Zagier and the first author. Our construction is motivated by the analysis of singular limits in Chern-Simons theory with gauge group $SL(2,C)$ at fixed level $k$, where $\zeta^k=1$.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Holography beyond conformal invariance and AdS isometry?

We suggest that the principle of holographic duality be extended beyond conformal invariance and AdS isometry. Such an extension is based on a special relation between functional determinants of the operators acting in the bulk and on its boundary, provided that the boundary operator represents the inverse propagators of the theory induced on the boundary by the Dirichlet boundary value problem in the bulk spacetime. This relation holds for operators of a general spin-tensor structure on generic manifolds with boundaries irrespective of their background geometry and conformal invariance, and it apparently underlies numerous O(N{sup 0}) tests of the AdS/CFT correspondence, based on direct calculation of the bulk and boundary partition functions, Casimir energies, and conformal anomalies. The generalized holographic duality is discussed within the concept of the “double-trace” deformation of the boundary theory, which is responsible in the case of large-N CFT coupled to the tower of higher-spin gauge fields for the renormalization group flow between infrared and ultraviolet fixed points. Potential extension of this method beyond the one-loop order is also briefly discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Asymmetric collapse by dissolution or melting in a uniform flow

An advection-diffusion-limited dissolution model of an object being eroded by a two-dimensional potential flow is presented. By taking advantage of the conformal invariance of the model, a numerical method is introduced that tracks the evolution of the object boundary in terms of a time-dependent Laurent series. Simulations of a variety of dissolving objects are shown, which shrink and collapse to a single point in finite time. The simulations reveal a surprising exact relationship, whereby the collapse point is the root of a non-Analytic function given in terms of the flow velocity and the Laurent series coefficients describing the initial shape. This result is subsequently derived using residue calculus. The structure of the non-Analytic function is examined for three different test cases, and a practical approach to determine the collapse point using a generalized Newton-Raphson root-finding algorithm is outlined. These examples also illustrate the possibility that the model breaks down in finite time prior to complete collapse, due to a topological singularity, as the dissolving boundary overlaps itself rather than breaking up into multiple domains (analogous to droplet pinch-off in fluid mechanics). In conclusion, the model raises fundamental mathematical questions about broken symmetries in finite-Time singularities of both continuous and stochastic dynamical systems.

42 ENGINEERING↗
NASA NTRS1995

The Mode 1 mechanical crack tip stress field in hyperelastic and incompressible materials

The finite deformation field of a plane strain Mode 1 crack in a hyperelastic and incompressible material was examined under the assumptions of small scale nonlinearity. Finite element analyses were performed for two different material laws, a Neo-Hookean material and a third order invariant of a Rivlin material. The numerical results for both materials were compared to the appropriate theoretical asymptotic solution. A local cavitation locus surrounding the crack tip was identified for the Neo-Hookean material. For the third order invariant Rivlin material, maximum values of the dominant stress component were found close to the surface of the crack, above and below the deformed crack tip.

Quigley, Claudia J.↗