Search NASA⌕ Search

SEARCH · Search NASA

Results for “equation solving”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

More About the Tetrahedral Unstructured Software System

TetrUSS is a comprehensive suite of computational fluid dynamics (CFD) programs that won the Software of the Year award in 1996 and has found increasing use in government, academia, and industry for solving realistic flow problems (especially in aerodynamics and aeroelastics of aircraft having complex shapes). TetrUSS includes not only programs for solving basic equations of flow but also programs that afford capabilities for efficient generation and utilization of computational grids and for graphical representation of computed flows (see figure). The 2004 version of the Tetrahedral Unstructured Software System (TetrUSS), which is one of two software systems reported in "NASA s 2004 Software of the Year," NASA Tech Briefs, Vol. 28, No. 10 (October 2004), page 18, has been improved greatly since 1996. These improvements include (1) capabilities to simulate viscous flow by solving the Navier-Stokes equations on unstructured grids, (2) portability to personal computers from diverse manufacturers, (3) advanced models of turbulence, (4) a parallel-processing version of one of the unstructured-grid Navier-Stokes-equation-solving programs, and (5) advanced programs for generating unstructured grids.

Abdol-Hamid, Khaled S.↗
NASA NTRS1988

Investigation of viscous/inviscid interaction in transonic flow over airfoils with suction

The viscous/inviscid interaction over transonic airfoils with and without suction is studied. The streamline angle at the edge of the boundary layer is used to couple the viscous and inviscid flows. The potential flow equations are solved for the inviscid flow field. In the shock region, the Euler equations are solved using the method of integral relations. For this, the potential flow solution is used as the initial and boundary conditions. An integral method is used to solve the laminar boundary-layer equations. Since both methods are integral methods, a continuous interaction is allowed between the outer inviscid flow region and the inner viscous flow region. To avoid the Goldstein singularity near the separation point the laminar boundary-layer equations are derived in an inverse form to obtain solution for the flows with small separations. The displacement thickness distribution is specified instead of the usual pressure distribution to solve the boundry-layer equations. The Euler equations are solved for the inviscid flow using the finite volume technique and the coupling is achieved by a surface transpiration model. A method is developed to apply a minimum amount of suction that is required to have an attached flow on the airfoil. The turbulent boundary layer equations are derived using the bi-logarithmic wall law for mass transfer. The results are found to be in good agreement with available experimental data and with the results of other computational methods.

Vemuru, C. S.↗
NASA NTRS1998

Computational Investigation and Validation of Twin-Tail Buffet Response Including Dynamics and Control

Multidisciplinary tools for prediction of single rectangular-tail buffet are extended to single swept-back-tail buffet in transonic-speed flow, and multidisciplinary tools for prediction and control of twin-tail buffet are developed and presented. The configuration model consists of a sharp-edged delta wing with single or twin tails that are oriented normal to the wing surface. The tails are treated as cantilevered beams fixed at the root and allowed to oscillate in both bending and torsion. This complex multidisciplinary problem is solved sequentially using three sets of equations on a dynamic single or multi-block grid structure. The first set is the unsteady, compressible, Reynolds-averaged Navier-Stokes equations which are used for obtaining the flow field vector and the aerodynamic loads on the tails. The Navier-Stokes equations are solved accurately in time using the implicit, upwind, flux-difference splitting, finite volume scheme. The second set is the coupled bending and torsion aeroelastic equations of cantilevered beams which are used for obtaining the bending and torsion deflections of the tails. The aeroelastic equations'are solved accurately in time using, a fifth-order-accurate Runge-Kutta scheme. The third set is the grid-displacement equations and the rigid-body dynamics equations, which are used for updating the grid coordinates due to the tail deflections and rigid-body motions. The tail-buffet phenomenon is predicted for highly-swept, single vertical tail placed at the plane of geometric symmetry, and for highly-swept, vertical twin tails placed at three different spanwise separation distances. The investigation demonstrates the effects of structural inertial coupling and uncoupling of the bending and torsion modes of vibration, spanwise positions of the twin-tail, angle of attack, and pitching and rolling dynamic motions of the configuration model on the tail buffet loading and response. The fundamental issue of twin-tail buffet alleviation is addressed using two active flow-control methods. These methods are the tangential leading-edge blowing and the flow suction from the leading-edge vortex cores along their paths. Qualitative and quantitative comparisons with the available experimental data are presented. The comparisons indicate that the present multidisciplinary aeroelastic analysis tools are robust, accurate and efficient.

Kandil, Osama A.↗
NASA NTRS2001

Investigation of Conjugate Heat Transfer in Turbine Blades and Vanes

We report on work carried out to develop a 3-D coupled Finite Volume/BEM-based temperature forward/flux back (TFFB) coupling algorithm to solve the conjugate heat transfer (CHT) which arises naturally in analysis of systems exposed to a convective environment. Here, heat conduction within a structure is coupled to heat transfer to the external fluid which is convecting heat into or out of the solid structure. There are two basic approaches to solving coupled fluid structural systems. The first is a direct coupling where the solution of the different fields is solved simultaneously in one large set of equations. The second approach is a loose coupling strategy where each set of field equations is solved to provide boundary conditions for the other. The equations are solved in turn until an iterated convergence criterion is met at the fluid-solid interface. The loose coupling strategy is particularly attractive when coupling auxiliary field equations to computational fluid dynamics codes. We adopt the latter method in which the BEM is used to solve heat conduction inside a structure which is exposed to a convective field which in turn is resolved by solving the NASA Glenn compressible Navier-Stokes finite volume code Glenn-HT. The BEM code features constant and bi-linear discontinuous elements and an ILU-preconditioned GMRES iterative solver for the resulting non-symmetric algebraic set arising in the conduction solution. Interface of flux and temperature is enforced at the solid/fluid interface, and a radial-basis function scheme is used to interpolated information between the CFD and BEM surface grids. Additionally, relaxation is implemented in passing the fluxes from the conduction solution to the fluid solution. Results from a simple test example are reported.

Kassab, A. J.↗
NASA NTRS1985

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a given delta wing (Butler wing) at zero incident in a uniform stream. The simulation is done by integrating a set of flow field equations. This set of equations governs the unsteady, viscous, compressible, heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. These equations are solved by the finite difference method, and results obtained for high-speed freestream conditions are compared with theoretical and experimental results. In this study, the Navier-Stokes equations are solved numerically. These equations are unsteady, compressible, viscous, and three-dimensional without neglecting any terms. The time dependency of the governing equations allows the solution to progress naturally for an arbitrary initial initial guess to an asymptotic steady state, if one exists. The equations are transformed from physical coordinates to the computational coordinates, allowing the solution of the governing equations in a rectangular parallel-piped domain. The equations are solved by the MacCormack time-split technique which is vectorized and programmed to run on the CDC VPS 32 computer.

Abolhassani, J. S.↗
DOE OSTI2022

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2009

Method and Apparatus for Predicting Unsteady Pressure and Flow Rate Distribution in a Fluid Network

A method and apparatus for analyzing steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics, external body forces such as gravity and centrifugal force and conjugate heat transfer. In some embodiments, a graphical user interface provides for the interactive development of a fluid network simulation having nodes and branches. In some embodiments, mass, energy, and specific conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. In some embodiments, contained herein are data objects for computing thermodynamic and thermophysical properties for fluids. In some embodiments, the systems of equations describing the fluid network are solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods.

Majumdar, Alok K.↗
NASA NTRS1995

Shape optimization of single- and two-element airfoils on multiblock grids

A multiblock, discrete sensitivity analysis method is used to couple a direct optimization method and a flow analysis method. The domain is divided into smaller subdomains for which the sensitivities are obtained separately. Then, an effective sensitivity equation is solved to complete the coupling of all the sensitivity information. The flow analysis is based on the thin-layer Navier-Stokes equations solved by an implicit, upwind-biased, finite-volume method. The method of feasible directions is used for the present gradient-based optimization approach. First, a transonic airfoil is optimized to investigate the behavior of the method in highly nonlinear flows as well as the effect of different blocking strategies on the procedure. A supercritical airfoil is produced from an initially symmetric airfoil with multiblocking affecting the path but not the final shape. Secondly, a two-element airfoil is shape optimized in subsonic flow to demonstrate the present method's capability of shaping aerodynamically interfering elements simultaneously. For a very low and a very high Reynolds number cases, the shape of the main airfoil and the flap are optimized to yield improved lift-to-drag ratios.

Lacasse, James M.↗
NASA NTRS2004

Numerical Modeling of Conjugate Heat Transfer in Fluid Network

Fluid network modeling with conjugate heat transfer has many applications in Aerospace engineering. In modeling unsteady flow with heat transfer, it is important to know the variation of wall temperature in time and space to calculate heat transfer between solid to fluid. Since wall temperature is a function of flow, a coupled analysis of temperature of solid and fluid is necessary. In cryogenic applications, modeling of conjugate heat transfer is of great importance to correctly predict boil-off rate in propellant tanks and chill down of transfer lines. In TFAWS 2003, the present author delivered a paper to describe a general-purpose computer program, GFSSP (Generalized Fluid System Simulation Program). GFSSP calculates flow distribution in complex flow circuit for compressible/incompressible, with or without heat transfer or phase change in all real fluids or mixtures. The flow circuit constitutes of fluid nodes and branches. The mass, energy and specie conservation equations are solved at the nodes where as momentum conservation equations are solved at the branches. The proposed paper describes the extension of GFSSP to model conjugate heat transfer. The network also includes solid nodes and conductors in addition to fluid nodes and branches. The energy conservation equations for solid nodes solves to determine the temperatures of the solid nodes simultaneously with all conservation equations governing fluid flow. The numerical scheme accounts for conduction, convection and radiation heat transfer. The paper will also describe the applications of the code to predict chill down of cryogenic transfer line and boil-off rate of cryogenic propellant storage tank.

Majumdar, Alok↗
NASA NTRS2004

Generalized Fluid System Simulation Program

A general purpose program implemented on a computer analyzes steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal force. A preprocessor provides for the inter- active development of a fluid network simulation having nodes and branches. Mass, energy, and specie conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. Contained herein are subroutines for computing "real fluid" thermodynamic and thermophysical properties for 12 fluids, and a number of different source options are provided for model- ing momentum sources or sinks in the branches. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. Application and verification of this invention are provided through an example problem, which demonstrates that the predictions of the present invention compare most reasonably with test data.

Alok Kumar Majumdar↗
DOE OSTI2023

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗
NASA NTRS1985

A hybrid numerical technique for predicting the aerodynamic and acoustic fields of advanced turboprops

A hybrid numerical procedure is presented for the prediction of the aerodynamic and acoustic performance of advanced turboprops. A hybrid scheme is proposed which in principle leads to a consistent simultaneous prediction of both fields. In the inner flow a finite difference method, the Approximate-Factorization Alternating-Direction-Implicit (ADI) scheme, is used to solve the nonlinear Euler equations. In the outer flow the linearized acoustic equations are solved via a Boundary-Integral Equation (BIE) method. The two solutions are iteratively matched across a fictitious interface in the flow so as to maintain continuity. At convergence the resulting aerodynamic load prediction will automatically satisfy the appropriate free-field boundary conditions at the edge of the finite difference grid, while the acoustic predictions will reflect the back-reaction of the radiated field on the magnitude of the loading source terms, as well as refractive effects in the inner flow. The equations and logic needed to match the two solutions are developed and the computer program implementing the procedure is described. Unfortunately, no converged solutions were obtained, due to unexpectedly large running times. The reasons for this are discussed and several means to alleviate the situation are suggested.

Homicz, G. F.↗
NASA NTRS2007

Receptivity of Hypersonic Boundary Layers Due to Acoustic Disturbances over Blunt Cone

The transition process induced by the interaction of acoustic disturbances in the free-stream with boundary layers over a 5-degree straight cone and a wedge with blunt tips is numerically investigated at a free-stream Mach number of 6.0. To compute the shock and the interaction of shock with the instability waves the Navier-Stokes equations are solved in axisymmetric coordinates. The governing equations are solved using the 5th -order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. After the mean flow field is computed, acoustic disturbances are introduced at the outer boundary of the computational domain and unsteady simulations are performed. Generation and evolution of instability waves and the receptivity of boundary layer to slow and fast acoustic waves are investigated. The mean flow data are compared with the experimental results. The results show that the instability waves are generated near the leading edge and the non-parallel effects are stronger near the nose region for the flow over the cone than that over a wedge. It is also found that the boundary layer is much more receptive to slow acoustic wave (by almost a factor of 67) as compared to the fast wave.

Kara, K.↗
NASA NTRS1991

Unsteady vortex-dominated flow around wings with oscillating leading-edge flaps

Unsteady vortex-dominated flow around delta wings with oscillating leading-edge flaps represents an important classs of problems for supermaneuverability and flow control of advanced aircraft. The problem is solved using time accurate integration of the unsteady, compressible, thin-layer Navier-Stokes equations in conjunction with the unsteady, linearized, Navier-displacement equations. Starting with an initial configuration of the wing and its flaps, the Navier-Stokes equations are solved on an initial structured grid for the steady flow. The forced oscillation of the flaps is then applied, and the problem is solved accurately in time. The Navier-displacement equations are solved for the grid deformation and the Navier-Stokes equations are solved for the flowfield. Symmetric and anti-symmetric flaps oscillations are presented to study the effect of the flaps oscillation on the leading-edge vortical flow.

Kandil, Osama A.↗
NASA NTRS1998

A Generalized Fluid System Simulation Program to Model Flow Distribution in Fluid Networks

This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.

Majumdar, Alok↗
NASA NTRS1969

Numerical solutions of differential equations

Various numerical methods for solving differential equations were analyzed and refined in an effort to develop a method which was adaptable to a large class of problems. The prime capabilities of the method included accuracy, numerical stability, and economic use of computer time. In multistep processes the corrector was changed at each step.

Wesson, J. R.↗
NASA NTRS1992

Neural network error correction for solving coupled ordinary differential equations

A neural network is presented to learn errors generated by a numerical algorithm for solving coupled nonlinear differential equations. The method is based on using a neural network to correctly learn the error generated by, for example, Runge-Kutta on a model molecular dynamics (MD) problem. The neural network programs used in this study were developed by NASA. Comparisons are made for training the neural network using backpropagation and a new method which was found to converge with fewer iterations. The neural net programs, the MD model and the calculations are discussed.

Shelton, R. O.↗