Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Taylor approximation”

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 145 records · Page 8

Decoupling the Role of Particle Inertia and Gravity on Particle Dispersion

Particle dispersion and the influence that particle momentum exchange has on the properties of a turbulent carrier flow in micro-gravity environments challenge present understanding and predictive schemes. The objective of this effort has been to develop and assess high-fidelity simulation tools for predicting particle transport within micro-gravity environments suspended in turbulent flows. The computational technique is based on Direct Numerical Simulation (DNS) of the incompressible Navier-Stokes equations. The particular focus of the present work is on the class of dilute flows in which particle volume fractions and inter-particle collisions are negligible. Particle motion is assumed to be governed by drag with particle relaxation times ranging from the Kolmogorov scale to the Eulerian timescale of the turbulence and particle mass loadings up to one. The velocity field was made statistically stationary by forcing the low wavenumbers of the flow. The calculations were performed using 96(exp 3) collocation points and the Taylor-scale Reynolds number for the stationary flow was 62. The effect of particles on the turbulence was included in the Navier-Stokes equations using the point-force approximation in which 96(exp 3) particles were used in the calculations. DNS results show that particles increasingly dissipate fluid kinetic energy with increased loading, with the reduction in kinetic energy being relatively independent of the particle relaxation time. Viscous dissipation in the fluid decreases with increased loading and is larger for particles with smaller relaxation times. Fluid energy spectra show that there is a non-uniform distortion of the turbulence with a relative increase in small-scale energy. The non-uniform distortion significantly affects the transport of the dissipation rate, with the production and destruction of dissipation exhibiting completely different behaviors. The spectrum of the fluid-particle energy exchange rate shows that the fluid drags particles at low wavenumbers while the converse is true at high wavenumbers for small particles. A spectral analysis shows that the increase of the high wavenumber portion of the fluid energy spectrum can be attributed to transfer of the fluid-particle covariance by the fluid turbulence. This in turn explains the relative increase of small-scale energy caused by small particles observed in the present simulations as well as those of others.

Squires, Kyle D.↗

A Method for Calculating Viscosity and Thermal Conductivity of a Helium-Xenon Gas Mixture

A method for calculating viscosity and thermal conductivity of a helium-xenon (He-Xe) gas mixture was employed, and results were compared to AiResearch (part of Honeywell) analytical data. The method of choice was that presented by Hirschfelder with Singh's third-order correction factor applied to thermal conductivity. Values for viscosity and thermal conductivity were calculated over a temperature range of 400 to 1200 K for He-Xe gas mixture molecular weights of 20.183, 39.94, and 83.8 kg/kmol. First-order values for both transport properties were in good agreement with AiResearch analytical data. Third-order-corrected thermal conductivity values were all greater than AiResearch data, but were considered to be a better approximation of thermal conductivity because higher-order effects of mass and temperature were taken into consideration. Viscosity, conductivity, and Prandtl number were then compared to experimental data presented by Taylor.

Johnson, Paul K.↗

Mathematical models of the open magnetosphere - Application to dayside auroras.

Two static mathematical models of the open or Dungey model of the magnetosphere are constructed. The process of construction is similar to that for early closed magnetosphere models, such as the Taylor-Hones model. The first model in fact is simply an addition of an interplanetary field in arbitrary direction to a Taylor-Hones image dipole model. In order to preserve the shape of the magnetosphere at high latitudes, and to partially exclude the exterior field, another model is constructed with the magnetopause approximated by a diamagnetic sphere. We find that there are some interplanetary field lines connected to the earth for all orientations of the interplanetary field other than strictly northward, and that the maximum number of connected field lines occurs with a due southward field. For an average spiral hose angle of the interplanetary field, the dayside neutral point occurs on the magnetopause at about 10 o'clock local time. Dayside auroras, convection patterns, and other phenomena may exhibit symmetry about this local time. For a positive (negative) interplanetary field sector, energetic, anisotropic particle fluxes should have direct access to the northern (southern) polar caps, as is supported by many recent observations.

Forbes, T. G.↗

The Dynamics of Miscible Interfaces: A Space Flight Experiment

Experiments as well as accompanying simulations are described that serve in preparation of a space flight experiment to study the dynamics of miscible interfaces. The investigation specifically addresses the importance of both nonsolenoidal effects as well as nonconventional Korteweg stresses in flows that give rise to steep but finite concentration gradients. The investigation focuses on the flow in which a less viscous fluid displaces one of higher viscosity and different density within a narrow capillary tube. The fluids are miscible in all proportions. An intruding finger forms that occupies a fraction of the total tube diameter. Depending on the flow conditions, as expressed by the Peclet number, a dimensionless viscosity ratio, and a gravity parameter, this fraction can vary between approximately 0.9 and 0.2. For large Pe values, a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the specific flow conditions, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip. For small Pe values, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). Comparisons between experiments and numerical simulations based on the 'conventional' assumption of solenoidal velocity fields and without Korteweg stresses yield poor agreement as far as the Pe value is concerned that distinguishes these two regimes. As one possibility, we attribute this lack of agreement to the disregard of these terms. An attempt is made to use scaling arguments in order to evaluate the importance of the Korteweg stresses and of the assumption of solenoidality. While these effects should be strongest in absolute terms when steep concentration fronts exist, i.e., at large Pe, they may be relatively most important at lower values of Pe. We subsequently compare these conventional simulations to more complete simulations that account for nonvanishing divergence as well as Korteweg stresses. While the exact value of the relevant stress coefficients are not known, ballpark numbers do exist, and their use in the simulations indicates that these stresses may indeed be important. We plan to evaluate these issues in detail by means of comparing a space experiment with corresponding simulations, in order to extract more accurate Korteweg stress coefficients, and to confirm or deny the importance of such stresses.

Maxworthy, Tony↗

Methods and approximations for the computation of transmission profiles in the v4 band of methane in the atmosphere of Jupiter.

This note discusses the validity of certain band models and scaling approximations for computing transmissions in the v4 band of methane along inhomogeneous paths in the atmosphere of Jupiter. It is shown that Goody's random band model approximates the results of a rigorous numerical line-by-line calculation of the transmission profile of a Jovian model atmosphere.

Taylor, F. W.↗

Comparison of predicted and measured drag for a single-engine airplane

Renewed interest in natural laminar flow (NLF) has rekindled designers' concerns that manufacturing deviations, (loss of surface contours or other surface imperfections) may destroy the effectiveness of NLF for an operational airplane. This paper reports on experimental research that compares predicted and measured boundary layer transition, total drag, and two-dimensional drag coefficients for three different wing surface conditions on an airplane typical of general aviation manufacturing technology. The three flight test phases included: (1) assessment of an unpainted airframe, (2) flight tests of the same airplane after painstakingly filling and sanding the wings to design contours, and (3) similar measurements after this airplane was painted. In each flight phase, transition locations were monitored using either sublimating chemicals or pigmented oil. As expected, total drag changes were difficult to measure. Two-dimensional drag coefficients were estimated using the Eppler-Somers code and measured with a wake rake in a method very similar to Jones' pitot traverse method. The net change in two-dimensional drag was approximately 20 counts between the unpainted airplane and the 'hand-smoothed' airplane for typical cruise flight conditions.

Ward, D. T.↗

Flight Controllability Limits and Related Human Transfer Functions as Determined from Simulator and Flight Tests

A simulator study and flight tests were performed to determine the levels of static stability and damping necessary to enable a pilot to control the longitudinal and lateral-directional dynamics of a vehicle for short periods of time. Although a basic set of aerodynamic characteristics was used, the study was conducted so that the results would be applicable to a wide range of flight conditions and configurations. Novel piloting techniques were found which enabled the pilot to control the vehicle at conditions that were otherwise uncontrollable. The influence of several critical factors in altering the controllability limits was also investigated. Several human transfer functions were used which gave fairly good representations of the controllability limits determined experimentally for the short-period longitudinal, directional, and lateral modes. A transfer function with approximately the same gain and phase angle as the pilot at the controlling frequencies along the controllability limits was also derived.

Taylor, Lawrence W., Jr.↗

X-ray digital imaging petrography of lunar mare soils: modal analyses of minerals and glasses

It is essential that accurate modal (i.e., volume) percentages of the various mineral and glass phases in lunar soils be used for addressing and resolving the effects of space weathering upon reflectance spectra, as well as for their calibration such data are also required for evaluating the resource potential of lunar minerals for use at a lunar base. However, these data are largely lacking. Particle-counting information for lunar soils, originally obtained to study formational processes, does not provide these necessary data, including the percentages of minerals locked in multi-phase lithic fragments and fused-soil particles, such as agglutinates. We have developed a technique for modal analyses, sensu stricto, of lunar soils, using digital imaging of X-ray maps obtained with an energy-dispersive spectrometer mounted on an electron microprobe. A suite of nine soils (90 to 150 micrometers size fraction) from the Apollo 11, 12, 15, and 17 mare sites was used for this study. This is the first collection of such modal data on soils from all Apollo mare sites. The abundances of free-mineral fragments in the mare soils are greater for immature and submature soils than for mature soils, largely because of the formation of agglutinitic glass as maturity progresses. In considerations of resource utilization at a lunar base, the best lunar soils to use for mineral beneficiation (i.e., most free-mineral fragments) have maturities near the immature/submature boundary (Is/FeO approximately or = 30), not the mature soils with their complications due to extensive agglutination. The particle data obtained from the nine mare soils confirm the generalizations for lunar soils predicted by L.A. Taylor and D.S. McKay (1992, Lunar Planet Sci. Conf. 23rd, pp. 1411-1412 [Abstract]).

NASA Discipline Exobiology↗

Augmenting subspace optimization methods with linear bandits

In this work, we consider the framework of methods for unconstrained minimization that are, in each iteration, restricted to a model that is only a valid approximation to the objective function on some affine subspace containing an incumbent point. These methods are of practical interest in computational settings where derivative information is either expensive or impossible to obtain. Recent attention has been paid in the literature to employing randomized matrix sketching for generating the affine subspaces within this framework. We consider a relatively straightforward, deterministic augmentation of such a generic subspace optimization method. In particular, we consider a sequential optimization framework where actions consist of one-dimensional linear subspaces and rewards consist of (approximations to) the magnitudes of directional derivatives computed in the direction of the action subspace. Reward maximization in this context is consistent with maximizing lower bounds on descent guaranteed by first-order Taylor models. This sequential optimization problem can be analysed through the lens of dynamic regret. We modify an existing linear upper confidence bound (UCB) bandit method and prove sublinear dynamic regret in the subspace optimization setting. We demonstrate the efficacy of employing this linear UCB method in a setting where forward-mode algorithmic differentiation can provide directional derivatives in arbitrary directions and in a derivative-free setting. For the derivative-free setting, we propose SS-POUNDers, an extension of the derivative-free optimization method POUNDers that employs the linear UCB mechanism to identify promising subspaces. Our numerical experiments suggest a preference, in either computational setting, for employing a linear UCB mechanism within a subspace optimization method.

97 MATHEMATICS AND COMPUTING↗

On incorporating damping and gravity effects in models of structural dynamics of the SCOLE configuration

The damping for structural dynamic models of flexible spacecraft is usually ignored and then added after modal frequencies and mode shapes are calculated. It is common practice to assume the same damping ratio for all modes, although it is known that damping due to bending and that due to torsion are sometimes ignored. Two methods of including damping in the modeling process from its onset are examined. First, the partial derivative equations of motion are analyzed for a pinned-pinned beam with damping. The end conditions are altered to handle bodies with mass and inertia for the Spacecraft Control Laboratory Experiment (SCOLE) configuration. Second, a massless beam approximation is used for the modes with low frequencies, and a clamped-clamped system is used to approximate the modes for arbitrarily high frequency. The model is then modified to include gravity effects and is compared with experimental results.

Taylor, Larry↗

Derivation of generalized transition/boundary conditions for planar multiple-layer structures

Infinite-order generalized impedance boundary conditions (GIBCs) and generalized sheet transition conditions (GSTCs) for planar multilayer configurations are developed via the Taylor series expansion method. The conditions are derived in a matrix product form where each matrix corresponds to a specific layer. An overall composite boundary/transition condition is obtained by making finite-order approximations to the elements of each matrix for the cases of 'low-contrast' and 'high-contrast' material layers. The accuracy of the truncated boundary conditions is examined by comparing their implied reflection and transmission coefficients with the corresponding exact coefficients. Design curves are also given which relate the maximum order of the conditions required to simulate a coating or layer of specific thickness and contrast. Expressions are then derived for the reflection and transmission coefficients of the GIBC/GSTC sheets, and these are compared to exact coefficients to demonstrate the validity of the derived GIBCs/GSTCs.

Ricoy, M. A.↗

Probing Venus's cloud structure with Galileo NIMS

The spectral image cubes obtained by the Near-Infrared Mapping Spectrometer (NIMS) on Galileo as it flew by Venus have been analyzed to constrain the vertical structure of the clouds, the nature of the aerosol particles, and the location and particle properties of the opacity variations responsible for high-contrast features observed in the near-infrared windows at 1.7 and 2.3 micrometers. A radiative transfer program was used to simulate mid-latitude curves of limb darkening at 3.7 micrometers. Best-fit models to these curves demonstrate that the upper clouds are dominated by mode 2 particles (r-bar = 1.0 micrometers), with a contribution of approximately 15% of opacity from mode 1 particles (r-bar = 0.3 micrometers). The low-latitude upper cloud is well represented by a dual scale-height model, with a particle scale height of approximately 1 km from an altitude of 61-63 km, and a scale height of approximately 6 km above this, up to the level where tau = 1 at approximately 71 km. This model also successfully simulates limb-darkening curves at 11.5 micrometers from the Pioneer Venus Orbiter Infrared Radiometer. Successful simulations of correlation plots of 1.7 vs 2.3 micrometers intensities reveal that mode 3 particles (r-bar = 3.65 micrometers) represent the dominant source of opacity in the lower and middle clouds, and that variation in total cloud opacity reflects chiefly the addition and removal of mode 3 particles near the cloud base. We find that the full spectrum of brightnesses at 1.7 and 2.3 micrometers implies that the total cloud optical depth varies from approximately 25 to approximately 40.

Grinspoon, D. H.↗

Paleo-Pole Positions from Martian Magnetic Anomaly Data

Magnetic component anomaly maps were made from five mapping cycles of the Mars Global Surveyor s magnetometer data. Our goal was to find and isolate positive and negative anomaly pairs which would indicate magnetization of a single source body. From these anomalies we could compute the direction of the magnetizing vector and subsequently the location of the magnetic pole existing at the time of magnetization. We found nine suitable anomaly pairs and from these we computed four North and 3 South poles with two at approximately 60 degrees north latitude. These results suggest that during the existence of the Martian main magnetic field it experienced several reversals.

Taylor, Patrick T.↗

Approach for Input Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

An implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for quasi 3-D Euler CFD code is presented. Given uncertainties in statistically independent, random, normally distributed input variables, first- and second-order statistical moment procedures are performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, these moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

Polynominal Interpolation Methods for Viscous Flow Calculations

Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and by Hermitian (Taylor series) finite-difference discretization. The similarities and special features of these different developments are discussed. The governing systems apply for both uniform and variable meshes. Hybrid schemes resulting from two different polynomial approximations for the first and second derivatives lead to a nonuniform mesh extension of the so-called compact or Pad? difference technique (Hermite 4). A variety of fourth-order methods are described and the Hermitian approach is extended to sixth-order (Hermite 6). The appropriate spline boundary conditions are derived for all procedures. For central finite differences, this leads to a two-point, second-order accurate generalization of the commonly used three-point end-difference formula. Solutions with several spline and Hermite procedures are presented for the boundary layer equations, with and without mass transfer, and for the incompressible viscous flow in a driven cavity. Divergence and nondivergence equations are considered for the cavity. Among the fourth-order techniques, it is shown that spline 4 has the smallest truncation error. The spline 4 procedure generally requires one-quarter the number of mesh points in a given coordinate direction as a central finite-difference calculation of equal accuracy. The Hermite 6 procedure leads to remarkably accurate boundary layer solutions.

Rubin, S. G.↗

Solution of the linear shallow water equations by the fourth-order leapfrog scheme

Numerical schemes of the first and second order of approximation introduce numerical distortion when the wave propagation over a long distance is investigated. To alleviate this problem, the fourth-order leapfrog scheme is constructed. The standard leapfrog method is based on the truncated Taylor series expansion which depicts an error proportional to the second-order terms. In the proposed method the numerical solution is corrected for these terms. The space and time corrections work well in diminishing numerical dispersion and dissipation.

Kowalik, Z.↗

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗

Formation of a rotational accretion column

We consider the problem of the flow produced by the accretion of gas from a rotating cloud onto a point gravitational source. It is shown that, for a cloud initially in uniform rotation, the effect of accretion is to produce a Taylor column in the cloud parallel to the rotation axis. Only gas within the column is accreted; the gas outside the column undergoes oscillatory flow, which is probably subject to shear instability. The column forms in approximately one-half rotation period from the time at which accretion begins. Its initial radius is determined by the accretion rate and the angular velocity of the cloud. Application of the results to the formation of Jupiter by accretion is discussed.

Cassen, P.↗