Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Generalised Primal-Dual Grids for Unstructured Co-Volume Schemes

The generation of high-quality staggered unstructured grids is considered, leading to the development of a new optimisation-based strategy designed to construct weighted 'Regular-Power' tessellations appropriate for co-volume type numerical discretisation techniques. This new framework aims to extend the conventional Delaunay-Voronoi primal-dual structure; seeking to assemble generalised orthogonal tessellations with enhanced geometric quality. The construction of these grids is motivated by the desire to improve the performance and accuracy of numerical methods based on unstructured co-volume type schemes, including various staggered grid techniques for the simulation of fluid dynamics and hyperbolic transport. In this study, a new hybrid optimisation strategy is proposed; seeking to optimise the geometry, topology and weights associated with general, two-dimensional Regular-Power tessellations using a combination of gradient-ascent and energy-based techniques. The performance of this new method is tested experimentally, with a range of complex, multi-resolution primal-dual grids generated for various coastal and regional ocean modelling applications.

power diagrams↗
DOE OSTI2024

Evaluation of Radiography for TRISO Buffer Layer Density Measurement

Tristructural isotropic (TRISO) fuel particles consist of a central uranium-bearing kernel and a series of coating layers designed to retain fission products and to ensure fuel performance. Several parameters such as thickness and density must be measured for these coating layers to show that they conform with fuel specifications. Current methods for measuring the density of pyrolytic carbon and silicon carbide layers (liquid gradient density column) and the buffer layer (mercury porosimetry) generate Resource Conservation and Recovery Act (RCRA) radiological-mixed waste. In addition, measurement of buffer and inner pyrolytic carbon layer densities require hot sampling or interrupted coating runs and the mercury porosimetry method used for buffer density measurement only measures the mean buffer density, not the interparticle distribution. A new approach has been evaluated to measure the density of coating layers in TRISO particles based on the dependence of x-ray attenuation in radiographs on material density. This method does not generate RCRA mixed waste, measures density on a particle-by-particle basis, and in principle is capable of measuring the density of all coating layers in a single process. Initial results using thinned TRISO particle sections to evaluate radiography measurement of density as a quality control characterization method are reported herein. In this work, the primary focus is on measurement of the density of the buffer layer; however, with appropriate calibration the method should be applicable to other coating layers. Improvements to the initial method and a full demonstration of the method on the remaining coating layers may be pursued as a future effort.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1987

A finite-difference outer layer and integral inner layer method for the solution of the turbulent boundary layer equations

A new turbulent boundary-layer method is developed which models the inner region with the law of the wall while the outer region uses Clauser's eddy viscosity in Matsuno's finite-difference method. The match point between the inner and outer regions as well as the wall shear stress are determined at each marching step during the computation. Results obtained for incompressible, two-dimensional flow over flat plates and ellipses are compared with solutions from a baseline method which uses a finite-difference method for the entire boundary layer. Since the present method used the finite-difference method in the outer region only, the number of grid points required was about half that needed for the baseline method. Accurate displacement and momentum thicknesses were predicted for all cases. Skin friction was predicted well for the flat plate, but the accuracy decreased significantly for the ellipses. Adding a wake functions to the law of the wall allows some of the pressure gradient effect to be taken into account thereby increasing the accuracy of the method.

Barnwell, R. W.↗
DOE OSTI2022

Minimum feature size control in level set topology optimization via density fields

A level set topology optimization approach that uses an auxiliary density field to nucleate holes during the optimization process and achieves minimum feature size control in optimized designs is explored. The level set field determines the solid-void interface and the density field describes the distribution of a fictitious porous material using the solid isotropic material with penalization. These fields are governed by two sets of independent optimization variables which are initially coupled using a penalty for hole nucleation. The strength of the density field penalization and projection is gradually increased during the optimization process to promote a 0-1 density distribution. In addition, a second penalty regulates the evolution of the density field in the void phase. The treatment of the density field combined with the second penalty mitigate the appearance of small design features. The minimum feature size of optimized designs is controlled by the radius of the linear filter applied to the density optimization variables. The structural response is predicted by the extended finite element method, the sensitivities by the adjoint method, and the optimization variables are updated by a gradient-based optimization algorithm. Numerical examples investigate the robustness of this approach with respect to algorithmic parameters and mesh refinement. The results show the applicability of the combined density level set topology optimization approach for both optimal hole nucleation and for minimum feature size control in 2D and 3D. This comes, however, at the cost of a more complex problem formulation and additional computational cost due to an increased number of optimization variables.

42 ENGINEERING↗
DOE OSTI2024

Topological Optimization with Big Steps

Using persistent homology to guide optimization has emerged as a novel application of topological data analysis. Existing methods treat persistence calculation as a black box and backpropagate gradients only onto the simplices involved in particular pairs. We show how the cycles and chains used in the persistence calculation can be used to prescribe gradients to larger subsets of the domain. In particular, we show that in a special case, which serves as a building block for general losses, the problem can be solved exactly in linear time. This relies on another contribution of this paper, which eliminates the need to examine a factorial number of permutations of simplices with the same value. Here, we present empirical experiments that show the practical benefits of our algorithm: the number of steps required for the optimization is reduced by an order of magnitude.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2023

Generalized moving least squares vs. radial basis function finite difference methods for approximating surface derivatives

Approximating differential operators defined on two-dimensional surfaces is an important problem that arises in many areas of science and engineering. Over the past ten years, localized meshfree methods based on generalized moving least squares (GMLS) and radial basis function finite differences (RBF-FD) have been shown to be effective for this task as they can give high orders of accuracy at low computational cost, and they can be applied to surfaces defined only by point clouds. However, there have yet to be any studies that perform a direct comparison of these methods for approximating surface differential operators (SDOs). The first purpose of this work is to fill that gap. For this comparison, we focus on an RBF-FD method based on polyharmonic spline kernels and polynomials (PHS+Poly) since they are most closely related to the GMLS method. Additionally, we use a relatively new technique for approximating SDOs with RBF-FD called the tangent plane method since it is simpler than previous techniques and natural to use with PHS+Poly RBF-FD. Further, the second purpose of this work is to relate the tangent plane formulation of SDOs to the local coordinate formulation used in GMLS and to show that they are equivalent when the tangent space to the surface is known exactly. The final purpose is to use ideas from the GMLS SDO formulation to derive a new RBF-FD method for approximating the tangent space for a point cloud surface when it is unknown. For the numerical comparisons of the methods, we examine their convergence rates for approximating the surface gradient, divergence, and Laplacian as the point clouds are refined for various parameter choices. We also compare their efficiency in terms of accuracy per computational cost, both when including and excluding setup costs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Simulation of multi-shell fullerenes using Machine-Learning Gaussian Approximation Potential

Multi-shell fullerenes ”buckyonions ” were simulated, starting from initially random configurations, using a density-functional-theory (DFT)-trained machine-learning carbon potential within the Gaussian Approximation Potential (GAP) Framework [Volker L. Deringer and Gábor Csányi, Phys. Rev. B 95, 094203 (2017)]. Fullerenes formed from seven different system sizes, ranging from 60 ~ 3774 atoms, were considered. The buckyonions are formed by clustering and layering starting from the outermost shell and proceeding inward. Inter-shell cohesion is partly due to interaction between delocalized π electrons protruding into the gallery. The energies of the models were validated ex post facto using density functional codes, VASP and SIESTA , revealing an energy difference within the range of 0.02 - 0.08 eV/atom after conjugate gradient energy convergence of the models was achieved with both methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Dynamics of Glyceline and Interactions of Constituents: A Multitechnique NMR Study

In this work, the dynamics of the organic components of the deep eutectic solvent (DES) glyceline are analyzed using an array of complementary nuclear magnetic resonance (NMR) methods. Fast-Field Cycling 1 H relaxometry, Pulsed Field Gradient diffusion, Nuclear Overhauser Effect Spectroscopy (NOESY), 13 C NMR relaxation and pressure dependent NMR experiments are deployed to sample a range of frequencies and modes of motion of the glycerol and choline components of the DES. Generally, translational and rotational diffusion of glycerol are more rapid than those of choline while short range rotational motions observed from 13 C relaxation indicate slow local motion of glycerol at low choline chloride (ChCl) content. The rates of glycerol and choline local motions become more similar at higher ChCl. This result taken together with pressure dependent NMR studies show that the addition of ChCl makes it easier to disrupt glycerol packing. Finally, a relatively slow hydroxyl H-exchange process between glycerol and choline protons is deduced from the data. Consistent with this, NOESY results indicate relatively little direct H-bonding between glycerol and choline. These results suggest that the glycerol H-bonding network is disrupted as choline is added, but primarily in regions where there is intimate mixing of the two components. Thus, the local dynamics of most of the glycerol, resembles that of pure glycerol until substantial choline chloride is present.

13C NMR relaxation↗
DOE OSTI2024

Spin Decoherence in VOPc@graphene Nanoribbon Complexes

Carbon nanoribbons or nanographene qubit arrays can facilitate quantum-to-quantum transduction between light, charge, and spin, making them an excellent testbed for fundamental science in quantum coherent systems and for the construction of higher-level qubit circuits. Here, in this work, we study spin decoherence due to coupling with a surrounding nuclear spin bath of an electronic molecular spin of a vanadyl phthalocyanine (VOPc) molecule integrated onto an armchair-edged graphene nanoribbon (GNR). Density functional theory (DFT) is used to obtain ground-state atomic configurations. Decay of spin coherence in Hahn echo experiments is then simulated using the cluster correlation expansion method with a spin Hamiltonian involving hyperfine and electric field gradient tensors calculated from DFT. We find that the decoherence time T 2 is anisotropic with respect to magnetic field orientation and determined only by the hydrogen nuclear spins on both VOPc and GNR. Large electron spin echo envelope modulation (ESEEM) due to nitrogen and vanadium nuclear spins is present at specific field ranges and can be completely suppressed by tuning the magnetic field. The relation between these field ranges and the hyperfine interactions is analyzed. The effects of interactions with the nuclear quadrupole moments are also studied, validating the applicability and limitations of the spin Hamiltonian when they are disregarded.

Hamiltonians↗
DOE OSTI2026

Explainable machine learning for incipient anomaly detection in compact molten salt heat exchanger with overlapping feature distributions

High-temperature molten salt-cooled reactors (MSCRs) are a promising next-generation nuclear technology option, offering efficient power conversion and inherent safety features. However, the reliability of these systems depends on the robust operation of heat exchangers (HXs), which are susceptible to failure due to temperature gradients and channel plugging caused by fluid freezing. Conventional monitoring methods, relying on inlet and outlet measurements, lack the spatial resolution needed to detect early-stage faults. We propose a novel design of a compact salt-to-salt matrix-type HX design consisting of interleaved arrays of parallel tubes, with integrated synthetic fiber optic distributed temperature sensing (DTS) to enable localized detection of incipient faults. To evaluate performance of this design, we generate high-fidelity synthetic data using heat transfer computational modeling to simulate channel plugging, and introduce sensor noise for realistic modeling of measurements. The dataset comprises of 97% normal operation and 3% anomaly cases, with each anomaly class representing 1% of the data. These early anomalies result in overlapping temperature profiles between normal and faulty channels, producing a non-separable dataset that challenges traditional classification techniques. We benchmark eight supervised machine learning (ML) models and demonstrate that XGBoost achieves the highest performance. To improve transparency, we develop an explainability framework combining Shapley values and partially ordered sets (POSETs) to quantify and structurally analyze feature importance. This approach identifies both dominant predictors and ambiguous feature relationships, enhancing trust and interpretability. Our results highlight the potential of combining DTS and explainable ML with intelligent feature selection to improve predictive maintenance and ensure operational resilience in advanced nuclear systems.

Prantikos, Konstantinos [Argonne National Laborato↗
DOE OSTI2021

Lagrangian analysis for turbulent transport in variable-density turbulence

In this report Lagrangian analysis of materially conserved scalars is applied to the problem of turbulent transport in variable-density flows. The consequences of an additional material conserved quantity, the density, is generally not acknowledged and leads to significant and meaningfully different expressions for turbulent transport in the moment equations. The formal Lagrangian analysis produces gradient transport expressions substantially different from those obtained by the physically intuitive “argument by analogy” method used in computational models. Various intuitive arguments, in Favre and Reynolds averaged settings, are contrasted to the formal Lagrangian results. Using expressions from the formal analysis, we derive consistent gradient transport closures for the turbulent transport terms that appear in the first- and second-order Favre moment equations. Results for coupled multispecies turbulent transport are given. The analysis is limited to variable-density turbulence in which the dilatation of the fluctuating velocity is small. The results are applicable to turbulent combustion and to stellar convection problems in which the density fluctuations are on the order of the mean density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

First Measurement of Z Opacity Sample Evolution near Solar Interior Conditions Using Time-Resolved Spectroscopy

Opacity model differences with Fe opacity measurements at high temperature (T>1⁢8⁢0 eV ) and high electron density (𝑛 𝑒 >3×10 22 cm −3 ) at nearly solar interior conditions have remained unresolved [Bailey et al., Nature 517, 56 (2015) and Nagayama et al., Phys. Rev. Lett. 122, 235001 (2019)]. Systematic errors from temporal gradients are one hypothesis for the discrepancy. Past data recorded on x-ray film provided spectral measurements over a time determined by the 2.8-ns backlighter duration. Here, we present the first measurements of opacity sample temporal evolution using novel hCMOS ultrafast x-ray camera technology. The measured conditions, measured backlighter time history, and modeled opacities are used to show that temporal gradients do not resolve the model-data discrepancy. The methods demonstrated provide potential advantages, including opacities at more extreme conditions, spectral line shift measurements, and improved measurements at other facilities.

atomic spectra↗
DOE OSTI2023

Topological Guided Detection of Extreme Wind Phenomena: Implications for Wind Energy

Extreme wind phenomena play a crucial role in the efficient operation of wind farms for renewable energy generation. However, existing detection methods are computationally expensive, limited to specific coordinate. In real-world scenarios, understanding the occurrence of these phenomena over a large area is essential. Therefore, there is a significant demand for a fast and accurate approach to forecast such events. In this paper, we propose a novel method for detecting wind phenomena using topological analysis, leveraging the gradient of wind speed or critical points in a topological framework. By extracting topological features from the wind speed profile within a defined region, we employ topological distance to identify extreme wind phenomena. Our results demonstrate the effectiveness of utilizing topological features derived from regional wind speed profiles. We validate our approach using high-resolution simulations with the Weather Research and Forecasting model (WRF) over a month in the US East Coast.

MATHEMATICS AND COMPUTING,WIND ENERGY↗
DOE OSTI2023

Physics-Informed Deep Learning-Based Proof-of-Concept Study of a Novel Elastohydrodynamic Seal for Supercritical CO 2 Turbomachinery

Supercritical carbon dioxide (sCO 2 ) power cycles show promising potential of higher plant efficiencies and power densities for a wide range of power generation applications such as fossil fuel power plants, nuclear power production, solar power, and geothermal power generation. sCO 2 leakage through the turbomachinery has been one of the main concerns in such applications. To offer a potential solution, we propose an elastohydrodynamic (EHD) seal that can work at elevated pressures and temperatures with low leakage and minimal wear. The EHD seal has a very simple, sleeve-like structure, wrapping on the rotor with minimal initial clearance at micron levels. Here, in this work, a proof-of-concept study for the proposed EHD seal was presented by using the simplified Reynolds equation and Lame’s formula for the fluid flow in the clearance and for seal deformation, respectively. The set of nonlinear equations was solved by using both the conventional Prediction–Correction (PC) method and modern Physics-Informed Neural Network (PINN). It was shown that the physics-informed deep learning method provided good computational efficiency in resolving the steep pressure gradient in the clearance with good accuracy. The results showed that the leakage rates increased quadratically with working pressures and reached a steady-state at high-pressure values of 15~20 MPa, where Q = 300 g/s at 20 MPa for an initial seal clearance of 255 μm. This indicates that the EHD seal could be tailored to become a potential solution to minimize the sCO 2 discharge in power plants.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1972

Random response of rectangular panels to the pressure field beneath a turbulent boundary layer in subsonic flows

The response of a rectangular panel under the excitation of a turbulent boundary layer with a zero longitudinal mean pressure gradient in a subsonic flow was studied in detail. The method of normal mode was used together with the technique of spectral analysis. Both simply supported and clamped edge conditions of a panel were considered, and the displacement power spectral density of the panel response was computed. The results for the clamped edge panel compare favorably with existing data. Charts of structural acceptance, which provide a framework for estimating the response of other rectangular panels, are presented, and the physical significances of structural acceptances discussed.

Chyu, W. J.↗
NASA NTRS1983

Moire deflectometry with deferred analysis

A simplified version of moire deflectometry for flow visualization is presented. The modification comprises replacement of the second grating, the ground glass screen, and the camera with a plate film holder to obtain a shadow picture. Postanalysis proceeds by placing a grid transparency of the same periodicity of the remaining glass grating on the negative for projection viewing. The rotational angle between the grating and the plate alters the moire fringe periodicity. Use of the method in studying flows is predicted to reveal weak density gradient areas and shocks.

Rhodes, D. B.↗
NASA NTRS1985

A space-marching method for incompressible Navier-Stokes equations

This paper deals with the development of a space-marching method for incompressible flows. The method solves the continuity and momentum equations as a coupled system at each streamwise station. The character of the system of equations has been changed from elliptic to hyperbolic/parabolic in order to enable the equations to be marched in space. The present method has many advantages compared to the existing parabolic or space-marching methods for incompressible flow: (1) it avoids the solution of Poisson equations, (2) it conserves the mass flow with no additional computation, (3) it does not require the specification of an assumed pressure field when used in the prediction of duct flows. The present method can capture strong secondary velocities and strong transverse pressure gradients. Predictions of the flow through straight and curved ducts are in good agreement with analytical and experimental results.

Pouagare, M.↗
NASA NTRS1985

An improved viscid/inviscid interaction procedure for transonic flow over airfoils

A new interacting boundary layer approach for computing the viscous transonic flow over airfoils is described. The theory includes a complete treatment of viscous interaction effects induced by the wake and accounts for normal pressure gradient effects across the boundary layer near trailing edges. The method is based on systematic expansions of the full Reynolds equation of turbulent flow in the limit of Reynolds numbers, Reynolds infinity. Procedures are developed for incorporating the local trailing edge solution into the numerical solution of the coupled full potential and integral boundary layer equations. Although the theory is strictly applicable to airfoils with cusped or nearly cusped trailing edges and to turbulent boundary layers that remain fully attached to the airfoil surface, the method was successfully applied to more general airfoils and to flows with small separation zones. Comparisons of theoretical solutions with wind tunnel data indicate the present method can accurately predict the section characteristics of airfoils including the absolute levels of drag.

Melnik, R. E.↗