Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical accuracy”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

ON THE ANALATICAL SOLUTIONS AND NUMERICAL VERIFICATIONS OF THE TWO-PHASE WATER FAUCET PROBLEM

The two-phase water faucet problem is one of the benchmark problems originally proposed by Ransom. The test consists of a vertical pipe, which is initially filled with a uniform column of liquid moving with an initial velocity and an annulus of gas sitting still. At the top boundary, liquid is supplied at the same velocity as the initial velocity, and the bottom of the pipe opens to the ambient. Due to the gravity effect, the liquid column accelerates and becomes thinner. With simplifications, such as massless gas phase and no wall and interfacial frictions, an analytical solution could be obtained for the transient liquid velocity and void fraction distribution. The test has been widely used for the purposes of code assessment, benchmark and numerical verifications. This exact solution was used in this work to verify the numerical accuracy of a high-resolution spatial discretization scheme. The high-resolution spatial discretization scheme is based on the commonly used staggered grid mesh, and slope limiter was introduced to achieve high resolutions in space. It was found that, at the steady state, an expected second-order spatial accuracy could not be achieved. A further investigation showed that this exact solution does not actually satisfy the commonly used two-fluid single pressure two-phase flow equations. More importantly, we have also found that a steady-state exact solution could be obtained for the water faucet problem. In this work, we have demonstrated the derivation of the steady-state exact solution for the water faucet problem. The exact solution was then used to verify a first-order and a high-resolution spatial discretization scheme, from which expected order of accuracy was achieved by doing mesh refinement.

42 ENGINEERING↗
DOE OSTI2011

Stochastic weighted particle methods for population balance equations

Highlights: {yields} Weight transfer functions for Monte Carlo simulation of coagulation. {yields} Efficient support for single-particle growth processes. {yields} Comparisons to analytic solutions and soot formation problems. {yields} Better numerical accuracy for less common particles. - Abstract: A class of coagulation weight transfer functions is constructed, each member of which leads to a stochastic particle algorithm for the numerical treatment of population balance equations. These algorithms are based on systems of weighted computational particles and the weight transfer functions are constructed such that the number of computational particles does not change during coagulation events. The algorithms also facilitate the simulation of physical processes that change single particles, such as growth, or other surface reactions. Four members of the algorithm family have been numerically validated by comparison to analytic solutions to simple problems. Numerical experiments have been performed for complex laminar premixed flame systems in which members of the class of stochastic weighted particle methods were compared to each other and to a direct simulation algorithm. Two of the weighted algorithms have been shown to offer performance advantages over the direct simulation algorithm in situations where interest is focused on the larger particles in a system. The extent of this advantage depends on the particular system and on the quantities of interest.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Numerical analysis of laminar and turbulent incompressible flows using the finite element Fluid Dynamics Analysis Package (FIDAP)

The purpose of the study is the evaluation of the numerical accuracy of FIDAP (Fluid Dynamics Analysis Package). Accordingly, four test problems in laminar and turbulent incompressible flows are selected and the computational results of these problems compared with other numerical solutions and/or experimental data. These problems include: (1) 2-D laminar flow inside a wall-driven cavity; (2) 2-D laminar flow over a backward-facing step; (3) 2-D turbulent flow over a backward-facing step; and (4) 2-D turbulent flow through a turn-around duct.

Sohn, Jeong L.↗
DOE OSTI2020

Scalability of OpenFOAM Density-Based Solver with Runge–Kutta Temporal Discretization Scheme

Compressible density-based solvers are widely used in OpenFOAM, and the parallel scalability of these solvers is crucial for large-scale simulations. In this paper, we report our experiences with the scalability of OpenFOAM’s native rhoCentralFoam solver, and by making a small number of modifications to it, we show the degree to which the scalability of the solver can be improved. The main modification made is to replace the first-order accurate Euler scheme in rhoCentralFoam with a third-order accurate, four-stage Runge-Kutta or RK4 scheme for the time integration. The scaling test we used is the transonic flow over the ONERA M6 wing. This is a common validation test for compressible flows solvers in aerospace and other engineering applications. Numerical experiments show that our modified solver, referred to as rhoCentralRK4Foam, for the same spatial discretization, achieves as much as a 123.2% improvement in scalability over the rhoCentralFoam solver. As expected, the better time resolution of the Runge–Kutta scheme makes it more suitable for unsteady problems such as the Taylor–Green vortex decay where the new solver showed a 50% decrease in the overall time-to-solution compared to rhoCentralFoam to get to the final solution with the same numerical accuracy. Finally, the improved scalability can be traced to the improvement of the computation to communication ratio obtained by substituting the RK4 scheme in place of the Euler scheme. All numerical tests were conducted on a Cray XC40 parallel system, Theta, at Argonne National Laboratory.

Li, Sibo↗
DOE OSTI2013

Guiding Center Equations of High Accuracy

Guiding center simulations are an important means of predicting the effect of resistive and ideal magnetohydrodynamic instabilities on particle distributions in toroidal magnetically confined thermonuclear fusion research devices. Because saturated instabilities typically have amplitudes of δ B/B of a few times 10 -4 numerical accuracy is of concern in discovering the effect of mode particle resonances. We develop a means of following guiding center orbits which is greatly superior to the methods currently in use. In the presence of ripple or time dependent magnetic perturbations both energy and canonical momentum are conserved to better than one part in 10 14 , and the relation between changes in canonical momentum and energy is also conserved to very high order.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Memory-efficient emulation of physical tabular data using quadtree decomposition

Computationally expensive functions are sometimes replaced in simulations with an emulator that approx-imates the true function (e.g., equations of state, wavelength-dependent opacity, or composition-dependent materials properties). For functions that have a constrained domain of interest, this can be done by discretizing the domain and performing a local interpolation on the tabulated function values of each local domain. For these so-called tabular data methods, the method of discretizing the domain and mapping the input space to each subdomain can drastically influence the memory and computational costs of the emulator. This is especially true for functions that vary drastically in different regions. We present a method for domain discretization and mapping that utilizes quadtrees, which results in significant reductions in the size of the emulator with minimal increases to computational costs or loss of global accuracy. We apply our method to the electron-positron Helmholtz free energy equation of state and show over an order of magnitude reduction in memory costs for reasonable levels of numerical accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Nonlinear Dynamic Characteristics and Control of Giant Magnetostrictive Ultrasonic Transducer

The nonlinear dynamic characteristics of giant magnetostrictive ultrasonic transducer (GMUT) are researched in this paper, and the optimal control strategy is proposed. New nonlinear differential terms are introduced to explain the hysteretic phenomena of giant magnetostrictive materials (GMM), and the nonlinear dynamic model of a giant magnetostrictive ultrasonic transducer subjected to harmonic magnetic field and stochastic disturbance is developed. The system’s nonlinear dynamic responses are discussed, and the optimal control strategy is proposed to improve the system’s accuracy. Numerical simulation and experiment results show that there are abundant nonlinear dynamic behaviors in the system; the stability of the system varies with the parameters, which causes stochastic bifurcation; the system’s accuracy is improved through stochastic optimal control. The results of this study are helpful for the application of giant magnetostrictive ultrasonic transducer in engineering fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Accuracy of Lagrange-sinc functions as a basis set for electronic structure calculations of atoms and molecules

We developed a self-consistent field program based on Kohn-Sham density functional theory using Lagrange-sinc functions as a basis set and examined its numerical accuracy for atoms and molecules through comparison with the results of Gaussian basis sets. The result of the Kohn-Sham inversion formula from the Lagrange-sinc basis set manifests that the pseudopotential method is essential for cost-effective calculations. The Lagrange-sinc basis set shows faster convergence of the kinetic and correlation energies of benzene as its size increases than the finite difference method does, though both share the same uniform grid. Using a scaling factor smaller than or equal to 0.226 bohr and pseudopotentials with nonlinear core correction, its accuracy for the atomization energies of the G2-1 set is comparable to all-electron complete basis set limits (mean absolute deviation ≤1 kcal/mol). The same basis set also shows small mean absolute deviations in the ionization energies, electron affinities, and static polarizabilities of atoms in the G2-1 set. In particular, the Lagrange-sinc basis set shows high accuracy with rapid convergence in describing density or orbital changes by an external electric field. Moreover, the Lagrange-sinc basis set can readily improve its accuracy toward a complete basis set limit by simply decreasing the scaling factor regardless of systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1969

Numerical solutions of differential equations

Various numerical methods for solving differential equations were analyzed and refined in an effort to develop a method which was adaptable to a large class of problems. The prime capabilities of the method included accuracy, numerical stability, and economic use of computer time. In multistep processes the corrector was changed at each step.

Wesson, J. R.↗
NASA NTRS1980

Implicit Finite-Difference Simulations of Three-Dimensional Compressible Flow

An implicit finite-difference procedure for unsteady three-dimensional flow capable of handling arbitrary geometry through the use of general coordinate transformations is described. Viscous effects are optionally incorporated with a "thin-layer" approximation of the Navier-Stokes equations. An implicit approximate factorization technique is employed so that the small grid sizes required for spatial accuracy and viscous resolution do not impose stringent stability limitations. Results obtained from the program include transonic inviscid or viscous solutions about simple body configurations. Comparisons with existing theories and experiments are made. Numerical accuracy and the effect of three-dimensional coordinate singularities are also discussed.

Pulliam, Thomas H.↗
NASA NTRS2002

Application of the generalized vertical coordinate ocean model for better representing satellite data

It is found that two adaptive parametric functions can be introduced into the basic ocean equations for utilizing the optimal or hybrid features of commonly used z-level, terrain- following, isopycnal, and pressure coordinates in numerical ocean models. The two parametric functions are formulated by combining three techniques: the arbitrary vertical coordinate system of Kasahara (1 974), the Jacobian pressure gradient formulation of Song (1 998), and a newly developed metric factor that permits both compressible (non-Boussinesq) and incompressible (Boussinesq) approximations. Based on the new formulation, an adaptive modeling strategy is proposed and a staggered finite volume method is designed to ensure conservation of important physical properties and numerical accuracy. Implementation of the combined techniques to SCRUM (Song and Haidvogel1994) shows that the adaptive modeling strategy can be applied to any existing ocean model without incurring computational expense or altering the original numerical schemes. Such a generalized coordinate model is expected to benefit diverse ocean modelers for easily choosing optimal vertical structures and sharing modeling resources based on a common model platform. Several representing oceanographic problems with different scales and characteristics, such as coastal canyons, basin-scale circulation, and global ocean circulation, are used to demonstrate the model's capability for multiple applications. New results show that the model is capable of simultaneously resolving both Boussinesq and non-Boussinesq, and both small- and large-scale processes well. This talk will focus on its applications of multiple satellite sensing data in eddy-resolving simulations of Asian Marginal Sea and Kurosio. Attention will be given to how Topex/Poseidon SSH, TRMM SST; and GRACE ocean bottom pressure can be correctly represented in a non- Boussinesq model.

numerical ocean model generalized coordinate syste↗
DOE OSTI2025

Enhancing the accuracy of XPS calculations: Exploring hybrid basis set schemes for CVS-EOMIP-CCSD calculations

Reliable computational methodologies and basis sets for modeling x-ray spectra are essential for extracting and interpreting electronic and structural information from experimental x-ray spectra. In particular, the trade-off between numerical accuracy and computational cost due to the size of the basis set is a major challenge, since molecular orbitals undergo extreme relaxation in the core-hole state. To gain clarity on the changes in electronic structure induced by the formation of a core-hole, the use of sufficiently flexible basis for expanding the orbitals, particularly for the core region, has been shown to be essential. This work focuses on the refinement of core-hole ionized state calculations using the equation-of-motion coupled cluster family of methods through an extensive analysis on the effectiveness of “hybrid” and mixed basis sets. In this investigation, we utilize the CVS-EOMIP-CCSD method in combination and construct hybrid basis sets piecewise from readily available Dunning’s correlation consistent basis sets in order to calculate x-ray ionization energies (IEs) for a set of small gas phase molecules. Our results provide insights into the impact of basis sets on the CVS-EOMIP-CCSD calculations of K-edge IEs of first-row p-block elements. Furthermore, these insights enable us to understand more about the basis set dependence of the core IEs computed and allow us to establish a protocol for deriving reliable and cost-effective theoretical estimates for computing IEs of small molecules containing such elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1978

On implicit finite-difference simulations of three-dimensional flow

An implicit finite-difference procedure for unsteady three-dimensional flow capable of handling arbitrary geometry through the use of general coordinate transformations is described. Viscous effects are optionally incorporated with a 'thin layer' approximation of the Navier-Stokes equations. An implicit approximate factorization technique is employed so that the small grid sizes required for spatial accuracy and viscous resolution do not impose stringent stability limitations. Results obtained from the program include transonic inviscid and laminar-turbulent solutions about simple body configurations. Comparisons with existing theories and experiments are made. Numerical accuracy and the effect of three-dimensional coordinate singularities are also discussed.

Pulliam, T. H.↗
NASA NTRS1992

Initial weather regimes as predictors of numerical 30-day mean forecast accuracy

Thirty 30-day mean 500-mb-height anomaly forecasts generated by the NCAR Community Climate Model (CCM) for the year 1978 are examined in order to determine if the forecast accuracy can be estimated with the initial conditions. The initial weather regimes were defined in such a way that the regimes could discriminate between the best and the worst 30-day mean forecasts run from the initial fields in this data set. On the basis of the CCM experiments, it is suggested that the accuracy of numerical 30-day mean forecasts may depend upon the accuracy with which the cyclones and their interactions with the planetary scale are predicted early in the forecast cycle, and that this accuracy may depend upon the initial conditions.

Colucci, Stephen J.↗
DOE OSTI2024

Bi-fidelity variational auto-encoder for uncertainty quantification

Quantifying the uncertainty of quantities of interest (QoIs) from physical systems is a primary objective in model validation. However, achieving this goal entails balancing the need for computational efficiency with the requirement for numerical accuracy. To address this trade-off, we propose a novel bi-fidelity formulation of variational auto-encoders (BF-VAE) designed to estimate the uncertainty associated with a QoI from low-fidelity (LF) and high-fidelity (HF) samples of the QoI. Here, this model allows for the approximation of the statistics of the HF QoI by leveraging information derived from its LF counterpart. Specifically, we design a bi-fidelity auto-regressive model in the latent space which is integrated within the VAE’s probabilistic encoder–decoder structure. An effective algorithm is proposed to maximize the variational lower bound of the HF log-likelihood in the presence of limited HF data, resulting in the synthesis of HF realizations with a reduced computational cost. Additionally, we introduce the concept of the bi-fidelity information bottleneck (BF-IB) to provide an information-theoretic interpretation of the proposed BF-VAE model. Our numerical results demonstrate that the BF-VAE leads to considerably improved accuracy, as compared to a VAE trained using only HF data, when limited HF data is available.

42 ENGINEERING↗
DOE OSTI2018

A high-order multiscale finite-element method for time-domain acoustic-wave modeling

Accurate and efficient wave equation modeling is vital for many applications in such as acoustics, electromagnetics, and seismology. However, solving the wave equation in large-scale and highly heterogeneous models is usually computationally expensive because the computational cost is directly proportional to the number of grids in the model. We develop a novel high-order multiscale finite-element method to reduce the computational cost of time-domain acoustic-wave equation numerical modeling by solving the wave equation on a coarse mesh based on the multiscale finite-element theory. In contrast to existing multiscale finite-element methods that use only first-order multiscale basis functions, our new method constructs high-order multiscale basis functions from local elliptic problems which are closely related to the Gauss–Lobatto–Legendre quadrature points in a coarse element. Essentially, these basis functions are not only determined by the order of Legendre polynomials, but also by local medium properties, and therefore can effectively convey the fine-scale information to the coarse-scale solution with high-order accuracy. Numerical tests show that our method can significantly reduce the computation time while maintain high accuracy for wave equation modeling in highly heterogeneous media by solving the corresponding discrete system only on the coarse mesh with the new high-order multiscale basis functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A General Framework for Progressive Data Compression and Retrieval

In scientific simulations, observations, and experiments, the transfer of data to and from disk and across networks has become a major bottleneck for data analysis and visualization. Compression techniques have been employed to tackle this challenge, but traditional lossy methods often demand conservative error tolerances to meet the numerical accuracy requirements of both anticipated and unknown data analysis tasks. Progressive data compression and retrieval has emerged as a promising solution, where each analysis task dictates its own accuracy needs. However, few analysis algorithms inherently support progressive data processing, and adapting compression techniques, file formats, client/server frameworks, and APIs to support progressivity can be challenging. Here, this paper presents a framework that enables progressive-precision data queries for any data compressor or numerical representation. Our strategy hinges on a multi-component representation that successively reduces the error between the original and compressed field, allowing each field in the progressive sequence to be expressed as a partial sum of components. We have implemented this approach with four established scientific data compressors and assessed its effectiveness using real-world data sets from the SDRBench collection. The results show that our framework competes in accuracy with the standalone compressors it is based upon. Additionally, (de)compression time is proportional to the number of components requested by the user. Finally, our framework allows for fully lossless compression using lossy compressors when a sufficient number of components are employed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Modeling turbulent, reacting flow

Several of the approximations or models involved in the development of a numerical combustor flow code are examined. In the first section, the importance of numerical accuracy is illustrated, and the impact that improved-accuracy schemes have on slowing convergence is demonstrated. Solution algorithms that can speed convergence are discussed and some performance features of these algorithms are illustrated. A sample calculation displaying the importance of boundary conditions on a three-dimensional numerical prediction is presented. The inaccuracy of a current turbulence model in highly turbulent (nonequilibrium) regions is described. Finally, the surprisingly good performance of a six-flux model in describing radiation heat transfer is displayed. In all the areas examined, continued research is still needed, but valuable engineering tools are available today.

Claus, Russell W.↗