Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical solver”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling the TetraSpar Floating Offshore Wind Turbine Foundation as a Flexible Structure in OrcaFlex and OpenFAST

Floating offshore wind turbine technology has seen an increasing and continuous development in recent years. When designing the floating platforms, both experimental and numerical tools are applied, with the latter often using time-domain solvers based on hydro-load estimation from a Morison approach or a boundary element method. Commercial software packages such as OrcaFlex, or open-source software such as OpenFAST, are often used where the floater is modeled as a rigid six degree-of-freedom body with loads applied at the center of gravity. However, for final structural design, it is necessary to have information on the distribution of loads over the entire body and to know local internal loads in each component. This paper uses the TetraSpar floating offshore wind turbine design as a case study to examine new modeling approaches in OrcaFlex and OpenFAST that provide this information. The study proves the possibility of applying the approach and the extraction of internal loads, while also presenting an initial code-to-code verification between OrcaFlex and OpenFAST. As can be expected, comparing the flexible model to a rigid-body model proves how motion and loads are affected by the flexibility of the structure. OrcaFlex and OpenFAST generally agree, but there are some differences in results due to different modeling approaches. Since no experimental data are available in the study, this paper only forms a baseline for future studies but still proves and describes the possibilities of the approach and codes.

17 WIND ENERGY↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2017

Thermodynamical effects and high resolution methods for compressible fluid flows

One of the fundamental differences of compressible fluid flows from incompressible fluid flows is the involvement of thermodynamics. This difference should be manifested in the design of numerical schemes. Unfortunately, the role of entropy, expressing irreversibility, is often neglected even though the entropy inequality, as a conceptual derivative, is verified for some first order schemes. In this paper, we refine the GRP solver to illustrate how the thermodynamical variation is integrated into the design of high resolution methods for compressible fluid flows and demonstrate numerically the importance of thermodynamic effects in the resolution of strong waves. As a by-product, we show that the GRP solver works for generic equations of state, and is independent of technical arguments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗
DOE OSTI2016

solveME: fast and reliable solution of nonlinear ME models

Background Genome-scale models of metabolism and macromolecular expression (ME) significantly expand the scope and predictive capabilities of constraint-based modeling. ME models present considerable computational challenges: they are much (>30 times) larger than corresponding metabolic reconstructions (M models), are multiscale, and growth maximization is a nonlinear programming (NLP) problem, mainly due to macromolecule dilution constraints. Results Here, we address these computational challenges. We develop a fast and numerically reliable solution method for growth maximization in ME models using a quad-precision NLP solver (Quad MINOS). Our method was up to 45 % faster than binary search for six significant digits in growth rate. We also develop a fast, quad-precision flux variability analysis that is accelerated (up to 60× speedup) via solver warm-starts. Finally, we employ the tools developed to investigate growth-coupled succinate overproduction, accounting for proteome constraints. Conclusions Just as genome-scale metabolic reconstructions have become an invaluable tool for computational and systems biologists, we anticipate that these fast and numerically reliable ME solution methods will accelerate the wide-spread adoption of ME models for researchers in these fields.

Biochemistry & Molecular Biology↗
NASA NTRS1990

ULTRA-SHARP nonoscillatory convection schemes for high-speed steady multidimensional flow

For convection-dominated flows, classical second-order methods are notoriously oscillatory and often unstable. For this reason, many computational fluid dynamicists have adopted various forms of (inherently stable) first-order upwinding over the past few decades. Although it is now well known that first-order convection schemes suffer from serious inaccuracies attributable to artificial viscosity or numerical diffusion under high convection conditions, these methods continue to enjoy widespread popularity for numerical heat transfer calculations, apparently due to a perceived lack of viable high accuracy alternatives. But alternatives are available. For example, nonoscillatory methods used in gasdynamics, including currently popular TVD schemes, can be easily adapted to multidimensional incompressible flow and convective transport. This, in itself, would be a major advance for numerical convective heat transfer, for example. But, as is shown, second-order TVD schemes form only a small, overly restrictive, subclass of a much more universal, and extremely simple, nonoscillatory flux-limiting strategy which can be applied to convection schemes of arbitrarily high order accuracy, while requiring only a simple tridiagonal ADI line-solver, as used in the majority of general purpose iterative codes for incompressible flow and numerical heat transfer. The new universal limiter and associated solution procedures form the so-called ULTRA-SHARP alternative for high resolution nonoscillatory multidimensional steady state high speed convective modelling.

Leonard, B. P.↗
DOE OSTI2021

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

New developments in the method of space-time conservation element and solution element: Applications to the Euler and Navier-Stokes equations

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods--i.e., finite difference, finite volume, finite element, and spectral methods. It is conceptually simple and designed to avoid several key limitations to the above traditional methods. An explicit model scheme for solving a simple 1-D unsteady convection-diffusion equation is constructed and used to illuminate major differences between the current method and those mentioned above. Unexpectedly, its amplification factors for the pure convection and pure diffusion cases are identical to those of the Leapfrog and the DuFort-Frankel schemes, respectively. Also, this explicit scheme and its Navier-Stokes extension have the unusual property that their stabilities are limited only by the CFL condition. Moreover, despite the fact that it does not use any flux-limiter or slope-limiter, the Navier-Stokes solver is capable of generating highly accurate shock tube solutions with shock discontinuities being resolved within one mesh interval. An accurate Euler solver also is constructed through another extension. It has many unusual properties, e.g., numerical diffusion at all mesh points can be controlled by a set of local parameters.

Chang, Sin-Chung↗
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗
DOE OSTI2021

Toward performance-portable PETSc for GPU-based exascale systems

The Portable Extensible Toolkit for Scientific computation (PETSc) library delivers scalable solvers for nonlinear time-dependent differential and algebraic equations and for numerical optimization. The PETSc design for performance portability addresses fundamental GPU accelerator challenges and stresses flexibility and extensibility by separating the programming model used by the application from that used by the library, and it enables application developers to use their preferred programming model, such as Kokkos, RAJA, SYCL, HIP, CUDA, or OpenCL, on upcoming exascale systems. Furthermore, a blueprint for using GPUs from PETSc-based codes is provided, and case studies emphasize the flexibility and high performance achieved on current GPU-based systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A high-frequency, low-power resonant radio-frequency neutron spin flipper for high-resolution spectroscopy

Here, we present a resonant-mode, transverse-field, radio-frequency (rf) neutron spin flipper design that uses high-temperature superconducting films to ensure sharp transitions between uniform magnetic field regions. Resonant mode allows for low-power, high-frequency operation but requires strict homogeneity of the magnetic fields inside the device. This design was found to efficiently flip neutrons at 96.6 ± 0.6% at an effective frequency of 4 MHz in bootstrap configuration with a beam size of 2.4 × 2.5 cm 2 and a wavelength of 0.4 nm. The high frequency and efficiency enable this device to perform high-resolution neutron spectroscopy with comparable performance with currently implemented rf flipper designs. The limitation of the maximum frequency was found due to the field homogeneity of the device. We numerically analyze the maximum possible efficiency of this design using a Bloch solver simulation with magnetic fields generated from finite-element simulations. We also discuss future improvements of the efficiency and frequency to the design based on the experimental and simulation results.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Butterfly Factorization Via Randomized Matrix-Vector Multiplications

This paper presents an adaptive randomized algorithm for computing the butterfly factorization of an m × n matrix with m ≈ n provided that both the matrix and its transpose can be rapidly applied to arbitrary vectors. The resulting factorization is composed of O(log n) sparse factors, each containing O(n) nonzero entries. The factorization can be attained using O(n 3/2 log n) computation and O(n log n) memory resources. Furthermore, the proposed algorithm can be implemented in parallel and can apply to matrices with strong or weak admissibility conditions arising from surface integral equation solvers as well as multi-frontal-based finite-difference, finite-element, or finite-volume solvers. A distributed-memory parallel implementation of the algorithm demonstrates excellent scaling behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

The General Atomics Center for Magnetohydrodynamic Modeling (CEMM) Final Report

The Center for Extended Magnetohydrodynamic Modeling (CEMM) is built around the further development, verification, validation, and other application of two major extended magnetohydrodynamics (MHD) codes NIMROD and M3D-C1. The General Atomics (GA) CEMM Program focuses on three main research areas: extension of the M3D-C1 MHD code plasma response calculations to new regimes, coupling of the DK4D drift kinetic equation solver to a temperature evolution equation, and simulations of disruption mitigation for unstable plasmas with pre-existing magnetic islands using massive gas injection (MGI) with the NIMROD MHD code. A systematic scan of the zero-crossing of the rotation profile in a DIII-D ITER-similar shape equilibrium using linear, time-independent modeling with M3D-C1 was performed. A new Python utility, autoC1, has been developed to automate the running of M3D-C1 for linear simulations. A number of ELM suppression experiments were analyzed including a validation and cross-code verification with the MARS-F MHD code of magnetic signal phase shifts across ELM-suppression bifurcations. A linearized, time-dependent temperature equation was tightly coupled to the electron and ion drift-kinetic equations (DKEs) in the DK4D drift-kinetic equation solver. Various explicit and implicit implementations were attempted, but all methods proved numerically unstable, even with the addition of numerical dissipation. Several simulations were completed investigating the effects of pre-existing islands on disruption mitigation by MGI in DIII-D using NIMROD.

3D MHD Simulation↗
DOE OSTI2019

Enhancement of EM Pump Performance through Modeling and Testing

Liquid sodium has seen development as a nuclear reactor coolant in the United States for over 50 years. While sodium has many desirable characteristics as a heat transfer fluid, its reactivity poses numerous engineering challenges for even the most basic plant components like pumps. Despite the challenges, mechanical pumps have seen success in sodium service albeit at the expense of simplicity and flexibility. In fact, significant failures of these pumps have been recorded which can be attributed to their complexity. Notwithstanding the complexity, mechanical pumps are still widely proposed for the main circuit pump of sodium cooled reactors. However, an alternative to mechanical pumps are Electromagnetic Induction Pumps (EMIPs) which use the principle of the Lorentz Force. EMIPs avoid mechanically moving parts in contact with sodium which simplifies the design by avoiding mechanical wear and seals. Therefore, these pumps are generally safer and more reliable than their mechanical counterparts. That being said, EMIPs suffer from relatively low efficiencies that range from 5% to 45%. Some of the inefficiencies are inherent to the EMIPs pumps; namely, inefficiencies arise in the form of resistive heating losses in the fluid itself. However, other inefficiencies arise from Edge Effects which include Finite Length and Finite Width effects. Due to the significant safety advantages EMIPs have compared to mechanical pumps, there is a strong motivation to integrate these technologies into the design of Generation IV liquid-metal-cooled nuclear power plants. However, EMIPs have low efficiencies partly attributed to Edge Effects. Much of the literature focusing on Edge Effects is dated or simply does not exist in the public domain. More recently, numerical models of EMIPs have shown promising results. However, these models are often custom solvers or use complex software suits which limit their application. This report describes the work completed to understand Edge Effects in two EMIPs designs. These designs are the Annular Linear Induction Pump (ALIP) and Permanent Magnet Induction Pump (PMIP) also known as a Moving Magnet Pump (MMP). The methodology developed during this project aim at measuring key performance parameters (pressure-flowrate and efficiency-slip performance) using high temperature liquid-sodium. Baseline experimental measurements were then used to verify the accuracy of numerical models developed in the software suit ANSYS/FLUENT. On the experimental side, instrumentation was developed to characterize the pressure-flowrate and efficiency-slip performance of both the ALIP and PMIP in high-temperature liquid-sodium. Additionally, instrumentation was developed to measure potential output pressure pulsations and the magnetic field distributions. Moreover, a specialized ALIP was constructed to study the impacts magnetic field tapering, also known as coil grading, on the Finite Length Effect. Lastly, experimental performance curves of the PMIP were collected at 300 oC and 600 oC and compared to analytical theory which accounted for Finite Width Effects. These results showed that the Finite Width Effect and non-uniform magnetic field in a PMIP impact the performance. On the numerical side, models were developed using the software suit FLUENT/ANSYS. These methods were initially compared against a simple MHD flow in a Permanent Magnet Electromagnet Flowmeter (PMEMFM). The results were significant in that they showed agreeance between the numerical model and demonstrated that PMEMFM’s could be calibrated numerically. These modeling methods were then applied to a prototypic PMIP. Results showed that this simple model provided reasonably accurate results of the PMIP p-Q performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

GentenMPI: Distributed Memory Sparse Tensor Decomposition

GentenMPl is a toolkit of sparse canonical polyadic (CP) tensor decomposition algorithms that is designed to run effectively on distributed-memory high-performance computers. Its use of distributed-memory parallelism enables it to efficiently decompose tensors that are too large for a single compute node's memory. GentenMPl leverages Sandia's decades-long investment in the Trilinos solver framework for much of its parallel-computation capability. Trilinos contains numerical algorithms and linear algebra classes that have been optimized for parallel simulation of complex physical phenomena. This work applies these tools to the data science problem of sparse tensor decomposition. In this report, we describe the use of Trilinos in GentenMPl, extensions needed for sparse tensor decomposition, and implementations of the CP-ALS (CP via alternating least squares) and GCP-SGD (generalized CP via stochastic gradient descent) sparse tensor decomposition algorithms. We show that GentenMPl can decompose sparse tensors of extreme size, e.g., a 12.6-terabyte tensor on 8192 computer cores. We demonstrate that the Trilinos backbone provides good strong and weak scaling of the tensor decomposition algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Computationally Efficient Control Co-Design Optimization Framework with Mixed-Fidelity Fluid and Structure Analysis: OpenTurbineCoDe (OTCD) - ARPA-E ATLANTIS Phase I Project (Final Technical Report)

The goal of the Department of Energy (DOE) Advanced Research Projects Agency-Energy (ARPA-E) Aerodynamic Turbines Lighter and Afloat with Nautical Technologies and Integrated Servo-control (ATLANTIS) Program is to develop new technologies for floating offshore wind turbines, or wind farms, using the discipline of control co-design (CCD.) In this context, our goal is to develop a computationally efficient optimization framework for design of floating offshore wind turbines. Our specific aim is to utilize high-fidelity structural, aerodynamic, aero-structural tools, and to derive control-oriented reduced- or low-order models directly from the high-fidelity tools. We proposed a mixed-fidelity modeling approach which means that we are also using low- and mid-fidelity tools when necessary. This research is conducted by a multidisciplinary team consisting of Rutgers University, University of Michigan, Brigham Young University, and the National Renewable Energy Laboratory (NREL). The computational framework, called OpenTurbineCoDe, is designed to integrate, where possible, traditional structural, aerodynamic, aeroelastic models (e.g., OpenFAST) and advanced control algorithms with higher fidelity simulation tools including Reynolds-averaged Navier–Stokes (RANS) solvers, and three-dimensional structural finite element solvers. All the high-fidelity tools used in this research provide numerically exact gradients to facilitate both efficient optimization and local linearization for control implementation.

17 WIND ENERGY↗
NASA NTRS1987

Flow solution on a dual-block grid around an airplane

The compressible flow around a complex fighter-aircraft configuration (fuselage, cranked delta wing, canard, and inlet) is simulated numerically using a novel grid scheme and a finite-volume Euler solver. The patched dual-block grid is generated by an algebraic procedure based on transfinite interpolation, and the explicit Runge-Kutta time-stepping Euler solver is implemented with a high degree of vectorization on a Cyber 205 processor. Results are presented in extensive graphs and diagrams and characterized in detail. The concentration of grid points near the wing apex in the present scheme is shown to facilitate capture of the vortex generated by the leading edge at high angles of attack and modeling of its interaction with the canard wake.

Eriksson, Lars-Erik↗