Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “General relativity equations & solutions”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

On the tesseral-harmonics resonance problem in artificial-satellite theory

The longitude-dependent part of the geopotential usually gives rise only to short-period effects in the motion of an artificial satellite. However, when the motion of the satellite is commensurable with that of the earth, the path of the satellite repeats itself relative to the earth and perturbations build up at each passage of the satellite in the same spot, so that there can be important long-period effects. In order to take these effects into account in deriving a theoretical solution to the equations of motion of an artificial satellite, it is necessary to select terms in the longitude-dependent part of the geopotential that will contribute significantly to the perturbations. Attempts made to obtain a selection that is valid in a general case, regardless of the initial eccentricity of the orbit and of the order of the resonance, are reported. The solution to the equations of motion of an artificial satellite, in a geopotential thus determined, is then derived by using Hori's method by Lie series, which, by its properties regarding canonical invariance, has proved advantageous in the classical theory.

Romanowicz, B. A.↗
NASA NTRS1981

Evolution of magnetic field and atmospheric response. I - Three-dimensional formulation by the method of projected characteristics. II - Formulation of proper boundary equations

The method described as the method of nearcharacteristics by Nakagawa (1980) is renamed the method of projected characteristics. Making full use of properties of the projected characteristics, a new and simpler formulation is developed. As a result, the formulation for the examination of the general three-dimensional problems is presented. It is noted that since in practice numerical solutions must be obtained, the final formulation is given in the form of difference equations. The possibility of including effects of viscous and ohmic dissipations in the formulation is considered, and the physical interpretation is discussed. A systematic manner is then presented for deriving physically self-consistent, time-dependent boundary equations for MHD initial boundary problems. It is demonstrated that the full use of the compatibility equations (differential equations relating variations at two spatial locations and times) is required in determining the time-dependent boundary conditions. In order to provide a clear physical picture as an example, the evolution of axisymmetric global magnetic field by photospheric differential rotation is considered.

Nakagawa, Y.↗
NASA NTRS1999

High Order Approximations for Compressible Fluid Dynamics on Unstructured and Cartesian Meshes

The development of high-order accurate numerical discretization techniques for irregular domains and meshes is often cited as one of the remaining challenges facing the field of computational fluid dynamics. In structural mechanics, the advantages of high-order finite element approximation are widely recognized. This is especially true when high-order element approximation is combined with element refinement (h-p refinement). In computational fluid dynamics, high-order discretization methods are infrequently used in the computation of compressible fluid flow. The hyperbolic nature of the governing equations and the presence of solution discontinuities makes high-order accuracy difficult to achieve. Consequently, second-order accurate methods are still predominately used in industrial applications even though evidence suggests that high-order methods may offer a way to significantly improve the resolution and accuracy for these calculations. To address this important topic, a special course was jointly organized by the Applied Vehicle Technology Panel of NATO's Research and Technology Organization (RTO), the von Karman Institute for Fluid Dynamics, and the Numerical Aerospace Simulation Division at the NASA Ames Research Center. The NATO RTO sponsored course entitled "Higher Order Discretization Methods in Computational Fluid Dynamics" was held September 14-18, 1998 at the von Karman Institute for Fluid Dynamics in Belgium and September 21-25, 1998 at the NASA Ames Research Center in the United States. During this special course, lecturers from Europe and the United States gave a series of comprehensive lectures on advanced topics related to the high-order numerical discretization of partial differential equations with primary emphasis given to computational fluid dynamics (CFD). Additional consideration was given to topics in computational physics such as the high-order discretization of the Hamilton-Jacobi, Helmholtz, and elasticity equations. This volume consists of five articles prepared by the special course lecturers. These articles should be of particular relevance to those readers with an interest in numerical discretization techniques which generalize to very high-order accuracy. The articles of Professors Abgrall and Shu consider the mathematical formulation of high-order accurate finite volume schemes utilizing essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) reconstruction together with upwind flux evaluation. These formulations are particularly effective in computing numerical solutions of conservation laws containing solution discontinuities. Careful attention is given by the authors to implementational issues and techniques for improving the overall efficiency of these methods. The article of Professor Cockburn discusses the discontinuous Galerkin finite element method. This method naturally extends to high-order accuracy and has an interpretation as a finite volume method. Cockburn addresses two important issues associated with the discontinuous Galerkin method: controlling spurious extrema near solution discontinuities via "limiting" and the extension to second order advective-diffusive equations (joint work with Shu). The articles of Dr. Henderson and Professor Schwab consider the mathematical formulation and implementation of the h-p finite element methods using hierarchical basis functions and adaptive mesh refinement. These methods are particularly useful in computing high-order accurate solutions containing perturbative layers and corner singularities. Additional flexibility is obtained using a mortar FEM technique whereby nonconforming elements are interfaced together. Numerous examples are given by Henderson applying the h-p FEM method to the simulation of turbulence and turbulence transition.

Barth, Timothy↗
NASA NTRS1977

Application of the method of integral relations to laminar boundary layers in three dimensions

The method of integral relations is extended to general three-dimensional compressible laminar boundary layer flows. The transformation employed to transform the basic three-dimensional compressible boundary layer equations into quasi-incompressible form is an extension of the Howarth transformation. The resulting system of differential equations is integrated numerically by the method of integral relations as proposed by Dorodnitsyn. To demonstrate the accuracy of the method, it is applied to calculation of the parabolic flow over a flat plate and the boundary flow over an infinite yawed cylinder, for which solutions are known. It is then applied to the flow over a flat plate disturbed by a cylinder normal to the plate, for which a finite-difference solution is available for comparison. It is finally applied to calculating the crossflow velocity variation for supersonic flow over a swept wedge.

Holt, M.↗
NASA NTRS1989

Rotordynamic coefficients for stepped labyrinth gas seals

The basic equations are derived for compressible flow in a stepped labyrinth gas seal. The flow is assumed to be completely turbulent in the circumferential direction where the friction factor is determined by the Blasius relation. Linearized zeroth and first-order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth-order pressure distribution is found by satisfying the leakage equation while the circumferential velocity distribution is determined by satisfying the momentum equations. The first order equations are solved by a separation of variables solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients. The results of this analysis are presented in the form of a parametric study, since there are no known experimental data for the rotordynamic coefficients of stepped labyrinth gas seals. The parametric study investigates the relative rotordynamic stability of convergent, straight and divergent stepped labyrinth gas seals. The results show that, generally, the divergent seal is more stable, rotordynamically, than the straight or convergent seals. The results also show that the teeth-on-stator seals are not always more stable, rotordynamically, then the teeth-on-rotor seals as was shown by experiment by Childs and Scharrer (1986b) for a 15 tooth seal.

Scharrer, Joseph K.↗
NASA NTRS1989

Transition regime droplet growth and evaporation - An integrodifferential variational approach

A variational principle based on the integrodifferential form of the Boltzmann equation, which allows very general forms of the collision integral and of the boundary operator, has been used to compute the mass flux to a spherical particle. For simplicity the problem has been restricted to the one-speed, constant cross-section approximation, the black sphere problem of neutron transport theory. Using the Hilbert expansion of the solution far from the sphere as a trial function leads to a rational expression in the inverse Knudsen number for the extrapolation distance, which is simply related to the number flux to the sphere. A five-term expansion gives results that are no more than 6 percent in error when compared to an accurate numerical solution. The technique can be applied to other forms of the collision model and other boundary conditions.

Cipolla, J. W., Jr.↗
DOE OSTI2023

Scattering in black hole backgrounds and higher-spin amplitudes. Part II

We continue to investigate correspondences between, on the one hand, scattering amplitudes for massive higher-spin particles and gravitons in appropriate quantum-to-classical limits, and on the other hand, classical gravitational interactions of spinning black holes according to general relativity. We first construct an ansatz for a gravitational Compton amplitude, at tree level, constrained only by locality, crossing symmetry, unitarity and consistency with the linearized-Kerr 3-point amplitude, to all orders in the black hole’s spin. We then explore the extent to which a unique classical Compton amplitude can be identified by comparing with the results of the classical process of scattering long-wavelength gravitational waves off an exact Kerr black hole, determined by appropriate solutions of the Teukolsky equation. Up to fourth order in spin, we find complete agreement with a previously conjectured exponential form of the tree-level Compton amplitude. At higher orders, we extract tree-level contributions from the Teukolsky amplitude by an analytic continuation from a physical (a/GM < 1) to a particle-like (a/GM > 1) regime. Up to the sixth order in spin, we identify a unique conservative part of the amplitude which is insensitive both to the choice of boundary conditions at the black hole horizon and to branch choices in the analytic continuation. The remainder of the amplitude is determined modulo an overall sign from a branch choice, with the sign flipping under exchanging purely ingoing and purely outgoing boundary conditions at the horizon. Along the way, we make contact with novel applications of massive spinor-helicity variables pertaining to their relation to EFT operators and (spinning) partial amplitudes.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1983

A short cut integration scheme to determine the dynamic response of a launch vehicle with several payloads

The coupled base motion equations of a system composed of a launch vehicle and several payloads are derived. In previous work it was shown how the special form of these equations allow for a quick and accurate solution both for the response and the loads of the system. In this paper a short-cut version of this direct integration technique is presented and evaluated. The result is a method which, depending on the nature of the structure at hand, vasilates between a full-scale coupled base motion approach and an open loop base drive method. Although completely general, the presented technique is most effective when the number of interface degrees of freedom is relatively small compared to the overall number of degrees of freedom of the system. The S.T.S.-S.T.-O.M.S. Kit system is used as an illustrative example.

Engels, R. C.↗
NASA NTRS1982

A formulation of rotor-airframe coupling for design analysis of vibrations of helicopter airframes

A linear formulation of rotor airframe coupling intended for vibration analysis in airframe structural design is presented. The airframe is represented by a finite element analysis model; the rotor is represented by a general set of linear differential equations with periodic coefficients; and the connections between the rotor and airframe are specified through general linear equations of constraint. Coupling equations are applied to the rotor and airframe equations to produce one set of linear differential equations governing vibrations of the combined rotor airframe system. These equations are solved by the harmonic balance method for the system steady state vibrations. A feature of the solution process is the representation of the airframe in terms of forced responses calculated at the rotor harmonics of interest. A method based on matrix partitioning is worked out for quick recalculations of vibrations in design studies when only relatively few airframe members are varied. All relations are presented in forms suitable for direct computer implementation.

Kvaternik, R. G.↗
NASA NTRS2022

Simultaneous Stoquasticity

Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.

Jacob Bringewatt↗
NASA NTRSDate not supplied

Simultaneous Stoquasticity

Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.

Monte Carlo↗
NASA NTRS1978

New formulation of de Sitter's theory of motion for Jupiter I-IV. I - Equations of motion and the disturbing function

Elliptic orbits are substituted for circular orbits in the first approximation, in an analysis of the common retrograde motion of Jupiter's satellites. A modification of the de Sitter theory, made possible by extended observations of the satellites, is presented with attention to that aspect of the theory which eliminates small divisors at all stages of the solution. The convergence problem is circumvented by use of Poincare's canonical relative coordinates. In addition, modified Delaunay variables and their associated Poincare variables are applied to the disturbing function, which is expanded by means of generalized Newcomb operators.

Aksnes, K.↗
DOE OSTI2022

Simultaneous stoquasticity

Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Furthermore, solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Conservative-variable average states for equilibrium gas multi-dimensional fluxes

Modern split component evaluations of the flux vector Jacobians are thoroughly analyzed for equilibrium-gas average-state determinations. It is shown that all such derivations satisfy a fundamental eigenvalue consistency theorem. A conservative-variable average state is then developed for arbitrary equilibrium-gas equations of state and curvilinear-coordinate fluxes. Original expressions for eigenvalues, sound speed, Mach number, and eigenvectors are then determined for a general average Jacobian, and it is shown that the average eigenvalues, Mach number, and eigenvectors may not coincide with their classical pointwise counterparts. A general equilibrium-gas equation of state is then discussed for conservative-variable computational fluid dynamics (CFD) Euler formulations. The associated derivations lead to unique compatibility relations that constrain the pressure Jacobian derivatives. Thereafter, alternative forms for the pressure variation and average sound speed are developed in terms of two average pressure Jacobian derivatives. Significantly, no additional degree of freedom exists in the determination of these two average partial derivatives of pressure. Therefore, they are simultaneously computed exactly without any auxiliary relation, hence without any geometric solution projection or arbitrary scale factors. Several alternative formulations are then compared and key differences highlighted with emphasis on the determination of the pressure variation and average sound speed. The relevant underlying assumptions are identified, including some subtle approximations that are inherently employed in published average-state procedures. Finally, a representative test case is discussed for which an intrinsically exact average state is determined. This exact state is then compared with the predictions of recent methods, and their inherent approximations are appropriately quantified.

Iannelli, G. S.↗
NASA NTRS2002

Analytical Modal Analysis for Thin-Film Flat Lenses

Due to strong potential applications and more demanding requirements imposed upon thin-film structures for space deployable, there has been increasing research and development activities during recent years in the field of vibration analysis of these types of structures. Moreover, interests in employing these structural components have received renewed emphasis in recent years within NASA and the Air Force. This is due to their inherent lightweight, low packaging and launch volume, and relative simplicity of deployment. Among the potential mission concepts for which these structural elements are included, one can mention solar sails, space solar power generation systems, solar thermal propulsion vehicles, large space telescopes, and inflatable communication antennas. This paper presents analytical procedures to determine vibration and physical characteristics of thin film lenses with circular and elliptical shapes membranes considered in design of a solar concentrator. In general, three methods are used to obtain approximate solutions of Helmholtz boundary value problems. One method requires that solution satisfy the differential equation exactly and the boundary condition approximately. Another method demands a solution that satisfies the boundary conditions exactly and the governing equations approximately. The third method sees a solution that satisfies both the governing equation and boundary conditions approximately. Extensive reviews of vibrations of membrane and plates are provided by Leissa and Mazumdar.

Hamid R. Hamidzadeh↗
DOE OSTI2024

Simplified Universal Equations for Ionic Conductivity and Transference Number

Nernst-Einstein equation can provide a reasonable estimate of the ionic conductivity of dilute solutions. For concentrated solutions, alternate methods such as Green–Kubo relations and Einstein relations are more suitable to account for ion-ion interactions. Such computations can be expensive for multicomponent systems. Simplified mathematical expressions like the Nernst-Einstein equation do not exist for concentrated multicomponent mixtures. Newman’s treatment of multicomponent concentrated solutions yields a conductivity relation in terms of species concentration and Onsager phenomenological coefficients. However, the estimation of these phenomenological coefficients is not straightforward. Here, mathematical formulations that relate the phenomenological coefficients with the friction coefficients are developed, leading to simplified, ready-to-use expressions of conductivity and transference numbers that can be used for a wide range of ionic mixtures. This approach involves spectral decomposition of the matrix of Onsager phenomenological coefficients. The general analytical expressions for conductivity and transference number are simplified for binary electrolytes, and numerical solutions are provided for ternary and quaternary mixtures with ion dissociation.

Electrochemistry↗
NASA NTRS1977

Scale-covariant theory of gravitation and astrophysical applications

A scale-covariant theory of gravitation is presented which is characterized by a set of equations that are complete only after a choice of the scale function is made. Special attention is given to gauge conditions and units which allow gravitational phenomena to be described in atomic units. The generalized gravitational-field equations are derived by performing a direct scale transformation, by extending Riemannian geometry to Weyl geometry through the introduction of the notion of cotensors, and from a variation principle. Modified conservation laws are provided, a set of dynamical equations is obtained, and astrophysical consequences are considered. The theory is applied to examine certain homogeneous cosmological solutions, perihelion shifts, light deflections, secular variations of planetary orbital elements, stellar structure equations for a star in quasi-static equilibrium, and the past thermal history of earth. The possible relation of the scale-covariant theory to gauge field theories and their predictions of cosmological constants is discussed.

Canuto, V.↗
NASA NTRS1991

A fractional step solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

The time-dependent, three-dimensional incompressible Navier-Stokes equations are presently solved in generalized coordinate systems by means of a fractional-step method whose primitive variable formulation uses as dependent variables, in place of the Cartesian components of the velocity: (1) pressure (defined at the center of the computational cell), and (2) volume fluxes across the faces of the cells. The momentum equations are solved by means of an approximate factorization method. A novel 'ZEBRA' scheme incorporating four-color ordering efficiently solves the Poisson equation. Illustrative two- and three-dimensional laminar flow test cases are computed and evaluated relative to extant numerical and experimental results, and good agreement is obtained.

Rosenfeld, Moshe↗