Hypersonic laminar boundary layers around slender bodies.
Nonsimilar compressible laminar boundary layer equations for hypersonic flow around slender bodies solved by integration based on steepest descent method
Jeng, D. R.↗
SEARCH · Search NASA
Results for “boundary integral equation”
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.
Experimental and theoretical investigation of three-dimensional turbulent boundary layers and turbulence characteristics inside an axial flow inducer passage
Analytical and experimental investigations of the characteristics of three dimensional turbulent boundary layers in a rotating helical passage of an inducer rotor are reported. Expressions are developed for the velocity profiles in the inner layer, where the viscous effects dominate, in the outer layer, where the viscous effects are small, and in the interference layer, where the end walls influence the flow. The prediction of boundary layer growth is based on the momentum integral technique. The equations derived are general enough to be valid for all turbomachinery rotors with arbitrary pressure gradients. The experimental investigations are carried out in a flat plate inducer 3 feet in diameter. The mean velocity profiles, turbulence intensities and shear stresses, wall shear stress, and limiting streamline angles are measured at various radial and chordwise locations by using rotating probes. The measurements are in general agreement with the predictions. The radial flows are well represented by an expression which includes the effect of stagger angle and radial pressure gradient. The radial flows in the rotor channel are higher than those on a single blade. The collateral region exists only very near the blade surface. The radial component of turbulence intensity is higher than the streamwise component because of the effect of rotation.
Anand, A. K.↗
Computation of transonic boattail flow with separation
The relaxation procedure of South and Jameson for the full potential transonic flow equation was coupled to a modified Reshotko-Tucker integral boundary-layer technique with an empirical model for separated flow. The viscous and inviscid flows were solved iteratively until convergence was obtained. This iterative method was then applied to the subsonic and transonic flow over a series of axisymmetric circular-arc boattails with solid jet plume simulators. Comparisons of theoretical and experimental surface pressures and boattail drag are presented over a free-stream Mach numbers below 0.90. The qualitative variation of boattail drag with free-stream Mach number and boattail angle well into the region of transonic drag rise was correctly predicted; however, the absolute drag levels were significantly underpredicted. For separated flows, the empirical discriminating streamline model gives good results up to a free-stream Mach number of about 0.90 and allows reasonable predictions for shock-induced separation if the proper separation location and separation turning angle are known.
Wilmoth, R. G.↗
Generation of C-type cascade grids for viscous flow computation
A rapid procedure for generating C-type cascade grids suitable for viscous flow computations in turbomachinery blade rows is presented. The resulting mesh is periodic from one blade passage to the next, nearly orthogonal, and continuous across the wake downstream of a blade. The procedure employs a pair of conformal mappings that take the exterior of the cascade into the interior of an infinite strip with curved boundaries. The final transformation to a rectangular computational domain is accomplished numerically. The boundary values are obtained from a panel solution of an integral equation and the interior values by a rapid ADI solution of Laplace's equation. Examples of C-type grids are presented for both compressor and turbine blades and the extension of the procedure to three dimensions is briefly outlined.
Sockol, P. M.↗
Generalized multi-point inverse airfoil design
In a rather general sense, inverse airfoil design can be taken to mean the problem of specifying a desired set of airfoil characteristics, such as the airfoil maximum thickness ratio, pitching moment, part of the velocity distribution or boundary-layer development, etc., then from this information determine the corresponding airfoil shape. This paper presents a method which approaches the design problem from this perspective. In particular, the airfoil is divided into segments along which, together with the design conditions, either the velocity distribution or boundary-layer development may be prescribed. In addition to these local desired distributions, single parameters like the airfoil thickness can be specified. The problem of finding the airfoil shape is determined by coupling an incompressible, inviscid, inverse airfoil design method with a direct integral boundary-layer analysis method and solving the resulting nonlinear equations via a multidimensional Newton iteration technique. The approach is fast and easily allows for interactive design. It is also flexible and could be adapted to solving compressible, inverse airfoil design problems.
Selig, Michael S.↗
Theoretical study of the incompressible Navier-Stokes equations by the least-squares method
Usually the theoretical analysis of the Navier-Stokes equations is conducted via the Galerkin method which leads to difficult saddle-point problems. This paper demonstrates that the least-squares method is a useful alternative tool for the theoretical study of partial differential equations since it leads to minimization problems which can often be treated by an elementary technique. The principal part of the Navier-Stokes equations in the first-order velocity-pressure-vorticity formulation consists of two div-curl systems, so the three-dimensional div-curl system is thoroughly studied at first. By introducing a dummy variable and by using the least-squares method, this paper shows that the div-curl system is properly determined and elliptic, and has a unique solution. The same technique then is employed to prove that the Stokes equations are properly determined and elliptic, and that four boundary conditions on a fixed boundary are required for three-dimensional problems. This paper also shows that under four combinations of non-standard boundary conditions the solution of the Stokes equations is unique. This paper emphasizes the application of the least-squares method and the div-curl method to derive a high-order version of differential equations and additional boundary conditions. In this paper, an elementary method (integration by parts) is used to prove Friedrichs' inequalities related to the div and curl operators which play an essential role in the analysis.
Jiang, Bo-Nan↗
Integrating matrix solution of the hybrid state vector equations for beam vibration
A simple, versatile, and efficient computational technique has been developed for dynamic analysis of linear elastic beam and rod type of structures. Moreover, the method provides a rather general solution approach for two-point boundary value problems that are described by a single independent spatial variable. For structural problems, the method is implemented by a mixed state vector formulation of the differential equations, combined with an integrating matrix solution procedure. Highly accurate solutions are easily achieved with this approach. Example solutions are given for beam vibration problems including discontinuous stiffness and mass parameters, elastic restraint boundary conditions, concentrated inertia loading, and rigid body modes
Lehman, L. L.↗
Computations of internal turbulent flow with large separated flow regions
An implicit two-equation turbulence solver in generalized coordinates has been developed and is used in conjunction with the three-dimensional incompressible Navier-Stokes solver, INS3D, to calculate the internal flow in one channel and in an additional channel with a sudden 2:3 expansion. A new and consistent boundary procedure for a low Reynolds number form of the kappa-epsilon turbulence model is chosen to integrate the equations up to the wall. The high Reynolds number form of the equations are integrated using wall functions. The latter approach yields a faster convergence to the steady state solution than the former. For the case of channel flow, both the wall function and wall boundary condition approaches yield results in good agreement with the experimental data. the back-step (sudden expansion) flow is calculated using wall function approach. The predictions are in reasonable agreement with the experimental data.
Kaul, U. K.↗
Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules
There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.
Papadopoulos, Dimitrios↗
Approximate method of integration of laminar boundary layer in incompressible fluid
A method is given for the approximate solution of the equations of the two-dimensional laminar boundary layer in an incompressible fluid. The method is based on the use of a system of equations of successive moments that is easily solved for simple supplementary assumptions. The solution obtained is given in closed form by simple formulas and is claimed to be no less accurate than the complicated solutions previously obtained, which were based on the use of special classes of flows.
BOUNDARY LAYER - WING SECTIONS↗
Reduction of the Dirichlet problem to an initial value problem.
Although the derivation is concerned with solutions for plane regions with prescribed boundary values, the approach presented could by easily generalized to higher dimensions. The initial-value method is derived by a combination of invariant imbedding techniques and the Fredholm integral equation method of representation of the potential as a function of a dilayer distribution on the boundary of the region in question.
Kalaba, R.↗
Analytical evaluation of a surface integral expressing the coupling between interior and exterior volumes in a FE-IE approach
Recently an approach which combines the finite element technique and an integral equation to determine the fields scattered by inhomogeneous bodies of complicated shape has been proposed. Basically, a mathematical surface which encloses the scatterers is introduced, thus dividing the space into an interior and an exterior volume, in which the finite element technique and an integral equation for EM scattering, respectively, are applied. The integral equation is set up for the tangential components of the fields at the surface, while the interior volume the unknowns are the total fields. Continuity of the tangential fields at the boundary, as required by Maxwell's equations, is imposed, thus coupling the two methods to obtain a consistent solution. The coupling term is expressed by a surface integral formed by the dot product of a FE basis function and an IE testing function, or viceversa. By choosing the boundary to be a surface of revolution and by making a convenient selection of IE basis (testing) functions, it is possible to evaluate the integrals analytically on surfaces such as curved triangles, curved quadrilaterals and curved pentagons. We will illustrate the salient steps involved in setting up and carrying out these integrals and discuss what class of basis (testing) functions and analytic surfaces of revolution they are applicable to. Analytic calculations offer the advantage of better accuracy than purely numerical ones, and, when combined with them, often shed light on issues of numerical convergence and limiting values. Furthermore, they may reduce computation time and storage requirements.
Zuffada, C.↗
Turbulent atmospheric flow over a backward-facing step
The phenomenon of atmospheric shear layer separation over a man-made structure such as a building (modeled as a backward-facing step) has been analyzed theoretically by (1) solving the two-dimensional equations of motion in the two variables, stream function and vorticity, and by (2) employing an approximate integral technique. Boundary conditions for the undisturbed flow are that of the turbulent atmospheric shear flow over a rough terrain. In the first approach a two-equation model of turbulence was used. In the second approach an approximate technique was utilized in an attempt to describe the details of the flow in the recirculation zone behind the step. The results predict velocity profiles in sufficient detail that the presence of the corner eddy in the region of negative surface pressure gradient is evident. The magnitude of the reversed flow velocity in the recirculation eddy has been found to agree with that found from experiments. Also, a surface eddy viscosity distribution has been an outgrowth of the method which realistically follows the magnitude of the surface pressure gradient distribution as found experimentally.
Kaul, U. K.↗
Using EIGER for Antenna Design and Analysis
EIGER (Electromagnetic Interactions GenERalized) is a frequency-domain electromagnetics software package that is built upon a flexible framework, designed using object-oriented techniques. The analysis methods used include moment method solutions of integral equations, finite element solutions of partial differential equations, and combinations thereof. The framework design permits new analysis techniques (boundary conditions, Green#s functions, etc.) to be added to the software suite with a sensible effort. The code has been designed to execute (in serial or parallel) on a wide variety of platforms from Intel-based PCs and Unix-based workstations. Recently, new potential integration scheme s that avoid singularity extraction techniques have been added for integral equation analysis. These new integration schemes are required for facilitating the use of higher-order elements and basis functions. Higher-order elements are better able to model geometrical curvature using fewer elements than when using linear elements. Higher-order basis functions are beneficial for simulating structures with rapidly varying fields or currents. Results presented here will demonstrate curren t and future capabilities of EIGER with respect to analysis of installed antenna system performance in support of NASA#s mission of exploration. Examples include antenna coupling within an enclosed environment and antenna analysis on electrically large manned space vehicles.
Champagne, Nathan J.↗
Finite element analysis of unsteady incompressible flow around an oscillating obstacle of arbitrary shape.
An analytical procedure based on Navier-Stokes equations was developed for representing unsteady flow patterns around oscillating obstacles. A variational formulation of the Helmholtz vorticity equation was discretized in finite element form and integrated numerically. At each step of the numerical integration the velocity field around the obstacle was determined from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints at each time step of the numerical integration. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.
Bratanow, T.↗
Feasibility study of the numerical integration of shell equations using the field method
The field method is developed for arbitrary open branch domains subjected to general linear boundary conditions. Although closed branches are within the scope of the method, they are not treated here. The numerical feasibility of the method has been demonstrated by implementing it in a computer program for the linear static analysis of open branch shells of revolution under asymmetric loads. For such problems the field method eliminates the well-known numerical problem of long subintervals associated with the rapid growth of extraneous solutions. Also, the method appears to execute significantly faster than other numerical integration methods.
Cohen, G. A.↗
Application of the method of integral relations to the calculation of incompressible turbulent boundary layers
Integral relations method efficiency and accuracy in solving incompressible turbulent boundary layer equations, noting acceptability for turbulent flows with small wake components
Murphy, J. D.↗
Numerical solutions in three dimensional elastostatics
Three dimensional elastostatics problems solution method using singular integral equations solvable numerically for surface tractions and displacements of fully mixed boundary value problem
Cruse, T. A.↗