Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “general relativity equations and 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 163 records · Page 9

Remote sensing of chlorophyll concentration: State-of-the-art, 1975

Remote measurement of chlorophyll concentration of the world's oceans from satellite observations could potentially be extremely useful for assessments of productivity in large areas for which measurements by other means would be impractical. The basis of these measurements rests with the physics of the interaction of light with material dissolved and suspended in the water. It is theoretically possible to predict the nature of light upwelled from the ocean surface from a solution to the radiative transfer equation. Practically, however, this is difficult. Monte-Carlo methods presently are thought to be the most viable method to treat the general theoretical problem. With restrictive assumptions of the nature of scattering, it is possible to construct simpler models. Algorithms developed to relate chlorophyll concentration (or some other parameter, i.e., seechi depth) to the upwelled light spectrum are discussed.

Atwell, B. H.↗

The origins of cosmic rays and quantum effects on gravity

The energy spectrum of primary cosmic rays is explained by particles emitted during a thermal expansion of explosive objects inside and near the galaxy, remnants of which may be supernova and/or active talaxies, or even stars or galaxies that disappeared from our sight after the explosion. A power law energy spectrum for cosmic rays, E to the (-alpha -1, is obtained from an expansion rate T is proportional to R to the alpha. Using the solution of the Einstein equation, we obtain a spectrum which agrees very well with experimental data. The implication of an inflationary early universe on the cosmic ray spectrum is also discussed. It is also suggested that the conflict between this model and the singularity theorem in classical general relativity may be eliminated by quantum effects.

Tomozawa, Y.↗

Numerical methods for turbulent flow

It has generally become accepted that the Navier-Strokes equations predict the dynamic behavior of turbulent as well as laminar flows of a fluid at a point in space away form a discontinuity such as a shock wave. Turbulence is also closely related to the phenomena of non-uniqueness of solutions of the Navier-Strokes equations. These second order, nonlinear partial differential equations can be solved analytically for only a few simple flows. Turbulent flow fields are much to complex to lend themselves to these few analytical methods. Numerical methods, therefore, offer the only possibility of achieving a solution of turbulent flow equations. In spite of recent advances in computer technology, the direct solution, by discrete methods, of the Navier-Strokes equations for turbulent flow fields is today, and in the foreseeable future, impossible. Thus the only economically feasible way to solve practical turbulent flow problems numerically is to use statistically averaged equations governing mean-flow quantities. The objective is to study some recent developments relating to the use of numerical methods to study turbulent flow.

Turner, James C., Jr.↗

Solving the Grad–Shafranov equation using spectral elements for tokamak equilibrium with toroidal rotation

The Grad–Shafranov equation is solved using spectral elements for tokamak equilibrium with toroidal rotation. The Grad–Shafranov solver builds upon and extends the NIMEQ code (Howell and Sovinec, 2014) previously developed for static tokamak equilibria. Both geometric and algebraic convergence are achieved as the polynomial degree of the spectral-element basis increases. A new analytical solution to the Grad–Shafranov equation is obtained for Solov’ev equilibrium in presence of rigid toroidal rotation, in addition to a previously obtained analytical solution for a different set of equilibrium and rotation profiles. The numerical solutions from the extended NIMEQ are benchmarked with the analytical solutions, with good agreements. Besides, the extended NIMEQ code is benchmarked with the FLOW code (Guazzotto et al., 2004). The modification of pressure profile induced by toroidal flow is investigated. The relative change of pressure profile is found significant around the edge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Purfication kinetics of beryllium during vacuum induction melting

The kinetics of evaporation in binary alloys were quantitatively treated. The formalism so developed works well for several systems studied. The kinetics of purification of beryllium was studied through evaporation data actually acquired during vacuum induction melting. Normal evaporation equations are shown to be generally valid and useful for understanding the kinetics of beryllium purification. The normal evaporation analysis has been extended to cover cases of limited liquid diffusion. It was shown that under steady-state evaporation, the solute concentration near the surface may be up to six orders of magnitude different from the bulk concentration. Corrections for limited liquid diffusion are definitely needed for the highly evaporative solute elements, such as Zn, Mg, and Na, for which the computed evaporation times are improved by five orders of magnitude. The commonly observed logarithmic relation between evaporation time and final concentration further supports the validity of the normal evaporation equations.

Mukherjee, J. L.↗

Drift-ordered fluid vorticity equation with energy consistency

Although drift-ordered fluid models are widely applied in tokamak edge turbulence simulations, the models used are acknowledged not to conserve energy or even electrical charge. The present paper aims to remove many of the existing pitfalls in drift-fluid models, however, with the objective of finding a solution simple enough to be implemented in numerical applications. Our main result is an improved version of the drift-Braginskii equations involving a generalized vorticity function. In the new drift-Braginskii system, the quasi-neutrality condition translates into a transport equation for a generalized vorticity, expressed in conservation form, and related to the total mass-weighted circulation. It is found that kinetic energy conservation can be achieved if the polarization flow is defined recursively. In conclusion, the resulting model conserves the kinetic energy associated with E×B and diamagnetic flows and retains the associated perpendicular kinetic energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Unsteady subsonic compressible flow around finite thickness wings.

A general formulation for the unsteady subsonic compressible potential flow around aircraft having arbitrary configurations is presented. An integral representation of the velocity potential is obtained. From this a linear integral equation relating the perturbation potential and its normal derivative (which is known from the boundary conditions) is derived. For the numerical solution of the integral equation, the surface of the aircraft is divided into small elements and the potential is assumed to be constant within each element. Numerical results are obtained for an oscillating finite-thickness wing and indicate good convergence and excellent agreement with existing lifting surface solutions.

Morino, L.↗

A general multiblock Euler code for propulsion integration. Volume 1: Theory document

A general multiblock Euler solver was developed for the analysis of flow fields over geometrically complex configurations either in free air or in a wind tunnel. In this approach, the external space around a complex configuration was divided into a number of topologically simple blocks, so that surface-fitted grids and an efficient flow solution algorithm could be easily applied in each block. The computational grid in each block is generated using a combination of algebraic and elliptic methods. A grid generation/flow solver interface program was developed to facilitate the establishment of block-to-block relations and the boundary conditions for each block. The flow solver utilizes a finite volume formulation and an explicit time stepping scheme to solve the Euler equations. A multiblock version of the multigrid method was developed to accelerate the convergence of the calculations. The generality of the method was demonstrated through the analysis of two complex configurations at various flow conditions. Results were compared to available test data. Two accompanying volumes, user manuals for the preparation of multi-block grids (vol. 2) and for the Euler flow solver (vol. 3), provide information on input data format and program execution.

Chen, H. C.↗

Gravitational wave propagation beyond general relativity: waveform distortions and echoes

We study the cosmological propagation of gravitational waves (GWs) beyond general relativity (GR) across homogeneous and isotropic backgrounds. We consider scenarios in which GWs interact with an additional tensor field and use a parametrized phenomenological approach that generically describes their coupled equations of motion. We analyze four distinct classes of derivative and non-derivative interactions: mass, friction, velocity, and chiral. We apply the WKB formalism to account for the cosmological evolution and obtain analytical solutions to these equations. We corroborate these results by analyzing numerically the propagation of a toy GW signal. We then proceed to use the analytical results to study the modified propagation of realistic GWs from merging compact binaries, assuming that the GW signal emitted is the same as in GR. We generically find that tensor interactions lead to copies of the originally emitted GW signal, each one with its own possibly modified dispersion relation. These copies can travel coherently and interfere with each other leading to a scrambled GW signal, or propagate decoherently and lead to echoes arriving at different times at the observer that could be misidentified as independent GW events. Depending on the type of tensor interaction, the detected GW signal may exhibit amplitude and phase distortions with respect to a GW waveform in GR, as well as birefringence effects. Furthermore, we discuss observational probes of these tensor interactions with both individual GW events, as well as population studies for both ground- and space-based detectors.

79 ASTRONOMY AND ASTROPHYSICS↗

Asymptotic Expansions and Estimates for the Capillary Problem

The asymptotic properties for the small Bond number B of the equilibrium capillary interface interior to a circular cylindrical tube vertically dipped in an infinite reservoir of liquid are discussed. (The Bond number B is a dimensionless parameter which is the ratio of gravitational to capillary forces.) The formal expansion in powers of B of the solution to the differential equation describing the equilibrium surface (as can be obtained by standard perturbation methods) is proved to be truly asymptotic -- to all orders and uniformly in the variable and parameter gamma, the contact angle. Sequences of general estimates, in closed form, from above and from below, are also given for the solution and related functions. The M-th term in these sequences are asymptotically exact to order m. An idiosyncrasy of the problem, crucial in obtaining these estimates, is the absolute monotonicity of the structural function of the system in integral form.

Brulois, F. P.↗

Statistical-mechanical selection of the shapes of disk galaxies

A new method is proposed for selecting steady state shapes of disk galaxies as 'most probable states' of a large number of stars, given only the total energy and total angular momentum. A partial differential equation is derived for the mean gravitational potential; it is closely related to the 'sinh-Poisson' equation for the mean-field description of a line vortex system or electrostatic guiding-center plasma. A 'water bag' approximation to the distribution function for bound stars renders the equation analytically tractable, but accurate solution of it may require numerical integration, in view of its general nonlinearity.

Montgomery, David↗

Application of two-point implicit central-difference methods to hyperbolic systems

This paper presents a general solution algorithm for the set of difference equations that arise when two-point central differences are used to approximate the flux difference terms in systems of hyperbolic differential equations. The general algorithm eliminates the weak points associated with the nonstandard algorithm reported by Wornom and Hafez (1986). The disadvantages of their algorithm relate to its implementation. It consists of separate algorithms for subsonic, supersonic, sonic and shock cells, applied individually, which presents a major bookkeeping problem when multiple sonic and shock cells are present. The general algorithm eliminates this problem and introduces an improved shock treatment which produces shocks with at most one interior shock point.

Wornom, Stephen F.↗

A note on relative motion in the general three-body problem.

It is shown that the equations of the general three-body problem take on a very symmetric form when one considers only their relative positions, rather than position vectors relative to some given coordinate system. From these equations one quickly surmises some well known classical properties of the three-body problem, such as the first integrals and the equilateral triangle solutions. Some new Lagrangians with relative coordinates are also obtained. Numerical integration of the new equations of motion is about 10% faster than with barycentric or heliocentric coordinates.

Broucke, R.↗

A general algorithm for relating ground trajectory distance, elapsed flight time, and aircraft airspeed and its application to 4-D guidance

A general solution using an elliptic integral approximation which relates flight time, aircraft airspeed, and ground distance on straight-line and circular-arc trajectory segments is developed. The solution procedure is applicable to both constant and accelerating aircraft flight. In addition, wind shear including both magnitude and heading change is incorporated in the solution. The solution equations are used in a four-dimensional (4-D) control algorithm where both flight time and final airspeed are specified. The results show that the algorithm converges rapidly and accurately.

Foudriat, E. C.↗

Nonlinear rotordynamics analysis

Three analytic consequences of the nonlinear Jeffcott equations are examined. The primary application of these analyses is directed toward understanding the excessive vibrations recorded in the Liquid Oxygen (LOX) pump of the Space Shuttle Main Engine (SSME) during hot firing ground testing. The first task is to provide bounds on the coefficients of the equations which delimit the two cases of numerical solution as a circle or an annulus. The second task examines the mathematical generalization to multiple forcing functions, which includes the special problems of mass imbalance, side force, rubbing, and combination of these forces. Finally, stability and boundedness of the steady-state solutions is discussed and related to the corresponding linear problem.

Day, W. B.↗

Surfactant Facilitated Spreading of Aqueous Drops on Hydrophobic Surfaces

Microgravity technologies often require aqueous phases to spread over nonwetting hydrophobic solid surfaces. Surfactants facilitate the wetting of water on hydrophobic surfaces by adsorbing on the water/air and hydrophobic solid/water interfaces and lowering the surface tensions of these interfaces. The tension reductions decrease the contact angle, which increases the equilibrium wetted area. Hydrocarbon surfactants; (i.e., amphiphiles with a hydrophobic moiety consisting of an extended chain of (aliphatic) methylene -CH2- groups attached to a large polar group to give aqueous solubility) are capable of reducing the contact angles on surfaces which are not very hydrophobic, but do not reduce significantly the contact angles of the very hydrophobic surfaces such as parafilm, polyethylene or self assembled monolayers. Trisiloxane surfactants (amphiphiles with a hydrophobe consisting of methyl groups linked to a trisiloxane backbone in the form of a disk ((CH3)3-Si-O-Si-O-Si(CH3)3) and an extended ethoxylate (-(OCH2CH2)a-) polar group in the form of a chain with four or eight units) can significantly reduce the contact angle of water on a very hydrophobic surface and cause rapid and complete (or nearly complete) spreading (termed superspreading). The overall goal of the research described in this proposal is to establish and verify a theory for how trisiloxanes cause superspreading, and then use this knowledge as a guide to developing more general hydrocarbon based surfactant systems which superspread. We propose that the trisiloxane surfactants superspread because their structure allows them to strongly lower the high hydrophobic solid/aqueous tension when they adsorb to the solid surface. When the siloxane adsorbs, the hydrophobic disk parts of the molecule adsorb onto the surface removing the surface water. Since the cross-sectional area of the disk is larger than that of the extended ethoxylate chain, the disks can form a space-filling mat on the surface which removes a significant amount of the surface water. In this presentation, we report the results of measurements of the molecular packing and rates of kinetic exchange of the trisiloxane surfactants at the air/water interface in order to confirm our picture of trisiloxane packing, and provide additional insight into the superspreading process. We used the pendant bubble technique as a Langmuir trough to measure the trisiloxane equation of state which relates the tension to the surface concentration. From these measurements we obtain accurate values for the maximum packing density. We find that trisiloxanes with 4 and 8 ethoxylate groups have the same maximum packing concentration, indicating that the maximum packing is controlled by the cross section of the head group. For trisiloxanes with larger than eight ethoxylates, the maximum packing increases with ethoxylate number, indicating that the disposition of the ethoxylate chain (i.e., its effective size) is controlling. This supports our picture of superspreading: The superspreading ability of trisiloxanes decreases considerably for trisiloxanes with larger than eight ethoxylates; the packing measurements indicate that for the higher ethoxylate number trisiloxanes, the compact nonpolar head groups are pushed apart by the ethoxylate chain. They leave spaces of surface water on adsorption and do not lower the solid tension as much as their lower chain analogues. Finally the report measurements of the dynamic tension reduction accompanying the adsorption of trisiloxanes onto an initially clean interface using the pendant bubble technique, and we obtain from these relaxations, the equation of state and a mass transfer model, the rate constants for kinetic exchange. We find that the rate constants for desorption of trisiloxanes are generally much slower than for analogous aliphatic polyethoxylate surfactants with identical ethoxylate chain lengths. When an aqueous drop of a superspreader solution is placed on a hydrophobic surface and begins to spread, the lower desorption rates allows the tension at the drop center to remain reduced relative to the tension of the expanding periphery, thus strengthening Marangoni forces which can assist the spreading. Marangoni forces can be especially significant in the case of superspreaders because their maximum equilibrium reduction in tension is several dynes/cm lower than for aliphatic surfactants.

Kumar, Nitin↗

Investigation of the three-dimensional flow field within a transonic fan rotor: Experiment and analysis

An experimental investigation of the three dimensional flow field through a low aspect ratio, transonic, axial flow fan rotor has been conducted using an advanced laser anemometer (LA) system. Laser velocimeter measurements of the rotor flow field at the design operating speed and over a range of through flow conditions are compared to analytical solutions. The numerical technique used herein yields the solution to the full, three dimensional, unsteady Euler equations using an explicit time marching, finite volume approach. The numerical analysis, when coupled with a simplified boundary layer calculation, generally yields good agreement with the experimental data. The test rotor has an aspect ratio of 1.56, a design total pressure ratio of 1.629 and a tip relative Mach number of 1.38. The high spatial resolution of the LA data matrix (9 radial by 30 axial by 50 blade to blade) permits details of the transonic flow field such as shock location, turning distribution and blade loading levels to be investigated and compared to analytical results.

Pierzga, M. J.↗

Interaction of transient radiation in nongray gaseous systems

A general formulation is presented to investigate the transient radiative interaction in nongray absorbing-emitting species between two parallel plates. Depending on the desired sophistication and accuracy, any nongray absorption model from line-by-line models to the wide band model correlations can be employed in the formulation to investigate the radiative interaction. Special attention is directed to investigate the radiative interaction in a system initially at a uniform reference temperature and suddenly the temperature of the bottom plate is reduced to a lower but constant temperature. The interaction is considered for the case of radiative equilibrium as well as for combined radiation and conduction. General as well as limiting forms of the governing equations are presented and solutions are obtained numerically by employing the method of variation of parameters. Specific results are obtained for CO, CO2, H2O, and OH. The information on species H2O and OH is of special interest for the proposed scramjet engine application. The results demonstrate the relative ability of different species for radiative interactions.

Tiwari, S. N.↗