SEARCH · Engineering Papers

Results for “Finite Volume”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1982

A vectorized, finite-volume, adaptive grid algorithm applied to planetary entry problems

An adaptive grid, finite-volume method has been applied to problems in planetary entry for computing complete flowfields. The adaption algorithm is implicit in nature and is keyed to resolve user specified gradients. The finite-volume algorithm is explicit, utilizing a maximum time step advancement at each grid point to accelerate convergence to the steady state. The present version of the code is for the laminar flow of a perfect gas. The role of the adaption algorithm in resolving various features of blunt body/wake flow for planetary entry conditions is emphasized.

Gnoffo, P. A.↗

DOE OSTI2023

Three relativistic neutrons in a finite volume

We generalize the relativistic field-theoretic (RFT) three-particle finite-volume formalism to systems of three identical, massive, spin-1/2 fermions, such as three neutrons. This allows, in principle, for the determination of the three-neutron interaction from the finite-volume spectrum of three-neutron states, which can be obtained from lattice QCD calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1983

A time-split finite-volume algorithm for three-dimensional flow-field simulation

A general finite-volume algorithm is developed for solving three-dimensional, time-dependent, compressible Navier-Stokes equations for high Reynolds number flows over an arbitrary geometry. This algorithm adapts MacCormack's (1982) explicit-implicit scheme to a time-split, three-dimensional finite-volume concept in a general coordinate system. It is shown that the thin-layer approximation in all three spatial directions significantly reduces the evaluation of viscous terms and allows the algorithm to solve more complicated geometries with all boundaries in two or all three directions. The calculated results using this method are found to be in good agreement with the experimental measurements of a blunt-fin induced shock wave and boundary-layer interaction problems. Observations of the existence of peak pressure, primary horseshoe and secondary vortices, and reversed supersonic zones show that computational fluid dynamics can effectively supplement the wind tunnel tests for aerodynamic design as well as for understanding basic fluid dynamics.

Hung, C. M.↗

NASA NTRS1986

Finite volume solution of the compressible boundary-layer equations

A box-type finite volume discretization is applied to the integral form of the compressible boundary layer equations. Boundary layer scaling is introduced through the grid construction: streamwise grid lines follow eta = y/h = const., where y is the normal coordinate and h(x) is a scale factor proportional to the boundary layer thickness. With this grid, similarity can be applied explicity to calculate initial conditions. The finite volume method preserves the physical transparency of the integral equations in the discrete approximation. The resulting scheme is accurate, efficient, and conceptually simple. Computations for similar and non-similar flows show excellent agreement with tabulated results, solutions computed with Keller's Box scheme, and experimental data.

Loyd, B.↗

DOE OSTI2024

Complex scaling in finite volume

Quantum resonances, i.e., metastable states with a finite lifetime, play an important role in nuclear physics and other domains. Describing this phenomenon theoretically is generally a challenging task. In this work, we combine two established techniques to address this challenge. Complex scaling makes it possible to calculate resonances with bound-state-like methods. Finite-volume simulations exploit the fact that the infinite-volume properties of quantum systems are encoded in how discrete energy levels change as one varies the size of the volume. Herein we apply complex scaling to systems in finite periodic boxes and derive the volume dependence of states in this scenario, demonstrating with explicit examples how one can use these relations to infer infinite-volume resonance energies and lifetimes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1986

Development of a finite volume time marching method

Progress made in the development of a time-marching finite volume method from September 1984 to December 1984 is reported. The objective of the work is to develop and demonstrate a Navier-Stokes approach for transonic flow which includes viscous terms in the finite volume method. The accuracy of the computational method will be verified using a transonic diffuser as a test case. The computational goal is to calculate the flow in sufficient detail and with sufficient accuracy that the loss generating mechanisms can be studied to assess the sources of inefficiency in the transonic diffuser.

Nicholson, Stephen↗

DOE OSTI2021

Deployment of the Finite Volume Method in Pronghorn for Gas and Salt cooled Pebble Bed Reactors

This report summarizes the activities related to ”Complete FHR and HTGR pebble bed simulator, including initial validation” funded by the NEAMS thermal-hydraulics focus area. The activity revolves around the coarse-mesh thermal-hydraulics code Pronghorn and its application to gas and salt cooled Pebble bed reactor (PBR). The main difference between gas and salt cooled PBR from a thermal-hydraulics perspective is the fluid. To address the difference in fluid behavior, two separate approaches are implemented in the MOOSE Navier Stokes module: 1) a Boussinesq approximation and 2) a fully compressible formulation. The developed finite volume method capabilities are used for improving pre-existing gas-cooled and salt-cooled pebble-bed reactor models. A steady-state, multiphysics gas-cooled pebble-bed reactor model is created that couples the equilibrium core depletion capability developed in previous work, and the finite volume method capability developed for this report. The salt-cooled pebble-bed reactor model is upgraded to use the incompressible finite volume method capability and then extended to three spatial dimensions. Finally, several verification-and-validation exercises performed with Pronghorn are documented using the verification-and-validation report of the MooseDoc system. The goal of this effort to document the verification-and-validation level of Pronghorn and improve stakeholder confidence in the results obtained with Pronghorn.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗

NASA NTRS2006

Finite Volume Numerical Methods for Aeroheating Rate Calculations from Infrared Thermographic Data

The use of multi-dimensional finite volume heat conduction techniques for calculating aeroheating rates from measured global surface temperatures on hypersonic wind tunnel models was investigated. Both direct and inverse finite volume techniques were investigated and compared with the standard one-dimensional semi-infinite technique. Global transient surface temperatures were measured using an infrared thermographic technique on a 0.333-scale model of the Hyper-X forebody in the NASA Langley Research Center 20-Inch Mach 6 Air tunnel. In these tests the effectiveness of vortices generated via gas injection for initiating hypersonic transition on the Hyper-X forebody was investigated. An array of streamwise-orientated heating striations was generated and visualized downstream of the gas injection sites. In regions without significant spatial temperature gradients, one-dimensional techniques provided accurate aeroheating rates. In regions with sharp temperature gradients caused by striation patterns multi-dimensional heat transfer techniques were necessary to obtain more accurate heating rates. The use of the one-dimensional technique resulted in differences of 20% in the calculated heating rates compared to 2-D analysis because it did not account for lateral heat conduction in the model.

Daryabeigi, Kamran↗

NASA NTRS2002

Application of Local Discretization Methods in the NASA Finite-Volume General Circulation Model

We present the basic ideas of the dynamics system of the finite-volume General Circulation Model developed at NASA Goddard Space Flight Center for climate simulations and other applications in meteorology. The dynamics of this model is designed with emphases on conservative and monotonic transport, where the property of Lagrangian conservation is used to maintain the physical consistency of the computational fluid for long-term simulations. As the model benefits from the noise-free solutions of monotonic finite-volume transport schemes, the property of Lagrangian conservation also partly compensates the accuracy of transport for the diffusion effects due to the treatment of monotonicity. By faithfully maintaining the fundamental laws of physics during the computation, this model is able to achieve sufficient accuracy for the global consistency of climate processes. Because the computing algorithms are based on local memory, this model has the advantage of efficiency in parallel computation with distributed memory. Further research is yet desirable to reduce the diffusion effects of monotonic transport for better accuracy, and to mitigate the limitation due to fast-moving gravity waves for better efficiency.

Yeh, Kao-San↗

DOE OSTI2026

Real-time estimators for scattering observables: A full account of finite-volume errors for quantum simulation

The real-time correlators of quantum field theories can be directly probed through new approaches to simulation, such as quantum computing and tensor networks. This provides a new framework for computing scattering observables in lattice formulations of strongly interacting theories, such as lattice quantum chromodynamics. In this paper, we prove that the proposal of real-time estimators of scattering observables is universally applicable to all scattering observables of gapped quantum field theories. All finite-volume errors are exponentially suppressed, and the rate of this suppression is controlled by the regulator considered, namely, a displacement of the spectrum of the theory into the complex plane. A partial restoration of Lorentz symmetry by averaging over different boosts gives an additional suppression of finite volume errors. Our results also apply to the simulation of wave packet scattering, where a similar averaging is performed to construct the wave packets that regulate the finite volume effects. This result represents a necessary key step toward determining a broad class of scattering observables via quantum computing that are currently inaccessible via classical computing. Such observables are relevant for various applications, including hadron spectroscopy, hadron structure, and precision tests of the Standard Model. We also comment on potential applications of our results to traditional computational schemes.

Burbano, Ivan M. [University of California, Berkel↗

NASA NTRS1987

Composite grid and finite-volume LU implicit scheme for turbine flow analysis

A composite grid was generated in an attempt to improve grid quality for a typical turbine blade with large camber in terms of mesh control, smoothness, and orthogonality. This composite grid consists of the C grid (or O grid) in the immediate vicinity of the blade and the H grid in the upstream region and in the middle of the blade passage between the C grids. It provides a good boundary layer resolution around the leading edge region for viscous calculation, has orthogonality at the blade surface and slope continuity at the C-H (or O-H) interface, and has flexibility in controlling the mesh distribution in the upstream region without using excessive grid points. This composite grid eliminates the undesirable qualities of a single grid when generated for a typical turbine geometry. A finite-volume lower-upper (LU) implicit scheme can be used in solving for the turbine flows on the composite grid. This grid has a special grid node that is connected to more than four neighboring nodes in two dimensions and to more than six nodes in three dimensions. But the finite-volume approach poses no problem at the special point because each interior cell has only four neighboring cells in two dimensions and only six cells in three dimensions. The finite-volume LU implicit scheme was demonstrated to be robust and efficient for both external and internal flows in a broad flow regime.

Choo, Yung K.↗

NASA NTRS1987

Composite grid and finite-volume LU implicit scheme for turbine flow analysis

A composite grid was generated in an attempt to improve grid quality for a typical turbine blade with large camber in terms of mesh control, smoothness, and orthogonality. This composite grid consists of the C grid (or O grid) in the immediate vicinity of the blade and the H grid in the upstream region and in the middle of the blade passage between the C grids. It provides a good boundary layer resolution around the leading edge region for viscous calculation, has orthogonality at the blade surface and slope continuity at the C-H (or O-H) interface, and has flexibility in controlling the mesh distribution in the upstream region without using excessive grid points. This composite grid eliminates the undesirable qualities of a single grid when generated for a typical turbine geometry. A finite-volume lower-upper (LU) implicit schemes can be used in solving for the turbine flows on the composite grid. This grid has a special grid node that is connected to more than four neighboring nodes in two dimensions and to more than six nodes in three dimensions. But the finite-volume approach poses no problem at the special point because each interior cell has only four neighboring cells in two dimensions and only six cells in three dimensions. The finite-volume LU implicit scheme was demonstrated to be robust and efficient for both external and internal flows in a broad flow regime.

Choo, Yung K.↗

DOE OSTI2022

Advanced Finite-Volume Numerics and Source Term Assumptions for Kernel and G-Equation Modelling of Propane/Air Flames

Here G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device. We tested the explicit (sub-cycled) vs. implicit formulation for the standard transport operators (advection, diffusion, compressibility). In addition to the accurate flame swept-volume method for chemistry and species source term, we developed a more accurate estimator for the burnt/unburnt split cell composition. Then, we developed a signed-distance-function (SDF) based method which provides a more stable reinitialization of the level-set field at every time-step. We found that simplifying assumptions common to several G-Equation implementations, for straightforward terms such as compressibility and advection, lead to large errors in predicting the propagation of even laminar flames, with deviations up to ~300% in simulated vs. formulated flame speed. Conversely, the enhanced numerics enabled through the SDF field reinitialization and improved chemistry source term improve simulation stability and smooth flame propagation even with significantly larger solver time-steps.

42 ENGINEERING↗

DOE OSTI2024

Generalized boost transformations in finite volumes and application to Hamiltonian methods

The investigation of hadron interactions within lattice QCD has been facilitated by the well-known quantisation condition, linking scattering phase shifts to finite-volume energies. Additionally, the ability to utilise systems at finite total boosts has been pivotal in smoothly charting the energy-dependent behaviour of these phase shifts. The existing implementations of the quantization condition at finite boosts rely on momentum transformations between rest and moving frames, defined directly in terms of the energy eigenvalues. This energy dependence is unsuitable in the formulation of a Hamiltonian. In this work, we introduce a novel approach to generalise the three-momentum boost prescription, enabling the incorporation of energy-independent finite-volume Hamiltonians within moving frames. We demonstrate the application of our method through numerical comparisons, employing a phenomenological ππ scattering example.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Decay amplitudes to three hadrons from finite-volume matrix elements

We derive relations between finite-volume matrix elements and infinite-volume decay amplitudes, for processes with three spinless, degenerate and either identical or non-identical particles in the final state. This generalizes the Lellouch-Lüscher relation for two-particle decays and provides a strategy for extracting three-hadron decay amplitudes using lattice QCD. Unlike for two particles, even in the simplest approximation, one must solve integral equations to obtain the physical decay amplitude, a consequence of the nontrivial finite-state interactions. We first derive the result in a simplified theory with three identical particles, and then present the generalizations needed to study phenomenologically relevant three-pion decays. The specific processes we discuss are the CP-violating K → 3π weak decay, the isospin-breaking η → 3π QCD transition, and the electromagnetic γ* → 3π amplitudes that enter the calculation of the hadronic vacuum polarization contribution to muonic g - 2.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Signs of nonmonotonic finite-volume corrections to 𝑔 𝐴

We study finite-volume (FV) corrections to determinations of 𝑔 𝐴 via lattice quantum chromodynamics (QCD) using analytic results and numerical analysis. We observe that 𝑆⁢𝑈⁡(2) heavy Baryon chiral perturbation theory does not provide an unambiguous prediction for the sign of the FV correction, which is not surprising when one also considers large-𝑁 𝑐 constraints on the axial couplings. We further show that nonmonotonic FV corrections are naturally allowed when one considers either including explicit Δ-resonance degrees of freedom or one works to higher orders in the chiral expansion. We investigate the potential impact of these FV corrections with a precision study of 𝑔 𝐴 using models of FV corrections that are monotonic and nonmonotonic. Using lattice QCD data that is approximately at the 1% level of precision, we do not see significant evidence of nonmonotonic corrections. Looking forward to the next phase of lattice QCD calculations, we estimate that calculations that are between the 0.1% and 1% level of precision may be sensitive to these FV artifacts. Finally, we present an update of the CalLat prediction of 𝑔 𝐴 in the isospin limit with subpercent precision, 𝑔$^{QCD}_{𝐴}$ = 1.2674⁢(96).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1989

Viscous analysis of high speed flows using an upwind finite volume technique

A numerical study using an upwind finite volume technique to analyze high speed viscous flows is described. Two test cases were considered: a Mach 4.0 shock wave/laminar-boundary-layer interaction, and Mach 14.1 compression corners (15-deg and 24-deg ramp angle). PNS, TLNS, and NS equation results were obtained on each test case. As expected, PNS solutions are not adequate to resolve the flow physics. For the grids used in this study, there is no noticeable difference between the TLNS and full NS solutions. Both the TLNS and full NS solutions compared favorably with experimental data.

Ng, W. F.↗