Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VISCOSITY”

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 1,099 records · Page 61

VISCOUS ACCRETION OF A POLYTROPIC SELF-GRAVITATING DISK IN THE PRESENCE OF WIND

Self-similar and semi-analytical solutions are found for the height-averaged equations governing the dynamical behavior of a polytropic, self-gravitating disk under the effects of winds around the nascent object. In order to describe the time evolution of the system, we adopt a radius-dependent mass loss rate, then highlight its importance on both the traditional {alpha} and innovative {beta} models of viscosity prescription. In agreement with some other studies, our solutions represent that the Toomre parameter is less than one in most regions on the {beta}-disk, which indicates that in such disks gravitational instabilities can occur at various distances from the central accretor. So, the {beta}-disk model might provide a good explanation of how the planetary systems form. The purpose of the present work is twofold: examining the structure of a disk with wind in comparison to a no-wind solution and seeing whether the adopted viscosity prescription significantly affects the dynamical behavior of the disk-wind system. We also considered the temperature distribution in our disk by a polytropic condition. The solutions imply that, under our boundary conditions, the radial velocity is larger for {alpha}-disks and increases as wind becomes stronger in both viscosity models. Also, we noticed that the disk thickness increases by amplifying the wind or adopting larger values for the polytropic exponent {gamma}. It also may globally decrease if one prescribes a {beta}-model for the viscosity. Moreover, in both viscosity models, the surface density and mass accretion rate diminish as the wind gets stronger or {gamma} increases.

79 ASTRONOMY AND ASTROPHYSICS↗

A SEMI-ANALYTICAL DESCRIPTION FOR THE FORMATION AND GRAVITATIONAL EVOLUTION OF PROTOPLANETARY DISKS

We investigate the formation process of self-gravitating protoplanetary disks in unmagnetized molecular clouds. The angular momentum is redistributed by the action of gravitational torques in the massive disk during its early formation. We develop a simplified one-dimensional accretion disk model that takes into account the infall of gas from the envelope onto the disk and the transfer of angular momentum in the disk with an effective viscosity. First we evaluate the gas accretion rate from the cloud core onto the disk by approximately estimating the effects of gas pressure and gravity acting on the cloud core. We formulate the effective viscosity as a function of the Toomre Q parameter that measures the local gravitational stability of the rotating thin disk. We use a function for viscosity that changes sensitively with Q when the disk is gravitationally unstable. We find a strong self-regulation mechanism in the disk evolution. During the formation stage of protoplanetary disks, the evolution of the surface density does not depend on the other details of the modeling of effective viscosity, such as the prefactor of the viscosity coefficient. Next, to verify our model, we compare the time evolution of the disk calculated with our formulation with that of three-dimensional hydrodynamical simulations. The structures of the resultant disks from the one-dimensional accretion disk model agree well with those of the three-dimensional simulations. Our model is a useful tool for the further modeling of chemistry, radiative transfer, and planet formation in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗

Origin of shear thickening in semidilute wormlike micellar solutions and evidence of elastic turbulence

The origin of shear thickening in an equimolar semidilute wormlike micellar solution of cetylpyridinium chloride and sodium salicylate was investigated in this work by using Couette rheometry, flow visualization, and capillary Rheo-particle image velocimetry. The use of the combined methods allowed the discovery of gradient shear banding flow occurring from a critical shear stress and consisting of two main bands, one isotropic (transparent) of high viscosity and one structured (turbid) of low viscosity. Mechanical rheometry indicated macroscopic shear thinning behavior in the shear banding regime. However, local velocimetry showed that the turbid band increased its viscosity along with the shear stress, even though barely reached the value of the viscosity of the isotropic phase. This shear band is the precursor of shear induced structures that subsequently give rise to the average increase in viscosity or apparent shear thickening of the solution. Further increase in the shear stress promoted the growing of the turbid band across the flow region and led to destabilization of the shear banding flow independently of the type of rheometer used, as well as to vorticity banding in Couette flow. At last, vorticity banding disappeared and the flow developed elastic turbulence with chaotic dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Connecting Structural and Transport Properties of Ionic Liquids with Cationic Oligoether Chains

We used X-ray diffraction and molecular dynamics simulations to probe the structures of two families of ionic liquids containing oligoether tails on the cations. Imidazolium and pyrrolidinium bis(trifluoromethylsulfonyl)amide ILs with side chains ranging from 4 to 10 atoms in length, including both linear alkyl and oligo-ethylene oxide tails, were prepared. Furthermore, their physical properties, such as viscosity, conductivity and thermal profile, were measured and compared for systematic trends. Consistent with earlier literature, a single ether substituent substantially decreases the viscosity of pyrrolidinium and imidazolium ILs compared to their alkyl congeners. Remarkably, as the number of ether units in the pyrrolidinium ILs increases there is hardly any increase in the viscosity, in contrast to alkylpyrrolidinium ILs where the viscosity increases steadily with chain length. Viscosities of imidazolium ether ILs increase with chain length but always remain well below their alkyl congeners. To complement the experimentally determined properties, molecular dynamics simulations were run on the two ILs with the longest ether chains. Our results point to specific aspects that could be useful for researchers designing ILs for specific applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Temperature Effects on Droplet Oscillation Decay with Application to Fuel Property Measurement

Observation of oscillation decay in droplets has been shown to be an effective approach in determining physical properties such as viscosity and surface tension for emerging biofuels using μl quantities. Herein this work extends the approach to higher temperatures relevant to fuel injection conditions for internal combustion engines. Experiments are conducted that use high-resolution strobed imaging of moving, heated, μm-sized, fuel droplets to capture shape oscillation decay through image processing and analysis. Two fuels are investigated, iso-butanol and a primary reference fuel (PRF 84), which is a mixture of iso-octane and n-heptane. A piezoelectric droplet generator is used to generate a continuous train of single droplets, which are given an initial perturbation and subsequently undergo damped oscillations, captured using strobed imaging. Surface tension and viscosity calculated at four different temperatures ranging from 30°C-55°C using the frequency and decay time associated with the fundamental mode are found to be within ~ 10% of reference values obtained from the literature. Complementary numerical simulations are performed that utilize a volume-of-fluid approach to track transient droplet oscillation phenomena along with heat and mass transfer in a 2D axisymmetric domain. Simulations, where droplet size and temperature can be independently varied, capture an expected decrease in oscillation frequency and increase in decay time, with increase in fuel temperature. The simulations are further used to study the relative contribution to deviations in surface tension and viscosity predictions due to heat and mass transfer from the droplet, as well as viscous effects violating the inviscid flow assumption in the droplet oscillation theory. Mass loss effects are found to be negligible. Temperature change due to heat transfer is found to have the next highest sensitivity, particularly for the more volatile PRF 84, which has a lower viscosity and associated Ohnesorge number. For isobutanol, which has a higher viscosity and Ohnesorge number, viscous effects contribute the most to deviation in fuel property predictions.

09 BIOMASS FUELS↗

Comparison of observed rheological properties of hard wheat flour dough with predictions of the Giesekus-Leonov, White-Metzner and Phan-Thien Tanner models

The measured rheological behavior of hard wheat flour dough was predicted using three nonlinear differential viscoelastic models. The Phan-Thien Tanner model gave good zero shear viscosity prediction, but overpredicted the shear viscosity at higher shear rates and the transient and extensional properties. The Giesekus-Leonov model gave similar predictions to the Phan-Thien Tanner model, but the extensional viscosity prediction showed extension thickening. Using high values of the mobility factor, extension thinning behavior was observed but the predictions were not satisfactory. The White-Metzner model gave good predictions of the steady shear viscosity and the first normal stress coefficient but it was unable to predict the uniaxial extensional viscosity as it exhibited asymptotic behavior in the tested extensional rates. It also predicted the transient shear properties with moderate accuracy in the transient phase, but very well at higher times, compared to the Phan-Thien Tanner model and the Giesekus-Leonov model. None of the models predicted all observed data consistently well. Overall the White-Metzner model appeared to make the best predictions of all the observed data.

Non-NASA Center↗

Drug marker absorption in relation to pellet size, gastric motility and viscous meals in humans

PURPOSE: The objective of this study was to evaluate drug marker absorption in relation to the gastric emptying (GE) of 0.7 mm and 3.6 mm enteric coated pellets as a function of viscosity and the underlying gastric motility. METHODS: Twelve subjects were evaluated in a 3-way crossover study. 0.7 mm caffeine and 3.6 mm acetaminophen enteric coated pellets were concurrently administered with a viscous caloric meal at the levels of 4000, 6000 and 8000 cP. Gastric motility was simultaneously measured with antral manometry and compared to time events in the plasma profiles of the drug markers. RESULTS: Caffeine, from the 0.7 mm pellets, was observed significantly earlier in the plasma than acetaminophen, from the 3.6 mm pellets, at all levels of viscosity. Motility related size differentiated GE was consistently observed at all viscosity levels, however, less variability was observed with the 4000 cP meal. Specifically, the onset of absorption from the of 3.6 mm pellets correlated with the onset of Phase II fasted state contractions (r = 0.929, p < 0.01). CONCLUSIONS: The timeframe of drug marker absorption and the onset of motility events were not altered within the range of viscosities evaluated. Rather, the differences in drug marker profiles from the non-digestible solids were most likely the result of the interaction between viscosity and motility influencing antral flow dynamics. The administration of the two sizes of pellets and a viscous caloric meal with subsequent monitoring of drug marker profiles is useful as a reference to assess the influence of motility patterns on the absorption profile of orally administered agents.

NASA Center JSC↗

Evaluation of the New Boron Model of the TRACE Code

System codes such as TRACE and RELAP5 have the option of adding boron to be dissolved in the water. In principle, the density and the viscosity of the water are not affected when the amounts of dissolved boron are very small. However, if the boron concentration is significant, water density and viscosity will be affected. Version 5.890 of the TRACE code has a new boron model with the option of modifying the water density and viscosity to account for the presence of dissolved boron. The input variable ISOLUT is used for this purpose. If ISOLUT=0, no boron is dissolved in the water. If ISOLUT=1, boron is dissolved in the water, and the liquid water density and viscosity are not modified by the dissolved boron. If ISOLUT=2, boron is dissolved in the water, and the liquid water density and viscosity are modified by the dissolved boron. This paper evaluates some features of this new boron model in the TRACE code.

42 ENGINEERING↗

Laser‐Sintered Silver Metallization for Silicon Heterojunction Photovoltaic Cells

Herein, a novel metallization technique is reported for crystalline silicon heterojunction (SHJ) solar cells in which silver (Ag) fingers are printed on the SHJ substrates by dispensing Ag nanoparticle‐based inks through a needle and then sintered with a continuous‐wave carbon dioxide (CO 2 ) laser. The impact of the Ag ink viscosity on the line quality and the line resistance is investigated on three Ag inks with different viscosities. Increasing ink viscosity yields higher Ag contact heights, larger aspect ratios, and lower line resistance values. The Ag line height increases from less than a micrometer to ≈18.62 ± 3.48 μm with the increasing viscosity. Photoluminescence imaging shows that the low‐resistance Ag metal contacts obtained do not result in any passivation damage of the SHJ substrate. This is because the wavelength of light emitted from the CO 2 laser (i.e., 10.6 μm) leads to optical absorption in the Ag, but this light is effectively transparent to the transparent conductive oxide film, amorphous silicon films, and crystalline silicon substrate. Bulk resistivity values as low as 6.5 μΩ cm are obtained for the laser‐sintered Ag contact and printed using the Ag ink with the highest viscosity in this work.

Mousumi, Jannatul Ferdous↗

Optimal Wetting Angles in Lattice Boltzmann Simulations of Viscous Fingering

We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios M=ν r /ν b ∈[0.001,100] and wetting angles θ w ∈[0°,180°]. As expected, at low-viscosity ratios M=ν r /ν b $\ll$1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios M$\gg$1, we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.

Engineering↗

Representing ductile damage with the dual domain material point method

In this study, we incorporate a ductile damage material model into a computational framework based on the Dual Domain Material Point (DDMP) method. As an example, simulations of a flyer plate experiment involving ductile void growth and material failure are performed. The results are compared with experiments performed on high purity tantalum. We also compare the numerical results obtained from the DDMP method with those obtained from the traditional Material Point Method (MPM). Effects of an overstress model, artificial viscosity, and physical viscosity are investigated. Our results show that a physical bulk viscosity and overstress model are important in this impact and failure problem, while physical shear viscosity and artificial shock viscosity have negligible effects. A simple numerical procedure with guaranteed convergence is introduced to solve for the equilibrium plastic state from the ductile damage model.

36 MATERIALS SCIENCE↗

Comprehensive study of the interactions between the critical dimensionless numbers associated with multiphase flow in 3D porous media

Multiphase flow in porous media is of signifcant interest in many engineering applications, such as geologic carbon sequestration, enhanced oil recovery, and groundwater contamination and remediation. In order to advance the fundamental understanding of multiphase flow in complex three-dimensional (3D) porous media, the interactions between the critical dimensionless numbers, including the contact angle, viscosity ratio, and capillary (Ca) number, were explored using X-ray micro-computed tomography (micro-CT) scanning and lattice Boltzmann (LB) modeling. In this study, the 3D pore structure information was extracted from micro-CT images and then used as interior boundary conditions of flow modeling in a pore-scale LB simulator to simulate multiphase flow within the pore space. A Berea sandstone sample was scanned and then two-phase flow LB simulations were performed based on the micro-CT images. The LB-simulated water/CO 2 distributions agreed well with the micro-CT scanned images. Simulation findings showed that a decreasing contact angle causes a decrease in wetting-fluid relative permeability and an increase in non-wetting fluid relative permeability. A rising Ca number increases both wetting and non-wetting fluid relative permeabilities. An increasing viscosity ratio (the ratio of non-wetting fluid viscosity to wetting fluid viscosity) facilitates the increase of non-wetting fluid relative permeability and mitigates the reduction of wetting fluid relative permeability, when the contact angle decreases continuously. The primary novel finding of this study is that the viscosity ratio affects the rate of change of the relative permeability curves for both phases when the contact angle changes continuously. To date, we believe it is the first time that comprehensive interactions between these dimensionless numbers are demonstrated in a sandstone sample based on real 3D structures. We also investigated the role of the changes of flow direction and sample location on relative permeability curves. Simulation results showed that the change in non-wetting fluid relative permeability was larger when the flow direction was switched from vertical to horizontal, which indicated that there was stronger anisotropy in larger pore networks that were primarily occupied by the non-wetting fluid.

42 ENGINEERING↗

Machine learning from RANS and LES to inform coarse grid simulations

Nuclear system thermal hydraulic analysis has historically relied on computationally inexpensive 1D codes. However, such tools are unable to capture multiscale multidimensional effects in large nuclear reactor enclosures. On the other hand, simulations with higher fidelity can be too expensive for such purposes. One of the ways to reduce computational cost is to perform simulations on a coarse grid, which, unfortunately, introduces large discretization errors. In this paper, two high-to-low data-driven approaches are investigated: (1) a coarse grid turbulence model to predict eddy viscosity and (2) correction of errors in coarse grid velocity fields. The approaches aim to reduce grid- and turbulence model-induced errors in coarse grid Reynolds-averaged Navier–Stokes (RANS) simulations. Two sources of high-fidelity data, RANS and large eddy simulations (LES), are explored. To extract the eddy viscosity from the LES data, an inverse optimization problem is solved. However, the LES eddy viscosity is shown to be comparable to the RANS eddy viscosity in terms of error reduction. Therefore, the directly available RANS eddy viscosity was used to develop a coarse grid data-driven turbulence model. Additionally, error correction in velocity is used to reduce the remaining uncertainties and bring the results closer to reality. In conclusion, the performance of the frameworks is demonstrated for a scaled upper plenum of a gas-cooled reactor facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The role of diffusion-driven pure climb creep on the rheology of bridgmanite under lower mantle conditions

The viscosity of Earth’s lower mantle is poorly constrained due to the lack of knowledge on some fundamental variables that affect the deformation behaviour of its main mineral phases. This study focuses on bridgmanite, the main lower mantle constituent, and assesses its rheology by developing an approach based on mineral physics. Following and revising the recent advances in this field, pure climb creep controlled by diffusion is identified as the key mechanism driving deformation in bridgmanite. The strain rates of this phase under lower mantle pressures, temperatures and stresses are thus calculated by constraining diffusion and implementing a creep theoretical model. The viscosity of MgSiO3 bridgmanite resulting from pure climb creep is consequently evaluated and compared with the viscosity profiles available from the literature. We show that the inferred variability of viscosity in these profiles can be fully accounted for with the chosen variables of our calculation, e.g., diffusion coefficients, vacancy concentrations and applied stresses. A refinement of these variables is advocated in order to further constrain viscosity and match the observables.

58 GEOSCIENCES↗

Sunlight can turn smoldering pine wood smoke into a glass

Wildfires inject biomass burning organic aerosols (BBOA) into the atmosphere. During their lifetimes of weeks to months, they are exposed to ultraviolet (UV) irradiation. Viscosity and phase behaviour are essential properties for understanding their chemical and climate impacts. Here, we quantify changes in viscosity and phase behavior after UV exposure. After an atmospheric equivalent of 8.7 days of boundary layer UV exposure, BBOA develop a highly viscous (glassy) outer phase with a viscosity at least five orders of magnitude higher than unaged BBOA, which persists up to at least 58% relative humidity. High-resolution mass spectrometry indicates that UV-aging increases oxidation and molecular weight. Using our viscosity results, we predict that UV-aged BBOA are frequently glassy above ~2.5 km and can have viscosities up to eight orders of magnitude higher than unaged particles in some atmospheric regions, which may influence lifetimes of pollutants and brown carbon, and impact stratospheric ozone chemistry.

Golay, Zoe↗

The effect of spin induced magnetization on Jeans instability of viscous and resistive quantum plasma

The effect of spin induced magnetization and electrical resistivity incorporating the viscosity of the medium is examined on the Jeans instability of quantum magnetoplasma. Formulation of the system is done by using the quantum magnetohydrodynamic model. The analysis of the problem is carried out by normal mode analysis theory. The general dispersion relation is derived from set of perturbed equations to analyse the growth rate and condition of self-gravitational Jeans instability. To discuss the influence of resistivity, magnetization, and viscosity parameters on Jeans instability, the general dispersion relation is reduced for both transverse and longitudinal mode of propagations. In the case of transverse propagation, the gravitating mode is found to be affected by the viscosity, magnetization, resistivity, and magnetic field strength whereas Jeans criterion of instability is modified by the magnetization and quantum parameter. In the longitudinal mode of propagation, the gravitating mode is found to be modified due to the viscosity and quantum correction in which the Jeans condition of instability is influenced only by quantum parameter. The other non-gravitating Alfven mode in this direction is affected by finite electrical resistivity, spin induced magnetization, and viscosity. The numerical study for the growth rate of Jeans instability is carried out for both in the transverse and longitudinal direction of propagation to the magnetic field. The effect of various parameters on the growth rate of Jeans instability in quantum plasma is analysed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the accuracy of compressibility transformations

This study highlights the importance of satisfying the eddy viscosity equivalence below the logarithmic layer, to deriving accurate compressibility transformations. First, we analyze the ability of known transformations to satisfy the eddy viscosity equivalence and show that the accuracy of these transformations is strongly dependent on this ability. Second, in a step-by-step manner, we devise new transformations that satisfy this hypothesis. An approach based on curve fitting of the incompressible Direct Numerical Simulation data for eddy viscosity profiles below the logarithmic layer provides an extremely accurate transformation, which motivates self-contained methods, making use of mixing length formulas in the inner region. It is shown that the accuracy of existing transformations can be significantly improved by applying these ideas, below the logarithmic layer. Motivated by the effectiveness of the formulations derived from eddy viscosity equivalence, we introduce a new integral transformation based on Reynolds number equivalence between compressible and incompressible flows. This approach is based on defining a new compressible velocity scale, which affects the accuracy of transformations. Several choices for the velocity scale are tested, and in each attempt, it is shown that the eddy viscosity equivalence plays a very important role for the accuracy of compressibility transformations.

42 ENGINEERING↗

Feasibility of observing small differences in friction mean effective pressure between different lubricating oil formations using small, single-cylinder motored engine rig

Here, the feasibility of using a motored single-cylinder 517 cc diesel engine to observe small frictional differences between oil formulations is investigated. Friction mean effective pressure (FMEP) is measured and compared for an SAE 10W-30 and an SAE 5W-20 oil in three stages of production: base oil, commercial oil without a friction and wear reducing additive, and fully formulated commercial oil. In addition, a commercial SAE 5W-30 engine oil is investigated. Friction mean effective pressure is plotted versus oil dynamic viscosity to compare the lubricant FMEP at a given viscosity. Linear regressions and average friction mean effective pressure are used as a secondary means of comparing FMEP for the various oil formulations. Differences between the oils are observed with the base oil having higher friction at a given viscosity but a lower average FMEP due to the temperature distribution of the test and lower viscosities reached by the base oil. The commercial oil is shown to have both a higher FMEP at a given viscosity and a higher average FMEP than the commercial oil without a friction and wear reducing additive. The increase in friction for the oil without a friction and wear reduction additive indicates that the operational regime of the engine may be out of the bounds of the optimal regime for the additive or that the additive is more optimized for wear reduction. Results show that it is feasible to observe small differences in FMEP between lubricating oil formulations using a small, single-cylinder motored engine.

33 ADVANCED PROPULSION SYSTEMS↗