Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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 1,027 records · Page 57

Improvement of the Earth's gravity field from terrestrial and satellite data

The determination of the Earth's gravitational potential can be done through the analysis of satellite perturbations, the analysis of surface gravity data, or both. The combination of the two data types yields a solution that combines the strength of each method: the longer wavelength strength in the satellite analysis with the better high frequency information from surface gravity data. Since 1972, Ohio State has carried out activities that have provided surface gravity data to a number of organizations who have developed combination potential coefficient models that describe the Earth's gravitational potential.

Rapp, Richard H.↗

Chemical vapor deposition fluid flow simulation modelling tool

Accurate numerical simulation of chemical vapor deposition (CVD) processes requires a general purpose computational fluid dynamics package combined with specialized capabilities for high temperature chemistry. In this report, we describe the implementation of these specialized capabilities in the spectral element code NEKTON. The thermal expansion of the gases involved is shown to be accurately approximated by the low Mach number perturbation expansion of the incompressible Navier-Stokes equations. The radiative heat transfer between multiple interacting radiating surfaces is shown to be tractable using the method of Gebhart. The disparate rates of reaction and diffusion in CVD processes are calculated via a point-implicit time integration scheme. We demonstrate the use above capabilities on prototypical CVD applications.

Bullister, Edward T.↗

Method for the calculation of spacecraft umbra and penumbra shadow terminator points

A method for calculating orbital shadow terminator points is presented. The current method employs the use of an iterative process which is used for an accurate determination of shadow points. This calculation methodology is required since orbital perturbation effects can introduce large errors when a spacecraft orbits a planet in a high altitude and/or highly elliptical orbit. To compensate for the required iteration methodology, all reference frame change definitions and calculations are performed with quaternions. Quaternion algebra significantly reduces the computational time required for the accurate determination of shadow terminator points.

Ortizlongo, Carlos R.↗

Interplanetary magnetic field orientation for transient events in the outer magnetosphere

It is generally believed that flux transfer events (FTEs) in the outer dayside magneosphere, usually identified by transient (approximately 1 min) bipolar magneitc field perturbations in the direction normal to the nominal magnetopause, occur when the magnetosheath magetic field has a southward component. We compare the results of three methods for determining the magnetosheath magnetic field orientationat the times of previously identified UKS/IRM events: (1) the average magnetosheath magnetic field orientation in the 30-min period adjacent to the nearest magnetopause crossing, (2) the magnetosheath magnetic field orientation observed just outside the magnetopause, and (3) the lagged interplanetary magnetic field (IMF) orientation at the time of the transient events. Whereas the results of method 2 indicate that the events tend to occur for a southward magnetosheath magnetic field, the results of methods 1 and 3 show no such tnedency. The fact that the three methods yield significantly diffeent results emphasizes the need for caution in future studies.

Sibeck, D. G.↗

Gauche Ethyl Alcohol: Laboratory Assignments and Interstellar Identification

Ethyl alcohol (ethanol) is known to possess a pair of closely spaced excited torsional substates (gauche+, gauche-) at an energy of approximately 57 K above the ground (trans) torsional substate. We report an extended analysis of some gauche - gauche+ Q-branch ((Delta)J = 0) transitions with a three-substate fixed frame axis method (FFAM) Hamiltonian. Our approach accounts for complex trans-gauche interactions for the first time. In addition, we are able to obtain intensities for perturbed rotational transitions, and to determine the trans to gauche+ separation to be 1185399.1 MHz. A complete ground state rotational-torsional partition function accounting for the previously neglected gauche substates is presented. Based on our analysis, a total of 14 U lines obtained towards Orion KL can now be assigned to gauche substates of ethanol. Analysis of these lines yields a rotational temperature of 223 K and a total (trans + gauche) column density of 7.0 x 10(exp 15)/sq cm. The column density is in reasonable agreement with the recent value of 2-3 x 10(exp 15)/sq cm based on observations of trans-ethanol by Ohishi et al., although there is some disparity in the rotational temperatures. Eight additional U lines in the literature are assigned to transitions of gauche ethanol.

Pearson, J. C.↗

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗

A Program to Improve the Triangulated Surface Mesh Quality Along Aircraft Component Intersections

A computer program has been developed for improving the quality of unstructured triangulated surface meshes in the vicinity of component intersections. The method relies solely on point removal and edge swapping for improving the triangulations. It can be applied to any lifting surface component such as a wing, canard or horizontal tail component intersected with a fuselage, or it can be applied to a pylon that is intersected with a wing, fuselage or nacelle. The lifting surfaces or pylon are assumed to be aligned in the axial direction with closed trailing edges. The method currently maintains salient edges only at leading and trailing edges of the wing or pylon component. This method should work well for any shape of fuselage that is free of salient edges at the intersection. The method has been successfully demonstrated on a total of 125 different test cases that include both blunt and sharp wing leading edges. The code is targeted for use in the automated environment of numerical optimization where geometric perturbations to individual components can be critical to the aerodynamic performance of a vehicle. Histograms of triangle aspect ratios are reported to assess the quality of the triangles attached to the intersection curves before and after application of the program. Large improvements to the quality of the triangulations were obtained for the 125 test cases; the quality was sufficient for use with an automated tetrahedral mesh generation program that is used as part of an aerodynamic shape optimization method.

Cliff, Susan E.↗

Initial Results of an MDO Method Evaluation Study

The NASA Langley MDO method evaluation study seeks to arrive at a set of guidelines for using promising MDO methods by accumulating and analyzing computational data for such methods. The data are collected by conducting a series of re- producible experiments. In the first phase of the study, three MDO methods were implemented in the SIGHT: framework and used to solve a set of ten relatively simple problems. In this paper, we comment on the general considerations for conducting method evaluation studies and report some initial results obtained to date. In particular, although the results are not conclusive because of the small initial test set, other formulations, optimality conditions, and sensitivity of solutions to various perturbations. Optimization algorithms are used to solve a particular MDO formulation. It is then appropriate to speak of local convergence rates and of global convergence properties of an optimization algorithm applied to a specific formulation. An analogous distinction exists in the field of partial differential equations. On the one hand, equations are analyzed in terms of regularity, well-posedness, and the existence and unique- ness of solutions. On the other, one considers numerous algorithms for solving differential equations. The area of MDO methods studies MDO formulations combined with optimization algorithms, although at times the distinction is blurred. It is important to

Alexandrov, Natalia M.↗

Black hole normal modes - A semianalytic approach

A new semianalytic technique for determining the complex normal mode frequencies of black holes is presented. The method is based on the WKB approximation. It yields a simple analytic formula that gives the real and imaginary parts of the frequency in terms of the parameters of the black hole and of the field whose perturbation is under study, and in terms of the quantity (n + 1/2), where n = 0, 1, 2,... and labels the fundamental mode, first overtone mode, and so on. In the case of the fundamental gravitational normal modes of the Schwarzschild black hole, the WKB estimates agree with numerical results to better than 7 percent in the real part of the frequency and 0.7 percent in the imaginary part, with the relative agreement improving with increasing angular harmonic. Carried to higher order the method may provide an accurate and systematic means to study black hole normal modes.

Schutz, B. F.↗

Monte Carlo solution of Boltzmann equation for a simple model of highly nonequilibrium diatomic gases - Translational rotational energy relaxation

Theoretical results pertaining to internally excited translational-rotational energy relaxation in a spatially uniform diatomic gas far removed from solid boundaries are obtained by solving the Boltzmann equation by means of the Monte Carlo direct simulation method. The analysis is based on calculations involving three different types of initial conditions: equilibrium, nonequilibrium-equipartition (i.e., equipartition is satisfied, but the distributions are perturbed), and nonequilibrium-nonequipartition (i.e., both equipartition and the distributions are perturbed). Results of monatomic-gas simulations are also included to facilitate comparisons with the coupled translational-rotational relaxation simulations, and some simulations for a normal shock-wave structure are briefly examined. The results show that: (1) single-step transitions are the significant mechanisms of intermodal energy transfer; (2) translational-rotational transitions are coupled most efficiently for low-lying states of rotationally excited molecules and least efficiently for highly rotationally excited molecules; and (3) relaxation occurs via a successive set of distributions that are not Maxwell-Boltzmann (nonlocal Maxwellian).

Yoshikawa, K. K.↗

Stability of a dual-spin spacecraft with spherical dampers

The present investigation is concerned with the stability characteristics of a specific dual-spin satellite configuration marked by a high degree of symmetry. The configuration includes a platform and a rotor. Both components contain arbitrarily located internal spherical dampers. The symmetry of the system configuration makes it possible to illustrate clearly the relationship between Routhian analysis, energy sink analysis, and digital simulation of the full nonlinear equations. Although the dual-spin spacecraft configuration contains energy dissipating devices on both platform and rotor, it is still possible to employ the rigorous, but relatively simple, Routh stability method. This method, unlike Floquet theory, has the potential of producing closed-form stability criteria. The energy sink method is capable of providing a closed-form stability criterion. Numerical simulation is a necessary requirement in the latter stages of design when the realistic perturbation environment must be considered.

Laskin, R. A.↗

Astrophysical materials science: Theory

A method of structural expansions for use in determining the equation of state of metallic hydrogen (and indeed other metals) up to the 4th order in the perturbation theory was developed. The electrical and thermal transport properties of the planetary interior of Jupiter were calculated. The nature of the interaction between molecules at short range and the importance of multicenter terms in arriving at an adequate description of the thermodynamic functions of condensed molecular hydrogen were also investigated.

Ashcroft, N. W.↗

Cut-and-connect of two antiparallel vortex tubes

Motivated by an early conjecture that vortex cut-and-connect plays a key role in mixing and production of turbulence, helicity and aerodynamic noise, the cross-linking of two antiparallel viscous vortex tubes via direct numerical simulation is studied. The Navier-Stokes equations are solved by a dealiased pseudo-spectral method with 64 cubed grid points in a periodic domain for initial Reynolds numbers Re up to 1000. The vortex tubes are given an initial sinusoidal perturbation to induce a collision and keep the two tubes pressed against each other as annihilation continues. Cross-sectional and wire plots of various properties depict three stages of evolution: (1) Inviscid induction causing vortex cores to first approach and form a contact zone with a dipole cross-section, and then to flatten and stretch; (2) Vorticity annihilation in the contact zone accompanied by bridging between the two vortices at both ends of the contact zone due to a collection of cross-linked vortex lines, now orthogonal to the initial vortex tubes. The direction of dipole advection in the contact zone reverses; and (3) Threading of the remnants of the original vortices in between the bridges as they pull apart. The crucial stage 2 is shown to be a simple consequence of vorticity annihilation in the contact zone, link-up of the un-annihilated parts of vortex lines, and stretching and advection by the vortex tube swirl of the cross-linked lines, which accumulate at stagnation points in front of the annihilating vortex dipole. It is claimed that bridging is the essence of any vorticity cross-linking and that annihilation is sustained by stretching of the dipole by the bridges. Vortex reconnection details are found to be insensitive to asymmetry. Modeling of the reconnection process is briefly examined. The 3D spatial details of scalar transport (at unity Schmidt number), enstrophy production, dissipation and helicity are also examined.

Melander, Mogens V.↗

An asymptotic description of transient settling and ultrafiltration of colloidal dispersions

The present model for the sedimentation of colloidal systems includes a diffusion term in the governing equation which, in the regions above the sediment, acts as small perturbation to the Kynch (1952) theory. Within the sediment, diffusion is comparable to convection due to the the high solid volume fraction. An application of the method of matched asymptotic expansions to the conservation equation allows a complete description of the settling process to be formulated, with specific attention to volume-fraction evolution in the sediment.

Davis, K. E.↗

Development of non-conventional control methods for high angle of attack flight using vortex manipulation

Low-speed wind tunnel experiments have been conducted to define methods for the manipulation of a generic fighter forebody's vortices in order to generate controlled yawing moments at high alpha. Attention is given to two vortex-control schemes: one involving independently controlled forebody tip strakes, and the other individually controlled blowing ports on the forebody. The most effective methods for forebody yawing-moment control were found to be the minimization of natural asymmetry with a pair of symmetrically mounted tip strakes, and the perturbation of the vortex system away from the symmetric condition with asymmetric blowing.

Malcolm, Gerald N.↗

Development of non-conventional control methods for high angle of attack flight using vortex manipulation

Methods of manipulating the forebody vortices on a generic fighter model to produce controlled yawing moments at high angles of attack were investigated. Forces and moments were measured on the complete configuration and on the isolated forebody. Two schemes for vortex control on the forebody were evaluated: individually-controlled tip strakes and individually-controlled blowing ports. The effectiveness of the forebody strakes in controlling forebody side forces was strongly dependent on the size and location of the strakes. A yawing moment can be produced by deploying the forebody strakes asymmetrically, or can be eliminated by deploying the strakes symmetrically. The most effective strake position was found to be near the primary separation point, between 105 and 120 deg from windward. Blowing on the surface of the model was also shown to have a strong effect on the yawing moment. Blowing either forward or aft tangential to the surface appears to be more effective than blowing normal to the surface. The most effective method to control the yawing moment on the forebody was to minimize the natural asymmetry with a pair of small symmetrically mounted tip strakes and to perturb the vortex system away from the symmetric condition with asymmetric blowing.

Malcolm, Gerald N.↗

Instabilities of two-dimensional inviscid compressible vortices

An investigation of the linear stability and subsequent nonlinear evolution and acoustic radiation of a planar inviscid compressible vortex is presented. The effects of the entropy gradient are investigated, and for the particular entropy profile chosen, the positive average entropy in the vortex core is destabilizing, while the opposite is true for the negative average entropy gradient. Finite-difference methods are used to study the linear initial value problem. These methods are found to be capable of accurately computing the frequencies and weak growth rates of the normal modes. When the initial condition consists of random perturbations, the long-time behavior is found to correspond to the most unstable normal mode in all cases. The numerical solution of the Euler equations is used to study the nonlinear development of an elliptic-mode perturbation.

Chan, W. M.↗