Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing techniques”

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 451 records · Page 25

An adaptive grid, unsteady model for two-dimensional magnetohydrodynamic (MHD) flow

An adaptive grid finite difference method for solving multi-dimensional, time-dependent, magnetohydrodynamic (MHD) equations is developed. The method is capable of solving problems that include high gradients due to geometry, propagation of shock waves and unsteady boundary conditions. The grid generation technique is based on variational principles with direct control over grid concentration, smoothness and skewness. An example for a two-dimensional MHD simulation of the propagation of a solar-flare-generated shock wave in solar wind flow in the heliographic equatorial plane is selected for illustration of this method.

Panitchob, Supat↗

A generalized gradient algorithm for dynamic optimization

A gradient algorithm is developed that determines optimal trajectories with path equality constraints and terminal constraints. A generalized gradient is formed which improves both the performance index and the path equality constraints simultaneously. The algorithm is extended to treat terminal constraints by using Bryson's impulse response technique. The main features of this algorithm are its numerical stability and smooth convergence near the optimum.

Zhao, Yiyuan↗

A multigrid nonoscillatory method for computing high speed flows

A multigrid method using different smoothers has been developed to solve the Euler equations discretized by a nonoscillatory scheme up to fourth order accuracy. The best smoothing property is provided by a five-stage Runge-Kutta technique with optimized coefficients, yet the most efficient smoother is a backward Euler technique in factored and diagonalized form. The singlegrid solution for a hypersonic, viscous conic flow is in excellent agreement with the solution obtained by the third order MUSCL and Roe's method. Mach 8 inviscid flow computations for a complete entry probe have shown that the accuracy is at least as good as the symmetric TVD scheme of Yee and Harten. The implicit multigrid method is four times more efficient than the explicit multigrid technique and 3.5 times faster than the single-grid implicit technique. For a Mach 8.7 inviscid flow over a blunt delta wing at 30 deg incidence, the CPU reduction factor from the three-level multigrid computation is 2.2 on a grid of 37 x 41 x 73 nodes.

Li, C. P.↗

Constrained optimization using design of experiment surfaces

An algorithm for solving constrained optimization problems is presented. First, design of experiment techniques are used to survey the design space. After evaluating the objective and constraint functions, as specified by Taguchi orthogonal arrays, analytical models of these functions are generated using a least-squares regression analysis. Next, a nonlinear programming package is used to optimize the analytical model. Based on the optimization information, the design space is reduced so as to close in around the minimum, and the entire procedure is repeated until convergence. An important feature of the algorithm is that function gradients are not required; therefore, for problems in which gradients would have to be estimated using finite-differences the number of function evaluations required for the optimization is significantly reduced, when compared with traditional nonlinear programming techniques. In addition, there is no requirement that the gradients must be smooth and continuous.

Bolt, Marvin Vance↗

A Numerical Simulation of Scattering from One-Dimensional Inhomogeneous Dielectric Random Surfaces

In this paper, an efficient numerical solution for the scattering problem of inhomogeneous dielectric rough surfaces is presented. The inhomogeneous dielectric random surface represents a bare soil surface and is considered to be comprised of a large number of randomly positioned dielectric humps of different sizes, shapes, and dielectric constants above an impedance surface. Clods with nonuniform moisture content and rocks are modeled by inhomogeneous dielectric humps and the underlying smooth wet soil surface is modeled by an impedance surface. In this technique, an efficient numerical solution for the constituent dielectric humps over an impedance surface is obtained using Green's function derived by the exact image theory in conjunction with the method of moments. The scattered field from a sample of the rough surface is obtained by summing the scattered fields from all the individual humps of the surface coherently ignoring the effect of multiple scattering between the humps. The statistical behavior of the scattering coefficient sigma(sup 0) is obtained from the calculation of scattered fields of many different realizations of the surface. Numerical results are presented for several different roughnesses and dielectric constants of the random surfaces. The numerical technique is verified by comparing the numerical solution with the solution based on the small perturbation method and the physical optics model for homogeneous rough surfaces. This technique can be used to study the behavior of scattering coefficient and phase difference statistics of rough soil surfaces for which no analytical solution exists.

Sarabandi, Kamal↗

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

Residue-free fabrication of high-performance graphene devices by patterned PMMA stencil mask

Two-dimensional (2D) atomic crystals and their hybrid structures have recently attracted much attention due to their potential applications. The fabrication of metallic contacts or nanostructures on 2D materials is very common and generally achieved by performing electron-beam (e-beam) lithography. However, e-beam lithography is not applicable in certain situations, e.g., cases in which the e-beam resist does not adhere to the substrates or the intrinsic properties of the 2D materials are greatly altered and degraded. Here, we present a residue-free approach for fabricating high-performance graphene devices by patterning a thin film of e-beam resist as a stencil mask. This technique can be generally applied to substrates with varying surface conditions, while causing negligible residues on graphene. The technique also preserves the design flexibility offered by e-beam lithography and therefore allows us to fabricate multi-probe metallic contacts. The graphene field-effect transistors fabricated by this method exhibit smooth surfaces, high mobility, and distinct magnetotransport properties, confirming the advantages and versatility of the presented residue-free technique for the fabrication of devices composed of 2D materials.

36 MATERIALS SCIENCE↗

The Center for Momentum Transport and Flow Organization in Plasmas - Final Scientific Report

Overview of University of Colorado Efforts: The University of Colorado group has focused on two primary fronts during the grant period: development of a variety of multi-point diagnostic and/or imaging analysis techniques, and momentum-transport related experiments on a variety of devices (NSTX at PPPL, CSDX at UCSD, LAPD at UCLA, DIII-D at GA). Experimental work has taken advantage of several diagnostic instruments, including fast-framing cameras for imaging of electron density fluctuations (either directly or using injected gas puffs), ECEI for imaging of electron temperature fluctuations, and multi-tipped Langmuir and magnetic probes for corroborating measurements of Reynolds and Maxwell stresses. Mode Characterization in CSDX: We have performed a series of experiments at the CSDX linear device at UCSD, in collaboration with Center PI G. Tynan's group. The experiments included a detailed study of velocity estimation techniques, including direct comparisons between Langmuir probes and image-based velocimetry from fast-framing camera data. We used the camera data in a second set of studies to identify the spatial and spectral structure of coherent modes, which illuminates wave behavior to a level of detail previously unavailable, and enables direct comparison of dispersion curves to theoretical estimates. In another CSDX study, similar techniques were used to demonstrate a controlled transition from nonlinearly coupled discrete eigenmodes to fully developed broadband turbulence. The axial magnetic field was varied from 40-240 mT, which drove the transition. At low magnetic fields, the plasma is dominated by drift waves. As the magnetic field is increased, a strong potential gradient at the edge introduces an ExB shear-driven instability. At the transition, another mode with signatures of a rotation-induced Rayleigh–Taylor instability appears at the central plasma region. Concurrently, large axial velocities were found in the plasma core. For larger magnetic fields, all of the instabilities co-exist, leading to rich plasma dynamics and fully developed broadband turbulence. Edge-Turbulence and Flow Experiments in NSTX: A series of Gas Puff Imaging (GPI) observations on NSTX revealed a quasi-periodic oscillation in the plasma edge preceding the L-H transition in a limited set of neutral beam heated plasmas. These ~3 kHz flow oscillations exhibit both long wavelength and long correlation lengths, suggesting they are zonal-flow-like. The flow oscillations are strongly correlated with modulations of the level of edge turbulence, thus the system appears to undergo a predator--prey-type limit-cycle preceding the L-H transition. However, a clear trigger for the L-H transition was not observed. Reynolds stress profiles were obtained directly from image velocimetry for L-mode periods ELM-Precursor Studies in NSTX: A separate study based on NSTX-GPI data captured the two-dimensional evolution of edge-localized mode (ELM) precursors. Precursor events were observed preceding ELMs and ELM-induced H–L back-transitions in radio-frequency heated H-mode plasmas, and the growth of the precursor mode through the ELM filamentation was imaged in the plane perpendicular to the local B-field. Strong edge intensity modulations appeared to propagate in the electron diamagnetic direction while steadily drifting radially outwards. Intensity fluctuations were observed at frequencies around 20 kHz and wavenumbers of 0.05-0.2 cm -1 . Upon growing to a trigger point, precursor fluctuations were seen to form filamentary structures and move into the scrape-off layer (SOL) explosively with radial velocities peaking at 8 km/s. Once in the SOL, filaments reverse their propagation direction and travel in the ion diamagnetic direction. Edge intensity fluctuations were strongly correlated with magnetic signals from Mirnov coils, and toroidally distributed coils estimated toroidal mode numbers of n=5-10. Quantitatively similar precursors have been observed in ohmic H-mode plasmas as well, though significantly fewer events are seen in the ohmic cases and none were observed in the near-threshold NBI H-modes studied. Development of Velocimetry and Image Analysis Techniques: Along with the experiments listed above, the Colorado group has continued to explore various velocimetry techniques and their range of validity. We have developed a “linear optical flow” code, which calculates smooth velocity maps while accurately assessing local regions of high curl. This is critical for separating spatial scales of velocity behavior, and thus transport. This code has recently come on-line, and we are currently using it to revisit a number of older datasets. Additionally, we have worked on developing pattern-recognition techniques for imaging diagnostics, based on established digital image compression algorithms. This has the potential to open the analysis of turbulent plasma behavior beyond the well-trodden Fourier and wavelet approaches. Finally, we have extended several of these image-analysis routines to multiple other diagnostic sets, including GPI datasets in NSTX and ECEI imaging on DIII-D. In one study, correlation maps were used to measure the full 2-D mode structure of drift-wave level structures using ECEI for the first time. This enabled direct comparison to gyrokinetic simulations using the GEM code, from which is was determined that TEM modes were being measured in DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ideal-observer detectability in photon-counting differential phase-contrast imaging using a linear-systems approach

Purpose: To provide a cascaded-systems framework based on the noise-power spectrum (NPS), modulation transfer function (MTF), and noise-equivalent number of quanta (NEQ) for quantitative evaluation of differential phase-contrast imaging (Talbot interferometry) in relation to conventional absorption contrast under equal-dose, equal-geometry, and, to some extent, equal-photon-economy constraints. The focus is a geometry for photon-counting mammography. Methods: Phase-contrast imaging is a promising technology that may emerge as an alternative or adjunct to conventional absorption contrast. In particular, phase contrast may increase the signal-difference-to-noise ratio compared to absorption contrast because the difference in phase shift between soft-tissue structures is often substantially larger than the absorption difference. We have developed a comprehensive cascaded-systems framework to investigate Talbot interferometry, which is a technique for differential phase-contrast imaging. Analytical expressions for the MTF and NPS were derived to calculate the NEQ and a task-specific ideal-observer detectability index under assumptions of linearity and shift invariance. Talbot interferometry was compared to absorption contrast at equal dose, and using either a plane wave or a spherical wave in a conceivable mammography geometry. The impact of source size and spectrum bandwidth was included in the framework, and the trade-off with photon economy was investigated in some detail. Wave-propagation simulations were used to verify the analytical expressions and to generate example images. Results: Talbot interferometry inherently detects the differential of the phase, which led to a maximum in NEQ at high spatial frequencies, whereas the absorption-contrast NEQ decreased monotonically with frequency. Further, phase contrast detects differences in density rather than atomic number, and the optimal imaging energy was found to be a factor of 1.7 higher than for absorption contrast. Talbot interferometry with a plane wave increased detectability for 0.1-mm tumor and glandular structures by a factor of 3-4 at equal dose, whereas absorption contrast was the preferred method for structures larger than {approx}0.5 mm. Microcalcifications are small, but differ from soft tissue in atomic number more than density, which is favored by absorption contrast, and Talbot interferometry was barely beneficial at all within the resolution limit of the system. Further, Talbot interferometry favored detection of 'sharp' as opposed to 'smooth' structures, and discrimination tasks by about 50% compared to detection tasks. The technique was relatively insensitive to spectrum bandwidth, whereas the projected source size was more important. If equal photon economy was added as a restriction, phase-contrast efficiency was reduced so that the benefit for detection tasks almost vanished compared to absorption contrast, but discrimination tasks were still improved close to a factor of 2 at the resolution limit. Conclusions: Cascaded-systems analysis enables comprehensive and intuitive evaluation of phase-contrast efficiency in relation to absorption contrast under requirements of equal dose, equal geometry, and equal photon economy. The benefit of Talbot interferometry was highly dependent on task, in particular detection versus discrimination tasks, and target size, shape, and material. Requiring equal photon economy weakened the benefit of Talbot interferometry in mammography.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

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

Two-step fabrication technique of gold tips for use in point-contact spectroscopy

For a successful point-contact spectroscopy (PCS) measurement, metallic tips of proper shape and smoothness are essential to ensure the ballistic nature of a point-contact junction. Until recently, the fabrication of Au tips suitable for use in point-contact spectroscopy has remained more of an art involving a trial and error method rather than an automated scientific process. To address these issues, we have developed a technique with which one can prepare high quality Au tips reproducibly and systematically. It involves an electronic control of the driving voltages used for an electrochemical etching of a gold wire in a HCl-glycerol mixture or a HCl solution. We find that a stopping current, below which the circuit is set to shut off, is a single very important parameter to produce an Au tip of desired shape. We present detailed descriptions for a two-step etching process for Au tips and also test results from PCS measurements using them.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

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

Simple Hydrothermal Synthesis for Mesoporous Pebble-Like NiCo{sub 2}O{sub 4} Nanostructures and Their Attractive Magnetic Properties

The mesoporous pebble-like NiCo{sub 2}O{sub 4} nanostructures were smoothly obtained via a simple hydrothermal route and the succeeding annealing treatment. Various techniques, including XRD, SEM, TEM, Raman spectra, FTIR, and MH curves have been used to study the obtained products. The obtained results exhibited that the morphology of spinel NiCo{sub 2}O{sub 4} could be controlled by adjusting the content of cetyltrimethylammonium bromide (CTAB). The morphology of NiCo{sub 2}O{sub 4} nanostructures transformed from a 300-nm popcorn-like to 30–40-nm pebble-like nanostructure by changing the content of CTAB. Compared with the popcorn-like NiCo{sub 2}O{sub 4} nanostructures, the pebble-like samples showed effective blue shifting of Raman spectroscopy, which might be due to the nano-sized effect of NiCo{sub 2}O{sub 4} particles. Finally, vibration sample magnetometer (VSM) measurements indicated the samples transformed from weak ferromagnetic (popcorn-like) to paramagnetic (pebble-like) due to small particle size.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A reduced basis approach for PDEs on parametrized geometries based on the shifted boundary finite element method and application to a Stokes flow

Here we propose a model order reduction technique integrating the Shifted Boundary Method (SBM) with a POD-Galerkin strategy. This approach allows to deal with complex parametrized domains in an efficient and straightforward way. The impact of the proposed approach is threefold. First, problems involving parametrizations of complex geometrical shapes and/or large domain deformations can be efficiently solved at full-order by means of the SBM. This unfitted boundary method permits to avoid remeshing and the tedious handling of cut cells by introducing an approximate surrogate boundary. Second, the computational effort is reduced by the development of a Reduced Order Model (ROM) technique based on a POD-Galerkin approach. Third, the SBM provides a smooth mapping from the true to the surrogate domain, and for this reason, the stability and performance of the reduced order basis are enhanced. This feature is the net result of the combination of the proposed ROM approach and the SBM. Similarly, the combination of the SBM with a projection-based ROM gives the great advantage of an easy and fast to implement algorithm considering geometrical parametrization with large deformations. The transformation of each geometry to a reference geometry (morphing) is in fact not required. These combined advantages will allow the solution of PDE problems more efficiently. We illustrate the performance of this approach on a number of two-dimensional Stokes flow problems.

42 ENGINEERING↗

Connecting the dots: Tracking galaxy evolution using constant cumulative number density at 3 ≤ z ≤ 7

Using the cosmological smoothed particle hydrodynamical code GADGET-3, we make a realistic assessment of the technique of using constant cumulative number density as a tracer of galaxy evolution at high redshift. We find that over a redshift range of 3 ≤ z ≤ 7 one can on average track the growth of the stellar mass of a population of galaxies selected from the same cumulative number density bin to within ∼0.20 dex. Over the stellar mass range that we probe (10{sup 10.4}⩽M{sub s}/M{sub ⊙}⩽10{sup 10.8} at z = 3 and 10{sup 8.5}⩽M{sub s}/M{sub ⊙}⩽10{sup 9.6} at z = 7), one can reduce this bias by selecting galaxies based on an evolving cumulative number density. We find that this cumulative number density evolution exhibits a trend toward higher values which can be quantified by simple linear formulations of −0.10Δz for descendants and 0.12Δz for progenitors. Utilizing such an evolving cumulative number density increases the accuracy of descendant/progenitor tracking by a factor of ∼2. This result is in excellent agreement, within 0.10 dex, with abundance matching results over the same redshift range. However, we find that our more physically realistic cosmological hydrodynamic simulations produce a much larger scatter in descendant/progenitor stellar masses than previous studies, particularly when tracking progenitors. This large scatter makes the application of either the constant cumulative number density or evolving cumulative number density technique limited to average stellar masses of populations only, as the diverse mass assembly histories caused by stochastic physical processes such as gas accretion or mergers lead to an even larger scatter in other physical properties such as metallicity and star formation rate.

79 ASTRONOMY AND ASTROPHYSICS↗