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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗

NASA NTRS2018

Characterization of Titanium Alloys Produced by Electron Beam Directed Energy Deposition

Functionally graded materials offer the potential to improve structural efficiency by allowing the material composition and/or microstructural features to spatially vary within a component. Additive manufacturing techniques enable the fabrication of such graded materials and structures. While examining several titanium alloys, this paper focuses on Ti-8Al-1Er as it has a unique microstructure that is only feasible when produced by rapid solidification methods like electron beam directed energy deposition, an additive manufacturing process. The results show that, when mixed, Ti-8Al-1Er and commercially-pure titanium uniformly mix at various ratios and the resultant static tensile properties of the mixed alloys behave according to rule-of-mixtures. At discontinuous interfaces between Ti-8Al-1Er and commercially-pure titanium, the crack growth behavior progresses smoothly across the discontinuity as the crack transitions from one crack growth regime into another. Studies on monolithic samples shows the mechanisms of damage in the Ti-8Al-1Er; specifically, that strain localization occurs near grain boundaries of high mis-orientation on the microscale and that twinning and dislocation density is concentrated near erbia-strengthening particles (Er2O3) on the nanoscale.

Newman, John A.↗

NASA NTRS1995

A 3D High-Order Unstructured Finite-Volume Algorithm for Solving Maxwell's Equations

A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Liu, Yen↗

NASA NTRS1980

Supersonic flow development in slotted wind tunnels

The development of test section slot shapes for achieving smooth supersonic Mach number distribution without overexpansion or waviness has, in the past, been largely an experimentally iterative or 'cutand try' procedure for each wind tunnel. To overcome the obvious disadvantages of time and expense involved in such an experimental approach, a simple analytical method has been developed to predict the supersonic flow development in a two-dimensional slotted tunnel given only the variation of open area ratio with downstream distance and the Mach number corresponding to the plenum static pressure. The well known method of characteristics is used with the constraint that it be compatible with the quadratic cross-flow pressure drop boundary condition at the slotted wall. The predicted results from this method agree remarkably well with the experimental calibration data available for some of the existing facilities. The flow mechanism responsible for causing overexpansion in the centerline Mach number distribution with some slot shapes has been brought to light.

Ramaswamy, M. A.↗

NASA NTRS2015

High-Order Hyperbolic Residual-Distribution Schemes on Arbitrary Triangular Grids

In this paper, we construct high-order hyperbolic residual-distribution schemes for general advection-diffusion problems on arbitrary triangular grids. We demonstrate that the second-order accuracy of the hyperbolic schemes can be greatly improved by requiring the scheme to preserve exact quadratic solutions. We also show that the improved second-order scheme can be easily extended to third-order by further requiring the exactness for cubic solutions. We construct these schemes based on the LDA and the SUPG methodology formulated in the framework of the residual-distribution method. For both second- and third-order-schemes, we construct a fully implicit solver by the exact residual Jacobian of the second-order scheme, and demonstrate rapid convergence of 10-15 iterations to reduce the residuals by 10 orders of magnitude. We demonstrate also that these schemes can be constructed based on a separate treatment of the advective and diffusive terms, which paves the way for the construction of hyperbolic residual-distribution schemes for the compressible Navier-Stokes equations. Numerical results show that these schemes produce exceptionally accurate and smooth solution gradients on highly skewed and anisotropic triangular grids, including curved boundary problems, using linear elements. We also present Fourier analysis performed on the constructed linear system and show that an under-relaxation parameter is needed for stabilization of Gauss-Seidel relaxation.

Mazaheri, Alireza↗

NASA NTRS1985

Numerical approach for the aerodynamic analysis if airfoils with laminar separation

A numerical method for simultaneously and efficiently coupling an external subsonic potential flow and an interior viscous flow such that the two flows match at an interfacing boundary is discussed. Both a panel method and a simple point compressible vortex model are used for the outer potential field. The interior flow solvers which were used are the Navier-Stokes and Euler codes of T. J. Coakley and the Euler code of A. Verhoff. In order to test compatibility, the panel method is coupled to the less expensive Euler codes since the coupling procedure is identical with the Navier-Stokes code. The results show significant efficiency improvements can be obtained over the uncoupled approach. Results also indicate the outer potential flow is best represented by the simple point compressible vortex model. The panel method couples smoothly to Coakley's implicit code but is numerically incompatible as coupled with the explicit Euler code. An improved Navier-Stokes code is under initial development which extends the Euler code to include the necessary viscous terms. Results are shown for all infinite length channel with one wavy periodic wall with and without laminar separation.

Halt, D. W.↗

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↗