Search NASA⌕ Search

SEARCH · Search NASA

Results for “Field boundaries”

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 127 records · Page 7

A numerical study of three-dimensional vortex breakdown

A numerical simulation of bubble-type vortex breakdown using a unique discrete form of the full 3-D, unsteady incompressible Navier-Stokes equations was performed. The Navier-Stokes equations were written in a vorticity-velocity form and the physical problem was not restricted to axisymmetric flow. The problem was parametized on a Rossby- Reynolds-number basis. Utilization of this parameter duo was shown to dictate the form of the free-field boundary condition specification and allowed control of axial breakdown location within the computational domain. The structure of the breakdown bubble was studied through time evolution plots of planar projected velocity vectors as well as through plots of particle traces and vortex lines. These results compared favorably with previous experimental studies. In addition, profiles of all three velocity components are presented at various axial stations and a Fourier analysis was performed to identify the dominant circumferential modes. The dynamics of the breakdown process were studied through plots of axial variation of rate of change of integrated total energy and rate of change of integrated enstrophy, as well as through contour plots of velocity, vorticity and pressure.

Spall, Robert E.↗

Advances in the computation of transonic separated flows over finite wings

The transonic flow around a low-aspect-ratio infinite wing in a wind tunnel is investigated by means of numerical simulations, with a focus on shock-induced separated flows. A coarse global grid with far-field boundaries matching those of the test section is subdivided into zones: the flow in the clustered zones near the wing is analyzed by solving the Reynolds-averaged Navier-Stokes equations, while that farther from the wing is modeled with the Euler equations. The results are presented graphically, and it is shown that a mushroomlike separated flow with two counterrotating vortices can be simulated when the correct shock strength is imposed (by careful selection of the artificial dissipation, the boundary conditions, the grid refinement, the algebraic turbulence model, and the geometry representation).

Kaynak, Unver↗

Unsteady transonic algorithm improvements for realistic aircraft applications

Improvements to a time-accurate approximate factorization (AF) algorithm were implemented for steady and unsteady transonic analysis of realistic aircraft configurations. These algorithm improvements were made to the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code developed at the Langley Research Center. The code permits the aeroelastic analysis of complete aircraft in the flutter critical transonic speed range. The AF algorithm of the CAP-TSD code solves the unsteady transonic small-disturbance equation. The algorithm improvements include: an Engquist-Osher (E-O) type-dependent switch to more accurately and efficiently treat regions of supersonic flow; extension of the E-O switch for second-order spatial accuracy in these regions; nonreflecting far field boundary conditions for more accurate unsteady applications; and several modifications which accelerate convergence to steady-state. Calculations are presented for several configurations including the General Dynamics one-ninth scale F-16C aircraft model to evaluate the algorithm modifications. The modifications have significantly improved the stability of the AF algorithm and hence the reliability of the CAP-TSD code in general.

Batina, John T.↗

Steady and unsteady transonic small disturbance analysis of realistic aircraft configurations

A transonic unsteady aerodynamic and aeroelasticity code CAP-TSD (Computational Aeroelasticity Program-Transonic Small Disturbance) has been developed for application to realistic aircraft configurations. The code uses a time-accurate approximate factorization (AF) algorithm for solution of the unsteady transonic small-disturbance potential equation. The paper gives an overview of the CAP-TSD code development effort and reports on recent algorithm modifications. The algorithm modifications include: an Engquist-Osher (E-O) type-dependent switch to treat regions of supersonic flow, extension of the E-O switch for second-order spatial accuracy, nonisentropic effects to treat strong-shock cases, nonreflecting far field boundary conditions for unsteady applications, and several modifications to accelerate convergence to steady state. The modifications have significantly improved the stability of the AF algorithm and hence the reliability of the CAP-TSD code in general. Calculations are also presented from a flutter analysis of a 45-deg sweptback wing which agree well with experimental data. The paper present descriptions of the CAP-TSD code and algorithm details along with results.

Batina, John T.↗

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for unsteady and steady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux-vector-splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flows results are made as well as with experimental data where available.

Anderson, W. Kyle↗

Unsteady transonic algorithm improvements for realistic aircraft applications

Improvements to a time-accurate approximate factorization (AF) algorithm have been implemented for steady and unsteady transonic analysis of realistic aircraft configurations. The algorithm improvements include: an Engquist-Osher (E-O) type-dependent switch to more accurately and efficiently treat regions of supersonic flow, extension of the E-O switch for second-order spatial accuracy in these regions, nonreflecting far field boundary conditions for more accurate unsteady applications, and several modifications which accelerate convergence to steady-state. Calculations are presented for several configurations to evaluate the algorithm modifications. The modifications have significantly improved the stability of the AF algorithm and hence the reliability of the CAP-TSD code in general.

Batina, John T.↗

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for steady and unsteady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux vector splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flow results are made, as well as with experimental data where available.

Anderson, W. Kyle↗

Numerical solution of periodic vortical flows about a thin airfoil

A numerical method is developed for computing periodic, three-dimensional, vortical flows around isolated airfoils. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Solutions for thin airfoils at zero degrees incidence to the mean flow are presented in this paper. Using an elliptic coordinate transformation, the computational domain is transformed into a rectangle. The Sommerfeld radiation condition is applied to the unsteady pressure on the grid line corresponding to the far field boundary. The results are compared with a Possio solver, and it is shown that for maximum accuracy the grid should depend on both the Mach number and reduced frequency. Finally, in order to assess the range of validity of the classical thin airfoil approximation, results for airfoils with zero thickness are compared with results for airfoils with small thickness.

Scott, James R.↗

Numerical solution of periodic vortical flows about a thin airfoil

A numerical method is developed for computing periodic, three-dimensional, vortical flows around isolated airfoils. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Solutions for thin airfoils at zero degrees incidence to the mean flow are presented in this paper. Using an elliptic coordinate transformation, the computational domain is transformed into a rectangle. The Sommerfield radiation condition is applied to the unsteady pressure on the grid line corresponding to the far field boundary. The results are compared with a Possion solver, and it is shown that for maximum accuracy the grid should depend on both the Mach number and reduced frequency. Finally, in order to assess the range of validity of the classical thin airfoil approximation, results for airfoils with zero thickness are compared with results for airfoils with small thickness.

Scott, James R.↗

A two-dimensional adaptive-wall test section with ventilated walls in the Ames 2- by 2-foot transonic wind tunnel

The first tests conducted in the adaptive-wall test section of the Ames Research Center's 2- by 2-Foot Transonic Wind Tunnel are described. A procedure was demonstrated for reducing wall interference in transonic flow past a two-dimensional airfoil by actively controlling flow through the slotted walls of the test section. Flow through the walls was controlled by adjusting pressures in compartments of plenums above and below the test section. Wall interference was assessed by measuring (with a laser velocimeter) velocity distributions along a contour surrounding the model, and then checking those measurements for their compatibility with free-air far-field boundary conditions. Plenum pressures for minimum wall interference were determined from empirical influence coefficients. An NACA 0012 airfoil was tested at angles of attach of 0 and 2, and at Mach numbers between 0.70 and 0.85. In all cases the wall-setting procedure greatly reduced wall interference. Wall interference, however, was never completely eliminated, primarily because the effect of plenum pressure changes on the velocities along the contour could not be accurately predicted.

Schairer, Edward T.↗

Numerical investigation of a plasma beam entering transverse magnetic fields

Two-dimensional electrostatic and electromagnetic particle-in-cell codes are used to investigate plasma beam propagation across transverse magnetic fields. Both large- and small-ion-gyroradius beams are studied in the case of small beam momentum or energy (low drift kinetic beta). Large-ion-gyroradius beams with a large dielectric constant epsilon much greater than (M/m) to the 1/2 (where M/m is the ion-to-electron mass ratio) are found to propagate across the magnetic field via E x B drifts at nearly the initial injection velocity. Beam degradation and undulations are observed, in agreement with previous experimental and analytical results. When epsilon is of order (M/m) to the 1/2 the plasma beam propagates across field lines at only half its initial velocity and loses its coherent structure. When epsilon is much less than (M/m) to the 1/2 the beam particles decouple at the magnetic field boundary, scattering the electrons and deflecting the ions slightly.

Koga, J.↗

Grid generation and inviscid flow computation about aircraft geometries

Grid generation and Euler flow about fighter aircraft are described. A fighter aircraft geometry is specified by an area ruled fuselage with an internal duct, cranked delta wing or strake/wing combinations, canard and/or horizontal tail surfaces, and vertical tail surfaces. The initial step before grid generation and flow computation is the determination of a suitable grid topology. The external grid topology that has been applied is called a dual-block topology which is a patched C (exp 1) continuous multiple-block system where inner blocks cover the highly-swept part of a cranked wing or strake, rearward inner-part of the wing, and tail components. Outer-blocks cover the remainder of the fuselage, outer-part of the wing, canards and extend to the far field boundaries. The grid generation is based on transfinite interpolation with Lagrangian blending functions. This procedure has been applied to the Langley experimental fighter configuration and a modified F-18 configuration. Supersonic flow between Mach 1.3 and 2.5 and angles of attack between 0 degrees and 10 degrees have been computed with associated Euler solvers based on the finite-volume approach. When coupling geometric details such as boundary layer diverter regions, duct regions with inlets and outlets, or slots with the general external grid, imposing C (exp 1) continuity can be extremely tedious. The approach taken here is to patch blocks together at common interfaces where there is no grid continuity, but enforce conservation in the finite-volume solution. The key to this technique is how to obtain the information required for a conservative interface. The Ramshaw technique which automates the computation of proportional areas of two overlapping grids on a planar surface and is suitable for coding was used. Researchers generated internal duct grids for the Langley experimental fighter configuration independent of the external grid topology, with a conservative interface at the inlet and outlet.

Smith, Robert E.↗

Utility of remotely sensed data for identification of soil conservation practices

Discussed are a variety of remotely sensed data sources that may have utility in the identification of conservation practices and related linear features. Test sites were evaluated in Alabama, Kansas, Mississippi, and Oklahoma using one or more of a variety of remotely sensed data sources, including color infrared photography (CIR), LANDSAT Thematic Mapper (TM) data, and aircraft-acquired Thermal Infrared Multispectral Scanner (TIMS) data. Both visual examination and computer-implemented enhancement procedures were used to identify conservation practices and other linear features. For the Kansas, Mississippi, and Oklahoma test sites, photo interpretations of CIR identified up to 24 of the 109 conservation practices from a matrix derived from the SCS National Handbook of Conservation Practices. The conservation practice matrix was modified to predict the possibility of identifying the 109 practices at various photographic scales based on the observed results as well as photo interpreter experience. Some practices were successfully identified in TM data through visual identification, but a number of existing practices were of such size and shape that the resolution of the TM could not detect them accurately. A series of computer-automated decorrelation and filtering procedures served to enhance the conservation practices in TM data with only fair success. However, features such as field boundaries, roads, water bodies, and the Urban/Ag interface were easily differentiated. Similar enhancement techniques applied to 5 and 10 meter TIMS data proved much more useful in delineating terraces, grass waterways, and drainage ditches as well as the features mentioned above, due partly to improved resolution and partly to thermally influenced moisture conditions. Spatially oriented data such as those derived from remotely sensed data offer some promise in the inventory and monitoring of conservation practices as well as in supplying parameter data for a variety of computer-implemented agricultural models.

Pelletier, R. E.↗

Electromagnetic scattering from arbitrary shaped conducting bodies coated with lossy materials of arbitrary thickness

A simple and efficient numerical technique is presented to solve the electromagnetic scattering problem of coated conducting bodies of arbitrary shape. The surface equivalence principle is used to formulate the problem in terms of a set of coupled integral equations involving equivalent electric and magnetic surface currents which represent boundary fields. The conducting structures and the dielectric materials are modeled by planar triangular patches, and the method of moments is used to solve the integral equations. Numerical results for scattering cross sections are given for various structures and compared with other available data. These results are proved accurate by a number of representative examples.

Rao, Sadasiva M.↗

One-way nesting for a primitive equation ocean model

Prognostic numerical models for atmospheric and oceanic circulations require initial fields, boundary conditions, and forcing functions in addition to a consistent set of partial differential equations, including a state relation and equations expressing conservation of mass, momentum, and energy. Depending on the horizontal domain to be modeled, the horizontal boundary conditions are either physically obvious or extremely difficult to specify consistently. If the entire atmosphere is modeled, periodic horizontal boundary conditions are appropriate. On the other hand, the physical horizontal boundaries on the entire ocean are solid walls. Obviously, the normal velocity at a solid wall is zero while the specification of the tangential velocity depends on the mathematical treatment of the horizontal viscous terms. Limitations imposed by computer capacity and cost, as well as research interests, have led to the use of limited area models to study flows in the atmosphere and ocean. The limited area models do not have physical horizontal boundaries, merely numerical ones. Correctly determining these open boundary conditions for limited-area numerical models has both intrigued and frustrated numerical modelers for decades. One common approach is to use the closed or solid wall boundary conditions for a limited-area model. The argument given for this approach is that the boundary conditions affect flow near the walls but that none of these effects are propagated into the interior. Therefore, one chooses a big enough domain that the central region of interest is not corrupted by the boundary flow. Research in progress to model the North Atlantic circulation vividly illustrates the pitfalls of this approach. Two model runs are compared: (1) the southern boundary at 20S between latitudes 0 and 40W is artificially closed; and (2) the same boundary is specified as open with an inward transport of 15 Sv (determined from a global model with the same physics) uniformly spread across the boundary. A comparison of both runs is presented.

Blake, D. W.↗

Two-dimensional adaptive-wall tests in the NASA Ames two- by two-foot transonic wind tunnel

The first tests conducted in the adaptive-wall test section of the Ames Research Center's 2- by 2-Foot Transonic Wind Tunnel are described. A procedure was demonstrated for reducing wall interference in transonic flow past a two-dimensional airfoil by actively controlling flow through the slotted walls of the test section. Flow through the walls was controlled by adjusting pressures in compartments of plenums above and below the test section. Wall interference was assessed by measuring (with a laser velocimeter) velocity distributions along a contour surrounding the model, and then checking those measurements for their compatibility with free-air far-field boundary conditions. Plenum pressures for minimum wall interference were determined from empirical influence coefficients. An NACA 0012 airfoil was tested at angles of attack of 0 and 2, and at Mach numbers between 0.70 and 0.85. In all cases the wall-setting procedure greatly reduced wall interference. Wall interference, however, was never completely eliminated, primarily because the effect of plenum pressure changes on the velocities along the contour could not be accurately predicted.

Schairer, Edward T.↗

Evidence for mass outflow in the low solar corona over a large sunspot

Spatially resolved EUV coronal emission-line profiles have been obtained in a solar active region, including a large sunspot, using an EUV imaging spectrograph. Relative Doppler velocities were measured in the lines of Mg IX, Fe XV, and Fe XVI with a sensitivity of 2-3 km/s at 350 A. The only significant Doppler shift occurred over the umbra of the large sunspot, in the emission line of Mg IX (at Te of about 1.1 x 10 exp 6 K). The maximum shift corresponded to a peak velocity toward the observer of 14 +/- 3 km/s relative to the mean of measurements in this emission line made elsewhere over the active region. The magnetic field in the low corona was aligned to within 10 deg of the line of sight at the location of maximum Doppler shift. Depending on the closure of the field, such a mass flow could either contribute to the solar wind or reappear as a downflow of material in distant regions on the solar surface. The site of the source, near a major photospheric field boundary, was consistent with origins of low-speed solar wind typically inferred from interplanetary plasma observations.

Neupert, Werner M.↗

The 24 GHz measurements of 2.2 lambda conical horn antennas illuminating a conducting sheet

Monostatic reflection-coefficient magnitude, absolute value of Gamma, measurements occurring between a radiating horn and a metal reflecting plate are presented for a family of three 2.2 lambda diameter conical horn antennas. The three horns have different aperture phase deviations: 6 deg, 22.5 deg, and 125 deg. Measurements of the magnitude of absolute value of Gamma as a function of horn-plate separation (d) extend from an effective antenna aperture short (d = O) to beyond the far-field boundary (d = 2D(sup 2)/lambda, where D is the antenna diameter). Measurement data are presented with various physical environments for each of the horns. Measured scalar data are compared with theoretical data from two models, a numerical model for a circular waveguide aperture in a ground plane and a scalar diffraction theory model. This work was conducted in support of the development effort for a spaceborne multifrequency microwave reflectometer designed to accurately determine the distance from a space vehicle's surface to a reflecting plasma boundary. The metal reflecting plate was used to simulate the RF reflectivity of a critically dense plasma. The resulting configuration, a ground plane mounted aperture facing a reflecting plane in close proximity, produces a strong interaction between the ground plane and the reflecting plate, especially at integral half-wavelength separations. The transition coefficient is characterized by large amplitude variations.

Cross, A. E.↗