Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical mixing”

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 127 records · Page 7

Improved methods of vibration analysis of pretwisted, airfoil blades

Vibration analysis of pretwisted blades of asymmetric airfoil cross section is performed by using two mixed variational approaches. Numerical results obtained from these two methods are compared to those obtained from an improved finite difference method and also to those given by the ordinary finite difference method. The relative merits, convergence properties and accuracies of all four methods are studied and discussed. The effects of asymmetry and pretwist on natural frequencies and mode shapes are investigated. The improved finite difference method is shown to be far superior to the conventional finite difference method in several respects. Close lower bound solutions are provided by the improved finite difference method for untwisted blades with a relatively coarse mesh while the mixed methods have not indicated any specific bound.

Subrahmanyam, K. B.↗

Fully-coupled analysis of jet mixing problems. Three-dimensional PNS model, SCIP3D

Numerical procedures formulated for the analysis of 3D jet mixing problems, as incorporated in the computer model, SCIP3D, are described. The overall methodology closely parallels that developed in the earlier 2D axisymmetric jet mixing model, SCIPVIS. SCIP3D integrates the 3D parabolized Navier-Stokes (PNS) jet mixing equations, cast in mapped cartesian or cylindrical coordinates, employing the explicit MacCormack Algorithm. A pressure split variant of this algorithm is employed in subsonic regions with a sublayer approximation utilized for treating the streamwise pressure component. SCIP3D contains both the ks and kW turbulence models, and employs a two component mixture approach to treat jet exhausts of arbitrary composition. Specialized grid procedures are used to adjust the grid growth in accordance with the growth of the jet, including a hybrid cartesian/cylindrical grid procedure for rectangular jets which moves the hybrid coordinate origin towards the flow origin as the jet transitions from a rectangular to circular shape. Numerous calculations are presented for rectangular mixing problems, as well as for a variety of basic unit problems exhibiting overall capabilities of SCIP3D.

Wolf, D. E.↗

Remote detection of ocean features in the Lesser Antilles using ERTS-1 data

The author has identified the following significant results. Photographic data received from the ERTS-1 satellite over the Lesser Antilles Islands shows distinct ocean features on the leeward side of each island. Attempts to relate these features to ocean eddy formations with the aid of ground truth data proved unsuccessful. However, analysis of surface and upper air wind data correlate extremely well with the size, shape, and downwind extent of the ocean features. Studies to date indicate strongly that these features result from horizontal differences in sea surface roughness due to the wind shadow effect of the islands. The results suggest that horizontal variations in the reflectance of the sea surface will make remote sensing of the ocean mixed layer more difficult than previously anticipated. The surface reflection seems to be large enough to mask the smaller variations in backscattered energy from the mixed layer. Efforts to limit the effect of surface reflectance by photographic differencing of two MSS bands were unsuccessful. A supplementary study to enhance the energy reflected from the mixed layer through numerical differencing of ERTS-1 digital data is planned.

Hanson, K. J.↗

Magnetohydrodynamic and gasdynamic aspects of solar-wind flow around terrestrial planets: A critical review

Problems in space physics are discussed whose models, in simplified form, reduce to a supersonic flow scheme with a detached shock wave, namely: (1) solar wind interaction with an intrinsic planetary magnetic field; (2) solar wind interaction with the ionized component of the atmosphere of a comet; and (3) solar wind interaction with the ionosphere of a planet which does not possess its own magnetic field. The numerical study of the above problems is performed with the use of magnetogasdynamic equations for an ideal single fluid model. From the physical veiwpoint, the problems are solved in terms of as simple phenomena as possible; the principal objective is to make recently developed methods of numerical analysis of mixed flows applicable to space physics problems.

Spreiter, J. R.↗

Error analysis of finite difference schemes applied to hyperbolic initial boundary value problems

Finite difference methods for the numerical solution of mixed initial boundary value problems for hyperbolic equations are studied. The reported investigation has the objective to develop a technique for the total error analysis of a finite difference scheme, taking into account initial approximations, boundary conditions, and interior approximation. Attention is given to the Cauchy problem and the initial approximation, the homogeneous problem in an infinite strip with inhomogeneous boundary data, the reflection of errors in the boundaries, and two different boundary approximations for the leapfrog scheme with a fourth order accurate difference operator in space.

Skollermo, G.↗

Equatorial waves in the NCAR stratospheric general circulation model

Equatorially trapped wave modes are very important in the tropical stratospheric momentum balance. Kelvin waves and mixed Rossby-gravity waves are believed to be responsible for the quasi-biennial oscillation of the zonal winds in the equatorial lower stratosphere. Both Kelvin and mixed Rossby-gravity waves have been identified in observations and in numerical models. Kelvin and mixed Rossby-gravity waves are identified in a general circulation model extending from the surface into the mesosphere and looks at the effect on the waves of lowering the top of the model.

Boville, B. A.↗

Implementation of equivalent domain integral method in the two-dimensional analysis of mixed mode problems

An equivalent domain integral (EDI) method for calculating J-intergrals for two-dimensional cracked elastic bodies is presented. The details of the method and its implementation are presented for isoparametric elements. The total and product integrals consist of the sum of an area of domain integral and line integrals on the crack faces. The line integrals vanish only when the crack faces are traction free and the loading is either pure mode 1 or pure mode 2 or a combination of both with only the square-root singular term in the stress field. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all problems analyzed. The EDI method when applied to a problem of an interface crack in two different materials showed that the mode 1 and mode 2 components are domain dependent while the total integral is not. This behavior is caused by the presence of the oscillatory part of the singularity in bimaterial crack problems. The EDI method, thus, shows behavior similar to the virtual crack closure method for bimaterial problems.

Raju, I. S.↗

An equivalent domain integral for analysis of two-dimensional mixed mode problems

An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies subjected to mixed mode loading is presented. The total and product integrals consist of the sum of an area or domain integral and line integrals on the crack faces. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all the problems analyzed.

Raju, I. S.↗

Vapor condensation on liquid surface due to laminar jet-induced mixing: The effects of system parameters

The effects of system parameters on the interface condensation rate in a laminar jet induced mixing tank are numerically studied. The physical system consists of a partially filled cylindrical tank with a slightly subcooled jet discharged from the center of the tank bottom toward the liquid-vapor interface which is at a saturation temperature corresponding to the constant tank pressure. Liquid is also withdrawn from the outer part of the tank bottom to maintain the constant liquid level. The jet velocity is selected to be low enough such that the free surface is approximately flat. The effect of vapor superheat is assumed to be negligible. Therefore, the interface condensation rate can be determined from the resulting temperature field in the liquid region alone. The nondimensional form of the steady state conservation equations are solved by a finite difference method for various system parameters including liquid height to tank diameter ratio, tank to jet diameter ratio, liquid inflow to outflow area ratio, and a heat leak parameter which characterizes the uniform wall heat flux. Detailed analyses based on the numerical solutions are performed and simplified equations are suggested for the prediction of condensation rate.

Lin, Chin-Shun↗

Vapor condensation on liquid surface due to laminar jet-induced mixing - The effects of system parameters

The effects of system parameters on the interface condensation rate in a laminar jet induced mixing tank are numerically studied. The physical system consists of a partially filled cylindrical tank with a slightly subcooled jet discharged from the center of the tank bottom toward the liquid-vapor interface which is at a saturation temperature corresponding to the constant tank pressure. Liquid is also withdrawn from the outer part of the tank bottom to maintain the constant liquid level. The jet velocity is selected to be low enough such that the free surface is approximately flat. The effect of vapor superheat is assumed to be negligible. Therefore, the interface condensation rate can be determined from the resulting temperature field in the liquid region alone. The nondimensional form of the steady state conservation equations are solved by a finite difference method for various system parameters including liquid height to tank diameter ratio, tank to jet diameter ratio, liquid inflow to outflow area ratio, and a heat leak parameter which characterizes the uniform wall heat flux. Detailed analyses based on the numerical solutions are performed and simplified equations are suggested for the prediction of condensation rate.

Lin, Chin-Shun↗

An equivalent domain integral method in the two-dimensional analysis of mixed mode crack problems

An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies is presented. The details of the method and its implementation are presented for isoparametric elements. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented.

Raju, I. S.↗

Source conductance scaling for high frequency superconducting quasiparticle receivers

It has been suggested that the optimum source conductance G(sub s) for the superconductor-insulator-superconductor (SIS) quasiparticle mixer should have a l/f dependence. This would imply that the critical current density of SIS junctions used for mixing should increase as frequency squared, a stringent constraint on the design of submillimeter SIS mixers, rather than in simple proportion to frequency as previously believed. We have used Tucker's quantum theory of mixing for extensive numerical calculations to determine G(sub s) for an optimized SIS receiver. We find that G(sub s) is very roughly independent of frequency (except for the best junctions at low frequency), and discuss the implications of our results for the design of submillimeter SIS mixers.

Ke, Qing↗

Development of homotopy algorithms for fixed-order mixed H2/H(infinity) controller synthesis

A major difficulty associated with H-infinity and mu-synthesis methods is the order of the resulting compensator. Whereas model and/or controller reduction techniques are sometimes applied, performance and robustness properties are not preserved. By directly constraining compensator order during the optimization process, these properties are better preserved, albeit at the expense of computational complexity. This paper presents a novel homotopy algorithm to synthesize fixed-order mixed H2/H-infinity compensators. Numerical results are presented for a four-disk flexible structure to evaluate the efficiency of the algorithm.

Whorton, M.↗

A mapping method for the gravitational few-body problem with dissipation

Recently a new class of numerical integration methods - 'mixed variable symplectic integrators' - has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of magnitude faster than conventional ordinary differential equations (ODE) integration methods. Here we present a simple modification of this method to include small non-gravitational forces. The new scheme provides a similar advantage of computational speed for a larger class of problems in Solar System dynamics.

Malhotra, Renu↗

Low gravity transfer line chilldown

A code has been developed that solves for the transfer line chilldown time and flow and heat transfer characteristics in one-g environment. The code solves the transient, one dimensional, space averaged mass, momentum, and energy conservation equations for liquid-vapor two-phase flow in tubes. The physical configuration solved is that appropriate for bottom coolant injection in a vertically supported heated tube. Four distinct regions are considered consecutively: fully liquid, inverted annular, dispersed, and fully vapor flow. The conservation equations for both the liquid and the vapor are solved in each region separately. Also, in each region the mass and energy transport between each phase as well as the energy and momentum transport between the tube wall and the fluid are accounted for. A finite wall thickness is also considered. The model described above was solved numerically through a mixed finite difference scheme with forward time marching. The inverted annular regime was resolved using a semi-implicit finite differencing while the dispersed regime was solved explicitly. Also, a staggered mesh was used in which the velocity was resolved at mesh boundaries while all other field variables were resolved at the mesh centroids. Different mesh sizes were used depending on the region of solution. A coarse mesh was used in the dispersed flow region while a much finer mesh was used in both the inverted annular flow region and the tube walls.

B. N. Antar↗

The Robustness of Acoustic Analogies

Acoustic analogies for the prediction of flow noise are exact rearrangements of the flow equations N(right arrow q) = 0 into a nominal sound source S(right arrow q) and sound propagation operator L such that L(right arrow q) = S(right arrow q). In practice, the sound source is typically modeled and the propagation operator inverted to make predictions. Since the rearrangement is exact, any sufficiently accurate model of the source will yield the correct sound, so other factors must determine the merits of any particular formulation. Using data from a two-dimensional mixing layer direct numerical simulation (DNS), we evaluate the robustness of two analogy formulations to different errors intentionally introduced into the source. The motivation is that since S can not be perfectly modeled, analogies that are less sensitive to errors in S are preferable. Our assessment is made within the framework of Goldstein's generalized acoustic analogy, in which different choices of a base flow used in constructing L give different sources S and thus different analogies. A uniform base flow yields a Lighthill-like analogy, which we evaluate against a formulation in which the base flow is the actual mean flow of the DNS. The more complex mean flow formulation is found to be significantly more robust to errors in the energetic turbulent fluctuations, but its advantage is less pronounced when errors are made in the smaller scales.

Freund, J. B.↗

Integrated lunar materials manufacturing process

A manufacturing plant and process for production of oxygen on the moon uses lunar minerals as feed and a minimum of earth-imported, process materials. Lunar feed stocks are hydrogen-reducible minerals, ilmenite and lunar agglutinates occurring in numerous, explored locations mixed with other minerals in the pulverized surface layer of lunar soil known as regolith. Ilmenite (FeTiO.sub.3) and agglutinates contain ferrous (Fe.sup.+2) iron reducible by hydrogen to yield H.sub.2 O and metallic Fe at about 700.degree.-1,200.degree. C. The H.sub.2 O is electrolyzed in gas phase to yield H.sub.2 for recycle and O.sub.2 for storage and use. Hydrogen losses to lunar vacuum are minimized, with no net hydrogen (or any other earth-derived reagent) consumption except for small leaks. Feed minerals are surface-mined by front shovels and transported in trucks to the processing area. The machines are manned or robotic. Ilmenite and agglutinates occur mixed with silicate minerals which are not hydrogen-reducible at 700.degree.-1,200.degree. C. and consequently are separated and concentrated before feeding to the oxygen generation process. Solids rejected from the separation step and reduced solids from the oxygen process are returned to the mine area. The plant is powered by nuclear or solar power generators. Vapor-phase water electrolysis, a staged, countercurrent, fluidized bed reduction reactor and a radio-frequency-driven ceramic gas heater are used to improve thermal efficiency.

Gibson, Michael A.↗