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 613 records · Page 34

Neural Network Solver for Small Quantum Clusters

Machine learning approaches have recently been applied to the study of various problems in physics. Most of these studies are focused on interpreting the data generated by conventional numerical methods or the data on an existing experimental database. An interesting question is whether it is possible to use a machine learning approach, in particular a neural network, for solving the many-body problem. In this paper, we present a neural network solver for the single impurity Anderson model, the paradigm of an interacting quantum problem in small clusters. We demonstrate that the neural-network-based solver provides quantitative accurate results for the spectral function as compared to the exact diagonalization method. This opens the possibility of utilizing the neural network approach as an impurity solver for other many-body numerical approaches, such as the dynamical mean field theory.

36 MATERIALS SCIENCE↗

Focused wave interactions with floating structures: a blind comparative study

The paper presents results from the Collaborative Computational Project in Wave Structure Interaction (CCP-WSI) Blind Test Series 2. Without prior access to the physical data, participants, with numerical methods ranging from low-fidelity linear models to fully non-linear Navier-Stokes (NS) solvers, simulate the interaction between focused wave events and two separate, taut-moored, floating structures: a hemispherical-bottomed cylinder and a cylinder with a moonpool. The ‘blind’ numerical predictions for heave, surge, pitch and mooring load, are compared against physical measurements. Dynamic time warping is used to quantify the predictive capability of participating methods. In general, NS solvers and hybrid methods give more accurate predictions; however, heave amplitude is predicted reasonably well by all methods; and a WEC-Sim implementation, with CFD-informed viscous terms, demonstrates comparable predictive capability to even the stronger NS solvers. Large variations in the solutions are observed (even among similar methods), highlighting a need for standardisation in the numerical modelling of WSI problems.

42 ENGINEERING↗

On the Direct Construction of the Steady Traveling Solution to High Explosive Sandwich, Cylinder and Aquarium Tests via a Streamline Finite Volume Approximation

The cylinder, sandwich and aquarium tests are routine experiments used to investigate the performance of high explosives. Here, we present a solution technique which simulates these experiments as two dimensional problems in steady traveling coordinates. Using a coordinate system aligned with the streamlines of the flow, an efficient numerical scheme is developed which uses an approximate Riemann solver to model the change in state across streamlines. The technique can also accommodate arbitrary equations of state for the materials. Two verification exercises show that the algorithm converges at second order. Validation studies using the aquarium experiment agreed well with simulations, though the omission of a material strength model leads to some discrepancies between simulation and experiments for the cylinder test. This computational technique is able to predict the behavior of these common high explosive experiments on the order of minutes on a single processor with very fine spatial resolution, on the order of 10 μm.

36 MATERIALS SCIENCE↗

Numerical Simulation and Experimental Study of the Gas–Solid Flow Behavior Inside a Full-Loop Circulating Fluidized Bed: Evaluation of Different Drag Models

Both experimental and computational studies of the fluidization of high-density polyethylene (HDPE) particles in a small-scale full-loop circulating fluidized bed (CFB) are conducted. Experimental measurements of pressure drop are taken at various locations along the bed. The solids circulation rate is measured with an advanced particle image velocimetry (PIV) technique. Bed height of the quasi-static region in the standpipe is also measured. Comparative numerical simulations are performed with a computational fluid dynamics solver utilizing a discrete element method (CFD-DEM). Here, this paper examines the effect of different drag laws used in the CFD simulations through a detailed and direct comparison with experimental data from a small-scale, full-loop circulating fluidized bed. The Hill–Koch–Ladd drag correlation was shown to have good agreement with respect to system component pressure drop and inventory height in the standpipe.

42 ENGINEERING↗

Application of mesh refinement to relativistic magnetic reconnection

During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Orbital differentiation in Hund metals

Orbital differentiation is a common theme in multiorbital systems, yet a complete understanding of it is still missing. Here, we consider a minimal model for orbital differentiation in Hund metals with a highly accurate method: We use the numerical renormalization group as a real-frequency impurity solver for a dynamical mean-field study of three-orbital Hubbard models, where a crystal field shifts one orbital in energy. The individual phases are characterized with dynamic correlation functions and their relation to diverse Kondo temperatures. Upon approaching the orbital-selective Mott transition, we find a strongly suppressed spin coherence scale and uncover the emergence of a singular Fermi liquid and interband doublon-holon excitations. Our theory describes the diverse polarization-driven phenomena in the $t$ 2 g bands of materials such as ruthenates and iron-based superconductors, and our methodological advances pave the way toward real-frequency analyses of strongly correlated materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scientific Core Library Stack (SCLS) v2026

SCLS (Scientific Core Library Stack) is an opinionated build and packaging system for scientific computing libraries developed at Lawrence Berkeley National Laboratory. It produces a coherent, reproducible stack of numerical libraries — including BLAS/LAPACK, MPI, sparse direct and iterative solvers, graph partitioners, and parallel I/O libraries (e.g., PETSc, SLEPc, HDF5, NetCDF, MUMPS, OpenBLAS) — that work together without manual repair by downstream scientific software. From a single recipe-and-flavor model, SCLS produces native RPM packages for RHEL-family Linux, DEB packages for Debian/Ubuntu, direct Unix-style prefix installs for HPC and locked-down environments, and native macOS builds. Multiple build "flavors" (e.g., GCC+OpenBLAS, GCC+MKL, Intel+MKL, debug) coexist in distinct prefixes on the same host. Compared to general-purpose meta-build frameworks, SCLS is deliberately curated rather than infinitely configurable. It enforces deterministic, audit-friendly behavior: explicit build dependencies, no silent feature autodetection, a clear open-source license policy, and rpath-based runtime linkage so installs integrate cleanly with standard package-manager workflows.

Messe, Christian [Lawrence Berkeley National Labor↗

CUDO: closed-form universal dwell-time optimization for computer-controlled optical surfacing

Precision optical figuring demands fast and accurate dwell time optimization to reach nanometer- and sub-nanometer-level accuracy in next-generation optical systems. We introduce CUDO (closed-form universal dwell-time optimization), the first, to the best of our knowledge, unified closed-form analytical framework that supports both function-form and matrix-form dwell time models in computer-controlled optical surfacing (CCOS). In contrast to traditional methods, which rely on iterative optimization and hyperparameter tuning, our framework derives direct analytical solutions with no adjustable parameters. This approach unifies the solution principles of existing methods within a single mathematical model, delivering three key advantages: (1) accuracy on par with, or superior to, iterative solvers, (2) substantial reduction in computation time, and (3) numerical robustness. Comparative studies with prior art confirm that closed-form solutions achieve equivalent residual error while removing runtime bottlenecks. By simplifying the implementation and enabling real-time, scalable deployment, CUDO establishes a practical foundation for future deterministic fabrication of large-aperture and high-performance optics.

36 MATERIALS SCIENCE↗

Status Report on NEAMS System Analysis Module Development

Under the Reactor Product Line (RPL) of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an advanced SFR System Analysis Module (SAM) is being developed at Argonne National Laboratory. The goal of the SAM development is to provide fast-running, improved-fidelity, whole-plant transient analyses capabilities. SAM utilizes an object-oriented application framework MOOSE), and its underlying meshing and finite-element library libMesh, as well as linear and non-linear solvers PETSc, to leverage modern advanced software environments and numerical methods. It also incorporates advances in physical and empirical models and seeks closure models based on information from high-fidelity simulations and experiments. This report provides an update on the SAM development, and summarizes the activities performed in FY15 and the first quarter of FY16. The tasks include: (1) implement the support of 2nd-order finite elements in SAM components for improved accuracy and computational efficiency; (2) improve the conjugate heat transfer modeling and develop pseudo 3-D full-core reactor heat transfer capabilities; (3) perform verification and validation tests as well as demonstration simulations; (4) develop the coupling requirements for SAS4A/SASSYS-1 and SAM integration.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Numerical simulation of a full-loop circulating fluidized bed under different operating conditions

Both experimental and computational studies of the fluidization of high-density polyethylene (HDPE) particles in a small-scale full-loop circulating fluidized bed are conducted. Experimental measurements of pressure drop are taken at different locations along the bed. The solids circulation rate is measured with an advanced Particle Image Velocimetry (PIV) technique. The bed height of the quasi-static region in the standpipe is also measured. Comparative numerical simulations are performed with a Computational Fluid Dynamics solver utilizing a Discrete Element Method (CFD-DEM). This paper reports a detailed and direct comparison between CFD-DEM results and experimental data for realistic gas-solid fluidization in a full-loop circulating fluidized bed system. The comparison reveals good agreement with respect to system component pressure drop and inventory height in the standpipe. In addition, the effect of different drag laws applied within the CFD simulation is examined and compared with experimental results.

01 COAL, LIGNITE, AND PEAT↗

SAM Plug-in Development (Phase I Final Report)

The DOE Office of Nuclear Energy (NE) has created an extensive set of advanced modeling and simulation tools for nuclear engineering analysis. The advanced capabilities of these newer analysis codes require more in-depth training, skills, and knowledge in order to effectively utilize them for the design, analysis, and licensing of advanced nuclear systems and experiments. A high learning curve for inexperienced users may deter organizations from incorporating these tools into their internal processes. This project involved development of a plug-in to the Symbolic Nuclear Analysis Package (SNAP) for the System Analysis Module (SAM) tool. SAM is an advanced system analysis tool for reactor transient analyses being developed at Argonne National Laboratory under the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. SAM utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage modern advanced software environments and numerical methods. SNAP provides a highly flexible framework for creating, modifying and documenting input for engineering analysis codes such as SAM as well as extensive functionality for submitting, monitoring, and interacting with the codes through an intuitive graphical user interface (GUI). The common user interface provided by SNAP minimizes the learning curve for engineers starting with a new analysis code and provides an intuitive framework for transitioning between different analysis codes. SNAP provides a powerful but intuitive interface to facilitate access to advanced modeling and simulation tools for inexperienced users. Unlike many “form based” GUI’s, SNAP maps each engineering code’s component input to an internal database which manages all component input parameters along with component interconnections. This level of abstraction permits SNAP to support several advanced capabilities such as renodalization, model validation and consistency checks, embedded documentation, model notebook generation, data ownership and reviewer tracking, and variable assignment for inputs to name a few. SNAP includes a built-in Python interpreter and is interfaced to several commercial and open source packages including CPython, MATLAB/OCTAVE, Microsoft Office, Open Office, and SANDIA’s DAKOTA package which provides Uncertainty Quantification analysis through the SNAP plug-ins. Phase I of this project involved development a fully functional basic SAM plug-in to SNAP. This plug-in provides the ability to import existing models, graphically construct, edit and submit models using SNAP’s extensive functionality.

99 GENERAL AND MISCELLANEOUS↗

Exponential Time Differencing Gauge Method for Incompressible Viscous Flows

Abstract In this paper, we study an exponential time differencing method for solving the gauge system of incompressible viscous flows governed by Stokes or Navier-Stokes equations. The momentum equation is decoupled from the kinematic equation at a discrete level and is then solved by exponential time stepping multistep schemes in our approach. We analyze the stability of the proposed method and rigorously prove that the first order exponential time differencing scheme is unconditionally stable for the Stokes problem. We also present a compact representation of the algorithm for problems on rectangular domains, which makes FFT-based solvers available for the resulting fully discretized system. Various numerical experiments in two and three dimensional spaces are carried out to demonstrate the accuracy and stability of the proposed method.

Physics↗

Improved 3-D turbomachinery CFD algorithm

The building blocks of a computer algorithm developed for the time-accurate flow analysis of rotating machines are described. The flow model is a finite volume method utilizing a high resolution approximate Riemann solver for interface flux definitions. This block LU implicit numerical scheme possesses apparent unconditional stability. Multi-block composite gridding is used to orderly partition the field into a specified arrangement. Block interfaces, including dynamic interfaces, are treated such as to mimic interior block communication. Special attention is given to the reduction of in-core memory requirements by placing the burden on secondary storage media. Broad applicability is implied, although the results presented are restricted to that of an even blade count configuration. Several other configurations are presently under investigation, the results of which will appear in subsequent publications.

Janus, J. Mark↗

Navier-Stokes computation of compressible turbulent flows with a second order closure

The objective was the development of a complete second order closure for wall bounded flows, including all components of the dissipation rate tensor and a numerical solution procedure for the resulting system of equations. The main topics discussed are the closure of the pressure correlations and the viscous destruction terms in the dissipation rate equations and the numerical solution scheme based on a block-tridiagonal solver for the nine equations required for the prediction of plane or axisymmetric flows.

Dingus, C.↗

Non-linear k-epsilon-v(sup 2)(bar) modeling with application to high-lift

The k-epsilon-v(sup 2)(bar) model has been investigated to quantify its predictive performance on two high-lift configurations: 2D flow over a single-element aerofoil, involving closed-type separation; 3D flow over a prolate spheroid, involving open-type separation. A 'code-friendly' modification has been proposed which enhances the numerical stability, in particular, for explicit and uncoupled flow solvers. As a result of introducing Reynolds-number dependence into a coefficient of the s-equation, the skin-friction distribution for the by-pass transitional flow over a flat plate is better predicted. In order to improve deficiencies arising from the Boussinesq approximation, a nonlinear stress-strain constitutive relation was adopted, in which the only one free constant is calibrated on the basis of DNS data, and the Reynolds-stress anisotropy near the wall is fairly well represented.

Lien, F. S.↗

Navier-Stokes Flowfield Simulation of Boeing 747-200 as Platform for SOFIA

Steady and unsteady viscous, three-dimensional flowfields are calculated using a thin layer approximation of Navier-Stokes equations in conjunction with Chimera overset grids. The finite-difference numerical scheme uses structured grids and a pentadiagonal flow solver called "OVERFLOW". The configuration of Boeing 747-200 has been chosen as one of configurations to be used as a platform for the SOFIA (Stratospheric Observatory For Infrared Astronomy). Initially, the steady flowfield of the full aircraft is calculated for the clean configuration (without a cavity to house telescope). This solution is then used to start the unsteady flowfield of a configuration containing cavity housing the observation telescope and its peripheral units. Analysis of unsteady flowfield in the cavity and its influence on the tail empennage, as well as the noise due to turbulence and optical quality of the flow are the main focus of this study. For the configuration considered here, the telescope housing cavity is located slightly downstream of the portwing. The entire flow-field is carefully constructed using 45 overset grids and consists of nearly 4 million grid points. All the computations axe done at one freestream flow condition of M(sub infinity) = 0.85, alpha = 2.5deg, and a Reynolds of Re = 1.85x10deg

Srinivasan, G.R.↗

Near-Field Noise Computation for a Supersonic Circular Jet

A fully expanded, high-Reynolds-number, supersonic circular jet of Mach number 1.4 is simulated, using a 3-D finite-volume Navier-Stokes solver, with emphasis on the near field noise. The numerical results are generally in good agreement with existing experimental findings.

Loh, Ching Y.↗

Newly-Developed 3D GRMHD Code and its Application to Jet Formation

We have developed a new three-dimensional general relativistic magnetohydrodynamic code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-interpolated constrained transport scheme is used to maintain a divergence-free magnetic field. We have performed various 1-dimensional test problems in both special and general relativity by using several reconstruction methods and found that the new 3D GRMHD code shows substantial improvements over our previous model. The . preliminary results show the jet formations from a geometrically thin accretion disk near a non-rotating and a rotating black hole. We will discuss the jet properties depended on the rotation of a black hole and the magnetic field strength.

Mizuno, Y.↗