Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite difference scheme”

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 523 records · Page 29

A comparison of two formulations for high-order accurate essentially non-oscillatory schemes

The finite-volume and finite-difference implementations of high-order accurate essentially nonoscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted.

Casper, J.↗

Prediction of asymmetric vortical flows around slender bodies using Navier-Stokes equations

Steady and unsteady asymmetric vortical flows around slender bodies at high angles of attack are solved using the unsteady, compressible, thin-layer Navier-Stokes equations. An implicit, upwind-biased, flux-difference splitting, finite-volume scheme is used for the numerical computations. For supersonic flows past point cones, the locally conical flow assumption have been used for efficient computational studies of this phenomenon. Asymmetric flows past a 5-deg semiapex-angle circular cone at different angles of attack, free-stream Mach numbers, and Reynolds numbers have been studied in responses to different sources of disturbances. The effects of grid fineness and computational domain size have also been investigated. Next, the responses of three-dimensional supersonic asymmetric flow around a 5-deg circular cone at different angles of attack and Reynolds numbers to short-duration sideslip disturbances are presented. The results show that flow asymmetry becomes stronger as the Reynolds number and angles of attack are increased. One of the cases of flow over a cone-cylinder configuration is validated fairly well by experimental data.

Liu, C. H.↗

Effect of Spatial Filtering in Implicit Large-Eddy Simulations of Separated Flows

The relatively high Reynolds number of turbulent flows encountered in various applications puts these problems well beyond the reach of direct numerical simulation (DNS) at present. Meanwhile, lower-fidelity Reynolds-averaged Navier-Stokes (RANS) calculations are known to be not accurate enough in complex problems, such as smooth-body flow separation and other flows involving highly-unsteady phenomena. Hence, given the current infeasibility of DNS and the unsatisfactory performance of RANS, intermediate techniques such as large-eddy simulation (LES) and hybrid RANS-LES, whose fidelity lie between RANS and DNS, have received much attention for application to various problems of practical importance. Modeling of the effect of missing scales on resolved scales, also known as subgrid-scale (SGS) modeling, is an important subject for LES. SGS models can be broadly categorized as explicit or implicit approaches. The explicit approach is based on an SGS model that explicitly appears in the governing equations expressed in the form of so-called “filtered Navier-Stokes equations”, which describe the evolution of the turbulence scales resolved by the LES grid. The effect of the scales unresolved by the grid is represented by the SGS model. The implicit modeling approach, on the other hand, does not employ an explicit model but instead treats the intrinsic dissipation of the numerical discretization scheme as an implicit SGS model. An LES without an explicit SGS model is commonly termed as an implicit LES (ILES). The relative merits of one SGS modeling approach over another is a subject of ongoing debate. We have opted to employ an ILES methodology, based on high-order compact finite-difference and spatial filtering schemes, in our recent investigations of separated flow problems [1, 2]. Further discussion of our preference of ILES over explicit LES is provided in Uzun and Malik [2]. The spatial filtering operation, described in the next section, is treated as an implicit SGS model for the ILES. Some observations made during the course of our recent investigations, which pointed out to excessive numerical dissipation in certain parts of the flowfield, prompted us to take a closer look at the potential effect of the spatial filter on ILES predictions. This technical note is therefore devoted to spatial filter effects in the context of a high Reynolds number, transonic shock-induced separated flow.

Ali Uzun↗

A comparison of two formulations for high-order accurate essentially non-oscillatory schemes

The finite-volume and finite-difference implementations of high-order accurate essentially non-oscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted.

Casper, Jay↗

Recent advances in numerical simulation and control of asymmetric flows around slender bodies

The problems of asymmetric flow around slender bodies and its control are formulated using the unsteady, compressible, thin-layer or full Navier-Stokes equations which are solved using an implicit, flux-difference splitting, finite-volume scheme. The problem is numerically simulated for both locally-conical and three-dimensional flows. The numerical applications include studies of the effects of relative incidence, Mach number and Reynolds number on the flow asymmetry. For the control of flow asymmetry, the numerical simulation cover passive and active control methods. For the passive control, the effectiveness of vertical fins placed in the leeward plane of geometric symmetry and side strakes with different orientations is studied. For the active control, the effectiveness of normal and tangential flow injection and surface heating and a combination of these methods is studied.

Kandil, Osama A.↗

Supercomputing of supersonic flows using upwind relaxation and MacCormack schemes

The performance of two numerical solution schemes, (1) an implicit upwind relaxation with a finite-volume discretization (Thomas and Walters, 1985) and (2) an explicit-implicit MacCormack (1981) scheme with a finite-difference discretization, is compared in two-dimensional simulations of supersonic flow past a flat plate with leading edge, a rearward-facing step, a 10-deg compression corner, a NACA 0012 airfoil at high angle of attack, and a cavity. The algebraic turbulence model, the solution methods, and the boundary conditions and SIMD coding are explained, and the results are presented in tables and graphs and characterized with reference to published experimental data. Scheme (1) is found to converge more rapidly and to give more accurate results than (2) in a wide range of problem types.

Baysal, O.↗

Active control of asymmetric conical flow using spinning and rotatory oscillations

The effectiveness of active control on asymmetric flows around circular cones is investigated computationally using cone spinning and rotatory oscillation around its axis. The investigation uses the time-accurate solution of the unsteady, compressible, full Navier-Stokes equations with the implicit, upwind, flux-difference splitting, finite-volume scheme. The present solutions are obtained under the locally-conical-flow assumption in order to understand the flow physics using very fine grids for reasonable flow resolution at low computational cost. For all the computational solutions, a grid of 241 x 81 x 2 points in the wrap-around, normal and axial directions, respectively, is used. The grid is spinning or oscillating rigidly with the cone according to its motion and the kinematical and dynamical boundary conditions are modified accordingly. The computational applications include the effects of uniform spinning rates and periodic rotatory oscillations at different amplitudes and frequencies on the flow asymmetry.

Kandil, Osama A.↗

Applications of an adaptive unstructured solution algorithm to the analysis of high speed flows

An upwind cell-centered scheme for the solution of steady laminar viscous high-speed flows is implemented on unstructured two-dimensional meshes. The first-order implementation employs Roe's (1981) approximate Riemann solver, and a higher-order extension is produced by using linear reconstruction with limiting. The procedure is applied to the solution of inviscid subsonic flow over an airfoil, inviscid supersonic flow past a cylinder, and viscous hypersonic flow past a double ellipse. A detailed study is then made of a hypersonic laminar viscous flow on a 24-deg compression corner. It is shown that good agreement is achieved with previous predictions using finite-difference and finite-volume schemes. However, these predictions do not agree with experimental observations. With refinement of the structured grid at the leading edge, good agreement with experimental observations for the distributions of wall pressure, heating rate and skin friction is obtained.

Thareja, R. R.↗

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗

Combustion of hydrogen injected into a supersonic airstream (a guide to the HISS computer program)

A computer program based on a finite-difference, implicit numerical integration scheme is described for the prediction of hydrogen injected into a supersonic airstream at an angle ranging from normal to parallel to the airstream main flow direction. Results of calculations for flow and thermal property distributions were compared with 'cold flow data' taken by NASA/Langley and show excellent correlation. Typical results for equilibrium combustion are presented and exhibit qualitatively plausible behavior. Computer time required for a given case is approximately one minute on a CDC 7600. A discussion of the assumption of parabolic flow in the injection region is given which demonstrates that improvement in calculation in this region could be obtained by a partially-parabolic procedure which has been developed. It is concluded that the technique described provides an efficient and reliable means for analyzing hydrogen injection into supersonic airstreams and the subsequent combustion.

Dyer, D. F.↗

Equatorial wave-mean flow interaction - A numerical study of the role of latitudinal shear

A time-dependent primitive equation model for an equatorial channel is used to assess the interaction of equatorial Kelvin and mixed Rossby-gravity waves with the mean flow. The proposed model involves a semiimplicit time-differencing scheme and a finite-difference grid in the meridional plane. It is shown that forced equatorial waves interact with mean flow to produce equatorial jets characterized by downward-moving westerly (Kelvin wave forcing) and easterly (mixed Rossby-gravity wave forcing) shear zones, respectively. For parameters characteristic of the observed waves in the equatorial stratosphere, the wave-mean flow interaction process always reduces the amplitude of any initial cross-equatorial mean wind shear. The mean flow profile tends to become symmetric about the equator as the interaction process continues.

Holton, J. R.↗

Navier-Stokes computations of symmetric and asymmetric vortex shedding around slender bodies

The unsteady, compressible Reynolds-averaged thin-layer Navier-Stokes equations are used to solve for symmetric and asymmetric vortical flows around slender pointed bodies of revolution. The modified Baldwin and Lomax algebraic two-layer turbulent model is used for the eddy viscosity calculation in prescribed turbulent regions. The implicit, upwind flux-difference splitting finite-volume scheme is used to obtain the solutions. Computational results are presented for a low-speed symmetrical vortical flow around a 3.5-caliber tangent-ogive cylinder for fully laminar flow. Computed results are presented for a high-speed asymmetric vortex-shedding flow around a cone at 3.286 relative incidence, using a fine computational grid. The asymmetric vortex-shedding flow has been produced through a small asymmetric transitional perturbation in the cross-flow plane.

Kandil, Osama A.↗

Integral solution of unsteady full-potential equation for a transonic pitching airfoil

The unsteady full-potential equation formulation in a moving frame of reference has been has been developed and used to solve unsteady transonic flow problems. An unsteady integral-equation shock-capturing (IE-SC) scheme has been developed. The resulting unsteady IE-SC scheme is applied to a NACA 0012 airfoil undergoing a pitching oscillation. The numerical results are compared with those of an implicit, approximately factored, finite-volume Euler scheme. The present scheme is efficient in terms of the number of iterations as compared to the other existing schemes, which use finite-difference or finite-volume methods.

Kandil, Osama A.↗

Direct Finite-Difference Simulations Of Turbulent Flow

Report discusses use of upwind-biased finite-difference numerical-integration scheme to simulate evolution of small disturbances and fully developed turbulence in three-dimensional flow of viscous, incompressible fluid in channel. Involves use of computational grid sufficiently fine to resolve motion of fluid at all relevant length scales.

Rai, Man Mohan↗

Asymmetric supersonic flow around cones with noncircular sections

The unsteady, compressible, single and double thin-layer, Navier-Stokes equations are used to solve for steady and unsteady, asymmetric, supersonic flow around pointed bodies with noncircular sections at high incidence and zero side slip. The equations are solved by using an implicit, upwind, flux-difference splitting finite-volume scheme. Since the flow is locally conical, the solutions are presented on a cross-flow plane at the axial station of unit. The grid is generated by using a modified Joukowski transformation. The computational applications cover noncircular sections with elliptic and diamond shapes. Unsteady asymmetric vortex shedding has been captured at large angles of attack. It is shown that for the same flow conditions and same cross-section fineness ratio, the diamond-section cones develop less flow asymmetry than the elliptic-section cones. Passive control of flow asymmetry was also demonstrated computationally.

Kandil, Osama A.↗

Computation of steady and unsteady compressible quasi-axisymmetric vortex flow and breakdown

The unsteady, compressible Navier-Stokes equations are used to compute and analyze compressible quasi-axisymmetric isolated vortices. The Navier-Stokes equations are solved using an implicit, upwind, flux-difference splitting finite-volume scheme. The developed three-dimensional solver has been verified by comparing its solution profiles with those of a slender, quasi-axisymmetric vortex solver for a subsonic, isolated quasi-axisymmetric vortex in an unbounded domain. The Navier-Stokes solver is then used to solve for a supersonic quasi-axisymmetric vortex flow in a configured circular duct. Steady and unsteady vortex-shock interactions and breakdown have been captured. The problem has also been calculated using the Euler solver of the same code and the results are compared with those of the Navier-Stokes solver. The effect of the initial swirl has been tentatively studied.

Kandil, Osama A.↗

Parallelized solvers for heat conduction formulations

Based on multilevel partitioning, this paper develops a structural parallelizable solution methodology that enables a significant reduction in computational effort and memory requirements for very large scale linear and nonlinear steady and transient thermal (heat conduction) models. Due to the generality of the formulation of the scheme, both finite element and finite difference simulations can be treated. Diverse model topologies can thus be handled, including both simply and multiply connected (branched/perforated) geometries. To verify the methodology, analytical and numerical benchmark trends are verified in both sequential and parallel computer environments.

Padovan, Joe↗

Numerical simulation of steady and unsteady asymmetric vortical flow

The unsteady, compressible, thin-layer, Navier-Stokes equations are solved to simulate steady and unsteady, asymmetric, vortical laminar flow around cones at high incidences and supersonic Mach numbers. The equations are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme. Locally conical flows are assumed and the solutions are obtained by forcing the conserved components of the flowfield vector to be equal at two axial stations located at 0.95 and 1.0. Computational examples cover steady and unsteady asymmetric flows around a circular cone and its control using side strakes. Unsteady asymmetric flows are also presented for elliptic- and diamond-section cones, which model asymmetric vortex shedding around round and sharp-edged delta wing flows.

Kandil, O. A.↗