Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “low viscous fluids”

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 91 records · Page 5

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .

Close-clearance impellers↗

Effects of viscosity and constraints on the dispersion and dissipation of waves in large blood vessels. I.

Investigation of the effects of blood viscosity on dissipation as well as dispersion of small waves in arteries and veins by means of a parametric study. A linearized analysis of axisymmetric waves in a cylindrical membrane that contains a viscous fluid indicates that there are two families of waves: a family of slow waves and one of fast waves. The faster waves are shown to be more sensitive to variations in the elastic properties of the medium surrounding the blood vessels and at high values of the frequency parameter alpha. At low values of alpha the effects of viscosity on attenuation are reversed.

Jones, E.↗

Application of computer generated color graphic techniques to the processing and display of three dimensional fluid dynamic data

Color coding techniques used in the processing of remote sensing imagery were adapted and applied to the fluid dynamics problems associated with turbofan mixer nozzles. The computer generated color graphics were found to be useful in reconstructing the measured flow field from low resolution experimental data to give more physical meaning to this information and in scanning and interpreting the large volume of computer generated data from the three dimensional viscous computer code used in the analysis.

Anderson, B. H.↗

Mantle flow pressure and the angle of subduction - Non-Newtonian corner flows

Corner flows of Newtonian and non-Newtonian fluids are used to model the flow in a subduction zone which is viscously driven by the motions of the converging plates and the descending slab. The pressures induced by the flow tend to lift the slab up beneath the overriding plate thereby offsetting the tendency of gravity to align the slab with the vertical. The low angles of subduction observed in Peru and Central Chile may be the result of strong dynamic pressures forcing the slab up against the overriding plate. Viscous coupling between the overriding plate and the downgoing slab is essential if the nonvertical dips of slabs are a consequence of the balance between gravitational and pressure torques. For a Newtonian mantle, shear stresses and pressures on the top of the slab are comparable. If the mantle is non-Newtonian, however, the pressures greatly exceed the shear stresses, for most acute dip angles. Thus frictional forces on the top and bottom surfaces of slabs are less important in resisting slab descent into a non-Newtonian mantle than they are in resisting penetration into a Newtonian mantle.

Tovish, A.↗

CFD Applications in Propulsion

An overview of various applications of Computational Fluid Dynamics (CFD) algorithms to propulsion problems is given. Topics of interest include incompressible, low speed compressible, transonic, and supersonic problems. A common family of algorithms is used for all applications and emphasis is placed on maintaining accuracy and convergence efficiency for all problems. Specific problems include pump hydrodynamics, simultaneous combustion and mixing in rocket engines, viscous nozzle flow, and CFD applications to combustion stability.

Merkle, Charles L.↗

Ice-atmosphere interactions in the central Arctic: Remotely-sensed and simulated ice concentration and motion

The response of the Beaufort Sea ice pack to the passage of strong low pressure systems is studied using SAR (synthetic aperture radar), AVHRR (advanced very high resolution radiometer), and SSM/I (special sensor microwave/imager) data combined with model simulations. The SAR, AVHRR, and modeled concentrations concur generally in showing a 1 to 5 percent decrease in ice fraction during the passage of the lows. SSM/I derived concentrations appear to underestimate overall concentration and overestimate the opening of the pack relative to the other data sets. Ice motion derived from SAR and AVHRR and model simulations using two ice rheologies agree reasonably well in direction and magnitude, although the cavitating fluid approximation tends to overestimate ice motion compared to a viscous plastic rheology. Ice motions from SAR, AVHRR, and buoy observations can be merged to yield uniform grids particularly well suited to synoptic scale studies and model validation.

Maslanik, J. A.↗

Convective and morphological instability in vapor crystal growth

Theoretical and experimental work on the fluid dynamics of physical vapor transport is reported. It is shown that diffusion in viscous interaction with container walls leads to concentration gradients normal to the main transport direction. Consequently any convection threshold is removed. This coupling with convective instabilities makes the theoretical interpretation of earthbound (anisotropic) morphological stability studies on vapor-solid interfaces intractable. Hence low-gravity experiments are suggested that will allow for the establishment of morphological stability criteria under diffusion controlled conditions.

Rosenberger, F.↗

An explicit Runge-Kutta method for turbulent reacting flow calculations

The paper presents a numerical method for the solution of the conservation equations governing steady, reacting, turbulent viscous flow in two-dimensional geometries, in both Cartesian and axisymmetric coordinates. These equations are written in Favre-averaged form and closed with a first order model. A two-equation K-epsilon model, where low Reynolds number and compressibility effects are included, and a modified eddy-break up model are used to simulate fluid mechanics turbulence, chemistry and turbulence-combustion interaction. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multi-stage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. The method is applied to the computation of both diffusion and premixed turbulent reacting flows. The computed temperature distributions compare favorably with experimental data.

Boretti, A. A.↗

Progress Toward Generation of a Navier-Stokes Database for a Harrier in Ground Effect

The Harrier YAV-8B aircraft is capable of vertical and short-field take-off and landing (V/STOL) by directing its four exhaust nozzles toward the ground, or conventional flight by rotating its nozzles into a horizontal position. The British Royal Air Force and the United States Marine Corps have used this aircraft for more than 30 years to provide a quick reaction time for troop support, and reduce the need for long runways. The success of this powered-lift (PL) vehicle has also prompted the more recent design of the Joint Strike Fighter (JSF). However there are significant safety issues that must be addressed when operating a PL vehicle in close proximity to the ground. Hot Gas Ingestion (HGI) by the inlets can result in a rapid loss of powered lift; and high-speed jet flows along the ground plane can induce low pressures underneath the vehicle, causing a 'suck-down' effect. Under these conditions, departure from controlled flight may occur. Moreover, unsteady ground vortices and jet fountains can affect the aircraft,s controllability and its proximity to ground troops. The viscous, time-dependent flow fields of PL vehicles are difficult to accurately and efficiently predict using Computational Fluid Dynamics (CFD). A number of researchers have used the time-dependent Reynolds-averaged Navier-Stokes (RANS) equations to compute flows for single and multiple jets in a cross-flow. A few have added some geometric complexity to the problem by computing flows for jet-augmented delta wings near a ground plane. Smith et.al. computed for the first time a single RANS solution about a simplified Harrier. This geometry included a fuselage, wing, leading edge root extension (LERX), inlets, and exhaust nozzles. All of these investigations cite two practical problems with computing these flows: 1) the need for improved solution accuracy; and, 2) the need for faster solution methods. We view the need for faster solution methods as key to improving the solution accuracy and making this class of computation more routine. One can hardly refine grids, explore the use of advanced turbulence models, and generate databases when it takes weeks of dedicated computer time for a single solution. Chaderjian, Ahmad, Pandya, and Murman have focused on reducing the time-to-solution for this very difficult and complex problem through process automation and exploitation of parallel computing. They began with the Harrier geometry reported, and added a deflected wing flap and empennage for greater realism. To date more than 80 solutions have been carried out. This paper will describe this process and progress made in reducing the time required to generate a simple longitudinal force and moment database for a Harrier in ground effect. It shows a typical snap-shot from an unsteady streakline animation, where fluid particles are colored by temperature. The ground vortex and a jet-fountain vortex are highlighted. It also shows a similar streakline image, where HGI occurs due to the vehicle in close proximity to the ground. It is show the mean lift coefficient as a function of angle of attack and height. The angle of attack range was 4 deg less than or = alpha less than or = 10 deg with an increment of 1 degree, and the height range was 10 ft less than or = h less than or = 30ft with an increment of 5 feet. This 35 solution database was extended to over 2500 cases using a monotone cubic-spline interpolation procedure. The suck-down effect (reduction of lift near the ground) is highlighted in the figure. The "cushion effect," the conventional reduction of lift as the vehicle moves out of ground effect, is also indicated. All 35 RANS solutions were obtained using 952 Silicon Graphics Origin 2000 and 3000 processors in dedicated mode for one week. Typically, 112 processors were assigned to each case. Some other cases used fewer processors to utilize all available CPUS. The final paper will report on the automation of the solution process, including: grid generation, job monitoring, solution completion criteria, and post processing. Moreover, improvements in parallel efficiency for a dual time-step algorithm for the RANS equations will also be presented. Results will be discussed in detail using unsteady streakline flow visualization to correlate unsteady flow structures with dominant aerodynamic frequencies. The stability derivatives, CL, and CL, will also be presented.

Chaderjian, Neal M.↗

Wall-Modeled Large-Eddy Simulation of Supersonic Parachute Inflation

Supersonic parachutes have been used in nearly every robotic mission to another planetary body with an atmosphere because they are one of the most mass-efficient ways to decelerate a payload to land on the surface. Short of performing flight tests in the upper Earth atmosphere, we currently cannot reliably predict a novel parachute system’s performance or potential failure modes, or even confidently explain it after the fact, as was the case with the Low Density Supersonic Decelerator flight tests and subsequent investigations. Fluid-structure interaction (FSI) simulations have the potential to bridge this gap. To this end, we present several improvements to the state-of-the-art for simulating supersonic parachutes using FSI: a higher effective resolution convective flux, an immersed boundary turbulent wall layer modeling approach to capture viscous effects, and a novel method to obtain a more realistic initial parachute shape. A recent supersonic parachute flight test is simulated and compared to measurements for the purposes of model validation.

ESM↗

The Iodine Satellite (iSAT) Propellant Feed System - Design and Development

The development, modeling, and testing of components and subsystems required to feed iodine propellant to a 200-W Hall thruster and cathode are described. This work aims to address design deficiencies and issues associated with the propellant feed system that were revealed by an integrated thruster-cathode-feed system test. The feed system design is modified to use materials that are more resistant to the highly-reactive nature of iodine propellant. Dynamic modeling indicates that the inclusion of additional constraints on feed system tubing will reduce the vibrationally-induced stresses that occur during launch. Full spacecraft thermal modeling show that the feed system heater power levels are sufficient to heat the tank and propellant lines to operating temperatures, where iodine in the tank is sublimed to supply propellant for operation and the tubing is elevated in temperature to keep propellant from redepositing to block the flow. Experiments are conducted to demonstrate that is it possible through the application of heating to clear an iodine deposit blocking the flow. Deposits in the low-pressure portion of the system near the exit to vacuum are shown to be relatively easy to remove in this manner while blockages forming upstream nearer to the higher-pressure propellant tank require significantly more effort to remove. Fluid flow modeling of the feed system is performed, exhibiting some qualitative agreement with experimental data. However, the highly viscous nature of the fluid flow and the dependence of the component flow coefficients on the Reynolds number are likely causes of the generally-poor quantitative agreement between the modeling results and experimentally-measured fluid flow properties.

Polzin, Kurt A.↗

Dissipative heating effects and end corrections for viscous Newtonian flow in high shear capillary tube viscometry

The effect of dissipation heating on the apparent viscosity measured in capillary tube viscometry is described in this paper. Conditions of low Reynolds number and high shear are assumed. End corrections to the tube flow are incorporated. The flow curves show decreasing apparent viscosity when the shear stress increases. The configuration of the flow curves plotted in logarithmic presentation are found to be identical for fluids with Newtonian behavior. Convection is the predominant mechanism in removal of the heat in short capillary tube. The estimated upper bound for the shear stress obtainable in short length capillary tubes appears to be of the order of magnitude of 10 MPa limited primarily by the pressure drop associated with the constant end correction from the flat ended inlet and exit of the tube.

Jakobsen, J.↗

Finite element analysis of low speed viscous and inviscid aerodynamic flows

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

Baker, A. J.↗

Computational Challenges of Viscous Incompressible Flows

Over the past thirty years, numerical methods and simulation tools for incompressible flows have been advanced as a subset of the computational fluid dynamics (CFD) discipline. Although incompressible flows are encountered in many areas of engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to the rather stringent requirements for predicting aerodynamic performance characteristics of flight vehicles, while flow devices involving low-speed or incompressible flow could be reasonably well designed without resorting to accurate numerical simulations. As flow devices are required to be more sophisticated and highly efficient CFD took become increasingly important in fluid engineering for incompressible and low-speed flow. This paper reviews some of the successes made possible by advances in computational technologies during the same period, and discusses some of the current challenges faced in computing incompressible flows.

Kwak, Dochan↗

Unusual Contact-Line Dynamics of Thick Films and Drops

We report several novel phenomena In contact-line and fingering dynamics of macroscopic spinning drops and gravity-driven films with dimensions larger than the capillary length. It is shown through experimental and theoretical analysis that such macroscopic films can exhibit various interfacial shapes, including multi valued ones, near the contact line due to a balance between the external body forces with capillarity. This rich variety of front shapes couples with the usual capillary, viscous, and intermolecular forces at the contact line to produce a rich and unexpected spectrum of contact-line dynamics. A single finger develops when part of the front becomes multivalued on a partially wetting macroscopic spinning drop in contrast to a different mechanism for microscopic drops of completely wetting fluids. Contrary to general expectation, we observe that, at high viscosity and low frequencies of rotation, the speed of a glycerine finger increases with increasing viscosity. Completely wetting Dow Corning 200 Fluid spreads faster over a dry inclined plane than a prewetted one. The presence of a thin prewetted film suppresses fingering both for gravity-driven flow and for spin coating. We analyze some of these unique phenomena in detail and offer qualitative physical explanations for the others.

Veretennikov, Igor↗

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗

Analysis of Unsteady Flow Phenomena: Shock - Vortex and Shock - Boundary Layer Interactions

The interaction of shock waves with vortices has received much attention in the past, mainly because shock-vortex interaction closely models the interaction of a shock wave with the coherent structures of a turbulent flow-field, and is a key feature in the broad-band shock noise for supersonic jets in off-project conditions. Chu and Kovasznay have shown that a weak disturbance in a viscous heat conducting fluid can be decomposed as the sum of three basic modes, namely acoustic, vortical and entropy mode; the interaction of any of these modes with a shock wave gives rise to all three disturbance modes downstream of the shock. The vortical mode is important since it constitutes the basis of the coherent structures that have been observed to dominate turbulence for low- to moderate-flow speed. Hollingsworth et al. have experimentally investigated the interaction of a cylindrical shock-induced starting vortex with a plane normal shock, and have shown that the interaction generates a cylindrical acoustic pulse that exhibits a quadrupolar structure consisting of four alternate compression and expansion regions centered around the transmitted vortex. The investigations of Hollingsworth and Richards have been extended by Dosanjh and Weeks that have analyzed the interaction of a columnar spiral vortex with a normal shock wave, thus obtaining quantitative measurements and confirming the generation of a progressive cylindrical wavefront of alternate compression-expansion nature. Naumann and Hermanns' have experimentally addressed the non-linear aspects of shock-vortex interaction, and have shown that the interaction causes both a diffraction and a reflection of the shock with a pattern consisting of either a regular-or a Mach-reflection depending on the shock and the vortex strengths. An attempt to theoretically explain the production of sound from the shock-vortex interaction was carried out by Ribner. Pao and Salas have numerically studied two-dimensional shock-vortex interactions.

Grasso, Francesco↗

Computational fluid dynamics studies of nuclear rocket performance

A CFD analysis of a low pressure nuclear rocket concept is presented with the use of an advanced chemical kinetics, Navier-Stokes code. The computations describe the flow field in detail,including gas dynamic, thermodynamic and chemical properties, as well as global performance quantities such as specific impulse. Computational studies of several rocket nozzle shapes are conducted in an attempt to maximize hydrogen recombination. These Navier-Stokes calculations, which include real gas and viscous effects, predict lower performance values than have been reported heretofore.

Stubbs, Robert M.↗