Engineering Papers⌕ Search

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.

At least 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Experimental and Theoretical Studies of Pulsating Turbulent Flow

The objective of this investigation was to study the effects of small amplitude sinusoidal pulsations on fully developed turbulent flow in a tube from both experimental and theoretical viewpoints. Theoretical models for the macroscopic behavior of pulsating turbulent tube flow were developed for the two cases of very low and very high pulsation frequencies. The models are based on assumptions of quasi-steady and frozen eddy viscosity flow behavior, respectively. The models successfully predict unsteady velocity profiles, thereby supporting the currently proposed definitions of frequency regimes in pulsating turbulent flow. Experimental measurements were made of the time-dependent pressure drop and velocity profiles over the range of frequency-to-Reynolds number ratios from 0.0095 to 0.24. The two macroscopic models developed in this study predict unsteady velocity profiles which are in moderately good agreement with the experiments in their respective frequency regimes, and a previously developed quasi-steady model is found to predict experimental velocity profiles well in both the quasisteady and the frozen eddy viscosity frequency regimes. The effect of flow pulsations on the dissipation of turbulence energy in the vicinity of the wall was measured in the lower transition frequency regime. The long-time averaged dissipation was observed to be unchanged from the steady flow dissipation, within the accuracy of the experiment. A theoretical model of the periodic viscous sublayer was also developed and applied to pulsating flow in a tube, in order to investigate the effects of flow pulsations on the rate of production of turbulence in the region of the wall. The periodic viscous sublayer model predicts sublayer growth periods in steady flow which agree with the published experimental data. When the model is applied to pulsating flow, the response of the sublayer growth period falls into three frequency regimes, the parameters of which are in approximate agreement with the frequency regimes which are defined on the basis of macroscopic flow behavior. The sublayer renewal cycle exhibits quasi-steady flow behavior when the sublayer growth period is much less than the pulsation period, transition behavior when these two periods are approximately equal, and frozen eddy viscosity behavior when the sublayer period is much longer than the pulsation period. The effect of the sublayer growth and renewal cycle on the level of turbulence was investigated by two methods. The velocity fluctuations seen by a point velocity probe located close to the wall were predicted from the model in one method and the rate of turbulence production was estimated from the frequency of sublayer renewal events in the other.

Kingston, G. C.↗
NASA NTRS1993

A numerical study of the effect of geometry variation, turbulence models, and dissipation on the flow past control surfaces

The steady 3D viscous flow past the ONERA M6 wing and a slender delta wing-body with trailing edge control surfaces has been computed. A cell-centered finite-volume Navier-Stokes patched zonal method has been used for the numerical simulation. Both diagonalized and LUSGS schemes have been implemented. Besides the standard nonplanar zonal interfacing techniques, a new virtual zone capability has been employed. For code validation, the transonic flow past the ONERA M5 wing is calculated for angles-of-attack of 3.06 deg and 5.06 deg and compared with the available experiments. The wing-body computational results are compared with experimental data for both trailing-edge flaps deflected. The experimental flow conditions are M subinfinity = 0.4, a turbulent Reynolds number of 5.41 million based on a mean aerodynamic chord of 25.959 inches, adiabatic wall, and angles-of-attack varying from 0 deg to 23.85 deg. The computational results are presented for the 23.85 deg angle-of-attack case. The effects of the base flow due to a model sting, the varying second and fourth order numerical dissipation, and the turbulence model are all considered.

Chaussee, Denny S.↗
NASA NTRS2017

Baseline Experimental Results on the Effect of Oil Temperature on Shrouded Meshed Spur Gear Windage Power Loss

Rotorcraft gearbox efficiencies are reduced at increased surface speeds due to viscous and impingement drag on the gear teeth. This windage power loss can affect overall mission range, payload, and frequency of transmission maintenance. Experimental and analytical studies on shrouding for single gears have shown it be potentially effective in mitigating windage power loss. Efficiency studies on unshrouded meshed gears have shown the effect of speed, oil viscosity, temperature, load, lubrication scheme, etc. on gear windage power loss. The open literature does not cite data on shrouded meshed spur gears. Gear windage power loss test results are presented on shrouded meshed spur gears at elevated oil inlet temperatures and constant oil pressure both with and without shrouding. Shroud effectiveness is compared at four oil inlet temperatures. The results are compared to the available literature and follow-up work is outlined.

windage↗
NASA NTRS2017

Baseline Experimental Results on the Effect of Oil Temperature on Shrouded Meshed Spur Gear Windage Power Loss

Rotorcraft gearbox efficiencies are reduced at increased surface speeds due to viscous and impingement drag on the gear teeth. This windage power loss can affect overall mission range, payload, and frequency of transmission maintenance. Experimental and analytical studies on shrouding for single gears have shown it to be potentially effective in mitigating windage power loss. Efficiency studies on unshrouded meshed gears have shown the effect of speed, oil viscosity, temperature, load, lubrication scheme, etc. on gear windage power loss. The open literature does not contain experimental test data on shrouded meshed spur gears. Gear windage power loss test results are presented on shrouded meshed spur gears at elevated oil inlet temperatures and constant oil pressure both with and without shrouding. Shroud effectiveness is compared at four oil inlet temperatures. The results are compared to the available literature and follow-up work is outlined.

gear↗
NASA NTRS1985

Tidal disruption of dissipative planetesimals

A self-consistent numerical model is developed for the tidal disruption of a solid planetesimal. The planetesimal is treated as a highly viscous, slightly compressible fluid whose disturbed parts are an inviscid, pressureless fluid undergoing distortion and disruption. The distortions were constrained to being symmetrical above and below the equatorial plane. The tidal potential is expanded in terms of Legendre polynomials, which eliminates the center of mass acceleration effects, permitting definition of equations of motion in a noninertial frame. Consideration is given to viscous dissipation and to characteristics of the solid-atmosphere boundary. The model is applied to sample cases in one, two and three dimensions.

Mizuno, H.↗
NASA NTRS1994

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗
NASA NTRS1985

Solution of viscous transonic flow over wings

Since the calculations of plane steady transonic flows conducted by Murman and Cole (1971), steady progress has been made with respect to the computation of inviscid transonic flows. It has been found that it is inadequate to consider practical wing design at transonic speeds without considering the effects of viscosity and turbulence. The present study has the objective to develop a zonal method for viscous transonic flow over three-dimensional (3-D) wings. The employed approach follows closely the viscous/inviscid interaction techniques discussed by Melnik et al. (1983) for viscous flow about airfoils. Attention is given to inviscid flow equations and boundary conditions, the solution of 3-D boundary layer and wake, an iterative solution to viscid-inviscid interaction analysis, and results obtained for a transonic cruise wing of transport type.

Chow, R. R.↗
NASA NTRS1974

Reacting viscous-shock-layer solutions with multicomponent diffusion and mass injection

Numerical solutions of the viscous-shock-layer equation where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium are presented. Also the effects of the diffusion model, surface catalysis, and mass injection on surface transport and flow parameters are considered. The flow is treated as a mixture of five inert and thermally perfect species. The viscous-shock-layer equations are solved by using an implicit-difference scheme. All calculations are for hyperboloids with included angles of 20 and 45. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented to show the effects of the chemical models; diffusion models; surface catalysis; and mass injection of air on heat transfer; skin friction; shock standoff distance; wall pressure distribution; and tangential velocity, temperature, and species profiles. The results show that an equilibrium analysis can substantially overpredict the heat-transfer rates for flow conditions experienced by earth-orbital entry vehicles. Moreover, at such conditions surface catalysis significantly influences heat-transfer and flow-field properties. If a binary rather than a multicomponent diffusion model is assumed, negligible errors in most flow properties result. Quantitative results are presented that show the effect of mass injection on flow properties within and downstream of the injection region.

Moss, J. N.↗
NASA NTRS1990

An investigation of design optimization using a 2-D viscous flow code with multigrid

Computational fluid dynamics (CFD) codes have advanced to the point where they are effective analytical tools for solving flow fields around complex geometries. There is also a need for their use as a design tool to find optimum aerodynamic shapes. In the area of design, however, a difficulty arises due to the large amount of computer resources required by these codes. It is desired to streamline the design process so that a large number of design options and constraints can be investigated without overloading the system. There are several techniques which have been proposed to help streamline the design process. The feasibility of one of these techniques is investigated. The technique under consideration is the interaction of the geometry change with the flow calculation. The problem of finding the value of camber which maximizes the ratio of lift over drag for a particular airfoil is considered. In order to test out this technique, a particular optimization problem was tried. A NACA 0012 airfoil was considered at free stream Mach number of 0.5 with a zero angle of attack. Camber was added to the mean line of the airfoil. The goal was to find the value of camber for which the ratio of lift over drag is a maximum. The flow code used was FLOMGE which is a two dimensional viscous flow solver which uses multigrid to speed up convergence. A hyperbolic grid generation program was used to construct the grid for each value of camber.

Doria, Michael L.↗
NASA NTRS1979

Lubricant limiting shear stress effect on EHD film thickness

A Grubin-like EHD inlet analysis utilizing a non-linear viscous fluid model with a limiting shear stress is reported. The shear rheological equation requirs only a low shear stress viscosity and the limiting shear stress both functions of pressure. Values employed for these properties are taken from measurements on typical lubricants. Reductions of EHD film thickness are found to be up to 40 percent compared with the standard Grubin prediction for typical operating conditions. Slide-roll ratio, limiting shear stress dependence on pressure, and atmospheric pressure value of limiting shear stress are new variables required to determine film thickness with the first two being more important than the last. The EHD film thickness is reduced by increasing slide-roll ratio and/or decreasing the pressure dependence of the limiting shear stress.

Gecim, B.↗
NASA NTRS1984

Impact experiments in viscous fluid media

Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. The resulting models, for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-veocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studied. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SHSV in the vicinity of the low-velocity zone. Near the East Pacific Rise, however, SVSH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of about 50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures A new model for Average Earth is obtained which includes a thin LID. This model extends the fit of a P.R.E.M. type model to shorter period surface waves.

Greeley, R.↗
NASA NTRS1984

Host turbine heat transfer overview

Improved methods of predicting airfoil local metal temperatures require advances in the understanding of the physics and methods of analytically predicting the following four aerothermal loads: hot gas flow over airfoils, heat transfer rates on the gas-side of airfoils, cooling air flow inside airfoils, and heat transfer rates on the coolant-side of airfoils. A systematic building block research approach is being pursued to investigate these four areas of concern from both the experimental and analytical sides. Experimental approaches being pursued start with fundamental experiments using simple shapes and flat plates in wind tunnels, progress to more realistic cold and hot cascade tests using airfoils, continue to progress in large low-speed rigs and turbines and warm turbines, and finally, combine all the interactive effects in tests using real engines or real engine type turbine rigs. Analytical approaches being pursued also build from relatively simple steady two dimensional inviscid flow and boundary layer heat transfer codes to more advanced steady two and three dimensional viscous flow and heat transfer codes. These advanced codes provide more physics to model better the interactive effects and the true real-engine environment.

Rohde, J. E.↗
DOE OSTI2025

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗
NASA NTRS1983

Micromechanically Based Constitutive Relations for Polycrystalline Solids

A basic method to estimate the overall mechanical response of solids which contain periodically distributed defects is presented. The method estimates the shape and growth pattern of voids periodically distributed over the grain boundaries in a viscous matrix. The relaxed moduli are obtained for a polycrytalline solid that undergoes relaxation by grain boundary sliding which accounts for the interaction effects. The overall inelastic nonlinear response at elevated temperatures in terms of a model which considers nonlinear power law creep within the grains, and linear viscous flow in the grain boundaries is discussed.

Nemat-Nasser, S.↗
DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗
DOE OSTI2020

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗
DOE OSTI2022

On the Injection Scale of the Turbulence in the Partially Ionized Very Local Interstellar Medium

The cascade of magnetohydrodynamic (MHD) turbulence is subject to ion–neutral collisional damping and neutral viscous damping in the partially ionized interstellar medium. By examining the damping effects in the warm and partially ionized local interstellar medium, we find that the interstellar turbulence is damped by neutral viscosity at ~261 au and cannot account for the turbulent magnetic fluctuations detected by Voyager 1 and 2. The MHD turbulence measured by Voyager in the very local interstellar medium (VLISM) should be locally injected in the regime where ions are decoupled from neutrals for its cascade to survive the damping effects. With the imposed ion–neutral decoupling condition and the strong turbulence condition for the observed Kolmogorov magnetic energy spectrum, we find that the turbulence in the VLISM is sub-Alfvénic, and its largest possible injection scale is ~194 au.

79 ASTRONOMY AND ASTROPHYSICS↗