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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Rough-Wall Channel Analysis Using Suboptimal Control Theory

The original aim of this work was to shed some light on the physics of turbulence over rough walls using large-eddy simulations and the suboptimal-control wall boundary conditions introduced by Nicoud et al. It was hoped that, if that algorithm was used to fit the mean velocity profile of the simulations to that of a rough-walled channel, instead of to a smooth one, the wall stresses introduced by the control algorithm would give some indication of what aspects of rough walls are most responsible for the modification of the flow in real turbulence. It was similarly expected that the structure of the resulting velocity fluctuations would share some of the characteristics of rough-walled flows, thus again suggesting what is intrinsic and what is accidental in the effect of geometric wall roughness. A secondary goal was to study the effect of 'unphysical' boundary conditions on the outside flow by observing how a relatively major change of the target velocity profile, and therefore presumably of the applied wall stresses, modifies properties such as the dominant length scales of the velocity fluctuations away from the wall. As will be seen below, this secondary goal grew more important during the course of the study, which was carried out during a short summer visit of the first two authors to the CTR. It became clear that there are open questions about the way in which the control algorithm models the boundary conditions, even for smooth walls, and that these questions make the physical interpretation of the results difficult. Considerable more work in that area seems to be needed before even relatively advanced large-eddy simulations, such as these, can be used to draw conclusions about the physics of wall-bounded turbulent flows. The numerical method is the same as in Nicoud et al. The modifications introduced in the original code are briefly described in section 2, but the original paper should be consulted for a full description of the algorithm. The results are presented in section 3 and summarized in section 4. The elementary properties of turbulence over rough walls which are used in the text have been taken from recent reviews such as Raupach et al. or Jimenez.

Flores, O.↗

NASA NTRS2012

Images of the Extended Outer Regions of the Debris Ring around HR 4796 A

We present high-contrast images of HR 4796 A taken with Subaru/HiCIAO in H-band, resolving the debris disk in scattered light. The application of specialized angular differential imaging methods (ADI) allows us to trace the inner edge of the disk with high precision, and reveals a pair of "streamers" extending radially outwards from the ansae. Using a simple disk model with a power-law surface brightness profile, we demonstrate that the observed streamers can be understood as part of the smoothly tapered outer boundary of the debris disk, which is most visible at the ansae. Our observations are consistent with the expected result of a narrow planetesimal ring being ground up in a collisional cascade, yielding dust with a wide range of grain sizes. Radiation forces leave large grains in the ring and push smaller grains onto elliptical, or even hyperbolic trajectories. We measure and characterize the disk's surface brightness profile, and confirm the previously suspected offset of the disk's center from the star's position along the ring's major axis. Furthermore, we present first evidence for an offset along the minor axis. Such offsets are commonly viewed as signposts for the presence of unseen planets within a disk's cavity. Our images also offer new constraints on the presence of companions down to the planetary mass regime (approx 9 M(sub Jup) at 0".5, approx 3 M(sub Jup) at 1").

Thalmann, C.↗

NASA NTRS1987

Frequency-shaped large-angle maneuvers

The paper considers the problem of maneuvering a flexible spacecraft through large angles in finite time. The basic control problem is divided into two parts. The first part consists of generating a frequency-shaped open-loop solution for the nonlinear rigid body as the nominal solution. The resulting two-point boundary-value problem is solved by introducing a continuation method for altering the mass distribution and boundary conditions for the spacecraft. For the second part, a feedback control is designed by linearizing the flexible body response about several points along the rigid body nominal solution. The perturbation gains are designed by using a frequency-shaped cost functional approach. The gains are linearly interpolated to produce smooth control time-histories as the linear piecewise constant plant models change during the maneuver.

Chun, Hon M.↗

NASA NTRS1986

Wall turbulence control

A variety of wall turbulence control devices which were experimentally investigated are discussed; these include devices for burst control, alteration of outer flow structures, large eddy substitution, increased heat transfer efficiency, and reduction of wall pressure fluctuations. Control of pre-burst flow was demonstrated with a single, traveling surface depression which is phase-locked to elements of the burst production process. Another approach to wall turbulence control is to interfere with the outer layer coherent structures. A device in the outer part of a boundary layer was shown to suppress turbulence and reduce drag by opposing both the mean and unsteady vorticity in the boundary layer. Large eddy substitution is a method in which streamline curvature is introduced into the boundary layer in the form of streamwise vortices. Riblets, which were already shown to reduce turbulent drag, were also shown to exhibit superior heat transfer characteristics. Heat transfer efficiency as measured by the Reynolds Analogy Factor was shown to be as much as 36 percent greater than a smooth flat plate in a turbulent boundary layer. Large Eddy Break-Up (LEBU) which are also known to reduce turbulent drag were shown to reduce turbulent wall pressure fluctuation.

Wilkinson, Stephen P.↗

NASA NTRS2014

An Experimental and Computational Study on Soot Formation in a Coflow Jet Flame Under Microgravity and Normal Gravity

Upon the completion of the Structure and Liftoff in Combustion Experiment (SLICE) in March 2012, a comprehensive and unique set of microgravity coflow diffusion flame data was obtained. This data covers a range of conditions from weak flames near extinction to strong, highly sooting flames, and enabled the study of gravitational effects on phenomena such as liftoff, blowout and soot formation. The microgravity experiment was carried out in the Microgravity Science Glovebox (MSG) on board the International Space Station (ISS), while the normal gravity experiment was performed at Yale utilizing a copy of the flight hardware. Computational simulations of microgravity and normal gravity flames were also carried out to facilitate understanding of the experimental observations. This paper focuses on the different sooting behaviors of CH4 coflow jet flames in microgravity and normal gravity. The unique set of data serves as an excellent test case for developing more accurate computational models.Experimentally, the flame shape and size, lift-off height, and soot temperature were determined from line-of-sight flame emission images taken with a color digital camera. Soot volume fraction was determined by performing an absolute light calibration using the incandescence from a flame-heated thermocouple. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the chemically reacting flow, and the soot evolution was modeled by the sectional aerosol equations. The governing equations and boundary conditions were discretized on an axisymmetric computational domain by finite differences, and the resulting system of fully coupled, highly nonlinear equations was solved by a damped, modified Newtons method. The microgravity sooting flames were found to have lower soot temperatures and higher volume fraction than their normal gravity counterparts. The soot distribution tends to shift from the centerline of the flame to the wings from normal gravity to microgravity.

microgravity↗

NASA NTRS2014

Applying Workspace Limitations in a Velocity-Controlled Robotic Mechanism

A robotic system includes a robotic mechanism responsive to velocity control signals, and a permissible workspace defined by a convex-polygon boundary. A host machine determines a position of a reference point on the mechanism with respect to the boundary, and includes an algorithm for enforcing the boundary by automatically shaping the velocity control signals as a function of the position, thereby providing smooth and unperturbed operation of the mechanism along the edges and corners of the boundary. The algorithm is suited for application with higher speeds and/or external forces. A host machine includes an algorithm for enforcing the boundary by shaping the velocity control signals as a function of the reference point position, and a hardware module for executing the algorithm. A method for enforcing the convex-polygon boundary is also provided that shapes a velocity control signal via a host machine as a function of the reference point position.

Abdallah, Muhammad E.↗

NASA NTRS1990

Multigrid for hypersonic inviscid flows

The use of multigrid methods to solve the Euler equations for hypersonic flow is discussed. The steady state equations are considered with a Runge-Kutta smoother based on the time accurate equations together with local time stepping and residual smoothing. The effect of the Runge-Kutta coefficients on the convergence rate was examined considering both damping characteristics and convection properties. The importance of boundary conditions on the convergence rate for hypersonic flow is discussed. Also of importance are the switch between the second and fourth difference viscosity. Solutions are given for flow around the bump in a channel and flow around a biconic section.

Decker, Naomi H.↗

NASA NTRS1991

A multiple-block multigrid method for the solution of the three-dimensional Euler and Navier-Stokes equations

A multiple block multigrid method for the solution of the three dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell vertex method which employs central difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques and variable coefficient numerical dissipation results in an efficient and robust scheme is discussed. The multiblock strategy places the block loop within the Runge-Kutta Loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one and two block calculations in which the two block grid is generated by splitting the one block grid. Results are presented for both Euler and Navier-Stokes computations of wing/fuselage combinations.

Atkins, Harold↗

NASA NTRS2008

Aeroheating Test of CEV Entry Vehicle at Turbulent Conditions

An investigation of the aeroheating environment of the Project Orion Crew Entry Vehicle has been performed in the Arnold Engineering Development Center Tunnel 9. Data were measured on a approx. 3.5% scale model (0.1778m/7-inch diam.) of the vehicle using coaxial thermocouples in the Mach 8 and Mach 10 nozzles of Tunnel 9. Runs were performed at free stream Reynolds numbers of 1 106/ft to 20 10(exp 6)/ft in the Mach 10 nozzle and 8 10(exp 6)/ft to 48 10(exp 6)/ft in the Mach 8 nozzle. The test gas in Tunnel 9 is pure N2, which at these operating conditions remains un-dissociated and may be treated as a perfect gas. At these conditions, laminar, transitional, and turbulent flow was produced on the model at Mach 10, and transitional and turbulent conditions were produced on the model at Mach 8. The majority of runs were made on a clean, smooth-surface model configuration and a limited number of runs were made in which inserts with varying boundary-layer trips configurations were used to force the occurrence of transition. Laminar and turbulent predictions were generated for all wind tunnel test conditions and comparisons were performed with the data for the purpose of helping to define uncertainty margins for the computational method. Data from both the wind tunnel test and the computations are presented herein. Figure 1 shows a schematic of the thermocouple locations on the model and figures 2 and 3 show a photo and schematic of the AEDC Hypervelocity Tunnel 9. Figure 4 shows a typical grid used in the computations. From the comparisons shown in figures 5 through 8 it was concluded that for perfect-gas conditions, the computations could predict either fully-laminar or full-turbulent flow to within +/-10% of the experimental data. The experimental data showed that transition began on the leeside of the heatshield at a free stream Reynolds number of 9 10(exp 6)/ft in the Mach 10 nozzle and fully-developed turbulent flow was produced at 20 10(exp 6)/ft. In the Mach 8 nozzle, transition on the leeside of the heat-shield was observed for all test conditions, and full-developed turbulent flow occurred at a free stream Reynolds number of 18 10(exp 6)/ft. On the aftbody of the vehicle no evidence of turbulence was detected at Mach 10 conditions, and at Mach 8 conditions, transition appeared to begin on the windside of the aftbody at free stream Reynolds number of 18x10(exp 6)/ft with fully-developed turbulent flow occurring only at the highest test condition of 48x10(exp 6)/ft.

Hollis, Brian R.↗

NASA NTRS1991

A multi-block multigrid method for the solution of the Euler and Navier-Stokes equations for three-dimensional flows

A multi-block multigrid method for the solution of the three-dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell-vertex method which employs central-difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques, and variable-coefficient numerical smoothing results in an efficient and robust scheme. The multi-block strategy places the block loop within the Runge-Kutta loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one- and two-block calculations in which the two-block grid is generated by splitting the one-block grid. Results are presented for both Euler and Navier-Stokes computations of wings and wing-fuselage combinations.

Atkins, H. L.↗

NASA NTRS2014

A Computational and Experimental Study of Coflow Laminar Methane/Air Diffusion Flames: Effects of Fuel Dilution, Inlet Velocity, and Gravity

The influences of fuel dilution, inlet velocity, and gravity on the shape and structure of laminar coflow CH4-air diffusion flames were investigated computationally and experimentally. A series of nitrogen-diluted flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) on board the International Space Station was assessed numerically under microgravity (mu g) and normal gravity (1g) conditions with CH4 mole fraction ranging from 0.4 to 1.0 and average inlet velocity ranging from 23 to 90 cm/s. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the reactive gaseous mixture, and soot evolution was modeled by sectional aerosol equations. The governing equations and boundary conditions were discretized on a two-dimensional computational domain by finite differences, and the resulting set of fully coupled, strongly nonlinear equations was solved simultaneously at all points using a damped, modified Newton's method. Experimentally, flame shape and soot temperature were determined by flame emission images recorded by a digital color camera. Very good agreement between computation and measurement was obtained, and the conclusions were as follows. (1) Buoyant and nonbuoyant luminous flame lengths are proportional to the mass flow rate of the fuel mixture; computed and measured nonbuoyant flames are noticeably longer than their 1g counterparts; the effect of fuel dilution on flame shape (i.e., flame length and flame radius) is negligible when the flame shape is normalized by the methane flow rate. (2) Buoyancy-induced reduction of the flame radius through radially inward convection near the flame front is demonstrated. (3) Buoyant and nonbuoyant flame structure is mainly controlled by the fuel mass flow rate, and the effects from fuel dilution and inlet velocity are secondary.

flames↗

NASA NTRSDate not supplied

Evaluation of Near Singular Integrals for Computational Electromagnetics by Dimensionality Reduction

With the need for ever faster codes, a limiting factor that must be dealt with is the accurate yet efficient evaluation of interaction integrals between the more problematic near-field elements. Several recent works have together shown that all evaluations of source potential integrals and their derivatives for the most common bases and elements can be reduced to the evaluation of boundary line integrals; these can be evaluated by Gauss-Legendre quadrature, though integrand-smoothing transforms are often needed to accelerate their computation. In this paper, we modify the reported approach to eliminate cancellation errors in the line integral integrand, reinterpret the integral as a vertex function, and study the scalar potential integral form under the sinh transform and static subtraction acceleration methods.

D R Wilton↗

NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Olsson & Oliger, Gerritsen & Olsson, and Yee et al. (15, 7, 37) is revisited. The objective of this paper is to prove for the first time that the entropy split scheme is an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for nonlinear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. The proof is to replace the spatial derivatives by summation-by-parts (SBP) difference operators in the entropy split form of the equations using the physical entropy of the Euler equations. The numerical boundary closure follows directly from the SBP operator. No additional numerical boundary procedure is required. In contrast, Tadmor-type entropy conserving schemes (31) using mathematical entropies and more recently in (35], do not naturally come with a numerical boundary closure and a generalized SBP operator has to be developed (18). Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grubber splitting.

Sjogreen, Bjorn↗

NASA NTRS2006

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗

NASA NTRS1980

Numerical generation of two-dimensional orthogonal curvilinear coordinates in an Euclidean space

A noniterative method for the numerical generation of orthogonal curvilinear coordinates for plane annular regions between two arbitrary smooth closed curves was developed. The basic generating equation is the Gaussian equation for an Euclidean space which is solved analytically. The method is applied in many cases and these test results demonstrate that the proposed method can be readily applied to a wide variety of problems. The method can also be used for simply connected regions only by obtaining the solution of the linear equation under the changed boundary conditions.

Warsi, Z. U. A.↗

NASA NTRS1987

Solution of the surface Euler equations for accurate three-dimensional boundary-layer analysis of aerodynamic configurations

The three-dimensional boundary-layer equations in the limit as the normal coordinate tends to infinity are called the surface Euler equations. The present paper describes an accurate method for generating edge conditions for three-dimensional boundary-layer codes using these equations. The inviscid pressure distribution is first interpolated to the boundary-layer grid. The surface Euler equations are then solved with this pressure field and a prescribed set of initial and boundary conditions to yield the velocities along the two surface coordinate directions. Results for typical wing and fuselage geometries are presented. The smoothness and accuracy of the edge conditions obtained are found to be superior to the conventional interpolation procedures.

Iyer, V.↗

NASA NTRS1986

An efficient method for solving the steady Euler equations

An efficient numerical procedure for solving a set of nonlinear partial differential equations is given, specifically for the steady Euler equations. Solutions of the equations were obtained by Newton's linearization procedure, commonly used to solve the roots of nonlinear algebraic equations. In application of the same procedure for solving a set of differential equations we give a theorem showing that a quadratic convergence rate can be achieved. While the domain of quadratic convergence depends on the problems studied and is unknown a priori, we show that firstand second-order derivatives of flux vectors determine whether the condition for quadratic convergence is satisfied. The first derivatives enter as an implicit operator for yielding new iterates and the second derivatives indicates smoothness of the flows considered. Consequently flows involving shocks are expected to require larger number of iterations. First-order upwind discretization in conjunction with the Steger-Warming flux-vector splitting is employed on the implicit operator and a diagonal dominant matrix results. However the explicit operator is represented by first- and seond-order upwind differencings, using both Steger-Warming's and van Leer's splittings. We discuss treatment of boundary conditions and solution procedures for solving the resulting block matrix system. With a set of test problems for one- and two-dimensional flows, we show detailed study as to the efficiency, accuracy, and convergence of the present method.

Liou, M. S.↗

DOE OSTI2026

Electropolishing-induced topographic defects in niobium: Insights and Implications for SRF

Electropolishing (EP) is the premier surface preparation method for high-Q, high-gradient superconducting RF cavities made of Nb. This leaves behind an apparently smooth surface, yet the achievable peak magnetic fields fall well below the superheating field of Nb, in most cases. In this work, the ultimate surface finish of EP was investigated by studying its effect on highly polished Nb samples. EP introduces high slope angle sloped-steps at grain boundaries. The magnetic field enhancement and superheating field suppression factors associated with such a geometry are calculated in the London theory. Despite the by-eye smoothness of electropolished Nb, such defects compromise the stability of the low-loss Meissner state, likely limiting the achievable peak accelerating fields in superconducting RF cavities. Finally, the impact of surface roughness on impurity diffusion is investigated which can link surface roughness to the effectiveness of heat treatments like low-temperature baking or nitrogen infusion in the vortex nucleation or hydride hypotheses. Surface roughness tends to decrease the effective dose of impurities as a result of the expansion of impurities into regions with greater internal angle. The effective dose of impurities can be protected by minimizing slope angles and step heights, ensuring uniformity.

Lechner, Eric [Thomas Jefferson National Accelerat↗