SEARCH · Engineering Papers

Results for “discretization error”

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 OSTI2021

Physics-Informed Neural Networks for Heat Transfer Problems

Abstract Physics-informed neural networks (PINNs) have gained popularity across different engineering fields due to their effectiveness in solving realistic problems with noisy data and often partially missing physics. In PINNs, automatic differentiation is leveraged to evaluate differential operators without discretization errors, and a multitask learning problem is defined in order to simultaneously fit observed data while respecting the underlying governing laws of physics. Here, we present applications of PINNs to various prototype heat transfer problems, targeting in particular realistic conditions not readily tackled with traditional computational methods. To this end, we first consider forced and mixed convection with unknown thermal boundary conditions on the heated surfaces and aim to obtain the temperature and velocity fields everywhere in the domain, including the boundaries, given some sparse temperature measurements. We also consider the prototype Stefan problem for two-phase flow, aiming to infer the moving interface, the velocity and temperature fields everywhere as well as the different conductivities of a solid and a liquid phase, given a few temperature measurements inside the domain. Finally, we present some realistic industrial applications related to power electronics to highlight the practicality of PINNs as well as the effective use of neural networks in solving general heat transfer problems of industrial complexity. Taken together, the results presented herein demonstrate that PINNs not only can solve ill-posed problems, which are beyond the reach of traditional computational methods, but they can also bridge the gap between computational and experimental heat transfer.

Engineering↗

DOE OSTI2021

Laplacian Smoothing Stochastic Gradient Markov Chain Monte Carlo

As an important Markov chain Monte Carlo (MCMC) method, the stochastic gradient Langevin dynamics (SGLD) algorithm has achieved great success in Bayesian learning and posterior sampling. Furthermore, SGLD typically suffers from a slow convergence rate due to its large variance caused by the stochastic gradient. In order to alleviate these drawbacks, we leverage the recently developed Laplacian smoothing technique and propose a Laplacian smoothing stochastic gradient Langevin dynamics (LS-SGLD) algorithm. We prove that for sampling from both log-concave and non-log-concave densities, LS-SGLD achieves strictly smaller discretization error in 2-Wasserstein distance, although its mixing rate can be slightly slower. Experiments on both synthetic and real datasets verify our theoretical results and demonstrate the superior performance of LS-SGLD on different machine learning tasks including posterior sampling, Bayesian logistic regression, and training Bayesian convolutional neural networks.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

A Fast Solver for the Fractional Helmholtz Equation

The purpose of this paper is to study a Helmholtz problem with a spectral fractional Laplacian, instead of the standard Laplacian. Recently, it has been established that such a fractional Helmholtz problem better captures the underlying behavior in Geophysical Electromagnetics. In this work, we establish the well-posedness and regularity of this problem. We introduce a hybrid spectral-finite element approach to discretize it and show well-posedness of the discrete system. In addition, we derive a priori discretization error estimates. Finally, we introduce an efficient solver that scales aswell as the best possible solver for the classical integer-order Helmholtz equation. We conclude withseveral illustrative examples that confirm our theoretical findings.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

CFDverify

Estimating the discretization error of computational fluid dynamics (CFD) or other scientific codes as part of solution verification is often a non-trivial part of the analysis process. Methods can be complicated and may include assumptions/qualifications that need to be checked during analysis. CFD analysts, therefore, are likely to make errors when trying to conduct this necessary analysis on their own in not knowing about the best method for their problem, not correctly implementing a method, or in not have diagnostic tools to determine if the method was correctly applied.

Weinmeister, Justin [Oak Ridge National Laboratory↗

DOE OSTI2020

Evaluating Implementations of the Immersed Boundary Method in the Weather Research and Forecasting Model

The terrain-following coordinate system used by many atmospheric models can cause numerical instabilities due to discretization errors as resolved terrain slopes increase and the grid becomes highly skewed. The immersed boundary (IB) method, which does not require the grid to conform to the terrain, has been shown to alleviate these errors, and has been used successfully for high-resolution atmospheric simulations over steep terrain, including vertical building surfaces. Since many previous applications of IB methods to atmospheric models have used very fine grid resolution (5 m or less), the present study seeks to evaluate IB method performance over a range of grid resolutions and aspect ratios. Two classes of IB algorithms, velocity reconstruction and shear stress reconstruction, are tested within the common framework of the Weather Research and Forecasting (WRF) Model. Performance is evaluated in two test cases, one with flat terrain and the other with the topography of Askervein Hill, both under neutrally stratified conditions. WRF-IB results are compared to similarity theory, observations, and native WRF results. Despite sensitivity to the location at which the IB intersects the model grid, the velocity reconstruction IB method shows consistent performance when used with a hybrid RANS/LES surface scheme. The shear stress reconstruction IB method is not sensitive to the grid intersection, but is less consistent and near-surface velocity errors can occur at coarse resolutions. This study represents an initial investigation of IB method variability across grid resolutions in WRF. Future work will focus on improving IB method performance at intermediate to coarse resolutions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Nek5000: improvements in the available RANS models, meshing, tutorials, and training

This year, the Nuclear Energy Advanced Modeling Simulation program (NEAMS) thermal-hydraulics report for Nek5000 NRC- and verification and validation (V&V)-driven development focuses on following areas of code application and improvement. First we have continued improvements of RANS modeling capabilities in Nek5000 including improved k-tau model focusing mostly on wallfunction initial implementation with spectral element method (SEM) and initiating investigation of an alternative approach XSEM that greatly reduces discretization errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

CFD Benchmark of Pressure Drop in a 61-Pin Wire-Wrapped Assembly with Blocked Channels Using NekRS

Thermal hydraulic behaviors of wire-wrapped rod bundles continue to be a subject of intense research. One of the leading next-generation designs, the sodium fast reactor, relies on a hexagonal assembly of wire-wrapped fuel pins. The issue of potential channel blockage has been raised as a safety concern due to the tightly packed arrangement of the fuel pins. This has led to several recent experimental and computational studies working to quantify the potential impact on the fluid flow and heat transfer behaviors of such blockages. The objective of the present study is to benchmark the high-fidelity NekRS CFD solver in predicting pressure drop for large blockages against available experimental data. A 61-pin wire-wrapped fuel assembly with two flow blockage configurations has been simulated and investigated at various low to moderate Reynolds numbers. The NekRS solver has been shown to yield exponentially decreasing spatial discretization errors with increasing polynomial order. All simulated results agreed well with measured data, which indicates that the overall methodology is adequate and consistent. The results of this benchmark study demonstrate the accuracy of NekRS for sodium fast reactor hydrodynamic simulations, increasing the confidence in its use for design, licensing, and analysis activities.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Center for Tokamak Transient Simulations (RPI Unstructured Mesh Developments for FES SciDAC4 Partnerships) (Final Report)

The goal of this project was the development of unstructured mesh technologies for fusion simulation codes” for FES SciDAC partnerships and to integrate those technologies into the simulation workflows of those partnerships. Specific developments were executed in support of the following FES SciDAC4 partnerships: Partnership Center for High‐fidelity Boundary Plasma Simulation (HBPS), Center for Integrated Simulation of Fusion Relevant RF (RF‐SciDAC), Center for Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces (PSI2), and Center for Tokamak Transient Simulations (CTTS). The key unstructured mesh development areas addressed in this project include (i) methods to most effectively perform PIC calculations on unstructured meshes; (ii) creating meshes for fusion systems accounting for any desired level of geometric complexity and providing physics aware mesh configurations, (iii) adapting unstructured meshes to control the discretization errors, (iv) executing unstructured mesh calculations on GPU accelerated systems, (v) supporting physics/application‐specific PIC operations including surface/wall interactions of particles, (vi) providing infrastructure tools to support the interactions of solvers with unstructured meshes, and (vii) providing advanced methods for coupling plasma physics codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Partnership Center for High-Fidelity Boundary Plasma Simulation (Final Report)

Within the Partnership Center for High-Fidelity Boundary Plasma Simulation (HBPS), work at UT-Austin was aimed at improved verification, validation, and uncertainty quantification (VVUQ) for edge plasma simulations and on performing gyrokinetics simulations of pedestal instabilities and turbulence in order to expand foundational understanding of pedestal transport. Regarding VVUQ, the accomplishments can be summarized as follows. First, it was shown that the Moment Preserving Constrained Resampling technique, when applied periodically in particle-in-cell simulations in the XGC code, can dramatically improve the accuracy of the simulation at essentially equivalent computational cost. Second, a technique for estimating model correlations, which are required to solve the model selection and sample allocation problem in multifidelity UQ techniques, without sampling the highest fidelity, most computationally expensive model, was developed and demonstrated. Third, previously developed methods for estimating statistical and discretization errors were applied to numerical methods relevant to edge plasma simulations, namely in particle-in-cell-based approaches, and shown to work. Finally, benchmark studies for comparing gyrokinetic codes were developed and performed, leading to reasonable agreement between four commonly used codes. Regarding physics studies, gyrokinetic simulations to investigate microtearing modes in the DIII-D pedestal were performed using the GENE code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Comparison of Full-Field and Integrated CFD Convergence Based on Richardson Extrapolation

This work investigated the usefulness of Richardson extrapolation--based discretization error estimates across all points in a solution field to produce a spatial convergence field for a computational fluid dynamics (CFD) simulation. The presented work used previously developed methods for Richardson extrapolation to compute the convergence orders of a CFD simulation at all points of the base (coarsest) mesh solution. Three test cases of increasing complexity were considered: Poiseuille flow, incompressible flow around a sharp corner, and transonic flow over an RAE 2822 airfoil. These test cases highlighted the potential of the proposed method to identify error sources and their relation to the model system-response-quantity convergence orders. However, these test cases also revealed the immaturity of the proposed method stemming from the unreliability of computing observed convergence orders at single points. Nonetheless, the test cases highlighted that the observed convergence orders allow for a more accurate diagnosis of constructive and destructive error transport than mesh pair error estimates. In the long run, the proposed method can be a tool for developing efficient and advanced error management strategies like adaptive mesh refinement.

Weinmeister, Justin↗

DOE OSTI2025

The Fluid Dynamics Uncertainty Quantification Challenge Problem: XFOIL vs. MFOIL

Uncertainty quantification (UQ) has become more critical in aerospace engineering due to the growing dependence on computational tools for design optimization and performance analyses of aerospace vehicles. Even though the significance of UQ in assessing the credibility of computational analyses is well recognized, its costs and complexity impede its integration into standard practices, particularly in computational fluid dynamics (CFD) and other fluid analyses. This paper presents a UQ study for low-fidelity computational aerodynamics analyses with XFOIL and mfoil (i.e., the MATLAB version of XFOIL with several implementation modifications); these tools are utilized widely in both research and education. The main contributions of this paper are as follows: 1) improved precision in quantifying the uncertainty of the baseline Monte Carlo results used to benchmark surrogate modeling techniques for UQ, 2) quantification of the effect of the implementation differences between XFOIL and mfoil on solution quantities of interest (QoIs), such as lift and pitching moment coefficients, and 3) development of an open-source UQ library for use with XFOIL and mfoil, which has educational values and helps promote UQ for fluid analyses with aerospace applications. Results and discussions revolve around cases 1-4 of the challenge problem posed by the AIAA Fluid Dynamics Technical Committee’s Uncertainty Quantification Discussion Group (UQDG). In case 3, this work employs CFDverify, an open-source solution verification software, to quantify the discretization error and evaluate the extrapolated QoIs based on the grid convergence index (GCI). This UQ study differentiates itself from previous studies in the rigor of handling baseline Monte Carlo uncertainty and in including mfoil, which is a more accessible alternative to XFOIL. Finally, despite the growing computing power, low-fidelity computational tools remain valuable, such as for aerodynamic shape optimization at Mach numbers below 0.65 and low-to-mid Reynolds numbers.

Lay, Aidan S [University of Tennessee, Knoxville (↗

DOE OSTI2020

Uncertainty Assessment of CFD Investigation of the Nonlinear Difference-Frequency Wave Loads on a Semisubmersible FOWT Platform

Current mid-fidelity modeling approaches for floating offshore wind turbines (FOWTs) have been found to underpredict the nonlinear, low-frequency wave excitation and the response of semisubmersible FOWTs. To examine the cause of this underprediction, the OC6 project is using computational fluid dynamics (CFD) tools to investigate the wave loads on the OC5-DeepCwind semisubmersible, with a focus on the nonlinear difference-frequency excitation. This paper focuses on assessing the uncertainty of the CFD predictions from simulations of the semisubmersible in a fixed condition under bichromatic wave loading and on establishing confidence in the results for use in improving mid-fidelity models. The uncertainty for the nonlinear wave excitation is found to be acceptable but larger than that for the wave-frequency excitation, with the spatial discretization error being the dominant contributor. Further, unwanted free waves at the difference frequency have been identified in the CFD solution. A wave-splitting and wave load-correction procedure are presented to remove the contamination from the free waves in the results. A preliminary comparison to second-order potential-flow theory shows that the CFD model predicted significantly higher difference-frequency wave excitations, especially in surge, suggesting that the CFD results can be used to better calibrate the mid-fidelity tools.

17 WIND ENERGY↗

NASA NTRS1973

A variable step incremental procedure

Description of a variable step incremental procedure for the solution of nonlinear equations in finite element structural analysis. The proposed procedure is effective in improving the accuracy of the basic incremental technique and in providing, in addition, an accurate estimate of the discretization error. The proposed approach is highly appropriate for solving problems for which the user has no a priori estimate of the step size to use.

Thomas, G. R.↗

NASA NTRS1979

A vector-continuous loading concept for aerodynamic panel methods

An approach to the reduction of discretization errors in aerodynamic panel methods is presented. The approach is based on preventing the occurence of induced velocity singularities at panel slope discontinuities by maintaining continuity of the velocity jump vector across the panels. The approach was implemented in a two-dimensional incompressible panel method formulation and evaluated by application to several external and internal flow problems. The method is shown to exhibit a second order accuracy trend and to produce smaller errors with velocity component boundary conditions imposed on the real flow than with equipotential boundary conditions imposed on the imaginary flow behind the panels. For flows around airfoil sections with either sharp or blunt trailing edges, the method gives excellent agreement with results from a well developed finite difference method. The method is well behaved and is insensitive to irregularities in panel size distribution.

Kemp, W. B., Jr.↗

NASA NTRS1980

Vortex methods for two- and three-dimensional flow simulations

The point vortex and vortex blob methods for two dimensional flows are presented. Several results are discussed concerning the numerical analysis of the latter scheme, e.g., the preservation of globally conserved quantities and the analysis of the spatial discretization error resulting from the convection of fixed blobs of vorticity. An application to the two dimensional mixing layer is briefly described. The contour dynamics method is also discussed. The simulation of three dimensional flows with vortex methods is discussed. A natural way to represent the vorticity is in the form of closed tubes of filaments of vorticity, although other schemes are examined. Applications to aircraft trailing vortices and to a turbulent spot in a laminar boundary layer are presented. Hybrid schemes that use an Eulerian mesh to solve the Poisson equation for the velocity field are discussed. The goal of these schemes is to avoid the high cost of the Biot-Savart integration if many vortex elements are used while enjoying most of the advantages of pure Lagrangian schemes.

Leonard, A.↗

NASA NTRS1981

The use of adaptive grids in conjunction with shock capturing methods

The use of shock capturing finite-difference techniques in computing flow fields containing shocks results in a smeared or oscillatory solution in the vicinity of the shocks. This smearing or oscillatory behavior is due to the discretized form of the governing differential equations used to compute the solution. The discretization error can be reduced by a proper clustering of mesh points in the region of the shock and by using shock aligned grids. This paper uses a simple method that was developed earlier to cluster points near the shocks and serves to introduce a new method of generating a shock aligned mesh. Applications to the one-dimensional inviscid Burgers' equation, supersonic flow over a wedge with the associated straight oblique shock, one- and two-dimensional inviscid flows through an expanding duct and the problem of a straight oblique shock in a uniform supersonic freestream are presented. Significant reduction in the oscillatory behavior of the solution is demonstrated.

Rai, M. M.↗

NASA NTRS1981

The nonlinear modified equation approach to analyzing finite difference schemes

The nonlinear modified equation approach is taken in this paper to analyze the generalized Lax-Wendroff explicit scheme approximation to the unsteady one- and two-dimensional equations of gas dynamics. Three important applications of the method are demonstrated. The nonlinear modified equation analysis is used to (1) generate higher order accurate schemes, (2) obtain more accurate estimates of the discretization error for nonlinear systems of partial differential equations, and (3) generate an adaptive mesh procedure for the unsteady gas dynamic equations. Results are obtained for all three areas. For the adaptive mesh procedure, mesh point requirements for equal resolution of discontinuities were reduced by a factor of five for a 1-D shock tube problem solved by the explicit MacCormack scheme.

Klopfer, G. H.↗