Viscous-transonic flow in the accretion and stellar-wind problems.
Closed form solution expressing transonic flow in idealized stellar wind and spherical accretion problems including viscosity and thermoconductivity effects
TRANSONIC FLOW↗
SEARCH · Engineering Papers
Results for “viscous effects”
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.
Experimental study of nozzle wall boundary layers at Mach numbers 20 to 47
The nozzle wall boundary layer of an M-50 helium tunnel was investigated with pitot pressure, total temperature, skin friction, and wall heat transfer measurements at five stations and hot wire measurements at two stations. The results indicated that the boundary layer was turbulent with a thick viscous sublayer. Pressure gradients were observed across the boundary layer; the effect of these gradients on the equations of motion are discussed. The direct skin friction measurements were higher than expected from empirical predictions; the Reynolds analogy factors however were lower than expected. Hot wire measurements indicated mass flow fluctuations as large as 80 percent of the local mean mass flow at the edge of the viscous sublayer with a maximum value relative to the edge mass flow of about 15 percent at gamma about equal to 0.8.
Kemp, J. H., Jr.↗
Impacts of Mixed-Wettability on Brine Drainage and Supercritical CO 2 Storage Efficiency in a 2.5-D Heterogeneous Micromodel
Geological carbon storage (GCS) involves unstable drainage processes, the formation of patterns in a morphologically unstable interface between two fluids in a porous medium during drainage. The unstable drainage processes affect CO 2 storage efficiency and plume distribution and can be greatly complicated by the mixed-wet nature of rock surfaces common in hydrocarbon reservoirs where supercritical CO 2 (scCO 2 ) is used in enhanced oil recovery. Here, we performed scCO 2 injection (brine drainage) experiments at 8.5 MPa and 45°C in heterogeneous micromodels, two mixed-wet with varying water- and intermediate-wet patches, and one water-wet. The flow regime changes from capillary fingering through crossover to viscous fingering in the micromodels of the same pore geometry but different wetting surfaces at displacement rates with logCa (capillary number) increasing from –8.1 to –4.4. While the mixed-wet micromodel with uniformly distributed intermediate-wet patches yields ~0.15 scCO 2 saturation increase at both capillary fingering and crossover flow regimes (–8.1 ≤ logCa ≤ –6.1), the one heterogeneous wetting to scCO 2 results in ~0.09 saturation increase only at the crossover flow regime (–7.1 ≤ logCa ≤ –6.1). The interconnected flow paths in the former are quantified and compared to the channelized scCO 2 flow through intermediate-wet patches in the latter by topological analysis. At logCa > –6.1 (near well), the effects of wettability and pore geometry are suppressed by strong viscous force. Both scCO 2 saturation and distribution suggest the importance of wettability on CO 2 storage efficiency and plume shape in reservoirs and capillary leakage through caprock at GCS conditions.
58 GEOSCIENCES↗
Beyond the Stokes approximation: shallow visco-elastic ice-sheet models
Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.
58 GEOSCIENCES↗
The effect of wall compliance on the Goertler vortex instability
The stability of the flow of a viscous incompressible fluid over a curved compliant wall to longitudinal Goertler vortices is investigated. The compliant wall is modeled by a particularly simple equation relating the induced wall displacement to the pressure in the overlying fluid. Attention is restricted to the large Goertler number regime; this regime being appropriate to the most unstable Goertler mode. The effect of wall compliance on this most unstable mode is investigated.
Denier, J. P.↗
Tidal Modulation of Ice-shelf Flow: a Viscous Model of the Ross Ice Shelf
Three stations near the calving front of the Ross Ice Shelf, Antarctica, recorded GPS data through a full spring-neap tidal cycle in November 2005. The data revealed a diurnal horizontal motion that varied both along and transverse to the long-term average velocity direction, similar to tidal signals observed in other ice shelves and ice streams. Based on its periodicity, it was hypothesized that the signal represents a flow response of the Ross Ice Shelf to the diurnal tides of the Ross Sea. To assess the influence of the tide on the ice-shelf motion, two hypotheses were developed. The first addressed the direct response of the ice shelf to tidal forcing, such as forces due to sea-surface slopes or forces due to sub-ice-shelf currents. The second involved the indirect response of ice-shelf flow to the tidal signals observed in the ice streams that source the ice shelf. A finite-element model, based on viscous creep flow, was developed to test these hypotheses, but succeeded only in falsifying both hypotheses, i.e. showing that direct tidal effects produce too small a response, and indirect tidal effects produce a response that is not smooth in time. This nullification suggests that a combination of viscous and elastic deformation is required to explain the observations.
Tidal↗
A multi-grid method for the computation of viscid/inviscid interactions on airfoils
An improved version of the 'GRUMFOIL' code has been developed for the computation of airfoil flows. The method employs a conservative difference scheme for the potential equation, integral methods for the boundary layer, and viscous coupling conditions that fully account for the wake and strong interaction effects at trailing edges. The improved version uses Jameson's 'MAD' scheme to accelerate convergence of the inviscid solution, an improved 2nd order artificial viscosity and far field 'TARE' correction to reduce spatial truncation errors and Carter's semi-inverse method for the viscous solution. Results are presented which demonstrate a factor of ten reduction in computing cost.
Melnik, R. E.↗
On higher eigenstates of Goertler instability in compressible flow
The effect of secondary instabilities of Goertler vortices with respect to their growth and amplification was investigated. Instability of viscous flow along a concave surface is in the form of counterrotating streamwise vortices commonly known as Goertler vortices. The effects of suction on higher modes of Goertler vortices are also considered.
Goglia, G. L.↗
A three-dimensional viscous flow analysis for the helicopter tip vortex generation problem
The tip vortex flow field occurring in the vicinity of the tip region of a a helicopter rotor blade is a very complicated three-dimensional, viscous flow phenomenon. The details of the flow in the tip region can have a major effect in determining the generated rotor noise and can significantly affect the performance and dynamic loading of the rotor blade. The three-dimensional viscous subsonic tip vortex generation processes is investigated by a numerical procedure which allows spatial forward-marching integration, utilizing flow approximations from the velocity-decomposition approach of Briley and McDonald. The approach has been applied to compute the laminar and turbulent tip vortex flows for a constant thickness slab airfoil with a square tip, a constant thickness slab airfoil with a half round tip and a NACA 0012 airfoil with a half round tip. The basic mechanism of the tip vortex generation process as well as the prediction of vortex appearance, strength and secondary flow shown by the calculations are in qualitative agreement with experimental results.
Lin, S. J.↗
Boundary layer effects on the flowfield about flow-through nacelles
A numerical study has been made of boundary layer effects on the flowfield about flow-through nacelles. A viscous-inviscid interacting computational model for investigating the problem was constructed by coupling a three-dimensional explicit Euler solution procedure with a compressible, 'lag-entrainment' integral boundary layer solution technique. Solutions obtained with the interacting model for long-duct turbofan engine nacelles at a free stream Mach number of 0.8 and an angle of attack of 0 deg are compared with experimental data. The numerical experiments have revealed insight into the boundary layer effects on the interactions between the internal and external flows.
Compton, W. B.↗
Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface
Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.
Palaparthi, Ravi↗
Long discontinuous fiber composite structure: Forming and structural mechanics
Cost effective composite structure has motivated the investigation of several new approaches to develop composite structure from innovative material forms. Among the promising new approaches is the conversion of planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. In the present study, the authors have established a framework which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. The initial study focuses upon the establishment of micromechanics models for prediction of the effective anisotropic viscosities of the oriented fiber assembly in a viscous matrix. Next, the developed constitutive relation is employed through an analogy with incompressible elasticity to exercise the finite element technique for determination of local fiber orientation and laminate thickness after forming. Results are presented for the stretch bending of a curved beam from an arbitrary composite laminate and the bulging of a clamped sheet. Structural analyses are conducted to determine the effect of microstructure on the performance of curved beams manufactured from long discontinuous fiber composites. For the purposes of this study, several curved beams with ideal and non-ideal microstructures are compared for response under pure bending. Material parameters are determined from a separate microstructural analysis.
Pipes, R. B.↗
Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields
The primary accomplishments of the project were as follows: (1) From an overall standpoint, the primary accomplishment of this research was the development of a complete gasdynamic-radiatively coupled nonequilibrium viscous shock layer solution method for axisymmetric blunt bodies. This method can be used for rapid engineering modeling of nonequilibrium re-entry flowfields over a wide range of conditions. (2) Another significant accomplishment was the development of an air radiation model that included local thermodynamic nonequilibrium (LTNE) phenomena. (3) As part of this research, three electron-electronic energy models were developed. The first was a quasi-equilibrium electron (QEE) model which determined an effective free electron temperature and assumed that the electronic states were in equilibrium with the free electrons. The second was a quasi-equilibrium electron-electronic (QEEE) model which computed an effective electron-electronic temperature. The third model was a full electron-electronic (FEE) differential equation model which included convective, collisional, viscous, conductive, vibrational coupling, and chemical effects on electron-electronic energy. (4) Since vibration-dissociation coupling phenomena as well as vibrational thermal nonequilibrium phenomena are important in the nonequilibrium zone behind a shock front, a vibrational energy and vibration-dissociation coupling model was developed and included in the flowfield model. This model was a modified coupled vibrational dissociation vibrational (MCVDV) model and also included electron-vibrational coupling. (5) Another accomplishment of the project was the usage of the developed models to investigate radiative heating. (6) A multi-component diffusion model which properly models the multi-component nature of diffusion in complex gas mixtures such as air, was developed and incorporated into the blunt body model. (7) A model was developed to predict the magnitude and characteristics of the shock wave precursor ahead of vehicles entering the Earth's atmosphere. (8) Since considerable data exists for radiating nonequilibrium flow behind normal shock waves, a normal shock wave version of the blunt body code was developed. (9) By comparing predictions from the models and codes with available normal shock data and the flight data of Fire II, it is believed that the developed flowfield and nonequilibrium radiation models have been essentially validated for engineering applications.
Carlson, Leland A.↗
Seismic Anisotropy of Temperate Ice in Polar Ice Sheets
Here we present a series of simple shear numerical simulations of dynamic recrystallization of two-phase nonlinear viscous materials that represent temperate ice. First, we investigate the effect of the presence of water on the resulting microstructures and, second, how water influences on P wave (V p ) and fast S wave (V s ) velocities. Regardless the water percentage, all simulations evolve from a random fabric to a vertical single maximum. For a purely solid aggregate, the highest V p quickly aligns with the maximum c-axis orientation. At the same time, the maximum c-axis development reduces V s in this orientation. When water is present, the developed maximum c-axis orientation is less intense, which results in lower V p and V s . At high percentage of water, V p does not align with the maximum c-axis orientation. If the bulk modulus of ice is assumed for the water phase (i.e., implying that water is at high pressure), we find a remarkable decrease of V s while V p remains close to the value for purely solid ice. These results suggest that the decrease in V s observed at the base of the ice sheets could be explained by the presence of water at elevated pressure, which would reside in isolated pockets at grain triple junctions. Under these conditions water would not favor sliding between ice grains. However, if we consider that deformation dominates over recrystallization, water pockets get continuously stretched, allowing water films to be located at grain boundaries. This configuration would modify and even overprint the maximum c-axis-dependent orientation and the magnitude of seismic anisotropy.
58 GEOSCIENCES↗
The landscape of disc outflows from black hole–neutron star mergers
ABSTRACT We investigate mass ejection from accretion discs formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH–NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disc configurations expected from BH–NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion discs that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower limit to the fraction of the initial disc mass ejected. We find a nearly linear inverse dependence of this fraction on disc compactness (BH mass over initial disc radius). The dependence is related to the fraction of the disc mass accreted before the ouflow is launched, which depends on the disc position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disc mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disc masses. We estimate the radioactive luminosity from the disc outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH–NS parameter space, the disc outflow contribution is well below the kilonova mass upper limits for GW190814.
79 ASTRONOMY AND ASTROPHYSICS↗
Literature Review of Thermomechanical Analysis, Its Instrument Components and Applications
Thermomechanical analysis (TMA) is used to measure the dimensional changes in a sample as a function of temperature. TMA allows for several modes of analysis, with expansion and compression being the most common. TMA is based on the theory of dilatometry and relies on the fact that polymers are viscoelastic materials, possessing both viscous and elastic properties. TMA yields information about the thermal effects on a polymer such as shrinkage, expansion, softening, melting, glass transition temperature and allows for the calculation of the coefficient of thermal expansion (CTE), which can be useful for many applications of polymers.
36 MATERIALS SCIENCE↗
Literature Review of Thermomechanical Analysis, Its Instrument Components and Applications
Thermomechanical analysis (TMA) is used to measure the dimensional changes in a sample as a function of temperature. TMA allows for several modes of analysis, with expansion and compression being the most common. TMA is based on the theory of dilatometry and relies on the fact that polymers are viscoelastic materials, possessing both viscous and elastic properties. TMA yields information about the thermal effects on a polymer such as shrinkage, expansion, softening, melting, glass transition temperature and allows for the calculation of the coefficient of thermal expansion (CTE), which can be useful for many applications of polymers.
36 MATERIALS SCIENCE↗
Theoretical investigations of plasma processes
System analyses are presented for electrically sustained, collision dominated plasma centrifuges, in which the plasma rotates under the influence of the Lorentz forces resulting from the interaction of the current density fields with an external magnetic field. It is shown that gas discharge centrifuges are technically feasible in which the plasma rotates at speeds up to 1 million cm/sec. The associated centrifugal forces produce a significant spatial isotope separation, which is somewhat perturbed in the viscous boundary layers at the centrifuge walls. The isotope separation effect is the more pronounced. The induced magnetic fields have negligible influence on the plasma rotation if the Hall coefficient is small. In the technical realization of collision dominated plasma centrifuges, a trade-off has to be made between power density and speeds of rotation. The diffusion of sputtered atoms to system surfaces of ion propulsion systems and the deposition of the atoms are treated theoretically by means of a simple model which permits an analytical solution. The problem leads to an inhomogeneous integral equation.
Wilhelm, H. E.↗