Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing technique”

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 361 records · Page 20

Open Specy 1.0: Automated (Hyper)spectroscopy for Microplastics

Microplastic spectral analysis is one of the most time-consuming processes in studying microplastic pollution, often requiring days per sample. Researchers are transitioning to automated batch and hyperspectral image analysis techniques to enhance efficiency. Open Specy, initially aimed at manual single-spectrum analysis, has now integrated automated methods. This updated version, Open Specy 1.0, introduces several new features, including two algorithms for automated processing (smoothing and particle compression), an extensive library containing over 40,000 open-source Raman and FTIR spectra, and two machine learning classifiers (logistic regression and k medoids) developed from this library. Furthermore, it includes a revamped user interface, an R package, and a benchmark data set for testing future advancements in automated techniques. Researchers evaluated various configurations for hyperspectral smoothing, particle identification, compression, and splitting, to achieve combined recovery rates between 50 and 150% particle counts, identities, and sizes with a coefficient of variation (CV) of less than 40% (the accredited standard). Mean absorbance times the standard deviation provided a consistent particle identification. Hyperspectral smoothing led to a 96% combined recovery rate and reduced variability (CV = 38%) compared to the 86% recovery (CV = 83%) of nonsmoothed controls. Additionally, compressing spectra for particles was significantly faster (>3x) and showed similar accuracy but with reduced variability than processing each pixel individually. Key challenges persist in automating spectral analysis, particularly in refining particle splitting algorithms, and improving identification routines to minimize false positives and negatives. In conclusion, new methods in sample preparation for better stabilization and dispersion of particles could overcome some of these issues.

13 HYDRO ENERGY↗

Multi-grid solution of Neumann pressure problem for viscous flows using primitive variables

The multi-grid (MG) technique has been advanced for use with the Neumann boundary-value problem in clustered curvilinear orthogonal coordinates. This comprises an important step in the analysis of viscous flows using the velocity-pressure formulation of the Navier-Stokes equations. With successive over-relaxation (SOR) as the smoothing operator and with suitably formulated restriction and coarse-grid-correction operators, a 4-grid procedure enhances the efficiency of fine-grid solutions of the Neumann problem by a factor of 3 to 5, depending on the problem parameters. Thy influence of the smoothing operator is also examined by employing the alternating-direction implicit and the strongly implicit techniques instead of SOR.

Ghia, U.↗

Development of computer models for the prediction of large distorted antenna characteristics

A program for calculating the radiation patterns of analytic reflector antennas, using geometric optics and aperture plane integration techniques, was modified and extended to include nonanalytic surfaces. The original program, while capable of predicting patterns for a variety of reflectors, both smooth and multipanelled, was subject to the restriction that the surface had to be expressed analytically, i.e., described by a set of mathematical equations. This restriction was removed by adapting a technique which allows the reflector configuration to be interpolated from a finite set of points measured on the surface. Criteria for choosing both the number and distribution of these measured surface points were developed by comparing the predicted antenna radiation patterns with measured laboratory results. Test cases involved relatively smooth as well as grossly distorted arbitrary surfaces. Finally, surface point accuracies were established by applying random error analysis to the measurement process.

Moskowitz, S. M.↗

Detecting hidden transient events in noisy nonlinear time-series

The information impulse function (IIF), running Variance, and local Hölder Exponent are three conceptually different time-series evaluation techniques. These techniques examine time-series for local changes in information content, statistical variation, and point-wise smoothness, respectively. Using simulated data emulating a randomly excited nonlinear dynamical system, this study interrogates the utility of each method to correctly differentiate a transient event from the background while simultaneously locating it in time. Computational experiments are designed and conducted to evaluate the efficacy of each technique by varying pulse size, time location, and noise level in time-series. Our findings reveal that, in most cases, the first instance of a transient event is more easily observed with the information-based approach of IIF than with the Variance and local Hölder Exponent methods. While our study highlights the unique strengths of each technique, the results suggest that very robust and reliable event detection for nonlinear systems producing noisy time-series data can be obtained by incorporating the IIF into the analysis.

42 ENGINEERING↗

COUPLING SMOOTHED PARTICLE HYDRODYNAMICS WITH FINITE ELEMENT METHOD TO SIMULATE RESIDUAL STRESSES FROM FRICTION STIR PROCESSING

Friction stir processing (FSP) is a solid-state material processing technique that locally modifies the microstructure but also induces undesirable residual stresses. A robust numerical model for the FSP can help in mitigating these residual stresses. Heat source models within a finite element method (FEM) framework suffer from inaccuracies. In contrast, smoothed particle hydrodynamics (SPH) model that explicitly captures the material flow near the tool and the associated heat generation are accurate. However, the computational expense of SPH simulations can be prohibitive. In this work, we propose a coupled SPH-FEM framework. SPH is used to model the heat generation accurately near the tool and which is then inserted into to the FEM model as a heat source. To verify this proposed coupling approach, a test case is set up with typical FSP conditions and it is modeled in both SPH and SPH-FEM. The temperatures profiles were compared after the simulations have reached steady-state temperatures. The similarity of the temperature profiles from SPH-FEA and SPH validated the proposed coupling approach. This proposed approach achieves the accuracy of the SPH method while potentially retaining the low computational expense of FEM.

smoothed particle hydrodynamics, Finite Element Me↗

Efficient matrix partitioning for optical computing

Techniques for partitioning optical linear algebra problems to make them amenable to solution using optical processors programmed with simple algorithms are explored. Generalized methods for splitting a linear algebra matrix into a series of submatrices are reviewed, showing that simple forms can be pipelined smoothly and that parallel accumulation can be achieved by beam combining on detectors or by summing electronically. The techniques offer simplified bookkeeping, algorithmic independence, and high efficiency. The computational speed will depend on the number of multiplier-accumulators devoted to the task.

Caulfield, H. J.↗

Landscaper v1

Understanding the inner workings of machine learning models through their loss landscapes offers crucial insights into model properties, optimization dynamics, and generalizability. However, accessing these insights has traditionally required specialized mathematical expertise, limiting broader adoption. Landscaper is an open-source Python package designed to bridge this gap. Landscaper seamlessly integrates a suite of multi-dimensional loss landscape analyses with cutting-edge topological data analysis (TDA) methods. This powerful combination makes both fundamental loss landscape analysis and advanced TDA techniques accessible to the broader scientific ML community, without requiring deep pre-existing mathematical knowledge. Landscaper offers three key functionalities: * Construction: Builds detailed loss landscape representations through versatile low and high-dimensional sampling techniques. * Quantification: Applies advanced metrics, including a novel topological data analysis (TDA) based smoothness metric, enabling new perspectives on model behavior. * Visualization: Offers intuitive tools to visualize and interpret loss landscapes, providing actionable insights beyond traditional performance metrics.

Weber, Gunther [Lawrence Berkeley National Laborat↗

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.↗

Minimal-resource computer program for automatic generation of ocean wave ray or crest diagrams in shoaling waters

A computer program for studying linear ocean wave refraction is described. The program features random-access modular bathymetry data storage. Three bottom topography approximation techniques are available in the program which provide varying degrees of bathymetry data smoothing. Refraction diagrams are generated automatically and can be displayed graphically in three forms: Ray patterns with specified uniform deepwater ray density, ray patterns with controlled nearshore ray density, or crest patterns constructed by using a cubic polynomial to approximate crest segments between adjacent rays.

Poole, L. R.↗

Construction of composite three dimensional grids from subregion grids generated by elliptic systems

A technique is presented for constructing a boundary-conforming grid throughout a general three-dimensional flow region as a composite of subregion grids. Each subregion grid is generated numerically by solving a nonlinear system of elliptic equations. The boundary values represent nodal points in a quasi-two-dimensional grid that covers the curved surface bounding the subregion, and are generated numerically by a modified elliptic system. A primary feature of the technique is that the composite three-dimensional grid remains both continuous and smooth across the surface of juncture between any two adjoining subregions. Sample results are given for a wing-body combination.

Thomas, P. D.↗

Analysis of combustion spectra containing organ pipe tone by cepstral techniques

Signal reinforcements and cancellations due to standing waves may distort constant bandwidth combustion spectra. Cepstral techniques previously applied to the ground reflection echo problem are used to obtain smooth broadband data and information on combustion noise propagation. Internal fluctuating pressure measurements made using a J47 combustor attached to a 6.44 m long duct are analyzed. Measurements made with Jet A and hydrogen fuels are compared. The acoustic power levels inferred from the measurements are presented for a range of low heat release rate operating conditions near atmospheric pressure. For these cases, the variation with operating condition of the overall acoustic broadband power level for both hydrogen and Jet A fuels is consistent with previous results showing it was proportional to the square of the heat release rate. However, the overall acoustic broadband power level generally is greater for hydrogen than for Jet A.

Miles, J. H.↗

System Identification for Nonlinear Control Using Neural Networks

An approach to incorporating artificial neural networks in nonlinear, adaptive control systems is described. The controller contains three principal elements: a nonlinear inverse dynamic control law whose coefficients depend on a comprehensive model of the plant, a neural network that models system dynamics, and a state estimator whose outputs drive the control law and train the neural network. Attention is focused on the system identification task, which combines an extended Kalman filter with generalized spline function approximation. Continual learning is possible during normal operation, without taking the system off line for specialized training. Nonlinear inverse dynamic control requires smooth derivatives as well as function estimates, imposing stringent goals on the approximating technique.

Stengel, Robert F.↗

Full resolution hologram-like autostereoscopic display

Under this program, Dimension Technologies Inc. (DTI) developed a prototype display that uses a proprietary illumination technique to create autostereoscopic hologram-like full resolution images on an LCD operating at 180 fps. The resulting 3D image possesses a resolution equal to that of the LCD along with properties normally associated with holograms, including change of perspective with observer position and lack of viewing position restrictions. Furthermore, this autostereoscopic technique eliminates the need to wear special glasses to achieve the parallax effect. Under the program a prototype display was developed which demonstrates the hologram-like full resolution concept. To implement such a system, DTI explored various concept designs and enabling technologies required to support those designs. Specifically required were: a parallax illumination system with sufficient brightness and control; an LCD with rapid address and pixel response; and an interface to an image generation system for creation of computer graphics. Of the possible parallax illumination system designs, we chose a design which utilizes an array of fluorescent lamps. This system creates six sets of illumination areas to be imaged behind an LCD. This controlled illumination array is interfaced to a lenticular lens assembly which images the light segments into thin vertical light lines to achieve the parallax effect. This light line formation is the foundation of DTI's autostereoscopic technique. The David Sarnoff Research Center (Sarnoff) was subcontracted to develop an LCD that would operate with a fast scan rate and pixel response. Sarnoff chose a surface mode cell technique and produced the world's first large area pi-cell active matrix TFT LCD. The device provided adequate performance to evaluate five different perspective stereo viewing zones. A Silicon Graphics' Iris Indigo system was used for image generation which allowed for static and dynamic multiple perspective image rendering. During the development of the prototype display, we identified many critical issues associated with implementing such a technology. Testing and evaluation enabled us to prove that this illumination technique provides autostereoscopic 3D multi perspective images with a wide range of view, smooth transition, and flickerless operation given suitable enabling technologies.

Eichenlaub, Jesse B.↗

Improving Global Mass Flux Solutions from Gravity Recovery and Climate Experiment (GRACE) Through Forward Modeling and Continuous Time Correlation

We describe Earth's mass flux from April 2003 through November 2008 by deriving a time series of mas cons on a global 2deg x 2deg equal-area grid at 10 day intervals. We estimate the mass flux directly from K band range rate (KBRR) data provided by the Gravity Recovery and Climate Experiment (GRACE) mission. Using regularized least squares, we take into account the underlying process dynamics through continuous space and time-correlated constraints. In addition, we place the mascon approach in the context of other filtering techniques, showing its equivalence to anisotropic, nonsymmetric filtering, least squares collocation, and Kalman smoothing. We produce mascon time series from KBRR data that have and have not been corrected (forward modeled) for hydrological processes and fmd that the former produce superior results in oceanic areas by minimizing signal leakage from strong sources on land. By exploiting the structure of the spatiotemporal constraints, we are able to use a much more efficient (in storage and computation) inversion algorithm based upon the conjugate gradient method. This allows us to apply continuous rather than piecewise continuous time-correlated constraints, which we show via global maps and comparisons with ocean-bottom pressure gauges, to produce time series with reduced random variance and full systematic signal. Finally, we present a preferred global model, a hybrid whose oceanic portions are derived using forward modeling of hydrology but whose land portions are not, and thus represent a pure GRACE-derived signal.

Sabaka, T. J.↗

Waves in Radial Gravity Using Magnetic Fluid

Terrestrial laboratory experiments studying various fluid dynamical processes are constrained, by being in an Earth laboratory, to have a gravitational body force which is uniform and unidirectional. Therefore fluid free-surfaces are horizontal and flat. Such free surfaces must have a vertical solid boundary to keep the fluid from spreading horizontally along a gravitational potential surface. In atmospheric, oceanic, or stellar fluid flows that have a horizontal scale of about one-tenth the body radius or larger, sphericity is important in the dynamics. Further, fluids in spherical geometry can cover an entire domain without any sidewall effects, i.e. have truly periodic boundary conditions. We describe spherical body-force laboratory experiments using ferrofluid. Ferrofluids are dilute suspensions of magnetic dipoles, for example magnetite particles of order 10 nm diameter, suspended in a carrier fluid. Ferrofluids are subject to an additional body force in the presence of an applied magnetic field gradient. We use this body force to conduct laboratory experiments in spherical geometry. The present study is a laboratory technique improvement. The apparatus is cylindrically axisymmetric. A cylindrical ceramic magnet is embedded in a smooth, solid, spherical PVC ball. The geopotential field and its gradient, the body force, were made nearly spherical by careful choice of magnet height-to-diameter ratio and magnet size relative to the PVC ball size. Terrestrial gravity is eliminated from the dynamics by immersing the "planet" and its ferrofluid "ocean" in an immiscible silicone oil/freon mixture of the same density. Thus the earth gravity is removed from the dynamics of the ferrofluid/oil interface and the only dynamically active force there is the radial magnetic gravity. The entire apparatus can rotate, and waves are forced on the ferrofluid surface by exterior magnets. The biggest improvement in technique is in the wave visualization. Fluorescing dye is added to the oil/freon mixture and an argon ion laser generates a horizontal light that can be scanned vertically. Viewed from above, the experiment is a black circle with wave deformations surrounded by a light background. A contour of the image intensity at any light sheet position gives the surface of the ferrofluid "ocean" at that "latitude". Radial displacements of the waves as a function of longitude are obtained by subtracting the contour line positions from a no-motion contour at that laser sheet latitude. The experiments are run by traversing the forcing magnet with the laser sheet height fixed and images are frame grabbed to obtain a time-series at one latitude. The experiment is then re-run with another laser-sheet height to generate a full picture of the three-dimensional wave structure in the upper hemisphere of the ball as a function of time. We concentrate here on results of laboratory studies of waves that are important in Earth's atmosphere and especially the ocean. To get oceanic scaling in the laboratory, the experiment must rotate rapidly (4-second rotation period) so that the wave speed is slow compared to the planetary rotation speed as in the ocean. In the Pacific Ocean, eastward propagating Kelvin waves eventually run into the South American coast. Theory predicts that some of the wave energy should scatter into coastal-trapped Kelvin waves that propagate north and south along the coast. Some of this coastal wave energy might then scatter into mid-latitude Rossby waves that propagate back westward. Satellite observations of the Pacific Ocean sea-surface temperature and height seem to show signatures of westward propagating mid-latitude Rossby waves, 5 to 10 years after the 1982-83 El Nino. The observational data is difficult to interpret unambiguously owing to the large range of motions that fill the ocean at shorter timescales. This series of reflections giving eastward, north- ward, and then westward traveling waves is observed cleanly in the laboratory experiments, confirming the theoretical expectations.

Ohlsen, D. R.↗

Optimization of self-acting herringbone journal bearings for maximum stability

The stability of herringbone grooved journal bearings was increased by optimizing the groove depth, width, length, and angle. Optimization was performed on the basis of a small-eccentricity infinite-groove analysis in combination with a previously developed Newton-Raphson procedure for bearings with the smooth member rotating or with the grooved member rotating at low compressibility numbers. A newly-developed vector technique for bearings with the grooved member rotating at high compressibility numbers also was used in the analysis. The design curves obtained are suitable for choosing the optimum bearing for a wide range of operating conditions (compressibility numbers from 0 to 80, length to diameter ratios from 1/4 to 2, and smooth or grooved member rotating). Compared to bearings optimized for maximum load, bearings optimized for maximum stability are shown to provide a drastic increase in bearing-supported mass. Bearings with grooved member rotating are substantially more stable than bearing with smooth member rotating, particularly at high compressibility numbers.

Fleming, D. P.↗

The Surface-Topography Challenge: A Multi-Laboratory Benchmark Study to Advance the Characterization of Topography

Surface performance is critically influenced by topography in virtually all real-world applications. The current standard practice is to describe topography using one of a few industry-standard parameters. The most commonly reported number is Ra, the average absolute deviation of the height from the mean line (at some, not necessarily known or specified, lateral length scale). However, other parameters, particularly those that are scale-dependent, influence surface and interfacial properties; for example the local surface slope is critical for visual appearance, friction, and wear. The present Surface-Topography Challenge was launched to raise awareness for the need of a multi-scale description, but also to assess the reliability of different metrology techniques. In the resulting international collaborative effort, 153 scientists and engineers from 64 research groups and companies across 20 countries characterized statistically equivalent samples from two different surfaces: a “rough” and a “smooth” surface. The results of the 2088 measurements constitute the most comprehensive surface description ever compiled. We find wide disagreement across measurements and techniques when the lateral scale of the measurement is ignored. Consensus is established through scale-dependent parameters while removing data that violates an established resolution criterion and deviates from the majority measurements at each length scale. Our findings suggest best practices for characterizing and specifying topography. The public release of the accumulated data and presented analyses enables global reuse for further scientific investigation and benchmarking.

42 ENGINEERING↗

Correlated Trajectory Uncertainty for Adaptive Sequential Decision Making

One of the great challenges with decision making tasks on real world systems is the fact that data is sparse and acquiring additional data is expensive. In these cases, it is often crucial to make a model of the environment to assist in making decisions. At the same time, limited data means that learned models are erroneous, making it just as important to equip the model with good predictive uncertainties. In the context of learning sequential decision making policies, these uncertainties can prove useful for informing which data to collect for the greatest improvement in policy performance \citep{mehta2021experimental, mehta2022exploration} or informing the policy about unsure regions of state and action space to avoid during test time \citep{yu2020mopo}. Additionally, assuming that realistic samples of the environment can be drawn, an adaptable policy can be trained that attempts to make optimal decisions for any given possible instance of the environment \citep{ghosh2022offline, chen2021offline}. In this work, we examine the so-called ``probabilistic neural network'' (PNN) model that is ubiquitous in model-based reinforcement learning (MBRL) works. We argue that while PNN models may have good marginal uncertainties, they form a distribution of non-smooth transition functions. Not only are these samples unrealistic and may hamper adaptability, but we also assert that this leads to poor uncertainty estimates when predicting multiple step trajectory estimates. To address this issue, we propose a simple sampling method that can be implemented on top of pre-existing models.We evaluate our sampling technique on a number of environments, including a realistic nuclear fusion task, and find that, not only do smooth transition function samples produce more calibrated uncertainties, but they also lead to better downstream performance for an adaptive policy.

Offline Reinforcement Learning↗