Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical schemes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Reliability of Complex Nonlinear Numerical Simulations

This work describes some of the procedure to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Yee, H. C.↗
NASA NTRS2005

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations.

Yee, H. C.↗
NASA NTRS1981

Numerical boundary condition procedure for the transonic axisymmetric inverse problem

Two types of boundary condition procedures for the axisymmetric inverse problem are described. One is a Neumann type boundary condition (analogous to the analysis problem) and the other is a Dirichlet type boundary conditon, both requiring special treatments to make the inverse scheme numerically stable. The dummy point concept is utilized in implementing both. Results indicate the Dirichlet type inverse boundary condition is more robust and conceptually simpler to implement than the Neumann type procedure. A few results demonstrating the powerful capability of the newly developed inverse method that can handle both shocked as well as shockless body design are included.

Shankar, V.↗
NASA NTRS1981

Mixed spectral/finite difference approximations for slightly viscous flows

A numerical scheme is presented, which employs a standard finite difference approximation for the viscous terms in high Reynolds number flows, and resorts to pseudo-spectral methods for the advection terms to greatly reduce the major source of numerical error without dramatically increasing computational cost. The spectral/finite difference (SFD) method evaluates the advection term and second-order differences to evaluate the diffusion term. The fully finite difference (FFD) method with second-order central differences on both terms is also used. The SFD method can handle strong shocks and can outperform the FFD method at moderate viscosity.

Zang, T. A.↗
NASA NTRS1983

Three-dimensional unsteady Euler equation solutions using flux vector splitting

A method for numerically solving the three dimensional unsteady Euler equations using flux vector splitting is developed. The equations are cast in curvilinear coordinates and a finite volume discretization is used. An explicit upwind second-order predictor-corrector scheme is used to solve the discretized equations. The scheme is stable for a CFL number of two and local time stepping is used to accelerate convergence for steady-state problems. Characteristic variable boundary conditions are developed and used in the far field and at surfaces. No additional dissipation terms are included in the scheme. Numerical results are compared with results from an existing three dimensional Euler code and experimental data.

Whitfield, D. L.↗
NASA NTRS1992

Unsteady flowfield simulation of ducted prop-fan configurations

A technique for the simulation of unsteady flows in and around complex rotating machinery is presented. Additional domain decomposition mechanisms are introduced which extend the range of applicability of software developed for the time-accurate simulation of rotating machinery flowfields. The flow models uses the unsteady 3D Euler equations, discretized as a finite-volume method, utilizing a high-resolution approximate Riemann solver for cell interface flux definitions. Multiblock domain decomposition is used to partition the field radially, axially, as well as circumferentially into an ordered arrangement of blocks which exhibit varying degrees of similarity. A general high-order numerical scheme is applied to satisfy the geometric conservation law. Two configurations are presented - ducted single rotation prop-fan and a rotor-deswirl vane combination which form a single stage fan. Comparisons are made to other numerical solutions for these geometries and to available experimental data.

Janus, J. M.↗
NASA NTRS1983

A rotationally biased upwind difference scheme for the Euler equations

The upwind difference schemes of Godunov, Osher, Roe and van Leer are able to resolve one dimensional steady shocks for the Euler equations within one or two mesh intervals. Unfortunately, this resolution is lost in two dimensions when the shock crosses the computing grid at an oblique angle. To correct this problem, a numerical scheme was developed which automatically locates the angle at which a shock might be expected to cross the computing grid and then constructs separate finite difference formulas for the flux components normal and tangential to this direction. Numerical results which illustrate the ability of this method to resolve steady oblique shocks are presented.

Davis, S. F.↗
NASA NTRS1984

A rotationally biased upwind difference scheme for the Euler equations

The upwind difference schemes of Godunov, Osher, Roe and van Leer are able to resolve one dimensional steady shocks for the Euler equations within one or two mesh intervals. Unfortunately, this resolution is lost in two dimensions when the shock crosses the computing grid at an oblique angle. To correct this problem, a numerical scheme was developed which automatically locates the angle at which a shock might be expected to cross the computing grid and then constructs separate finite difference formulas for the flux components normal and tangential to this direction. Numerical results which illustrate the ability of this method to resolve steady oblique shocks are presented.

Davis, S. F.↗
NASA NTRS1984

Three-dimensional unsteady Euler equations solution using flux vector splitting

A method for numerically solving the three-dimensional unsteady Euler equations using flux vector splitting is developed. The equations are cast in curvilinear coordinates and a finite volume discretization is used. An explicit upwind second-order predictor-corrector scheme is used to solve the discretized equations. The scheme is stable for a CFL number of 2 and local time stepping is used to accelerate convergence for steady-state problems. Characteristic variable boundary conditions are developed and used in the far-field and at surfaces. No additional dissipation terms are included in the scheme. Numerical results are compared with results from an existing three-dimensional Euler code and experimental data.

Whitfield, D. L.↗
NASA NTRS1987

Incompressible Navier-Stokes computations of vortical flows over double-delta wings

An implicit flux-difference splitting scheme is used to compute incompressible laminar vortical flows about a thin round-edged double-delta wing with 80 deg and 60 deg leading-edged sweep for the strake and the main wing respectively. The numerical scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. It is second-order accurate spatially by applying TVD-like upwind discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. The interaction between the two primary vortices emanating from the apex and kink leading-edges is successfully simulated. The computed trajectories of vortical cores compare well with experimental data.

Hsu, Chung-Hao↗
NASA NTRS1994

Investigation of supersonic jet plumes using an improved two-equation turbulence model

Supersonic jet plumes were studied using a two-equation turbulence model employing corrections for compressible dissipation and pressure-dilatation. A space-marching procedure based on an upwind numerical scheme was used to solve the governing equations and turbulence transport equations. The computed results indicate that two-equation models employing corrections for compressible dissipation and pressure-dilatation yield improved agreement with the experimental data. In addition, the numerical study demonstrates that the computed results are sensitive to the effect of grid refinement and insensitive to the type of velocity profiles used at the inflow boundary for the cases considered in the present study.

Lakshmanan, B.↗
DOE OSTI2024

Investigation of two-dimensional radio-frequency sheath properties using a microscale fluid model

Abstract In previous work (Kohno and Myra 2023 Comput. Phys. Commun. 291 108841), we developed a numerical scheme based on a two-dimensional microscale radio-frequency (RF) sheath model with periodically curved wall boundaries. Here, we expand the capability of this scheme through modification of the boundary conditions (BCs) on the conducting walls, which allows the ion flow to turn back to the plasma at locations on the walls where the electromagnetic force on the ions is reversed from its usual direction. Numerical simulations are carried out to investigate the dependences of the surface-integrated admittances on the wall bump height, ion magnetization, ion mobility, and the magnetic field angle, and to visualize the sheath structures in several cases. One of the main results is the ion cyclotron admittance resonance observed under the condition of low ion mobility (high normalized frequency). It is shown that the amplitude of the resonance peak depends on the wall bump height and the ion velocity is reversed on the sides of the bump in an RF cycle for the resonance cases. Furthermore, the differences in the admittances between the one- and two-dimensional microscale models are assessed for the purpose of understanding non-locality of the sheath near the wall surface for the parameters considered in this study. This information will be essential for improving the sheath BC for macroscale calculations in the future.

Kohno, H. (ORCID:0000000332690010)↗
NASA NTRS1973

Dynamic loads due to moving liquid.

The marker-and-cell method has been extended to simulate the flow and velocity fields of a moving liquid in an axisymmetric container. Dynamic loads induced by the moving liquid are calculated from the pressure distribution along the container wall as functions of time. Using viscosity and equivalent gravitational acceleration as parameters, several cases of the transient flow of a liquid in a container are studied. Subsequently, a three-dimensional marker-and-cell formulation is obtained and a numerical scheme to implement the formulation is developed. The resulting computer program is used to simulate the three-dimensional flow of a liquid in a rectangular container. The study shows that dynamic loads induced by a moving liquid can be computed numerically.

Feng, G. C.↗
NASA NTRS1984

Enthalpy damping for the steady Euler equations

For inviscid steady flow problems where the enthalpy is constant at steady state, it was previously proposed to use the difference between the local enthalpy and the steady state enthalpy as a driving term to accelerate convergence of iterative schemes. This idea is analyzed, both on the level of the partial differential equation and on the level of a particular finite difference scheme. It is shown that for the two-dimensional unsteady Euler equations, a hyperbolic system with eigenvalues on the imaginary axis, there is no enthalpy damping strategy which moves all the eigenvalues into the open left half plane. For the numerical scheme, however, the analysis shows and examples verify that enthalpy damping is potentially effective in accelerating convergence to steady state.

Jespersen, D. C.↗
NASA NTRS1985

Enthalpy damping for the steady Euler equations

For inviscid steady flow problems where the enthalpy is constant at steady state, it was previously proposed to use the difference between the local enthalpy and the steady state enthalpy as a driving term to accelerate convergence of iterative schemes. This idea is analyzed, both on the level of the partial differential equation and on the level of a particular finite difference scheme. It is shown that for the two-dimensional unsteady Euler equations, a hyperbolic system with eigenvalues on the imaginary axis, there is no enthalpy damping strategy which moves all the eigenvalues into the open left half plane. For the numerical scheme, however, the analysis shows and examples verify that enthalpy damping is potentially effective in accelerating convergence to steady state.

Jespersen, D. C.↗
NASA NTRS1986

Interaction of decaying trailing vortices in ground shear

The drift of trailing vortices in a cross wind near to the ground is modeled by an unsteady, two-dimensional, rotational flowfield with a concentration of large vorticity in spots having finite total strength. The problem is analyzed by a combination of the method of matched asymptotic analyses for the decay of the vortical spots and the Euler solution for the unsteady rotational flow. A primary feature of the numerical scheme is that the grid size and time step depend only the length and velocity scales of the background flow and is independent of the effective core size of a vortical spot which can be much smaller than the grid size. Numerical results are presented to demonstrated the strong interaction between the trajectjories of the vortical spots and the redistribution of vorticity in the background flowfield.

Liu, C. H.↗
NASA NTRS1993

Effects of spatial order of accuracy on the computation of vortical flowfields

The effect of the order-of-accuracy, used for the spatial discretization, on the resolution of the leading edge vortices over sharp-edged delta wings is investigated. The flowfield is computed using a viscous/inviscid zonal approach. The viscous flow in the vicinity of the wing is computed using the conservative formulation of the compressible, thin-layer Navier-Stokes equations. The leeward-side vortical flowfield and the other flow regions away from the surface are computed as inviscid. The time integration is performed with both an explicit fourth-order Runge-Kutta scheme and an implicit, factorized, iterative scheme. High-order-accurate inviscid fluxes are computed using both a conservative and a non-conservative (primitive variable) formulation. The nonlinear, inviscid terms of the primitive variable form of the governing equations are evaluated with a finite-difference numerical scheme based on the sign of the eigenvalues. High-order, upwind-biased, finite difference formulas are used to evaluate the derivatives of the nonlinear convective terms. Computed results are compared with available experimental data, and comparisons of the flowfield in the vicinity of the vortex cores are presented.

Ekaterinaris, J. A.↗
NASA NTRS2013

Numerical Modeling of Ablation Heat Transfer

A unique numerical method has been developed for solving one-dimensional ablation heat transfer problems. This paper provides a comprehensive description of the method, along with detailed derivations of the governing equations. This methodology supports solutions for traditional ablation modeling including such effects as heat transfer, material decomposition, pyrolysis gas permeation and heat exchange, and thermochemical surface erosion. The numerical scheme utilizes a control-volume approach with a variable grid to account for surface movement. This method directly supports implementation of nontraditional models such as material swelling and mechanical erosion, extending capabilities for modeling complex ablation phenomena. Verifications of the numerical implementation are provided using analytical solutions, code comparisons, and the method of manufactured solutions. These verifications are used to demonstrate solution accuracy and proper error convergence rates. A simple demonstration of a mechanical erosion (spallation) model is also provided to illustrate the unique capabilities of the method.

Ewing, Mark E.↗