Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORTHOGONAL FUNCTION”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Spherical balloon response to three-dimensional time-dependent flows

The concept of the Lagrangian displacement of a balloon is introduced. It is shown that the general balloon response problem is extremely complicated because the wind-forcing functions in the balloon equations of motion are functions of the wind velocity vector and its Eulerian first derivatives evaluated at the location of the balloon. The linear perturbation equations for a spherical balloon are derived by perturbing the components of velocity of the balloon about a terminal velocity state which is in equilibrium with a space-time invariant mean horizontal flow. The atmospheric flow is also perturbed such that the resulting equations can be used to analyze the responses of spherical balloons to three-dimensional time-dependent flows. The wind field is represented in terms of a four-fold Fourier integral that involves three orthogonal wave numbers and a frequency, while the balloon components of velocity are represented as Fourier integrals involving a frequency which, in turn, is a function of the wind field wave numbers and frequency and the unperturbed flow components of velocity.

Fichtl, G. H.↗
NASA NTRS1994

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗
DOE OSTI2025

Designing Antifouling and Antimicrobial Interfaces: Structural Characterization using CryoEM, Automated Microscopy, and AI Image Segmentation

The design of functionalized surfaces for interactions with biological systems is critical across sectors such as healthcare, energy, and agriculture. Tailoring materials for specific applications, such as antifouling and antimicrobial surfaces, demands a comprehensive understanding of topology and chemistry across multiple length and time scales on both biological and materials systems. This work presents the development and characterization of nanostructured surfaces with controlled topographies and chemistries that enhance bacterial membrane disruption, reduce biofilm formation, and improve antimicrobial and antifouling capabilities. Two specific use cases will be presented - the use of cellulose nanocrystals (CNCs) for bacterial growth inhibition and the development of antifouling surfaces to prevent protein and bacterial adsorption [1-4]. By leveraging large language models (LLMs) for image segmentation and training [5], we enable automated analysis of terabyte-scale cryogenic electron microscopy (cryoEM) datasets. This analysis provides statistical insights into the biotic/abiotic interface and facilitates automated electron microscopy experiments to mitigate time and dose. The integration of cryogenic electron tomography (cryoET) and cryogenic focused ion beam (cryoFIB) milling enables high-resolution, near-native-state imaging and 3D reconstructions of bio/material interfaces [6]. Orthogonal characterization techniques and computational modeling further enhances our understanding, offering a robust platform for the design and optimization of next-generation functional surfaces [7].

Williams, Alexis [ORNL] (ORCID:0000000252835822)↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Communication: Global minimum search of Ag{sub 10}{sup +} with molecular beam optical spectroscopy

The present study is focused on the optical properties of the Ag{sub 10}{sup +} cluster in the photon energy range ℏω = 1.9–4.4 eV. Absorption spectra are recorded by longitudinal molecular beam depletion spectroscopy and compared to optical response calculations using time-dependent density functional theory. Several cluster isomers obtained by the new pool-based parallel implementation of the Birmingham Cluster Genetic Algorithm, coupled with density functional theory, are used in excited state calculations. The experimental observations, together with additional simulations of ion mobilities for the several geometries found within this work using different models, clearly identify the ground state isomer of Ag{sub 10}{sup +} to be composed of two orthogonal interpenetrating pentagonal bipyramids, having overall D{sub 2d} symmetry.

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

Polarization contributions to intermolecular interactions revisited with fragment electric-field response functions

The polarization energy in intermolecular interactions treated by self-consistent field electronic structure theory is often evaluated using a constraint that the atomic orbital (AO) to molecular orbital transformation is blocked by fragments. This approach is tied to AO basis sets, overestimates polarization energies in the overlapping regime, particularly in large AO basis sets, and lacks a useful complete basis set limit. These problems are addressed by the construction of polarization subspaces based on the responses of isolated fragments to weak electric fields. These subspaces are spanned by fragment electric-field response functions, which can capture effects up to the dipole (D), or quadrupole (DQ) level, or beyond. Schemes are presented for the creation of both non-orthogonal and orthogonal fragment subspaces, and the basis set convergence of the polarization energies computed using these spaces is assessed. Numerical calculations for the water dimer, water–Na{sup +}, water–Mg{sup 2+}, water–F{sup −}, and water–Cl{sup −} show that the non-orthogonal DQ model is very satisfactory, with small differences relative to the orthogonalized model. Additionally, we prove a fundamental difference between the polarization degrees of freedom in the fragment-blocked approaches and in constrained density schemes. Only the former are capable of properly prohibiting charge delocalization during polarization.

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

Modeling and characterizing anisotropic inclusion orientation in heterogeneous material via directional cluster functions and stochastic microstructure reconstruction

We present a framework to model and characterize the microstructure of heterogeneous materials with anisotropic inclusions of secondary phases based on the directional correlation functions of the inclusions. Specifically, we have devised an efficient method to incorporate both directional two-point correlation functions S{sub 2} and directional two-point cluster functions C{sub 2} that contain non-trivial topological connectedness information into the simulated annealing microstructure reconstruction procedure. Our framework is applied to model an anisotropic aluminum alloy and the accuracy of the reconstructed structural models is assessed by quantitative comparison with the actual microstructure obtained via x-ray tomography. We show that incorporation of directional clustering information via C{sub 2} significantly improves the accuracy of the reconstruction. In addition, a set of analytical “basis” correlation functions are introduced to approximate the actual S{sub 2} and C{sub 2} of the material. With the proper choice of basis functions, the anisotropic microstructure can be represented by a handful of parameters including the effective linear sizes of the iron-rich and silicon-rich inclusions along three orthogonal directions. This provides a general and efficient means for heterogeneous material modeling that enables one to significantly reduce the data set required to characterize the anisotropic microstructure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Quantum mechanics without potential function

In the standard formulation of quantum mechanics, one starts by proposing a potential function that models the physical system. The potential is then inserted into the Schrödinger equation, which is solved for the wavefunction, bound states energy spectrum, and/or scattering phase shift. In this work, however, we propose an alternative formulation in which the potential function does not appear. The aim is to obtain a set of analytically realizable systems, which is larger than in the standard formulation and may or may not be associated with any given or previously known potential functions. We start with the wavefunction, which is written as a bounded infinite sum of elements of a complete basis with polynomial coefficients that are orthogonal on an appropriate domain in the energy space. Using the asymptotic properties of these polynomials, we obtain the scattering phase shift, bound states, and resonances. This formulation enables one to handle not only the well-known quantum systems but also previously untreated ones. Illustrative examples are given for two- and three-parameter systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Chebyshev polynomial filtered subspace iteration in the discontinuous Galerkin method for large-scale electronic structure calculations

The Discontinuous Galerkin (DG) electronic structure method employs an adaptive local basis (ALB) set to solve the Kohn-Sham equations of density functional theory in a discontinuous Galerkin framework. The adaptive local basis is generated on-the-fly to capture the local material physics and can systematically attain chemical accuracy with only a few tens of degrees of freedom per atom. A central issue for large-scale calculations, however, is the computation of the electron density (and subsequently, ground state properties) from the discretized Hamiltonian in an efficient and scalable manner. We show in this work how Chebyshev polynomial filtered subspace iteration (CheFSI) can be used to address this issue and push the envelope in large-scale materials simulations in a discontinuous Galerkin framework. We describe how the subspace filtering steps can be performed in an efficient and scalable manner using a two-dimensional parallelization scheme, thanks to the orthogonality of the DG basis set and block-sparse structure of the DG Hamiltonian matrix. The on-the-fly nature of the ALB functions requires additional care in carrying out the subspace iterations. We demonstrate the parallel scalability of the DG-CheFSI approach in calculations of large-scale twodimensional graphene sheets and bulk three-dimensional lithium-ion electrolyte systems. In conclusion, employing 55 296 computational cores, the time per self-consistent field iteration for a sample of the bulk 3D electrolyte containing 8586 atoms is 90 s, and the time for a graphene sheet containing 11 520 atoms is 75 s.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Nonadiabatic Dynamics of Photoexcited cis-Stilbene Using Ab Initio Multiple Spawning

The photochemistry of cis-stilbene proceeds through two pathways: cis-trans isomerization and ring closure to 4a,4b-dihydrophenanthrene (DHP). Despite serving for many decades as a model system for photoisomerization, the photodynamics of cis-stilbene is still not fully understood. In this work, we use ab initio multiple spawning on a SA-2-CASSCF(2,2) potential energy surface to simulate the nonadiabatic dynamics of isolated cis-stilbene. We find the cyclization (to DHP and cis-stilbene) and isomerization (to trans- and cis-stilbene) reaction coordinates to be orthogonal; branching between the two pathways is determined on the S1 excited state within 150 fs of photoexcitation. Furthermore ,trajectory basis functions (TBFs) undergoing cyclization decay rapidly to the ground state in 250 fs, while TBFs moving along the isomerization coordinate remain on the excited state longer, with the majority decaying between 300 and 500 fs. We observe three avoided crossing regions in the dynamics: two along the isomerization coordinate (displaying pyramidalization and migration of an ethylenic hydrogen or phenyl group), and one DHP-like conical intersection along the cyclization coordinate. The isomeric form of the vibrationally hot photoproducts (as determined by measurement 2 ps after photoexcitation) is determined within less than 50 fs of decay to the ground state mediated by passage through a conical intersection. Excess vibrational energy of ground state cis- and trans-stilbene is channelled into phenyl torsions (with mostly opposing directionality). Finally, our simulations are validated by direct comparison to experiment for the absorption spectrum, branching ratio of the three photoproducts (44:52:4 cis-stilbene:trans-stilbene:DHP), and excited state lifetime (520 ± 40 fs).

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

Excited state calculations in solids by auxiliary-field quantum Monte Carlo

Here, we present an approach for ab initio many-body calculations of excited states in solids. Using auxiliary-field quantum Monte Carlo, we introduce an orthogonalization constraint with virtual orbitals to prevent collapse of the stochastic Slater determinants in the imaginary-time propagation. Trial wave functions from density-functional calculations are used for the constraints. Detailed band structures can be calculated. Results for standard semiconductors are in good agreement with experiments; comparisons are also made with GW calculations and the connections and differences are discussed. For the challenging ZnO wurtzite structure, we obtain a fundamental band gap of 3.26(16) eV, consistent with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Origins of anisotropic transport in the electrically switchable antiferromagnet Fe 1/3 NbS 2

Recent experiments on the antiferromagnetic intercalated transition metal dichalcogenide Fe 1/3 NbS 2 have demonstrated reversible resistivity switching by application of orthogonal current pulses below its magnetic ordering temperature, making Fe 1/3 NbS 2 promising for spintronics applications. Here, we perform density functional theory calculations with Hubbard U corrections of the magnetic order, electronic structure, and transport properties of crystalline Fe 1/3 NbS 2 , clarifying the origin of the different resistance states. The two experimentally proposed antiferromagnetic ground states, corresponding to in-plane stripe and zigzag ordering, are computed to be nearly degenerate. In-plane cross sections of the calculated Fermi surfaces are anisotropic for both magnetic orderings, with the degree of anisotropy sensitive to the Hubbard U value. The in-plane resistance, computed within the Kubo linear response formalism using a constant relaxation time approximation, is also anisotropic, supporting a hypothesis that the current-induced resistance changes are due to a repopulating of antiferromagnetic domains. Our calculations indicate that the transport anisotropy of Fe 1/3 NbS 2 in the zigzag phase is reduced relative to stripe, consistent with the relative magnitudes of resistivity changes in experiment. Finally, our calculations reveal the likely directionality of the current-domain response, specifically, which domains are energetically stabilized for a given current direction.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Parabolic sturmians approach to the three-body continuum Coulomb problem

The three-body continuum Coulomb problem is treated in terms of the generalized parabolic coordinates. Approximate solutions are expressed in the form of a Lippmann-Schwinger-type equation, where the Green's function includes the leading term of the kinetic energy and the total potential energy, whereas the potential contains the non-orthogonal part of the kinetic energy operator. As a test of this approach, the integral equation for the (e{sup -}, e{sup -}, He{sup ++}) system has been solved numerically by using the parabolic Sturmian basis representation of the (approximate) potential. Convergence of the expansion coefficients of the solution has been obtained as the basis set used to describe the potential is enlarged.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Raynal–Revai coefficients for a general kinematic rotation

In a three-body system, transitions between different sets of normalized Jacobi coordinates are described as general kinematic transformations that include an orthogonal or a pseudoorthogonal rotation. For such rotations, the Raynal–Revai coefficients execute a unitary transformation between three-body hyperspherical functions. Recurrence relations that make it possible to calculate the Raynal–Revai coefficients for arbitrary angular momenta are derived on the basis of linearized representations of products of hyperspherical functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Large amplitude flexural vibration of thin elastic flat plates and shells

The general equations governing the large amplitude flexural vibration of any thin elastic shell using curvilinear orthogonal coordinates are derived and consist of two coupled, nonlinear, partial differential equations in the normal displacement w and the stress function F. From these equations, the governing equations for the case of shells of revolution or flat plates can be readily obtained as special cases. The material of the shell or plate is isotropic and homogeneous and Hooke's law for the two-dimensional case is valid. It is suggested that the difference between the hardening type of nonlinearity in the case of flat plates and straight beams and the softening type of nonlinearity in the case of shells and rings can, in general, be traced to the amount of curvature present in the underformed median surface of the structure concerned.

Pandalia, K. A. V.↗
NASA NTRS2012

Earthquake Damage Assessment Using Objective Image Segmentation: A Case Study of 2010 Haiti Earthquake

In this study, we perform a case study on imagery from the Haiti earthquake that evaluates a novel object-based approach for characterizing earthquake induced surface effects of liquefaction against a traditional pixel based change technique. Our technique, which combines object-oriented change detection with discriminant/categorical functions, shows the power of distinguishing earthquake-induced surface effects from changes in buildings using the object properties concavity, convexity, orthogonality and rectangularity. Our results suggest that object-based analysis holds promise in automatically extracting earthquake-induced damages from high-resolution aerial/satellite imagery.

image segmentation↗
DOE OSTI2018

Algebraic calculations for spectrum of superintegrable system from exceptional orthogonal polynomials

We introduce an extended Kepler–Coulomb quantum model in spherical coordinates. The Schrödinger equation of this Hamiltonian is solved in these coordinates and it is shown that the wave functions of the system can be expressed in terms of Laguerre, Legendre and exceptional Jacobi polynomials (of hypergeometric type). We construct ladder and shift operators based on the corresponding wave functions and obtain their recurrence formulas. These recurrence relations are used to construct higher-order, algebraically independent integrals of motion to prove superintegrability of the Hamiltonian. The integrals form a higher rank polynomial algebra. By constructing the structure functions of the associated deformed oscillator algebras we derive the degeneracy of energy spectrum of the superintegrable system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗