Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximate boundary”

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 865 records · Page 48

Driven, steady-state RFP computations

The pseudospectral three-dimensional MHD code of Dahlburg et al. (1986 and 1987) is used to compute the dynamical behavior of a channel of magnetofluid carrying an axial current and magnetic flux. This situation contains the essential MHD behavior of the reversed-field pinch (RFP). An externally imposed electric field is applied to an initially current-free magnetofluid and drives currents that rise and eventually fluctuate about values corresponding to pinch ratios Theta of about 1.3, 2.2, and 4.5. A period of violent turbulence leads to an approximately force-free core, surrounded by an active MHD boundary layer that is not force-free. A steady state is reached that can apparently be sustained indefinitely (for several hundred Alfven transit times or longer). The turbulence level and time variability in the steady state increase with increasing Theta. The average toroidal magnetic field at the wall reverses for Theta = 2.2 and 4.5, but not for Theta = 1.3. Negative toroidal current filaments are observed. The Lundquist numbers are of the order of a few hundred.

Dahlburg, J. P.↗

Heat transfer with very high free-stream turbulence and streamwise vortices

Results are presented for two experimental programs related to augmentation of heat transfer by complex flow characteristics. In one program, high free stream turbulence (up to 63 percent) was shown to increase the Stanton number by more than a factor of 5, compared with the normally expected value based on x-Reynolds number. These experiments are being conducted in a free-jet facility, near the margins of the jet. To a limited extent, the mean velocity, turbulence intensity, and integral length scale can be separately varied. The results show that scale is a very important factor in determining the augmentation. Detailed studies of the turbulence structure are being carried out using an orthogonal triple hot-wire anemometer equipped with a fourth wire for measuring temperature. The v' component of turbulence appears to be distributed differently from u' or w'. In the second program, the velocity distributions and boundary layer thicknesses associated with a pair of counter-rotating, streamwise vortices were measured. There is a region of considerably thinned boundary layer between the two vortices when they are of approximately the same strength. If one vortex is much stronger than the other, the weaker vortex may be lifted off the surface and absorbed into the stronger.

Moffat, Robert J.↗

Boundary conditions and the simulation of low Mach number flows

The problem of accurately computing low Mach number flows, with the specific intent of studying the interaction of sound waves with incompressible flow structures, such as concentrations of vorticity is considered. This is a multiple time (and/or space) scales problem, leading to various difficulties in the design of numerical methods. Concentration is on one of these difficulties - the development of boundary conditions at artificial boundaries which allow sound waves and vortices to radiate to the far field. Nonlinear model equations are derived based on assumptions about the scaling of the variables. Then these are linearized about a uniform flow and exact boundary conditions are systematically derived using transform methods. Finally, useful approximations to the exact conditions which are valid for small Mach number and small viscosity are computed.

Hagstrom, Thomas↗

Navier-Stokes analysis of an oxidizer turbine blade with tip clearance with and without a mini-shroud

The Gas Generator Oxidizer Turbine (GGOT) Blade is being analyzed by various investigators under the NASA MSFC-sponsored Turbine Stage Technology Team design effort. The present work concentrates on the tip clearance region flow and associated losses; however, flow details for the passage region are also obtained in the simulations. The present calculations simulate the rotor blade row in a rotating reference frame with the appropriate coriolis and centrifugal acceleration term included in the momentum equations. The upstream computational boundary is located about one axial chord from the blade leading edge. The boundary conditions at this location have been determined by Pratt & Whitney using an Euler analysis without the vanes to obtain approximately the same flow profiles at the rotor as were obtained with the Euler stage analysis including the vanes. Inflow boundary layer profiles are then constructed assuming the skin friction coefficient at both the hub and the casing. The downstream computational boundary is located about one axial chord from the blade trailing edge, and the circumferentially averaged static pressure at this location was also obtained from the P&W Euler analysis. Results obtained for the 3-D baseline GGOT geometry at the full scale design Reynolds number show a region of high loss in the region near the casing. Particle traces in the near tip region show vortical flow behavior of the fluid which passes through the clearance region and exits at the downstream edge of the gap. In an effort to reduce clearance flow losses, the mini-shroud concept was proposed by the Pratt & Whitney design team. Calculations were performed on the GGO geometry with the mini-shroud. Results of these calculations indicate that the mini-shroud does not significantly affect the flow in the passage region, and although the tip clearance flow is different, the mini-shroud does not seem to prevent the above-mentioned vortical flow behavior. Since both flow distortion and total pressure losses are similar for both geometries, the addition of the mini-shroud does not seem to reduce the tip clearance flow effects.

Chan, Tony↗

Flight Measurements of Average Skin-Friction Coefficients on a Parabolic Body of Revolution (NACA RM-10) at Mach Numbers from 1.0 to 3.7

Measurement of average skin-friction coefficients have been made on six rocket-powered free-flight models by using the boundary-layer rake technique. The model configuration was the NACA RM-10, a 12.2-fineness-ratio parabolic body of revolution with a flat base. Measurements were made over a Mach number range from 1 to 3.7, a Reynolds number range 40 x 10(exp 6) to 170 x 10(exp 6) based on length to the measurement station, and with aerodynamic heating conditions varying from strong skin heating to strong skin cooling. The measurements show the same trends over the test ranges as Van Driest's theory for turbulent boundary layer on a flat plate. The measured values are approximately 7 percent higher than the values of the flat-plate theory. A comparison which takes into account the differences in Reynolds number is made between the present results and skin-friction measurements obtained on NACA RM-10 scale models in the Langley 4- by 4-foot supersonic pressure tunnel, the Lewis 8- by 6-foot supersonic tunnel, and the Langley 9-inch supersonic tunnel. Good agreement is shown at all but the lowest tunnel Reynolds number conditions. A simple empirical equation is developed which represents the measurements over the range of the tests.

Loposer, J. Dan↗

BASS Code Development

The BASS computational aeroacoustic code solves the fully nonlinear Euler equations in the time domain in two-dimensions. The acoustic response of the stator is determined simultaneously for the first three harmonics of the convected vortical gust of the rotor. The spatial mode generation, propagation and decay characteristics are predicted by assuming the acoustic field away from the stator can be represented as a uniform flow with small harmonic perturbations superimposed. The computed field is then decomposed using a joint temporal-spatial transform to determine the wave amplitudes as a function of rotor harmonic and spatial mode order. This report details the following technical aspects of the computations and analysis. 1) the BASS computational technique; 2) the application of periodic time shifted boundary conditions; 3) the linear theory aspects unique to rotor-stator interactions; and 4) the joint spatial-temporal transform. The computational results presented herein are twofold. In each case, the acoustic response of the stator is determined simultaneously for the first three harmonics of the convected vortical gust of the rotor. The fan under consideration here like modern fans is cut-off at +, and propagating acoustic waves are only expected at 2BPF and 3BPF. In the first case, the computations showed excellent agreement with linear theory predictions. The frequency and spatial mode order of acoustic field was computed and found consistent with linear theory. Further, the propagation of the generated modes was also correctly predicted. The upstream going waves propagated from the domain without reflection from the in ow boundary. However, reflections from the out ow boundary were noticed. The amplitude of the reflected wave was approximately 5% of the incident wave. The second set of computations were used to determine the influence of steady loading on the generated noise. Toward this end, the acoustic response was determined with three steady loading conditions: design, low-loading, high-loading. The overall trend showed significant (approximately 10 dB) increases in the generated noise for the highly loaded stator.

Sawyer, Scott↗

Modeling the Geologic History of Mt. Sharp

Gale is an approximately 155 km diameter crater located on the martian dichotomy boundary (5 deg S 138 deg E). Gale is estimated to have formed 3.8 - 3.5 Gya, in the late Noachian or early Hesperian. Mt. Sharp, at the center of Gale Crater, is a crescent shaped sedimentary mound that rises 5.2 km above the crater floor. Gale is one of the few craters that has a peak reaching higher than the rim of the crater wall. The Curiosity rover is currently fighting to find its way across a dune field at the northwest base of the mound searching for evidence of habitability. This study used orbital images and topographic data to refine models for the geologic history of Mt. Sharp by analyzing its morphological features. In addition, it assessed the possibility of a peak ring in Gale. The presence of a peak ring can offer important information to how Mt. Sharp was formed and eroded early in Gale's history.

Pascuzzo, A.↗

On the Convergence of Numerical Computations for Both Exact and Approximate Solutions for Electromagnetic Scattering by Nonspherical Dielectric Particles

We summarize the size parameter range of the applicability of four lightscattering computational methods for nonspheric dielectric particles. These methods include two exact methods - the extended boundary condition method (EBCM) and the invariant imbedding T-matrix method (II-TM) and two approximate approaches - the physical-geometric optics method (PGOM) and the improved geometric optics method (IGOM). For spheroids, the single-scattering properties computed by EBCM and II-TM agree for size parameters up to 150, and the comparison gives us confidence in using IITM as a benchmark for size parameters up to 150 for other geometries (e.g., hexagonal columns) because the applicability of II-TM with respect to particle shape is generic, as demonstrated in our previous studies involving a complex aggregate. This study demonstrates the convergence of the exact II-TM and approximate PGOM solutions for the complete set of single-scattering properties of a nonspherical shape other than spheroids and circular cylinders with particle sizes of ~48λ, specifically a hexagonal column with a size parameter of length as kL=300, where k=2π/λ and L is the column length. IGOM is also quite accurate except near the exact 180º backscattering direction. This study demonstrates that a synergetic combination of the numerically-exact II-TM and the approximate PGOM can seamlessly cover the entire size parameter range of practical interest. To demonstrate the applicability of the approach, we compute the optical properties of dust particles with a downstream application to the retrieval of dust aerosol optical thickness and effective particle size from polarimetric observations.

Ping Yang↗

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗

Investigation of Blade-row Flow Distributions in Axial-flow-compressor Stage Consisting of Guide Vanes and Rotor-blade Row

A 30-inch tip-diameter axial-flow compressor stage was investigated with and without rotor to determine individual blade-row performance, interblade-row effects, and outer-wall boundary-layer conditions. Velocity gradients at guide-vane outlet without rotor approximated design assumptions, when the measured variation of leaving angle was considered. With rotor in operation, Mach number and rotor-blade effects changed flow distribution leaving guide vanes and invalidated design assumption of radial equilibrium. Rotor-blade performance correlated interpolated two-dimensional results within 2 degrees, although tip stall was indicated in experimental and not two-dimensional results. Boundary-displacement thickness was less than 1.0 and 1.5 percent of passage height after guide vanes and after rotor, respectively, but increased rapidly after rotor when tip stall occurred.

COMPRESSORS - AXIAL FLOW↗

A Model for Turbulence Spectra in the Equilibrium Range of the Stable Atmospheric Boundary Layer

Stratification can cause turbulence spectra to deviate from Kolmogorov's isotropic -5/3 power law scaling in the universal equilibrium range at high Reynolds numbers. However, a consensus has not been reached with regard to the exact shape of the spectra. Here we propose a shape of the turbulent kinetic energy and temperature spectra in horizontal wavenumber for the equilibrium range that consists of three regimes at small Froude number: the buoyancy subrange, a transition region, and the isotropic inertial subrange through dimensional analysis and substantial revision of previous theoretical approximation. These spectral regimes are confirmed by various observations in the atmospheric boundary layer. The representation of the transition region in direct numerical simulations will require large-scale separation between the Dougherty-Ozmidov scale and the Kolmogorov scale for strongly stratified turbulence at high Reynolds numbers, which is still challenging computationally. In addition, we suggest that the failure of Monin-Obukhov similarity theory in the very stable atmospheric boundary layer is due to the fact that it does not consider the buoyancy scale that characterizes the transition region.

54 ENVIRONMENTAL SCIENCES↗

Boundaries of quantum supremacy via random circuit sampling

Abstract Google’s quantum supremacy experiment heralded a transition point where quantum computers can evaluate a computational task, random circuit sampling, faster than classical supercomputers. We examine the constraints on the region of quantum advantage for quantum circuits with a larger number of qubits and gates than experimentally implemented. At near-term gate fidelities, we demonstrate that quantum supremacy is limited to circuits with a qubit count and circuit depth of a few hundred. Larger circuits encounter two distinct boundaries: a return of a classical advantage and practically infeasible quantum runtimes. Decreasing error rates cause the region of a quantum advantage to grow rapidly. At error rates required for early implementations of the surface code, the largest circuit size within the quantum supremacy regime coincides approximately with the smallest circuit size needed to implement error correction. Thus, the boundaries of quantum supremacy may fortuitously coincide with the advent of scalable, error-corrected quantum computing.

97 MATHEMATICS AND COMPUTING↗

Numerical computations on one-dimensional inverse scattering problems

An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.

Dunn, M. H.↗

Numerical computations on one-dimensional inverse scattering problems

An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.

Dunn, M. H.↗

A Parametric Study of the Aerodynamic Characteristics of Nose-Cylinder-Flare Bodies at a Mach Number of 6.0

Force tests were conducted at a Mach number of 6.0 on nose-cylinder-flare bodies to determine the effect of nose shape, cylinder length, flare angle, and flare length on the longitudinal aerodynamic characteristics. A particular investigation was conducted to determine the effect of flare angle for constant flare length, surface area, and diameter. Results indicated that at a Reynolds number of approximately 0.92 x l0 (exp 6) (based on body diameter), the boundary-layer separation effects were significant only with respect to the slope of the normal-force and pitching-moment curve at low angles of attack. The variations of the aerodynamic characteristics with the various parameters were, in general, similar to those predicted by Newtonian theory below a flare angle of 30 degrees and a ratio of flare base diameter to cylinder diameter of less than approximately 2.2. The limiting diameter ratio is consistent with the extent of the low-constant dynamic-pressure region near the body caused by the bow-shock influences as predicted by axisymmetric characteristic theory. The effects of the various parameters for the flares that exceeded the limiting diameter ratio follow the trends predicted by the computed flow-field properties. The axial force for these flare configurations at zero angle of attack was, in general, computed within 10 percent by using these properties. For a constant flare length and surface area the flare effectiveness increased with increasing flare angle; however, for constant flare diameter only the axial-force coefficient was affected by flare angle.

Ashby, George C., Jr.↗

Approximate solutions to the Navier-Stokes initial value problem

A Galerkin-Ritz procedure for any arbitrary system of field equations is shown to follow generically from 'two-functional' variational conditions. The initial-value problem for boundary-free incompressible Navier-Stokes flow is solved analytically in the two-parameter-function approximation.

Rosen, G.↗

The interaction between a strong longitudinal vortex and a turbulent boundary layer

This paper presents the preliminary results of a long-term investigation of the interaction of strong longitudinal vortices with turbulent boundary layers. A pair of trailing vortices is generated by a delta wing, which is mounted ahead of a flat plate so that the trailing vortices merge with the boundary layer on the upper surface of the flat plate, while the non-rolled-up part of the delta-wing wake passes below the plate. This approximates an aircraft wing with a close-coupled canard. The measurements show how the boundary-layer fluid is first subjected to large cross flow and is then transported away from the plate to circulate around the vortex core.

Cutler, A. D.↗

Numerical Study of Pressure Fluctuations due to High-Speed Turbulent Boundary Layers

Direct numerical simulations (DNS) are used to examine the pressure fluctuations generated by fully developed turbulence in supersonic turbulent boundary layers with an emphasis on both pressure fluctuations at the wall and the acoustic fluctuations radiated into the freestream. The wall and freestream pressure fields are first analyzed for a zero pressure gradient boundary layer with Mach 2.5 and Reynolds number based on momentum thickness of approximately 2835. The single and multi-point statistics reported include the wall pressure fluctuation intensities, frequency spectra, space-time correlations, and convection velocities. Single and multi-point statistics of surface pressure fluctuations show good agreement with measured data and previously published simulations of turbulent boundary layers under similar flow conditions. Spectral analysis shows that the acoustic fluctuations outside the boundary layer region have much lower energy content within the high-frequency region. The space-time correlations reflect the convective nature of the pressure field both at the wall and in the freestream, which is characterized by the downstream propagation of pressure-carrying eddies. Relative to those at the wall, the pressure-carrying eddies associated with the freestream signal are larger and convect at a significantly lower speed. The preliminary DNS results of a Mach 6 boundary layer show that the pressure rms in the freestream region is significantly higher than that of the lower Mach number case.

Duan, Lian↗