SEARCH · Search NASA

Results for “free boundary problem”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Variational formulation of hybrid problems for fully 3-D transonic flow with shocks in rotor

Based on previous research, the unified variable domain variational theory of hybrid problems for rotor flow is extended to fully 3-D transonic rotor flow with shocks, unifying and generalizing the direct and inverse problems. Three variational principles (VP) families were established. All unknown boundaries and flow discontinuities (such as shocks, free trailing vortex sheets) are successfully handled via functional variations with variable domain, converting almost all boundary and interface conditions, including the Rankine Hugoniot shock relations, into natural ones. This theory provides a series of novel ways for blade design or modification and a rigorous theoretical basis for finite element applications and also constitutes an important part of the optimal design theory of rotor bladings. Numerical solutions to subsonic flow by finite elements with self-adapting nodes given in Refs., show good agreement with experimental results.

Liu, Gao-Lian↗

NASA NTRS1992

Bypass transition in compressible boundary layers

Transition to turbulence in aerospace applications usually occurs in a strongly disturbed environment. For instance, the effects of free-stream turbulence, roughness and obstacles in the boundary layer strongly influence transition. Proper understanding of the mechanisms leading to transition is crucial in the design of aircraft wings and gas turbine blades, because lift, drag and heat transfer strongly depend on the state of the boundary layer, laminar or turbulent. Unfortunately, most of the transition research, both theoretical and experimental, has focused on natural transition. Many practical flows, however, defy any theoretical analysis and are extremely difficult to measure. Morkovin introduced in his review paper the concept of bypass transition as those forms of transition which bypass the known mechanisms of linear and non-linear transition theories and are currently not understood by experiments. In an effort to better understand the mechanisms leading to transition in a disturbed environment, experiments are conducted studying simpler cases, viz. the effects of free stream turbulence on transition on a flat plate. It turns out that these experiments are very difficult to conduct, because generation of free stream turbulence with sufficiently high fluctuation levels and reasonable homogeneity is non trivial. For a discussion see Morkovin. Serious problems also appear due to the fact that at high Reynolds numbers the boundary layers are very thin, especially in the nose region of the plate where the transition occurs, which makes the use of very small probes necessary. The effects of free-stream turbulence on transition are the subject of this research and are especially important in a gas turbine environment, where turbulence intensities are measured between 5 and 20 percent, Wang et al. Due to the fact that the Reynolds number for turbine blades is considerably lower than for aircraft wings, generally a larger portion of the blade will be in a laminar transitional state. The effects of large free stream turbulence in compressible boundary layers at Mach numbers are examined both in the subsonic and transonic regime using direct numerical simulations. The flow is computed over a flat plate and curved surface. while many applications operate in the transonic regime. Due the nature of their numerical scheme, a non-conservation formulation of the Navier-Stokes equations, it is a non-trivial extension to compute flow fields in the transonic regime. This project aims at better understanding the effects of large free-stream turbulence in compressible boundary layers at mach number both in the subsonic and transonic regime using direct numerical simulations. The present project aims at computing the flow over a flat plate and curved surface. This research will provide data which can be used to clarify mechanisms leading to transition in an environment with high free stream turbulence. This information is useful for the development of turbulence models, which are of great importance for CFD applications, and are currently unreliable for more complex flows, such as transitional flows.

Vandervegt, J. J.↗

NASA NTRS1991

Linear and nonlinear dynamic analysis by boundary element method

An advanced implementation of the direct boundary element method (BEM) applicable to free-vibration, periodic (steady-state) vibration and linear and nonlinear transient dynamic problems involving two and three-dimensional isotropic solids of arbitrary shape is presented. Interior, exterior, and half-space problems can all be solved by the present formulation. For the free-vibration analysis, a new real variable BEM formulation is presented which solves the free-vibration problem in the form of algebraic equations (formed from the static kernels) and needs only surface discretization. In the area of time-domain transient analysis, the BEM is well suited because it gives an implicit formulation. Although the integral formulations are elegant, because of the complexity of the formulation it has never been implemented in exact form. In the present work, linear and nonlinear time domain transient analysis for three-dimensional solids has been implemented in a general and complete manner. The formulation and implementation of the nonlinear, transient, dynamic analysis presented here is the first ever in the field of boundary element analysis. Almost all the existing formulation of BEM in dynamics use the constant variation of the variables in space and time which is very unrealistic for engineering problems and, in some cases, it leads to unacceptably inaccurate results. In the present work, linear and quadratic isoparametric boundary elements are used for discretization of geometry and functional variations in space. In addition, higher order variations in time are used. These methods of analysis are applicable to piecewise-homogeneous materials, such that not only problems of the layered media and the soil-structure interaction can be analyzed but also a large problem can be solved by the usual sub-structuring technique. The analyses have been incorporated in a versatile, general-purpose computer program. Some numerical problems are solved and, through comparisons with available analytical and numerical results, the stability and high accuracy of these dynamic analysis techniques are established.

Ahmad, Shahid↗

NASA NTRS1987

Nonlinear force-free magnetic fields: Calculation and applicatin to astrophysics

The problem concerned in this work is that of calculating magnetic field configurations in which the Lorentz force (vector)j x (vector)B is everywhere zero, subject to specified boundary conditions. The magnetic field is represented in terms of Clebsch variables in the form (vector)B = del alpha x del beta. These variables are constant on any field line. The most appropriate choice of boundary conditions is to specify the values of alpha and beta on the bounding surface. It is proposed that the field lines move in the direction of local Lorentz force and relax towards a force-free field configuration. This concept leads to an iteration procedure for modifying the variables alpha and beta that tends asymptotically towards the force-free state. This method is first applied to a simple problem in two rectangular dimensions; the calculation shows that the convergence of magnetic field energy to a minimum state (force-free) is close to exponential. This method is then applied to study some astrophysical force-free magnetic fields, such as the structures and evolution of magnetic fields of rotating sunspots and accretion disks. The implication of the results, as related to the mechanisms of solar flares, extragalactic radio sources and radio jets, are discussed.

Yang, Wei-Hong↗

NASA NTRS1982

Computer prediction of three-dimensional potential flow fields in which aircraft propellers operate

A computer program was developed to calculate the three dimensional, steady, incompressible, inviscid, irrotational flow field at the propeller plane (propeller removed) located upstream of an arbitrary airframe geometry. The program uses a horseshoe vortex of known strength to model the wing. All other airframe surfaces are modeled by a network source panels of unknown strength which is exposed to a uniform free stream and the wing-induced velocity field. By satisfying boundary conditions on each panel (the Neumann problem), relaxed boundary conditions being used on certain panels to simulate inlet inflow, the source strengths are determined. From the known source and wing vortex strengths, the resulting velocity fields on the airframe surface and at the propeller plane are obtained. All program equations are derived in detail, and a brief description of the program structure is presented. A user's manual which fully documents the program is cited. Computer predictions of the flow on the surface of a sphere and at a propeller plane upstream of the sphere are compared with the exact mathematical solutions. Agreement is good, and correct program operation is verified.

Jumper, S. J.↗

NASA NTRS1985

An internal variable assisted formulation of moisture sorption models for polymers and composites

The purpose of the paper is to demonstrate that constitutive equations for moisture sorption and swelling can be obtained with an internal variable method which is rooted in the theory of irreversible thermodynamics. Realism is injected into the model through the concept of free volume. The power of the method is illustrated by the solution of the formidable nonlinear boundary value problem for diffusion with swelling-stress dependent diffusivity.

Hiel, C. C.↗

NASA NTRS1987

An improved method of Nusselt number calculation

A novel method for calculating the Nusselt number, Nu, in a steady-state Rayleigh-Benard convection problem is presented, in which calculations are done for a square box with constant temperature, free-slip boundary conditions at the top and bottom, and a reflection symmetry along the side walls. The element heat flux is obtained by averaging over the entire element; element heat fluxes are then projected to the adjacent nodes. Compared with previous methods, the approach reduces the calculated depth variation in horizontally averaged flux by more than a factor of 10 and shows more rapid convergence of Nu as a function of grid size.

Ho-Liu, Phyllis↗

NASA NTRS1990

Dynamic substructuring by the boundary flexibility vector method of component mode synthesis

Component mode synthesis (CMS) is a method of dynamic analysis, for structures having a large number of degrees of freedom (DOF). These structures often required lengthy computer CPU time and large computer memory resources, if solved directly by the finite-element method (FEM). In CMS, the structure is divided into independent components in which the DOF are defined by a set of generalized coordinates defined by displacement shapes. The number of the generalized coordinates are much less than the original number of physical DOF, in the component. The displacement shapes are used to transform the component property matrices and any applied external loads, to a reduced system of coordinates. Reduced system property matrices are assembled, and any type of dynamic analysis is carried out in the reduced coordinate system. Any obtained results are back transformed to the original component coordinate systems. In all conventional methods of CMS, the mode shapes used for components are dynamic mode shapes, supplemented by static deflected shapes. Historically, all the dynamic mode shapes used in conventional CMS are the natural modes (eigenvectors) of components. A new method of CMS, namely the boundary flexibility vector method of CMS, is presented. The method provides for the incorporation of a set of static Ritz vectors, referred to as boundary flexibility vectors, as a replacement and/or supplement to conventional eigenvectors, as displacement shapes for components. The generation of these vectors does not require the solution of a costly eigenvalue problem, as in the case of natural modes in conventional CMS, and hence a substantial saving in CPU time can be achieved. The boundary flexibility vectors are generated from flexibility (or stiffness) properties of components. The formulation presented is for both free and fixed-interface components, and for both the free and forced vibration problems. Free and forced vibration numerical examples are presented to verify the accuracy of the method and the saving in CPU time. Compared to conventional methods of CMS, the results indicate that by using the new method, more accurate results can be obtained with a substantial saving in CPU time.

Abdallah, Ayman Ahmed↗

NASA NTRS1977

Stress-free end problem in layered materials

In this paper the plane elastostatic problem for a medium which consists of periodically arranged two sets of bonded dissimilar layers or strips is considered. First it is assumed that one set of strips contains a crack which crosses the bimaterial interfaces. Then, by letting the collinear cracks join, the stress-free end problem is formulated. The singular behavior of the solutions at the point on intersection of the stress-free boundary and the interfaces is examined and appropriate stress intensity factors are defined. The results of some numerical examples are then presented which include the cases of both plane stress and plane strain.

Erdogan, F.↗

NASA NTRS1981

Eigenvalues of the Rayleigh-Benard and Marangoni problems

The eigenvalues of the linear Benard-Marangoni stability problem are discussed. Pearson and Nield boundary conditions, which correspond to a rigid, isothermal lower boundary and a stress-free conducting upper boundary are considered. It is shown that although a critical value of the Marangoni number can be determined, the number is not, strictly speaking, an eigenvalue and cannot be used as an eigenvalue parameter for the determination of an eigenvector set.

Rosenblat, S.↗

NASA NTRS1986

The construction and use of divergence free vector expansions for incompressible fluid flow calculations

For incompressible fluids the law of mass conservation reduces to a constraint on the velocity vector, namely that it be divergence free. This constraint has long been a source of great difficulty to the numericist seeking to discretize the Navier-Stokes and Euler equations. A spectral method is discussed which overcomes this difficulty. Its efficacy is demonstrated on some simple problems. The velocity is approximated by a finite sum of divergence free vectors, each of which satisfies the same boundary conditions as the velocity. Projecting the governing equation onto the space of inviscid vector fields eliminates the pressure term and produces a set of ordinary differential equations that must be solved for the coefficents in the velocity. The pressure can then be recovered if it is needed.

Mhuiris, N. M. G.↗

NASA NTRS1984

Probabilistic radiative transfer

A computationally efficient method has been developed for highly nonlinear problems in which radiative transfer is an important aspect of the heating and cooling of the medium. An approximate probabilistic radiative transfer equation is derived for one-dimensional plane-parallel atmospheres of finite or semi-infinite extent, for both spectral lines and bound-free continua. Boundary conditions, accuracy, escape probabilities, and practical aspects of complete linearization, are discussed. The method is accurate to a few tens of percent of a wide variety of realistic problems in which frequency redistribution of scattered photons dominates the transfer and escape of radiation.

Canfield, R. C.↗

NASA NTRS1984

Surface cracks in a plate of finite width under extension or bending

In this paper the problem of a finite plate containing collinear surface cracks is considered. The problem is solved by using the line spring model with plane elasticity and Reissner's plate theory. The main purpose of the study is to investigate the effect of interaction between two cracks or between cracks and stress-free plate boundaries on the stress intensity factors and to provide extensive numerical results which may be useful in applications. First, some sample results are obtained and are compared with the existing finite element results. Then the problem is solved for a single (internal) crack, two collinear cracks and two corner cracks for wide range of relative dimensions. Particularly in corner cracks the agreement with the finite element solution is surprisingly very good. The results are obtained for semielliptic and rectangular crack profiles which may, in practice, correspond to two limiting cases of the actual profile of a subcritically growing surface crack.

Erdogan, F.↗

NASA NTRS1956

An Engineering Approach to the Variable Fluid Property Problem in Free Convection

An analysis is made for the variable fluid property problem for laminar free convection on an isothermal vertical flat plate. For a number of specific cases, solutions of the boundary layer equations appropriate to the variable property situation were carried out for gases and liquid mercury. Utilizing these findings, a simple and accurate shorthand procedure is presented for calculating free convection heat transfer under variable property conditions. This calculation method is well established in the heat transfer field. It involves the use of results which have been derived for constant property fluids, and of a set of rules (called reference temperatures) for extending these constant property results to variable property situations. For gases, the constant property heat transfer results are generalized to the variable property situation by replacing beta (expansion coefficient) by one over T sub infinity and evaluating the other properties at T sub r equals T sub w minus zero point thirty-eight (T sub w minus T sub infinity). For liquid mercury, the generalization may be accomplished by evaluating all the properties (including beta) at this same T sub r. It is worthwhile noting that for these fluids, the film temperature (with beta equals one over T sub infinity for gases) appears to serve as an adequate reference temperature for most applications. Results are also presented for boundary layer thickness and velocity parameters.

Gregg, J. L.↗

NASA NTRS1990

On the instability of hypersonic flow past a flat plate

The instability of hypersonic boundary-layer flows over flat plates is considered. The viscosity of the fluid is taken to be governed by Sutherland's law, which gives a much more accurate representation of the temperature dependence of fluid viscosity at hypersonic speeds than Chapman's approximate linear law; although at lower speeds the temperature variation of the mean state is less pronounced so that the Chapman law can be used with some confidence. Attention is focussed on the so-called (vorticity) mode of instability of the viscous hypersonic boundary layer. This is thought to be the fastest growing inviscid disturbance at hypersonic speeds; it is also believed to have an asymptotically larger growth rate than any viscous or centrifugal instability. As a starting point the instability of the hypersonic boundary layer which exists far downstream from the leading edge of the plate is investigated. In this regime the shock that is attached to the leading edge of the plate plays no role, so that the basic boundary layer is non-interactive. It is shown that the vorticity mode of instability of this flow operates on a significantly different lengthscale than that obtained if a Chapman viscosity law is assumed. In particular, it is found that the growth rate predicted by a linear viscosity law overestimates the size of the growth rate by O(M(exp 2). Next, the development of the vorticity mode as the wavenumber decreases is described, and it is shown that acoustic modes emerge when the wavenumber has decreased from it's O(1) initial value to O(M (exp -3/2). Finally, the inviscid instability of the boundary layer near the leading edge in the interaction zone is discussed and particular attention is focussed on the strong interaction region which occurs sufficiently close to the leading edge. It is found that the vorticity mode in this regime is again unstable, and that it is concentrated in the transition layer at the edge of the boundary layer where the temperature adjusts from its large, O(M(exp 2), value in the viscous boundary layer, to its O(1) free stream value. The existence of the shock indirectly, but significantly, influences the instability problem by modifying the basic flow structure in this layer.

Blackaby, Nicholas↗

NASA NTRS1991

The unsteady laminar boundary layer on an axisymmetric body subject to small-amplitude fluctuations in the free-stream velocity

The effect of small-amplitude, time-periodic, free-stream disturbances on an otherwise steady axisymmetric boundary layer on a circular cylinder is considered. Numerical solutions to the problem are presented, and an asymptotic solution to the flow, valid far downstream along the axis of the cylinder is detailed. Particular emphasis is placed on the unsteady eigensolutions that occur far downstream, which turn out to be very different from the analogous planar eigensolutions. These axisymmetric eigensolutions are computed numerically and also are described by asymptotic analyses valid for low and high frequencies of oscillation.

Duck, Peter W.↗

NASA NTRS1994

Role of pressure diffusion in non-homogeneous shear flows

A non-local model is presented for approximating the pressure diffusion in calculations of turbulent free shear and boundary layer flows. It is based on the solution of an elliptic relaxation equation which enables local diffusion sources to be distributed over lengths of the order of the integral scale. The pressure diffusion model was implemented in a boundary layer code within the framework of turbulence models based on both the kappa-epsilon-(bar)upsilon(exp 2) system of equations and the full Reynolds stress equations. Model computations were performed for mixing layers and boundary layer flows. In each case, the pressure diffusion model enabled the well-known free-stream edge singularity problem to be eliminated. There was little effect on near-wall properties. Computed results agreed very well with experimental and DNS data for the mean flow velocity, the turbulent kinetic energy, and the skin-friction coefficient.

Demuren, A. O.↗

NASA NTRS1963

INVESTIGATION OF PROBLEMS IN THERMAL CONVECTION

Numerical solution of thermal convection equations for velocity and temperature fluctuations as well as the stability analysis and heat transport for free boundary conditions

THERMAL CONVECTION↗