Results for “finite volume discretization” — Page 18

DOE OSTI2026

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗

DOE OSTI2021

NRC Multiphysics Analysis Capability Deployment FY21: Part 3

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The deliverables completed for this report are: Deliverable 1c: the capability to model gas mixtures was added to Pronghorn. A test problem mimicking the conditions achieved in a depressurized loss of forced cooling (DLOFC) event was solved with both RELAP-5 and Pronghorn. Pronghorn employed a finite vol ume method with the Kurganov-Tadmor discretization. The comparison between the mass fraction spatial profiles computed with RELAP-5 and Pronghorn clearly shows the presence of numerical artifacts (i.e., overly diffusive behavior at low Mach numbers). We confirmed that the problem disappears at higher Mach numbers. We recommend future work on the implementation of a low Mach finite volume formulation to better treat low Mach number problems. Deliverable 2a: we demonstrated two approaches to model the radiation/conduction/natural convection heat transfer across a stagnant gas for the PBMR-400 design using Pronghorn. The first approach is based on the net radiation method, which relies on the computation of view factors with the Multiphysics Object-Oriented Simulation Environment (MOOSE) ray tracing capability. The second method is a traditional thermal resistance approach. The test problems include both 2D and 3D geometries. In all cases, the results show very good agreement during a DLOFC transient. This confirms that the faster thermal resistance method produces solutions that are equivalent to the net radiation method for this geometry. Deliverable 3d: we demonstrated the use of the advection kernel for the delayed neutron precursor equation in Griffin with a 2D MSFR model. The results appear physical but further verification is recommended. We also recommend the addition of conjugate heat transfer to compute the temperatures and model the thermomechanic behavior of the reflectors and other structures. Significant memory and performance issues were encountered in the 3D axisymmetric model. Future work is recommended in this area. Task 8g: this task allows multidimensional MOOSE applications to be coupled to system codes (RELAP-7 and SAM). We implemented a faster multiphysics iteration coupling algorithm, which provides an overall 6× acceleration of the 3D-1D coupling of the core multidi- mensional fluid flow solver and the 1D primary and secondary loop model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1997

Complex Volume Grid Generation Through the Use of Grid Reusability

This paper presents a set of surface and volume grid generation techniques which reuse existing surface and volume grids. These methods use combinations of data manipulations to reduce grid generation time, improve grid characteristics, and increase the capabilities of existing domain discretization software. The manipulation techniques utilize physical and computational domains to produce basis function on which to operate and modify grid character and smooth grids using Trans-Finite Interpolation, a vector interpolation method and parametric re-mapping technique. With these new techniques, inviscid grids can be converted to viscous grids, multiple zone grid adaption can be performed to improve CFD solver efficiency, and topological changes to improve modeling of flow fields can be done simply and quickly. Examples of these capabilities are illustrated as applied to various configurations.

Alter, Stephen J.↗

DOE OSTI2022

A Moving Embedded Boundary Approach for the Compressible Navier-Stokes Equations in a Block-Structured Adaptive Refinement Framework

A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement (SAMR) framework of AMReX [1], [91], [92]. We leverage the SAMR capability to obtain quantitatively accurate results whilst using robust, second-order finite volume schemes. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. A third-order least-squares formulation has been developed to compute the wall velocity gradients, and was found to significantly improve the performance of the solver in terms of the quantitative comparison of surface quantities such as the skin friction coefficient. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds - acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries - an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating that the numerical scheme is conservative with grid refinement, but is not discretely conservative. Viscous test cases involve that of a horizontally moving cylinder at Re = 40, an inline oscillating cylinder at Re = 100, and a transversely oscillating cylinder at Re = 185. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.

adaptive refinement↗

DOE OSTI2025

Parallel diffusion operator for magnetized plasmas with improved spectral fidelity

Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge for describing non-grid-aligned magnetic structures in Eulerian (grid-based) simulations. Here, the present paper describes and validates a new method for parallel diffusion in magnetized plasmas based on the anti-symmetry representation [Halpern and Waltz, Phys. Plasmas 25, 060703 (2018)]. In the anti-symmetry formalism, diffusion manifests as a flow operator involving the logarithmic derivative of the transported quantity. Qualitative plane wave analysis shows that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart. Numerical simulations comparing the new method against existing finite difference methods are carried out, showing significant improvement. In particular, we find that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods.

Anisotropic diffusion↗

NASA NTRS1972

Thermal elastoplastic structural analysis of non-metallic thermal protection systems

An incremental theory and numerical procedure to analyze a three-dimensional thermoelastoplastic structure subjected to high temperature, surface heat flux, and volume heat supply as well as mechanical loadings are presented. Heat conduction equations and equilibrium equations are derived by assuming a specific form of incremental free energy, entropy, stresses and heat flux together with the first and second laws of thermodynamics, von Mises yield criteria and Prandtl-Reuss flow rule. The finite element discretization using the linear isotropic three-dimensional element for the space domain and a difference operator corresponding to a linear variation of temperature within a small time increment for the time domain lead to systematic solutions of temperature distribution and displacement and stress fields. Various boundary conditions such as insulated surfaces and convection through uninsulated surface can be easily treated. To demonstrate effectiveness of the present formulation a number of example problems are presented.

Chung, T. J.↗

DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1993

The 3-D numerical study of airflow in the compressor/combustor prediffuser and dump diffuser of an industrial gas turbine

This paper describes the 3D computations of a flow field in the compressor/combustor diffusers of an industrial gas turbine. The geometry considered includes components such as the combustor support strut, the transition piece and the impingement sleeve with discrete cooling air holes on its surface. Because the geometry was complex and 3D, the airflow path was divided into two computational domains sharing an interface region. The body-fitted grid was generated independently in each of the two domains. The governing equations for incompressible Navier-Stokes equations were solved using the finite volume approach. The results show that the flow in the prediffuser is strongly coupled with the flow in the dump diffuser and vice versa. The computations also revealed that the flow in the dump diffuser is highly nonuniform.

Agrawal, Ajay K.↗

NASA NTRS2017

Numerical Modeling of Saturated Boiling in a Heated Tube

This paper describes a mathematical formulation and numerical solution of boiling in a heated tube. The mathematical formulation involves a discretization of the tube into a flow network consisting of fluid nodes and branches and a thermal network consisting of solid nodes and conductors. In the fluid network, the mass, momentum and energy conservation equations are solved and in the thermal network, the energy conservation equation of solids is solved. A pressure-based, finite-volume formulation has been used to solve the equations in the fluid network. The system of equations is solved by a hybrid numerical scheme which solves the mass and momentum conservation equations by a simultaneous Newton-Raphson method and the energy conservation equation by a successive substitution method. The fluid network and thermal network are coupled through heat transfer between the solid and fluid nodes which is computed by Chen's correlation of saturated boiling heat transfer. The computer model is developed using the Generalized Fluid System Simulation Program and the numerical predictions are compared with test data.

Majumdar, Alok↗

DOE OSTI2024

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1976

The Relation of Finite Element and Finite Difference Methods

Finite element and finite difference methods are examined in order to bring out their relationship. It is shown that both methods use two types of discrete representations of continuous functions. They differ in that finite difference methods emphasize the discretization of independent variable, while finite element methods emphasize the discretization of dependent variable (referred to as functional approximations). An important point is that finite element methods use global piecewise functional approximations, while finite difference methods normally use local functional approximations. A general conclusion is that finite element methods are best designed to handle complex boundaries, while finite difference methods are superior for complex equations. It is also shown that finite volume difference methods possess many of the advantages attributed to finite element methods.

Vinokur, M.↗

NASA NTRSDate not supplied

Aerothermal Shape Optimization of Actively-Cooled Battery Packs using Conjugate Heat Transfer

Thermal management for battery is important for electric aircraft because battery temperature is critically important to vehicle safety, and it also has direct impact on the efficiency of the battery system. Because ambient air is a readily available resource for aircraft, this paper considers an active cooling concept with forced convection of ambient air through the battery pack. Conjugate heat transfer analysis is used to solve the coupled aero-thermal problem, which consists of a finite-volume computational fluid dynamics solver for the fluid domain, and a conduction heat transfer solver for the solid domain. A mixed Neumann and Dirichlet boundary condition is developed for the fluid-solid interface, which allows the solid domain to completely submerge in the fluid domain. A gradient-based optimization method is adopted, and the discrete adjoint approach implemented in DAFoam is used to efficiently compute the gradients. The aero-thermal coupling for primal analysis and gradient computation is handled using the OpenMDAO-based MPhys framework. A constant heat source is prescribed for the battery cells, and the battery shape (design variable) is optimized to minimize cooling pump power and battery weight (composite objective function) while keeping the battery temperature below a threshold (constraint). The optimized design achieves a 44.6% and 1.5% reduction in the cooling pump power and battery weight, respectively, and the maximal temperature constraint is satisfied. This work has the potential to reduce battery-pack weight, improve performance, and reduce the weight of thermal management systems for electric vertical take-off and landing aircraft.

thermal management↗

NASA NTRSDate not supplied

Nodal Modeling of Tank Pressurization and Draining using a Multi-Node-Ullage Approach

Pressurized expulsion tests of liquid methane were modeled by a nodal code using a multinode ullage approach. Generalized Fluid System Simulation Program (GFSSP), a finite volume based nodal code was used to model the expulsion of liquid methane from a 1.52-meter (5 ft) diameter spherical tank by pressurizing with helium and gaseous methane. The purpose of the model was to estimate the amount of pressurant required and amount of condensation and evaporation of methane during the operation. The ullage was discretized into multiple nodes and each fluid node was connected to solid nodes. Both the fluid and solid nodes grow as the tank drains. The heat and mass transfer between fluid and solid node in the tank ullage was computed. The heat and mass transfer between the bottom ullage node and the liquid methane was also computed. The model predictions were compared with data from twelve test cases. The predicted pressurant consumption for helium pressurization compares with test data with a Mean Absolute Percent Error (MAPE) of 3 %. For autogenous pressurization, where gaseous methane was used to pressurize liquid methane, the predicted pressurant consumption compares with a MAPE of 6%.

Cryogenic Propellant↗

NASA NTRSDate not supplied

Nodal Modeling of Helium Pressurization and Autogenous Pressurization and Draining using a Multi-Node-Ullage Approach

Pressurized expulsion tests of liquid methane were modeled by a nodal code using a multinode ullage approach. Generalized Fluid System Simulation Program (GFSSP), a finite volume based nodal code was used to model the expulsion of liquid methane from a 1.52-meter (5 ft) diameter spherical tank by pressurizing with helium and gaseous methane. The purpose of the model was to estimate the amount of pressurant required and amount of condensation and evaporation of methane during the operation. The ullage was discretized into multiple nodes and each fluid node was connected to solid nodes. Both the fluid and solid nodes grow as the tank drains. The heat and mass transfer between fluid and solid node in the tank ullage was computed. The heat and mass transfer between the bottom ullage node and the liquid methane was also computed. The model predictions were compared with data from twelve test cases. The predicted pressurant consumption for helium pressurization compares with test data with a Mean Absolute Percent Error (MAPE) of 3%. For autogenous pressurization, where gaseous methane was used to pressurize liquid methane, the predicted pressurant consumption compares with a MAPE of 6%.

Nodal Model↗

NASA NTRSDate not supplied

Aerothermal Shape Optimization of Actively-Cooled Battery Packs Using Conjugate Heat Transfer

Thermal management for battery is important for electric aircraft because battery temperature is critically important to vehicle safety, and it also has direct impact on the efficiency of the battery system. Because ambient air is a readily available resource for aircraft, this paper considers an active cooling concept with forced convection of ambient air through the battery pack. Conjugate heat transfer analysis is used to solve the coupled aero-thermal problem, which consists of a finite-volume computational fluid dynamics solver for the fluid domain, and a conduction heat transfer solver for the solid domain. A mixed Neumann and Dirichlet boundary condition is developed for the fluid-solid interface, which allows the solid domain to completely submerge in the fluid domain. A gradient-based optimization method is adopted, and the discrete adjoint approach implemented in DAFoam is used to efficiently compute the gradients. The aero-thermal coupling for primal analysis and gradient computation is handled using the OpenMDAO-based MPhys framework. A constant heat source is prescribed for the battery cells, and the battery shape (design variable) is optimized to minimize cooling pump power and battery weight (composite objective function) while keeping the battery temperature below a threshold (constraint). The optimized design achieves a 44.6% and 1.5% reduction in the cooling pump power and battery weight, respectively, and the maximal temperature constraint is satisfied. This work has the potential to reduce battery-pack weight, improve performance, and reduce the weight of thermal management systems for electric vertical take-off and landing aircraft.

heat transfer↗

NASA NTRSDate not supplied

Edge Based Viscous Method for Node-Centered Formulations

This paper presents a novel, efficient, conservative, edge-based method for evaluation of mean flow viscous fluxes and turbulence-model diffusion terms of the Reynolds-averaged Navier-Stokes equations on tetrahedral grids. The new method is implemented in a practical, node-centered, finite-volume computational fluid dynamics solver. The baseline finite-volume scheme that is equivalent to a second-order accurate finite-element Galerkin approximation of viscous stresses is reformulated. The order of operations to compute the cell-based Green-Gauss gradients is changed to combine the operations by edge, which leads to an equivalent formulation on tetrahedral grids, improves efficiency, and preserves the compact discretization stencil based on the nearest neighbors. The computational results presented in this paper verify the implementation of this edge-based method by comparing its accuracy and iterative convergence with those of the well verified and validated baseline formulation. Efficiency gains for residual and Jacobian evaluations result in significant reduction of time to solution. This novel edge-based formulation on tetrahedra can be seamlessly combined with the baseline formulation on cells of other types for computing solutions on mixed-element grids.

Edge Based↗

NASA NTRS2023

A Hybrid Transfer Function Procedure for Broadband Auralization Within Small Flight Vehicle Interiors

Auralization of the sound fields within flight vehicle interiors is of interest to aircraft designers considering crew and passenger ride quality. Computer models can be used to predict the filtering effects that the fuselage structure and interior volume have on the exterior sound field in order to simulate the resulting cabin interior noise. However, this can become expensive when a wide frequency range of analysis suitable for passenger auralization is desired. This is particularly the case when using deterministic modeling methods such as finite element analysis (FEA) wherein prohibitively high levels of detail and finite element discretizations are required to resolve the wavelengths at high frequencies. On the contrary, Statistical Energy Analysis (SEA), while not capturing important modal behaviors at low frequencies, is better suited for higher frequency analyses where structural and interior dynamics exhibit high modal overlap and spatially averaged physical quantities suffice. For this reason, a hybrid procedure is considered here that utilizes FE modal solutions at lower frequencies and SEA results at higher frequencies. A transfer function dataset containing the exterior to interior vibroacoustic filtering effect is calculated accordingly and subsequently stored to be used for a variety of exterior dynamic load cases. At low frequencies, the transfer function dataset detail is related to the level of FE model refinement whereas at high frequencies it functions as an equalizer array with a level of detail corresponding to the SEA subsystem partitioning scheme. The above procedure is described herein and demonstrated on a six passenger flight vehicle with eight propulsors in a hover flight condition.

auralization↗

NASA NTRS2023

A Hybrid Transfer Function Procedure for Broadband Auralization Within Small Flight Vehicle Interiors

Auralization of the sound fields within flight vehicle interiors is of interest to aircraft designers considering crew and passenger ride quality. Computer models can be used to predict the filtering effects that the fuselage structure and interior volume have on the exterior sound field in order to simulate the resulting cabin interior noise. However, this can become expensive when a wide frequency range of analysis suitable for passenger auralization is desired. This is particularly the case when using deterministic modeling methods such as finite element analysis (FEA) wherein prohibitively high levels of detail and finite element discretizations are required to resolve the wavelengths at high frequencies. On the contrary, Statistical Energy Analysis (SEA), while not capturing important modal behaviors at low frequencies, is better suited for higher frequency analyses where structural and interior dynamics exhibit high modal overlap and spatially averaged physical quantities suffice. For this reason, a hybrid procedure is considered here that utilizes FE modal solutions at lower frequencies and SEA results at higher frequencies. A transfer function dataset containing the exterior to interior vibroacoustic filtering effect is calculated accordingly and subsequently stored to be used for a variety of exterior dynamic load cases. At low frequencies, the transfer function dataset detail is related to the level of FE model refinement whereas at high frequencies it functions as an equalizer array with a level of detail corresponding to the SEA subsystem partitioning scheme. The above procedure is described herein and demonstrated on a six passenger flight vehicle with eight propulsors in a hover flight condition.

auralization↗