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

Boundary-fitted coordinate systems for arbitrary computational regions

The method of Smith and Wiegel was used to generate meshes for mixer lobes and subsonic inlets that are compatible with flow analysis codes requiring a boundary fitted coordinate system. Successful application of this mesh generator required development of procedures to distribute the mesh points along the boundaries, to regulate the dependence of the connecting function to the local boundary slope, to concentrate the mesh into regions of special interest, and to modify the mesh grid so that it possessed a smooth progression of cell metrics and cell volumes in all directions. The method of Smith and Wiegel when coupled with the extensions mentioned above has proven to be easy to use and control for the inlet and mixer lobe geometries investigated.

Kowalski, E. J.↗

Vortex methods for separated flows

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented and includes the relationship with traditional point-vortex studies, convergence to smooth solutions of the Euler equations, and the essential differences between two and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. Two-dimensional flows around bluff bodies are emphasized. Robustness of the method and the assessment of accuracy, vortex-core profiles, time-marching schemes, numerical dissipation, and efficient programming are treated. Operation counts for unbounded and periodic flows are given, and two algorithms designed to speed up the calculations are described.

Spalart, Philippe R.↗

Application and sensitivity investigation of Fourier transforms for microwave radiometric inversions

Existing microwave radiometer technology now provides a suitable method for remote determination of the ocean surface's absolute brightness temperature. To extract the brightness temperature of the water from the antenna temperature equation, an unstable Fredholm integral equation of the first kind was solved. Fast Fourier Transform techniques were used to invert the integral after it is placed into a cross-correlation form. Application and verification of the methods to a two-dimensional modeling of a laboratory wave tank system were included. The instability of the Fredholm equation was then demonstrated and a restoration procedure was included which smooths the resulting oscillations. With the recent availability and advances of Fast Fourier Transform techniques, the method presented becomes very attractive in the evaluation of large quantities of data. Actual radiometric measurements of sea water are inverted using the restoration method, incorporating the advantages of the Fast Fourier Transform algorithm for computations.

Holmes, J. J.↗

Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method

In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.

58 GEOSCIENCES↗

Uniformly high order accurate essentially non-oscillatory schemes 3

In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear.

Harten, A.↗

Euler solutions for blunt bodies using triangular meshes - Artificial viscosity forms and numerical boundary conditions

A finite volume method is used to calculate compressible inviscid flows over blunt bodies using, in general, unstructured grids. Artificial viscosity forms are derived based on a simplified least squares procedure. The extra second order terms are consistent with the governing equations, hence a systematic treatment of the numerical boundary conditions can be easily implemented. A special treatment of blunt bodies may be required. The discrete equations are linearized and the resulting system is solved by a relaxation method. Preliminary results indicate that the effect of the numerical dissipation is minimal. For subsonic flows over smooth bodies, the solution is practically vorticity-free and the total pressure loss is of the same order as the truncation error. Finally, some extensions of the present method are briefly discussed.

Winterstein, R.↗

Optimization of La 2–x Sr x CuO 4 Single Crystal Film Growth via Molecular Beam Epitaxy

Atomic layer-by-layer molecular beam epitaxy (ALL-MBE) combined with ozone is one of the best methods to fabricate single-crystal thin films of complex oxides. Cuprate such as La 2-x Sr x CuO 4 (LSCO) is a representative complex-oxide high-temperature superconductor (HTS) material. Our group utilizes this method to produce high-quality single-crystal HTS films with atomically smooth surfaces and interfaces. In addition, ALL-MBE enables us to engineer multilayer heterostructures with atomic precision. This allows the fabrication of tunnel junctions, various nanostructures, and other HTS devices of interest for superconducting electronics. We have synthesized over three thousand LSCO thin films in the past two decades. These films’ structural and electronic properties have been studied and characterized by various methods. Here, we distill the extensive experience we accumulated into a step-by-step protocol to fabricate atomically perfect LSCO films. The recipe includes substrate preparation, ozone generation and distillation, source calibration, the in situ monitoring of the film synthesis, post-growth annealing, and ex situ characterization. It discloses a reproducible way to fabricate single-crystal LSCO films for basic research and HTS electronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Alternative Attitude Commanding and Control for Precise Spacecraft Landing

A report proposes an alternative method of control for precision landing on a remote planet. In the traditional method, the attitude of a spacecraft is required to track a commanded translational acceleration vector, which is generated at each time step by solving a two-point boundary value problem. No requirement of continuity is imposed on the acceleration. The translational acceleration does not necessarily vary smoothly. Tracking of a non-smooth acceleration causes the vehicle attitude to exhibit undesirable transients and poor pointing stability behavior. In the alternative method, the two-point boundary value problem is not solved at each time step. A smooth reference position profile is computed. The profile is recomputed only when the control errors get sufficiently large. The nominal attitude is still required to track the smooth reference acceleration command. A steering logic is proposed that controls the position and velocity errors about the reference profile by perturbing the attitude slightly about the nominal attitude. The overall pointing behavior is therefore smooth, greatly reducing the degree of pointing instability.

Singh, Gurkirpal↗

Method for error correction in scanning probe microscopy

Disclosed here is a scanning probe microscope system and method for operating the same for producing scanning probe microscope images at fast scan rates and reducing oscillation artifacts. In some embodiments, an inverse consistent image registration method is used to align forward and backward scan traces for each line of the scanning microscope image. In some embodiments, the aligned forward and backward scan traces are combined using a weighting factor favoring the scan trace with higher smoothness. In some embodiments, the scanning probe microscope image is a potentiometry map and a method is provided to extract from the potentiometry map a conductivity map.

47 OTHER INSTRUMENTATION↗

Efficient High Order Central Schemes for Multi-Dimensional Hamilton-Jacobi Equations: Talk Slides

This viewgraph presentation presents information on the attempt to produce high-order, efficient, central methods that scale well to high dimension. The central philosophy is that the equations should evolve to the point where the data is smooth. This is accomplished by a cyclic pattern of reconstruction, evolution, and re-projection. One dimensional and two dimensional representational methods are detailed, as well.

Bryson, Steve↗

High precision control of laser energy for laser–matter interaction studies

Precise, highly reproducible control of the laser energy is required for high confidence laser–matter interaction research such as in dynamic compression science and high energy density physics. The energy must be adjustable without affecting the pulse shape (time varying intensity) or beam smoothness. We have developed a convenient two-stage energy tuning method for a nominal 100 J, 351 nm (UV) laser. The energy is adjusted in 10 J (10%) increments by operating the laser at full energy and inserting a beam splitter in the laser output. As the splitter is located after the final frequency tripling optics, the UV pulse shape is unchanged. The energy is varied by substituting a splitter of different reflectivity. For finer 3 J (3%) increments, the infrared pulse is attenuated inside the laser before the final amplifier. This requires modest tuning to preserve the pulse shape. The demonstrated variation in shot-to-shot reproducibility is less than +/–2.5 J (5% of the full energy), irrespective of the laser output energy. These approaches can be adapted to most ~100 J class lasers. We describe these techniques and show two examples where they have elucidated the underlying physics in laser shock compression experiments. One used only the beam splitters to establish the pressure for melting in iron. The other combined both techniques to finely increment the peak stress (~2 GPa steps) in germanium to precisely determine the onset and completion of melting—including the melting kinetics. Furthermore, these unambiguous results would not be possible without the developments described here.

47 OTHER INSTRUMENTATION↗

A General Spatiotemporal Imputation Framework for Missing Sensor Data

Many applications from precision agriculture, environmental monitoring and transportation networks rely on data collected across space and time over a large geographic area. Missing data poses a significant challenge for any data-driven inference and control tasks. Data imputation or the estimation of missing data can help fill these gaps by utilizing inherent spatial relationships and temporal patterns. A variety of spatiotemporal imputation models have been developed to address missing data in spatiotemporal datasets. However, these classical methods rely on the assumption that the underlying data follows a smooth trend and fail to provide accurate estimates when there is a large number of missing points in the data. Even though there are machine learning driven tensor completion approaches such as convolutional neural network based tensor completion (CoSTCo) that capture the non-linear relationships in the dataset, the transductive nature makes the algorithm less scalable. Thus, existing approaches for estimating the missing information do not effectively capture all dimensions of the spatiotemporal data structure, resulting in erroneous predictions and poor performance. The main contributions of this paper are: (1) We propose a novel inductive framework (G-LSTM) for missing data imputation that integrates a graph neural network with LSTMs to effectively capture both spatial and temporal dependencies. (2) Experimental results on a traffic dataset demonstrate that the proposed GNN integrated with an LSTM framework achieves improved imputation and maintains steady performance even when there are extreme missing conditions in comparison with the state-of-the-art imputation framework (i.e, CoSTCo). (3) The simulation results on a traffic network show up to 69% reduction in mean absolute error and 61% reduction in root mean square error when compared to CoSTCo.

Tharzeen, Aabila↗

The Squeezed Bispectrum from CHIME H I Emission and Planck Cosmic Microwave Background Lensing: Current Sensitivity and Forecasts

Line intensity mapping using atomic hydrogen (H I) has the potential to efficiently map large volumes of the Universe if the signal can be successfully separated from overwhelmingly bright radio foreground emission. This motivates cross correlations, to ascertain the cosmological nature of measured H I fluctuations, and to study their connections with galaxies and the underlying matter density field. However, these same foregrounds render the cross correlation with projected fields such as the lensing of the cosmic microwave background (CMB) difficult. Indeed, the correlated Fourier modes vary slowly along the line of sight and are thus most contaminated by the smooth-spectrum radio continuum foregrounds. In this paper, we implement a method that avoids this issue by attempting to measure the nonlinear gravitational coupling of the small-scale 21 cm power from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) with large-scale Planck CMB lensing. This measurement is a position-dependent power spectrum, i.e., a squeezed integrated bispectrum. Using 94 nights of CHIME data between 1.0 < z < 1.3 and aggressive foreground filtering, we find that the expected signal is 5 times smaller than the current noise. We forecast that incorporating the additional nights of CHIME data already collected would enable a signal-to-noise ratio of 3, without any further improvements in filtering for foreground cleaning.

21cm↗

Inversion of solar extinction data from the Apollo-Soyuz Test Project Stratospheric Aerosol Measurement (ASTP/SAM) experiment

The inversion methods are reported that have been used to determine the vertical profile of the extinction coefficient due to the stratospheric aerosols from data measured during the ASTP/SAM solar occultation experiment. Inversion methods include the onion skin peel technique and methods of solving the Fredholm equation for the problem subject to smoothing constraints. The latter of these approaches involves a double inversion scheme. Comparisons are made between the inverted results from the SAM experiment and near simultaneous measurements made by lidar and balloon born dustsonde. The results are used to demonstrate the assumptions required to perform the inversions for aerosols.

Pepin, T. J.↗

Uniform high order spectral methods for one and two dimensional Euler equations

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two dimensional shock wave and a rotating vortex is simulated.

Cai, Wei↗

Surface roughness of flat and curved optical surfaces

Surface roughness measurement has several applications. Even a few A roughness will cause scattered light in optical systems. Smooth surfaces are required in a wide variety of instruments. For example, the outputs of the high power lasers are limited by the surface roughness of mirrors and windows. Similarly, the information storage capacity of magnetic media is limited by the roughness of the surface. Roughness reduces the resolving power of optics and distorts images. The performance of certain thin film components in electronic industries is affected by the roughness on the film surface. X-ray astronomical telescopes require smooth curved surfaces. To improve the surface quality, super sensitive detection methods are required. Wide ranging measurement techniques are developed based on interferometry, electron microscopy, C-rays, ellipsometry, light scattering, and using mechanical stylus, etc. Though there are several techniques available for measurement and evaluation of the surfaces, no single technique is fully adequate. Also, the technique used should be nondestructive and highly sensitive. So, an optical heterodyne profilometer was fabricated. Its current sensitivity is much better than 10A rms. It is a noncontact and nondestructive technique. The instrument can be operated even by unskilled personnel for routine measurements.

George, M. C.↗

Uniform high-order spectral methods for one- and two-dimensional Euler equations

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one-dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two-dimensional shock wave and a rotating vortex is simulated.

Cai, Wei↗