Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence rate”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Multigrid acceleration and turbulence models for computations of 3D turbulent jets in crossflow

A multigrid method is presented for the calculation of three-dimensional turbulent jets in crossflow. Turbulence closure is achieved with either the standard k-epsilon model or a Reynolds Stress Model (RSM). Multigrid acceleration enables convergence rates which are far superior to that for a single grid method. With the k-epsilon model the rate approaches that for laminar flow, but with RSM it is somewhat slower. The increased stiffness of the system of equations in the latter may be responsible. Computed results with both turbulence models are compared with experimental data for a pair of opposed jets in crossflow. Both models yield reasonable agreement with mean flow velocity but RSM yields better prediction of the Reynolds stresses.

Demuren, A. O.↗
NASA NTRS1992

Development of a multigrid transonic potential flow code for cascades

Finite-volume methods for discretizing transonic potential flow equations have proven to be very flexible and accurate for both two and three dimensional problems. Since they only use local properties of the mapping, they allow decoupling of the grid generation from the rest of the problem. A very effective method for solving the discretized equations and converging to a solution is the multigrid-ADI technique. It has been successfully applied to airfoil problems where O type, C type and slit mappings have been used. Convergence rates for these cases are more than an order of magnitude faster than with relaxation techniques. In this report, we describe a method to extend the above methods, with the C type mappings, to airfoil cascade problems.

Steinhoff, John↗
NASA NTRS1990

Acceleration, transport of and radiation by electrons in impulsive phase of flares

Recent work carried out in connection with various manifestations of the impulsive phase of solar flares is reviewed. Characteristics of accelerated electrons are determined, taking into account various possible interactions that electrons undergo as they traverse the flare plasma, and the expected radiation is evaluated via several mechanisms. Constraints are set on the characteristics of the accelerated electrons and on some properties of the flare plasma. Results obtained from comparisons with SMM and other observations of flares during the past solar cycle are described. It is shown that limits can be set on the range and distribution of the accelerated electron spectrum and pitch angle distribution, and on the column depth and field convergence rate of the flaring loop.

Petrosian, Vahe↗
NASA NTRS1992

Improved method for solving the viscous shock layer equations

An improved method for solving the viscous shock layer equations for supersonic/hypersonic flows past blunt-nosed bodies is presented. The method is capable of handling slender to thick bodies. The solution is obtained by solving a coupled set of five equations, built of the four basic viscous shock layer equations and an additional equation for the standoff distance. The coupling of the equations prevents the local iterations divergence problems encountered by previous methods of solution far downstream on slender bodies. It also eliminates the need for local iterations, which were required by previous methods of solution, for a first-order scheme in the streamwise direction. A new global iteration procedure is employed to impose the shock boundary conditions. The procedure prevents the global iteration instability encountered by the basic method of solution and improves the convergence rate of the global iteration procedure of later methods devised to overcome this difficulty. The new technique reduces the computation time by 65-95 percent as compared to previous methods of solution. The method can efficiently be implemented in vector/parallel computers.

Gordon, Rachel↗
NASA NTRS1992

Evaluation of a multigrid-based Navier-Stokes solver for aerothermodynamic computations

A multigrid acceleration technique developed for solving the three-dimensional Navier-Stokes equations is used for computing high Mach number flows over configurations of practical interest. An explicit multistage Runge-Kutta type of time-stepping scheme is used as the basic algorithm. Solutions are presented for a spherically blunted cone at Mach 10 and a modified shuttle orbiter at Mach 6. The computed surface heat-transfer distributions are shown to compare favorably with the experimental data. Effect of grid-refinement on computed heat-transfer distributions is also examined to assess the numerical accuracy of the computed solutions. The rapid convergence rate associated with multigrid schemes in previous applications at transonic speeds is observed at the higher Mach number flows investigated here.

Vatsa, Veer N.↗
NASA NTRS1993

Development of an unstructured solution adaptive method for the quasi-three-dimensional Euler and Navier-Stokes equations

A general solution adaptive scheme based on a remeshing technique is developed for solving the two-dimensional and quasi-three-dimensional Euler and Favre-averaged Navier-Stokes equations. The numerical scheme is formulated on an unstructured triangular mesh utilizing an edge-based pointer system which defines the edge connectivity of the mesh structure. Jameson's four-stage hybrid Runge-Kutta scheme is used to march the solution in time. The convergence rate is enhanced through the use of local time stepping and implicit residual averaging. As the solution evolves, the mesh is regenerated adaptively using flow field information. Mesh adaptation parameters are evaluated such that an estimated local numerical error is equally distributed over the whole domain. For inviscid flows, the present approach generates a complete unstructured triangular mesh using the advancing front method. For turbulent flows, the approach combines a local highly stretched structured triangular mesh in the boundary layer region with an unstructured mesh in the remaining regions to efficiently resolve the important flow features. One-equation and two-equation turbulence models are incorporated into the present unstructured approach. Results are presented for a wide range of flow problems including two-dimensional multi-element airfoils, two-dimensional cascades, and quasi-three-dimensional cascades. This approach is shown to gain flow resolution in the refined regions while achieving a great reduction in the computational effort and storage requirements since solution points are not wasted in regions where they are not required.

Jiang, Yi-Tsann↗
NASA NTRS1992

A conjugate gradients/trust regions algorithms for training multilayer perceptrons for nonlinear mapping

This paper addresses the issue of applying a globally convergent optimization algorithm to the training of multilayer perceptrons, a class of Artificial Neural Networks. The multilayer perceptrons are trained towards the solution of two highly nonlinear problems: (1) signal detection in a multi-user communication network, and (2) solving the inverse kinematics for a robotic manipulator. The research is motivated by the fact that a multilayer perceptron is theoretically capable of approximating any nonlinear function to within a specified accuracy. The algorithm that has been employed in this study combines the merits of two well known optimization algorithms, the Conjugate Gradients and the Trust Regions Algorithms. The performance is compared to a widely used algorithm, the Backpropagation Algorithm, that is basically a gradient-based algorithm, and hence, slow in converging. The performances of the two algorithms are compared with the convergence rate. Furthermore, in the case of the signal detection problem, performances are also benchmarked by the decision boundaries drawn as well as the probability of error obtained in either case.

Madyastha, Raghavendra K.↗
NASA NTRS1993

On waveform multigrid method

Waveform multigrid method is an efficient method for solving certain classes of time dependent PDEs. This paper studies the relationship between this method and the analogous multigrid method for steady-state problems. Using a Fourier-Laplace analysis, practical convergence rate estimates of the waveform multigrid iterations are obtained. Experimental results show that the analysis yields accurate performance prediction.

Taasan, Shlomo↗
NASA NTRS1993

Mapping robust parallel multigrid algorithms to scalable memory architectures

The convergence rate of standard multigrid algorithms degenerates on problems with stretched grids or anisotropic operators. The usual cure for this is the use of line or plane relaxation. However, multigrid algorithms based on line and plane relaxation have limited and awkward parallelism and are quite difficult to map effectively to highly parallel architectures. Newer multigrid algorithms that overcome anisotropy through the use of multiple coarse grids rather than line relaxation are better suited to massively parallel architectures because they require only simple point-relaxation smoothers. The parallel implementation of a V-cycle multiple semi-coarsened grid (MSG) algorithm or distributed-memory architectures such as the Intel iPSC/860 and Paragon computers is addressed. The MSG algorithms provide two levels of parallelism: parallelism within the relaxation or interpolation on each grid and across the grids on each multigrid level. Both levels of parallelism must be exploited to map these algorithms effectively to parallel architectures. A mapping of an MSG algorithm to distributed-memory architectures that demonstrate how both levels of parallelism can be exploited is described. The results is a robust and effective multigrid algorithm for distributed-memory machines.

Overman, Andrea↗
NASA NTRS1992

Multigrid acceleration and turbulence models for computations of 3D turbulent jets in crossflow

A multigrid method is presented for the calculation of three-dimensional turbulent jets in crossflow. Turbulence closure is achieved with either the standard k-epsilon model or a Reynolds stress model (RSM). Multigrid acceleration enables convergence rates which are far superior to that for a single grid method to be obtained with both turbulence models. With the k-epsilon model the rate approaches that for laminar flow, but with RSM it is somewhat slower. The increased stiffness of the system of equation in the latter may be responsible. Computed results with both turbulence models are compared to experimental data for a pair of opposed jets in crossflow. Both models yield reasonable agreement for the mean flow velocity, but RSM yields better predictions of the Reynolds stresses.

Demuren, A. O.↗
NASA NTRS1993

Distributed parallel processing applied to an implicit multigrid Euler/Navier-Stokes algorithm

An implicit multigrid algorithm for the solution of the Euler and Navier-Stokes equations has been implemented within the framework of multiple block-structured grids in which the physical domain is spatially decomposed into several blocks and the solution is advanced in parallel on each block. Utilities have been developed to implement such a scheme in a distributed computing environment. The multi-block algorithm is designed so that the explicit residual calculation is identical to that of single-block scheme, and therefore converged solutions for both schemes must be the same. To accelerate convergence, synchronous and asynchronous multigrid strategies are implemented. Significant speedups have been achieved in a multiple processor environment, while convergence rates similar to those of the single-block scheme are observed.

Tysinger, T. L.↗
NASA NTRS1993

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗
NASA NTRS1993

The application of preconditioning in viscous flows

The present time-derivative preconditioning algorithm is effective in flow conditions ranging from inviscid to very diffusive flows, as well as low subsonic to supersonic flow velocities. By means of a preconditioning matrix that (1) introduces well-conditioned eigenvalues and (2) avoids nonphysical time reversals for viscous flows, a mechanism is obtained which controls the inviscid and viscous time-step parameters at very diffusive flows. These capabilities are demonstrated for a variety of sample problems; convergence rates of solutions that are indistinguishable from those obtained without preconditioning are shown to be accelerated by as much as two orders of magnitude.

Choi, Y.-H.↗
NASA NTRS1993

Laplace's equation and the Dirichlet-Neumann map in multiply connected domains

A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.

Greenbaum, A.↗
NASA NTRS1991

The L sub 1 finite element method for pure convection problems

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

Jiang, Bo-Nan↗
NASA NTRS1993

Characteristics of three-dimensional turbulent jets in crossflow

Three-dimensional turbulent jets in crossflow at low to medium jet-to-crossflow velocity ratios are computed with a finite volume numerical procedure which utilizes a second-moment closure model to approximate the Reynolds stresses. A multigrid method is used to accelerate the convergence rate of the procedure. Comparison of the computations to measured data show good qualitative agreement. All trends are correctly predicted, though there is some uncertainty on the height of penetration of the jet. The evolution of the vorticity field is used to explore the jet-crossflow interaction.

Demuren, A. O.↗
NASA NTRS1993

Quasi-optical grids with thin rectangular patch/aperture elements

Theoretical analysis is presented for an efficient and accurate performance evaluation of quasi-optical grids comprised of thin rectangular patch/aperture elements with/without a dielectric substrate/superstrate. The convergence rate of this efficient technique is improved by an order of magnitude with the approximate edge conditions incorporated in the basis functions of the integral equation solution. Also presented are the interesting applications of this efficient analytical technique to the design and performance evaluation of the coupling grids and beam splitters in the optical systems as well as thermal protection sunshields used in the communication systems of satellites and spacecrafts.

Wu, Te-Kao↗
NASA NTRS1993

Implicit multigrid Euler solutions with symmetric Total-Variation-Diminishing dissipation

A symmetric Total-Variation-Diminishing (TVD) formulation of the numerical dissipation terms has been incorporated into a diagonalized alternating direction implicit multigrid algorithm to solve the Euler equations of inviscid compressible flow. The new treatment of the dissipation makes is possible to capture both very strong and very weak shocks, virtually without oscillation for the steady flows of interest here. In addition, the TVD constraint fixes one of the two previously arbitrary constants in the formulation of the dissipation, and results in both converged solutions and convergence rates which are relatively insensitive to the choice of the remaining dissipation parameter.

Caughey, David A.↗