Engineering topics
Publications and source records attributed to Flores, J..
Report describes study aimed at accelerating rates of convergence of iterative schemes for numerical integration of equations of hypersonic flow of viscous and inviscid fluids. Richardson-type overrelaxation method applied.
Explore the source record for details and available documents.
The convergence of inviscid and viscous hypersonic flow calculations using a two-dimensional flux-splitting code is accelerated by applying a Richardson-type overrelaxation method. Successful results are presented for various cases; and a 50 percent savings in computer time is usually achieved. An analytical formula for the overrelaxation factor is derived, and the performance of this scheme is confirmed numerically. Moreover, application of this overrelaxation scheme produces a favorable preconditioning for Wynn's epsilon-algorithm. Both techniques have been extended to viscous three-dimensional flows and applied to accelerate the convergence of the compressible Navier-Stokes code. A savings of 40 percent in computer time is achieved in this case.
The objective of our reactor design is to supply a lunar-based research facility with 20 MW(e). The fundamental layout of this lunar-based system includes the reactor, power conversion devices, and a radiator. The additional aim of this reactor is a longevity of 12 to 15 years. The reactor is a liquid metal fast breeder that has a breeding ratio very close to 1.0. The geometry of the core is cylindrical. The metallic fuel rods are of beryllium oxide enriched with varying degrees of uranium, with a beryllium core reflector. The liquid metal coolant chosen was natural lithium. After the liquid metal coolant leaves the reactor, it goes directly into the power conversion devices. The power conversion devices are Stirling engines. The heated coolant acts as a hot reservoir to the device. It then enters the radiator to be cooled and reenters the Stirling engine acting as a cold reservoir. The engines' operating fluid is helium, a highly conductive gas. These Stirling engines are hermetically sealed. Although natural lithium produces a lower breeding ratio, it does have a larger temperature range than sodium. It is also corrosive to steel. This is why the container material must be carefully chosen. One option is to use an expensive alloy of cerbium and zirconium. The radiator must be made of a highly conductive material whose melting point temperature is not exceeded in the reactor and whose structural strength can withstand meteor showers.
A fast diagonalized Beam-Warming algorithm is coupled with a zonal approach to solve the three dimensional Euler/Navier-Stokes equations. The computer code, called Transonic Navier-Stokes (TNS), uses a total of four zones for wing configurations (or can be extended to complete aircraft configurations by adding zones). In the inner blocks near the wing surface, the thin-layer Navier-Stokes equations are solved, while in the outer two blocks the Euler equations are solved. The diagonal algorithm yields a speedup of as much as a factor of 40 over the original algorithm/zonal method code. The TNS code, in addition, has the capability to model wind tunnel walls. Transonic viscous solutions are obtained on a 150,000-point mesh for a NACA 0012 wing. A three-order-of-magnitude drop in the L2-norm of the residual requires approximately 500 iterations, which takes about 45 min of CPU time on a Cray-XMP processor. Simulations are also conducted for a different geometrical wing called WING C. All cases show good agreement with experimental data.
A fast diagonalized Beam-Warming algorithm is coupled with a zonal approach to solve the three-dimensional Euler/Navier-Stokes equations. The computer code, called Transonic Navier-Stokes (TNS), uses a total of four zones for wing configurations (or can be extended to complete aircraft configurations by adding zones). In the inner blocks near the wing surface, the thin-layer Navier-Stokes equations are solved, while in the outer two blocks the Euler equations are solved. The diagonal algorithm yields a speedup of as much as a factor of 40 over the original algorithm/zonal method code. The TNS code, in addition, has the capability to model wind tunnel walls. Transonic viscous solutions are obtained on a 150,000-point mesh for a NACA 0012 wing. A three-order-of-magnitude drop in the L2-norm of the residual requires approximately 500 iterations, which takes about 45 min of CPU time on a Cray-XMP processor. Simulations are also conducted for a different geometrical wing called WING C. All cases show good agreement with experimental data.
Photolysis experiments were performed on the H2-CD4-NH3 and He-CD4-NH3 systems. The photolysis (1849 A) involved only NH3. Mixtures of H2:CD4:NH3 included all combinations of the ratios (200,400,800):(10,20,40):4. Two He:CD4:NH3 mixtures were examined where the ratios equalled the combinations 100:(10,20):4. Abstraction of a D from CD4 by the photolytically produced hot hydrogen from ammonia was monitored by mass spectrometric determination of HD. Both experiment and semiempirical hot-atom theory show that H2 is a very poor thermalizer of hot hydrogens with excess kinetic energy of about 2 eV. Applications of the hard-sphere collision model to the H2-CD4-NH3 system resulted in predicted ratios of net HD production to NH3 decomposition that were two orders of magnitude smaller than the experimental ratios. On the other hand, helium is found to be a very efficient thermalizer; here, the classical model yields reasonable agreement with experiments. Application of a semiempirical hot-atom program gave quantitative agreement with experiment for either system.
Transonic flow fields about wing geometries are computed using an Euler/Navier-Stokes approach in which the flow field is divided into several zones. The flow field immediately adjacent to the wing surface is resolved with fine grid zones and solved using a Navier-Stokes algorithm. Flow field regions removed from the wing are resolved with less finely clustered grid zones and are solved with an Euler algorithm. Computational issues associated with this zonal approach, including data base management aspects, are discussed. Solutions are obtained that are in good agreement with experiment, including cases with significant wind tunnel wall effects. Additional cases with significant shock induced separation on the upper wing surface are also presented.
The viscous transonic flow about a multielement airfoil (augmentor-wing) is simulated by coupling full-potential and direct/inverse differential boundary-layer algorithms. Solutions have been obtained for a variety of conditions and are in fair agreement with available experimental data. Typical results from this transonic augmentor-wing code (TAUG-V) require approximately three minutes of CRAY-XMP CPU time. Since this viscous transonic code accounts for most of the important flow physics, yet is still economical, it is a practical tool for the design aerodynamicist.
Transonic flow fields about wing geometries are computed using an Euler/Navier-Stokes approach in which the flow field is divided into several zones. The grid zones immediately adjacent to the wing surface are suitably clustered and solved with the Navier-Stokes equations. Grid zones removed from the wing are less finely clustered and are solved with the Euler equations. Wind tunnel wall effects are easily and accurately modeled with the new grid-zoning algorithm because the wind tunnel grid is constructed as an exact subset of the corresponding free-air grid. Solutions are obtained that are in good agreement with experiment, including cases with significant wind tunnel wall effects and shock-induced separation on the upper wing surface.
Explore the source record for details and available documents.
A fast diagonal algorithm is coupled with a zonal approach to solve the three-dimensional Euler/Navier-Stokes equations. Transonic viscous solutions are obtained on a 150,000 point mesh for a NACA 0012 wing. The new computational approach yields a speedup by as much as a factor of 40 over the standard Beam-Warming algorithm/zonal method originally coded. A three-order-of-magnitude drop in the L2-norm of the residual requires approximately 500 iterations, which takes about 45 min of CPU time on a Cray-XMP. The numerically computed solutions are in good agreement with experimental results. Effects on convergence rate owing to increasing the zonal boundary overlap regions, different stretching distributions in the viscous regions, and different CFL values are also explored.
Previously cited in issue 5, p. 585, Accession no. A83-16678
A quantitative comparison between the Euler and full potential formulations with respect to speed and accuracy is presented. The robustness of the codes used is tested by a number of transonic airfoil cases. The computed results are from four transonic airfoil computer codes. The full potential codes use fully implicit iteration algorithms. The first Euler code uses a fully implicit ADI iteration scheme. The second Euler code uses an explicit Runge Kutta time stepping algorithm which is enhanced by a multigrid convergence acceleration scheme. Quantitative comparisons are made using various plots of lift coefficient versus the average mesh spacing along the airfoil. Besides yielding an asymptotic limit to the lift coefficient, these results also demonstrate the truncation error behavior of the various codes. Quantitative conclusions regarding the full potential and Euler formulations with respect to accuracy, speed, and robustness can be presented.
Results are presented which strengthen the hypothesis of heterogeneity among the carbon- and nitrogen-bearing phases of the Allende meteorite. These data also highlight the possibility of performing physical separations yielding samples in which some of the noble gas- and carbon-bearing phases are extraordinarily predominant over others. The conclusion, based on mass and isotope balance arguments, that a significant portion of the carbonaceous matter in Allende is likely to be gas-poor or gas-free need not weaken the case for carbonaceous carriers for the major noble gas components. The concept that acid-soluble carbonaceous phases contain a multiplicity of components, each of which may have formed under a multiplicity of different physical-chemical conditions, is reemphasized by the results of the present study.
Transonic flow solutions are obtained over a multielement airfoil (augmentor-wing) using the full-potential equation. Solutions obtained for a subcritical case and a strong shock case show good quantitative agreement with experiment in regions not dominated by viscous effects. In those regions where viscous effects are dominant, the results are still in good qualitative agreement. For the strong shock case, Mach number and angle-of-attack corrections were necessary to match experimental coefficient of lift. Typical results from the transonic augmentor-wing Potential Code on the Cray-1S computer require about 10 sec of CPU time for a three-order-of-magnitude drop in the maximum residual. The speed with which solutions can be generated, and the associated low cost, will make this code a practical tool for the design aerodynamicist.
A new spatial differencing scheme for the transonic full-potential equation in conservative form has been developed. This scheme guarantees zero truncation error on any curvilinear mesh for freestream flows in either two- or three-space dimensions. Solutions obtained with this new differencing scheme, away from freestream regions, exhibit greatly improved accuracy, especially for nonsmooth or singular meshes.
Features taken from various models of Titan's atmosphere are combined in a working composite model that provides environmental constraints within which different pathways for organic chemical synthesis are determined. Experimental results and theoretical modeling suggest that the organic chemistry of the satellite is dominated by two processes: photochemistry and energetic particle bombardment. Photochemical reactions of CH4 in the upper atmosphere can account for the presence of C2 hydrocarbons. Reactions initiated at various levels of the atmosphere by cosmic rays, Saturn 'wind', and solar wind particle bombardment of a CH4-N2 atmospheric mixture can account for the UV-visible absorbing stratospheric haze, the reddish appearance of the satellite, and some of the C2 hydrocarbons. In the lower atmosphere photochemical processes will be important if surface temperatures are sufficiently high for gaseous NH3 to exist. It is concluded that the surface of Titan may contain ancient or recent organic matter (or both) produced in the atmosphere.