Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical 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.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Hybrid iterative approach for simulation of radio-frequency fields in plasma

In this work, a novel iterative approach for solving discretized linear wave equations in a frequency domain, which combines time evolution with iterative relaxation schemes, is presented. In this hybrid approach, each iteration cycle consists of evolution of electromagnetic (EM) fields in time over a specified number of field periods followed by several iterative relaxations. Provided that there is sufficient dissipation, both the time evolution and the iterative relaxations contribute to the convergence of the EM fields to the solution of the formulated full wave boundary value problem. Time evolution rapidly distributes EM fields, propagating with group velocity, over the simulation domain, while the iterative relaxations smooth the fields, reducing the numerical errors such that iteration cycles converge to a steady state solution, approximating the solution of the formulated problem. This approach is intended for large scale simulations which are beyond the capabilities of direct solvers presently used for solving wave equations in the frequency domain. Lastly, the technique is demonstrated for solving wave equations on a regular grid using a cold plasma dielectric model with collisions for 2D modeling of EM fields in tokamak in an electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Empirical correlations between the function of entropy ( Z S ) and net artificial viscous work in a shock physics hydrocode

Entropy is a state variable that may be obtained from any thermodynamically complete equation of state (EOS). However, hydrocode calculations that output the entropy often contain numerical errors; this is not because of the EOS, but rather the solution techniques that are used in hydrocodes (especially Eulerian) such as convection, remapping, and artificial viscosity. Here, in this work, empirical correlations are investigated to reduce the errors in entropy without altering the solution techniques for the conservation of mass, momentum, and energy. Specifically, these correlations are developed for the function of entropy Z S , and they depend upon the net artificial viscous work, as determined via Sandia National Laboratories’ shock physics hydrocode CTH. These results are a continuation of a prior effort to implement the entropy-based CREST reactive burn model in CTH, and they are presented here to stimulate further interest from the shock physics community. Future work is planned to study higher-dimensional shock waves, shock wave interactions, and possible ties between the empirical correlations and a physical law.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

EXTRAPOLATION OF THE SOLAR CORONAL MAGNETIC FIELD FROM SDO/HMI MAGNETOGRAM BY A CESE-MHD-NLFFF CODE

Due to the absence of direct measurement, the magnetic field in the solar corona is usually extrapolated from the photosphere in a numerical way. At the moment, the nonlinear force-free field (NLFFF) model dominates the physical models for field extrapolation in the low corona. Recently, we have developed a new NLFFF model with MHD relaxation to reconstruct the coronal magnetic field. This method is based on CESE-MHD model with the conservation-element/solution-element (CESE) spacetime scheme. In this paper, we report the application of the CESE-MHD-NLFFF code to Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) data with magnetograms sampled for two active regions (ARs), NOAA AR 11158 and 11283, both of which were very non-potential, producing X-class flares and eruptions. The raw magnetograms are preprocessed to remove the force and then inputted into the extrapolation code. Qualitative comparison of the results with the SDO/AIA images shows that our code can reconstruct magnetic field lines resembling the EUV-observed coronal loops. Most important structures of the ARs are reproduced excellently, like the highly sheared field lines that suspend filaments in AR 11158 and twisted flux rope which corresponds to a sigmoid in AR 11283. Quantitative assessment of the results shows that the force-free constraint is fulfilled very well in the strong-field regions but apparently not that well in the weak-field regions because of data noise and numerical errors in the small currents.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Nonlinear force-free extrapolation of the coronal magnetic field based on the magnetohydrodynamic relaxation method

We develop a nonlinear force-free field (NLFFF) extrapolation code based on the magnetohydrodynamic (MHD) relaxation method. We extend the classical MHD relaxation method in two important ways. First, we introduce an algorithm initially proposed by Dedner et al. to effectively clean the numerical errors associated with ∇ · B . Second, the multigrid type method is implemented in our NLFFF to perform direct analysis of the high-resolution magnetogram data. As a result of these two implementations, we successfully extrapolated the high resolution force-free field introduced by Low and Lou with better accuracy in a drastically shorter time. We also applied our extrapolation method to the MHD solution obtained from the flux-emergence simulation by Magara. We found that NLFFF extrapolation may be less effective for reproducing areas higher than a half-domain, where some magnetic loops are found in a state of continuous upward expansion. However, an inverse S-shaped structure consisting of the sheared and twisted loops formed in the lower region can be captured well through our NLFFF extrapolation method. We further discuss how well these sheared and twisted fields are reconstructed by estimating the magnetic topology and twist quantitatively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

ACCRETION OF SUPERSONIC WINDS ONTO BLACK HOLES IN 3D: STABILITY OF THE SHOCK CONE

Using numerical simulations we present the accretion of supersonic winds onto a rotating black hole in three dimensions. We study five representative directions of the wind with respect to the axis of rotation of the black hole and focus on the evolution and stability of the high-density shock cone that is formed during the process. We explore both the regime in which the shock cone is expected to be stable in order to confirm previous results obtained with two-dimensional simulations, and the regime in which the shock cone is expected to show a flip–flop (FF) type of instability. The methods used to attempt a triggering of the instability were (i) the accumulation of numerical errors and (ii) the explicit application of a perturbation on the velocity field after the shock cone was formed. The result is negative, that is, we did not find the FF instability within the parameter space we explored, including cases that are expected to be unstable.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

On the reach of perturbative descriptions for dark matter displacement fields

We study Lagrangian Perturbation Theory (LPT) and its regularization in the Effective Field Theory (EFT) approach. We evaluate the LPT displacement with the same phases as a corresponding N-body simulation, which allows us to compare perturbation theory to the non-linear simulation with significantly reduced cosmic variance, and provides a more stringent test than simply comparing power spectra. We reliably detect a non-vanishing leading order EFT coefficient and a stochastic displacement term, uncorrelated with the LPT terms. This stochastic term is expected in the EFT framework, and, to the best of our understanding, is not an artifact of numerical errors or transients in our simulations. This term constitutes a limit to the accuracy of perturbative descriptions of the displacement field and its phases, corresponding to a 1% error on the non-linear power spectrum at k = 0.2 h{sup −1}Mpc at z = 0. Predicting the displacement power spectrum to higher accuracy or larger wavenumbers thus requires a model for the stochastic displacement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

On the regularity of the covariance matrix of a discretized scalar field on the sphere

We present a comprehensive study of the regularity of the covariance matrix of a discretized field on the sphere. In a particular situation, the rank of the matrix depends on the number of pixels, the number of spherical harmonics, the symmetries of the pixelization scheme and the presence of a mask. Taking into account the above mentioned components, we provide analytical expressions that constrain the rank of the matrix. They are obtained by expanding the determinant of the covariance matrix as a sum of determinants of matrices made up of spherical harmonics. We investigate these constraints for five different pixelizations that have been used in the context of Cosmic Microwave Background (CMB) data analysis: Cube, Icosahedron, Igloo, GLESP and HEALPix, finding that, at least in the considered cases, the HEALPix pixelization tends to provide a covariance matrix with a rank closer to the maximum expected theoretical value than the other pixelizations. The effect of the propagation of numerical errors in the regularity of the covariance matrix is also studied for different computational precisions, as well as the effect of adding a certain level of noise in order to regularize the matrix. In addition, we investigate the application of the previous results to a particular example that requires the inversion of the covariance matrix: the estimation of the CMB temperature power spectrum through the Quadratic Maximum Likelihood algorithm. Finally, some general considerations in order to achieve a regular covariance matrix are also presented.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Multiscale convergence properties for spectral approximations of a model kinetic equation

Here in this work, we prove rigorous convergence properties for a semi-discrete, moment-based approximation of a model kinetic equation in one dimension. This approximation is equivalent to a standard spectral method in the velocity variable of the kinetic distribution and, as such, is accompanied by standard algebraic estimates of the form N -q , where N is the number of modes and q > 0 depends on the regularity of the solution. However, in the multiscale setting, the error estimate can be expressed in terms of the scaling parameter ϵ, which measures the ratio of the mean-free-path to the characteristic domain length. We show that, for isotropic initial conditions, the error in the spectral approximation is $\mathcal{O}$($ϵ^{N+1}$). More surprisingly, the coefficients of the expansion satisfy super convergence properties. In particular, the error of the ℓ th coefficient of the expansion scales like $\mathcal{O}$($ϵ^{2N}$ ) when ℓ = 0 and $\mathcal{O}$($ϵ^{2N+2-ℓ}$ ) for all 1 ≤ ℓ ≤ N. This result is significant, because the low-order coefficients correspond to physically relevant quantities of the underlying system. All the above estimates involve constants depending on N, the time t, and the initial condition. We investigate specifically the dependence on N, in order to assess whether increasing N actually yields an additional factor of ϵ in the error. Numerical tests will also be presented to support the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

The Lyman α forest in optically thin hydrodynamical simulations

We study the statistics of the Lyα forest in a flat Λ cold dark matter cosmology with the N-body + Eulerian hydrodynamics code NYX.We produce a suite of simulations, covering the observationally relevant redshift range 2 ≤ z ≤ 4. We find that a grid resolution of 20 h-1kpc is required to produce 1 per cent convergence of Lyα forest flux statistics, up to k = 10 h -1 Mpc. In addition to establishing resolution requirements, we study the effects of missing modes in these simulations, and find that box sizes of L > 40h -1 Mpc are needed to suppress numerical errors to a sub-per cent level. Our optically thin simulations with the ionizing background prescription of Haardt & Madau reproduce an intergalactic medium density-temperature relation with T 0 ≈ 10 4 K and γ ≈ 1.55 at z = 2, with a mean transmitted flux close to the observed values. When using the ionizing background prescription of Faucher-Giguère et al., the mean flux is 10-15 per cent below observed values at z = 2, and a factor of 2 too small at z = 4. We show the effects of the common practice of rescaling optical depths to the observed mean flux and how it affects convergence rates. We also investigate the practice of 'splicing' results from a number of different simulations to estimate the 1D flux power spectrum and show it is accurate at the 10 per cent level. Finally, we find that collisional heating of the gas from dark matter particles is negligible in modern cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Ab-initio Calculations of Absorption Spectra of Semiconducting Nanowires Within Many Body Perturbation Theory

We present a study of the optical absorption spectra of thin silicon nanowires using many-body perturbation theory. We solve the Bethe-Salpeter equation in the static approximation using a technique that avoids explicit calculation of empty electronic states, as well as storage and inversion of the dielectric matrix. We provide a detailed assessment of the numerical accuracy of this technique, when using plane wave basis sets and periodically repeated supercells. Our calculations show that establishing numerical error bars of computed spectra is critical, in order to draw meaningful comparisons with experiments and between results obtained within different algorithms. We also discuss the influence of surface structure on the absorption spectra of nanowires with {approx_equal}1-nm diameter. Finally, we compare our calculations with those obtained within time-dependent density functional theory and find substantial differences, more pronounced than in the case of Si nanoparticles with the same diameter.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Role of Collins effect in the single spin asymmetry A{sub N} in p(transv. Polarized) p --> h X processes

The much debated issue of the transverse single spin asymmetry A{sub N} observed in the inclusive large P{sub T} production of a single hadron in pp interactions, p(transv. polarized) p --> {pi} X, is considered in a Transverse Momentum Dependent factorization scheme. A previous result stating that the maximum contribution of the Collins effect is strongly suppressed, is revisited, correcting a numerical error. New estimates are given, adopting the Collins functions recently extracted from Semi-Inclusive Deep Inelastic Scattering and e{sup +}e{sup -} data, and phenomenological consequences are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Accurate and Scalable O(N) Algorithm for First-Principles Molecular-Dynamics Computations on Large Parallel Computers

We present the first truly scalable first-principles molecular dynamics algorithm with O(N) complexity and controllable accuracy, capable of simulating systems with finite band gaps of sizes that were previously impossible with this degree of accuracy. By avoiding global communications, we provide a practical computational scheme capable of extreme scalability. Accuracy is controlled by the mesh spacing of the finite difference discretization, the size of the localization regions in which the electronic wave functions are confined, and a cutoff beyond which the components of the overlap matrix can be omitted when computing selected elements of its inverse. We demonstrate the algorithm's excellent parallel scaling for up to 101 952 atoms on 23 328 processors, with a wall-clock time of the order of 1 min per molecular dynamics time step and numerical error on the forces of less than 7x10 -4 Ha/Bohr.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Robust Solution Verification Experiments on Nonuniform Meshes

The activities of verification, validation, and uncertainty quantification (VVUQ) provide a comprehensive means to assess the credibility of computational models. Within VVUQ, solution verification assesses numerical errors and evaluates whether the simulation is sufficiently accurate for its intended applications. As computational modeling gains traction in the development of complex, high-consequence systems, the need for robust solution verification intensifies, particularly because experimental data for these systems are often limited. This work examines improvements in the robustness of Richardson extrapolation (RE), a method commonly used in solution verification to study the discretization error of computational models using a power law. Nonuniform mesh refinement is discussed alongside other pollutants that affect the robustness of the power law model. Maximum likelihood estimation (MLE) is proposed as a robust strategy to address the uncertainty generated by nonuniform mesh refinement. An exploratory computational fluid dynamics (CFD) study of a 2D planar Poiseuille flow is conducted to determine if nonuniform mesh noise can be modeled with this MLE approach for more robust RE.

Weinmeister, Justin [ORNL] (ORCID:0000000160090237↗
DOE OSTI2017

Assessment of Vertical Mesh Refinement in Concurrently Nested Large-Eddy Simulations Using the Weather Research and Forecasting Model

To facilitate multiscale simulation using the Weather Research and Forecasting Model, vertical mesh refinement for one-way concurrent nested simulation was recently introduced. Grid refinement in the vertical dimension removes the requirement of different grid aspect ratios on the bounding versus the nested domain, such that results from refinement are in the horizontal directions only, and thereby can also reduce numerical errors on the bounding domain arising from large aspect ratios in the presence of complex terrain. Herein, the impacts of vertical grid refinement on the evolving downstream flow in nested large-eddy simulations are evaluated in relation to other model configuration choices, including turbulence subfilter-scale (SFS) stress models, mesh configuration, and alternative methods for calculating several near-surface flow parameters. Although vertical nesting requires coarsening of the vertical grid on the bounding domain, leading to a smaller range of resolved turbulence scales in the nest’s lateral boundary conditions, parameter values within the nested domains are generally only minimally impacted, relative to nesting using the same vertical grid on each domain. Two dynamic SFS models examined herein generally improved the simulated mean wind speed, turbulence kinetic energy, stresses and spectra, on both domains, and accelerated equilibration rates within nested domains, relative to two constant coefficient models. A new method of extrapolating horizontal velocity components to near-surface locations at nested domain lateral boundaries, and a correction to the calculation of deformation elements near the surface, are each shown to slightly alter the mean parameter values, yet only minimally impact equilibration rates within the nested domain.

14 SOLAR ENERGY↗
DOE OSTI2018

Topographic Effects on Radiation in the WRF Model with the Immersed Boundary Method: Implementation, Validation, and Application to Complex Terrain

Abstract Topographic effects on radiation, including both topographic shading and slope effects, are included in the Weather Research and Forecasting (WRF) Model, and here they are made compatible with the immersed boundary method (IBM). IBM is an alternative method for representing complex terrain that reduces numerical errors over sloped terrain, thus extending the range of slopes that can be represented in WRF simulations. The implementation of topographic effects on radiation is validated by comparing land surface fluxes, as well as temperature and velocity fields, between idealized WRF simulations both with and without IBM. Following validation, the topographic shading implementation is tested in a semirealistic simulation of flow over Granite Mountain, Utah, where topographic shading is known to affect downslope flow development in the evening. The horizontal grid spacing is 50 m and the vertical grid spacing is approximately 8–27 m near the surface. Such a case would fail to run in WRF with its native terrain-following coordinates because of large local slope values reaching up to 55°. Good agreement is found between modeled surface energy budget components and observations from the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) program at a location on the east slope of Granite Mountain. In addition, the model captures large spatiotemporal inhomogeneities in the surface sensible heat flux that are important for the development of thermally driven flows over complex terrain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Mesoscale to Microscale Simulations over Complex Terrain with the Immersed Boundary Method in the Weather Research and Forecasting Model

Improvements to the Weather Research and Forecasting (WRF) Model are made to enable multiscale simulations over highly complex terrain with dynamically downscaled boundary conditions from the mesoscale to the microscale. Over steep terrain, the WRF Model develops numerical errors that are due to grid deformation of the terrain-following coordinates. An alternative coordinate system, the immersed boundary method (IBM), has been implemented into WRF, allowing for simulations over highly complex terrain; however, the new coordinate system precluded nesting within mesoscale simulations using WRF’s native terrain-following coordinates. Here, the immersed boundary method and WRF’s grid-nesting framework are modified to seamlessly work together. This improved framework for the first time allows for large-eddy simulation over complex (urban) terrain with IBM to be nested within a typical mesoscale WRF simulation. Simulations of the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma, are performed using a multiscale five-domain nested configuration, spanning horizontal grid resolutions from 6 km to 2 m. These are compared with microscale-only simulations with idealized lateral boundary conditions and with observations of wind speed/direction and SF 6 concentrations from a controlled release from intensive observation period 3. The multiscale simulation, which is configured independent of local observations, shows similar model skill predicting wind speed/direction and improved skill predicting SF 6 concentrations when compared with the idealized simulations, which require use of observations to set mean flow conditions. Use of this improved multiscale framework shows promise for enabling large-eddy simulation over highly complex terrain with dynamically downscaled boundary conditions from mesoscale models.

42 ENGINEERING↗
DOE OSTI2013

Statistical Mechanics Modeling of Mesoscale Deformation in Metals

The research under this project focused on a theoretical and computational modeling of dislocation dynamics of mesoscale deformation of metal single crystals. Specifically, the work aimed to implement a continuum statistical theory of dislocations to understand strain hardening and cell structure formation under monotonic loading. These aspects of crystal deformation are manifestations of the evolution of the underlying dislocation system under mechanical loading. The project had three research tasks: 1) Investigating the statistical characteristics of dislocation systems in deformed crystals. 2) Formulating kinetic equations of dislocations and coupling these kinetics equations and crystal mechanics. 3) Computational solution of coupled crystal mechanics and dislocation kinetics. Comparison of dislocation dynamics predictions with experimental results in the area of statistical properties of dislocations and their field was also a part of the proposed effort. In the first research task, the dislocation dynamics simulation method was used to investigate the spatial, orientation, velocity, and temporal statistics of dynamical dislocation systems, and on the use of the results from this investigation to complete the kinetic description of dislocations. The second task focused on completing the formulation of a kinetic theory of dislocations that respects the discrete nature of crystallographic slip and the physics of dislocation motion and dislocation interaction in the crystal. Part of this effort also targeted the theoretical basis for establishing the connection between discrete and continuum representation of dislocations and the analysis of discrete dislocation simulation results within the continuum framework. This part of the research enables the enrichment of the kinetic description with information representing the discrete dislocation systems behavior. The third task focused on the development of physics-inspired numerical methods of solution of the coupled dislocation kinetics and crystal mechanics framework. To a large extent, this task has also been successfully started. We have developed a custom finite-element approach with mesh points being a subset of the underlying crystal structure. When used to predict the evolution of the dislocation system, the planar motion of dislocations is naturally captured for all slip systems, thus minimizing numerical errors and providing simple ways to investigate cross slip and dislocation reactions. Preliminary results in this direction show that we are closer than ever in building a predictive framework for dislocation dynamics and mesoscale plasticity based on the first principles of dislocation dynamics. The rest of the report gives and overview of the research performed under this project and highlights the key results and open questions left for future investigations.

36 MATERIALS SCIENCE↗