Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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 235 records · Page 13

Velocity-split Navier-Stokes solution procedure for incompressible high Reynolds number external flows

A method for solving the Navier-Stokes equations based on splitting the velocity vector into its rotational and irrotational parts was sucessfully applied to internal flow computations. The applicability of the method to external flows is examined by studying several model problems. The model problems are those of laminar and turbulent incompressible flow past a semi-infinite flat plate and laminar incompressible flow past a finite flat plate. For these problems, the procedure accurately reproduces the known solutions and is computationally very efficient even at high Reynolds numbers. Computational aspects of the method are discussed along with the possibility of using the procedure to retrofit a viscous capability into existing potential flow codes.

Dwoyer, D. L.↗

Application of a self-adaptive grid method to complex flows

A directional-split, modular, user-friendly grid point distribution code is applied to several test problems. The code is self-adaptive in the sense that grid point spacing is determined by user-specified constants denoting maximum and minimum grid spacings and constants relating the relative influence of smoothness and orthogonality. Estimates of truncation error, in terms of flow-field gradients and/or geometric features, are used to determine the point distribution. Points are redistributed along grid lines in a specified direction in an elliptic manner over a user-specified subdomain, while orthogonality and smoothness are controlled in a parabolic (marching) manner in the remaining directions. Multidirectional adaption is achieved by sequential application of the method in each coordinate direction. The flow-field solution is redistributed onto the newly distributed grid points after each unidirectional adaption by a simple one-dimensional interpolation scheme. For time-accurate schemes such interpolation is not necessary and time-dependent metrics are carried in the fluid dynamic equations to account for grid movement.

Deiwert, G. S.↗

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements.

Rogers, Stuart E.↗

A multigrid method for steady Euler equations on unstructured adaptive grids

A flux-difference splitting type algorithm is formulated for the steady Euler equations on unstructured grids. The polynomial flux-difference splitting technique is used. A vertex-centered finite volume method is employed on a triangular mesh. The multigrid method is in defect-correction form. A relaxation procedure with a first order accurate inner iteration and a second-order correction performed only on the finest grid, is used. A multi-stage Jacobi relaxation method is employed as a smoother. Since the grid is unstructured a Jacobi type is chosen. The multi-staging is necessary to provide sufficient smoothing properties. The domain is discretized using a Delaunay triangular mesh generator. Three grids with more or less uniform distribution of nodes but with different resolution are generated by successive refinement of the coarsest grid. Nodes of coarser grids appear in the finer grids. The multigrid method is started on these grids. As soon as the residual drops below a threshold value, an adaptive refinement is started. The solution on the adaptively refined grid is accelerated by a multigrid procedure. The coarser multigrid grids are generated by successive coarsening through point removement. The adaption cycle is repeated a few times. Results are given for the transonic flow over a NACA-0012 airfoil.

Riemslagh, Kris↗

A hybridized method for computing high-Reynolds-number hypersonic flow about blunt bodies

A hybridized method for computing the flow about blunt bodies is presented. In this method the flow field is split into its viscid and inviscid parts. The forebody flow field about a parabolic body is computed. For the viscous solution, the Navier-Stokes equations are solved on orthogonal parabolic coordinates using explicit finite differencing. The inviscid flow is determined by using a Moretti type scheme in which the Euler equations are solved, using explicit finite differences, on a nonorthogonal coordinate system which uses the bow shock as an outer boundary. The two solutions are coupled along a common data line and are marched together in time until a converged solution is obtained. Computed results, when compared with experimental and analytical results, indicate the method works well over a wide range of Reynolds numbers and Mach numbers.

Weilmuenster, K. J.↗

Improved Phase-Mask Fabrication of Fiber Bragg Gratings

An improved method of fabrication of Bragg gratings in optical fibers combines the best features of two prior methods: one that involves the use of a phase mask and one that involves interference between the two coherent laser beams. The improved method affords flexibility for tailoring Bragg wavelengths and bandwidths over wide ranges. A Bragg grating in an optical fiber is a periodic longitudinal variation in the index of refraction of the fiber core. The spatial period (Bragg wavelength) is chosen to obtain enhanced reflection of light of a given wavelength that would otherwise propagate relatively unimpeded along the core. Optionally, the spatial period of the index modulation can be made to vary gradually along the grating (such a grating is said to be chirped ) in order to obtain enhanced reflection across a wavelength band, the width of which is determined by the difference between the maximum and minimum Bragg wavelengths. In the present method as in both prior methods, a Bragg grating is formed by exposing an optical fiber to an ultraviolet-light interference field. The Bragg grating coincides with the pattern of exposure of the fiber core to ultraviolet light; in other words, the Bragg grating coincides with the interference fringes. Hence, the problem of tailoring the Bragg wavelength and bandwidth is largely one of tailoring the interference pattern and the placement of the fiber in the interference pattern. In the prior two-beam interferometric method, a single laser beam is split into two beams, which are subsequently recombined to produce an interference pattern at the location of an optical fiber. In the prior phase-mask method, a phase mask is used to diffract a laser beam mainly into two first orders, the interference between which creates the pattern to which an optical fiber is exposed. The prior two-beam interferometric method offers the advantage that the period of the interference pattern can be adjusted to produce gratings over a wide range of Bragg wavelengths, but offers the disadvantage that success depends on precise alignment and high mechanical stability. The prior phase-mask method affords the advantages of compactness of equipment and relative insensitivity to both misalignment and vibration, but does not afford adjustability of the Bragg wavelength. The present method affords both the flexibility of the prior two-beam interferometric method and the compactness and stability of the prior phase-mask method. In this method (see figure), a laser beam propagating along the x axis is normally incident on a phase mask that lies in the (y,z) plane. The phase of light propagating through the mask is modulated with a spatial periodicity, p, along the y axis chosen to diffract the laser light primarily to first order at the angle . (The zero-order laser light propagating along the x axis can be used for alignment and thereafter suppressed during exposure of the fiber.) The diffracted light passes through a concave cylindrical lens, which converts the flat diffracted wave fronts to cylindrical ones, as though the light emanated from a line source. Then two parallel flat mirrors recombine the diffracted beams to form an interference field equivalent to that of two coherent line sources at positions A and B (virtual sources). The interference pattern is a known function of the parameters of the apparatus and of position (x,y) in the interference field. Hence, the tilt, wavelength, and chirp of the Bragg grating can be chosen through suitable adjustments of the apparatus and/or of the position and orientation of the optical fiber. In particular, the Bragg wavelength can be adjusted by moving the fiber along the x axis, and the bandwidth can be modified over a wide range by changing the fiber tilt angle or by moving the phase mask and/or the fiber. Alignment is easy because the zero-order beam defines the x axis. The interference is relatively stable and insensitive to the mechanical vibration because of the gh symmetry and compactness of the apparatus, the fixed positions of the mirrors and lens, and the consequent fixed positions of the two virtual line sources, which are independent of the translations of the phase mask and the laser relative to the lens.

Grant, Joseph↗

Application of a self-adaptive grid method to complex flows

A directional-split modular user-friendly grid-point distribution code is applied to several test problems. Grid-point spacing is determined by user-specified constants denoting maximum and minimum grid spacings and relating the relative influence of smoothness and orthogonality. Estimates of truncation error, in terms of flowfield gradients and/or geometric features, are used to determine the point distribution.

Deiwert, G. S.↗

Calculation of the chordwise load distribution over airfoil sections with plain, split, or serially hinged trailing-edge flaps

A method is presented for the rapid calculation of the incremental chordwise normal-force distribution over an airfoil section due to the deflection of a plain flap or tab, a split flap, or a serially hinged flap. This report is intended as a supplement to NACA Report no. 631, wherein a method is presented for the calculation of the chordwise normal-force distribution over an airfoil without a flap or, as it may be considered, an airfoil with flap (or flaps) neutral. The method enables the determination of the form and magnitude of the incremental normal-force distribution to be made for an airfoil-flap combination for which the section characteristics have been determined. A method is included for the calculation of the flap normal-force and hinge-moment coefficients without necessitating a determination of the normal-force distribution.

Allen, H Julian↗

Choice of implicit and explicit operators for the upwind differencing method

The flux-vector and flux-difference splittings of Steger-Warming, van Leer and Roe are tested in all possible combinations on the implicit and explicit operators that can be distinguished in implicit relaxation methods for the steady Euler and Navier-Stokes equations. The tests include one-dimensional inviscid nozzle flow, and two-dimensional inviscid and viscous shock reflection. Roe's splitting, as anticipated, is found to uniformly yield the most accurate results. On the other hand, an approximate Roe splitting of the implicit operator (the complete Roe splitting is too complicated for practical use) proves to be the least robust with regard to convergence to the steady state. In this respect, the Steger-Warming splitting is the most robust; it leads to convergence when combined with any of the splittings in the explicit operator, although not necessarily in the most efficient way.

Liou, Meng-Sing↗

Choice of implicit and explicit operators for the upwind differencing method

The flux-vector and flux-difference splittings of Steger-Warming, Van Leer and Roe are tested in all possible combinations in the implicit and explicit operators that can be distinguished in implicit relaxation methods for the steady Euler and Navier-Stokes equations. The tests include one-dimensional inviscid nozzle flow, and two-dimensional inviscid and viscous shock reflection. Roe's splitting, as anticipated, is found to uniformly yield the most accurate results. On the other hand, an approximate Roe splitting of the implicit operator (the complete Roe splitting is too complicated for practical use) proves to be the least robust with regard to convergence to the steady state. In this respect, the Steger-Warming splitting is the most robust: it leads to convergence when combined with any of the splittings in the explicit operator, although not necessarily in the most efficient way.

Liou, Meng-Sing↗

The CLASSY clustering algorithm: Description, evaluation, and comparison with the iterative self-organizing clustering system (ISOCLS)

A clustering method, CLASSY, was developed, which alternates maximum likelihood iteration with a procedure for splitting, combining, and eliminating the resulting statistics. The method maximizes the fit of a mixture of normal distributions to the observed first through fourth central moments of the data and produces an estimate of the proportions, means, and covariances in this mixture. The mathematical model which is the basic for CLASSY and the actual operation of the algorithm is described. Data comparing the performances of CLASSY and ISOCLS on simulated and actual LACIE data are presented.

Lennington, R. K.↗

An implicit flux-split algorithm for the compressible Navier-Stokes equations

An implicit upwind scheme for the compressible Navier-Stokes equations is described and applied to the internal flow in a dual-throat nozzle. The method is second-order accurate spatially and naturally dissipative. A spatially-split approximate factorization method is used to obtain efficient steady-state solutions on the NASA Langley VPS-32 (CYBER 205) supercomputer.

Thomas, James L.↗

A viscous flow analysis for the tip vortex generation process

A three dimensional, forward-marching, viscous flow analysis is applied to the tip vortex generation problem. The equations include a streamwise momentum equation, a streamwise vorticity equation, a continuity equation, and a secondary flow stream function equation. The numerical method used combines a consistently split linearized scheme for parabolic equations with a scalar iterative ADI scheme for elliptic equations. The analysis is used to identify the source of the tip vortex generation process, as well as to obtain detailed flow results for a rectangular planform wing immersed in a high Reynolds number free stream at 6 degree incidence.

Shamroth, S. J.↗

A viscous flow analysis of the tip vortex generation process

A viscous primary-secondary flow analysis designed for numerical solution by forward marching integration is applied to the tip vortex generation problem. Equations governing streamwise momentum, streamwise vorticity, and a gas law are solved in conjunction with equations governing scalar and vector surface potentials for the secondary velocities. The numerical method used combines a consistently split linearized block implicit (LBI) scheme for parabolic equations and a scalar iterative ADI scheme for elliptic equations. Computed results are presented for a wing of rectangular planform immersed in a high Reynolds number stream at 6 degrees incidence. The physical mechanism for generation of the tip vortex, as revealed by the analysis and computations, is discussed. Detailed flow results for the case of 6 degree incidence are included.

Shamroth, S. J.↗

A strongly implicit procedure for the compressible Navier-Stokes equations

The strongly implicit procedure of Stone is used to obtain numerical solutions of the compressible Navier-Stokes equations in conservative form. In contrast to the spatially split Douglas-Gunn type methods, the method is shown to be numerically stable for the three-dimensional wave equation. The method is applied to a variety of external and internal two-dimensional flow problems involving shock wave boundary-layer interaction for both laminar and turbulent flows. The results are in good agreement with other methods and/or experiments. The storage penalty associated with the method is discussed and a simple, yet effective, means of minimizing the problem is presented.

Walters, R. W.↗

Computation of turbulent flows on a CAST 10 wing using an upwind scheme

Two and three dimensional turbulent calculations around RAE2822, CAST 10 airfoils and a wing side-wall configuration are presented. A finite volume, cell-center method based on flux difference splitting schemes is used to solve the steady state thin-layer Navier-Stokes equations. Comparisons between the performances of different Riemann solvers and reconstruction schemes are also presented. The influence of a side-wall on the symmetry plane of a cylindrical CAST 10 wing is found to be highly sensitive to the turbulence modeling in the vicinity of the wing side-wall juncture.

Marx, Yves P.↗

Comparative study of computational efficiency of two LU schemes for non-equilibrium reacting flows

The systems of equations governing chemically nonequilibrium flows are solved numerically by two lower-upper(LU) algorithms, namely, LU-SSOR (Yoon and Jameson, 1987) and LU-SW (Steger and Warming, 1981). The methods of flux-Jacobian splitting differentiate the two LU algorithms. The LU-SSOR gives a simpler implicit operator with less temporal damping, while the LU-SW converges faster but requires more computation for each iteration. Relative efficiencies and mutual merits of these two LU algorithms are compared. The results show that LU-SW with upwind-differenced right-hand-side gives the best convergence.

Tsai, Y.-L. Peter↗

Effects of bleed-hole geometry and plenum pressure on three-dimensional shock-wave/boundary-layer/bleed interactions

A numerical study was performed to investigate 3D shock-wave/boundary-layer interactions on a flat plate with bleed through one or more circular holes that vent into a plenum. This study was focused on how bleed-hole geometry and pressure ratio across bleed holes affect the bleed rate and the physics of the flow in the vicinity of the holes. The aspects of the bleed-hole geometry investigated include angle of bleed hole and the number of bleed holes. The plenum/freestream pressure ratios investigated range from 0.3 to 1.7. This study is based on the ensemble-averaged, 'full compressible' Navier-Stokes (N-S) equations closed by the Baldwin-Lomax algebraic turbulence model. Solutions to the ensemble-averaged N-S equations were obtained by an implicit finite-volume method using the partially-split, two-factored algorithm of Steger on an overlapping Chimera grid.

Chyu, Wei J.↗