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 757 records · Page 42

Interplanetary stream magnetism: Kinematic effects

The particle density, and the magnetic field intensity and direction are calculated in corotating streams of the solar wind, assuming that the solar wind velocity is constant and radial and that its azimuthal variations are not two rapid. The effects of the radial velocity profile in corotating streams on the magnetic fields were examined using kinematic approximation and a variety of field configurations on the inner boundary. Kinematic and dynamic effects are discussed.

Burlaga, L. F.↗

Three-dimensional elastic stress and displacement analysis of finite geometry solids containing cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement distributions in various bodies containing cracks. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The results obtained show a considerable potential for using this method in the three-dimensional analysis of finite geometry solids and suggest a possible extension of this technique to nonlinear material behavior.

Gyekenyesi, J. P.↗

The response of the dayside aurora to sharp northward and southward transitions of the interplanetary magnetic field and to magnetospheric substorms

Latitudinal shifts of the discrete dayside aurora are examined with respect to variations in the interplanetary magnetic field (IMF) and to magnetospheric substorms. Within 10-15 min after steplike southward (northward) transitions in the IMF the dayside auroral oval moves equatorward (poleward). In at least some cases the auroral shift is similar to an exponential relaxation in latitude from an initial to a final steady state value; for these cases the average exponential time constant is estimated to be 17 min. Substorm features in the dayside aurora include (1) an equatorward shift by 1-3 deg of the equatorward discrete auroral boundary, (2) a brightening of the discrete aurora near substorm onset, (3) the formation of multiple auroral bands, and (4) poleward motion of short-lived individual auroral arcs occurring nearly coincident in time with the equatorward boundary shifts, the average poleward velocity of these arcs being approximately 800 m/s.

Horwitz, J. L.↗

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗

Application of singular perturbation techniques /SPT/ and continuation methods for on-line aircraft trajectory optimization

The singular perturbation theory (SPT) approximation of optimal feedback control laws is presented and methods for on-line application of these approximations are discussed. It is demonstrated that SPT control laws break down when the current state is near the terminal target state. The use of continuation methods to improve the accuracy of the SPT approximation and to obtain global solutions of two-point boundary value problems is also discussed. As an illustration, consideration is given to the minimum-time control of a supersonic aircraft for a three-dimensional intercept problem.

Washburn, R. B.↗

A kinetic cross-field streaming instability

In a high-beta plasma the so-called modified-two-stream instability, which results from strongly magnetized electrons drifting relative to unmagnetized ions across a homogeneous magnetic field, is misnamed because the mode is highly kinetic, particularly when the relative streaming velocity exceeds the Alfven speed of the plasma. This kinetic cross-field streaming instability is investigated in detail, examining the effect of the electromagnetic terms and the stability boundaries in both low- and high-beta plasmas. An approximate dispersion relation showing the relation of this mode to the whistler is derived and solutions of it are compared with those obtained from the exact dispersion relation. The kinetic mode, unlike the usual modified-two-stream instability, is not stabilized by electromagnetic effects when the relative electron-ion drift speed exceeds the Alfven speed.

Wu, C. S.↗

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↗

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.↗