New implicit boundary procedures - Theory and applications
Previously cited in issue 5, p. 579, Accession no. A83-16536
Rai, M. M.↗
Engineering topics
Rai, M. M.
Publications and source records attributed to Rai, M. M..
Applications of a conservative zonal scheme to transient and geometrically complex problems
A conservative zoning technique, wherein the flow field for a finite-difference calculation is divided into several regions to simplify grid generation, is discussed and is applied in the solution of a two-dimensional problem of complex topology. Calculations are performed on two zonal, or patched, grid systems for the supersonic flow over a double-airfoil configuration. The solution is smooth and continuous across the zonal interfaces, and shock waves pass through the boundaries without distortion. In addition, the time-accuracy of the zonal-boundary method is verified by a two-zone cyclinder calculation with a stationary inner and a rotating outer mesh.The feasibility of the zonal approach for use in the solution of geometrically complex and unsteady problems is thus demonstrated.
Hessenius, K. A.↗
An implicit form for the Osher upwind scheme
Conservative upwind schemes for the Euler equations, such as the Osher scheme, accurately resolve flow discontinuities and correctly model the physics of the problem. However, these schemes require many more arithmetic operations per integration step than simple central-difference schemes and hence result in large computing times. An implicit version of the first-order- and second-order-accurate Osher schemes in two spatial dimensions and generalized coordinates is developed in this study. Because implicit schemes permit the use of large integration steps, in many cases they require fewer integration steps to reach steady-state (especially in calculations on grids with widely varying mesh-cell sizes). The implicit scheme developed in this study accelerated convergence speeds by almost an order of magnitude in the problems considered. Test cases include quasi-one-dimensional nozzle flow and supersonic flow past a cylinder.
Rai, M. M.↗
A conservative treatment of zonal boundaries for Euler equation calculations
Finite-difference calculations require the generation of a grid for the region of interest. A zonal approach, wherein the given region is subdivided into zones and the grid for each zone is generated independently, makes the grid-generation process for complicated topologies and for regions requiring selective grid refinement a fairly simple task. This approach results in new boundaries within the given region, that is, zonal boundaries at the interfaces of the various zones. The zonal-boundary scheme (the integration scheme used to update the points on the zonal boundary) for the Euler equations must be conservative, accurate, stable, and applicable to general curvilinear coordinate systems. A zonal-boundary scheme with these desirable properties is developed in this study. The scheme is designed for explicit, first-order-accurate integration schemes but can be modified to accommodate second-order-accurate explicit and implicit integration schemes. Results for inviscid flow, including supersonic flow over a cylinder, blast-wave diffraction by a ramp, and one-dimensional shock-tube flow are obtained on zonal grids. The conservative nature of the zonal-boundary scheme permits the smooth transition of the discontinuities associated with these flows from one zone to another. The calculations also demonstrate the continuity of contour lines across zonal boundaries that can be achieved with the present zonal scheme.
Rai, M. M.↗
Metric-discontinuous zonal grid calculations using the Osher scheme
Computations on zonal grids - in particular, grids with metric discontinuities resulting from the interspersion of highly clustered regions with coarse regions - are possible using a fully conservative form of the Osher upwind scheme. These zonal grids can result from an abrupt clustering of points near solution discontinuities or near other flow features that require improved resolution. The zonal approach is shown to capture shocks with almost 'shock-fitting' quality but with minimal effort. Results for inviscid flow, including quasi-one-dimensional nozzle flow, supersonic flow over a cylinder, and blast-wave diffraction by a ramp, are presented. These calculations demonstrate the powerful capabilities of the Osher scheme used in conjunction with zonal grids in simulating flow fields with complex shock patterns.
Rai, M. M.↗
Calculation of viscous supersonic flows over finned bodies
The parabolized Navier-Stokes (PNS) equations are used to calculate the viscous, supersonic flow fields about a six-finned projectile and a generic four-finned missile at angles of attack. Since current computer speeds and storage preclude a fully three-dimensional calculation using the unsteady, Reynolds-averaged, Navier-Stokes equations, the applicability of the PNS equations to the above flow fields is of considerable interest. Two important aspects of the calculation are grid generation and the type of smoothing used to prevent nonphysical solutions. This paper includes a description of the grid-generation process. Results in the form of density contours and velocity vector plots are presented for the two configurations. The applicability of the PNS equations to the complicated flow fields considered is successfully demonstrated.
Rai, M. M.↗
New implicit boundary procedures - Theory and applications
Analytical techniques for the application of implicit boundary conditions for inviscid flows to shock and body boundary layer conditions involving the Euler equations are presented. The theory of characteristics is used to update boundary points with spatial second order accuracy. The method is useful for implicit schemes which feature approximate factorization, and as such is incorporated into an existing PNS code. Examples are provided in terms of flows over a cone, over a maneuverable reentry vehicle, and over a finned vehicle. Improvements in the convergence rate are demonstrated for the conical flow solutions.
Rai, M. M.↗
The use of solution adaptive grids in solving partial differential equations
The grid point distribution used in solving a partial differential equation using a numerical method has a substantial influence on the quality of the solution. An adaptive grid which adjusts as the solution changes provides the best results when the number of grid points available for use during the calculation is fixed. Basic concepts used in generating and applying adaptive grids are reviewed in this paper, and examples illustrating applications of these concepts are presented.
Anderson, D. A.↗
Grid evolution in time asymptotic problems
The selection of the proper coordinate system in solving any fluid flow or heat transfer problem is a very important consideration. A new technique of moving mesh points in physical space is introduced so as to reduce the error in a computed asymptotic solution relative to that obtained using a fixed mesh. The technique has been used to solve the simple viscous Burgers' equation in one and two dimensions. Substantial error reductions are demonstrated. The treatment of boundary points and the effect of using different error criteria in generating grids are discussed.
Rai, M. M.↗
Application of adaptive grids to fluid-flow problems with asymptotic solutions
Coordinate system selection is an important consideration in the asymptotic numerical solution of any fluid-flow or heat transfer problem. This paper uses a new technique that provides a simple way of moving the mesh points in physical space in order to reduce the error in the computed asymptotic solution relative to that obtained using a fixed mesh. Applications to fluid-flow problems are presented, including boundary layer flow and inviscid supersonic flow over cylinders, and wedges with associated detached shocks. The treatment of curved boundaries, stationary and nonstationary boundaries, and systems of PDE's is discussed. Significant error reductions are demonstrated.
Rai, M. M.↗
The use of adaptive grids in conjunction with shock capturing methods
The use of shock capturing finite-difference techniques in computing flow fields containing shocks results in a smeared or oscillatory solution in the vicinity of the shocks. This smearing or oscillatory behavior is due to the discretized form of the governing differential equations used to compute the solution. The discretization error can be reduced by a proper clustering of mesh points in the region of the shock and by using shock aligned grids. This paper uses a simple method that was developed earlier to cluster points near the shocks and serves to introduce a new method of generating a shock aligned mesh. Applications to the one-dimensional inviscid Burgers' equation, supersonic flow over a wedge with the associated straight oblique shock, one- and two-dimensional inviscid flows through an expanding duct and the problem of a straight oblique shock in a uniform supersonic freestream are presented. Significant reduction in the oscillatory behavior of the solution is demonstrated.
Rai, M. M.↗
Grid evolution in time asymptotic problems
A technique for generating systems of coordinates for solving time asymptotic problems is described which provides a simple way of moving the mesh points in physical space and reduces the error in the solution relative to that obtained using a fixed mesh. First order partial differential equations are formulated for the grid point velocity in transient problems. Local flow information and boundary motion are used to determine the interior grid point motion.
Rai, M. M.↗