Variable mesh methods for differential equations
Variable mesh multistep methods for numerical solution of nonlinear differential equations
SEARCH · Search NASA
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.
Variable mesh multistep methods for numerical solution of nonlinear differential equations
Variable mesh multistep predictor-corrector method for iterative solution of ordinary differential equations, considering numerical stability and algorithm efficiency
Performance of stator with wire mesh shell blading during cold air investigation of turbine with transpiration cooled stator blades - Vol. 3
Cold-air tests to determine performance characteristics of single-stage turbine with stator blades employing transpiration coolant ejection through wire mesh shell
Capacitor field calculations for capacitive displacement measurement of open mesh structures such as radio telescopes and radar receivers
Computer program is described which computes and plots coordinates for two-dimensional orthogonal mesh for channel containing solid body, about which flow passes and which spans channel from one wall to the other.
An algorithm is presented to relabel automatically the nodes of an arbitrary finite-element mesh. The purpose of such relabeling is to reduce the bandwidth of the master stiffness matrix produced by the finite-element method. The algorithm uses a random process for the relabeling. Computing time is reduced substantially, compared to systematic methods.
Discontinuous, or weak, solutions of the wave equation, the inviscid form of Burgers equation, and the time-dependent, two-dimensional Euler equations are studied. A numerical method of second-order accuracy in two forms, differential and integral, is used to calculate the weak solutions of these equations for several initial value problems, including supersonic flow past a wedge, a double symmetric wedge, and a sphere. The effect of the computational mesh on the accuracy of computed weak solutions including shock waves and expansion phenomena is studied. Modifications to the finite-difference method are presented which aid in obtaining desired solutions for initial value problems in which the solutions are nonunique.
The computer program DVMESH and the use of the Tektronix DVST graphics terminal were described for applications of preparing mesh data for use in various two-dimensional axisymmetric finite element stress analysis and heat transfer codes.
A model (187 x 187 grid with ten layers) was used to produce a 24 hour forecast using initial conditions for 1200Z, 20 May 1976. This forecast was compared to its five layer counterpart on the 187 x 187 grid and to its 63 x 63 coarse mesh counterpart having the same number of layers. Increases in horizontal resolution lead to significant differences in a one day forecast. Many of these differences represent improvements. Increases in vertical resolution tend to produce smaller impacts on the forecast, except in the region near and above the tropopause. The effect on forecast precipitation is in the 10-20% range, as opposed to the 100% for increases in the horizontal resolution. With respect to model energetics (kinetic energy; square vorticity; square divergence), the model tend to group according to horizontal resolution. The time variations of these parameters show that dynamic initialization is needed to: (1) minimize initialization shock; and (2) stabilize the model context to prevent wash-out of small scale information during the adjustment period (first 6-12 forecast hours).
A new, implicit approximate factorization (AF) algorithm designed to solve the conservative full-potential equation for the transonic flow past arbitrary airfoils has been developed. The new algorithm uses an upwind bias of the density coefficient to provide stability in supersonic regions. This allows the simple two- and three-banded matrix form of the AF scheme to be retained over the entire flow field, even in regions of supersonic flow. A numerical transformation is used to establish an arbitrary body-fitted finite-difference mesh. Airfoil pressure distributions have been computed and are in good agreement with independent results.
The concept of equalizing energy levels was shown to be a viable additional criterion in laying out finite element grids according to the isoenergetic discretization technique. Similar problems specifically with respect to mesh refinement in piecewise approximation theory are being researched. Common criteria in both areas are developed in an effort to cope with the question of discretization for improved piecewise approximations.
The solution of partial differential equations which describe physical phenomena by the use of coordinate transformations is described. The constraints of the problems are stated in geometric terms and include boundary constraints, uniformity constraints, and internal constraints. Algebraic mesh generations are very satisfactory for these constraints.
Several aspects connected with the notion of computation with flow oriented mesh systems are presented. Simple, effective approaches to the ideas discussed are demonstrated in current applications to blown forebody shock layer flow and full bluff body shock layer flow including the massively separated wake region.
A multigrid method for the acceleration of transonic potential flow calculations based on a Galerkin finite element approach is described. In order to allow the use of arbitrary body fitted meshes it is necessary to introduce nonuniform interpolation and residual weighting. Emphasis is put on the construction of these operators consistent with the finite element approximation, while standard successive line overrelaxation is used as a smoothing step. Substantial convergence acceleration is obtained and results are presented for different transonic flow configurations including shocks.
An analysis was conducted for into mesh oil jet lubrication with an arbitrary offset and inclination angle from the pitch point for the case where the oil jet velocity is equal to or less than pitch line velocity. The analysis includes the case for the oil jet offset from the pitch point in the direction of the pinion and where the oil jet is inclined to intersect the common pitch point. Equations were developed for the minimum oil jet velocity required to impinge on the pinion or gear and the optimum oil jet velocity to obtain the maximum impingement depth.
An analysis was conducted for into mesh oil jet lubrication with an arbitrary offset and inclination angle from the pitch point for the case where the oil jet velocity is equal to or greater than gear pitch line velocity. Equations were developed for minimum and maximum oil jet impingement depth. The analysis also included the minimum oil jet velocity required to impinge on the gear or pinion and the optimum oil jet velocity required to obtain the best lubrication condition of maximum impingement depth and gear cooling. It was shown that the optimum oil jet velocity for best lubrication and cooling is when the oil jet velocity equals the gear pitch line velocity. When the oil jet velocity is slightly greater than the pitch line velocity the loaded side of the driven gear and the unloaded side of the pinion receive the best lubrication and cooling with slightly less impingement depth. As the jet velocity becomes much greater than the pitch line velocity the impingement depth is considerably reduced and may completely miss the pinion.
The penetration depth onto the tooth flank of a jet of oil at different velocities pointed at the pitch line on the outgoing side of mesh was determined. The analysis determines the impingement depth for both the gear and the pinion. It includes the cases for speed increasers and decreasers as well as for one to one gear ratio. In some cases the jet will strike the loaded side of the teeth, and in others it will strike the unloaded side of the teeth. In nearly all cases the top land will be cooled regardless of the penetration depth, and postimpingement oil spray will usually provide adequate amounts of oil for lubrication but is marginal or inadequate for cooling.