Search NASA⌕ Search

SEARCH · Search NASA

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 505 records · Page 28

Accelerated basis adaptation in homogeneous chaos spaces

Polynomial chaos expansions (PCE) provide an efficient approach to uncertainty quantification (UQ) and have been adapted to diverse applications across the spectrum of science and engineering. For situations involving large stochastic parameterizations, the curse of dimensionality renders PCE-based methods computationally prohibitive. A basis adaptation approach for PCE was proposed by Tipireddy and Ghanem (2014) which transforms the input random variables through an isometry such that sufficient probabilistic characterization of specific quantities of interest (QoI) is concentrated in an algebraic manifold embedded in the linear span of the dominant transformed variables. While quite versatile, that original version of the PCE basis adaptation exhibited slow convergence for a number of problems of practical significance. In the present paper, we propose two novel methods to accelerate the convergence of the original basis adaptation approach, thus expanding its reach while also providing insight into its performance. In the first method, information gained from a pilot PCE representation is used to correct the mean and Gaussian coefficients in the adapted space. By taking advantages of probabilistic information in higher dimensional adaptation gleaned from an initial adaptation, the second method updates the rotation matrix used to identify the dominant transformed variables. In this manner, the new rotation matrix concentrates even more probabilistic information in its first few dimensions. These two method can be combined to achieve even better performance, the combined method is referred to as sequentially optimized adaptation method. Here, the methods are demonstrated on an analytical test function and a model of a space structure with several sub-components and a non-smooth quantity of interest representing the maximum acceleration over time. Both methods achieve accelerated convergence of the basis adaptation approach with negligible additional costs.

Accelerated convergence↗

Hybrid PN-SN Calculations with SAAF for the Multiscale Transport Capability in Rattlesnake

Two interface conditions, the Lagrange multiplier method and the upwinding method, for hybrid \pn-\sn calculations is proposed for the self-adjoint angular flux (SAAF) formulation of the transport equation using the continuous finite element method (FEM) for spatial discretization. These interface conditions are implemented in Rattlesnake, the radiation transport application built on MOOSE, for the on-going multiscale transport simulation effort at INL. For smoothing the solution at the interface for the Lagrange multiplier method, a method based on \sn Lagrange interpolation on the sphere is proposed. Numerical results indicate that the interface conditions give the expected convergence.

97 MATHEMATICS AND COMPUTING↗

Airplane Mesh Development with Grid Density Studies

Automatic Grid Generation Wish List Geometry handling, including CAD clean up and mesh generation, remains a major bottleneck in the application of CFD methods. There is a pressing need for greater automation in several aspects of the geometry preparation in order to reduce set up time and eliminate user intervention as much as possible. Starting from the CAD representation of a configuration, there may be holes or overlapping surfaces which require an intensive effort to establish cleanly abutting surface patches, and collections of many patches may need to be combined for more efficient use of the geometrical representation. Obtaining an accurate and suitable body conforming grid with an adequate distribution of points throughout the flow-field, for the flow conditions of interest, is often the most time consuming task for complex CFD applications. There is a need for a clean unambiguous definition of the CAD geometry. Ideally this would be carried out automatically by smart CAD clean up software. One could also define a standard piece-wise smooth surface representation suitable for use by computational methods and then create software to translate between the various CAD descriptions and the standard representation. Surface meshing remains a time consuming, user intensive procedure. There is a need for automated surface meshing, requiring only minimal user intervention to define the overall density of mesh points. The surface mesher should produce well shaped elements (triangles or quadrilaterals) whose size is determined initially according to the surface curvature with a minimum size for flat pieces, and later refined by the user in other regions if necessary. Present techniques for volume meshing all require some degree of user intervention. There is a need for fully automated and reliable volume mesh generation. In addition, it should be possible to create both surface and volume meshes that meet guaranteed measures of mesh quality (e.g. minimum and maximum angle, stretching ratios, etc.).

Cliff, Susan E.↗

Automatic cloud tracking applied to GOES and Meteosat observations

An improved automatic processing method for the tracking of cloud motions as revealed by satellite imagery is presented and applications of the method to GOES observations of Hurricane Eloise and Meteosat water vapor and infrared data are presented. The method is shown to involve steps of picture smoothing, target selection and the calculation of cloud motion vectors by the matching of a group at a given time with its best likeness at a later time, or by a cross-correlation computation. Cloud motion computations can be made in as many as four separate layers simultaneously. For data of 4 and 8 km resolution in the eye of Hurricane Eloise, the automatic system is found to provide results comparable in accuracy and coverage to those obtained by NASA analysts using the Atmospheric and Oceanographic Information Processing System, with results obtained by the pattern recognition and cross correlation computations differing by only fractions of a pixel. For Meteosat water vapor data from the tropics and midlatitudes, the automatic motion computations are found to be reliable only in areas where the water vapor fields contained small-scale structure, although excellent results are obtained using Meteosat IR data in the same regions. The automatic method thus appears to be competitive in accuracy and coverage with motion determination by human analysts.

Endlich, R. M.↗

A frozen Gaussian approximation-based multi-level particle swarm optimization for seismic inversion

In this paper, we propose a frozen Gaussian approximation (FGA)-based multi-level particle swarm optimization (MLPSO) method for seismic inversion of high-frequency wave data. The method addresses two challenges in it: First, the optimization problem is highly non-convex, which makes hard for gradient-based methods to reach global minima. This is tackled by MLPSO which can escape from undesired local minima. Second, the character of high-frequency of seismic waves requires a large number of grid points in direct computational methods, and thus renders an extremely high computational demand on the simulation of each sample in MLPSO. We overcome this difficulty by three steps: First, we use FGA to compute high-frequency wave propagation based on asymptotic analysis on phase plane; Then we design a constrained full waveform inversion problem to prevent the optimization search getting into regions of velocity where FGA is not accurate; Last, we solve the constrained optimization problem by MLPSO that employs FGA solvers with different fidelity. The performance of the proposed method is demonstrated by a two-dimensional full-waveform inversion example of the smoothed Marmousi model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Transonic flow calculations using a flux vector splitting method for the Euler equations

A study of the flux vector splitting method of Steger and Warming for the solution of the time dependent Euler equations in strong conservation law form for arbitrary two-dimensional geometries is presented. The procedure employed here differs from that of Buning and Steger in that it uses a different algorithm and employs implicit boundary conditions. Moreover, the method, as implemented here, does not contain any explicit smoothing or any adjustable parameters. Calculations were carried out for an NACA 0012 airfoil at various Mach numbers and angles of attack, and cylinders. Steady symmetric solutions were obtained for the full cylinder at a freestream Mach number of .5 without imposing a symmetry condition. In general, good agreement with other methods was obtained.

Seaford, C. M.↗

Tank Tests of the Effect of Rivet Heads, etc., on the Water Performance of a Seaplane Float, Special Report

A 1/3.5 full-size model of the Mark V float of the Bureau of Aeronautics, Navy Department, was tested in the NACA tank both with smooth painted bottom surfaces and with roundhead rivets, plate laps, and keel plates fitted to simulate the actual bottom of a metal float. The augmentation in water resistance due to the added roughness was found to be from 10-12% at the hum speed and from 12-14% at high speeds. The effect of the roughness of the afterbody was found to be negligible except at high trims. The model data were extrapolated to full size by the usual method which assumes the forces to vary according to Froude's law, and in the case of the smooth model by a method of separation that takes into account the effect of scale on the frictional resistance. It was concluded that the effect of rivet heads on the takeoff performance of a relatively high-powered float seaplane is of little consequence but that it may be of greater importance in the case of more moderately powered flying boats.

Parkinson, J. B.↗

Spectral Methods for Computational Fluid Dynamics

As a tool for large-scale computations in fluid dynamics, spectral methods were prophesized in 1944, born in 1954, virtually buried in the mid-1960's, resurrected in 1969, evangalized in the 1970's, and catholicized in the 1980's. The use of spectral methods for meteorological problems was proposed by Blinova in 1944 and the first numerical computations were conducted by Silberman (1954). By the early 1960's computers had achieved sufficient power to permit calculations with hundreds of degrees of freedom. For problems of this size the traditional way of computing the nonlinear terms in spectral methods was expensive compared with finite-difference methods. Consequently, spectral methods fell out of favor. The expense of computing nonlinear terms remained a severe drawback until Orszag (1969) and Eliasen, Machenauer, and Rasmussen (1970) developed the transform methods that still form the backbone of many large-scale spectral computations. The original proselytes of spectral methods were meteorologists involved in global weather modeling and fluid dynamicists investigating isotropic turbulence. The converts who were inspired by the successes of these pioneers remained, for the most part, confined to these and closely related fields throughout the 1970's. During that decade spectral methods appeared to be well-suited only for problems governed by ordinary diSerential eqllations or by partial differential equations with periodic boundary conditions. And, of course, the solution itself needed to be smooth. Some of the obstacles to wider application of spectral methods were: (1) poor resolution of discontinuous solutions; (2) inefficient implementation of implicit methods; and (3) drastic geometric constraints. All of these barriers have undergone some erosion during the 1980's, particularly the latter two. As a result, the applicability and appeal of spectral methods for computational fluid dynamics has broadened considerably. The motivation for the use of spectral methods in numerical calculations stems from the attractive approximation properties of orthogonal polynomial expansions.

Zang, T. A.↗

Extrusion of Complex Surface Meshes Utilizing Face Offsetting and Mean Curvature Smoothing

Grids for three-dimensional Computational Fluid Dynamics problems frequently require a prismatic layer of cells, typically extruded in the off-body direction from a two-dimensional surface mesh to properly resolve boundary layers. When the surface geometry is complex, the extrusion process can fail, resulting in the physical boundary layer being thicker than the prismatic layer, leading to under-resolved boundary layers. To address the shortcomings of existing grid generation tools, a new tool has been developed as part of the Mesh Tools suite, that is capable of extruding complex surface meshes. The new tool uses a face offsetting method to preserve surface curvature and a mean curvature smoothing algorithm to prevent the cells from self-intersecting in concave regions. Testing of the new tool found that not only were complex surface meshes able to be extruded to the desired thickness, but extrusion of simple surface meshes was also improved, due to the face offsetting method.

mesh extrusion↗

Objective scatterometer wind ambiguity removal using smoothness and dynamical constraints

In the present investigation, a variational analysis method (VAM) is used to remove the ambiguity of the Seasat-A Satellite Scatterometer (SASS) winds. At each SASS data point, two, three, or four wind vectors (termed ambiguities) are retrieved. It is pointed out that the VAM is basically a least squares method for fitting data. The problem may be nonlinear. The best fit to the data and constraints is obtained on the basis of a minimization of the objective function. The VAM was tested and tuned at 12 h GMT Sept. 10, 1978. Attention is given to a case study involving an intense cyclone centered south of Japan at 138 deg E.

Hoffman, R. N.↗

Near atomically smooth alkali antimonide photocathode thin films

Nano-roughness is one of the major factors degrading the emittance of electron beams that can be generated by high efficiency photocathodes, such as the thermally reacted alkali antimonide thin films. In this paper, we demonstrate a co-deposition based method for producing alkali antimonide cathodes that produce near atomic smoothness with high reproducibility. Here, we calculate the effect of the surface roughness on the emittance and show that such smooth cathode surfaces are essential for operation of alkali antimonide cathodes in high field, low emittance radio frequency electron guns and to obtain ultracold electrons for ultrafast electron diffraction applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Tuning the Interpolation Basis in a Multigrid Decomposition for Local Error Control

In the compression of scientific data, error-controlled compressors enable to considerably decrease the size of the dataset while maintaining adequate levels of accuracy. In this paper, we note that multi-level refactoring scheme such as MGARD i) rely on an approximation of the data based on the interpolation of coefficients, ii) estimate the resulting error with global metrics on the dataset. To improve on these two aspects, we propose a method that aims to divide the original dataset into blocks based on their smoothness and refactors each block separately with the most relevant interpolation order. We show the relevance of such a method on tailored datasets and the benefits and challenges when applying it to large scientific data.

Vidal, Nicolas [ORNL]↗

Tank Tests to Show the Effect Rivet Heads on the Water Performance of a Seaplane-Float

A 1/3.5 full-sized model of a seaplane float constructed from lines supplied by the Bureau of Aeronautics, Navy Department, was tested first with smooth painted bottom surfaces and then with round-head rivets, plate laps, and keel plates fitted to simulate the actual bottom of a metal float. A percentage increase in water resistance caused by the added roughness was found to be from 5 to 20 percent at the hump speed and from 15 to 40 percent at high speeds. The effect of the roughness of the afterbody was found to be negligible except at high trims. The model data were extrapolated to full size by the usual method that assumes the forces to vary according to Froude's law and, in the case of the smooth model, by a method of separation that takes into account the effect of scale on the frictional resistance. It was concluded that the effect of rivet heads on the take-off performance of a relatively high-powered float seaplane is of little consequence, but it may be of greater importance in the case of more moderately powered flying boats.

FLOATS - RESISTANCE EFFECT OF SURFACE ROUGHNESS↗

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sr{sub 1.98}Eu{sub 0.02}SiO{sub 4} luminescence whisker based on vapor-phase deposition: Facile synthesis, uniform morphology and enhanced luminescence properties

Highlights: • For the first time, it is possible to obtain Sr{sub 1.98}Eu{sub 0.02}SiO{sub 4} whisker. • The whiskers are smooth and uniform with L/D ratio over 50. • Durability and thermal stability of the whisker are enhanced. - Abstract: A high performance strontium silicate phosphor has been successfully synthesized though a facile vapor-phase deposition method. The product consists of single crystal whiskers which are smooth and uniform, and with a sectional equivalent diameter of around 5 μm; the aspect ratio is over 50 and no agglomeration can be observed. X-ray diffraction result confirmed that the crystal structure of the whisker was α’-Sr{sub 2}SiO{sub 4}. The exact chemical composition was Sr{sub 1.98}Eu{sub 0.02}SiO{sub 4} which was analyzed by energy dispersive spectrometer and inductively coupled plasma-mass spectrometer. The whisker shows broad green emission with peak at 523 nm ranging from 470 to 600 nm (excited at 370 nm). Compared with traditional Sr{sub 2}SiO{sub 4}:Eu phosphor, durability (at 85% humidity and 85 °C) and thermal stability of the whisker are obviously improved. Moreover, growth mechanism of the Sr{sub 1.98}Eu{sub 0.02}SiO{sub 4} whiskers is Vapor–Liquid–Solid. On a macro-scale, the product is still powder which makes it suitable for the current packaging process of WLEDs.

36 MATERIALS SCIENCE↗

A Facile Degumming Method of Kenaf Fibers Using Deep Eutectic Solution

An eco-friendly and effective degumming method needed to be developed for producing kenaf bast fibers. As such, this study devised a novel deep eutectic solvent (DES) method coupled with microwave and alkaline-ultrasonic treatment. The novel method could effectively remove the gummy matters, providing a smooth and clean fiber surface. The properties of the fibers were assessed including chemical compositions, surface structure, crystallinity index (66.68%) and thermal properties. The residual gum content (9.419%), fiber fineness (4.125tex), breaking tenacity (13.650 cn/tex) of the refined dry fibers produced by the novel treatment were comparable with the fibers produced by the traditional two-step alkali boiling method. Besides, the novel method could reduce the usage of chemical, water and time by 48.9%, by 66.7% and 66.8%, respectively. These results revealed that the novel combined DES pretreatment is a practical and feasible pretreatment method for kenaf bast degumming, demonstrating its facile, green energy-saving and fast properties in the degumming process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Three-dimensional self-adaptive grid method for complex flows

A self-adaptive grid procedure for efficient computation of three-dimensional complex flow fields is described. The method is based on variational principles to minimize the energy of a spring system analogy which redistributes the grid points. Grid control parameters are determined by specifying maximum and minimum grid spacing. Multidirectional adaptation is achieved by splitting the procedure into a sequence of successive applications of a unidirectional adaptation. One-sided, two-directional constraints for orthogonality and smoothness are used to enhance the efficiency of the method. Feasibility of the scheme is demonstrated by application to a multinozzle, afterbody, plume flow field. Application of the algorithm for initial grid generation is illustrated by constructing a three-dimensional grid about a bump-like geometry.

Djomehri, M. Jahed↗