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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Computational Simulations of a Mach 0.745 Transonic Truss-Braced Wing Design

A joint effort between the NASA Ames and Langley Research Centers was undertaken to analyze the Mach 0.745 variant of the Boeing Transonic Truss-Braced Wing (TTBW) Design. Two different flow solvers, LAVA and USM3D, were used to predict the TTBW flight performance. Sensitivity studies related to mesh resolution and numerical schemes were conducted to define best practices for this type of geometry and flow regime. Validation efforts compared the numerical simulation results of various modeling methods against experimental data taken from the NASA Ames 11-foot Unitary Wind Tunnel experimental data. The fidelity of the computational representation of the wind tunnel experiment, such as utilizing a porous wall boundary condition to model the ventilated test section, was varied to examine how different tunnel effects influence CFD predictions. LAVA and USM3D results both show an approximate 0.5o angle of attack shift from experimental lift curve data. This drove an investigation that revealed that the trailing edge of the experimental model was rounded in comparison to the CAD model, due to manufacturing tolerances, which had not been accounted for in the initial simulations of the experiment. Simulating the TTBW with an approximation of this rounded trailing-edge reduces error by approximately 60%. An accurate representation of the tested TTBW geometry, ideally including any wing twists and deflections experienced during the test under various loading conditions, will be necessary for more thorough validation of the CFD.

Maldonado, Daniel↗
NASA NTRS2002

Development of Computational Aeroacoustics Code for Jet Noise and Flow Prediction

Accurate prediction of jet fan and exhaust plume flow and noise generation and propagation is very important in developing advanced aircraft engines that will pass current and future noise regulations. In jet fan flows as well as exhaust plumes, two major sources of noise are present: large-scale, coherent instabilities and small-scale turbulent eddies. In previous work for the NASA Glenn Research Center, three strategies have been explored in an effort to computationally predict the noise radiation from supersonic jet exhaust plumes. In order from the least expensive computationally to the most expensive computationally, these are: 1) Linearized Euler equations (LEE). 2) Very Large Eddy Simulations (VLES). 3) Large Eddy Simulations (LES). The first method solves the linearized Euler equations (LEE). These equations are obtained by linearizing about a given mean flow and the neglecting viscous effects. In this way, the noise from large-scale instabilities can be found for a given mean flow. The linearized Euler equations are computationally inexpensive, and have produced good noise results for supersonic jets where the large-scale instability noise dominates, as well as for the tone noise from a jet engine blade row. However, these linear equations do not predict the absolute magnitude of the noise; instead, only the relative magnitude is predicted. Also, the predicted disturbances do not modify the mean flow, removing a physical mechanism by which the amplitude of the disturbance may be controlled. Recent research for isolated airfoils' indicates that this may not affect the solution greatly at low frequencies. The second method addresses some of the concerns raised by the LEE method. In this approach, called Very Large Eddy Simulation (VLES), the unsteady Reynolds averaged Navier-Stokes equations are solved directly using a high-accuracy computational aeroacoustics numerical scheme. With the addition of a two-equation turbulence model and the use of a relatively coarse grid, the numerical solution is effectively filtered into a directly calculated mean flow with the small-scale turbulence being modeled, and an unsteady large-scale component that is also being directly calculated. In this way, the unsteady disturbances are calculated in a nonlinear way, with a direct effect on the mean flow. This method is not as fast as the LEE approach, but does have many advantages to recommend it; however, like the LEE approach, only the effect of the largest unsteady structures will be captured. An initial calculation was performed on a supersonic jet exhaust plume, with promising results, but the calculation was hampered by the explicit time marching scheme that was employed. This explicit scheme required a very small time step to resolve the nozzle boundary layer, which caused a long run time. Current work is focused on testing a lower-order implicit time marching method to combat this problem.

Keith, Theo G., Jr.↗
NASA NTRS1984

Using exact solutions to develop an implicit scheme for the baroclinic primitive equations

The exact solutions presently obtained by means of a novel method for nonlinear initial value problems are used in the development of numerical schemes for the computer solution of these problems. The method is applied to a new, fully implicit scheme on a vertical slice of the isentropic baroclinic equations. It was not possible to find a global scale phenomenon that could be simulated by the baroclinic primitive equations on a vertical slice.

Marchesin, D.↗
NASA NTRS1977

Supersonic separated turbulent boundary - layer over a wavy wall

A prediction method is developed for calculating distributions of surface heating rates, pressure and skin friction over a wavy wall in a two-dimensional supersonic flow. Of particular interest is the flow of thick turbulent boundary layers. The surface geometry and the flow conditions considered are such that there exists a strong interaction between the viscous and inviscid flow. First, using the interacting turbulent boundary layer equations, the problem is formulated in physical coordinates and then a reformulation of the governing equations in terms of Levy-Lees variables is given. Next, a numerical scheme for solving interacting boundary layer equations is adapted. A number of modifications which led to the improvement of the numerical algorithm are discussed. Finally, results are presented for flow over a train of up to six waves at various flow conditions.

Polak, A.↗
NASA NTRS1993

Time-accurate Navier-Stokes calculations with multigrid acceleration

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is implemented by modifying the multigrid-multiblock version of the steady Navier-Stokes equations solver, TLNS3D. The scheme is fully implicit in time and uses TLNS3D to iteratively invert the equations at each physical time step. The design objective of the scheme is unconditional stability (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability which is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. An analysis of the iterative procedure and the implementation of this procedure in TLNS3D are discussed. Numerical results are presented to show both the capabilities of the scheme and its speed up relative to the use of global minimum time stepping. Reductions in computational times of an order of magnitude are demonstrated.

Melson, N. Duane↗
NASA NTRS1992

A numerical study of three-dimensional surface tension driven convection with fre surface deformation

The steady three-dimensional thermocapillary motion with a deformable free surface is studied numerically in both normal and zero gravity environments. Flow configurations consist of a square cavity heated from the side. In the analysis, the free surface is allowed to deform and the grid distribution is adapted to the surface deformation. The divergence-free condition is satisfied by using a dual time-stepping approach in the numerical scheme. Convective flux derivatives are evaluated using a third-order accurate upwind-biased flux-split differencing technique. The numerical solutions at the midplane of the square cavity are compared with the results from two-dimensional calculations. In addition, numerial results for cases under zero and normal gravity conditions are compared. Significantly different flow structures and surface deformation have been observed. The comparison of calculated results will be compared with experimental data in the updated version of this paper.

Hsieh, Kwang-Chung↗
NASA NTRS1996

Multigrid Acceleration of Time-Accurate Navier-Stokes Calculations

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is fully implicit in time and is unconditionally stable (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability. This is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. A multigrid-multiblock, steady-state, three-dimensional Navier-Stokes solver, TLNS3D, was modified to iteratively invert the equations at each physical time step. The implementation of this procedure in TLNS3D is discussed. The implications of applying several popular turbulence models to unsteady flow are also considered. Numerical results are presented to show the application of the scheme to various two-dimensional turbulent flows. The results of a three-dimensional laminar flow calculation are also given.

Melson, N. Duane↗
DOE OSTI2021

A consistent and conservative model and its scheme for N -phase- M -component incompressible flows

Here, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Development and application of computational aerothermodynamics flowfield computer codes

Research was performed in the area of computational modeling and application of hypersonic, high-enthalpy, thermo-chemical nonequilibrium flow (Aerothermodynamics) problems. A number of computational fluid dynamic (CFD) codes were developed and applied to simulate high altitude rocket-plume, the Aeroassist Flight Experiment (AFE), hypersonic base flow for planetary probes, the single expansion ramp model (SERN) connected with the National Aerospace Plane, hypersonic drag devices, hypersonic ramp flows, ballistic range models, shock tunnel facility nozzles, transient and steady flows in the shock tunnel facility, arc-jet flows, thermochemical nonequilibrium flows around simple and complex bodies, axisymmetric ionized flows of interest to re-entry, unsteady shock induced combustion phenomena, high enthalpy pulsed facility simulations, and unsteady shock boundary layer interactions in shock tunnels. Computational modeling involved developing appropriate numerical schemes for the flows on interest and developing, applying, and validating appropriate thermochemical processes. As part of improving the accuracy of the numerical predictions, adaptive grid algorithms were explored, and a user-friendly, self-adaptive code (SAGE) was developed. Aerothermodynamic flows of interest included energy transfer due to strong radiation, and a significant level of effort was spent in developing computational codes for calculating radiation and radiation modeling. In addition, computational tools were developed and applied to predict the radiative heat flux and spectra that reach the model surface.

Venkatapathy, Ethiraj↗
NASA NTRS1982

Thermal radiation view factor: Methods, accuracy and computer-aided procedures

The computer aided thermal analysis programs which predicts the result of predetermined acceptable temperature range prior to stationing of these orbiting equipment in various attitudes with respect to the Sun and the Earth was examined. Complexity of the surface geometries suggests the use of numerical schemes for the determination of these viewfactors. Basic definitions and standard methods which form the basis for various digital computer methods and various numerical methods are presented. The physical model and the mathematical methods on which a number of available programs are built are summarized. The strength and the weaknesses of the methods employed, the accuracy of the calculations and the time required for computations are evaluated. The situations where accuracies are important for energy calculations are identified and methods to save computational times are proposed. Guide to best use of the available programs at several centers and the future choices for efficient use of digital computers are included in the recommendations.

Kadaba, P. V.↗
NASA NTRS1983

Application of a viscous-inviscid interaction procedure to predict separated flows with heat transfer

A viscous-inviscid interaction procedure is described for predicting heat transfer in separated flows. The separating flow in a rearward-facing step/asymmetric channel expansion is considered. For viscous regions, the boundary-layer momentum and continuity equations are solved inversely in a coupled manner by a finite-difference numerical scheme. The streamwise convective term is altered to permit marching the solution through regions of reversed flow. The inviscid flow is computed by numerically solving the Laplace equation for stream-function in the region bounded by the displacement surfaces used in the inverse boundary-layer solution. The viscous and inviscid solutions are repeated iteratively until the edge velocities obtained from both solutions are in agreement. Predictions using this method compare favorably with experimental data and other predictions.

Hall, E. J.↗
NASA NTRS2000

Internal Designs Application for Inlet and Nozzle Aeroperformance Improvement

The following research results are based on development of an approach previously proposed by the authors for optimum nozzle design to obtain maximum thrust. The design was denoted a Telescope nozzle. A Telescope nozzle contains one or several internal designs of certain location, which are inserted at certain locations into a divergent conical or planar main nozzle near its exit. Such a design provides additional thrust augmentation over 20% by comparison with the optimum single nozzle of equivalent lateral area. What is more, recent experimental acoustic tests have discovered an essential noise reduction due to Telescope nozzles application. In this paper, some additional theoretical results are presented for Telescope nozzles and a similar approach is applied for aeroperformance improvement of a supersonic inlet. In addition, a classic gas dynamics problem of a similar supersonic flow into a plate has been analyzed. In some particular cases, new exact analytical solutions are obtained for a flow into a wedge with an oblique shock wave. Numerical simulations were conducted for supersonic flow into a divergent portion of a 2D or axisymmetric nozzle with several plane or conuical designs as well as into a 2D or axisymmetric supersonic inlet with a forebody. The 1st order Kryko-Godunov marching numerical scheme for inviscid supersonic flows was used. Several cases were tested using the NASA CFL3d code based on full Navier-Stokes equations. Numerical simulation results have confirmed essential benefits of Telescope design applications in propulsion systems.

Gilinsky, M.↗
NASA NTRS2000

Internal Designs Application for Inlet and Nozzle Aeroperformance Improvement

The following research results are based on development of an approach previously proposed by the authors for optimum nozzle design to obtain maximum thrust. The design was denoted a Telescope nozzle. A Telescope nozzle contains one or several internal designs of certain location, which are inserted at certain locations into a divergent conical or planar main nozzle near its exit. Such a design provides additional thrust augmentation over 20% by comparison with the optimum single nozzle of equivalent lateral area. What is more, recent experimental acoustic tests have discovered an essential noise reduction due to Telescope nozzles application. In this paper, some additional theoretical results are presented for Telescope nozzles and a similar approach is applied for aeroperformance improvement of a supersonic inlet. In addition, a classic gas dynamics problem of a similar supersonic flow into a plate has been analyzed. In some particular cases, new exact analytical solutions are obtained for a flow into a wedge with an oblique shock wave. Numerical simulations were conducted for supersonic flow into a divergent portion of a 2D or axisymmetric nozzle with several plane or conical designs as well as into a 2D or axisymmetric supersonic inlet with a forebody. The 1st order Kryko-Godunov marching numerical scheme for inviscid supersonic flows was used. Several cases were tested using the NASA CFL3d code based on full Navier-Stokes equations. Numerical simulation results have confirmed essential benefits of Telescope design applications in propulsion systems.

Gilinsky, M.↗
NASA NTRS2000

Internal Designs Application for Inlet and Nozzle Aeroperformance Improvement

The following research results are based on development of an approach previously proposed by the authors for optimum nozzle design to obtain maximum thrust. The design was denoted a Telescope nozzle. A Telescope nozzle contains one or several internal designs of certain location, which are inserted at certain locations into a divergent conical or planar main nozzle near its exit. Such a design provides additional thrust augmentation over 20% by comparison with the optimum single nozzle of equivalent lateral area. What is more, recent experimental acoustic tests have discovered an essential noise reduction due to Telescope nozzles application. In this paper, some additional theoretical results are presented for Telescope nozzles and a similar approach is applied for aeroperformance improvement of a supersonic inlet. In addition, a classic gas dynamics problem of a similar supersonic flow into a plate has been analyzed. In some particular cases, new exact analytical solutions are obtained for a flow into a wedge with an oblique shock wave. Numerical simulations were conducted for supersonic flow into a divergent portion of a 2D or axisymmetric nozzle with several plane or conical designs as well as into a 2D or axisymmetric supersonic inlet with a forebody. The 1st order Kryko-Godunov march- ing numerical scheme for inviscid supersonic flows was used. Several cases were tested using the NASA CFL3d code based on full Navier-Stokes equations. Numerical simulation results have confirmed essential benefits of Telescope design applications in propulsion systems.

Gilinsky, M.↗
NASA NTRS1991

Compressible flows with periodic vortical disturbances around lifting airfoils

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method that is presented fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. The numerical scheme is validated through extensive comparisons with known solutions to unsteady vortical flow problems. In general, it is seen that the agreement between the numerical and analytical results is very good for reduced frequencies ranging from 0 to 4, and for Mach numbers ranging from .1 to .8. Numerical results are also presented for a wide variety of flow configurations for the purpose of determining the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. It is seen that each of these parameters can have a significant effect on the unsteady airfoil response to the incident disturbances, and that the effect depends strongly upon the reduced frequency and the dimensionality of the gust. For a one-dimensional (transverse) or two-dimensional (transverse and longitudinal) gust, the results indicate that airfoil thickness increases the unsteady lift and moment at the low reduced frequencies but decreases it at the high reduced frequencies. The results show that an increase in airfoil Mach number leads to a significant increase in the unsteady lift and moment for the low reduced frequencies, but a significant decrease for the high reduced frequencies.

Scott, James R.↗
NASA NTRS2003

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Apte, S. V.↗
NASA NTRS1998

A Second Law Based Unstructured Finite Volume Procedure for Generalized Flow Simulation

An unstructured finite volume procedure has been developed for steady and transient thermo-fluid dynamic analysis of fluid systems and components. The procedure is applicable for a flow network consisting of pipes and various fittings where flow is assumed to be one dimensional. It can also be used to simulate flow in a component by modeling a multi-dimensional flow using the same numerical scheme. The flow domain is discretized into a number of interconnected control volumes located arbitrarily in space. The conservation equations for each control volume account for the transport of mass, momentum and entropy from the neighboring control volumes. In addition, they also include the sources of each conserved variable and time dependent terms. The source term of entropy equation contains entropy generation due to heat transfer and fluid friction. Thermodynamic properties are computed from the equation of state of a real fluid. The system of equations is solved by a hybrid numerical method which is a combination of simultaneous Newton-Raphson and successive substitution schemes. The paper also describes the application and verification of the procedure by comparing its predictions with the analytical and numerical solution of several benchmark problems.

Majumdar, Alok↗
NASA NTRS1991

Navier-Stokes analysis of turbomachinery blade external heat transfer

The two-dimensional, compressible, thin-layer Navier-Stokes and energy equations were solved numerically to obtain heat transfer rates on turbomachinery blades. The Baldwin-Lomax algebraic model and the q-omega low Reynolds number two-equation model were used for modeling of turbulence. For the numerical solution of the governing equations a four-stage Runge-Kutta solver was employed. The turbulence model equations were solved using an implicit scheme. Numerical solutions are presented for two-dimensional flow within two vane cascades. The heat transfer results and the pressure distributions were compared with published experimental data. The agreement between the numerical calculations and the experimental values were found to be generally favorable. The position of transition from laminar to turbulent flow was also predicted accurately.

Gorla, Rama S. R.↗