SEARCH · Engineering Papers

Results for “Smoothed boundary method”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Effect of leading-edge geometry on boundary-layer receptivity to freestream sound

The receptivity to freestream sound of the laminar boundary layer over a semi-infinite flat plate with an elliptic leading edge is simulated numerically. The incompressible flow past the flat plate is computed by solving the full Navier-Stokes equations in general curvilinear coordinates. A finite-difference method which is second-order accurate in space and time is used. Spatial and temporal developments of the Tollmien-Schlichting wave in the boundary layer, due to small-amplitude time-harmonic oscillations of the freestream velocity that closely simulate a sound wave travelling parallel to the plate, are observed. The effect of leading-edge curvature is studied by varying the aspect ratio of the ellipse. The boundary layer over the flat plate with a sharper leading edge is found to be less receptive. The relative contribution of the discontinuity in curvature at the ellipse-flat-plate juncture to receptivity is investigated by smoothing the juncture with a polynomial. Continuous curvature leads to less receptivity. A new geometry of the leading edge, a modified super ellipse, which provides continuous curvature at the juncture with the flat plate, is used to study the effect of continuous curvature and inherent pressure gradient on receptivity.

Lin, Nay↗

NASA NTRS2018

Scale-Resolving Simulations of Low-Pressure Turbine Cascades with Wall Roughness Using A Spectral-Element Method

The accurate prediction of wall-roughness effects in turbomachinery is becoming critical as turbine designers address airfoil surface quality and degradation concerns arising from the shift to advanced ceramic matrix composite (CMC) or additively-manufactured airfoils operating in higher temperature environments. In this paper, a recently developed computational capability for accurate and efficient scale-resolving simulations of turbomachinery is extended to analyze the boundary- layer separation and transition characteristics in a rough-wall low-pressure turbine (LPT) cascade. The computational capability is based on an entropy-stable discontinuous-Galerkin spectral-element approach that extends to arbitrarily high orders of spatial and temporal accuracy, and is implemented in an efficient manner for a modern high performance computer architecture. Results from the scale-resolving simulations of both smooth and rough airfoil cascades are presented and compared to previous experiments and numerical simulations. The results show that the suction surface boundary layer undergoes laminar separation, transition, and turbulent reattachment for the smooth airfoil cascade, while in the presence of roughness the separation and transition behavior of the suction surface boundary layer is substantially modified. The differences between the smooth and rough airfoil cascades are then highlighted by a detailed analysis of their respective turbulent flow fields.

Spectral-Element↗

NASA NTRS1988

Entrance and exit region friction factor models for annular seal analysis

The Mach number definition and boundary conditions in Nelson's nominally-centered, annular gas seal analysis are revised. A method is described for determining the wall shear stress characteristics of an annular gas seal experimentally. Two friction factor models are developed for annular seal analysis; one model is based on flat-plate flow theory; the other uses empirical entrance and exit region friction factors. The friction factor predictions of the models are compared to experimental results. Each friction model is used in an annular gas seal analysis. The seal characteristics predicted by the two seal analyses are compared to experimental results and to the predictions of Nelson's analysis. The comparisons are for smooth-rotor seals with smooth and honeycomb stators. The comparisons show that the analysis which uses empirical entrance and exit region shear stress models predicts the static and stability characteristics of annular gas seals better than the other analyses. The analyses predict direct stiffness poorly.

Elrod, David Alan↗

NASA NTRS1982

Feasibility of SPF/DB titanium sandwich for LFC wings

The feasibility of fabricating SPF/DB titanium structures of sufficient smoothness to be used for laminar flow wing surfaces was demonstrated. Two methods of fabricating panels which meet the surface smoothness criteria were demonstrated. The first consists of superplasticity forming/diffusion bonding a panel using steel dies, then machining the surface to the required flatness and finish after forming. This approach, however, was estimated to be more costly than the second approach, in which the panel is formed against ceramic platens which produce the desired surface smoothness without subsequent finishing. Four panels were fabricated, and their surface quality was evaluated. The specified criteria for waviness compared with actual measured values on the panels is shown. The acceptable surface quality as well as feasibility of the laminar flow control surface design, in which separate strips incorporating the boundary layer bleed provisions are bonded into slots on the surface, was also demonstrated.

Mcquilkin, F. T.↗

NASA NTRS2014

Toward Automatic Verification of Goal-Oriented Flow Simulations

We demonstrate the power of adaptive mesh refinement with adjoint-based error estimates in verification of simulations governed by the steady Euler equations. The flow equations are discretized using a finite volume scheme on a Cartesian mesh with cut cells at the wall boundaries. The discretization error in selected simulation outputs is estimated using the method of adjoint-weighted residuals. Practical aspects of the implementation are emphasized, particularly in the formulation of the refinement criterion and the mesh adaptation strategy. Following a thorough code verification example, we demonstrate simulation verification of two- and three-dimensional problems. These involve an airfoil performance database, a pressure signature of a body in supersonic flow and a launch abort with strong jet interactions. The results show reliable estimates and automatic control of discretization error in all simulations at an affordable computational cost. Moreover, the approach remains effective even when theoretical assumptions, e.g., steady-state and solution smoothness, are relaxed.

Simulations↗

NASA NTRS1984

Surface grid generation for wing-fuselage bodies

In the application of finite-difference methods to obtain numerical solutions of viscous compressible fluid flow about wing-fuselage bodies, it is advantageous to transform the governing equations to an idealized boundary-fitted coordinate system. The advantages are reduced computational complexity and added accuracy in the application of boundary conditions. The solution process requires that a grid be superimposed on the physical solution domain which corresponds to a uniform grid on a rectangular computational domain (uniform rectangular parallel-epiped). Grid generation is the determination of a one to one relationship between grid points in the physical domain and grid points in the computational domain. A technique for computing wing-fuselage surface grids using the Harris geometry and software for smooth-surface representation is described. Grid spacing control concepts which govern the relationship between the wing-fuselage surface and the computational grid are also presented.

Smith, R. E.↗

NASA NTRS1985

Minimal norm constrained interpolation

In computational fluid dynamics and in CAD/CAM, a physical boundary is usually known only discreetly and most often must be approximated. An acceptable approximation preserves the salient features of the data such as convexity and concavity. In this dissertation, a smooth interpolant which is locally concave where the data are concave and is locally convex where the data are convex is described. The interpolant is found by posing and solving a minimization problem whose solution is a piecewise cubic polynomial. The problem is solved indirectly by using the Peano Kernal theorem to recast it into an equivalent minimization problem having the second derivative of the interpolant as the solution. This approach leads to the solution of a nonlinear system of equations. It is shown that Newton's method is an exceptionally attractive and efficient method for solving the nonlinear system of equations. Examples of shape-preserving interpolants, as well as convergence results obtained by using Newton's method are also shown. A FORTRAN program to compute these interpolants is listed. The problem of computing the interpolant of minimal norm from a convex cone in a normal dual space is also discussed. An extension of de Boor's work on minimal norm unconstrained interpolation is presented.

Irvine, L. D.↗

NASA NTRS2015

A Numerical and Experimental Study of Coflow Laminar Diffusion Flames: Effects of Gravity and Inlet Velocity

In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both computationally and experimentally. A series of flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity conditions with the fuel stream CH4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth vorticity-velocity formulation of the governing equations was employed to describe the reactive gaseous mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional computational domain was employed, where the grid on the axisymmetric domain was a nonuniform tensor product mesh. The governing equations and boundary conditions were discretized on the mesh by a nine-point finite difference stencil, with the convective terms approximated by a monotonic upwind scheme and all other derivatives approximated by centered differences. The resulting set of fully coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton's method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off height, and soot temperature were determined by flame emission images recorded by a digital camera, and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first compared to assess the accuracy of the numerical model. After its validity was established, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of laminar coflow methane-air diffusion flames was explored further by examining quantities derived from the computational results.

microgravity↗

NASA NTRS1987

Measurement of the effect of manufacturing deviations on natural laminar flow for a single engine general aviation airplane

Renewed interest in natural laminar flow (NLF) had rekindled designer concern that manufacuring deviations may destroy the effectiveness of NLF for an operational aircraft. Experiments are summarized that attemtped to measure total drag changes associated with three different wing surface conditions on an aircraft typical of current general aviation high performance singles. The speed power technique was first used in an attempt to quantify the changes in total drag. Predicted and measured boundary layer transition locations for three different wing surface conditions were also compared, using two different forms of flow visualization. The three flight test phases included: assessment of an unpainted airframe, flight tests of the same aircraft after painstakingly filling and sanding the wings to design contours, and similar measurement after this aricraft was painted. In each flight phase, transition locations were monitored using with sublimating chemicals or pigmented oil. Two-dimensional drag coefficients were estimated using the Eppler-Somers code and measured with a wake rake in a method very similar to Jones' pitot traverse method. The net change in two-dimensional drag coefficient was approximately 20 counts between the unpainted aircraft and the hand-smoothed aircraft for typical cruise flight conditions.

Source record↗

NASA NTRS1994

Numerical study of the effects of icing on viscous flow over wings

An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved using a diagonal form of an alternating-direction implicit (ADI) approximate factorization procedure. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic-based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are less than 60 percent of that required for a full-blown Navier-Stokes analysis for steady flow applications and about 60 percent of the Navier-Stokes CPU times for unsteady flows in non-vector processing machines. Applications of the method are presented for a rectangular NACA 0012 wing in low subsonic steady flow at moderate and high angles of attack, and for an F-5 wing in steady and unsteady subsonic and transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to Navier-Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier-Stokes equations can be retained with a significant savings in computational time.

Sankar, L. N.↗

NASA NTRS1987

New main reflector, subreflector and dual chamber concepts for compact range applications

A compact range is a facility used for the measurement of antenna radiation and target scattering problems. Most presently available parabolic reflectors do not produce ideal uniform plane waves in the target zone. Design improvements are suggested to reduce the amplitude taper, ripple and cross polarization errors. The ripple caused by diffractions from the reflector edges can be reduced by adding blended rolled edges and shaping the edge contour. Since the reflected edge continues smoothly from the parabola onto the rolled surface, rather than being abruptly terminated, the discontinuity in the reflected field is reduced which results in weaker diffracted fields. This is done by blending the rolled edges from the parabola into an ellipse. An algorithm which enables one to design optimum blended rolled edges was developed that is based on an analysis of the continuity of the surface radius of curvature and its derivatives across the junction. Futhermore, a concave edge contour results in a divergent diffracted ray pattern and hence less stray energy in the target zone. Design equations for three-dimensional reflectors are given. Various examples were analyzed using a new physical optics method which eliminates the effects of the false scattering centers on the incident shadow boundaries. A Gregorian subreflector system, in which both the subreflector and feed axes are tilted, results in a substantial reduction in the amplitude taper and cross polarization errors. A dual chamber configuration is proposed to eliminate the effects of diffraction from the subreflector and spillover from the feed. A computationally efficient technique, based on ray tracing and aperture integration, was developed to analyze the scattering from a lossy dielectric slab with a wedge termination.

Pistorius, C. W. I.↗

NASA NTRS2015

Inward-Turning Streamline-Traced Inlet Design Method for Low-Boom, Low-Drag Applications

A new design method for inward-turning, streamline-traced inlets is presented. Resulting designs are intended for moderate supersonic, low-drag, low-boom applications such as that required for NASA's proposed low-boom flight demonstration aircraft. A critical feature of these designs is the internal cowl lip angle that allows for little or no flow turning on the outer nacelle. Present methods using conical-flow Busemann parent flowfields have simply truncated, or otherwise modified the stream-traced contours to include this internal cowl angle. Such modifications disrupt the parent flowfield, reducing inlet performance and flow uniformity. The method presented herein merges a conical flowfield that includes a leading shock with a truncated Busemann flowfield in a manner that minimizes unwanted interactions. A leading internal cowl angle is now inherent in the parent flowfield, and inlet contours traced from this flowfield retain its high performance and good flow uniformity. CFD analysis of a candidate inlet design is presented that verifies the design technique, and reveals a starting issue with the basic geometry. A minor modification to the cowl lip region is shown to eliminate this phenomenon, thereby allowing starting and smooth transition to sub-critical operation as back-pressure is increased. An inlet critical-point total pressure recovery of 96 is achieved based on CFD results for a Mach 1.7 freestream design. Correction for boundary-layer displacement thickness, and sizing for a given engine airflow requirement are also discussed.

low sonic boom↗

NASA NTRS1981

Quasi-Lagrangian models of nascent thermals

Slab-symmetric and axisymmetric geometries are considered in the modeling of thermals by considering their boundaries as self-advecting vortex sheets, for the case of a neutral environment without turbulent processes. The initial accelerations for circular cross section thermals are obtained analytically, and a numerical method is developed to simulate the evolution of thermals by time-marching in which the vortex sheet is divided into finite segments whose positions are tracked in a quasi-Lagrangian fashion. Self-advection is considered, and a redistribution procedure is adopted to prevent the segments from bunching unmanageably. The induced field of motion is determined both inside and outside the thermal. It is shown that (1) the axisymmetric thermal rises more quickly than the slab-symmetric thermal, (2) vorticity maxima occur both within the concavities and on the trailing edges, and (3) the leading edges are very smooth.

Rambaldi, S.↗

NASA NTRS1985

Analysis of a finite composite plate with smooth rigid pin

An elastic continuum method of analysis developed for a smooth rigid pin in a finite composite plate by using complex potentials that satisfy the governing orthotropic differential equation is presented. The plate is subjected to uniaxial loading and the nonlinear problems of partial contact in the pin-plate are studied. Three types of pin fits are examined: interference, push, and clearance. The problems of stresses due to pure interference, load contact behavior, and stresses variations are addressed. The effects of the rectangular geometry and the orthotropic property of the plate are demonstrated through numerical results. It is found that the use of large hole diameter/plate width ratios have deleterious effects and the use of interference has beneficial effects. The solutions obtained will help in the development of special finite elements for fastener joints. Finally, the boundary conditions on the finite outer boundaries are found to be satisfied by a successive integration technique.

Naidu, A. C. B.↗

NASA NTRS1975

Remarks on the theory of aerodynamic noise

The accuracy of approximations employed in the Lighthill theory in current aerodynamic noise research is critically evaluated. Based on the method of matched asymptotic expansions, the full Navier-Stokes equations are expanded for small Mach numbers. To the first order, the near-field is generally nonisentropic. The pressure field (pseudo-sound) is generated by the incompressible Reynolds stresses in the turbulent flow and the velocity, pressure perturbation, and their derivatives on the boundaries. In the far-field, the first-order pressure (acoustic) field satisfying a linear wave equation is obtained by matching with the pseudo-sound field. A uniformly valid solution for the pressure field in a stationary or uniformly moving medium is obtained. The solution shows that the generation of the first-order aerodynamic noise does not depend on the viscous, thermal, or entropy effects in the adiabatic flow nor on the shear stress on a smooth rigid boundary.

Pan, Y. S.↗

NASA NTRS1984

Polycrystalline Silicon Sheets for Solar Cells by the Improved Spinning Method

Cost reduction of silicon materials in the photovoltaic program of materials was examined. The current process of producing silicon sheets is based entirely on the conventional Czochralski ingot growth and wafering used in the semiconductor industry. The current technology cannot meet the cost reduction demands for producing low cost silicon sheets. Alternative sheet production processes such as unconventional crystallization are needed. The production of polycrystalline silicon sheets by unconventional ingot technology is the casting technique. Though large grain sheets were obtained by this technique, silicon ribbon growth overcomes deficiencies of the casting process by obtaining the sheet directly from the melt. The need to solve difficulties of growth stability and impurity effects are examined. The direct formation process of polycrystalline silicon sheets with large grain size, smooth surface, and sharp edges from the melt with a high growth rate which will yield low cost silicon sheets for solar cells and the photovoltaic characteristics associated with this type of sheet to include an EBIC study of the grain boundaries are described.

Maeda, Y.↗

NASA NTRS1985

Comparison of predicted and measured drag for a single-engine airplane

Renewed interest in natural laminar flow (NLF) has rekindled designers' concerns that manufacturing deviations, (loss of surface contours or other surface imperfections) may destroy the effectiveness of NLF for an operational airplane. This paper reports on experimental research that compares predicted and measured boundary layer transition, total drag, and two-dimensional drag coefficients for three different wing surface conditions on an airplane typical of general aviation manufacturing technology. The three flight test phases included: (1) assessment of an unpainted airframe, (2) flight tests of the same airplane after painstakingly filling and sanding the wings to design contours, and (3) similar measurements after this airplane was painted. In each flight phase, transition locations were monitored using either sublimating chemicals or pigmented oil. As expected, total drag changes were difficult to measure. Two-dimensional drag coefficients were estimated using the Eppler-Somers code and measured with a wake rake in a method very similar to Jones' pitot traverse method. The net change in two-dimensional drag was approximately 20 counts between the unpainted airplane and the 'hand-smoothed' airplane for typical cruise flight conditions.

Ward, D. T.↗

NASA NTRS1991

The dependence of wind stress on wave height and wind speed

Three near-neutral boundary layer data sets were investigated with the aim of finding a dependence of wind stress on both wind speed and significant wave height. The data set most representative of open-ocean wave height, wind speed, and momentum flux conditions, was selected and analyzed by means of the least-squares method to produce a new parameterization for the wind stress as a function of both wind speed and significant wave height. This study shows that the wind stress, and consequently the drag coefficient, decreases with increasing wave height for a fixed wind speed. The study also shows that the curvature of the wind profile decreases with increasing wave height and that the C(DN) = A + BU-bar(10) form for the drag coefficient parameterization is inadequate. A drag coefficient that applies to both smooth and rough flows is proposed. These results are more applicable for open-ocean deep-water conditions and less applicable for sheltered, closed, shallow water sites.

Blake, Reginald A.↗