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 739 records · Page 41

Measurement of Convective Heat Transfer Coefficients with Supercritical CO2 in Novel Additively Manufactured Helically Patterned Pin Fin Tubes Using the Wilson Plot Technique

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO2 flowing in pin-fin patterned pipes in the Heat Exchange and Experimental Testing (HEET) facility at the US DoE’s National Energy Technology Laboratory (NETL) in Morgantown, WV. The measurement procedures in the HEET rig were validated by conducting benchmark tests with smooth stainless-steel tube and comparing the results with published correlations for Nusselt number (Nu) and friction factor. Over typical Reynolds number range in sCO2, the measured Nu and friction factors were within 7% of classical correlations for smooth tube flow.The candidate pin fin patterned pipes were additively manufactured (AM) at the Oak Ridge National Laboratory. The pins were circular or elliptical in cross-section. Pin length to diameter aspect ratios were 1.33 and 2, while the pin diameter to tube diameter ratio was 0.188 and 0.125. Tests were performed for ReD varying from 6.9×104 to 2.2×105 and at conditions equivalent to the low pressure outlet (8.69 MPa, 361 K) of the low temperature recuperator (LTR) in an indirect sCO2 power cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients.For the better performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.188), the local heat transfer coefficient increased by 112% relative to the Dittus-Boelter correlation at the LTR low pressure outlet. This corresponded to a 282% increase in the product of the heat transfer coefficient and the surface area. Large pressure drops across the test articles were observed.

Searle, Matthew↗

11-th order of accuracy for numerical solution of 3-D Poisson equation with irregular interfaces on unfitted Cartesian meshes

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10 6 times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

97 MATHEMATICS AND COMPUTING↗

Estimate the Electromechanical States Using Particle Filtering and Smoothing

Accurate knowledge of electromechanical states is critical for efficient and reliable control of a power system. This paper proposes a particle filtering approach to estimate the electromechanical states of power systems from Phasor Measurement Unit (PMU) data. Without having to go through laborious linearization procedure, the proposed particle filtering techniques can estimate states of a complex power system, which is often non-linear and has non-Gaussian noise. The proposed method is evaluated using a multi-machine system with both large and small disturbances. Sensitivity studies of the dynamic state estimation performance are also presented to show the robustness of the proposed method. The inherent decoupling properties of particle filtering make it highly scalable and the potential to reduce computational time through parallel implementation is very promising.

Auxiliary Particle Filter↗

Bootstrap AMG for spectral clustering

Graph Laplacian is a popular tool for analyzing graphs, particularly in graph partitioning and clustering. Given a notion of similarity (via an adjacency matrix), graph clustering refers to identifying different groups such that vertices in the same group are more similar compared to vertices across different groups. Data clustering can be reformulated in terms of a graph clustering problem when the given set of data is represented as a graph, also known as similarity graph. In this context, eigenvectors of the graph Laplacian are often used to obtain a new geometric representation of the original data set that generally enhances cluster properties and improves cluster detection. Here, we apply a bootstrap algebraic multigrid (AMG) method that constructs a set of vectors associated with the graph Laplacian. These vectors, referred to as algebraically smooth ones, span a low-dimensional Euclidean space, which we use to represent the data, enabling cluster detection both in synthetic and in realistic well-clustered graphs. We show that, in the case of a good quality bootstrap AMG, the computed smooth vectors employed in the construction of the final AMG operator, which by construction is spectrally equivalent to the originally given graph Laplacian, accurately approximate the space in the lower portion of the spectrum of the preconditioned operator. Thus, our approach can be viewed as a spectral clustering technique associated with the generalized spectral problem (Laplace operator versus the final AMG operator), and hence, it can be seen as an extension of the classical spectral clustering that employs a standard eigenvalue problem.

97 MATHEMATICS AND COMPUTING↗

A variant nonmonotone smoothing algorithm with improved numerical results for large-scale LWCPs

This paper considers the linear weighted complementarity problem (denoted by LWCP) introduced by Potra (SIAM J Optim 22:1634–1654, 2012). Based on two weighted smoothing functions, we propose a new nonmonotone smoothing algorithm for solving the LWCP and establish its global and local quadratic convergence without the strict complementarity assumption. Compared to existing nonmonotone smoothing algorithms, the proposed algorithm solves the linear system only approximately which can save the computation work when one solves large-scale LWCPs. Moreover, the nonmonotone line search technique adopted in this paper includes the usual monotone line search and some existing nonmonotone line searches as special cases. Numerical results show that our algorithm is considerably efficient for solving large-scale LWCPs.

97 MATHEMATICS AND COMPUTING↗

Electronic structure of epitaxial half-Heusler Co 1-x Ni x TiSb across the semiconductor to metal transition

In this paper, we report on the growth, electronic, and surface properties of the electron-doped half-Heusler series Co 1-x Ni x TiSb (001) grown by molecular beam epitaxy. High-quality epitaxial growth of thin films is achieved on InP (001) substrates using an InAlAs buffer layer for all nickel concentrations. The semiconductor to metal transition as a function of substitutional alloying was examined using electrical transport, Seebeck measurements, and angle-resolved photoemission spectroscopy (ARPES). Temperature-dependent electrical transport measurements of films with composition x ≤ 0.1 exhibit thermally activated behavior while x > 0.1 exhibit metallic behavior. Smooth, highly ordered film surfaces can be achieved following ex-situ transfer of the films and subsequent desorption of a sacrificial, protective antimony capping layer. Using this transfer technique, ARPES experiments were performed to investigate the effects of nickel alloying on the electronic band structure. Finally, an electron pocket is observed below the Fermi level at the bulk X point for compositions x > 0.1, in accordance with the crossover from semiconducting to metallic behavior observed in the transport measurements.

36 MATERIALS SCIENCE↗

Detected Changes in Precipitation Extremes at Their Native Scales Derived from In Situ Measurements

Abstract The gridding of daily accumulated precipitation—especially extremes—from ground-based station observations is problematic due to the fractal nature of precipitation, and therefore estimates of long period return values and their changes based on such gridded daily datasets are generally underestimated. In this paper, we characterize high-resolution changes in observed extreme precipitation from 1950 to 2017 for the contiguous United States (CONUS) based on in situ measurements only. Our analysis utilizes spatial statistical methods that allow us to derive gridded estimates that do not smooth extreme daily measurements and are consistent with statistics from the original station data while increasing the resulting signal-to-noise ratio. Furthermore, we use a robust statistical technique to identify significant pointwise changes in the climatology of extreme precipitation while carefully controlling the rate of false positives. We present and discuss seasonal changes in the statistics of extreme precipitation: the largest and most spatially coherent pointwise changes are in fall (SON), with approximately 33% of CONUS exhibiting significant changes (in an absolute sense). Other seasons display very few meaningful pointwise changes (in either a relative or absolute sense), illustrating the difficulty in detecting pointwise changes in extreme precipitation based on in situ measurements. While our main result involves seasonal changes, we also present and discuss annual changes in the statistics of extreme precipitation. In this paper we only seek to detect changes over time and leave attribution of the underlying causes of these changes for future work.

54 ENVIRONMENTAL SCIENCES↗

Titanium reinforced boron polyimide composite

Program involves development of process technique for boron-polyimide prepeg, lay-up and curing procedures for prepegs when processed under vacuum bag pressure, and development and evaluation of titanium hard points for smooth transition of loads from titanium attach points into boron reinforced body of structure.

Clark, G. A.↗

Pole-strength of the earth from Magsat and magnetic determination of the core radius

A model based on two days of Magsat data is used to numerically evaluate the unsigned magnetic flux linking the earth's surface, and a comparison of the 16.054 GWb value calculated with values from earlier geomagnetic field models reveals a smooth, monotonic, and recently-accelerating decrease in the earth's pole strength at a 50-year average rate of 8.3 MWb, or 0.052%/year. Hide's (1978) magnetic technique for determining the radius of the earth's electrically-conducting core is tested by (1) extrapolating main field models for 1960 and 1965 downward through the nearly-insulating mantle, and then separately comparing them to equivalent, extrapolated models of Magsat data. The two unsigned fluxes are found to equal the Magsat values at a radius which is within 2% of the core radius; and (2) the 1960 main field and secular variation and acceleration coefficients are used to derive models of 1930, 1940 and 1950. The same core magnetic radius value, within 2% of the seismic value, is obtained. It is concluded that the mantle is a nearly-perfect insulator, while the core is a perfect conductor, on the decade time scale.

Voorhies, G. V.↗

Energy functions for regularization algorithms

Regularization techniques are widely used for inverse problem solving in computer vision such as surface reconstruction, edge detection, or optical flow estimation. Energy functions used for regularization algorithms measure how smooth a curve or surface is, and to render acceptable solutions these energies must verify certain properties such as invariance with Euclidean transformations or invariance with parameterization. The notion of smoothness energy is extended here to the notion of a differential stabilizer, and it is shown that to void the systematic underestimation of undercurvature for planar curve fitting, it is necessary that circles be the curves of maximum smoothness. A set of stabilizers is proposed that meet this condition as well as invariance with rotation and parameterization.

Delingette, H.↗

Investigations of Remote Plasma Irregularites by Radio Sounding: Applications of the Radio Plasma Imager on IMAGE

The Radio Plasma Imager (RPI) on the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) mission operates like a radar by transmitting and receiving coherent electromagnetic pulses. Long-range echoes of electromagnetic sounder waves are reflected at remote plasma cutoffs. Thus, analyses of RPI observations will yield the plasma parameters and distances to the remote reflection points. These analyses assume that the reflecting plasma surfaces are cold and are sufficiently smooth that they effectively behave as plane mirrors to the incoming sounder waves, i.e., that geometric optics can be used. The RPI will employ pulse compression and spectral integration techniques, perfected in ground-based ionospheric digital sounders, in order to enhance the signal-to-noise ratio in long-range magnetospheric sounding. When plasma irregularities exist in the remote magnetospheric plasmas that are being probed by the sounder waves, echo signatures may become complicated. Ionospheric sounding experience indicates that while topside sounding echo strengths can actually be enhanced by the presence of irregularities, ground-based sounding indicates that coherent detection techniques can still be employed. In this paper we investigate the plasma conditions that will allow coherent signals to be detected by the RPI and the signatures to be expected, such as scattering and plasma resonances, in the presence of multi-scale irregularities, may possibly have on RPI signals. Sounding of irregular plasma structures in the plasmasphere, plasmapause and magnetopause are also discussed.

Fung, Shing F.↗

Investigation of the Attenuation of Plane Shock Waves Moving over very Rough Surfaces

Experimental measurements of the attenuation of plane shock waves moving over rough walls have been made in a shock tube. Measurements of the boundary-layer characteristics, including thickness and velocity distribution behind the shock, have also been made with the aid of new cal techniques which provide direct information on the local boundary-layer conditions at the rough walls. Measurements of shock speed and shock pressure ratio are presented for both smooth-wall and rough-wall flow over lengths of machined-smooth and rough strips which lined all four walls of the shock tube. A simplified theory based on Von Karman's expression for skin-friction coefficient for flow over rough walls, along with a wave-model concept and extensions to include time effects, is presented. In this theory, the shock-tube flow is assumed to be one-dimensional at all times and the wave-model concept is used to relate the local layer growth to decreases in shock strength. This concept assumes that local boundary-layer growths act as local mass-flow sinks, which give rise to expansion waves which, in turn, overtake the shock and lower its mass flow accordingly.

Huber, Paul W.↗

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗

Scattering Amplitudes in a Neighborhood of Limit Rays in Short-Wave Diffraction by Elongated Bodies of Revolution

In the paper, diffraction problems of a plane wave by smooth, convex, and elongated bodies of revolution are considered within the framework of short-wave approximation (axially symmetric cases). The scattering amplitudes are calculated in the direction of limit rays, and the influence of the elongation of the scatterers on the amplitudes behavior is investigated. Mathematical techniques of our approach are based on the Green’s formulas in the exterior of the scatterers and numerical calculations of the wave field current in the boundary layers in the vicinity of the light-shadow zone. It is established that the elongation of axially symmetric bodies relatively weakly affects the scattering amplitudes of the short-wave asymptotics. The main contribution to the amplitudes is made by the solution of the 2D diffraction problem by a convex, smooth curve in the cross section of the scatterers by a plane containing the rotation axis.

97 MATHEMATICS AND COMPUTING↗

Investigation by Schlieren Technique of Methods of Fixing Fully Turbulent Flow on Models at Supersonic Speeds

An investigation has been made of the effectiveness of two-dimensional wire and three-dimensional grandular-type roughnesses for fixing transition on a 27 deg cone at a free-stream Mach number of 2.20. One height of each type of roughness was investigated over a free-stream Reynolds number per foot range of about 1.0 x 10 (exp 6) to 7.0 x 10 (exp 6). The tests were made at an angle of attack of 0 deg and with essentially zero heat transfer. The results indicate that, within the Reynolds number per foot range investigated, both two-dimensional wire and three-dimensional granular types of roughness caused large forward movements of transition at constant Reynolds number per foot; however, the forward movement with increase in Reynolds number per foot was very slow for the two-dimensional wire roughness and transition never reached the roughness location as it did with the three-dimensional granular roughness. An extrapolated value of the roughness Reynolds number (based on roughness height and flow conditions at the outer edge of the roughness) for the first appearance of forward movement from the smooth-cone case for the three-dimensional granular roughness was in good agreement with the value of about 600 for the first appearance of turbulent bursts obtained on a similar cone by the hot-wire technique. The range of values of roughness Reynolds number for which a fully turbulent boundary layer existed in the vicinity of the roughness element was about 4,500 to 6,500.

Jackson, Mary W.↗

Search for dijet resonances with data scouting in proton-proton collisions at $\sqrt{s}=13$ TeV

A search is presented for narrow resonances, with a mass between 0.6 and 1.8 TeV, decaying to pairs of jets, in proton-proton collisions at $\sqrt{s}=13$ TeV. The search is performed using dijets that are reconstructed, selected, and recorded in a compact form by the high-level trigger in a technique referred to as “data scouting”, from data collected in 2016–2018 corresponding to an integrated luminosity of 117 fb −1 . The dijet mass spectra are well described by a smooth parameterization, and no significant evidence for the production of new particles is observed. Model-independent upper limits are presented on the product of the cross section, branching fraction, and acceptance for the individual cases of narrow quark-quark, quark-gluon, and gluon-gluon resonances, and are compared to the predictions from a variety of models of narrow dijet resonance production. The upper limit on the coupling of a dark matter mediator to quarks is presented as a function of the mediator mass. The sensitivity of this search goes beyond what is expected from statistical scaling with the integrated luminosity alone, as a consequence of the use of fewer parameters in the background function within a more robust statistical procedure.

beyond Standard Model↗