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 415 records · Page 23

Asynchronous numerical solution of PDEs on parallel computers

Asynchronous and partially-synchronous finite difference schemes are investigated for the one-dimensional steady-state and nonsteady-state heat equation, respectively. It is shown experimentally that the asynchronous scheme runs twice as fast as the parallel synchronous scheme. For relatively small values of time levels, however, the asynchronous scheme produces the least accurate results, while the partially-asynchronous scheme is almost as accurate as its synchronous counterpart.

Amitai, Dganit↗

Projection methods for the numerical solution of Markov chain models

Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.

Saad, Youcef↗

Numerical solution of a three-dimensional cubic cavity flow by using the Boltzmann equation

A three-dimensional cubic cavity flow has been analyzed for diatomic gases by using the Boltzmann equation with the Bhatnagar-Gross-Krook (B-G-K) model. The method of discrete ordinate was applied, and the diffuse reflection boundary condition was assumed. The results, which show a consistent trend toward the Navier-Stokes solution as the Knudson number is reduced, give us confidence to apply the method to a three-dimensional geometry for practical predictions of rarefied-flow characteristics. The CPU time and the main memory required for a three-dimensional geometry using this method seem reasonable.

Hwang, Danny P.↗

Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes

A new combination of a finite volume discretization in conjunction with carefully designed dissipative terms of third order, and a Runge Kutta time stepping scheme, is shown to yield an effective method for solving the Euler equations in arbitrary geometric domains. The method has been used to determine the steady transonic flow past an airfoil using an O mesh. Convergence to a steady state is accelerated by the use of a variable time step determined by the local Courant member, and the introduction of a forcing term proportional to the difference between the local total enthalpy and its free stream value.

Jameson, A.↗

Numerical solutions for heat flow in adhesive lap joints

The present formulation for the modeling of heat transfer in thin, adhesively bonded lap joints precludes difficulties associated with large aspect ratio grids required by standard FEM formulations. This quasi-static formulation also reduces the problem dimensionality (by one), thereby minimizing computational requirements. The solutions obtained are found to be in good agreement with both analytical solutions and solutions from standard FEM programs. The approach is noted to yield a more accurate representation of heat-flux changes between layers due to a disbond.

Howell, P. A.↗

Electromagnetically driven relativistic jets - A class of self-consistent numerical solutions

We explore the generalized Grad-Shafranov (GS) equation (Lovelace et al. 1986), which describes the cross-field balance of the magnetic flux surfaces in a stationary, axisymmetric, cold relativistic MHD wind. We construct a family of self-consistent, jet-type, non-self-similar solutions in which flux surfaces thread the equatorial plane vertically and are eventually collimated into nested cylinders. These solutions conserve the total energy and angular momentum along each flux surface and satisfy the relativistic GS equation in the cross-field direction. We find that the total specific energy in the flow depends primarily on a dimensionless combination of the poloidal magnetic flux, the rotation, and the mass loading rate, which reduces to Michel's (1969) magnetization parameter for radial winds. The final width of the flow and its ratio of kinetic energy flux to Poynting flux depend sensitively on the pressure distribution of the ambient medium, which is required to contain the jet at its outer edge when no jet-confirming external toroidal magnetic field is present. Furthermore, we show that the magnetic pressure in the flow can be much higher near the axis than close to the edge, reflecting the 'pinching' effect of the toroidal magnetic field.

Li, Zhi-Yun↗

Recommendations for numerical solution of reinforced-panel and fuselage-ring problems

Procedures are recommended for solving the equations of equilibrium of reinforced panels and isolated fuselage rings as represented by the external loads and the operations table established according to Southwell's method. From the solution of these equations the stress distribution can be easily determined. The method of systematic relaxations, the matrix-calculus method, and several other methods applicable in special cases are discussed. Definite recommendations are made for obtaining the solution of reinforced-panel problems which are generally designated as shear lag problems. The procedures recommended are demonstrated in the analysis of a number of panels. In the case of fuselage rings it is not possible to make definite recommendations for the solution of the equilibrium equations for all rings and loadings. However, suggestions based on the latest experience are made and demonstrated on several rings.

Hoff, N J↗

Numerical solutions of the complete Navier-Stokes equations

The objective of this study is to compare the use of assumed pdf (probability density function) approaches for modeling supersonic turbulent reacting flowfields with the more elaborate approach where the pdf evolution equation is solved. Assumed pdf approaches for averaging the chemical source terms require modest increases in CPU time typically of the order of 20 percent above treating the source terms as 'laminar.' However, it is difficult to assume a form for these pdf's a priori that correctly mimics the behavior of the actual pdf governing the flow. Solving the evolution equation for the pdf is a theoretically sound approach, but because of the large dimensionality of this function, its solution requires a Monte Carlo method which is computationally expensive and slow to coverage. Preliminary results show both pdf approaches to yield similar solutions for the mean flow variables.

Hassan, H. A.↗

Second-order numerical solution of time-dependent, first-order hyperbolic equations

A finite difference scheme is developed to find an approximate solution of two similar hyperbolic equations, namely a first-order plane wave and spherical wave problem. Finite difference approximations are made for both the space and time derivatives. The result is a conditionally stable equation yielding an exact solution when the Courant number is set to one.

Shah, Patricia L.↗

Numerical solution of differential equations by artificial neural networks

Conventionally programmed digital computers can process numbers with great speed and precision, but do not easily recognize patterns or imprecise or contradictory data. Instead of being programmed in the conventional sense, artificial neural networks (ANN's) are capable of self-learning through exposure to repeated examples. However, the training of an ANN can be a time consuming and unpredictable process. A general method is being developed by the author to mate the adaptability of the ANN with the speed and precision of the digital computer. This method has been successful in building feedforward networks that can approximate functions and their partial derivatives from examples in a single iteration. The general method also allows the formation of feedforward networks that can approximate the solution to nonlinear ordinary and partial differential equations to desired accuracy without the need of examples. It is believed that continued research will produce artificial neural networks that can be used with confidence in practical scientific computing and engineering applications.

Meade, Andrew J., Jr.↗

Effects of the computational time step on numerical solutions for turbulent flow

Effects of large computational time steps on the computed turbulence were investigated using a fully implicit method. In turbulent channel flow computations the largest computational time step in wall units which led to accurate prediction of turbulence statistics was determined. Turbulence fluctuations could not be sustained if the computational time step was near or larger than the Kolmogorov time scale.

Choi, Haecheon↗

Numerical solutions of the complete Navier-Stokes equations

A compressible turbulence model based on the kappa-zeta (enstrophy) model of Robinson et al. is developed. The model incorporates a new compressibility correction model that is quite different from previous models. Results are compared to the mixing layer data compiled by Settles and Dodson. It is shown that the predictions of the theory are in good agreement with experiment. Moreover, the present model is superior to the kappa-epsilon model with traditional compressibility corrections.

Hassan, H. A.↗