Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smoothing methods”

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.

At least 469 records · Page 26

Advances in Automation of Overset Structured Volume Mesh Generation and Domain Connectivity

Automation of overset structured surface mesh generation has recently been accomplished by the creation of face, edge, and node meshes based on Boundary Representation solids as the geometry input. The current work continues the automation effort in overset volume mesh generation and domain connectivity based on the auto-generated surface meshes. All near-body curvilinear volume meshes are automatically generated using hyperbolic methods. Automation of this step is enabled by appropriate surface grid point distribution, and selection of boundary-splay and smoothing parameters based on concave and convex surface features. The off-body domain is covered by two automatically generated grid systems. The first contains a single Cartesian mesh with a uniform core enclosing all near-body volume meshes and stretched layers that extend to the far field, while the second consists of a set of small stretched Cartesian grids covering pockets of off-body orphan points. With high quality mesh overlap mostly guaranteed by the surface meshing scheme, orphan points that need to be covered by the second Cartesian mesh system are located away from the fine grid spacing region near the wall. Using line-segment and ray-pierce tests against the surface grids, hole-cutting is accomplished on both near and off-body volume grids resulting in appropriate clearances from the wall. The complete mesh generation automation process is demonstrated on five test cases where flow solutions are also computed and compared with solutions obtained using other methods.

TTT↗

Method of making a rocket nozzle

A method is described for forming the interior of a nozzle having uneven walls so that a throat of smooth converging and diverging sides is provided for passing flow. A metallic insert material is placed within the flow passageway adjacent to the area where the sharper throat constriction is to be formed, so that the material will flow through the inlet into the throat space when liquefied.

Campbell, D. H.↗

Computation of aerodynamic interference effects on oscillating airfoils with controls in ventilated subsonic wind tunnels

Lift interference effects are discussed based on Bland's (1968) integral equation. A mathematical existence theory is utilized for which convergence of the numerical method has been proved for general (square-integrable) downwashes. Airloads are computed using orthogonal airfoil polynomial pairs in conjunction with a collocation method which is numerically equivalent to Galerkin's method and complex least squares. Convergence exhibits exponentially decreasing error with the number n of collocation points for smooth downwashes, whereas errors are proportional to 1/n for discontinuous downwashes. The latter can be reduced to 1/n to the m+1 power with mth-order Richardson extrapolation (by using m = 2, hundredfold error reductions were obtained with only a 13% increase of computer time). Numerical results are presented showing acoustic resonance, as well as the effect of Mach number, ventilation, height-to-chord ratio, and mode shape on wind-tunnel interference. Excellent agreement with experiment is obtained in steady flow, and good agreement is obtained for unsteady flow.

Fromme, J. A.↗

A computational fluid dynamic technique valid at the centerline for non-axisymmetric problems in cylindrical coordinates

A technique is described for the numerical solution of non-axisymmetric flow problems posed in cylindrical coordinates when the z-axis is included in the flowfield. The highlight of the technique is the manner in which the singularities at the centerline are handled. Specifically, the governing flowfield equations at r = 0 are put in a special form by applying L'Hospital's Rule. The required radial derivatives are evaluated using a one-sided, second-order accurate, first-difference. This leads to a smooth, convergent calculation of the flowfield at the centerline. This appears to be the first generally applicable numerical method for avoiding coordinate system singularities in the context of a finite-difference scheme, and could have application to many nonaxisymmetric flows. The technique is illustrated by specific results for the time-dependent flowfield inside an internal combustion engine.

Griffin, M. D.↗

Method of and apparatus for generating an interstitial point in a data stream having an even number of data points

Apparatus for doubling the data density rate of an analog to digital converter or doubling the data density storage capacity of a memory deviced is discussed. An interstitial data point midway between adjacent data points in a data stream having an even number of equal interval data points is generated by applying a set of predetermined one-dimensional convolute integer coefficients which can include a set of multiplier coefficients and a normalizer coefficient. Interpolator means apply the coefficients to the data points by weighting equally on each side of the center of the even number of equal interval data points to obtain an interstital point value at the center of the data points. A one-dimensional output data set, which is twice as dense as a one-dimensional equal interval input data set, can be generated where the output data set includes interstitial points interdigitated between adjacent data points in the input data set. The method for generating the set of interstital points is a weighted, nearest-neighbor, non-recursive, moving, smoothing averaging technique, equivalent to applying a polynomial regression calculation to the data set.

Edwards, T. R.↗

Global collocation methods for approximation and the solution of partial differential equations

Polynomial interpolation methods are applied both to the approximation of functions and to the numerical solutions of hyperbolic and elliptic partial differential equations. The derivative matrix for a general sequence of the collocation points is constructed. The approximate derivative is then found by a matrix times vector multiply. The effects of several factors on the performance of these methods including the effect of different collocation points are then explored. The resolution of the schemes for both smooth functions and functions with steep gradients or discontinuities in some derivative are also studied. The accuracy when the gradients occur both near the center of the region and in the vicinity of the boundary is investigated. The importance of the aliasing limit on the resolution of the approximation is investigated in detail. Also examined is the effect of boundary treatment on the stability and accuracy of the scheme.

Solomonoff, A.↗

Separated flow

A brief overview of flow separation phenomena is provided. Langley has many active research programs in flow separation related areas. Three cases are presented which describe specific examples of flow separation research. In each example, a description of the fundamental fluid physics and the complexity of the flow field is presented along with a method of either reducing or controlling the extent of separation. The following examples are discussed: flow over a smooth surface with an adverse pressure gradient; flow over a surface with a geometric discontinuity; and flow with shock-boundary layer interactions. These results will show that improvements are being made in the understanding of flow separation and its control.

Sellers, W. L., III↗

Direct numerical simulation of instabilities in parallel flow with spherical roughness elements

Results from a direct numerical simulation of laminar flow over a flat surface with spherical roughness elements using a spectral-element method are given. The numerical simulation approximates roughness as a cellular pattern of identical spheres protruding from a smooth wall. Periodic boundary conditions on the domain's horizontal faces simulate an infinite array of roughness elements extending in the streamwise and spanwise directions, which implies the parallel-flow assumption, and results in a closed domain. A body force, designed to yield the horizontal Blasius velocity in the absence of roughness, sustains the flow. Instabilities above a critical Reynolds number reveal negligible oscillations in the recirculation regions behind each sphere and in the free stream, high-amplitude oscillations in the layer directly above the spheres, and a mean profile with an inflection point near the sphere's crest. The inflection point yields an unstable layer above the roughness (where U''(y) is less than 0) and a stable region within the roughness (where U''(y) is greater than 0). Evidently, the instability begins when the low-momentum or wake region behind an element, being the region most affected by disturbances (purely numerical in this case), goes unstable and moves. In compressible flow with periodic boundaries, this motion sends disturbances to all regions of the domain. In the unstable layer just above the inflection point, the disturbances grow while being carried downstream with a propagation speed equal to the local mean velocity; they do not grow amid the low energy region near the roughness patch. The most amplified disturbance eventually arrives at the next roughness element downstream, perturbing its wake and inducing a global response at a frequency governed by the streamwise spacing between spheres and the mean velocity of the most amplified layer.

Deanna, R. G.↗

Mapping ERS-1 wind fields over north west Atlantic using a variational objective analysis

A variational method is implemented to produce five day mean gridded ERS-1 analyzed wind fields in the north west Atlantic with the aim of providing for wind forcing of basin scale ocean models. The method consists of minimizing a cost functional, designed to measure misfits to prescribed weighted constraints which express a smoothed behavior and the proximity to input data vectors and curl. The weights are empirically determined by comparison with independent ship and buoy data over four five day periods. Root mean square differences between analyzed winds and independent data are thus further decreased: they range from 0.8 up to 1.8 m/s. Some problems present in the initial data remain; they are principally due to incomplete data coverage (instrumental problems), and possibly unresolved ambiguities. The resulting curl fields are smoothed and show coherent patterns. A comparison with the European Center for Medium range Weather Forecasting (ECMWF) analysis is encouraging.

Siefridt, L.↗

Two-dimensional radiative transfer in cloudy atmospheres - The spherical harmonic spatial grid method

A new two-dimensional monochromatic method that computes the transfer of solar or thermal radiation through atmospheres with arbitrary optical properties is described. The model discretizes the radiative transfer equation by expanding the angular part of the radiance field in a spherical harmonic series and representing the spatial part with a discrete grid. The resulting sparse coupled system of equations is solved iteratively with the conjugate gradient method. A Monte Carlo model is used for extensive verification of outgoing flux and radiance values from both smooth and highly variable (multifractal) media. The spherical harmonic expansion naturally allows for different levels of approximation, but tests show that the 2D equivalent of the two-stream approximation is poor at approximating variations in the outgoing flux. The model developed here is shown to be highly efficient so that media with tens of thousands of grid points can be computed in minutes. The large improvement in efficiency will permit quick, accurate radiative transfer calculations of realistic cloud fields and improve our understanding of the effect of inhomogeneity on radiative transfer in cloudy atmospheres.

Evans, K. F.↗

Atmospheric Infrared Sounder

The microwave 'first-guess' algorithm was run on the cloudy test simulations. Eight datasets were considered in the cloudy test, comprising approximately 360 retrievals, of which one failed to converge. Retrievals were done on the AMSU-A grid. Examination of the true profiles (provided for the A, C and D tracks) showed numerous cases of very pronounced temperature inversion layers in the troposphere which the retrieval does not have enough vertical resolution to reproduce. A typical example with an inversion layer near 700 mbar is shown in Figure 1. The inversion layers also exhibit strong vertical gradients of water vapor which are not resolved in the retrieval. (Water vapor volume density is given in g/sq cm per layer. Layer thickness is 20 mb from 200 to 400 mb, and 25 mb from 400 to 1,000 mb.) The retrievals do reproduce the overall smoothed shape of the profiles, and therefore as a first guess should be within the range of linear methods for IR retrievals using AIRS.

Rosenkranz, Philip W.↗

Heat Transfer and Flow on the First Stage Blade Tip of a Power Generation Gas Turbine: Experimental Results - Part 1

A combined computational and experimental study has been performed to investigate the detailed distribution of convective heat transfer coefficients on the first stage blade tip surface for a geometry typical of large power generation turbines(>100MW). This paper is concerned with the design and execution of the experimental portion of the study. A stationary blade cascade experiment has been run consisting of three airfoils, the center airfoil having a variable tip gap clearance. The airfoil models the aerodynamic tip section of a high pressure turbine blade with inlet Mach number of 0.30, exit Mach number of 0.75, pressure ratio of 1.45, exit Reynolds number based on axial chord of 2.57 x 10(exp 6), and total turning of about 110 degrees. A hue detection based liquid crystal method is used to obtain the detailed heat transfer coefficient distribution on the blade tip surface for flat, smooth tip surfaces with both sharp and rounded edges. The cascade inlet turbulence intensity level took on values of either 5% or 9%. The cascade also models the casing recess in the shroud surface ahead of the blade. Experimental results are shown for the pressure distribution measurements on the airfoil near the tip gap, on the blade tip surface, and on the opposite shroud surface. Tip surface heat transfer coefficient distributions are shown for sharp-edge and rounded-edge tip geometries at each of the inlet turbulence intensity levels.

Bunker, Ronald S.↗

Quantification of rat retinal growth and vascular population changes after single and split doses of proton irradiation: translational study using stereology methods

This study quantified architectural and population changes in the rat retinal vasculature after proton irradiation using stereology. A 100 MeV conformal proton beam delivered 8, 14, 20 and 28 Gy as single and split doses to the whole eye. The vascular networks were prepared from retinal digests. Stereological methods were used to obtain the area of the retina and unbiased estimates of microvessel/artery/vein endothelial, pericyte and smooth muscle population, and vessel length. The retinal area increased progressively in the unirradiated, age-matched controls and in the retinas irradiated with 8 and 14 Gy, indicating uniform progressive retinal growth. No growth occurred after 20 and 28 Gy. Regression analysis of total endothelial cell number in all vessels (arteries, veins and capillaries) after irradiation documented a progressive time- and dose-dependent cell loss occurring over 15 to 24 months. The difference from controls was significant (P<0.01) after 28 Gy given in single and split doses and after 20 Gy given as a split dose (P<0.05). Total vessel length in microvessel was significantly shortened at 20 and 28 Gy compared to that of controls (P<0.05). No evident dose recovery was observed in the endothelial populations after split doses. At 10 Gy, the rate of endothelial cell loss, a dose parameter used to characterize the time- and dose-dependent loss of the endothelial population, was doubled.

NASA Discipline Radiation Health↗

Effects of Strain Rate and Surface Finish on Tensile Properties in High Pressure Hydrogen — A Case Study for LP-DED NASA HR-1

NASA has advanced NASA HR-1 as a solution for liquid rocket engine components that operate in hydrogen rich environments using various AM techniques. The NASA HR-1 for AM applications was specifically formulated for high ductility in high pressure gaseous hydrogen environment in addition to high strength and low cycle fatigue (LCF) resistance. NASA HR-1 meets materials requirements for liquid rocket engine components, including good hydrogen resistance, high conductivity, good LCF performance, and high elongation and strength for components in medium heat flux environments. NASA has completed fabrication of several subscale and full-scale channel wall nozzles in LP-DED NASA HR-1 and completed hot-fire testing. This includes refinement of the process to produce thin-walls and various channel geometries to meet the requirements for channel wall nozzle applications. The main objective of this investigation is to better understand the effects of tensile test strain rate and specimen surface finish on the amount of tensile ductility degradation in high pressure gaseous hydrogen environments. Performing smooth tensile tests in high pressure hydrogen and nitrogen environments can serve as a rapid and economical screening method to determine the relative susceptibility of materials to HEE. This paper presents the tensile test results for LP-DED NASA HR-1 in hydrogen and details the material evaluation on the effects of strain rate and surface finish condition on surface cracking behavior and fracture characteristics. The resistance to HEE for LP-DED NASA HR-1 and the wrought alloy is compared under the same testing conditions. Loss of tensile ductility is quantified in terms of changes in elongation at failure in hydrogen to compare HEE susceptibility. The most suitable strain rate and specimen surface finish condition for future HEE screening testing at MSFC are recommended.

Additive manufacturing↗

Unsteady transonic Euler solutions using finite elements

A finite element solution of the unsteady Euler equations is presented and demonstrated for 2D airfoil configurations oscillating in transonic flows. Computations are performed by spatially discretizing the conservation equations using the Galerkin weighted residual method and then employing a multistage Runge-Kutta scheme to march forward in time. A mesh deformation scheme has been developed to efficiently move interior points in a smooth fashion as the airfoil undergoes rigid body pitch and plunge motion. Both steady and unsteady results are presented, and a comparison is made with solutions obtained using finite-volume techniques. The effects of using either a lumped or consistent mass matrix are presented; the finite element method provides an accurate solution for unsteady transonic flows about isolated airfoils.

Davis, Gary A.↗

A Synthesis Method to Auralize Rotor Noise During Transitional Flight

This paper describes a new capability to synthesize Urban Air Mobility (UAM) vehicle rotor/propeller loading and thickness noise and broadband self noise for transitional flight auralization. A vehicle smoothly transitioning between flight states is approximated here as transitioning in time through a series of trimmed steady-state flight conditions. To generate information that helps synthesize loading and thickness noise for transitional flight, this paper describes the development of a tool that combines blade loading, motion, and geometry data from a series of separate trimmed flight states. This paper also describes a method to combine broadband self noise acoustic predictions, each representing a separate trimmed flight state, to synthesize broadband self noise for transitional flight. Different flight profiles of the six-passenger NASA UAM Quadrotor reference vehicle demonstrate these capabilities, with each profile specified by a series of trimmed flight states.

Auralization↗

Growth of Homoepitaxial ZnO Semiconducting Films

As a high temperature semiconductor, ZnO has been used for many applications such as wave-guide, solar cells, and surface acoustic wave devices. Since the ZnO material has an energy gap of 3.3 eV at room temperature and an excitonic binding energy (60 meV) that is possible to make excitonic lasering at room temperature a recent surge of interest is to synthesize ZnO films for electro-optical devices. These applications require films with a smooth surface, good crystal quality, and low defect density. Homoepitaxial films have been studied in terms of morphology, crystal structure, and electrical and optical properties. ZnO single crystals are grown by the hydrothermal method. Substrates are mechanically polished and annealed in air for four hours before deposited films. The annealing temperature-dependence of ZnO substrates is studied. Films are synthesized by the off-axis reactive sputtering deposition. The films have very smooth surface with a roughness </= 5 nm on a 5 micrometers x 5 micrometers area. The full width at half maximum of film theta rocking curves measured by the x-ray diffraction is slightly larger than the crystal substrate. The film quality is determined by measuring the film resistivity, the Hall mobility, carrier densities and the energy band gap. The properties of ZnO films grown of (0001) ZnO and (0001) sapphire substrates will be also compared and discussed in the presentation.

Zhu, Shen↗