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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Vector Field Decompositions Using Multiscale Poisson Kernel

Extraction of multiscale features using scale-space is one of the fundamental approaches to analyze scalar fields. However, similar techniques for vector fields are much less common, even though it is well known that, for example, turbulent flows contain cascades of nested vortices at different scales. The challenge is that the ideas related to scale-space are based upon iteratively smoothing the data to extract features at progressively larger scale, making it difficult to extract overlapping features. Instead, we consider spatial regions of influence in vector fields as scale, and introduce a new approach for the multiscale analysis of vector fields. Rather than smoothing the flow, we use the natural Helmholtz-Hodge decomposition to split it into small-scale and large-scale components using progressively larger neighborhoods. Our approach creates a natural separation of features by extracting local flow behavior, for example, a small vortex, from large-scale effects, for example, a background flow. We demonstrate our technique on large-scale, turbulent flows, and show multiscale features that cannot be extracted using state-of-the-art techniques.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Investigation of poly(phenylacetylene) derivatives for carbon precursor with high carbon yield and good solubility

This paper investigates a family of poly(phenylacetylene) derivatives with a p-electrons conjugated polymer backbone and side groups containing only sp 2 and sp carbons. The objective is to identify the suitable carbon precursor that is processible and can be transformed to carbonaceous material with high carbon yield by a simple (one-step) thermal transformation process (without any external reagent). Poly(phenylacetylene) with para-substituted acetylene group poly(PA-A) shows an exceptionally high C-yield (~90%) in one-step heating from ambient temperature to 1000 °C under N 2 atmosphere. Unfortunately, this poly(PA-A) polymer is quite sensitive to heat and light with very limited solubility. On the other hand, poly(phenylacetylene) with para-substituted phenylacetylene group, i.e. poly(PA-PA), offers a relatively high C-yield (~80%) and also good solubility in common organic solvents, such as toluene and tetrahydrofuran (THF). Several uniform dark-red poly(PA-PA) fibers with smooth surface and fiber diameter in the range of 3–6 μm were prepared from 30 wt% poly(PA-PA)/THF solution using electrospinning technique. Furthermore, the resulting precursor fibers were converted to the corresponding carbon structure in a one-step thermal heating process under N 2 atmosphere. Both X-ray and Raman spectra show the polymorphous carbon morphology with the graphene crystalline domains.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Mass measurements of neutron-rich nuclei near $N = 70$

The astrophysical origin for the chemical elements between the first and second r-process peaks is a matter of intense debate, with a number of nucleosynthesis processes at explosive stellar environments possibly contributing to their production. Reliable data on the trends of neutron separation energies of neutron-rich isotopes are required to model neutron-capture processes that would produce these elements. Masses of 104 Y, 106 Zr, 112 Mo, and 115 Tc have been measured with the time-of-flight-magnetic-rigidity (ToF–Bρ) technique at the National Superconducting Cyclotron Laboratory at Michigan State University. The experiment is the first application of the ToF–Bρ technique at the S800 spectrograph that reached the mass region relevant to heavy-element nucleosynthesis. Finally, the two-neutron separation energy deduced from the measured masses exhibits a smooth trend consistent with the theoretical predictions within the range of experimental uncertainty, indicating that there is no sudden shape transition in these isotopes as hinted at by previous data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Combined Dimensional and Topology Optimization of Synchronous Machine Rotors Using a Material Density Interpolation Method

This article presents a magneto-structural combined dimensional and topology optimization technique for interior permanent magnet synchronous machine (IPMSM) and wound field synchronous machine (WFSM) rotors. Dimensional changes to the permanent magnet (PM) and rotor winding location or size are accomplished by interpolating or projecting a smoothed Heaviside rectangular function representing the presence of PM material onto the IPMSM rotor design domain mesh, whereas the presence of copper material for WFSM rotor is interpolated in the same manner. A density based Solid Isotropic Material with Penalization (SIMP) topology optimization approach is then used to vary the presence of electrical steel in mesh elements to form flux barriers around the PM. The proposed method enforces a defined shape for the PM without requiring the mesh in the design domain to be deformed. Four examples are presented to demonstrate the technique: Two flat bar IPMSM, one V-shaped IPMSM, and a WFSM. Here, a comparative study is performed on one of the flat bar IPMSMs using Metamodel of Optimal Prognosis (MOP) based method.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2023

Temperature uncertainty modelling with proxy structural data as geostatistical constraints for well siting: an example applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural geology information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aims of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations. The data used to support the modelling are measured temperatures and structural proxies from 3D geologic modelling (distance to fault, distance to fault intersections and terminations, Coulomb stress change and dilation tendency), the latter considered ‘secondary’ data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally-likely temperature models (i.e. realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: V prior and entropy. V prior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because of the dipping ellipsoid shape created and that the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e. linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g. dilation v. distance to fault termination) should be used given the specific geothermal system. In general, V prior from local varying mean results identify locations that are close to high values for the structural proxies: areas with higher probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comparison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

15 GEOTHERMAL ENERGY↗

DOE OSTI2019

Suppression of hybridization by Cd doping in CeCoIn 5

The interaction between superconductivity and magnetism in heavy fermion compounds has remained a challenge ever since the discovery of unconventional superconductivity in such systems. The heavy fermion Ce T In 5 ( T = Co,Rh,Ir) family offers a good platform to study the relationship between superconductivity and magnetism, whose ground state can be easily tuned by chemical doping, making the ambient-pressure experimental techniques like angle-resolved photoemission spectroscopy possible to address this old but engaging question. By Cd doping, the superconducting state of CeCoIn 5 can be smoothly tuned to an antiferromagnetic ground state, which can be reversed by applying pressure. In this study, we present an electronic structure study of CeCo ( In 0.85 Cd 0.15 ) 5 with an antiferromagnetic ground state and the results are compared with CeCoIn 5 . We found that the hybridization strength between the f electrons and conduction electrons in CeCoIn 5 has been suppressed by Cd doping, but still the Fermi surface includes a significant contribution from the Ce f electrons. Our results find that the electronic tuning effect by Cd doping arises from the change in the hybridization strength between the Ce f electrons and conduction electrons caused by the Cd dopant, and the subtle interplay between itinerancy and localization of the f electrons determines the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Directional Crystallization of Conjugated Molecules during Coating Processes

The coating of organic molecules from the solution phase can result in directional crystal growth under certain conditions, even on a smooth isotropic surface and without the need of any kind of graphoexpitaxial preparation of the substrate. Based on reviewing the results from a variety of coating techniques and coating parameters, we identified that it is crucial for the coating speed to match the growth speed of the fastest growing crystal plane to achieve a high degree of directional crystallization.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Li-ion Battery Material phase prediction through Hierarchical Curriculum Learning

Li-ion Batteries (LIB), one of the most efficient energy storage devices, are widely adopted in many industrial applications. Imaging data of these battery electrodes obtained from X-ray tomography can explain the distribution of material constituents and allow reconstructions to study electron transport pathways. Therefore, it can eventually help quantify various associated properties of electrodes (e.g., volume-specific surface area, porosity) which determine the performance of batteries. However, these images often suffer from low image contrast between multiple material constituents , making it difficult for humans to distinguish and characterize these constituents through visualization. A minor error in detecting distributions among the material constituents can lead to a high error in the calculated parameters of material properties.We present a novel hierarchical curriculum learning framework to address the complex task of estimating material constituent distribution in battery electrodes. To provide spatially smooth prediction, our framework comprises three modules: (i) an uncertainty-aware model trained to yield inferences conditioned upon global knowledge of material distribution, (ii) a technique to capture relatively more fine-grained (local) distributional signals, (iii) an aggregator to appropriately fuse the local and global effects towards obtaining the final distribution.

Tabassum, Anika↗

DOE OSTI2024

Dynamically Learning Incentives for Load Control

As electrical generation becomes more distributed and volatile, and loads become more uncertain, controllability of distributed energy resources (DERs), regardless of their ownership status, will be necessary for grid reliability. Grid operators lack direct control over end-users' grid interactions, such as energy usage, but incentives can influence behavior -- for example, an end-user that receives a grid-driven incentive may adjust their consumption or expose relevant control variables in response. A key challenge in studying such incentives is the lack of data about human behavior, which usually motivates strong assumptions, such as distributional assumptions on compliance or rational utility-maximization. In this paper, we propose a general incentive mechanism in the form of a constrained optimization problem -- our approach is distinguished from prior work by modeling human behavior (e.g., reactions to an incentive) as an arbitrary unknown function. We propose feedback-based optimization algorithms to solve this problem that each leverage different amounts of information and/or measurements. We show that each converges to an asymptotically stable incentive with (near)-optimality guarantees given mild assumptions on the problem. Finally, we evaluate our proposed techniques in voltage regulation simulations on standard test beds. We test a variety of settings, including those that break assumptions required for theoretical convergence (e.g., convexity, smoothness) to capture realistic settings. In this evaluation, our proposed algorithms are able to find near-optimal incentives even when the reaction to an incentive is modeled by a theoretically difficult (yet realistic) function.

demand response↗

DOE OSTI2024

Smooth trends in fermium charge radii and the impact of shell effects

The quantum-mechanical nuclear-shell structure determines the stability and limits of the existence of the heaviest nuclides with large proton numbers Z ≳ 100. Shell effects also affect the sizes and shapes of atomic nuclei, as shown by laser spectroscopy studies in lighter nuclides. However, experimental information on the charge radii and the nuclear moments of the heavy actinide elements, which link the heaviest naturally abundant nuclides with artificially produced superheavy elements, is sparse. Here we present laser spectroscopy measurements along the fermium (Z = 100) isotopic chain and an extension of data in the nobelium isotopic chain (Z = 102) across a key region. Multiple production schemes and different advanced techniques were applied to determine the isotope shifts in atomic transitions, from which changes in the nuclear mean-square charge radii were extracted. A range of nuclear models based on energy density functionals reproduce well the observed smooth evolution of the nuclear size. Both the remarkable consistency of model prediction and the similarity of predictions for different isotopes suggest a transition to a regime in which shell effects have a diminished effect on the size compared with lighter nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Spatially quasi-periodic bifurcations from periodic traveling water waves and a method for detecting bifurcations using signed singular values

We present a method of detecting bifurcations by locating zeros of a signed version of the smallest singular value of the Jacobian. This enables the use of quadratically convergent root-bracketing techniques or Chebyshev interpolation to locate bifurcation points. Only positive singular values have to be computed, though the method relies on the existence of an analytic or smooth singular value decomposition (SVD). The sign of the determinant of the Jacobian, computed as part of the bidiagonal reduction in the SVD algorithm, eliminates slope discontinuities at the zeros of the smallest singular value. We use the method to search for spatially quasi-periodic traveling water waves that bifurcate from large-amplitude periodic waves. The water wave equations are formulated in a conformal mapping framework to facilitate the computation of the quasi-periodic Dirichlet-Neumann operator. We find examples of pure gravity waves with zero surface tension and overhanging gravity-capillary waves. In both cases, the waves have two spatial quasi-periods whose ratio is irrational. We follow the secondary branches via numerical continuation beyond the realm of linearization about solutions on the primary branch to obtain traveling water waves that extend over the real line with no two crests or troughs of exactly the same shape. The pure gravity wave problem is of relevance to ocean waves, where capillary effects can be neglected. Such waves can only exist through secondary bifurcation as they do not persist to zero amplitude. The gravity-capillary wave problem demonstrates the effectiveness of using the signed smallest singular value as a test function for multi-parameter bifurcation problems. This test function becomes mesh independent once the mesh is fine enough.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Data-Driven Reliability Assessment for Marine Renewable Energy Enabled Island Power Systems

Marine renewable energy (MRE) resources are highly predictable and persistent sources of energy, when compared to other renewable sources like wind and solar. These lend them favorably for potential grid applications, particularly for coastal/island power systems where their generation potential is high. Island power systems, on the other hand, are either supported by onsite generation or by transported energy from the mainland grid. Therefore, robustness of grid operations depend heavily on the diversity of onsite generation resources and the reliability of the power transportation medium. Issues relating to either of these two factors may lead to impediments in smooth and reliable operation of the power system. Analyzing and quantifying operational risks for such island power systems with diverse non-conventional generation portfolios through conventional techniques can also prove to be cumbersome, often requiring multiple different inputs. Therefore, in this paper, we firstly present a novel, purely data-driven formulation which quantifies the operational reliability of such island power systems through minimal input data. Specifically, our proposed methodology only relies on historical knowledge of typical hourly load and generation profiles to quantify associated operational risks. Subsequently, we use our proposed formulation to evaluate the effectiveness of MRE resources (over other renewable resources like wind and solar) in providing resilience benefits to island power systems. The proposed formulation is demonstrated with a case study for an island power system in Nantucket, MA.

Chalishazar, Vishvas H.↗

DOE OSTI2024

Dynamic Modeling and Analysis for Large-Scale Renewable Energy Integration

Lina He will present dynamic modeling techniques for the integration of large-scale renewable energy systems. The presentation will address the challenges associated with high renewable penetration and provide frameworks for ensuring smooth and reliable grid operations.

Dynamic modeling, large-scale renewable integratio↗

DOE OSTI2022

Molecular Reaction Imaging of a Surface Recombination Process Explains Performance Variation Among Smooth MoS 2 Photoelectrodes

Transition metal dichalcogenides (TMD) such as WSe 2 and MoS 2 are highly efficient and stable light absorbers in TMD∣I − ,I 3 − ∣Pt liquid junction solar cells. It is generally accepted that TMD crystals with a large fraction of exposed edge sites exhibit lower power conversion efficiencies (PCEs) than apparently smooth crystals. However, one open question is why does the PCE vary significantly from one crystal to another? Answering this critical question could lead to robust syntheses for high quality and uniform TMD samples. In this work, we apply nanoscale photoelectrochemical microscopy techniques to study n-type TMD nanoflake∣I − ,I 3 − ∣Pt cells. Using a combination of near-diffraction-limited photocurrent mapping and molecular reaction imaging techniques, we reveal a previously hidden surface recombination process: photogenerated holes in hidden p-type domains travel micron-scale distances parallel to the solid/liquid interface and preferentially react with iodide at step-edges. The overall efficiency of the nanoflake, as evidenced from whole nanoflake-level photoelectrochemical measurements, is dictated by the size, efficiency, and location of n- and p-type domains. Furthermore, these results provide a unifying view of efficiency losses in smooth TMD photoelectrodes and open the possibility to design electrode architectures that leverage the long-range lateral charge transport property for photoelectrocatalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Low-Cost Heliostat for High-Flux Small-Area Receivers (Final Technical Report)

This project analyzed a two-stage heliostat concept consisting of a tracking stage and a concentrating stage. The tracking stage uses mirrors mounted on a common drive that move to track the sun. The concentrating stage consists of stationary mirrors that each have a unique angle to direct rays towards a small-area, high-flux, point-focused receiver. By splitting the collection and concentrating process into two stages, multiple small, inexpensive mirrors can share a structure and be controlled by a single drive in the tracking stage. The project effort developed modeling techniques that were specifically relevant to this two-stage heliostat concept. Both field-level and unit-level models were developed. The field-level model does not explicitly consider unit-level losses which are predicted by the unit-level model and then integrated into the field-level model through a correlation referred to as an efficiency modifier. This approach is referred to as the two-model approach; the development and demonstration of this two-model approach for a multi-stage heliostat technology is a key outcome of this work. The field-level model is used to design a field that hits a specific design day power given a set of heliostat design parameters. An oversized field is simulated and then heliostat units are removed based on their annual energy production in order to generate the highest performing field. The field reduction procedure fits a smooth curve fit to annual energy production as a function of position in the field which has the effect of reducing the noise that is otherwise caused by the Monte Carlo ray tracing technique. This approach is referred to as the annual energy fit method and substantially reduces computational run time for a given field level modeling accuracy. The annual energy fit approach enables the selection of a properly sized, high-performing field using orders of magnitude fewer rays than would otherwise be possible and the development of this approach is a second key outcome of this work. These models are used within a genetic optimization algorithm in order to optimize the geometric parameters associated with a heliostat in order to achieve the lowest cost per unit of collected design day power. The cost modeling that underlies the optimization is a simple, scaling type analysis backed up by a much more detailed Design for Manufacture and Assembly (DFMA) analysis. Although the figure of merit used for optimization was not cost per mirror area, this metric is reasonable to use as a means of comparison. The optimally designed 500 kW design has a tracking mirror specific cost of $181.85/m 2 , which is significantly larger than the target value and also larger than the current state of the art. The cost of the torque-tube type linkages contributed substantially to the overall cost. Based on this observation, potentially attractive alternative design configuration utilizing a capstan type actuation system should be investigated. Finally, NREL compared the performance of the two-stage heliostat to the performance of a focused and different sized flat conventional heliostats and showed that, as expected, additional losses versus the convention heliostat caused by a worse cosine efficiency, two stages of reflection, and interstage interactions. The two-stage heliostat requires around 75% more reflective area than a flat 1x1 meter conventional heliostat (similar to a focused heliostat) and 40% more than a flat 2x2 meter conventional heliostat.

14 SOLAR ENERGY↗

DOE OSTI2022

Learning Optimal Multigrid Smoothers via Neural Networks

Multigrid methods are one of the most efficient techniques for solving large sparse linear systems arising from partial differential equations (PDEs) and graph Laplacians from machine learning applications. One of the key components of multigrid is smoothing, which aims at reducing high-frequency errors on each grid level. However, finding optimal smoothing algorithms is problem-dependent and can impose challenges for many problems. In this paper, we propose an efficient adaptive framework for learning optimized smoothers from operator stencils in the form of convolutional neural networks (CNNs). Here, the CNNs are trained on small-scale problems from a given type of PDEs based on a supervised loss function derived from multigrid convergence theories and can be applied to large-scale problems of the same class of PDEs. Numerical results on anisotropic rotated Laplacian problems and variable coefficient diffusion problems demonstrate improved convergence rates and solution time compared with classical hand-crafted relaxation methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

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↗

DOE OSTI2023

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↗