Search NASA⌕ Search

SEARCH · Search NASA

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Turbulent Prandtl number in the near-wall region of a turbulent channel flow

Calculations are presented of the turbulent Prandtl number Pr(T) in the near-wall region of a turbulent channel flow. It is shown that only the first-order terms in the Taylor series expansions for the eddy diffusivity and the Pr(T) are independent of the molecular Prandtl number Pr, at least to a first approximation. Also presented are calculations of the near-wall behavior of the correlations between temperature fluctuations and velocity fluctuations, as well as for their dependence on Pr.

Antonia, R. A.↗
NASA NTRS2002

Measurements of Turbulence Attenuation by a Dilute Dispersion of Solid Particles in Homogeneous Isotropic Turbulence

This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Eaton, John↗
NASA NTRS2007

Further Characterization of an Active Clearance Control Concept

A new test chamber and precision hydraulic actuation system were incorporated into an active clearance control (ACC) test rig at NASA Glenn Research Center. Using the improved system, a fast-acting, mechanically-actuated, ACC concept was evaluated at engine simulated temperatures and pressure differentials up to 1140 F and 120 psig, on the basis of secondary seal leakage and kinematic controllability. During testing, the ACC concept tracked a simulated flight clearance transient profile at 1140 F, 120 psig, with a maximum error of only 0.0012 in. Comparison of average dynamic leakage of the system with average static leakage did not show significant differences between the two operating conditions. Calculated effective clearance values for the rig were approximately 0.0002 in. at 120 psig, well below the industry specified effective clearance threshold of 0.001 in.

Taylor, Shawn C.↗
DOE OSTI2022

Fresnel zone plate point spread function approximation for zeroth order mitigation in millimetric field of view x-ray imaging

High spatial and temporal resolution x-ray radiography images are required at the National Ignition Facility (NIF) for high-energy density experiments. One technique that is in development to achieve the required resolution uses Fresnel zone plate (FZP) optics to image an object that is backlit by an x-ray source. The multiple FZP diffraction orders do not focus on the same plane, which increases the background and reduces the contrast. Understanding the point spread function of the different diffraction orders will allow the prediction of the expected background using simulations. We find that the two-dimensional point spread function of the FZP can be approximated by the addition of a sharp Gaussian with a disk. This allowed for the estimation of the background in NIF experimental images of Rayleigh–Taylor spikes and their interpretation. An alternative design of FZP is discussed to allow the inclusion of a zeroth order blocker to reduce the background.

Do, A. (ORCID:0000000174372147)↗
DOE OSTI2018

Subpercent-Scale Control of 3D Low Modes of Targets Imploded in Direct-Drive Configuration on OMEGA

In a series of direct-drive implosions on OMEGA, multiple time resolved x-ray images were used to tomographically measure their 3-D modes 1, 2, and 3 at a convergence ratio of ~3. Results show that the target modes vary linearly with the laser modes and are not affected by the Rayleigh–Taylor growth or lateral heat transport. This indicates that the residual modes (resulting from physical effects including beam mistiming, mispointing, and laser energy calibration) are approximately constant between shots. Lastly, this demonstrates that the low-mode amplitudes can be mitigated within by adjusting the laser-energy balance to compensate the residual target modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

A Fixed-Point Iteration Method with Quadratic Convergence

The fixed-point iteration algorithm is turned into a quadratically convergent scheme for a system of nonlinear equations. Most of the usual methods for obtaining the roots of a system of nonlinear equations rely on expanding the equation system about the roots in a Taylor series, and neglecting the higher order terms. Rearrangement of the resulting truncated system then results in the usual Newton-Raphson and Halley type approximations. In this paper the introduction of unit root functions avoids the direct expansion of the nonlinear system about the root, and relies, instead, on approximations which enable the unit root functions to considerably widen the radius of convergence of the iteration method. Methods for obtaining higher order rates of convergence and larger radii of convergence are discussed.

Walker, Kevin P.↗
NASA NTRS1979

The instability of a horizontal magnetic field in an atmosphere stable against convection

The theoretical problem posed by the buoyant escape of a magnetic field from the interior of a stably stratified body bears directly on the question of the present existence of primordial magnetic fields in stars. This paper treats the onset of the Rayleigh-Taylor instability of the upper boundary of a uniform horizontal magnetic field in a stably stratified atmosphere. The calculations are carried out in the Boussinesq approximation and show the rapid growth of the initial infinitesimal perturbation of the boundary. This result is in contrast to the extremely slow buoyant rise of a separate flux tube in the same atmosphere. Thus for instance, at a depth of 1/3 of a solar radius beneath the surface of the sun, a field of 100 G develops ripples over a scale of 1000 km in a characteristic time of 50 years, whereas the characteristic rise time of the same field in separate flux tubes with the same dimensions is 10 billion years. Thus, the development of irregularities proceeds quickly, soon slowing, however, to a very slow pace when the amplitude of the irregularities becomes significant. Altogether, the calculations show the complexity of the question of the existence of remnant primordial magnetic fields in stellar interiors.

Parker, E. N.↗
NASA NTRS1987

A full CI treatment of the 1A1, 1B1, and 3B1 states of SiH2

Full CI calculations are presented for the 1A1, 3B1, and 1B1 states of SiH2 at their respective equilibrium geometries and at geometries with the SiH bonds stretched. These results are compared with those obtained from single-reference and multireference CI calculations. The computed Te values agree well with the full CI results, provided that the effects of higher-than-double excitations are accounted for either by the Davidson correction or by a multireference approach. When the SiH bonds are stretched, the single-reference methods are not sufficiently flexible, and only CASSCF/MRCI achieves chemical accuracy (i.e., agrees with the full CI to 1 kcal/mol). Overall, the accuracy of the various approximate methods is very similar to that found for H2O, NH2, and CH2.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗
NASA NTRS1995

Interstellar clouds in high-speed, supersonic flows: Two-dimensional simulations

We present a series of gasdynamical simulations of the interaction of a dense, cool interstellar cloud with a high-speed, supersonic wind that confines and accelerates the embedded cloud. Our goal is to attempt to determine if such clouds can survive various potentially disruptive instabilities, that occur at their peripheries, long enough to be accelerated to speeds which are comparable to the wind velocity. These simulations are performed using two-dimensional, Eulerian gas dynamics on both an axisymmetric (about the cloud axis) and 'slab' geometric grid. The spatial and temporal resolutions of the simulations are varied over a wide range to investigate the effects of small-scale instabilities on the overall acceleration of clouds and the development of large-scale, disruptive instabilities. Also, we study the effects of wind/cloud Mach number variations by changing the wind speed constant at about 12 km/s (which corresponds to a cloud temperature of 10,000 K). The current simulations track the evolution of clouds as they are accelerated to speeds approximately 4-5 times greater than their internal sound speeds. Furthermore, the models with the highest resolution were extended far beyond quasi-linear Rayleigh-Taylor growth times reaching 6-7 Rayleigh-Taylor growth times for the largest scale instabilities before being terminated because of the accumulation of errors at the rear grid boundary.

Schiano, A. V. R.↗
DOE OSTI2019

Parametric investigation of the transition to turbulence in Rayleigh-Taylor mixing

In this study, a database is developed consisting of highly resolved large-eddy simulations of the transition process in a Rayleigh-Taylor mixing layer. Over 300 simulations are performed with parametric variations in Grashof number, Schmidt number, Atwood number, and initial perturbation spectrum. The database is then analyzed to explore how the factors that influence the transition-to-turbulence process impact mixing evolution and growth rate prior to the development of late-time, steady-state growth. The competition between entrainment-driven mixing and diffusion-driven mixing is observed, and two different mechanisms are suggested by which the homogeneity of mixing might increase, rather than decrease, at early time. The impact of parametric variation on growth rate is additionally examined, and all variations in initial conditions are found to converge to approximately the same growth parameter around α ≈ 0.025.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Roles of bulk viscosity on Rayleigh-Taylor instability: Non-equilibrium thermodynamics due to spatio-temporal pressure fronts

Direct numerical simulations of Rayleigh-Taylor instability (RTI) between two air masses with a temperature difference of 70 K is presented using compressible Navier-Stokes formulation in a non-equilibrium thermodynamic framework. The two-dimensional flow is studied in an isolated box with non-periodic walls in both vertical and horizontal directions. The non-conducting interface separating the two air masses is impulsively removed at t = 0 (depicting a heaviside function). No external perturbation has been used at the interface to instigate the instability at the onset. Computations have been carried out for rectangular and square cross sections. The formulation is free of Boussinesq approximation commonly used in many Navier-Stokes formulations for RTI. Effect of Stokes’ hypothesis is quantified, by using models from acoustic attenuation measurement for the second coefficient of viscosity from two experiments. Effects of Stokes’ hypothesis on growth of mixing layer and evolution of total entropy for the Rayleigh-Taylor system are reported. The initial rate of growth is observed to be independent of Stokes’ hypothesis and the geometry of the box. Following this stage, growth rate is dependent on the geometry of the box and is sensitive to the model used. As a consequence of compressible formulation, we capture pressure wave-packets with associated reflection and rarefaction from the non-periodic walls. The pattern and frequency of reflections of pressure waves noted specifically at the initial stages are reflected in entropy variation of the system.

42 ENGINEERING↗
NASA NTRS1975

On the analytic structure of the Taylor-Maccoll conical-flow solution

The conical flow solution for axisymmetric supersonic flow past cones has been found to be virtually independent of the ratio of specific heats when normalized in a certain way. A simple rational approximation to this flow is derived. The important singularities and the limiting behavior of the solution are also discussed.

Schwartz, L. W.↗
DOE OSTI2017

Optimal Bandwidth for Multitaper Spectrum Estimation

A systematic method for bandwidth parameter selection is desired for Thomson multitaper spectrum estimation. We give a method for determining the optimal bandwidth based on a mean squared error (MSE) criterion. When the true spectrum has a second-order Taylor series expansion, one can express quadratic local bias as a function of the curvature of the spectrum, which can be estimated by using a simple spline approximation. This is combined with a variance estimate, obtained by jackknifing over individual spectrum estimates, to produce an estimated MSE for the log spectrum estimate for each choice of time-bandwidth product. The bandwidth that minimizes the estimated MSE then gives the desired spectrum estimate. Additionally, the bandwidth obtained using our method is also optimal for cepstrum estimates. We give an example of a damped oscillatory (Lorentzian) process in which the approximate optimal bandwidth can be written as a function of the damping parameter. Furthermore, the true optimal bandwidth agrees well with that given by minimizing estimated the MSE in these examples.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

Some self starting integrators for x Prime equals f (x, t)

The integration is discussed of the vector differential equation X = F(x, t) from time t sub i to t sub (i = 1) where only the values of x sub i are available for the the integration. No previous values of x or x prime are used. Using an orbit integration problem, comparisons are made between Taylor series integrators and various types and orders of Runge-Kutta integrators. A fourth order Runge-Kutta type integrator for orbital work is presented, and approximate (there may be no exact) fifth order Runge-Kutta integrators are discussed. Also discussed and compared is a self starting integrator ising delta f/delta x. A numerical method for controlling the accuracy of integration is given, and the special equations for accurately integrating accelerometer data are shown.

Lear, W. M.↗
NASA NTRS1980

An atlas of 1977 and 1978 GEOS-3 radar altimeter data for tropical cyclone studies

All of the GEOS 3 satellite altimeter schedule information were collected with all of the available 1977 and 1978 tropical cyclone positional information. The time period covers from March 23, 1977 through Nov. 23, 1978. The geographical region includes all ocean area north of the equator divided into the following operational areas: the Atlantic area (which includes the Caribbean and Gulf of Mexico); the eastern Pacific area; the central and western Pacific area; and the Indian Ocean area. All available source material concerning tropical cyclones was collected. The date/time/location information was extracted for each disturbance. This information was compared with the GEOS 3 altimeter ON/OFF history information to determine the existence of any altimeter data close enough in both time and location to make the data potentially useful for further study (the very liberal criteria used was time less than 24 hours and location within 25 degrees). Geographic plots (cyclone versus GEOS 3 orbit track) were produced for all of the events found showing the approximate location of the cyclone and the GEOS 3 orbit traces for the full day.

Stanley, H. R.↗
NASA NTRS2013

Nonlinear Modeling by Assembling Piecewise Linear Models

To preserve nonlinearity of a full order system over a parameters range of interest, we propose a simple modeling approach by assembling a set of piecewise local solutions, including the first-order Taylor series terms expanded about some sampling states. The work by Rewienski and White inspired our use of piecewise linear local solutions. The assembly of these local approximations is accomplished by assigning nonlinear weights, through radial basis functions in this study. The efficacy of the proposed procedure is validated for a two-dimensional airfoil moving at different Mach numbers and pitching motions, under which the flow exhibits prominent nonlinear behaviors. All results confirm that our nonlinear model is accurate and stable for predicting not only aerodynamic forces but also detailed flowfields. Moreover, the model is robustness-accurate for inputs considerably different from the base trajectory in form and magnitude. This modeling preserves nonlinearity of the problems considered in a rather simple and accurate manner.

Yao, Weigang↗