Search NASA⌕ Search

SEARCH · Search NASA

Results for “recirculation”

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 793 records · Page 44

Protein Crystal Growth Under Forced Solution Flow: Experimental Setup and General Response of Lysozyme

We have experimentally studied the effects of solution flow on the growth kinetics of the protein lysozyme. To this end, we have expanded our interferometry setup by a novel crystallization cell and solution recirculation system. This combination permits monitoring of interface morphology and kinetics with a depth resolution of 200 A at bulk flow rates of up to 2000 micron/s. Particular attention was paid to the prevention of protein denaturation that is often associated with the pumping of protein solutions. We found that at bulk flow rates it less than 250 microns/s the average growth rate and step velocity, R(sub avg) and upsilon(sub avg) increase with increasing it. This can be quantitatively understood in terms of the enhanced, convective solute supply to the interface. With high-purity solutions, it u greater than 250 microns/s lead to growth deceleration, and, at low supersaturations sigma, to growth cessation. When solutions containing approx. 1% of other protein impurities were used, growth deceleration occurred at any u greater than 0 and cessation in the low sigma experiments was reached at about half the it causing cessation with pure solution. The flow-induced changes in R(sub avg) and upsilon(sub avg) including growth cessation, were reversible and reproducible, independent of the direction of the u-changes and solution purity. Hence, we attribute the deceleration to the convection-enhanced supply of impurities to the interface, which at higher flow rates overpowers the effects of enhanced interfacial solute concentration. Most importantly, we found that convective transport leads to a significant reduction in kinetics fluctuations, in agreement with our earlier expectations for the lysozyme system. This supports our hypothesis that these long-term fluctuations represent an intrinsic response feature of the coupled bulk transport-interfacial kinetics system in the mixed growth control regime.

Vekilov, P. G.↗

Atmospheric Transport and Photochemistry of Ozone Over Central Southern Africa During the Southern Africa Fire-Atmosphere Research Initiative

Vertically integrated back and forward trajectories for the 300-200, 700-500 and surface-800 hPa levels are calculated using Pretoria as point of origin for the Southern Africa Fire-Atmosphere Research Initiative (SAFARI) period September-October 1992. The transport fields are then combined to show both horizontal and vertical transport of air to and from Pretoria at the different levels. Air transport patterns in the vertical are linked to the occurrence of absolutely stable layers which are also evident in the 16 ozonesonde profiles recorded at Pretoria during SAFARI. The coherence of the stratification based on dynamical and ozone analysis permits the use of mean ozone profiles with air volume fluxes to interpret the ozone in terms of photochemistry and transport within stable layers. Extensive recirculation across the meridional plane at Pretoria implies that advection of ozone is slow and that photochemistry is responsible for the observed vertical structure over central southern Africa in September and October 1992. Requisite ozone formation rates are supported by model analysis of ozone and ozone precursors measured from SAFARI and Transport and Atmospheric Research Chemistry near the Equator-Atlantic aircraft.

Tyson, P. D.↗

Magnetic Damping of g-Jitter Driven Flows: 3-D Calculations

A 3-D numerical model is developed to represent the oscillating natural convection induced in a cylindrical cavity filled with Ga-doped germanium with and without the presence of an external magnetic field. The model is developed based on the penalty-finite element solution of the equations describing the transport of momentum, heat and solutal element as well as the electromagnetic field distribution in the melt pool. Automatic time step control is applied to help speed up the calculations. Numerical simulations are conducted to study the convection and magnetic damping effects as a function of frequency, directions and amplitudes of g-jitter and also the direction and magnitudes of the applied magnetic fields. The results show that the g-jitter driven flow is time dependent and exhibits a complex recirculating convection pattern in three dimensions and that an applied magnetic field can be employed to suppress this deleterious convective flow and both magnitude and orientation of the applied field are important in magnetic damping of the g-jitter induced convective flows.

Shang, D. Y.↗

Laminar Dust Flames: A Program of Microgravity and Ground Based Studies at McGill

Fundamental knowledge of heterogeneous combustion mechanisms is required to improve utilization of solid fuels (e.g. coal), safe handling of combustible dusts in industry, and solid propulsion systems. The objective of the McGill University research program on dust combustion is to obtain a reliable set of data on basic combustion parameters for dust suspensions (i.e. laminar burning velocity, flame structure, quenching distance, flammability limits, etc.) over a range of particle sizes, dust concentrations, and types of fuel. This set of data then permits theoretical models to be validated and, when necessary, new models to be developed to describe the detailed reaction mechanisms and transport processes. Microgravity is essential to the generation of a uniform dust suspension of arbitrary particle size and concentration. When particles with a characteristic size on the order of tens of microns are suspended, they rapidly settle in a gravitational field. To maintain a particulate in suspension for time duration adequate to carry out combustion experiments invariably requires continuous convective flow in excess of the gravitational settling velocity (which is comparable with and can even exceed the dust laminar burning velocity). This makes the experiments turbulent in nature and thus renders it impossible to study laminar dust flames. Even for small particle sizes on the order of microns, a stable laminar dust flow can be maintained only for relatively low dust concentrations at normal gravity conditions. High dust loading leads to gravitational instability of the dust cloud and to the formation of recirculation cells in the dust suspension in a confined volume, or to the rapid sedimentation of the dense dust cloud, as a whole, in an unconfined volume. Many important solid fuels such as carbon and boron also have low laminar flame speeds (of the order of several centimeters per second). Convection that occurs in combustion products due to buoyancy disrupts the low speed dust flames and makes observation of such flames at normal gravity difficult.

Goroshin, Sam↗

Solar disturbances and correlated geospace responses: Relativistic magnetospheric electron acceleration

The role of high-speed solar wind streams in driving relativistic electron acceleration within the earth's magnetosphere is discussed based on International Solar-Terrestrial Physics (ISTP) Observatory and related spacecraft observations. A 'recirculation' mechanism for electron acceleration and redistribution was invoked. Recently, an increase in the number of coronal mass ejections (CMEs) and related 'magnetic clouds' was seen at 1 AU. As these CME/cloud systems interact with the earth's magnetosphere, they are able to produce rapid enhancements in the magnetospheric electron population. The relativistic electron signatures observed by the POLAR, SAMPEX, and other spacecraft during recent magnetic cloud events, especially January 1997 and May 1997, were compared and contrasted. In these cases, there were large solar wind and IMF changes during the cloud passages and very rapid energetic electron acceleration was observed. The relative geoeffectiveness of these events is examined and 'space weather' predicatability is assessed.

Baker, D. N.↗

A Generalized Wall Function

The asymptotic solutions, described by Tennekes and Lumley (1972), for surface flows in a channel, pipe or boundary layer at large Reynolds numbers are revisited. These solutions can be extended to more complex flows such as the flows with various pressure gradients, zero wall stress and rough surfaces, etc. In computational fluid dynamics (CFD), these solutions can be used as the boundary conditions to bridge the near-wall region of turbulent flows so that there is no need to have the fine grids near the wall unless the near-wall flow structures are required to resolve. These solutions are referred to as the wall functions. Furthermore, a generalized and unified law of the wall which is valid for whole surface layer (including viscous sublayer, buffer layer and inertial sublayer) is analytically constructed. The generalized law of the wall shows that the effect of both adverse and favorable pressure gradients on the surface flow is very significant. Such as unified wall function will be useful not only in deriving analytic expressions for surface flow properties but also bringing a great convenience for CFD methods to place accurate boundary conditions at any location away from the wall. The extended wall functions introduced in this paper can be used for complex flows with acceleration, deceleration, separation, recirculation and rough surfaces.

Shih, Tsan-Hsing↗

Propellant Densification for Launch Vehicles: Simulation and Testing 1999

One of the many key technologies required to make single-stage to orbit an actuality, the technology of sub-cooling cryogenic propellants below their normal saturation temperatures and thereby making them more dense, is unquestionably on its way towards full and practical realization. The technology of Propellant Densification has been the subject of an extensive research and development program overseen by Lockheed Martin Michoud Space Systems and NASA Glenn Research Center over the past several years. This paper presents a status report of this research and development. Specifically examined within this paper is the status of the current and continuing efforts on the mathematical simulation of the in-tank propellant densification process currently baselined for the Lockheed Martin VentureStar Reusable Launch Vehicle (RLV). Keys to this modeling effort are an understanding and quantification of the effects of thermal stratification and the ability to capture the complex and unique multiple section tank geometries being proposed for future launch vehicles. A simulation that properly captures these phenomena has been developed by Lockheed Martin. Also discussed is the significant test program that has been undertaken in coordination with NASA Glenn Research Center. In this testing, the liquid hydrogen recirculation and densification process was simulated and the thermal stratification of the densified propellant was recorded throughout the tank. This testing marks the first time that such a process has been carried out within a multiple-lobe, flight-similar tank. The results from this testing have gone a long way towards grounding the mathematical models and towards demonstrating the readiness of the technology for near-term use. A further and even more ambitious test program examining the production and utilization of densified propellants is being planned for late-autumn 1999. An overview of these plans is presented.

Knowles, Timothy E.↗

Control of Meridional Flow by a Non-Uniform Rotational Magnetic Field

The diffusive mass transfer of species during crystal growth in vertical ampoules is significantly affected by fluid flow in the liquid mother phase (melt). For electrically conductive melts, an elegant way of remotely inducing and controlling this flow is by utilizing a uniform rotational magnetic field (RMF) in the transverse direction. It induces an azimuthal flow which tends to homogenize the thermal and solutal fields. The rotating field also reduces the diffusion boundary layer, stabilizes temperature fluctuations, and promotes better overall crystal growth. For moderate strengths of the applied magnetic field (2-20 m Tesla) with frequencies of up to 400 Hz, the induced secondary meridional flow becomes significant. It typically consists of one roll at the bottom of the liquid column and a second roll (vortex) at the top. The flow along the centerline (ampoule axis) is directed from the growing solid (interface) towards the liquid (melt). In case of convex interfaces (e.g. in floating zone crystal growth) such flow behavior is beneficial since it suppresses diffusion at the center. However, for concave interfaces (e.g. vertical Bridgman crystal growth) such a flow tends to exacerbate the situation in making the interface shape more concave. It would be beneficial to have some control of this meridional flow- for example, a single recirculating cell with controllable direction and flow magnitude will make this technique even more attractive for crystal growth. Such flow control is a possibility if a non-uniform PNE field is utilized for this purpose. Although this idea has been proposed earlier, it has not been conclusively demonstrated so far. In this work, we derive the governing equations for the fluid dynamics for such a system and obtain solutions for a few important cases. Results from parallel experimental measurements of fluid flow in a mercury column subjected to non-uniform RMF will also be presented.

Mazuruk, Konstantin↗

Unshrouded Centrifugal Turbopump Impeller

The ratio of rocket engine thrust to weight is a limiting constraint in placing more payload into orbit at a low cost. A key component of an engine's overall weight is the turbopump weight, Reducing the turbopump weight can result in significant engine weight reduction and hence, increased delivered payload. There are two main types of pumps: centrifugal and axial pumps. These types of pumps can be further sub-divided into those with shrouds and those without shrouds (unshrouded pumps). Centrifugal pumps can achieve the same pump discharge pressure as an axial pump and it requires fewer pump stages and lower pump weight than an axial pump. Also, with unshrouded centrifugal pumps (impeller), the number of stages and weight can be further reduced. However. there are several issues with regard to using an unshrouded impeller: 1) there is a pump performance penalty due to the front open face recirculation flow; 2) there is a potential pump axial thrust problem from the unbalanced front open face and the back shroud face; and, 3) since test data is very linu'ted for this configuration, there is uncertainty in the magnitude and phase of the rotordynamics forces due to the front impeller passage. The purpose of the paper is to discuss the design of an unshrouded impeller and to examine the design's hydrodynamic performance, axial thrust, and rotordynamics performance. The design methodology will also be discussed. This work will help provide some guidelines for unshrouded impeller design. In particular, the paper will discuss the design of three unshrouded impellers - one with 5 full and 5 partial blades (5+5). one with 6+6 blades and one with 8+8 blades. One of these designs will be selected for actual fabrication and flow test. Computational fluid dynamics (CFD) is used to help design and optimize the unshrouded impeller. The relative pump performance penalty is assessed by comparing the CFD results of the unshrouded impeller with the equivalent shrouded impeller for a particular design. Limited unshrouded - versus - shrouded impeller data from the J-2 pump is used to anchor the CFD. Since no detailed impeller blade force data is available, axial thrust and rotordynamic force predictions are based on the CFD model. For the axial thrust, the impeller front flow passage axial force is integrated from the CFD results and compared to the equivalent shrouded impeller axial force. For the rotordynamics forces, the fluid reaction forces are computed from unsteady flow CFD results using a moving boundary method; the rotor- shaft is moved at several whirl-to-speed frequency ratios to extract the rotordynamics coefficients.

Prueger, George↗

Effects of Radiative Emission and Absorption on the Propagation and Extinction of Premixed Gas Flames

Premixed gas flames in mixtures of CH4, O2, N2, and CO2 were studied numerically using detailed chemical and radiative emission-absorption models to establish the conditions for which radiatively induced extinction limits may exist independent of the system dimensions. It was found that reabsorption of emitted radiation led to substantially higher burning velocities and wider extinction limits than calculations using optically thin radiation models, particularly when CO2, a strong absorber, is present in the unburned gas, Two heat loss mechanisms that lead to flammability limits even with reabsorption were identified. One is that for dry hydrocarbon-air mixtures, because of the differences in the absorption spectra of H2O and CO2, most of the radiation from product H2O that is emitted in the upstream direction cannot be absorbed by the reactants. The second is that the emission spectrum Of CO2 is broader at flame temperatures than ambient temperature: thus, some radiation emitted near the flame front cannot be absorbed by the reactants even when they are seeded with CO2 Via both mechanisms, some net upstream heat loss due to radiation will always occur, leading to extinction of sufficiently weak mixtures. Downstream loss has practically no influence. Comparison with experiment demonstrates the importance of reabsorption in CO2 diluted mixtures. It is concluded that fundamental flammability limits can exist due to radiative heat loss, but these limits are strongly dependent on the emission-absorption spectra of the reactant and product -gases and their temperature dependence and cannot be predicted using gray-gas or optically thin model parameters. Applications to practical flames at high pressure, in large combustion chambers, and with exhaust-gas or flue-gas recirculation are discussed.

Ju, Yiguang↗

Finite Element Modeling of Magnetically-Damped Convection during Solidification

A fully 3-D, transient finite element model is developed to represent the magnetic damping effects on complex fluid flow, heat transfer and electromagnetic field distributions in a Sn- 35.5%Pb melt undergoing unidirectional solidification. The model is developed based on our in- house finite element code for the fluid flow, heat transfer and electromagnetic field calculations. The numerical model is tested against numerical and experimental results for water as reported in literature. Various numerical simulations are carried out for the melt convection and temperature distribution with and without the presence of a transverse magnetic field. Numerical results show that magnetic damping can be effectively applied to stabilize melt flow, reduce turbulence and flow levels in the melt and over a certain threshold value a higher magnetic field resulted in a greater reduction in velocity. Also, for the study of melt flow instability, a long enough running time is needed to ensure the final fluid flow recirculation pattern. Moreover, numerical results suggest that there seems to exist a threshold value of applied magnetic field, above which magnetic damping becomes possible and below which the 0 convection in the melt is actually enhanced.

deGroh, H. C.↗

A Seasonal Air Transport Climatology for Kenya

A climatology of air transport to and from Kenya has been developed using kinematic trajectory modeling. Significant months for trajectory analysis have been determined from a classification of synoptic circulation fields. Five-point back and forward trajectory clusters to and from Kenya reveal that the transport corridors to Kenya are clearly bounded and well defined. Air reaching the country originates mainly from the Saharan region and northwestern Indian Ocean of the Arabian Sea in the northern hemisphere and from the Madagascan region of the Indian Ocean in the southern hemisphere. Transport from each of these source regions show distinctive annual cycles related to the northeasterly Asian monsoon and the southeasterly trade wind maximum over Kenya in May. The Saharan transport in the lower troposphere is at a maximum when the subtropical high over northern Africa is strongly developed in the boreal winter. Air reaching Kenya between 700 and 500 hPa is mainly from Sahara and northwest India Ocean flows in the months of January and March, which gives way to southwest Indian Ocean flow in May and November. In contrast, air reaching Kenya at 400 hPa is mainly from southwest Indian Ocean in January and March, which is replaced by Saharan transport in May and November. Transport of air from Kenya is invariant, both spatially and temporally, in the tropical easterlies to the Congo Basin and Atlantic Ocean in comparison to the transport to the country. Recirculation of air has also been observed, but on a limited and often local scale and not to the extent reported in southern Africa.

Gatebe, C. K.↗

Processes Affecting Tropospheric Ozone Inferred from Ozonesonde and Other Tracer Data from the R/V R H Brown Atlantic Cruise (37N-34S) in January-February 1999

During the Aerosols-99 trans-Atlantic cruise from Norfolk, VA, to Cape Town, South Africa, 22 ozonesondes were launched from the NOAA R/V R H Brown between 17 Jan and 6 Feb 1999, with all sondes but one reaching 30 km. A composite of ozone profiles along the transect shows high free tropospheric ozone (up to 100 ppbv at 9 km) between 5N and 20S, a coherent feature straddling either side of the ITCZ. Latitudinal variations of tropospheric ozone are interpreted using correlative measurements of surface ozone, CO, water vapor, and aerosol optical thickness (column absorbance) measured from the ship. Elevated ozone in the lower troposphere results from photochemical reactions of precursors emitted by biomass burning north of the ITCZ. However, the greatest ozone mixing ratios are in the mid-troposphere south of the ITCZ, which gives evidence of interhemispheric transport. Column-integrated tropospheric ozone, 35 DU from 0-16 km, agrees with that derived from the TOMS satellite by the modified-residual method [Thompson and Hudson, 1999]. NCEP wind fields, ship-launched radiosondes and back trajectories are consistent with a picture of recirculating air parcels centered in the tropical Atlantic region which is identified with the maximum wave-one amplitude in total ozone seen in sondes and by satellite.

Thompson, Anne M.↗

Flow Simulation of Solid Rocket Motors: Injection Induced Water-Flow Tests from Porous Media - 1

Prior to selecting a proper porous material for use in simulating the internal port flow of a solid rocket motor (SRM), in cold-flow testing, the flow emerging from porous materials is experimentally investigated. The injection-flow emerging from a porous matrix always exhibits a lumpy velocity profile that is spatially stable and affects the development of the longitudinal port flow. This flow instability, termed pseudoturbulence, is an inherent signature of the porous matrix and is found to generally increase with the wall porosity and with the injection flow rate. Visualization studies further show that the flow from porous walls made from shaving-type material (sintered stainless-steel) exhibits strong recirculation zones that are conspicuously absent in walls made from nodular or spherical material (sintered bronze). Detailed flow visualization observations and hot-film measurements are reported from tests of injection-flow and a coupled cross-flow from different porous wall materials. Based on the experimental data, discussion is provided on the choice of suitable material for SRM model testing while addressing the consequences and shortcomings from such a test.

Ramachandran, N.↗

Incompressible Turbulent Wing-Body Junction Flow

The overall objective of this study is to contribute to the optimized design of fan bypass systems in advanced turbofan engines. Increasing the engine bypass ratios have provided a major boost in engine performance improvement over the last fifty years. An engine with high bypass ratio (11-16:1) such as the Advanced Ducted Propulsion (ADP) is being developed and is expected to provide an additional 25% improvement in overall efficiency over the early turbofans. Such significant improvements in overall efficiency would reduce the cost per seat mile, which is a major government and Industry challenge for the 21th century. The research is part of the Advanced Subsonic Technology (AST) program that involves a NASA, U.S. Industry and FAA partnership with the goal of a safe and highly productive global air transportation system. The immediate objective of the study is to perform numerical simulation of duct-strut interactions to elucidate the loss mechanisms associated with this configuration that is typical of advanced turbofan engines such as ADP. However, at present experimental data for a duct-strut configuration are not available. Thus, as a first step a wing-body junction flow would be studied and is the specific objective of the present study. At the outset it is to be recognized that while duct-strut interaction flow is similar to that of wing-body junction flows, there are some differences owing to the presence of a wall at both ends of the strut. Likewise, some differences are due to the sheared inflow (as opposed to a uniform inflow) velocity profile. It is however expected that some features of a wing-body junction flow would persist. Next, some of the salient aspects of the complex flow near a wing-body junction, as revealed by various studies reported in the literature will be reviewed. One of the principle characteristics of the juncture flow, is the presence of the mean flow components in a plane perpendicular to the direction of the oncoming free-stream flow. The lateral curvature of the wing/strat causes the oncoming turbulent layer to skew about am axis (x-axis) parallel to the plane (xz-plane) of the mean shear. This is the principle mechanism for the generation of secondary flow. Such skew-induced secondary flows are slow to be attenuated by Reynolds stresses. Additional contribution to the generation of secondary flow comes from anisotropies in Reynolds stresses. Upstream of the strut, the mean-vorticity is directed span wise (along the y-direction). The presence of secondary flow in the vicinity of the strut causes the vorticity to stretch around the obstacle in a horse-shoe shape, with each leg having a vorticity of the opposite sense. The blockage effect of the strut imposes a severe adverse pressure gradient on the oncoming turbulent shear layer, causing boundary layer separation ahead of the leading edge, resulting in a vortex that rolls up and flows downstream into the juncture region. The separation vortices trailing in the wake of the wing can alter the lift or drag characteristics of the surfaces downstream of the wing-body juncture. Likewise, on submarines, the wake flow behind the appendage can degrade the performance of the propeller located downstream. The complex nature of this flow is caused by the presence of all six components of Reynolds stresses. Devenport and Simpson report that in the vicinity of the horse-shoe vortex there is intense recirculation with turbulent stresses being much larger than those normally observed in turbulent flows. These features contribute to making this flow a challenge to predict numerically. Some of the past studies provide useful insights into this flow that would guide our numerical efforts. In measurements reported by Shabaka and Bradshaw, the eddy viscosity tensor is seen to be non-isotropic and has negative components in certain regions. In an effort to evaluate the closure assumptions of various turbulence models, Devenport and Simpson used their own extensive measurements in juncture flows around the nose of a wing-body junction. Measured values of mean-velocity and/or turbulence kinetic energy was used to predict the magnitude of the shear stress vector. Algebraic stress models performed the best followed by Cebeci-Smith eddy viscosity model. The flow is reported to be dominated by a pressure field produced by the wing and the velocity field generated by the horseshoe vortex that is wrapped around the junction between the wing and wall. Kubendran et al. conclude from an experimental study that the shape of leading edge of the wing as characterized by its slenderness ratio is a major factor in determining the flow fields in the juncture region. The more thinner the leading edge of the juncture, the weaker the horseshoe vortex is. Also, with a slender leading edge, the secondary flow in the juncture would be due mainly to the cross-stream gradients of Reynolds stresses rather than due to a lateral skewing of the shear layer.

Krishnamurthy, R.↗

Flight Control Using Distributed Shape-Change Effector Arrays

Recent discoveries in material science and fluidics have been used to create a variety of novel effector devices that offer great potential to enable new approaches to aerospace vehicle flight control. Examples include small inflatable blisters, shape-memory alloy diaphragms, and piezoelectric patches that may be used to produce distortions or bumps on the surface of an airfoil to generate control moments. Small jets have also been used to produce a virtual shape-change through fluidic means by creating a recirculation bubble on the surface of an airfoil. An advanced aerospace vehicle might use distributed arrays of hundreds of such devices to generate moments for stabilization and maneuver control, either augmenting or replacing conventional ailerons, flaps or rudders. This research demonstrates the design and use of shape-change device arrays for a tailless aircraft in a low-rate maneuvering application. A methodology for assessing the control authority of the device arrays is described, and a suite of arrays is used in a dynamic simulation to illustrate allocation and deployment methodologies. Although the authority of the preliminary shape-change array designs studied in this paper appeared quite low, the simulation results indicate that the effector suite possessed sufficient authority to stabilize and maneuver the vehicle in mild turbulence.

Raney, David L.↗

Generalized Wall Function for Complex Turbulent Flows

A generalized wall function was proposed by Shih et al., (1999). It accounts the effect of pressure gradients on the flow near the wall. Theory shows that the effect of pressure gradients on the flow in the inertial sublayer is very significant and the standard wall function should be replaced by a generalized wall function. Since the theory is also valid for boundary layer flows toward separation, the generalized wall function may be applied to complex turbulent flows with acceleration, deceleration, separation and recirculation. This paper is to verify the generalized wall function with numerical simulations for boundary layer flows with various adverse and favorable pressure gradients, including flows about to separate. Furthermore, a general procedure of implementation of the generalized wall function for National Combustion Code (NCC) is described, it can be applied to both structured and unstructured CFD codes.

Shih, Tsan-Hsing↗

A Design Tool for Liquid Rocket Engine Injectors

A practical design tool which emphasizes the analysis of flowfields near the injector face of liquid rocket engines has been developed and used to simulate preliminary configurations of NASA's Fastrac and vortex engines. This computational design tool is sufficiently detailed to predict the interactive effects of injector element impingement angles and points and the momenta of the individual orifice flows and the combusting flow which results. In order to simulate a significant number of individual orifices, a homogeneous computational fluid dynamics model was developed. To describe sub- and supercritical liquid and vapor flows, the model utilized thermal and caloric equations of state which were valid over a wide range of pressures and temperatures. The model was constructed such that the local quality of the flow was determined directly. Since both the Fastrac and vortex engines utilize RP-1/LOX propellants, a simplified hydrocarbon combustion model was devised in order to accomplish three-dimensional, multiphase flow simulations. Such a model does not identify drops or their distribution, but it does allow the recirculating flow along the injector face and into the acoustic cavity and the film coolant flow to be accurately predicted.

Farmer, R.↗