Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical solution”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Analysis of sound propagation in ducts using the wave envelope concept

A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow for plane wave input. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution does not oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude, and the number of grid points becomes independent of the sound frequency. Physically, the transformed pressure represents the amplitude of the conventional sound wave. Example solutions are presented for sound propagation in a one-dimensional straight hard-wall duct and in a two-dimensional straight soft-wall duct without steady flow. The numerical solutions show evidence of the existence along the duct wall of a developing acoustic pressure diffusion boundary layer which is similar in nature to the conventional viscous flow boundary layer. In order to better illustrate this concept, the wave equation and boundary conditions are written such that the frequency no longer appears explicitly in them. The frequency effects in duct propagation can be visualized solely as an expansion and stretching of the suppressor duct.

Baumeister, K. J.↗
NASA NTRS1977

Turbulent plane Couette flow using probability distribution functions

A numerical scheme employing a combination of the discrete ordinate method and finite differences is developed for solving the one-dimensional form of Lundgren's (1967) model equation for turbulent plane Couette flow. The approach used requires no a priori assumption about the form of the turbulent distribution function, and the numerical solution is obtained directly from the governing differential equations. Two different types of boundary conditions (zero-gradient and Chapman-Enskog) for the distribution function are evaluated by comparing the numerical results with experimental data. It is found that: (1) the present approach gives convergent and stable results over a wide range of Reynolds numbers; (2) Lundgren's equation yields results that compare well with experimental data for mean velocity and skin friction in the case of simple Couette flow; (3) the zero-gradient boundary condition leads to a logarithmic flow profile; and (4) the Chapman-Enskog boundary condition provides very good agreement with experimental data when applied within the near-wall region.

Srinivasan, R.↗
NASA NTRS1989

Radiation-driven winds of hot stars. VI - Analytical solutions for wind models including the finite cone angle effect

Analytical solutions for radiation-driven winds of hot stars including the important finite cone angle effect (see Pauldrach et al., 1986; Friend and Abbott, 1986) are derived which approximate the detailed numerical solutions of the exact wind equation of motion very well. They allow a detailed discussion of the finite cone angle effect and provide for given line force parameters k, alpha, delta definite formulas for mass-loss rate M and terminal velocity v-alpha as function of stellar parameters.

Kudritzki, R. P.↗
NASA NTRS1993

Application of an engineering inviscid-boundary layer method to slender three-dimensional vehicle forebodies

An engineering inviscid-boundary layer method has been modified for application to slender three-dimensional (3-D) forebodies which are characteristic of transatmospheric vehicles. An improved shock description in the nose region has been added to the inviscid technique which allows the calculation of a wider range of body geometries. The modified engineering method is applied to the perfect gas solution over a slender 3-D configuration at angle of attack. The method predicts surface pressures and laminar heating rates on the windward side of the vehicle that compare favorably with numerical solutions of the thin-layer Navier-Stokes equations. These improvements extend the 3-D capabilities of the engineering method and significantly increase its design applications.

Riley, Christopher J.↗
NASA NTRS1993

Application of an Engineering Inviscid-Boundary Layer Method to Slender Three-Dimensional Vehicle Forebodies

An engineering inviscid-boundary layer method has been modified for application to slender three-dimensional (3-D) forebodies which are characteristic of transatmospheric vehicles. An improved shock description in the nose region has been added to the inviscid technique which allows the calculation of a wider range of body geometries. The modified engineering method is applied to the perfect gas solution over a slender 3-D configuration at angle of attack. The method predicts surface pressures and laminar heating rates on the windward side of the vehicle that compare favorably with numerical solutions of the thin-layer Navier-Stokes equations. These improvements extend the 3-D capabilities of the engineering method and significantly increase its design applications.

Riley, Christopher J.↗
NASA NTRS1989

Transonic Navier-Stokes solutions of three-dimensional afterbody flows

The performance of a three-dimensional Navier-Stokes solution technique in predicting the transonic flow past a nonaxisymmetric nozzle was investigated. The investigation was conducted at free-stream Mach numbers ranging from 0.60 to 0.94 and an angle of attack of 0 degrees. The numerical solution procedure employs the three-dimensional, unsteady, Reynolds-averaged Navier-Stokes equations written in strong conservation form, a thin layer assumption, and the Baldwin-Lomax turbulence model. The equations are solved by using the finite-volume principle in conjunction with an approximately factored upwind-biased numerical algorithm. In the numerical procedure, the jet exhaust is represented by a solid sting. Wind-tunnel data with the jet exhaust simulated by high pressure air were also obtained to compare with the numerical calculations.

Compton, William B., III↗
NASA NTRS1982

Numerical analysis of confined turbulent flow

The considered investigation is concerned with the development of an efficient computational method for obtaining a physical understanding of an internal turbulent field. The employed approach makes use of a 'two equation' type model for the turbulence to obtain the numerical solution of a two-dimensional confined turbulent flow. The mean flow governing equations are considered along with the governing equation of the mean temperature and concentrations, and the boundary conditions. The numerical procedure for solving the turbulent flow is discussed, taking into account an approximation to the nonlinear terms, and the inner and outer coupling. Attention is given to a stability convergence analysis, the stability characteristics, and computational examples.

Lin, A.↗
NASA NTRS1984

Numerical simulation of reconnection in an emerging magnetic flux region

Numerical solutions in two dimensions are presented for the resistive MHD equations of an initial boundary value problem, simulating reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. This numerical model displays four main phases, which are interpreted in terms of (1) a slowly evolving quasi-steady phase during which most of the magnetic flux emerges, with approximate equilibrium between magnetic and pressure forces; (2) an impulsive phase in which either the reconnection or continuing emergence of the first phase comes to disrupt the equilibrium, and extensive acceleration occurs as the high pressure region within the emerging region drives the fluid upwards and outwards; (3) a second quasi-steady phase; and (4) a potential-static phase in which continuing reconnection, ohmic dissipation, and fluid transport through boundaries depletes the system of all currents and flows.

Forbes, T. G.↗
NASA NTRS1988

Numerical study of delta wing leading edge blowing

Spanwise and tangential leading edge blowing as a means of controlling the position and strength of the leading edge vortices are studied by numerical solution of the three-dimensional Navier-Stokes equations. The leading edge jet is simulated by defining a permeable boundary, corresponding to the jet slot, where suitable boundary conditions are implemented. Numerical results are shown to compare favorably with experimental measurements. It is found that the use of spanwise leading edge blowing at moderate angle of attack magnifies the size and strength of the leading edge vortices, and moves the vortex cores outboard and upward. The increase in lift primarily comes from the greater nonlinear vortex lift. However, spanwise blowing causes earlier vortex breakdown, thus decreasing the stall angle. The effects of tangential blowing at low to moderate angles of attack tend to reduce the pressure peaks associated with leading edge vortices and to increase the suction peak around the leading edge, so that the integrated value of the surface pressure remains about the same. Tangential leading edge blowing in post-stall conditions is shown to re-establish vortical flow and delay vortex bursting, thus increasing C sub L sub max and stall angle.

Yeh, David↗
NASA NTRS1992

On the numerical treatment of nonlinear source terms in reaction-convection equations

The objectives of this paper are to investigate how various numerical treatments of the nonlinear source term in a model reaction-convection equation can affect the stability of steady-state numerical solutions and to show under what conditions the conventional linearized analysis breaks down. The underlying goal is to provide part of the basic building blocks toward the ultimate goal of constructing suitable numerical schemes for hypersonic reacting flows, combustions and certain turbulence models in compressible Navier-Stokes computations. It can be shown that nonlinear analysis uncovers much of the nonlinear phenomena which linearized analysis is not capable of predicting in a model reaction-convection equation.

Lafon, A.↗
NASA NTRS1992

Acceleration of FDNS flow simulations using initial flowfields generated with a parabolized Navier-Stokes method

The capability of computational fluid dynamics (CFD) analysis to predict complex flowfields has been greatly advanced by the widespread use of flow simulation programs based on the Navier-Stokes (NS) equations. The flow physics are theoretically well represented, and with proper care the numerical solution should not introduce appreciable uncertainties. However, the computational cost of a typical simulation in terms of both memory and execution time are large by current standards. Therefore, the efficiency of the numerical algorithm in solving the set of model equations is one factor determining the usefulness of CFD tools. The objective of the present study is to reduce the execution time of the FDNS flow simulation code by using the parabolized Navier-Stokes (PNS) equations to provide a 'good' starting condition. The technique is not universal, however the PNS model can be applied to convection dominated flows with moderate deflection in geometries that are free of large obstructions.

Tenpas, Peter W.↗
NASA NTRS2021

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗
NASA NTRS2021

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

DAC↗
DOE OSTI2024

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗
NASA NTRS1974

Analysis of high-velocity and real-gas effects on the shock-interference pattern for delta-wing orbiters

Numerical solutions have been generated to study the two-dimensional flow field which results when supersonic flow encounters double-wedge configurations whose angles are such that a Type VI pattern occurs. Two numerical codes were developed: one which used the perfect-gas relations (including the ability to vary gamma to 'simulate' real-gas effects) and a second which incorporated a Mollier table to define equilibrium air properties. Theoretical surface-pressure and heat-transfer distributions in the interaction perturbed region along the simulated 'wing-leading edge' were obtained for velocities from 1167 m/sec to 7610 m/sec.

Bertin, J. J.↗
NASA NTRS1982

The CMC:3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 1: Theoretical

An order-of-magnitude analysis of the subsonic three dimensional steady time averaged Navier-Stokes equations, for semibounded aerodynamic juncture geometries, yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement of first order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. Numerical results document accuracy, convergence, and utility of the developed finite element algorithm, and the CMC:3DPNS computer code applied to an idealized wing-body juncture region. Additional results document accuracy aspects of the algorithm turbulence closure model.

Baker, A. J.↗
NASA NTRS1984

Impact of computers on aerodynamics research and development

Factors motivating the development of computational aerodynamics as a discipline are traced back to the limitations of the tools available to the aerodynamicist before the development of digital computers. Governing equations in exact and approximate forms are discussed together with approaches to their numerical solution. Example results obtained from the successively refined forms of the equations are presented and discussed, both in the context of levels of computer power required and the degree of the effect that their solution has on aerodynamic research and development. Factors pacing advances in computational aerodynamics are identified, including the amount of computational power required to take the next major step in the discipline. Finally, the Numerical Aerodynamic Simulation (NAS) Program - with its 1987 target of achieving a sustained computational rate of 1 billion floating-point operations per second operating on a memory of 240 million words - is briefly discussed in terms of its projected effect on the future of computational aerodynamics.

Peterson, V. L.↗