SEARCH · Engineering Papers

Results for “finite volume method”

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 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

NASA NTRS1999

A Pseubo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

Morrison, J. H.↗

NASA NTRS1999

A Pseudo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

White, J. A.↗

DOE OSTI2025

Lieb-Robinson Bounds with Exponential-in-Volume Tails

Lieb-Robinson bounds demonstrate the emergence of locality in many-body quantum systems. Intuitively, Lieb-Robinson bounds state that, with local or exponentially decaying interactions, the correlation that can be built up between two sites separated by distance 𝑟 after a time 𝑡 decays as exp (𝑣⁢𝑡 −𝑟), where 𝑣 is the emergent Lieb-Robinson velocity. In many problems, it is important to also capture how much of an operator grows to act on 𝑟 𝑑 sites in 𝑑 spatial dimensions. Perturbation theory and cluster expansion methods suggest that, at short times, these volume-filling operators are suppressed as exp (−𝑟 𝑑 ). We confirm this intuition, showing that, for 𝑟 >𝑣⁢𝑡, the volume-filling operator is suppressed by exp (−(𝑟−𝑣⁢𝑡) 𝑑 /(𝑣⁢𝑡) 𝑑−1 ). This closes a conceptual and practical gap between the cluster expansion and the Lieb-Robinson bound. We then present two very different applications of this new bound. Firstly, we obtain improved bounds on the classical computational resources necessary to simulate many-body dynamics with error tolerance 𝜀 for any finite time 𝑡: as 𝜀 becomes sufficiently small, only 𝜀 −O⁡(𝑡 𝑑−1 ) resources are needed. A protocol that likely saturates this bound is given. Secondly, we prove that disorder operators have volume-law suppression near the “solvable (Ising) point” in quantum phases with spontaneous symmetry breaking, which implies a new diagnostic for distinguishing many-body phases of quantum matter.

computational complexity↗

DOE OSTI2025

Resonance ultrasound prediction of residual stress within a hybrid layer for additively manufactured samples

Hybrid additive manufacturing (AM) involves secondary processes or energy sources to alter specified locations within the build volume. Each hybrid step can refine the grain size, increase dislocation density, or modify residual stresses. Typically, the changes in mechanical properties are not confined within a single layer but have a compounding effect on preceding layers. Existing methods of measuring AM residual stress are limited in terms of their sensitivity, or they are destructive measurements. We propose using resonant ultrasound spectroscopy (RUS) to measure the residual stress in hybrid-AM components noninvasively, based on changes to the resonances, compared to a stress-free component. In this paper, we use finite element models to simulate residual stress in hybrid-AM components and to examine the sensitivity of RUS measurements in terms of frequency shifts and mode shapes with respect to single hybrid layers. Then, the RUS results are used to predict stress for a layer at a known location with unknown stress. Here, the approach highlights the capabilities of RUS to address an AM characterization challenge.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Electronic Scanning Strategies in Adaptive Electrical Capacitance Volume Tomography: Tradeoffs and Prospects

Electrical Capacitance Volume Tomography (ECVT) has been applied for imaging of multiphase flows found in industrial applications. The ill-posed nature of the image reconstruction problem in ECVT and the consequent low resolution can be alleviated by employing electronic scanning enabled by electrode segmentation and reconfiguration during data acquisition. Here we study electronic scanning strategies that mimic physical rotation and shifting of the sensor along its symmetry axis. First, we study the feasibility of electronic scanning by analyzing the capacitance transducer circuit in SPICE. Then, we simulate electronic scanning using the finite element method for different electrode shapes, for which we compare the image reconstruction results and acquisition time. We find a noticeable improvement in image resolution for the scanning cases over conventional ECVT. Lastly, among the scanning cases, we emphasize a particular electrode shape that provides the best image resolution along with the minimal amount of acquisition time.

47 OTHER INSTRUMENTATION↗

NASA NTRS1995

Parallelized modelling and solution scheme for hierarchically scaled simulations

This two-part paper presents the results of a benchmarked analytical-numerical investigation into the operational characteristics of a unified parallel processing strategy for implicit fluid mechanics formulations. This hierarchical poly tree (HPT) strategy is based on multilevel substructural decomposition. The Tree morphology is chosen to minimize memory, communications and computational effort. The methodology is general enough to apply to existing finite difference (FD), finite element (FEM), finite volume (FV) or spectral element (SE) based computer programs without an extensive rewrite of code. In addition to finding large reductions in memory, communications, and computational effort associated with a parallel computing environment, substantial reductions are generated in the sequential mode of application. Such improvements grow with increasing problem size. Along with a theoretical development of general 2-D and 3-D HPT, several techniques for expanding the problem size that the current generation of computers are capable of solving, are presented and discussed. Among these techniques are several interpolative reduction methods. It was found that by combining several of these techniques that a relatively small interpolative reduction resulted in substantial performance gains. Several other unique features/benefits are discussed in this paper. Along with Part 1's theoretical development, Part 2 presents a numerical approach to the HPT along with four prototype CFD applications. These demonstrate the potential of the HPT strategy.

Padovan, Joe↗

NASA NTRS1999

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗

DOE OSTI2026

Broad frequency tuning of a Nb$_{3}$Sn superconducting microwave cavity for dark matter searches

We demonstrate a novel broad-frequency tuning mechanism for superconducting microwave cavities designed for dark matter searches. Using a Nb$_3$Sn-coated cigar-shaped cavity operating at approximately 9 GHz, we achieve continuous frequency tuning exceeding 1 GHz by mechanically separating the two cavity halves: a "tuning-by-opening" technique. Finite-element method simulations predict that radiative losses do not degrade the quality factor even for large openings, as a closed cavity with an intrinsic quality factor of $10^7$ maintains this value for apertures up to 9 mm, corresponding to a tuning range from 9.0 to 7.5 GHz. Experimental validation using both copper ring spacers and a continuous sliding mechanism confirms $Q_0$ values exceeding the dark matter quality factor across the entire explored frequency range, despite mechanical imperfections and film non-uniformities. This tuning approach avoids inserting elements into the resonant volume, making it particularly suitable for high-Q superconducting cavities in axion haloscope experiments and readily applicable to REBCO-based implementations capable of operating in multi-tesla magnetic fields.

Maiello, D. [Padua U.; INFN, Padua] (ORCID:0009000↗

DOE OSTI2023

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗

NASA NTRS1983

A finite element model for sound transmission through panels

A finite element method (FEM) is applied to predicting coupled frequencies and pressures within an acoustic cavity in order to characterize sound transmission through a panel. Structural equations of motion are defined and the FEM model is configured with four-noded plate elements, the lateral displacement and two slopes being the unknowns at every node. Each element then has 12 degrees of freedom (DOF) and the displacement variation is expressed by a 12-term nonconforming polynomial. A consistent mass matrix is used to represent the panel mass matrix and a wave equation governs the acoustic volume. Analysis of pressure and displacement over the panel yields a square coupling matrix, and an eigenanalysis leads to a solution of the forced vibration problem.

Ramakrishnan, J. V.↗

NASA NTRS2017

High-order Accurate Finite-Volume Formulations for the Pressure Gradient Force in Layered Ocean Models

Discretisation of the horizontal pressure gradient force in layered ocean models is a challenging task, with non-trivial interactions between the thermodynamics of the fluid and the geometry of the layers often leading to numerical difficulties. We present two new finite-volume schemes for the pressure gradient operator designed to address these issues. In each case, the horizontal acceleration is computed as an integration of the contact pressure force that acts along the perimeter of an associated momentum control-volume. A pair of new schemes are developed by exploring different control-volume geometries. Non-linearities in the underlying equation-of-state definitions and thermodynamic profiles are treated using a high-order accurate numerical integration framework, designed to preserve hydrostatic balance in a non-linear manner. Numerical experiments show that the new methods achieve high levels of consistency, maintaining hydrostatic and thermobaric equilibrium in the presence of strongly-sloping layer geometries, non-linear equations-of-state and non-uniform vertical stratification profiles. These results suggest that the new pressure gradient formulations may be appropriate for general circulation models that employ hybrid vertical coordinates andor terrain-following representations.

pressure gradient force↗

NASA NTRS2022

Laser Patterning Assisted Devitrification and Domain Engineering of Amorphous and Nanocrystalline Alloys

We have investigated laser-patterning and devitrification as a novel method to realize emergent properties in modified regions at the ribbon surface consisting of periodic localized laser spots thermally treated under rapid heating and cooling conditions. Fe-based amorphous ribbons were annealed for systematically varied laser power. Extremely high heating and cooling rates on the order of 10(exp 8)-10(exp 10)K/s have been estimated by finite element analysis simulations. Observations suggest surface melting followed by rapid solidification fast enough to quench and form an amorphous structure upon cooling. Diffusion of heat occurs from the laser irradiated surface region so that a finite volume of material surrounding this region rises in temperature above the crystallization temperature relevant for conventional isothermal annealing experiments. The underlying mechanism of laser annealing and consequences of heat transients involved are discussed in terms of impacts upon micro/nanostructure, residual stresses, and magnetic domain structure surrounding the laser irradiated region. The study illustrates the potential to exploit spatially optimized phase transformations in a scalable manufacturing process of amorphous and nanocrystalline alloys to locally access otherwise inaccessible extreme heating and cooling rates.

Amorphous ribbons↗

NASA NTRS1996

A Note on Multigrid Theory for Non-nested Grids and/or Quadrature

We provide a unified theory for multilevel and multigrid methods when the usual assumptions are not present. For example, we do not assume that the solution spaces or the grids are nested. Further, we do not assume that there is an algebraic relationship between the linear algebra problems on different levels. What we provide is a computationally useful theory for adaptively changing levels. Theory is provided for multilevel correction schemes, nested iteration schemes, and one way (i.e., coarse to fine grid with no correction iterations) schemes. We include examples showing the applicability of this theory: finite element examples using quadrature in the matrix assembly and finite volume examples with non-nested grids. Our theory applies directly to other discretizations as well.

Douglas, C. C.↗

DOE OSTI2024

Review of recent activities with MOOSE, an open-source finite element & finite volume multi-fidelity simulation framework

Modeling and simulation are an increasing part of engineering. This is undoubtedly driven by the high costs of constructing experimental facilities, but also enabled by the exponential increase in computing powers over the last decades, which allows computational models to be closer than ever to reality. One of the main drivers for the development of MOOSE is supporting advanced nuclear reactor simulations. A challenging aspect of modeling advanced nuclear reactors is the plurality of physics involved, including neutronics, thermal hydraulics and fuel performance. These physics are all coupled to some extent and are generally solved in a sequential but iterative fashion. The United States (U.S.) national laboratories have been developing MOOSE, an open source multiphysics framework since its inception at the Idaho National Laboratory (INL) in 2008. This framework enables seamless coupling of multiphysics simulations and facilitates the implementation of new physics and material governing laws. It is continuously expanded with novel numerical methods and new pre-implemented physics module. Numerous applications, developed within the Department of Energy (DOE) laboratories, academia, and industry, including outside of nuclear engineering, have been developed to study specialized physics problems. International collaborations are welcome on this open-source modeling and simulation project.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Three-body resonances in pionless effective field theory

We investigate the appearance of resonances in three-body systems using pionless effective field theory at leading order with two complementary methods. The Faddeev equation is analytically continued to the unphysical sheet adjacent to the positive real energy axis using a contour rotation. We consider both the three-boson system and the three-neutron system. For the former, we calculate the trajectory of Borromean three-body Efimov states turning into resonances as they cross the three-body threshold. For the latter, we find no sign of three-body resonances or virtual states at leading order. Furthermore, this result is validated by exploring the level structure of three-body states in a finite volume approach.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1980

Numerical calculation of transonic axial turbomachinery flows

A numerical method and the results of a computer program are presented for solving an exact, three dimensional, full potential equation that models rotating and nonrotating inviscid, absolutely irrotational, homentropic flows. Besides calculating the flows through an arbitrarily shaped rotor or stator blade row mounted on an axisymmetric hub and confined in an axisymmetric duct, the computer program is also capable of analyzing flow fields about arbitrarily shaped wing body combinations, propellers, helicopter rotors in hover, and wind turbine rotors. The governing equation is solved numerically in a fully conservative form by using an artificial time concept, a finite volume technique, rotated type dependent differencing, successive line overrelaxation, and sequential boundary conforming grid refinement. An artificial viscosity is added in fully conservative form, and an initial guess for the potential field is applied, as determined by a two dimensional cascade analysis.

Dulikravich, D. S.↗

NASA NTRS1993

High accuracy computation of fluid-structure interaction in transonic cascades

A coupling strategy for simulating fluid-structure interaction phenomena is formulated and applied to the prediction of flutter in transonic cascades. The flow is governed by the Euler equations and discretized using a finite volume flux-splitting scheme. The structure is modeled using an isoparametric finite element formulation. The coupling strategy successfully reconciles these two formulations at the fluid-structure interface by enforcing both kinematic and kinetic boundary conditions. In particular, the conservation laws applicable to the combined fluid-structure system are preserved across the interface. Since the primary mechanism driving aeroelastic phenomena involves energy exchange occurring at the interface, this highly accurate coupling mechanism is believed to improve the predictive capability of the scheme. The coupled equations are advanced simultaneously in time using an implicit time integration method. Results obtained using the coupling method are presented for cascade geometries operating in transonic flow.

Boschitsch, Alexander↗

NASA NTRS1995

Unstructured viscous flow solution using adaptive hybrid grids

A three dimensional finite volume scheme based on hybrid grids containing both tetrahedral and hexahedral cells is presented. The application to hybrid grids offers the possibility to combine the flexibility of tetrahedral meshes with the accuracy of hexahedral grids. An algorithm to compute a dual mesh for the entire computational domain was developed. The dual mesh technique guarantees conservation in the whole flow field even at interfaces between hexahedral and tetrahedral domains and enables the employment of an accurate upwind flow solver. The hybrid mesh can be adapted to the solution by dividing cells in areas of insufficient resolution. The method is tested on different viscous and inviscid cases for hypersonic, transonic and subsonic flows.

Galle, Martin↗