SEARCH · Search NASA

Results for “splitting methods”

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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1994

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗

NASA NTRS1993

Incompressible spectral-element method: Derivation of equations

A fractional-step splitting scheme breaks the full Navier-Stokes equations into explicit and implicit portions amenable to the calculus of variations. Beginning with the functional forms of the Poisson and Helmholtz equations, we substitute finite expansion series for the dependent variables and derive the matrix equations for the unknown expansion coefficients. This method employs a new splitting scheme which differs from conventional three-step (nonlinear, pressure, viscous) schemes. The nonlinear step appears in the conventional, explicit manner, the difference occurs in the pressure step. Instead of solving for the pressure gradient using the nonlinear velocity, we add the viscous portion of the Navier-Stokes equation from the previous time step to the velocity before solving for the pressure gradient. By combining this 'predicted' pressure gradient with the nonlinear velocity in an explicit term, and the Crank-Nicholson method for the viscous terms, we develop a Helmholtz equation for the final velocity.

Deanna, Russell G.↗

NASA NTRS1991

Pressure flux-split viscous solutions for subsonic diffusers

A method for predicting flow in diffusers with inlet swirl has been developed. Solutions of the pressure flux-split Reduced Navier Stokes (RNS) equations are obtained for flows in axisymmetric diffusers. Viscous-Inviscid interactions and flow fields with toroidal recirculation regions are efficiently captured. The computational model is verified by comparision with experimental data and other computations. The extreme sensitivity to grid, turbulent closure model and inlet profiles is discussed.

Pordal, H. S.↗

DOE OSTI2025

Holography for BCFTs with multiple boundaries: multi-splitting quenches

We elaborate on the method introduced in arXiv:2403.02165 for holographic duals of Boundary Conformal Field Theories (BCFTs) with multiple boundaries. Naïvely, as the number of boundaries grow large such a calculation becomes unmanageable. We show that is not always the case. Using these advances we calculate the entanglement entropy as a function of time for 1+1-dimensional CFTs that are split into N subsystems. We give explicit results for N = 4 and N = 17. We find that all qualitative differences that arise for larger N are present N = 4.

AdS-CFT correspondence↗

DOE OSTI2025

Informed total-error-minimizing priors: Interpretable cosmological parameter constraints despite complex nuisance effects

While Bayesian inference techniques are standard in cosmological analyses, it is common to interpret resulting parameter constraints with a frequentist intuition. This intuition can fail, for example, when marginalizing high-dimensional parameter spaces onto subsets of parameters, because of what has come to be known as projection effects or prior volume effects. We present the method of informed total-error-minimizing (ITEM) priors to address this problem. An ITEM prior is a prior distribution on a set of nuisance parameters, such as those describing astrophysical or calibration systematics, intended to enforce the validity of a frequentist interpretation of the posterior constraints derived for a set of target parameters (e.g., cosmological parameters). Our method works as follows. For a set of plausible nuisance realizations, we generate target parameter posteriors using several different candidate priors for the nuisance parameters. We reject candidate priors that do not accomplish the minimum requirements of bias (of point estimates) and coverage (of confidence regions among a set of noisy realizations of the data) for the target parameters on one or more of the plausible nuisance realizations. Of the priors that survive this cut, we select the ITEM prior as the one that minimizes the total error of the marginalized posteriors of the target parameters. As a proof of concept, we applied our method to the density split statistics measured in Dark Energy Survey Year 1 data. We demonstrate that the ITEM priors substantially reduce prior volume effects that otherwise arise and that they allow for sharpened yet robust constraints on the parameters of interest.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1981

Flux vector splitting of the inviscid gasdynamic equations with application to finite-difference methods

The conservation-law form of the inviscid gasdynamic equations has the remarkable property that the nonlinear flux vectors are homogeneous functions of degree one. This property readily permits the splitting of flux vectors into subvectors by similarity transformations so that each subvector has associated with it a specified eigenvalue spectrum. As a consequence of flux vector splitting, new explicit and implicit dissipative finite-difference schemes are developed for first-order hyperbolic systems of equations. Appropriate one-sided spatial differences for each split flux vector are used throughout the computational field even if the flow is locally subsonic. The results of some preliminary numerical computations are included.

Steger, J. L.↗

NASA NTRS1989

A class of high resolution explicit and implicit shock-capturing methods

An attempt is made to give a unified and generalized formulation of a class of high resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock wave computations. Included is a systematic review of the basic design principle of the various related numerical methods. Special emphasis is on the construction of the basis nonlinear, spatially second and third order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and the flux vector splitting approaches. Generalization of these methods to efficiently include equilibrium real gases and large systems of nonequilibrium flows are discussed. Some issues concerning the applicability of these methods that were designed for homogeneous hyperbolic conservation laws to problems containing stiff source terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for 1-, 2- and 3-dimensional gas dynamics problems.

Yee, H. C.↗

DOE OSTI2021

Preconditioners for multiphase poromechanics with strong capillarity

This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two-phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection-dominated to diffusion-dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.

42 ENGINEERING↗

DOE OSTI2020

Cobalt Metal–Cobalt Carbide Composite Microspheres for Water Reduction Electrocatalysis

Microspheres of cobalt metal–cobalt carbide (Co–Co x C, Co x C: Co 2 C and Co 3 C) composite with carbon shells were prepared via an OH – - and Cl – -assisted polyol method and explored for electrocatalytic activity and stability for the hydrogen evolution reaction (HER) in acidic media. From our transmission electron microscopy observations, the outermost surfaces of the as-prepared Co–Co x C composites were primarily covered with Co 2 C crystallites. Our best performing electrocatalyst exhibited superior HER activity with an overpotential of 78 mV to reach a current density of -10 mA·cm –2 , a Tafel slope of 87.8 mV·dec –1 , and 1 h of electrode durability. We show that this excellent HER performance is primarily due to the superior intrinsic activity of Co 2 C, as well as the high electrical conductivity resulting from the inclusion of cobalt metal and the presence of graphitic carbon shells in and on the composite, respectively. Using both computational and experimental approaches, we determine here that the carbon-rich cobalt carbide (Co 2 C) phase is more favorable for the HER than the carbon-poor phase (Co 3 C).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Analysis of the classical trajectory treatment of photon dynamics for polaritonic phenomena

Simulating photon dynamics in strong light–matter coupling situations via classical trajectories is proving to be powerful and practical. Here, in this work, we analyze the performance of the approach through the lens of the exact factorization approach. Since the exact factorization enables a rigorous definition of the potentials driving the photonic motion, it allows us to identify that the underestimation of photon number and intensities observed in earlier work is primarily due to an inadequate accounting of light–matter correlation in the classical Ehrenfest force rather than errors from treating the photons quasiclassically per se. The latter becomes problematic when the number of photons per mode begins to exceed a half.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Evaluation of Seismic Artificial Intelligence with Uncertainty

Artificial intelligence has transformed the seismic community with deep learning models (DLMs) that are trained to complete specific tasks within workflows. However, there is still a lack of robust evaluation frameworks for evaluating and comparing DLMs. Here, we address this gap by designing an evaluation framework that jointly incorporates two crucial aspects: performance uncertainty and learning efficiency. To target these aspects, we meticulously construct the training, validation, and test splits using a clustering method tailored to seismic data and enact an expansive training design to segregate performance uncertainty arising from stochastic training processes and random data sampling. The framework’s ability to guard against misleading declarations of model superiority is demonstrated through the evaluation of PhaseNet (Zhu and Beroza, 2018), a popular seismic phase picking DLM, under three training approaches. Our framework helps practitioners choose the best model for their problem and set performance expectations by explicitly analyzing model performance with uncertainty at varying budgets of training data.

58 GEOSCIENCES↗

NASA NTRS1974

Rotation in vibration, optimization, and aeroelastic stability problems

The effects of rotation in the areas of vibrations, dynamic stability, optimization, and aeroelasticity were studied. The governing equations of motion for the study of vibration and dynamic stability of a rapidly rotating deformable body were developed starting from the nonlinear theory of elasticity. Some common features such as the limitations of the classical theory of elasticity, the choice of axis system, the property of self-adjointness, the phenomenon of frequency splitting, shortcomings of stability methods as applied to gyroscopic systems, and the effect of internal and external damping on stability in gyroscopic systems are identified and discussed, and are then applied to three specific problems.

Kaza, K. R. V.↗

NASA NTRS1981

Analysis of turbulent flow about an isolated airfoil using a time dependent Navier-Stokes procedure

The procedure solves the Navier-Stokes equations by the consistently split linearized block implicit method of Briley and McDonald in a body fitted coordinate system. The procedure is described and results are presented for flow about an airfoil whose incidence changes from 6 degrees to 19 degrees at a Reynolds number of one million and Mach number of 0.2. In addition, the unsteady flow about an airfoil held at a constant 19 degree incidence is examined and compared to data.

Shamroth, S. J.↗

NASA NTRS1983

Navier-Stokes solution on the CYBER-203 by a pseudospectral technique

A three-level, time-split, mixed spectral/finite difference method for the numerical solution of the three-dimensional, compressible Navier-Stokes equations has been developed and implemented on the Control Data Corporation (CDC) CYBER-203. This method uses a spectral representation for the flow variables in the streamwise and spanwise coordinates, and central differences in the normal direction. The five dependent variables are interleaved one horizontal plane at a time and the array of their values at the grid points of each horizontal plane is a typical vector in the computation. The code is organized so as to require, per time step, a single forward-backward pass through the entire data base. The one-and two-dimensional Fast Fourier Transforms are performed using software especially developed for the CYBER-203.

Lambiotte, J. J.↗

NASA NTRS1984

Numerical calculations of complex Mach reflection

Numerical simulations of the interaction of a planar blast wave with a compression ramp are presented. The split coefficient matrix (SCM) method in conjunction with boundary shock and floating discontinuity-fitting procedures was employed to obtain the time-asymptotic solutions of the two-dimensional, unsteady Euler equations. The solutions were computed for the complex Mach reflection (CMR) regime of the shock diffraction problem in an attempt to explore the basic physical process governing the evolution of an incipient second Mach stem and the associated topological changes. Numerical results were obtained for shock diffraction over a 40 degree ramp with varying incident shock Mach numbers. The validity of the present approach has been substantiated by experimental observations and earlier numerical calculations.

Yamamoto, O.↗

NASA NTRS1985

Water quality monitor for recovered spacecraft water

A total organic carbon (TOC) analysis system based on ultraviolet absorption is described. The equation for measuring the intensity of the absorbed radiation of the organic substances, which is based on the Lambert-Beer law, is given; the intensity of the absorption is proportional to the concentration of the solution. The operation of the UV-Absorption analyzer, which utilizes a split beam, two wvaelength method, is studied. The influences of the cell path length and specific compounds in the solution flowing through the cell on absorbances is discussed. The performance and response of the analyzer is evaluated; good correlation is observed between the absorption value and TOC. The advantage of the UV-Absorption as compared with the UV-Oxidation are examined.

Ejzak, E. M.↗

NASA NTRS1987

A numerical simulation of finite-length Taylor-Couette flow

The processes leading to laminar-turbulent transition in finite-channel-length Taylor-Couette flow are investigated analytically, solving the unsteady incompressible Navier-Stokes equations by spectral-collocation methods. A time-split algorithm, implementable in both axisymmetric and fully three-dimensional time-accurate versions, and an algorithm based on the staggered-mesh discretization of Bernardi and Maday (1986) are described in detail, and results obtained by applying the axisymmetric version of the first algorithm and a steady-state version of the second are presented graphically and compared with published experimental data. The feasibility of full three-dimensional simulations of the progression through chaotic states to turbulence under the constraints of Taylor-Couette flow is demonstrated.

Streett, C. L.↗

NASA NTRS1988

Modulated-splitting-ratio fiber-optic temperature sensor

A fiber-optic temperature sensor is described, which uses a small silicon beamsplitter whose splitting ratio varies as a function of temperature. A four-beam technique is used to measure the sensor's temperature-indicating splitting ratio. This referencing method provides a measurement that is largely independent of the transmission properties of the sensor's optical fiber link. A significant advantage of this sensor, relative to other fiber-optic sensors, is its high stability, which permits the fiber-optic components to be readily substituted, thereby simplifying the sensor's installation and maintenance.

Beheim, Glenn↗