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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Natural convection in steady solidification - Finite element analysis of a two-phase Rayleigh-Benard problem

Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.

Chang, C. J.↗

NASA NTRS1984

On the stability of the decelerating laminar boundary layer

The stability of a decelerating boundary-layer flow is investigated experimentally and numerically. Experimentally, a flat plate having a Blasius boundary layer is decelerated in an 18 m towing tank. The boundary layer becomes unstable to two-dimensional waves, which break down into three-dimensional patterns, hairpin vortices, and finally turbulent bursts when the vortices lift off the wall. The unsteady boundary-layer equations are solved numerically to generate instantaneous velocity profiles for a range of boundary and initial conditions. A quasi-steady approximation is invoked and the stability of local velocity profiles is determined by solving the Orr-Sommerfeld equation using Chebyshev matrix methods. Comparisons are made between the numerical predictions and the experimentally observed instabilities.

Gad-El-hak, M.↗

NASA NTRS1983

Computation of viscous flows in turbomachinery cascades

Analytic and numerical experiments involving the use of the space-marching method for the prediction of the viscous flow field in turbomachinery cascades indicate that the grid lines in the physical domain must be nearly orthogonal in order to arrive at accurate predictions. For staggered cascades, a nonperiodic grid system was used and appropriate approximations were applied in place of the periodicity boundary conditions upstream and downstream of the cascade. Predicted drag coefficient, turning angle, boundary layer momentum thicknesses, and velocity profiles are compared with experimental data. Good agreement is obtained in most cases.

Pouagare, M.↗

NASA NTRS1985

Transient convection over the Amazon-Bolivia region and the dynamics of drought over Northeast Brazil

A two layer, nonlinear, equatorial beta-plane model, in p-coordinates is used to study the atmospheric response to a large scale prescribed heat source varying in time. The heat source is meant to represent a convective burst with total duration of approximately 48 hours over the Amazon/Bolivia region. The boundary conditions used are meridional velocity zero at 60 deg S, omega w = 0 at the top and zero geometric velocity at the lower boundary. Sensitivity study was done which includes initial state at rest, compared with realistic initial flow. The scale of the heat source is 1500 km in latitude and longitude and it is centered at 10 deg S. Special attention is paid to the distribution and intensity of the induced vertical motion. The model is integrated for two days and the preliminary results show agreement with the observed 200 mb flow. Of interest is the establishment of a trough and descending motion to the northeast of the heat source. A conjucture is thus made that the Amazon heat source and its fluctuations bear some relationship with the drought problem over Northeast Brazil.

Buchmann, J.↗

NASA NTRS1985

On the solution of integral equations with strong ly singular kernels

In this paper some useful formulas are developed to evaluate integrals having a singularity of the form (t-x) sup-m, m or = 1. Interpreting the integrals with strong singularities in Hadamard sense, the results are used to obtain approximate solutions of singular integral equations. A mixed boundary value problem from the theory of elasticity is considered as an example. Particularly for integral equations where the kernel contains, in addition to the dominant term (t,x) sup-m, terms which become unbounded at the end points, the present technique appears to be extremely effective to obtain rapidly converging numerical results.

Kaya, A. C.↗

NASA NTRS1986

On the solution of integral equations with strongly singular kernels

Some useful formulas are developed to evaluate integrals having a singularity of the form (t-x) sup-m ,m greater than or equal 1. Interpreting the integrals with strong singularities in Hadamard sense, the results are used to obtain approximate solutions of singular integral equations. A mixed boundary value problem from the theory of elasticity is considered as an example. Particularly for integral equations where the kernel contains, in addition to the dominant term (t-x) sup -m , terms which become unbounded at the end points, the present technique appears to be extremely effective to obtain rapidly converging numerical results.

Kaya, A. C.↗

NASA NTRS1986

Improved Euler analysis of advanced turboprop propeller flows

An implicit approximate factorization scheme in conjunction with a new boundary treatment was used to compute the inviscid flow field about an advanced high-speed propeller. The method of characteristics was used to apply impermeable boundary conditions. The convergence history of numerical calculations shows substantial decrease in the residual error decay. Euler solutions were computed for SR-3 propeller geometry and the results were compared with the experimental data and previous numerical results.

Yamamoto, O.↗

NASA NTRS1987

Temperature fields due to jet induced mixing in a typical OTV tank

The Eclipse Code is being developed as a general tool for analysis of cryogenic propellant behavior in spacecraft tankage. The focus of the work being reported is on prediction of temperature fields due to introduction of a cold jet along the centerline of a typical Orbit Transfer Vehicle tank. A brief description of the formulations used for modeling heat transfer and turbulent flow is presented. Code performance is verified through comparison to experimental data for mixing in small scale tanks. An unexpected difficulty in computing long duration flows is reviewed. Preliminary results for a partially filled full scale tank are obtained by approximating the free surface by a spherical solid boundary.

Hochstein, J. I.↗

NASA NTRS1987

The design and preliminary calibration of a boundary-layer flow channel

Design procedures for a new low-speed boundary-layer research channel are described. The channel is an open-circuit wind tunnel for the study of two-dimensional boundary layers under controlled pressure gradients, and follows design guidelines from published literature on blower tunnels with wide-angle diffusers. The contraction was arranged in a modular fashion permitting two different test sections of square and high-aspect-ratio cross section. A radical type of wide-angle diffuser was employed, and a stream-tube computer code (GE Streamtube Curvature Code) was used to check the contraction designs. The alternate test sections have the following specifications: 2- by 2-foot cross section with a fixed velocity of 23 ft/sec, and a boundary-layer section with a 0.5- by 2-foot cross section at a fixed velocity of approximately 89 ft/sec. Experimental techniques and data are described for the evaluation of diffuser effectiveness, boundary-layer channel characteristics, and overall performance of the facility.

Rao, Dhanvada M.↗

NASA NTRS1987

On the solution of integral equations with strongly singular kernels

Some useful formulas are developed to evaluate integrals having a singularity of the form (t-x) sup-m, m greater than or equal 1. Interpreting the integrals with strong singularities in Hadamard sense, the results are used to obtain approximate solutions of singular integral equations. A mixed boundary value problem from the theory of elasticity is considered as an example. Particularly for integral equations where the kernel contains, in addition to the dominant term (t-x) sup-m, terms which become unbounded at the end points, the present technique appears to be extremely effective to obtain rapidly converging numerical results.

Kaya, A. C.↗

NASA NTRS1989

Interactive boundary-layer calculations of a transonic wing flow

Results obtained from iterative solutions of inviscid and boundary-layer equations are presented and compared with experimental values. The calculated results were obtained with an Euler code and a transonic potential code in order to furnish solutions for the inviscid flow; they were interacted with solutions of two-dimensional boundary-layer equations having a strip-theory approximation. Euler code results are found to be in better agreement with the experimental data than with the full potential code, especially in the presence of shock waves, (with the sole exception of the near-tip region).

Kaups, Kalle↗

NASA NTRS1990

Dynamics of the outgoing turbulent boundary layer in a Mach 5 unswept compression ramp interaction

Wall pressure fluctuations have been measured under the unsteady separation shock and on the ramp face in an unswept Mach 5 compression ramp interaction. The freestream Reynolds number was 51.0 x 10 to the 6th/m, and the incoming turbulent boundary layer developed on the tunnel floor under approximately adiabatic wall temperature conditions. Standard data-acquisition methods, as well as real-time and posttest conditional sampling techniques were used. The results show that the mean and rms pressure levels are strong functions of separation shock position. At all stations on the ramp, from the corner to where the pressure reaches the theoretical inviscid value, the pressure signals have two dominant components: a low frequency component characteristic of the global unsteadiness, which correlates with the separation shock motion, and a higher frequency component associated with turbulence. The former is the major contributor to the overall signal variance.

Gramann, Richard A.↗

NASA NTRS1991

Doppler lidar for measurement of atmospheric wind fields

Measurements of wind fields in the earth's troposphere with daily global coverage is widely considered as a significant advance for forecasting and transport studies. For optimal use by NWP (Numerical Weather Prediction) models the horizontal and vertical resolutions should be approximately 100 km and 1 km, respectively. For boundary layer studies vertical resolution of a few hundred meters seems essential. Earth-orbiting Doppler lidar has a unique capability to measure global winds in the troposphere with the high vertical resolution required. The lidar approach depends on transmission of pulses with high spectral purity and backscattering from the atmospheric aerosol particles or layered clouds to provide a return signal. Recent field measurement campaigns using NASA research aircraft have resulted in collection of aerosol and cloud data which can be used to optimize the Doppler lidar instrument design and measurement strategy.

Menzies, Robert T.↗

NASA NTRS1993

Summation by parts, projections, and stability

We have derived stability results for high-order finite difference approximations of mixed hyperbolic-parabolic initial-boundary value problems (IBVP). The results are obtained using summation by parts and a new way of representing general linear boundary conditions as an orthogonal projection. By slightly rearranging the analytic equations, we can prove strict stability for hyperbolic-parabolic IBVP. Furthermore, we generalize our technique so as to yield strict stability on curvilinear non-smooth domains in two space dimensions. Finally, we show how to incorporate inhomogeneous boundary data while retaining strict stability. Using the same procedure one can prove strict stability in higher dimensions as well.

Olsson, Pelle↗

NASA NTRS1994

Radial diffusion models of energetic electrons and Jupiter's synchrotron radiation. 2: Time variability

We used a radial diffusion code for energetic electrons in Jupiter's magnetosphere to investigate variations in Jupiter's radio emission due to changes in the electron phase space density at L shells between 6 and 50, and due to changes in the radial diffusion parameters. We suggest that the observed variations in Jupiter's radio emission are likely caused by changes in the electron phase space density at some boundary L(sub 1) is greater than 6, if the primary mode of transport of energetic electrons is radial diffusion driven by fluctuating electric and/or magnetic fields induced by upper atmospheric turbulence. We noticed an excellent empirical correlation, both in phase and relative amplitude, between changes in the solar wind ram pressure and Jupiter's synchrotron radiation if the electron phase space density at the boundary L(sub 1) (L(sub 1) is approximately equal to 20-50) varies linearly with the square root of the solar wind ram pressure, f is approximately (N(sub s)nu(exp 2 sub s))(exp 1/2). The calculations were carried out with a diffusion coefficient D(sub LL) = D(sub n)L(exp n) with n = 3. The diffusion coefficient which best fit the observed variations in Jupiter's synchrotron radiation D(sub 3) = 1.3 +/- 0.2 x 10(exp -9)/s is approximately 0.041/yr, which corresponds to a lagtime of approximately 2 years. We further show that the observed short term (days-weeks) variations in Jupiter's radio emission cannot be explained adequately when radial diffusion is taken into account.

De Pater, Imke↗

NASA NTRS1994

MHD simulations of boundary layer formation along the dayside Venus ionopause due to mass loading

A two-dimensional magnetohydrodynamic (MHD) simulation of mass-loaded solar wind flow around the dayside of Venus is presented. For conditions appropriate to a low-altitude ionopause the simulations show that mass loading from the pickup of oxygen ions produces a boundary layer of finite thickness along the ionopause. Within this layer the temperatures exhibit strong gradients normal to and away from the ionopause. Furthermore, there is a shear in the bulk flow velocity across the boundary layer, such that the (predominantly tangential) flow decreases in speed as the ionopause is approached and remains small along the ionopause, consistent with Pioneer Venus observations. The total mass density increases significantly as the flow approaches the ionopause, where the contribution of O(+) to the total number density is a few percent. Numerical simulations are carried out for various mass addition rates and demonstrate that the boundary layer develops when oxygen ion production exceeds approximately 2 x 10(exp 5)/cu m/s. For the upstream solar wind parameters and mass loading rates chosen for these simulations, the results are consistent with observations made on the dayside of Venus for average ionopause conditions near 300 km.

Mcgary, J. E.↗

NASA NTRS1995

A new method of presentation the large-scale magnetic field structure on the Sun and solar corona

The large-scale photospheric magnetic field, measured at Stanford, has been analyzed in terms of surface harmonics. Changes of the photospheric field which occur within whole solar rotation period can be resolved by this analysis. For this reason we used daily magnetograms of the line-of-sight magnetic field component observed from Earth over solar disc. We have estimated the period during which day-to-day full disc magnetograms must be collected. An original algorithm was applied to resolve time variations of spherical harmonics that reflect time evolution of large-scale magnetic field within solar rotation period. This method of magnetic field presentation can be useful enough in lack of direct magnetograph observations due to sometimes bad weather conditions. We have used the calculated surface harmonics to reconstruct the large-scale magnetic field structure on the source surface near the sun - the origin of heliospheric current sheet and solar wind streams. The obtained results have been compared with spacecraft in situ observations and geomagnetic activity. We tried to show that proposed technique can trace shon-time variations of heliospheric current sheet and short-lived solar wind streams. We have compared also our results with those obtained traditionally from potential field approximation and extrapolation using synoptic charts as initial boundary conditions.

Ponyavin, D. I.↗

NASA NTRS1997

A Numerical Method for Incompressible Flow with Heat Transfer

A numerical method for the convective heat transfer problem is developed for low speed flow at mild temperatures. A simplified energy equation is added to the incompressible Navier-Stokes formulation by using Boussinesq approximation to account for the buoyancy force. A pseudocompressibility method is used to solve the resulting set of equations for steady-state solutions in conjunction with an approximate factorization scheme. A Neumann-type pressure boundary condition is devised to account for the interaction between pressure and temperature terms, especially near a heated or cooled solid boundary. It is shown that the present method is capable of predicting the temperature field in an incompressible flow.

Sa, Jong-Youb↗