Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Theoretical and numerical studies of inverse source problem for the linear parabolic equation with sparse boundary measurements

We consider the inverse source problem in the parabolic equation, where the unknown source possesses the semi-discrete formulation. Theoretically, we prove that the flux data from any nonempty open subset of the boundary can uniquely determine the semi-discrete source. This means the observed area can be extremely small, and that is the reason we call it sparse boundary data. For the numerical reconstruction, we formulate the problem from the Bayesian sequential prediction perspective and conduct the numerical examples which estimate the space-time-dependent source state by state. To better demonstrate the method’s performance, we solve two common multiscale problems from two models with a long source sequence. The numerical results illustrate that the inversion is accurate and efficient.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A numerical Poisson solver with improved radial solutions for a self-consistent locally scaled self-interaction correction method

Abstract The universal applicability of density functional approximations is limited by self-interaction error made by these functionals. Recently, a novel one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities was proposed. The locally scaled SIC (LSIC) method is exact for the one-electron densities, and unlike the well-known Perdew–Zunger SIC (PZSIC) method recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to PZSIC method as a special case when isoorbital indicator is set to the unity. Here, we present a numerical scheme that we have adopted to evaluate the Coulomb potential of the electron density scaled by the iso-orbital indicator required for the self-consistent LSIC calculations. After analyzing the behavior of the finite difference method (FDM) and the green function solution to the radial part of the Poisson equation, we adopt a hybrid approach that uses the FDM for the Coulomb potential due to the monopole and the GF for all higher-order terms. The performance of the resultant hybrid method is assessed using a variety of systems. The results show improved accuracy than earlier numerical schemes. We also find that, even with a generic set of radial grid parameters, accurate energy differences can be obtained using a numerical Coulomb solver in standard density functional studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A survey of numerical linear algebra methods utilizing mixed-precision arithmetic

The efficient utilization of mixed-precision numerical linear algebra algorithms can offer attractive acceleration to scientific computing applications. Especially with the hardware integration of low-precision special-function units designed for machine learning applications, the traditional numerical algorithms community urgently needs to reconsider the floating point formats used in the distinct operations to efficiently leverage the available compute power. In this study, we provide a comprehensive survey of mixed-precision numerical linear algebra routines, including the underlying concepts, theoretical background, and experimental results for both dense and sparse linear algebra problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Efficient Multi-Scale, Thermo-Hydro-Mechanical Numerical Model for Simulation of Long-Term Stability of Rooms for Nuclear Waste Disposal

This project demonstrated the feasibility of developing a numerical modeling tool that provides an efficient and unique approach to analyze the geomechanical performance of a deep geological repository (DGR) for high-level nuclear waste. The numerical modeling tool is physics-based and models the thermo-hydro-mechanical (THM) processes during the required time (up to one million years) for a DGR. The unique capability of the developed code is explicit representation of fracturing and a discrete fracture network (DFN) in the model that includes the entire repository layout for simulated time of one million years. Therefore, it can potentially address one of the challenges of analysis of a DGR, which is efficient numerical simulation of a multi-scale, physics-based models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Crack-Free Joining of High-Strength AA7055 Sheets by Friction Based Self-Piercing Riveting with the Aid of Numerical Design

Unique friction-based self-piercing riveting (F-SPR) was employed to join high-strength, low-ductility aluminum alloy 7055 for lightweight vehicle applications. This study aimed to maximize the joint strength of the AA7055 F-SPR joint while avoiding cracking issues due to low ductility at room temperature. A fully coupled Eulerian–Lagrangian (CEL) model was employed to predict the process temperature during F-SPR, and the temperature field was then mapped onto a 2D axisymmetric equivalent model for accelerated numerical analysis. The geometry, dimensions, and material strength of the rivet, as well as the depth of the die cavity and plunging depth, were investigated to enhance joint formation. Also, a static finite-element analysis model was developed to predict and analyze the stress distribution in the rivet under different mechanical testing loading conditions. Overall, the numerical model showed good agreement with the experiment results, such as joint formation and mechanical joint strength. With the aid of virtual fabrication through numerical modeling, the joint design iterations and process development time of F-SPR were greatly reduced regarding the goal of lightweight, high-strength aluminum joining.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessment of Numerical Diffusion in NRELAP5 Code for Density Wave Oscillations Applications

Density wave oscillation (DWO) in a boiling channel can be a delicate phenomenon, and the precision of numerical tools can influence predictions. In this paper, relevant numerical diffusion was analyzed for the finite difference schemes used in the NRELAP5 code. The numerical diffusion in NRELAP5 was investigated under both single-phase and two-phase conditions. A unique perturbation technique called the “V” ramping approach was introduced, and the NRELAP5-predicted DWO results were then compared to test data.

Numerical diffusion, NRELAP5, RELAP5-3D, density w↗
DOE OSTI2022

Removing numerical pathologies in a turbulence parameterization through convergence testing

Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e. clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, this noise can be amplified by an advective term in CLUBB's background diffusion. Smoothing the limiters and removing the advective component of the background diffusion reduces the noise and restores the expected first-order convergence in CLUBB's solutions. These model reformulations improve the results at coarser, near-operational grid spacing and time step in cumulus cloud and dry turbulence tests. In addition, convergence testing is proved to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

Mathematics and Computing, Geosciences↗
DOE OSTI2021

Propagation of Input Uncertainties in Numerical Simulations of Laser Powder Bed Fusion

Laser powder bed fusion has the potential of redefining state-of-the-art processing and production methods, but defect formation and inconsistent build quality have limited the implementation of this process. Numerical models are widely used to study this process and predict the formation of these defects. Presently, the uncertainties of model input parameters and thermophysical properties used by these numerical simulations have not been investigated. In the present study, the uncertainty in these input parameters and material properties are quantified for laser powder bed fusion, with and without a simulated powder bed, to determine their influence on the predictive accuracy of an experimentally validated numerical model. Accounting for all possible sources of uncertainty quickly becomes computationally expensive on account of the curse of dimensionality. Uncertainty in laser absorption, solid, and liquid specific heat of the metal were found to have the largest effect on model prediction reliability with or without the use of a powder bed. Results also illustrate that accounting for these three uncertain parameters still captures the majority of model prediction uncertainty. Furthermore, the methodology of this study may be used to understand the uncertainty in as-built microstructure through propagation to microstructure prediction models, or applied under processing conditions where high Péclet numbers are observed and the thermal convection and fluid flow within the molten pool are substantial.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Numerical and experimental analysis of precooled desiccant system

This paper presents a numerical and experimental study of four novel, desiccant air-conditioning systems based on pre-cooling of the ambient air. The purpose of the study is to: 1) establish which configuration is most profitable for the future studies through numerical and experimental methods; 2) verify the concept of precooled system and confirm that it is working in experimental conditions (for this purpose two prototypes of evaporative cooling units were developed and compared for different operation arrangements with desiccant wheel); 3) Analyze practical aspects of the systems, which includes water consumption analysis, economic analysis and establishing operational modes during a cooling season in temperate climate conditions (the proposed solution has a significant advantage during operation in a temperate climate, since it can utilize its first evaporative cooler the only source of cooling energy when the humidity of outdoor air allows to to meet the latent cooling load). Four arrangements of the system with different dew point evaporative air cooler prototypes were selected and developed for comparative analysis. The system was simulated as a whole, also, each component was individually tested on measurement station. develop the optimal arrangement of the system through numerical and experimental analysis. Four different performance indicators were considered: thermal COP (Coefficient Of Performance), the ERR (Energy Efficiency Ratio), coefficient a describing the potential in compensating indoor cooling loads, water consumption ratio. Also costs of the systems were compared through economic analysis. It was established that arrangement which is characterized by best-averaged performance in terms of all considered factors is a unit combined with two different types of evaporative coolers (marked as System B)- this is caused by the favorable arrangement of two units, where one is better suited to the pre-cooler role, while the other is better suited as the post-cooler.

42 ENGINEERING↗
DOE OSTI2021

On the competitive antagonism effect in combined chloride-sulfate attack: A numerical exploration

Highlights: • A chemo-physical-mechanical method based on thermodynamic modelling is implemented. • The contradicting findings from reported experiments are numerically explored. • The competitive antagonism effect is time-dependent rather than a conclusive term. • The competitive antagonism effect is quantifiable with a novel assessment criterion. It is of tremendous difficulty to draw a conclusive verdict on the complex process of mutual interference regarding the combined chloride-sulfate attack on concrete, where many have explored the topic, but contradicting results have been reported. In this paper, by taking advantages of the chemo-physical-mechanical method, based on the robust technique of thermodynamic modelling, the competitive antagonism effect in combined chloride-sulfate attack is investigated. Firstly, a range of reported tests are numerically modelled to demonstrate the effectiveness of the employed method in interpreting those opposing and diverse experimental outcomes. More importantly, a fresh conjecture incorporating a novel assessment criterion for the competitive antagonism effect is proposed and supported through a carefully designed illustrative example. Based on the numerical exploration, unique findings and implications for future engineering practices are revealed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A hybrid meshfree discretization to improve the numerical performance of peridynamic models

Efficient and accurate calculation of spatial integrals is of major interest in the numerical implementation of peridynamics (PD). The standard way to perform this calculation is a particle-based approach that discretizes the strong form of the PD governing equation. This approach has rapidly been adopted by the PD community since it offers some advantages. Additionally, it is computationally cheaper than other available schemes, can conveniently handle material separation, and effectively deals with nonlinear PD models. Nevertheless, PD models are still computationally very expensive compared with those based on the classical continuum mechanics theory, particularly for large-scale problems in three dimensions. This results from the nonlocal nature of the PD theory which leads to interactions of each node of a discretized body with multiple surrounding nodes. Here, we propose a new approach to significantly boost the numerical efficiency of PD models. We propose a discretization scheme that employs a simple collocation procedure and is truly meshfree; i.e., it does not depend on any background integration cells. In contrast to the standard scheme, the proposed scheme requires a much smaller set of neighboring nodes (keeping the same physical length scale) to achieve a specific accuracy and is thus computationally more efficient. Our new scheme is applicable to the case of linear PD models and within neighborhoods where the solution can be approximated by smooth basis functions. Therefore, to fully exploit the advantages of both the standard and the proposed schemes, a hybrid discretization is presented that combines both approaches within an adaptive framework. The high performance of the developed framework is illustrated by several numerical examples, including brittle fracture and corrosion problems in two and three dimensions.

42 ENGINEERING↗
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

MITNS: Multiple-Ion Transport Numerical Solver for magnetized plasmas

MITNS (Multiple-Ion Transport Numerical Solver) is a new numerical tool designed to perform 1D simulations of classical cross-field transport in magnetized plasmas. Furthermore, its detailed treatment of multi-species effects makes it a unique tool in the field. We describe the physical model it simulates, as well as its numerical implementation and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical modeling of hydrogen mixing in a direct-injection engine fueled with gaseous hydrogen

Hydrogen is considered as one of the most promising options to achieve effective decarbonization of the energy and transportation sectors. As such, it has recently been receiving increasing attention because of its promising potential as an energy carrier for advanced energy and propulsion systems. With a focus on internal combustion engines, direct injection (DI) of gaseous hydrogen during the compression stroke offers great potential for high engine efficiency and specific power while reducing the risk of backfiring and pre-ignition issues. Therefore, many experimental and numerical efforts have recently been dedicated to understanding the physical and chemical behaviors of hydrogen in engine during mixing and combustion. This study focuses on computational fluid dynamics (CFD) modeling of the hydrogen DI process in a hydrogen optical research engine. Under the conditions studied, gaseous hydrogen is injected into the combustion chamber via a centrally located single-hole injector at a pressure of 100 bar. Two configurations, namely low-and high-tumble, are investigated to understand the impact of different in-cylinder flow patterns on the fuel-air mixture preparation. Simulations are carried out using the commercial CFD software CONVERGE. Here, the in-cylinder turbulence is modeled with an unsteady Reynolds-averaged Navier-Stokes (URANS) formulation closed by the renormalization group (RNG) k-ε model. Several numerical methods and model constants, including but not limited to turbulent Schmidt number, are evaluated. The numerical results are systematically compared against experimental measurements of velocity and hydrogen concentration fields on the vertical center plane to assess the performance of the CFD model, unveil the physics of hydrogen mixing, and establish best practices for modeling hydrogen DI under relatively high injection pressure conditions.

33 ADVANCED PROPULSION SYSTEMS↗