Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

GPU acceleration of hybrid functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing hybrid functional calculations in generalized Kohn–Sham density functional theory. In particular, we develop a batch variant of the recently formulated Kronecker product-based linear solver for the simultaneous solution of multiple linear systems. We then develop a modular, math kernel based implementation for hybrid functionals on NVIDIA architectures, where computationally intensive operations are offloaded to the GPUs, while the remaining workload is handled by the central processing units (CPUs). Considering bulk and slab examples, we demonstrate that GPUs enable up to 8× speedup in node-hours and 80× in core-hours compared to CPU-only execution, reducing the time to solution on V100 GPUs to around 300 s for a metallic system with over 6000 electrons, and significantly reducing the computational resources required for a given wall time.

Kohn-Sham density functional theory↗
NASA NTRS1977

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.↗
NASA NTRS1989

Approximate formulation of redistribution in the Ly-alpha, Ly-beta, H-alpha system

Simple approximate formulas are given for the coupled redistribution of Ly-alpha, Ly-beta, and H-alpha, by using well-defined approximations to an essentially exact formulation. These formulas incorporate all the essential physics including Raman scattering, lower state radiative decay, and correlated terms representing emission during a collision which must be retained in order that the emission coefficients are properly behaved in the line wings. Approximate expressions for the appropriate line broadening parameters are collected. Finally, practical expressions for the source functions are given. These are formulated through newly introduced nonimpact redistribution functions, which are shown to be reasonably approximated by existing (ordinary and generalized) redistribution functions.

Cooper, J.↗
NASA NTRS1990

Spin bipolaron in the framework of emery model for high-T(sub c) copper oxide superconductors

The high-T(sub c) oxide compounds discovered recently exhibit a number of interesting physical properties. Two-dimensional antiferromagnetic spin order has been observed in these materials at the oxygen deficiency. This fact can be explained by strong correlation of the spins, situated on Cu sites in the conducting planes of the oxide superconductors. The doping or the oxygen deficiency lead to the occurrence of holes, occupying the oxygen p-orbitals according to the Emery model. At the small hole concentration they can move along the antiferromagnetic lattice of spins, localized on Cu sites. Researchers consider the two holes situation and describe in what way their behavior depends on the antiferromagnetic exchange interation J. It is known that in the framework of Hubbard model with strong on-site Coulomb repulsion, a single hole can form a spin polaron of the large radius. It is reasonable to admit that two holes with parallel spins (triplet) form the spin bipolaron complex owing to the hole excitations' capability to polarize Cu spin surroundings. Such an excitation was considered in the phenomenological way. Here the problem is discussed on the basis of the microscopic approach in the framework of the variational principle. A special kind of wave function is used for such a purpose. The wave function is constructed by generalizing the trial functions proposed in over two holes excitation situation (triplet) and then the region of spin bipolaron existance in the framework of Emery model is studied. In this model the Hamiltonian can be easily rewritten by forming the oxygen states transforming as the irreducible representations of the group D(sub 4).

Golub, A. A.↗
DOE OSTI2022

General gauge-Yukawa-quartic β-functions at 4-3-2-loop order

We determine the full set of coefficients for the completely general 4-loop gauge and 3-loop Yukawa β-functions for the most general renormalizable four-dimensional theories. Using a complete parametrization of the β-functions, we compare the general form to the specific β-functions of known theories to constrain the unknown coefficients. The Weyl consistency conditions provide additional constraints, completing the determination.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1972

Stability regions of large-scale systems.

Using an aggregated comparison system, methods are developed for estimating regions of asymptotic stability for large-scale systems composed of interconnected, exponentially stable subsystems; in general some of these subsystems are not globally stable. Three different forms of Liapunov functions are shown to be suitable for describing the stability behavior of the aggregated system, but one form in particular (heretofore not used in this context) is found to be the most natural. Using this function, some general results are obtained.

Weissenberger, S.↗
NASA NTRS1973

Stability regions of large-scale systems

Using an aggregated comparison system, methods are developed for estimating regions of asymptotic stability for large-scale systems composed of interconnected, exponentially stable subsystems; in general, some of these subsystems are not globally stable. Three different forms of Liapunov functions are shown to be suitable for describing the stability behavior of the aggregated system, but one form in particular, heretofore not used in this context, is found to be the most natural. Using this function, some general results are obtained.

Weissenberger, S.↗
NASA NTRS1975

Structural Properties and Estimation of Delay Systems

Two areas in the theory of delay systems were studied: structural properties and their applications to feedback control, and optimal linear and nonlinear estimation. The concepts of controllability, stabilizability, observability, and detectability were investigated. The property of pointwise degeneracy of linear time-invariant delay systems is considered. Necessary and sufficient conditions for three dimensional linear systems to be made pointwise degenerate by delay feedback were obtained, while sufficient conditions for this to be possible are given for higher dimensional linear systems. These results were applied to obtain solvability conditions for the minimum time output zeroing control problem by delay feedback. A representation theorem is given for conditional moment functionals of general nonlinear stochastic delay systems, and stochastic differential equations are derived for conditional moment functionals satisfying certain smoothness properties.

Kwong, R. H. S.↗
NASA NTRS1992

Loop transfer recovery for general nonminimum phase discrete time systems. II - Design

The authors consider the design of controllers for the recovery of target loop transfer function or sensitivity and complementary sensitivity functions for general nonminimum phase discrete time systems. The necessary design constraints and the available design freedom are reviewed. In view of the available freedom, possible specifications on the eigenstructure of the observer dynamic matrix are formulated. Three different types of controllers which are respectively based on prediction, current, and reduced-order estimators are considered. For each one of those controllers, three different design techniques are developed. The first one is an eigenstructure assignment scheme, while the other two are optimization-based designs. The eigenstructure assignment method yields a controller design which achieves any chosen recovery error matrix among a set of admissible recovery error matrices. On the other hand, one of the optimization-based design methods leads to a controller that achieves a recovery error matrix having the infimum H-infinity norm, while the other does the same except it achieves a recovery error matrix having the infimum H2 norm.

Chen, Ben M.↗
NASA NTRS1988

A generalization of Petterssen's frontogenesis function and its relation to the forcing of vertical motion

The classic Peterssen frontogenesis function, defined as the Lagrangian rate of change of the magnitude of the potential temperature gradient, is generalized to apply to the vector potential temperature gradient. The derivation of vector frontogenesis in natural coordinates is shown, and analytical examples of frontal evolution in nondivergent horizontal velocity fields are given which suggest that both the frontogenetic and rotational components of the vector frontogenesis function F may be comparable in developing frontal zones observed in nature. The relative importance of the magnitude and direction contributions to F is quantitatively investigated, and it is found that the Lagrangian rates of change of the magnitude and direction of the potential temperature gradient are comparable. The frontal circulation is found to be related to the magnitude component of the Q vector, whereas the background circulation is related to the direction component.

Keyser, Daniel↗
NASA NTRS1993

The quantum dynamics of electronically nonadiabatic chemical reactions

Considerable progress was achieved on the quantum mechanical treatment of electronically nonadiabatic collisions involving energy transfer and chemical reaction in the collision of an electronically excited atom with a molecule. In the first step, a new diabatic representation for the coupled potential energy surfaces was created. A two-state diabatic representation was developed which was designed to realistically reproduce the two lowest adiabatic states of the valence bond model and also to have the following three desirable features: (1) it is more economical to evaluate; (2) it is more portable; and (3) all spline fits are replaced by analytic functions. The new representation consists of a set of two coupled diabatic potential energy surfaces plus a coupling surface. It is suitable for dynamics calculations on both the electronic quenching and reaction processes in collisions of Na(3p2p) with H2. The new two-state representation was obtained by a three-step process from a modified eight-state diatomics-in-molecules (DIM) representation of Blais. The second step required the development of new dynamical methods. A formalism was developed for treating reactions with very general basis functions including electronically excited states. Our formalism is based on the generalized Newton, scattered wave, and outgoing wave variational principles that were used previously for reactive collisions on a single potential energy surface, and it incorporates three new features: (1) the basis functions include electronic degrees of freedom, as required to treat reactions involving electronic excitation and two or more coupled potential energy surfaces; (2) the primitive electronic basis is assumed to be diabatic, and it is not assumed that it diagonalizes the electronic Hamiltonian even asymptotically; and (3) contracted basis functions for vibrational-rotational-orbital degrees of freedom are included in a very general way, similar to previous prescriptions for locally adiabatic functions in various quantum scattering algorithms.

Truhlar, Donald G.↗
NASA NTRS1983

A theoretical study of plasmaspheric hiss generation

An investigation of the generation of plasmaspheric hiss by the whistler mode instability has been carried out. Using combined ray-tracing and growth rate calculations, the path-integrated growth rates of a large number of waves have been calculated using a distribution function. In general, the waves damp out before being magnetospherically reflected, or soon after. This is true of all the waves traced, independent of the distribution functions used for the energetic (resonant) electrons. Thus the growth rates are generally low, since the waves undergo few equatorial transits. A comparison of these results with observed spectra indicates that while some spectral characteristics are in qualitative agreement, there is a major discrepancy in the intensity levels.

Huang, C. Y.↗
NASA NTRS1975

Space Tug avionics definition study. Volume 2: Avionics functional requirements

Flight and ground operational phases of the tug/shuttle system are analyzed to determine the general avionics support functions that are needed during each of the mission phases and sub-phases. Each of these general support functions is then expanded into specific avionics system requirements, which are then allocated to the appropriate avionics subsystems. This process is then repeated at the next lower level of detail where these subsystem requirements are allocated to each of the major components that comprise a subsystem.

Source record↗
NASA NTRS2001

Viscoplastic Model Development to Account for Strength Differential: Application to Aged Inconel 718 at Elevated Temperature

The magnitude of yield and flow stresses in aged Inconel 718 are observed to be different in tension and compression. This phenomenon, called the Strength differential (SD), contradicts the metal plasticity axiom that the second deviatoric stress invariant alone is sufficient for representing yield and flow. Apparently, at least one of the other two stress invariants is also significant. A unified viscoplastic model was developed that is able to account for the SD effect in aged Inconel 718. Building this model involved both theory and experiments. First, a general threshold function was proposed that depends on all three stress invariants and then the flow and evolution laws were developed using a potential-based thermodynamic framework. Judiciously chosen shear and axial tests were conducted to characterize the material. Shear tests involved monotonic loading, relaxation, and creep tests with different loading rates and load levels. The axial tests were tension and compression tests that resulted in sufficiently large inelastic strains. All tests were performed at 650 C. The viscoplastic material parameters were determined by optimizing the fit to the shear tests, during which the first and the third stress invariants remained zero. The threshold surface parameters were then fit to the tension and compression test data. An experimental procedure was established to quantify the effect of each stress invariant on inelastic deformation. This requires conducting tests with nonproportional three-dimensional load paths. Validation of the model was done using biaxial tests on tubular specimens of aged Inconel 718 using proportional and nonproportional axial-torsion loading. These biaxial tests also helped to determine the most appropriate form of the threshold function; that is, how to combine the stress invariants. Of the set of trial threshold functions, the ones that incorporated the third stress invariant give the best predictions. However, inclusion of the first stress invariant does not significantly improve the model predictions. The model shows excellent predictive capability for nonproportional load paths. Additionally, it reduces to the well-known models of Mises Drucker and Drucker-Prager. The requisite experiments involve reasonably simple load paths in the axial-shear stress plane and hence can be performed on a variety of different materials: be they metallic, geological. polymeric, ceramic or granular. The general form of the threshold function allows representation of inelastic deformation in a range of materials.

Iyer, Saiganesh↗
NASA NTRS1996

The Acoustic Influence of Cell Depth on the Rotordynamic Characteristics of Smooth-Rotor/Honeycomb-Stator Annular Gas Seals

A two-control volume is employed for honeycomb-stator/smooth-rotor seals, with a conventional control-volume used for the through flow and a 'capacitance accumulator' model for the honeycomb cells. The control volume for the honeycomb cells is shown to cause a dramatic reduction in the effective acoustic velocity of the main flow, dropping the lowest acoustic frequency into the frequency range of interest for rotordynamics. In these circumstances, the impedance functions for the seals can not be modeled with conventional (frequency-independent) stiffness, damping, and mass coefficients. More general transfer functions are required to account for the reaction forces, and calculated here as a lead-lag term for the direct force function and a lag term for the cross-coupled function. These first order functions are simple compared to transfer functions for magnetic bearings or foundations, For synchronous response to imbalance, they can be approximated by running-speed-dependent stiffness and damping coefficients in conventional rotordynamic codes. Correct predictions for stability and transient response will require more general algorithms, pressumably using a state-space format.

Childs, Dara W.↗
NASA NTRS1975

On the Formulation of the Gravitational Potential in Terms of Equinoctial Variables

Analytical averaging techniques are used to expand the disturbing potential in the equinoctial coordinate frame by considering third body harmonics and zonal functions harmonics. General results are developed through applications of Legendre and associated Legendre polynomials and the Q sub nm functions for the gravitational potential.

Cefola, P. J.↗
DOE OSTI2023

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗