Energy flow and energy dissipation in a free surface.

NASA Astrophysics Data System (ADS)

Goldburg, Walter; Cressman, John

2005-11-01

Turbulent flows on a free surface are strongly compressible [1] and do not conserve energy in the absence of viscosity as bulk fluids do. Despite violation of assumptions essential to Kolmogorov's theory of 1941 (K41) [2, 3], surface flows show strong agreement with Kolmogorov scaling, though intermittency is larger there. Steady state turbulence is generated in a tank of water, and the spatially averaged energy flux is measured from the four-fifth's law at each instant of time. Likewise, the energy dissipation rate as measured from velocity gradients is also a random variable in this experiment. The energy flux - dissipation rate cross-correlation is measured to be correlated in incompressible bulk flows, but strongly anti-correlated on the surface. We argue that the reason for this discrepancy between surface and bulk flows is due to compressible effects present on the surface. [1] J. R. Cressman, J. Davoudi, W. I. Goldburg, and J. Schumacher, New Journal of Physics, 6, 53, 2004. [2] U. Frisch. Turbulence: The legacy of A. N. Kolmogorov, Cambridge University Press, Cambridge, 1995. [3] A. N. Kolmogorov, Doklady Akad. Nauk SSSR, 32, 16, 1941.

An experimental investigation of compressible three-dimensional boundary layer flow in annular diffusers

NASA Technical Reports Server (NTRS)

Om, Deepak; Childs, Morris E.

1987-01-01

An experimental study is described in which detailed wall pressure measurements have been obtained for compressible three-dimensional unseparated boundary layer flow in annular diffusers with and without normal shock waves. Detailed mean flow-field data were also obtained for the diffuser flow without a shock wave. Two diffuser flows with shock waves were investigated. In one case, the normal shock existed over the complete annulus whereas in the second case, the shock existed over a part of the annulus. The data obtained can be used to validate computational codes for predicting such flow fields. The details of the flow field without the shock wave show flow reversal in the circumferential direction on both inner and outer surfaces. However, there is a lag in the flow reversal between the inner nad the outer surfaces. This is an interesting feature of this flow and should be a good test for the computational codes.

Study of compressible flow through a rectangular-to-semiannular transition duct

NASA Technical Reports Server (NTRS)

Foster, Jeffry; Okiishi, Theodore H.; Wendt, Bruce J.; Reichert, Bruce A.

1995-01-01

Detailed flow field measurements are presented for compressible flow through a diffusing rectangular-to-semiannular transition duct. Comparisons are made with published computational results for flow through the duct. Three-dimensional velocity vectors and total pressures were measured at the exit plane of the diffuser model. The inlet flow was also measured. These measurements are made using calibrated five-hole probes. Surface oil flow visualization and surface static pressure data were also taken. The study was conducted with an inlet Mach number of 0.786. The diffuser Reynolds based on the inlet centerline velocity and the exit diameter of the diffuser was 3,200,000. Comparison of the measured data with previously published computational results are made. Data demonstrating the ability of vortex generators to reduce flow separation and circumferential distortion is also presented.

On predicting receptivity to surface roughness in a compressible infinite swept wing boundary layer

NASA Astrophysics Data System (ADS)

Thomas, Christian; Mughal, Shahid; Ashworth, Richard

2017-03-01

The receptivity of crossflow disturbances on an infinite swept wing is investigated using solutions of the adjoint linearised Navier-Stokes equations. The adjoint based method for predicting the magnitude of stationary disturbances generated by randomly distributed surface roughness is described, with the analysis extended to include both surface curvature and compressible flow effects. Receptivity is predicted for a broad spectrum of spanwise wavenumbers, variable freestream Reynolds numbers, and subsonic Mach numbers. Curvature is found to play a significant role in the receptivity calculations, while compressible flow effects are only found to marginally affect the initial size of the crossflow instability. A Monte Carlo type analysis is undertaken to establish the mean amplitude and variance of crossflow disturbances generated by the randomly distributed surface roughness. Mean amplitudes are determined for a range of flow parameters that are maximised for roughness distributions containing a broad spectrum of roughness wavelengths, including those that are most effective in generating stationary crossflow disturbances. A control mechanism is then developed where the short scale roughness wavelengths are damped, leading to significant reductions in the receptivity amplitude.

On the theory of oscillating airfoils of finite span in subsonic compressible flow

NASA Technical Reports Server (NTRS)

Reissner, Eric

1950-01-01

The problem of oscillating lifting surface of finite span in subsonic compressible flow is reduced to an integral equation. The kernel of the integral equation is approximated by a simpler expression, on the basis of the assumption of sufficiently large aspect ratio. With this approximation the double integral occurring in the formulation of the problem is reduced to two single integrals, one of which is taken over the chord and the other over the span of the lifting surface. On the basis of this reduction the three-dimensional problem appears separated into two two-dimensional problems, one of them being effectively the problem of two-dimensional flow and the other being the problem of spanwise circulation distribution. Earlier results concerning the oscillating lifting surface of finite span in incompressible flow are contained in the present more general results.

Further Investigations of Control Surface Seals for the X-38 Re-Entry Vehicle

NASA Technical Reports Server (NTRS)

Dunlap, Patrick H., Jr.; Steinetz, Bruce M.; Curry, Donald M.; Newquist, Charles W.; Verzemnieks, Juris

2001-01-01

NASA is currently developing the X-38 vehicle that will be used to demonstrate the technologies required for a potential crew return vehicle (CRV) for the International Space Station. This vehicle would serve both as an ambulance for medical emergencies and as an evacuation vehicle for the Space Station. Control surfaces on the X-38 (body flaps and rudder/fin assemblies) require high temperature seals to limit hot gas ingestion and transfer of heat to underlying low-temperature structures to prevent over-temperature of these structures and possible loss of the vehicle. NASAs Johnson Space Center (JSC) and Glenn Research Center (GRC) are working together to develop and evaluate seals for these control surfaces. This paper presents results for compression. flow, scrub, and arc jet tests conducted on the baseline X-38 rudder/fin seal design. Room temperature seal compression tests were performed at low compression levels to determine load versus linear compression, preload. contact area, stiffness. and resiliency characteristics under low load conditions. For all compression levels that were tested, unit loads and contact pressures for the seals were below the 5 lb/in. and 10 psi limits required to limit the loads on the adjoining Shuttle thermal tiles that the seals will contact. Flow rates through an unloaded (i.e. 0% compression) double arrangement were twice those of a double seal compressed to the 20% design compression level. The seals survived an ambient temperature 1000 cycle scrub test over relatively rough Shuttle tile surfaces. The seals were able to disengage and re-engage the edges of the rub surface tiles while being scrubbed over them. Arc jet tests were performed to experimentally determine anticipated seal temperatures for representative flow boundary conditions (pressures and temperatures) under simulated vehicle re-entry conditions. Installation of a single seat in the gap of the test fixture caused a large temperature drop (1710 F) across the seal location as compared to an open gap condition (140 F) confirming the need for seals in the rudder/fin gap location. The seal acted as an effective thermal barrier limiting heat convection through the seal gap and minimizing temperature increases downstream of the seal during maximum heating conditions.

Calculation of three-dimensional compressible laminar and turbulent boundary flows. Three-dimensional compressible boundary layers of reacting gases over realistic configurations

NASA Technical Reports Server (NTRS)

Kendall, R. M.; Bonnett, W. S.; Nardo, C. T.; Abbett, M. J.

1975-01-01

A three-dimensional boundary-layer code was developed for particular application to realistic hypersonic aircraft. It is very general and can be applied to a wide variety of boundary-layer flows. Laminar, transitional, and fully turbulent flows of compressible, reacting gases are efficiently calculated by use of the code. A body-oriented orthogonal coordinate system is used for the calculation and the user has complete freedom in specifying the coordinate system within the restrictions that one coordinate must be normal to the surface and the three coordinates must be mutually orthogonal.

Computational analysis of hypersonic flows past elliptic-cone waveriders

NASA Technical Reports Server (NTRS)

Yoon, Bok-Hyun; Rasmussen, Maurice L.

1991-01-01

A comprehensive study for the inviscid numerical calculation of the hypersonic flow past a class of elliptic-cone derived waveriders is presented. The theoretical background associated with hypersonic small-disturbance theory (HSDT) is reviewed. Several approximation formulas for the waverider compression surface are established. A CFD algorithm is used to calculate flow fields for the on-design case and a variety of off-design cases. The results are compared with HSDT, experiment, and other available CFD results. For the waverider shape used in previous investigations, the bow shock for the on-design condition stands off from the leading-edge tip of the waverider. It was found that this occurs because the tip was too thick according to the approximating shape formula that was used to describe the compression surface. When this was corrected, the bow shock became closer to attached as it should be. At Mach numbers greater than the design condition, a lambda-shock configuration develops near the tip of the compression surface. At negative angles of attack, other complicated shock patterns occur near the leading-edge tip. These heretofore unknown flow patterns show the power and utility of CFD for investigating novel hypersonic configurations such as waveriders.

Development of discrete gas kinetic scheme for simulation of 3D viscous incompressible and compressible flows

NASA Astrophysics Data System (ADS)

Yang, L. M.; Shu, C.; Wang, Y.; Sun, Y.

2016-08-01

The sphere function-based gas kinetic scheme (GKS), which was presented by Shu and his coworkers [23] for simulation of inviscid compressible flows, is extended to simulate 3D viscous incompressible and compressible flows in this work. Firstly, we use certain discrete points to represent the spherical surface in the phase velocity space. Then, integrals along the spherical surface for conservation forms of moments, which are needed to recover 3D Navier-Stokes equations, are approximated by integral quadrature. The basic requirement is that these conservation forms of moments can be exactly satisfied by weighted summation of distribution functions at discrete points. It was found that the integral quadrature by eight discrete points on the spherical surface, which forms the D3Q8 discrete velocity model, can exactly match the integral. In this way, the conservative variables and numerical fluxes can be computed by weighted summation of distribution functions at eight discrete points. That is, the application of complicated formulations resultant from integrals can be replaced by a simple solution process. Several numerical examples including laminar flat plate boundary layer, 3D lid-driven cavity flow, steady flow through a 90° bending square duct, transonic flow around DPW-W1 wing and supersonic flow around NACA0012 airfoil are chosen to validate the proposed scheme. Numerical results demonstrate that the present scheme can provide reasonable numerical results for 3D viscous flows.

Solution of plane cascade flow using improved surface singularity methods

NASA Technical Reports Server (NTRS)

Mcfarland, E. R.

1981-01-01

A solution method has been developed for calculating compressible inviscid flow through a linear cascade of arbitrary blade shapes. The method uses advanced surface singularity formulations which were adapted from those found in current external flow analyses. The resulting solution technique provides a fast flexible calculation for flows through turbomachinery blade rows. The solution method and some examples of the method's capabilities are presented.

Gas turbine engine with radial diffuser and shortened mid section

DOEpatents

Charron, Richard C.; Montgomery, Matthew D.

2015-09-08

An industrial gas turbine engine (10), including: a can annular combustion assembly (80), having a plurality of discrete flow ducts configured to receive combustion gas from respective combustors (82) and deliver the combustion gas along a straight flow path at a speed and orientation appropriate for delivery directly onto the first row (56) of turbine blades (62); and a compressor diffuser (32) having a redirecting surface (130, 140) configured to receive an axial flow of compressed air and redirect the axial flow of compressed air radially outward.

Two-dimensional compressible flow in centrifugal compressors with straight blades

NASA Technical Reports Server (NTRS)

Stanitz, John D; Ellis, Gaylord O

1950-01-01

Six numerical examples are presented for steady, two-dimensional, compressible, nonviscous flow in centrifugal compressors with thin straight blades, the center lines of which generate the surface of a right circular cone when rotated about the axis of the compressor. A seventh example is presented for incompressible flow. The solutions were obtained in a region of the compressors, including the impeller tip, that was considered to be unaffected by the diffuser vanes or by the impeller-inlet configuration. Each solution applies to radial and mixed flow compressors with various cone angles but with the same angle between blades on the conic flow surface. The solution also apply to radial and mixed flow turbines with the rotation and the flow direction reversed. The effects of variations in the following parameters were investigated: (1) flow rate, (2) impeller-tip speed, (3) variation of passage height with radius, and (4) angle between blades on conic flow surface. The numerical results are presented in plots of the streamlines and constant Mach number lines. Correlation equations are developed whereby the flow conditions in any impeller with straight blades can be determined (in the region investigated by this analysis) for all operating conditions.

Work Function of Oxide Ultrathin Films on the Ag(100) Surface.

PubMed

Sementa, Luca; Barcaro, Giovanni; Negreiros, Fabio R; Thomas, Iorwerth O; Netzer, Falko P; Ferrari, Anna Maria; Fortunelli, Alessandro

2012-02-14

Theoretical calculations of the work function of monolayer (ML) and bilayer (BL) oxide films on the Ag(100) surface are reported and analyzed as a function of the nature of the oxide for first-row transition metals. The contributions due to charge compression, charge transfer and rumpling are singled out. It is found that the presence of empty d-orbitals in the oxide metal can entail a charge flow from the Ag(100) surface to the oxide film which counteracts the decrease in the work function due to charge compression. This flow can also depend on the thickness of the film and be reduced in passing from ML to BL systems. A regular trend is observed along first-row transition metals, exhibiting a maximum for CuO, in which the charge flow to the oxide is so strong as to reverse the direction of rumpling. A simple protocol to estimate separately the contribution due to charge compression is discussed, and the difference between the work function of the bare metal surface and a Pauling-like electronegativity of the free oxide slabs is used as a descriptor quantity to predict the direction of charge transfer.

Measurement of corner separation zone response on a compression ramp to plasma actuation within the hypersonic boundary layer

NASA Astrophysics Data System (ADS)

Hedlund, Brock E.; Houpt, Alec W.; Gordeyev, Stanislav V.; Leonov, Sergey B.

2017-10-01

This study was performed to characterize the dominant frequencies present in the boundary layer uptsream of and in the corner separation zone of a compression surface in Mach 4.5 flow and to determine a control effect of transient plasma actuation on the boundary layer. Schlieren imaging was used to distinguish the corner separation zone for 20°, 25°, and 30° compression ramps mounted on flat plates. Spectra of the natural disturbances present in the boundary layer and separation zone were gathered using a high-speed Shack-Hartmann wavefront sensor and surface mounted PCBTM pressure sensors while varying flow parameters by adjusting total pressure, temperature, and ramp angle. Shallow cavity discharge plasma actuators were used as a high-frequency localized thermal forcing mechanism of the boundary layer. The plasma effect was negligible for forcing frequencies (50 kHz) below the natural dominant frequency (~55-80 kHz). High frequency perturbations that can promote the transition to turbulence were amplified when the plasma forcing frequency (100 kHz) was higher than the natural dominant frequency (~55-80 kHz). This technique can potentially be used for active control of hypersonic boundary layer transition and the supersonic flow structure on the compression surface.

Chemorheology of in-mold coating for compression molded SMC applications

NASA Astrophysics Data System (ADS)

Ko, Seunghyun; Straus, Elliott J.; Castro, Jose M.

2015-05-01

In-mold coating (IMC) is applied to compression molded sheet molding compound (SMC) exterior automotive or truck body panels as an environmentally friendly alternative to make the surface conductive for subsequent electrostatic painting operations. The coating is a thermosetting liquid that when injected onto the surface of the part cures and bonds to provide a smooth conductive surface. In order to optimize the IMC process, it is essential to predict the time available for flow, that is the time before the thermosetting reaction starts (inhibition time) as well as the time when the coating has enough structural integrity so that the mold can be opened without damaging the part surface (cure time). To predict both the inhibition time and the cure time, it is critical to study the chemorheology of IMC. In this paper, we study the chemorheology for a typical commercial IMC system, and show its relevance to both the flow and cure time for the IMC stage during SMC compression molding.

A model for the plastic flow of landslides

USGS Publications Warehouse

Savage, William Z.; Smith, William K.

1986-01-01

To further the understanding of the mechanics of landslide flow, we present a model that predicts many of the observed attributes of landslides. The model is based on an integration of the hyperbolic differential equations for stress and velocity fields in a two-dimensional, inclined, semi-infinite half-space of Coulomb plastic material under elevated pore pressure and gravity. Our landslide model predicts commonly observed features. For example, compressive (passive), plug, or extending (active) flow will occur under appropriate longitudinal strain rates. Also, the model predicts that longitudinal stresses increase elliptically with depth to the basal slide plane, and that stress and velocity characteristics, surfaces along which discontinuities in stress and velocity are propagated, are coincident. Finally, the model shows how thrust and normal faults develop at the landslide surface in compressive and extending flow.

Solution of plane cascade flow using improved surface singularity methods. [application of panel method to internal aerodynamics

NASA Technical Reports Server (NTRS)

Mcfarland, E. R.

1981-01-01

A solution method was developed for calculating compressible inviscid flow through a linear cascade of arbitrary blade shapes. The method uses advanced surface singularity formulations which were adapted from those in current external flow analyses. The resulting solution technique provides a fast flexible calculation for flows through turbomachinery blade rows. The solution method and some examples of the method's capabilities are presented.

Cones in Supersonic Flow

NASA Technical Reports Server (NTRS)

Hantzsche, W.; Wendt, H.

1947-01-01

In the case of cones in axially symmetric flow of supersonic velocity, adiabatic compression takes place between shock wave and surface of the cone. Interpolation curves betwen shock polars and the surface are therefore necessary for the complete understanding of this type of flow. They are given in the present report by graphical-numerical integration of the differential equation for all cone angles and airspeeds.

Viscous compressible flow direct and inverse computation and illustrations

NASA Technical Reports Server (NTRS)

Yang, T. T.; Ntone, F.

1986-01-01

An algorithm for laminar and turbulent viscous compressible two dimensional flows is presented. For the application of precise boundary conditions over an arbitrary body surface, a body-fitted coordinate system is used in the physical plane. A thin-layer approximation of tne Navier-Stokes equations is introduced to keep the viscous terms relatively simple. The flow field computation is performed in the transformed plane. A factorized, implicit scheme is used to facilitate the computation. Sample calculations, for Couette flow, developing pipe flow, an isolated airflow, two dimensional compressor cascade flow, and segmental compressor blade design are presented. To a certain extent, the effective use of the direct solver depends on the user's skill in setting up the gridwork, the time step size and the choice of the artificial viscosity. The design feature of the algorithm, an iterative scheme to correct geometry for a specified surface pressure distribution, works well for subsonic flows. A more elaborate correction scheme is required in treating transonic flows where local shock waves may be involved.

Cooled-Spool Piston Compressor

NASA Technical Reports Server (NTRS)

Morris, Brian G.

1994-01-01

Proposed cooled-spool piston compressor driven by hydraulic power and features internal cooling of piston by flowing hydraulic fluid to limit temperature of compressed gas. Provides sufficient cooling for higher compression ratios or reactive gases. Unlike conventional piston compressors, all parts of compressed gas lie at all times within relatively short distance of cooled surface so that gas cooled more effectively.

Method of analysis for compressible flow through mixed-flow centrifugal impellers of arbitrary design

NASA Technical Reports Server (NTRS)

Hamrick, Joseph T; Ginsburg, Ambrose; Osborn, Walter M

1952-01-01

A method is presented for analysis of the compressible flow between the hub and the shroud of mixed-flow impellers of arbitrary design. Axial symmetry was assumed, but the forces in the meridional (hub to shroud) plane, which are derived from tangential pressure gradients, were taken into account. The method was applied to an experimental mixed-flow impeller. The analysis of the flow in the meridional plane of the impeller showed that the rotational forces, the blade curvature, and the hub-shroud profile can introduce severe velocity gradients along the hub and the shroud surfaces. Choked flow at the impeller inlet as determined by the analysis was verified by experimental results.

Numerical calculation of the internal flow field in a centrifugal compressor impeller

NASA Technical Reports Server (NTRS)

Walitt, L.; Harp, J. L., Jr.; Liu, C. Y.

1975-01-01

An iterative numerical method has been developed for the calculation of steady, three-dimensional, viscous, compressible flow fields in centrifugal compressor impellers. The computer code, which embodies the method, solves the steady three dimensional, compressible Navier-Stokes equations in rotating, curvilinear coordinates. The solution takes place on blade-to-blade surfaces of revolution which move from the hub to the shroud during each iteration.

Physics of singularities in pressure-impulse theory

NASA Astrophysics Data System (ADS)

Krechetnikov, R.

2018-05-01

The classical solution in the pressure-impulse theory for the inviscid, incompressible, and zero-surface-tension water impact of a flat plate at zero dead-rise angle exhibits both singular-in-time initial fluid acceleration, ∂v /∂ t |t =0˜δ (t ) , and a near-plate-edge spatial singularity in the velocity distribution, v ˜r-1 /2 , where r is the distance from the plate edge. The latter velocity divergence also leads to the interface being stretched infinitely right after the impact, which is another nonphysical artifact. From the point of view of matched asymptotic analysis, this classical solution is a singular limit when three physical quantities achieve limiting values: sound speed c0→∞ , fluid kinematic viscosity ν →0 , and surface tension σ →0 . This leaves open a question on how to resolve these singularities mathematically by including the neglected physical effects—compressibility, viscosity, and surface tension—first one by one and then culminating in the local compressible viscous solution valid for t →0 and r →0 , demonstrating a nontrivial flow structure that changes with the degree of the bulk compressibility. In the course of this study, by starting with the general physically relevant formulation of compressible viscous flow, we clarify the parameter range(s) of validity of the key analytical solutions including classical ones (inviscid incompressible and compressible, etc.) and understand the solution structure, its intermediate asymptotics nature, characteristics influencing physical processes, and the role of potential and rotational flow components. In particular, it is pointed out that sufficiently close to the plate edge surface tension must be taken into account. Overall, the idea is to highlight the interesting physics behind the singularities in the pressure-impulse theory.

Analysis and calculation by integral methods of laminar compressible boundary-layer with heat transfer and with and without pressure gradient

NASA Technical Reports Server (NTRS)

Morduchow, Morris

1955-01-01

A survey of integral methods in laminar-boundary-layer analysis is first given. A simple and sufficiently accurate method for practical purposes of calculating the properties (including stability) of the laminar compressible boundary layer in an axial pressure gradient with heat transfer at the wall is presented. For flow over a flat plate, the method is applicable for an arbitrarily prescribed distribution of temperature along the surface and for any given constant Prandtl number close to unity. For flow in a pressure gradient, the method is based on a Prandtl number of unity and a uniform wall temperature. A simple and accurate method of determining the separation point in a compressible flow with an adverse pressure gradient over a surface at a given uniform wall temperature is developed. The analysis is based on an extension of the Karman-Pohlhausen method to the momentum and the thermal energy equations in conjunction with fourth- and especially higher degree velocity and stagnation-enthalpy profiles.

Laser driven supersonic flow over a compressible foam surface on the Nike lasera)

NASA Astrophysics Data System (ADS)

Harding, E. C.; Drake, R. P.; Aglitskiy, Y.; Plewa, T.; Velikovich, A. L.; Gillespie, R. S.; Weaver, J. L.; Visco, A.; Grosskopf, M. J.; Ditmar, J. R.

2010-05-01

A laser driven millimeter-scale target was used to generate a supersonic shear layer in an attempt to create a Kelvin-Helmholtz (KH) unstable interface in a high-energy-density (HED) plasma. The KH instability is a fundamental fluid instability that remains unexplored in HED plasmas, which are relevant to the inertial confinement fusion and astrophysical environments. In the experiment presented here the Nike laser [S. P. Obenschain et al., Phys. Plasmas 3, 2098 (1996)] was used to create and drive Al plasma over a rippled foam surface. In response to the supersonic Al flow (Mach=2.6±1.1) shocks should form in the Al flow near the perturbations. The experimental data were used to infer the existence and location of these shocks. In addition, the interface perturbations show growth that has possible contributions from both KH and Richtmyer-Meshkov instabilities. Since compressible shear layers exhibit smaller growth, it is important to use the KH growth rate derived from the compressible dispersion relation.

A Numerical Simulation of a Normal Sonic Jet into a Hypersonic Cross-Flow

NASA Technical Reports Server (NTRS)

Jeffries, Damon K.; Krishnamurthy, Ramesh; Chandra, Suresh

1997-01-01

This study involves numerical modeling of a normal sonic jet injection into a hypersonic cross-flow. The numerical code used for simulation is GASP (General Aerodynamic Simulation Program.) First the numerical predictions are compared with well established solutions for compressible laminar flow. Then comparisons are made with non-injection test case measurements of surface pressure distributions. Good agreement with the measurements is observed. Currently comparisons are underway with the injection case. All the experimental data were generated at the Southampton University Light Piston Isentropic Compression Tube.

Visualization and analysis of flow structures in an open cavity

NASA Astrophysics Data System (ADS)

Liu, Jun; Cai, Jinsheng; Yang, Dangguo; Wu, Junqiang; Wang, Xiansheng

2018-05-01

A numerical study is performed on the supersonic flow over an open cavity at Mach number of 1.5. A newly developed visualization method is employed to visualize the complicated flow structures, which provide an insight into major flow physics. Four types of shock/compressive waves which existed in experimental schlieren are observed in numerical visualization results. Furthermore, other flow structures such as multi-scale vortices are also obtained in the numerical results. And a new type of shocklet which is beneath large vortices is found. The shocklet beneath the vortex originates from leading edge, then, is strengthened by successive interactions between feedback compressive waves and its attached vortex. Finally, it collides against the trailing surface and generates a large number of feedback compressive waves and intensive pressure fluctuations. It is suggested that the shocklets beneath vortex play an important role of cavity self-sustained oscillation.

Progress in turbulence modeling for complex flow fields including effects of compressibility

NASA Technical Reports Server (NTRS)

Wilcox, D. C.; Rubesin, M. W.

1980-01-01

Two second-order-closure turbulence models were devised that are suitable for predicting properties of complex turbulent flow fields in both incompressible and compressible fluids. One model is of the "two-equation" variety in which closure is accomplished by introducing an eddy viscosity which depends on both a turbulent mixing energy and a dissipation rate per unit energy, that is, a specific dissipation rate. The other model is a "Reynolds stress equation" (RSE) formulation in which all components of the Reynolds stress tensor and turbulent heat-flux vector are computed directly and are scaled by the specific dissipation rate. Computations based on these models are compared with measurements for the following flow fields: (a) low speed, high Reynolds number channel flows with plane strain or uniform shear; (b) equilibrium turbulent boundary layers with and without pressure gradients or effects of compressibility; and (c) flow over a convex surface with and without a pressure gradient.

Shock wave-free interface interaction

NASA Astrophysics Data System (ADS)

Frolov, Roman; Minev, Peter; Krechetnikov, Rouslan

2016-11-01

The problem of shock wave-free interface interaction has been widely studied in the context of compressible two-fluid flows using analytical, experimental, and numerical techniques. While various physical effects and possible interaction patterns for various geometries have been identified in the literature, the effects of viscosity and surface tension are usually neglected in such models. In our study, we apply a novel numerical algorithm for simulation of viscous compressible two-fluid flows with surface tension to investigate the influence of these effects on the shock-interface interaction. The method combines together the ideas from Finite Volume adaptation of invariant domains preserving algorithm for systems of hyperbolic conservation laws by Guermond and Popov and ADI parallel solver for viscous incompressible NSEs by Guermond and Minev. This combination has been further extended to a two-fluid flow case, including surface tension effects. Here we report on a quantitative study of how surface tension and viscosity affect the structure of the shock wave-free interface interaction region.

Two-dimensional subsonic compressible flow past elliptic cylinders

NASA Technical Reports Server (NTRS)

Kaplan, Carl

1938-01-01

The method of Poggi is used to calculate, for perfect fluids, the effect of compressibility upon the flow on the surface of an elliptic cylinder at zero angle of attack and with no circulation. The result is expressed in a closed form and represents a rigorous determination of the velocity of the fluid at the surface of the obstacle insofar as the second approximation is concerned. Comparison is made with Hooker's treatment of the same problem according to the method of Janzen and Rayleight and it is found that, for thick elliptic cylinders, the two methods agree very well. The labor of computation is considerably reduced by the present solution.

An interface capturing scheme for modeling atomization in compressible flows

NASA Astrophysics Data System (ADS)

Garrick, Daniel P.; Hagen, Wyatt A.; Regele, Jonathan D.

2017-09-01

The study of atomization in supersonic flow is critical to ensuring reliable ignition of scramjet combustors under startup conditions. Numerical methods incorporating surface tension effects have largely focused on the incompressible regime as most atomization applications occur at low Mach numbers. Simulating surface tension effects in compressible flow requires robust numerical methods that can handle discontinuities caused by both shocks and material interfaces with high density ratios. In this work, a shock and interface capturing scheme is developed that uses the Harten-Lax-van Leer-Contact (HLLC) Riemann solver while a Tangent of Hyperbola for INterface Capturing (THINC) interface reconstruction scheme retains the fluid immiscibility condition in the volume fraction and phasic densities in the context of the five equation model. The approach includes the effects of compressibility, surface tension, and molecular viscosity. One and two-dimensional benchmark problems demonstrate the desirable interface sharpening and conservation properties of the approach. Simulations of secondary atomization of a cylindrical water column after its interaction with a shockwave show good qualitative agreement with experimentally observed behavior. Three-dimensional examples of primary atomization of a liquid jet in a Mach 2 crossflow demonstrate the robustness of the method.

An efficient user-oriented method for calculating compressible flow in an about three-dimensional inlets. [panel method

NASA Technical Reports Server (NTRS)

Hess, J. L.; Mack, D. P.; Stockman, N. O.

1979-01-01

A panel method is used to calculate incompressible flow about arbitrary three-dimensional inlets with or without centerbodies for four fundamental flow conditions: unit onset flows parallel to each of the coordinate axes plus static operation. The computing time is scarcely longer than for a single solution. A linear superposition of these solutions quite rigorously gives incompressible flow about the inlet for any angle of attack, angle of yaw, and mass flow rate. Compressibility is accounted for by applying a well-proven correction to the incompressible flow. Since the computing times for the combination and the compressibility correction are small, flows at a large number of inlet operating conditions are obtained rather cheaply. Geometric input is aided by an automatic generating program. A number of graphical output features are provided to aid the user, including surface streamline tracing and automatic generation of curves of curves of constant pressure, Mach number, and flow inclination at selected inlet cross sections. The inlet method and use of the program are described. Illustrative results are presented.

Computation of three-dimensional compressible boundary layers to fourth-order accuracy on wings and fuselages

NASA Technical Reports Server (NTRS)

Iyer, Venkit

1990-01-01

A solution method, fourth-order accurate in the body-normal direction and second-order accurate in the stream surface directions, to solve the compressible 3-D boundary layer equations is presented. The transformation used, the discretization details, and the solution procedure are described. Ten validation cases of varying complexity are presented and results of calculation given. The results range from subsonic flow to supersonic flow and involve 2-D or 3-D geometries. Applications to laminar flow past wing and fuselage-type bodies are discussed. An interface procedure is used to solve the surface Euler equations with the inviscid flow pressure field as the input to assure accurate boundary conditions at the boundary layer edge. Complete details of the computer program used and information necessary to run each of the test cases are given in the Appendix.

Incompressible material point method for free surface flow

NASA Astrophysics Data System (ADS)

Zhang, Fan; Zhang, Xiong; Sze, Kam Yim; Lian, Yanping; Liu, Yan

2017-02-01

To overcome the shortcomings of the weakly compressible material point method (WCMPM) for modeling the free surface flow problems, an incompressible material point method (iMPM) is proposed based on operator splitting technique which splits the solution of momentum equation into two steps. An intermediate velocity field is first obtained by solving the momentum equations ignoring the pressure gradient term, and then the intermediate velocity field is corrected by the pressure term to obtain a divergence-free velocity field. A level set function which represents the signed distance to free surface is used to track the free surface and apply the pressure boundary conditions. Moreover, an hourglass damping is introduced to suppress the spurious velocity modes which are caused by the discretization of the cell center velocity divergence from the grid vertexes velocities when solving pressure Poisson equations. Numerical examples including dam break, oscillation of a cubic liquid drop and a droplet impact into deep pool show that the proposed incompressible material point method is much more accurate and efficient than the weakly compressible material point method in solving free surface flow problems.

Compressible Boundary Layer Investigation for Ramjet/scramjet Inlets and Nozzles

NASA Astrophysics Data System (ADS)

Goldfeld, M. A.; Starov, A. V.; Semenova, Yu. V.

2005-02-01

The results of experimental investigation of a turbulent boundary layer on compression and expansion surfaces are presented. They include the study of the shock wave and/or expansion fan action upon the boundary layer, boundary layer separation and its relaxation. Complex events of paired interactions and the flow on compression convex-concave surfaces were studied [M. Goldfeld, 1993]. The possibility and conditions of the boundary layer relaminarization behind the expansion fan and its effect on the relaxation length are presented. Different model configurations for wide range conditions were investigated. Comparison of results for different interactions was carried out.

Development of a defect stream function, law of the wall/wake method for compressible turbulent boundary layers. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Wahls, Richard A.

1990-01-01

The method presented is designed to improve the accuracy and computational efficiency of existing numerical methods for the solution of flows with compressible turbulent boundary layers. A compressible defect stream function formulation of the governing equations assuming an arbitrary turbulence model is derived. This formulation is advantageous because it has a constrained zero-order approximation with respect to the wall shear stress and the tangential momentum equation has a first integral. Previous problems with this type of formulation near the wall are eliminated by using empirically based analytic expressions to define the flow near the wall. The van Driest law of the wall for velocity and the modified Crocco temperature-velocity relationship are used. The associated compressible law of the wake is determined and it extends the valid range of the analytical expressions beyond the logarithmic region of the boundary layer. The need for an inner-region eddy viscosity model is completely avoided. The near-wall analytic expressions are patched to numerically computed outer region solutions at a point determined during the computation. A new boundary condition on the normal derivative of the tangential velocity at the surface is presented; this condition replaces the no-slip condition and enables numerical integration to the surface with a relatively coarse grid using only an outer region turbulence model. The method was evaluated for incompressible and compressible equilibrium flows and was implemented into an existing Navier-Stokes code using the assumption of local equilibrium flow with respect to the patching. The method has proven to be accurate and efficient.

Linear stability of compressible Taylor-Couette flow

NASA Technical Reports Server (NTRS)

Kao, Kai-Hsiung; Chow, Chuen-Yen

1992-01-01

A temporal stability analysis of compressible Taylor-Couette flow is presented. The viscous flow studied in this paper is contained between two concentric cylinders of infinite length, which are rotating with different angular velocities and are kept at different surface temperatures. The effects of differential rotation and temperature difference on the stability of Taylor-Couette flow are contrasted for a range of Mach numbers ranging from incompressible to Mach 3.0. The relative motion of the cylinders dramatically affects the characteristics of the Couette flow at the onset of instability. The flow is stabilized or destabilized depending upon the temperature ratio and speeds of the two cylinders. Independent of Mach number and temperature ratio, increasing Reynolds number generally promotes a destabilizing effect, indicating the inviscid nature of the Taylor-Couette flow.

Wave phenomena in a high Reynolds number compressible boundary layer

NASA Technical Reports Server (NTRS)

Bayliss, A.; Maestrello, L.; Parikh, P.; Turkel, E.

1985-01-01

Growth of unstable disturbances in a high Reynolds number compressible boundary layer is numerically simulated. Localized periodic surface heating and cooling as a means of active control of these disturbances is studied. It is shown that compressibility in itself stabilizes the flow but at a lower Mach number, significant nonlinear distortions are produced. Phase cancellation is shown to be an effective mechanism for active boundary layer control.

Effect of an eigenstrain on slow viscous flow of compressible fluid films

DOE Office of Scientific and Technical Information (OSTI.GOV)

Murray, P.E.

We present a general formulation of the mechanics of slow viscous flow of slightly compressible fluid films in the presence of an eigenstrain. An eigenstrain represents a constrained volume change due to temperature, concentration of a dissolved species, or a chemical transformation. A silicon dioxide film grown on a silicon surface is an example of a viscous fluid film that is affected by a constrained volume change. We obtain a general expression for pressure in a fluid film produced by a surface chemical reaction accompanied by a volume change. This result is used to study the effect of an eigenstrainmore » on viscous stress relaxation in fluid films.« less

Compressible, unsteady lifting-surface theory for a helicopter rotor in forward flight

NASA Technical Reports Server (NTRS)

Runyan, H. L.; Tai, H.

1985-01-01

A lifting-surface theory has been developed for a helicopter rotor in forward flight for compressible and incompressible flow. The method utilizes the concept of the linearized acceleration potential and makes use of the doublet lattice procedure. Calculations demonstrating the application of the method are given in terms of the lift distribution on a one-bladed rotor, a two-bladed rotor, and a rotor with swept-forward and swept-back tips. Also, the lift on a rotor vibrating in a pitching mode at 4 per revolution is given. Compressibility effects and interference effects for a two-bladed rotor are discussed.

Compressible liquid flow in nano- or micro-sized circular tubes considering wall-liquid Lifshitz-van der Waals interaction

NASA Astrophysics Data System (ADS)

Zhang, Xueling; Zhu, Weiyao; Cai, Qiang; Shi, Yutao; Wu, Xuehong; Jin, Tingxiang; Yang, Lianzhi; Song, Hongqing

2018-06-01

Although nano- and micro-scale phenomena for fluid flows are ubiquitous in tight oil reservoirs or in nano- or micro-sized channels, the mechanisms behind them remain unclear. In this study, we consider the wall-liquid interaction to investigate the flow mechanisms behind a compressible liquid flow in nano- or micro-sized circular tubes. We assume that the liquid is attracted by the wall surface primarily by the Lifshitz-van der Waals (LW) force, whereas electrostatic forces are negligible. The long-range LW force is thus introduced into the Navier-Stokes equations. The nonlinear equations of motion are decoupled by using the hydrodynamic vorticity-stream functions, from which an approximate analytical perturbation solution is obtained. The proposed model considers the LW force and liquid compressibility to obtain the velocity and pressure fields, which are consistent with experimentally observed micro-size effects. A smaller tube radius implies smaller dimensionless velocity, and when the tube radius decreases to a certain radius Rm, a fluid no longer flows, where Rm is the lower limit of the movable-fluid radius. The radius Rm is calculated, and the results are consistent with previous experimental results. These results reveal that micro-size effects are caused by liquid compressibility and wall-liquid interactions, such as the LW force, for a liquid flowing in nano- or micro-sized channels or pores. The attractive LW force enhances the flow's radial resistance, and the liquid compressibility transmits the radial resistance to the streaming direction via volume deformation, thereby decreasing the streaming velocity.

Computer program for calculating laminar, transitional, and turbulent boundary layers for a compressible axisymmetric flow

NASA Technical Reports Server (NTRS)

Albers, J. A.; Gregg, J. L.

1974-01-01

A finite-difference program is described for calculating the viscous compressible boundary layer flow over either planar or axisymmetric surfaces. The flow may be initially laminar and progress through a transitional zone to fully turbulent flow, or it may remain laminar, depending on the imposed boundary conditions, laws of viscosity, and numerical solution of the momentum and energy equations. The flow may also be forced into a turbulent flow at a chosen spot by the data input. The input may contain the factors of arbitrary Reynolds number, free-stream Mach number, free-stream turbulence, wall heating or cooling, longitudinal wall curvature, wall suction or blowing, and wall roughness. The solution may start from an initial Falkner-Skan similarity profile, an approximate equilibrium turbulent profile, or an initial arbitrary input profile.

A compressible solution of the Navier-Stokes equations for turbulent flow about an airfoil

NASA Technical Reports Server (NTRS)

Shamroth, S. J.; Gibeling, H. J.

1979-01-01

A compressible time dependent solution of the Navier-Stokes equations including a transition turbulence model is obtained for the isolated airfoil flow field problem. The equations are solved by a consistently split linearized block implicit scheme. A nonorthogonal body-fitted coordinate system is used which has maximum resolution near the airfoil surface and in the region of the airfoil leading edge. The transition turbulence model is based upon the turbulence kinetic energy equation and predicts regions of laminar, transitional, and turbulent flow. Mean flow field and turbulence field results are presented for an NACA 0012 airfoil at zero and nonzero incidence angles of Reynolds number up to one million and low subsonic Mach numbers.

Numerical modeling of carrier gas flow in atomic layer deposition vacuum reactor: A comparative study of lattice Boltzmann models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pan, Dongqing; Chien Jen, Tien; Li, Tao

2014-01-15

This paper characterizes the carrier gas flow in the atomic layer deposition (ALD) vacuum reactor by introducing Lattice Boltzmann Method (LBM) to the ALD simulation through a comparative study of two LBM models. Numerical models of gas flow are constructed and implemented in two-dimensional geometry based on lattice Bhatnagar–Gross–Krook (LBGK)-D2Q9 model and two-relaxation-time (TRT) model. Both incompressible and compressible scenarios are simulated and the two models are compared in the aspects of flow features, stability, and efficiency. Our simulation outcome reveals that, for our specific ALD vacuum reactor, TRT model generates better steady laminar flow features all over the domainmore » with better stability and reliability than LBGK-D2Q9 model especially when considering the compressible effects of the gas flow. The LBM-TRT is verified indirectly by comparing the numerical result with conventional continuum-based computational fluid dynamics solvers, and it shows very good agreement with these conventional methods. The velocity field of carrier gas flow through ALD vacuum reactor was characterized by LBM-TRT model finally. The flow in ALD is in a laminar steady state with velocity concentrated at the corners and around the wafer. The effects of flow fields on precursor distributions, surface absorptions, and surface reactions are discussed in detail. Steady and evenly distributed velocity field contribute to higher precursor concentration near the wafer and relatively lower particle velocities help to achieve better surface adsorption and deposition. The ALD reactor geometry needs to be considered carefully if a steady and laminar flow field around the wafer and better surface deposition are desired.« less

A study of the compressible flow through a diffusing S-duct

NASA Technical Reports Server (NTRS)

Wellborn, Steven R.; Okiishi, Theodore H.; Reichert, Bruce A.

1993-01-01

Benchmark aerodynamic data are presented for compressible flow through a representative S-duct configuration. A numerical prediction of the S-duct flow field, obtained from a subsonic parabolized Navier-Stokes algorithm, is also shown. The experimental and numerical results are compared. Measurements of the three-dimensional velocity field, total pressures, and static pressures were obtained at five cross-sectional planes. Aerodynamic data were gathered with calibrated pneumatic probes. Surface static pressure and surface flow visualization data were also acquired. All reported tests were conducted with an inlet centerline Mach number of 0.6. The Reynolds number, based on the inlet centerline velocity and duct inlet diameter, was 2.6 x 10(exp 6). Thin inlet turbulent boundary layers existed. The collected data should be beneficial to aircraft inlet designers and the measurements are suitable for the validation of computational codes. The results show that a region of streamwise flow separation occurred within the duct. Details about the separated flow region, including mechanisms which drive this complicated flow phenomenon, are discussed. Results also indicate that the duct curvature induces strong pressure driven secondary flows. The cross flows evolve into counter-rotating vortices. These vortices convect low momentum fluid of the boundary layer toward the center of the duct, degrading both the uniformity and magnitude of the total pressure profile.

Experimental investigation of a Mach 6 fixed-geometry inlet featuring a swept external-internal compression flow field

NASA Technical Reports Server (NTRS)

Torrence, M. G.

1975-01-01

An investigation of a fixed-geometry, swept external-internal compression inlet was conducted at a Mach number of 6.0 and a test-section Reynolds number of 1.55 x 10 to the 7th power per meter. The test conditions was constant for all runs with stagnation pressure and temperature at 20 atmospheres and 500 K, respectively. Tests were made at angles of attack of -5 deg, 0 deg, 3 deg, and 5 deg. Measurements consisted of pitot- and static-pressure surveys in inlet throat, wall static pressures, and surface temperatures. Boundary-layer bleed was provided on the centerbody and on the cowl internal surface. The inlet performance was consistently high over the range of the angle of attack tested, with an overall average total pressure recovery of 78 percent and corresponding adiabatic kinetic-energy efficiency of 99 percent. The inlet throat flow distribution was uniform and the Mach number and pressure level were of the correct magnitude for efficient combustor design. The utilization of a swept compression field to meet the starting requirements of a fixed-geometry inlet produced neither flow instability nor a tendency to unstart.

Unsteady turbulent boundary layers in swimming rainbow trout.

PubMed

Yanase, Kazutaka; Saarenrinne, Pentti

2015-05-01

The boundary layers of rainbow trout, Oncorhynchus mykiss, swimming at 1.02±0.09 L s(-1) (mean±s.d., N=4), were measured by the particle image velocimetry (PIV) technique at a Reynolds number of 4×10(5). The boundary layer profile showed unsteadiness, oscillating above and beneath the classical logarithmic law of the wall with body motion. Across the entire surface regions that were measured, local Reynolds numbers based on momentum thickness, which is the distance that is perpendicular to the fish surface through which the boundary layer momentum flows at free-stream velocity, were greater than the critical value of 320 for the laminar-to-turbulent transition. The skin friction was dampened on the convex surface while the surface was moving towards a free-stream flow and increased on the concave surface while retreating. These observations contradict the result of a previous study using different species swimming by different methods. Boundary layer compression accompanied by an increase in local skin friction was not observed. Thus, the overall results may not support absolutely the Bone-Lighthill boundary layer thinning hypothesis that the undulatory motions of swimming fish cause a large increase in their friction drag because of the compression of the boundary layer. In some cases, marginal flow separation occurred on the convex surface in the relatively anterior surface region, but the separated flow reattached to the fish surface immediately downstream. Therefore, we believe that a severe impact due to induced drag components (i.e. pressure drag) on the swimming performance, an inevitable consequence of flow separation, was avoided. © 2015. Published by The Company of Biologists Ltd.

A finite difference model for free surface gravity drainage

DOE Office of Scientific and Technical Information (OSTI.GOV)

Couri, F.R.; Ramey, H.J. Jr.

1993-09-01

The unconfined gravity flow of liquid with a free surface into a well is a classical well test problem which has not been well understood by either hydrologists or petroleum engineers. Paradigms have led many authors to treat an incompressible flow as compressible flow to justify the delayed yield behavior of a time-drawdown test. A finite-difference model has been developed to simulate the free surface gravity flow of an unconfined single phase, infinitely large reservoir into a well. The model was verified with experimental results in sandbox models in the literature and with classical methods applied to observation wells inmore » the Groundwater literature. The simulator response was also compared with analytical Theis (1935) and Ramey et al. (1989) approaches for wellbore pressure at late producing times. The seepage face in the sandface and the delayed yield behavior were reproduced by the model considering a small liquid compressibility and incompressible porous medium. The potential buildup (recovery) simulated by the model evidenced a different- phenomenon from the drawdown, contrary to statements found in the Groundwater literature. Graphs of buildup potential vs time, buildup seepage face length vs time, and free surface head and sand bottom head radial profiles evidenced that the liquid refills the desaturating cone as a flat moving surface. The late time pseudo radial behavior was only approached after exaggerated long times.« less

Heat Transfer to Surfaces of Finite Catalytic Activity in Frozen Dissociated Hypersonic Flow

NASA Technical Reports Server (NTRS)

Chung, Paul M.; Anderson, Aemer D.

1961-01-01

The heat transfer due to catalytic recombination of a partially dissociated diatomic gas along the surfaces of two-dimensional and axisymmetric bodies with finite catalytic efficiencies is studied analytically. An integral method is employed resulting in simple yet relatively complete solutions for the particular configurations considered. A closed form solution is derived which enables one to calculate atom mass-fraction distribution, therefore catalytic heat transfer distribution, along the surface of a flat plate in frozen compressible flow with and without transpiration. Numerical calculations are made to determine the atom mass-fraction distribution along an axisymmetric conical body with spherical nose in frozen hypersonic compressible flow. A simple solution based on a local similarity concept is found to be in good agreement with these numerical calculations. The conditions are given for which the local similarity solution is expected to be satisfactory. The limitations on the practical application of the analysis to the flight of the blunt bodies in the atmosphere are discussed. The use of boundary-layer theory and the assumption of frozen flow restrict application of the analysis to altitudes between about 150,000 and 250,000 feet.

An improved weakly compressible SPH method for simulating free surface flows of viscous and viscoelastic fluids

NASA Astrophysics Data System (ADS)

Xu, Xiaoyang; Deng, Xiao-Long

2016-04-01

In this paper, an improved weakly compressible smoothed particle hydrodynamics (SPH) method is proposed to simulate transient free surface flows of viscous and viscoelastic fluids. The improved SPH algorithm includes the implementation of (i) the mixed symmetric correction of kernel gradient to improve the accuracy and stability of traditional SPH method and (ii) the Rusanov flux in the continuity equation for improving the computation of pressure distributions in the dynamics of liquids. To assess the effectiveness of the improved SPH algorithm, a number of numerical examples including the stretching of an initially circular water drop, dam breaking flow against a vertical wall, the impact of viscous and viscoelastic fluid drop with a rigid wall, and the extrudate swell of viscoelastic fluid have been presented and compared with available numerical and experimental data in literature. The convergent behavior of the improved SPH algorithm has also been studied by using different number of particles. All numerical results demonstrate that the improved SPH algorithm proposed here is capable of modeling free surface flows of viscous and viscoelastic fluids accurately and stably, and even more important, also computing an accurate and little oscillatory pressure field.

Goertler instability in compressible boundary layers along curved surfaces with suction and cooling

NASA Technical Reports Server (NTRS)

El-Hady, N.; Verma, A. K.

1982-01-01

The Goertler instability of the laminar compressible boundary layer flows along concave surfaces is investigated. The linearized disturbance equations for the three-dimensional, counter-rotating streamwise vortices in two-dimensional boundary layers are presented in an orthogonal curvilinear coordinate. The basic approximation of the disturbance equations, that includes the effect of the growth of the boundary layer, is considered and solved numerically. The effect of compressibility on critical stability limits, growth rates, and amplitude ratios of the vortices is evaluated for a range of Mach numbers for 0 to 5. The effect of wall cooling and suction of the boundary layer on the development of Goertler vortices is investigated for different Mach numbers.

Anelastic and Compressible Simulation of Moist Dynamics at Planetary Scales

NASA Astrophysics Data System (ADS)

Kurowski, M.; Smolarkiewicz, P. K.; Grabowski, W.

2015-12-01

Moist anelastic and compressible numerical solutions to the planetary baroclinic instability and climate benchmarks are compared. The solutions are obtained applying a consistent numerical framework for dis- crete integrations of the various nonhydrostatic flow equations. Moist extension of the baroclinic instability benchmark is formulated as an analog of the dry case. Flow patterns, surface vertical vorticity and pressure, total kinetic energy, power spectra, and total amount of condensed water are analyzed. The climate bench- mark extends the baroclinic instability study by addressing long-term statistics of an idealized planetary equilibrium and associated meridional transports. Short-term deterministic anelastic and compressible so- lutions differ significantly. In particular, anelastic baroclinic eddies propagate faster and develop slower owing to, respectively, modified dispersion relation and abbreviated baroclinic vorticity production. These eddies also carry less kinetic energy, and the onset of their rapid growth occurs later than for the compressible solutions. The observed differences between the two solutions are sensitive to initial conditions as they di- minish for large-amplitude excitations of the instability. In particular, on the climatic time scales, the anelastic and compressible solutions evince similar zonally averaged flow patterns with the matching meridional transports of entropy, momentum, and moisture.

A novel low profile wireless flow sensor to monitor hemodynamic changes in cerebral aneurysm

NASA Astrophysics Data System (ADS)

Chen, Yanfei; Jankowitz, Brian T.; Cho, Sung Kwon; Chun, Youngjae

2015-03-01

A proof of concept of low-profile flow sensor has been designed, fabricated, and subsequently tested to demonstrate its feasibility for monitoring hemodynamic changes in cerebral aneurysm. The prototype sensor contains three layers, i.e., a thin polyurethane layer was sandwiched between two sputter-deposited thin film nitinol layers (6μm thick). A novel superhydrophilic surface treatment was used to create hemocompatible surface of thin nitinol electrode layers. A finite element model was conducted using ANSYS Workbench 15.0 Static Structural to optimize the dimensions of flow sensor. A computational fluid dynamics calculations were performed using ANSYS Workbench Fluent to assess the flow velocity patterns within the aneurysm sac. We built a test platform with a z-axis translation stage and an S-beam load cell to compare the capacitance changes of the sensors with different parameters during deformation. Both LCR meter and oscilloscope were used to measure the capacitance and the resonant frequency shifts, respectively. The experimental compression tests demonstrated the linear relationship between the capacitance and applied compression force and decreasing the length, width and increasing the thickness improved the sensor sensitivity. The experimentally measured resonant frequency dropped from 12.7MHz to 12.48MHz, indicating a 0.22MHz shift with 200g ( 2N) compression force while the theoretical resonant frequency shifted 0.35MHz with 50g ( 0.5N). Our recent results demonstrated a feasibility of the low-profile flow sensor for monitoring haemodynamics in cerebral aneurysm region, as well as the efficacy of the use of the surface treated thin film nitinol for the low-profile sensor materials.

A composite velocity procedure for the compressible Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Khosla, P. K.; Rubin, S. G.

1982-01-01

A new boundary-layer relaxation procedure is presented. In the spirit of the theory of matched asymptotic expansions, a multiplicative composite of the appropriate velocity representations for the inviscid and viscous regions is prescribed. The resulting equations are structured so that far from the surface of the body the momentum equations lead to the Bernoulli relation for the pressure, while the continuity equation reduces to the familiar compressible potential equation. Close to the body surface, the governing equations and solution techniques are characteristic of those describing interacting boundary-layers; although, the full Navier-Stokes equations are considered here. Laminar flow calculations for the subsonic flow over an axisymmetric boattail simulator geometry are presented for a variety of Reynolds and Mach numbers. A strongly implicit solution method is applied for the coupled velocity components.

Contour forming of metals by laser peening

DOEpatents

Hackel, Lloyd; Harris, Fritz

2002-01-01

A method and apparatus are provided for forming shapes and contours in metal sections by generating laser induced compressive stress on the surface of the metal workpiece. The laser process can generate deep compressive stresses to shape even thick components without inducing unwanted tensile stress at the metal surface. The precision of the laser-induced stress enables exact prediction and subsequent contouring of parts. A light beam of 10 to 100 J/pulse is imaged to create an energy fluence of 60 to 200 J/cm.sup.2 on an absorptive layer applied over a metal surface. A tamping layer of water is flowed over the absorptive layer. The absorption of laser light causes a plasma to form and consequently creates a shock wave that induces a deep residual compressive stress into the metal. The metal responds to this residual stress by bending.

COMPRESSION SEAL AND SEALING MATERIAL THEREFOR

DOEpatents

Branin, T.G.

1962-05-29

This patent relates to compression seal and more particularly to a seaiing material therefor. The sealing surface is a coating consisting of alternate layers of gold and of a non-gold metal having similar plastic flow properties under pressure as gold. The coating is substantially free from oxidation effects when exposed to ambient atmosphere and does not become brittle when worked, as in a valve. (AEC)

Compression Ratio and Catalyst Aging Effects on Aqueous Ethanol Ignition (Year 2) : Part 2 Catalyst Aging and Effects of Water on Ignition

DOT National Transportation Integrated Search

2009-09-01

A tubular plug-flow reactor under low Reynolds Numbers Re flow regimes, along with a 127 um diameter coiled platinum (Pt) wire, were used to study catalytic surface reactions of nonflammable, fuel-lean mixtures of propane, oxygen, and water vapor dil...

Adaptive Methods for Compressible Flow

DTIC Science & Technology

1994-03-01

labor -intensive task of purpose of this work is to demonstrate the generating acceptable surface triangulations, advantages of integrating the CAD/CAM...sintilar results). L 1 (’-1)(2sn~p) boundary error (MUSCL) The flow variables wre then given by .04 .78% M=asOIne/i .02 AM% v= acosO /sintt .01 .0 p

Compressible Navier-Stokes Equations in a Polyhedral Cylinder with Inflow Boundary Condition

NASA Astrophysics Data System (ADS)

Kwon, Ohsung; Kweon, Jae Ryong

2018-06-01

In this paper our concern is with singularity and regularity of the compressible flows through a non-convex edge in R^3. The flows are governed by the compressible Navies-Stokes equations on the infinite cylinder that has the non-convex edge on the inflow boundary. We split the edge singularity by the Poisson problem from the velocity vector and show that the remainder is twice differentiable while the edge singularity is observed to be propagated into the interior of the cylinder by the transport character of the continuity equation. An interior surface layer starting at the edge is generated and not Lipshitz continuous due to the singularity. The density function shows a very steep change near the interface and its normal derivative has a jump discontinuity across there.

Hypersonic code efficiency and validation studies

NASA Technical Reports Server (NTRS)

Bennett, Bradford C.

1992-01-01

Renewed interest in hypersonic and supersonic flows spurred the development of the Compressible Navier-Stokes (CNS) code. Originally developed for external flows, CNS was modified to enable it to also be applied to internal high speed flows. In the initial phase of this study CNS was applied to both internal flow applications and fellow researchers were taught to run CNS. The second phase of this research was the development of surface grids over various aircraft configurations for the High Speed Research Program (HSRP). The complex nature of these configurations required the development of improved surface grid generation techniques. A significant portion of the grid generation effort was devoted to testing and recommending modifications to early versions of the S3D surface grid generation code.

30 CFR 7.506 - Breathable air components.

Code of Federal Regulations, 2010 CFR

2010-07-01

... cylinders, or boreholes with fans installed on the surface or compressors installed on the surface. Only....5 percent. (c) Breathable air supplied by compressed air from cylinders, fans, or compressors shall provide a minimum flow rate of 12.5 cubic feet per minute of breathable air for each person. (1) Fans or...

CFD Simulations of the IHF Arc-Jet Flow: Compression-Pad/Separation Bolt Wedge Tests

NASA Technical Reports Server (NTRS)

Gokcen, Tahir; Skokova, Kristina A.

2017-01-01

This paper reports computational analyses in support of two wedge tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using two different wedge models, each placed in a free jet downstream of a corresponding different conical nozzle in the Ames 60-MW Interaction Heating Facility. Panel test articles included a metallic separation bolt imbedded in the compression-pad and heat shield materials, resulting in a circular protuberance over a flat plate. As part of the test calibration runs, surface pressure and heat flux measurements on water-cooled calibration plates integrated with the wedge models were also obtained. Surface heating distributions on the test articles as well as arc-jet test environment parameters for each test configuration are obtained through computational fluid dynamics simulations, consistent with the facility and calibration measurements. The present analysis comprises simulations of the non-equilibrium flow field in the facility nozzle, test box, and flow field over test articles, and comparisons with the measured calibration data.

DEVELOPMENT OF PLASTICITY MODEL USING NON ASSOCIATED FLOW RULE FOR HCP MATERIALS INCLUDING ZIRCONIUM FOR NUCLEAR APPLICATIONS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Michael V. Glazoff; Jeong-Whan Yoon

2013-08-01

In this report (prepared in collaboration with Prof. Jeong Whan Yoon, Deakin University, Melbourne, Australia) a research effort was made to develop a non associated flow rule for zirconium. Since Zr is a hexagonally close packed (hcp) material, it is impossible to describe its plastic response under arbitrary loading conditions with any associated flow rule (e.g. von Mises). As a result of strong tension compression asymmetry of the yield stress and anisotropy, zirconium displays plastic behavior that requires a more sophisticated approach. Consequently, a new general asymmetric yield function has been developed which accommodates mathematically the four directional anisotropies alongmore » 0 degrees, 45 degrees, 90 degrees, and biaxial, under tension and compression. Stress anisotropy has been completely decoupled from the r value by using non associated flow plasticity, where yield function and plastic potential have been treated separately to take care of stress and r value directionalities, respectively. This theoretical development has been verified using Zr alloys at room temperature as an example as these materials have very strong SD (Strength Differential) effect. The proposed yield function reasonably well models the evolution of yield surfaces for a zirconium clock rolled plate during in plane and through thickness compression. It has been found that this function can predict both tension and compression asymmetry mathematically without any numerical tolerance and shows the significant improvement compared to any reported functions. Finally, in the end of the report, a program of further research is outlined aimed at constructing tensorial relationships for the temperature and fluence dependent creep surfaces for Zr, Zircaloy 2, and Zircaloy 4.« less

Diffuse interface method for a compressible binary fluid.

PubMed

Liu, Jiewei; Amberg, Gustav; Do-Quang, Minh

2016-01-01

Multicomponent, multiphase, compressible flows are very important in real life, as well as in scientific research, while their modeling is in an early stage. In this paper, we propose a diffuse interface model for compressible binary mixtures, based on the balance of mass, momentum, energy, and the second law of thermodynamics. We show both analytically and numerically that this model is able to describe the phase equilibrium for a real binary mixture (CO_{2} + ethanol is considered in this paper) very well by adjusting the parameter which measures the attraction force between molecules of the two components in the model. We also show that the calculated surface tension of the CO_{2} + ethanol mixture at different concentrations match measurements in the literature when the mixing capillary coefficient is taken to be the geometric mean of the capillary coefficient of each component. Three different cases of two droplets in a shear flow, with the same or different concentration, are simulated, showing that the higher concentration of CO_{2} the smaller the surface tension and the easier the drop deforms.

Studies on the compressibility of wax matrix granules of acetaminophen and their admixtures with various tableting bases.

PubMed

Uhumwangho, M U; Okor, R S

2006-04-01

Matrix granules of acetaminophen have been formed by a melt granulation process whereby the acetaminophen powder was triturated with the melted wax--goat wax, glyceryl monostearate or carnuba wax. The compressibility of the matrix granules and their admixture, with diluent granules (lactose, alpha-cellulose or microcrystalline cellulose) was investigated. The granules were compressed to tablets at a constant load (30 arbitrary units on the load scale) of a manesty single punch machine. Resulting tablets were evaluated for tensile strength (T) and disintegration times (DT). Granule flow was determined by measuring their angle of repose when allowed to fall freely on a level surface. Matrix granules prepared by melt granulation with goat wax or glyceryl monostearate were too sticky and therefore did not flow at all. They were also poorly compressible (T values = 0.20MN/m2). Inclusion of the diluent remarkably improved granule flow property and compressibility. The T values of the tablets (measure of compressibility) increased from about 0.24 to 0.65 MN/m2 during increase in diluent (lactose) content from 20 to 80 %w/w. Microcrystalline cellulose and alpha-cellulose were more effective than lactose in promoting compressibility of the granules. By contrast the matrix granules formed with carnuba wax were free flowing (angle of repose, 18.60). Addition of the diluent further improved flowability slightly. The matrix granules (without a diluent) were readily compressible (T value, 1.79MN/m2). Addition of the diluent (80%w/w) reduced T values (MN/m2) slightly to 1.32 (lactose), 1.48 (alpha-cellulose) and 1.74 (microcrystalline cellulose). Tablets of the matrix granules only, disintegrated rapidly within 3 minutes. DT was further reduced to <30 s by addition of any of the diluents. The indication is that the inclusion of the diluents studied can be used to improve the compressibility of the otherwise poorly compressible matrix granules. Based on the flowability, compressibility, and disintegration data, carnuba wax proved most promising in the melt granulation of the test drug for sustained release applications.

User's Manual for FOMOCO Utilities-Force and Moment Computation Tools for Overset Grids

NASA Technical Reports Server (NTRS)

Chan, William M.; Buning, Pieter G.

1996-01-01

In the numerical computations of flows around complex configurations, accurate calculations of force and moment coefficients for aerodynamic surfaces are required. When overset grid methods are used, the surfaces on which force and moment coefficients are sought typically consist of a collection of overlapping surface grids. Direct integration of flow quantities on the overlapping grids would result in the overlapped regions being counted more than once. The FOMOCO Utilities is a software package for computing flow coefficients (force, moment, and mass flow rate) on a collection of overset surfaces with accurate accounting of the overlapped zones. FOMOCO Utilities can be used in stand-alone mode or in conjunction with the Chimera overset grid compressible Navier-Stokes flow solver OVERFLOW. The software package consists of two modules corresponding to a two-step procedure: (1) hybrid surface grid generation (MIXSUR module), and (2) flow quantities integration (OVERINT module). Instructions on how to use this software package are described in this user's manual. Equations used in the flow coefficients calculation are given in Appendix A.

Further Investigations of Hypersonic Engine Seals

NASA Technical Reports Server (NTRS)

Dunlap, Patrick H., Jr.; Steinetz, Bruce M.; DeMange, Jeffrey J.

2004-01-01

Durable, flexible sliding seals are required in advanced hypersonic engines to seal the perimeters of movable engine ramps for efficient, safe operation in high heat flux environments at temperatures of 2000 to 2500 F. Current seal designs do not meet the demanding requirements for future engines, so NASA's Glenn Research Center is developing advanced seals and preloading devices to overcome these shortfalls. An advanced ceramic wafer seal design and two silicon nitride compression spring designs were evaluated in a series of compression, scrub, and flow tests. Silicon nitride wafer seals survived 2000 in. (50.8 m) of scrubbing at 2000 F against a silicon carbide rub surface with no chips or signs of damage. Flow rates measured for the wafers before and after scrubbing were almost identical and were up to 32 times lower than those recorded for the best braided rope seal flow blockers. Silicon nitride compression springs showed promise conceptually as potential seal preload devices to help maintain seal resiliency.

Sweep and Compressibility Effects on Active Separation Control at High Reynolds Numbers

NASA Technical Reports Server (NTRS)

Seifert, Avi; Pack, LaTunia G.

2000-01-01

This paper explores the effects of compressibility, sweep and excitation location on active separation control at high Reynolds numbers. The model, which was tested in a cryogenic pressurized wind tunnel, simulates the upper surface of a 20% thick GlauertGoldschmied type airfoil at zero angle of attack. The flow is fully turbulent since the tunnel sidewall boundary layer flows over the model. Without control, the flow separates at the highly convex area and a large turbulent separation bubble is formed. Periodic excitation is applied to gradually eliminate the separation bubble. Two alternative blowing slot locations as well as the effect of compressibility, sweep and steady suction or blowing were studied. During the test the Reynolds numbers ranged from 2 to 40 million and Mach numbers ranged from 0.2 to 0.7. Sweep angles were 0 and 30 deg. It was found that excitation must be introduced slightly upstream of the separation region regardless of the sweep angle at low Mach number. Introduction of excitation upstream of the shock wave is more effective than at its foot. Compressibility reduces the ability of steady mass transfer and periodic excitation to control the separation bubble but excitation has an effect on the integral parameters, which is similar to that observed in low Mach numbers. The conventional swept flow scaling is valid for fully and even partially attached flow, but different scaling is required for the separated 3D flow. The effectiveness of the active control is not reduced by sweep. Detailed flow field dynamics are described in the accompanying paper.

Sweep and Compressibility Effects on Active Separation Control at High Reynolds Numbers

NASA Technical Reports Server (NTRS)

Seifert, Avi; Pack, LaTunia G.

2000-01-01

This paper explores the effects of compressibility, sweep and excitation location on active separation control at high Reynolds numbers. The model, which was tested in a cryogenic pressurized wind tunnel, simulates the upper surface of a 20% thick Glauert Goldschmied type airfoil at zero angle of attack. The flow is fully turbulent since the tunnel sidewall boundary layer flows over the model. Without control, the flow separates at the highly convex area and a large turbulent separation bubble is formed. Periodic excitation is applied to gradually eliminate the separation bubble. Two alternative blowing slot locations as well as the effect of compressibility, sweep and steady suction or blowing were studied. During the test the Reynolds numbers ranged from 2 to 40 million and Mach numbers ranged from 0.2 to 0.7. Sweep angles were 0 and 30 deg. It was found that excitation must be introduced slightly upstream of the separation region regardless of the sweep angle at low Mach number. Introduction of excitation upstream of the shock wave is more effective than at its foot. Compressibility reduces the ability of steady mass transfer and periodic excitation to control the separation bubble but excitation has an effect on the integral parameters, which is similar to that observed in low Mach numbers. The conventional swept flow scaling is valid for fully and even partially attached flow, but different scaling is required for the separated 3D flow. The effectiveness of the active control is not reduced by sweep. Detailed flow field dynamics are described in the accompanying paper.

Numerical and Analytical Study of Nonlinear Effects of Transonic Flow Past a Wing Airfoil in Oscillation of its Surface Element

NASA Astrophysics Data System (ADS)

Aul'chenko, S. M.; Zamuraev, V. P.; Kalinina, A. P.

2014-05-01

The present work is devoted to a criterial analysis and mathematical modeling of the influence of forced oscillations of surface elements of a wing airfoil on the shock-wave structure of transonic flow past it. Parameters that govern the regimes of interaction of the oscillatory motion of airfoil sections with the breakdown compression shock have been established. The qualitative and quantitative influence of these parameters on the wave resistance of the airfoil has been investigated.

Methodology for the Design of Streamline-Traced External-Compression Supersonic Inlets

NASA Technical Reports Server (NTRS)

Slater, John W.

2014-01-01

A design methodology based on streamline-tracing is discussed for the design of external-compression, supersonic inlets for flight below Mach 2.0. The methodology establishes a supersonic compression surface and capture cross-section by tracing streamlines through an axisymmetric Busemann flowfield. The compression system of shock and Mach waves is altered through modifications to the leading edge and shoulder of the compression surface. An external terminal shock is established to create subsonic flow which is diffused in the subsonic diffuser. The design methodology was implemented into the SUPIN inlet design tool. SUPIN uses specified design factors to design the inlets and computes the inlet performance, which includes the flow rates, total pressure recovery, and wave drag. A design study was conducted using SUPIN and the Wind-US computational fluid dynamics code to design and analyze the properties of two streamline-traced, external-compression (STEX) supersonic inlets for Mach 1.6 freestream conditions. The STEX inlets were compared to axisymmetric pitot, two-dimensional, and axisymmetric spike inlets. The STEX inlets had slightly lower total pressure recovery and higher levels of total pressure distortion than the axisymmetric spike inlet. The cowl wave drag coefficients of the STEX inlets were 20% of those for the axisymmetric spike inlet. The STEX inlets had external sound pressures that were 37% of those of the axisymmetric spike inlet, which may result in lower adverse sonic boom characteristics. The flexibility of the shape of the capture cross-section may result in benefits for the integration of STEX inlets with aircraft.

Integrated LTCC pressure/flow/temperature multisensor for compressed air diagnostics.

PubMed

Fournier, Yannick; Maeder, Thomas; Boutinard-Rouelle, Grégoire; Barras, Aurélie; Craquelin, Nicolas; Ryser, Peter

2010-01-01

We present a multisensor designed for industrial compressed air diagnostics and combining the measurement of pressure, flow, and temperature, integrated with the corresponding signal conditioning electronics in a single low-temperature co-fired ceramic (LTCC) package. The developed sensor may be soldered onto an integrated electro-fluidic platform by using standard surface mount device (SMD) technology, e.g., as a standard electronic component would be on a printed circuit board, obviating the need for both wires and tubes and thus paving the road towards low-cost integrated electro-fluidic systems. Several performance aspects of this device are presented and discussed, together with electronics design issues.

Integrated LTCC Pressure/Flow/Temperature Multisensor for Compressed Air Diagnostics†

PubMed Central

Fournier, Yannick; Maeder, Thomas; Boutinard-Rouelle, Grégoire; Barras, Aurélie; Craquelin, Nicolas; Ryser, Peter

2010-01-01

We present a multisensor designed for industrial compressed air diagnostics and combining the measurement of pressure, flow, and temperature, integrated with the corresponding signal conditioning electronics in a single low-temperature co-fired ceramic (LTCC) package. The developed sensor may be soldered onto an integrated electro-fluidic platform by using standard surface mount device (SMD) technology, e.g., as a standard electronic component would be on a printed circuit board, obviating the need for both wires and tubes and thus paving the road towards low-cost integrated electro-fluidic systems. Several performance aspects of this device are presented and discussed, together with electronics design issues. PMID:22163518

User's manual for three dimensional boundary layer (BL3-D) code

NASA Technical Reports Server (NTRS)

Anderson, O. L.; Caplin, B.

1985-01-01

An assessment has been made of the applicability of a 3-D boundary layer analysis to the calculation of heat transfer, total pressure losses, and streamline flow patterns on the surface of both stationary and rotating turbine passages. In support of this effort, an analysis has been developed to calculate a general nonorthogonal surface coordinate system for arbitrary 3-D surfaces and also to calculate the boundary layer edge conditions for compressible flow using the surface Euler equations and experimental data to calibrate the method, calculations are presented for the pressure endwall, and suction surfaces of a stationary cascade and for the pressure surface of a rotating turbine blade. The results strongly indicate that the 3-D boundary layer analysis can give good predictions of the flow field, loss, and heat transfer on the pressure, suction, and endwall surface of a gas turbine passage.

Computation of subsonic flow around airfoil systems with multiple separation

NASA Technical Reports Server (NTRS)

Jacob, K.

1982-01-01

A numerical method for computing the subsonic flow around multi-element airfoil systems was developed, allowing for flow separation at one or more elements. Besides multiple rear separation also sort bubbles on the upper surface and cove bubbles can approximately be taken into account. Also, compressibility effects for pure subsonic flow are approximately accounted for. After presentation the method is applied to several examples and improved in some details. Finally, the present limitations and desirable extensions are discussed.

A Navier-Stokes solution of the three-dimensional viscous compressible flow in a centrifugal compressor impeller

NASA Technical Reports Server (NTRS)

Harp, J. L., Jr.

1977-01-01

A two-dimensional time-dependent computer code was utilized to calculate the three-dimensional steady flow within the impeller blading. The numerical method is an explicit time marching scheme in two spatial dimensions. Initially, an inviscid solution is generated on the hub blade-to-blade surface by the method of Katsanis and McNally (1973). Starting with the known inviscid solution, the viscous effects are calculated through iteration. The approach makes it possible to take into account principal impeller fluid-mechanical effects. It is pointed out that the second iterate provides a complete solution to the three-dimensional, compressible, Navier-Stokes equations for flow in a centrifugal impeller. The problems investigated are related to the study of a radial impeller and a backswept impeller.

Diffraction of a shock wave by a compression corner; regular and single Mach reflection

NASA Technical Reports Server (NTRS)

Vijayashankar, V. S.; Kutler, P.; Anderson, D.

1976-01-01

The two dimensional, time dependent Euler equations which govern the flow field resulting from the injection of a planar shock with a compression corner are solved with initial conditions that result in either regular reflection or single Mach reflection of the incident planar shock. The Euler equations which are hyperbolic are transformed to include the self similarity of the problem. A normalization procedure is employed to align the reflected shock and the Mach stem as computational boundaries to implement the shock fitting procedure. A special floating fitting scheme is developed in conjunction with the method of characteristics to fit the slip surface. The reflected shock, the Mach stem, and the slip surface are all treated as harp discontinuities, thus, resulting in a more accurate description of the inviscid flow field. The resulting numerical solutions are compared with available experimental data and existing first-order, shock-capturing numerical solutions.

A Unified View of Global Instabilities of Compressible Flow Over Open Cavities

DTIC Science & Technology

2005-06-30

the early work of Rossiter [3], have treated the shear-layer emanating from the upstream comer of the cavity in isolation ( using parallel flow... using a domain-decomposition method. The code has optional equation sets to solve either (i) nonlinear Navier-Stokes, (ii) Navier-Stokes equations...early experments of Maull and East [15]. They used oil flow visualization of surface streamlines on the cavity bottom to show the existence, under certain

Distribution of soap molecules in flowing soap films

NASA Astrophysics Data System (ADS)

Kim, Ildoo; Sane, Aakash; Mandre, Shreyas

Flowing soap films are useful tools to simulate two-dimensional flows. The Marangoni elasticity due to the presence of soap molecules not only stabilizes the soap film but also imparts it compressibility to the two-dimensional flow in the soap film. Therefore, it is desirable to measure the surface concentration cs of soap molecules to understand the physics flowing soap films. In this study, we present an indirect measurement of cs, by making a direct measurement of the surface tension and the Marangoni elasticity. Using a two-stage model for soap distribution in the flows, the range of cs is calculated for different thickness and the soap solution concentration. Our model shows that the soap film will have the same cs for the range of parameters in popular use and in agreements with experimental data.

Passive margins getting squeezed in the mantle convection vice

NASA Astrophysics Data System (ADS)

Husson, Laurent; Yamato, Philippe; Becker, Thorsten; Pedoja, Kevin

2013-04-01

Quaternary coastal geomorphology reveals that passive margins underwent wholesale uplift at least during the glacial cycle. In addition, these not-so-passive margins often exhibit long term exhumation and tectonic inversion, which suggest that compression and tectonic shortening could be the mechanism that triggers their overall uplift. We speculate that the compression in the lithosphere gradually increased during the Cenozoic. The many mountain belts at active margins that accompany this event readily witness this increase. Less clear is how that compression increase affects passive margins. In order to address this issue, we design minimalist 2D viscous models to quantify the impact of plate collision on the stress regime. In these models, a sluggish plate is disposed on a less viscous mantle. It is driven by a "mantle conveyor belt" alternatively excited by lateral shear stresses that represent a downwelling on one side, an upwelling on the other side, or both simultaneously. The lateral edges of the plate are either free or fixed, respectively representing the cases of free convergence and collision. In practice, it dramatically changes the upper boundary condition for mantle circulation and subsequently, for the stress field. The flow pattern transiently evolves almost between two end-members, starting from a situation close to a Couette flow to a pattern that looks like a Poiseuille flow with an almost null velocity at the surface (though in the models, the horizontal velocity at the surface is not strictly null, as the lithosphere deforms). In the second case, the lithosphere is highly stressed horizontally and deforms. For an equivalent bulk driving force, compression increases drastically at passive margins if upwellings are active because they push plates towards the collision. Conversely, if only downwellings are activated, compression occurs on one half of the plate and extension on the other half, because only the downwelling is pulling the plate. Thus, active upwellings underneath oceanic plates are required to explain compression at passive margins. This conclusion is corroborated by "real-Earth" 3D spherical models, wherein the flow is alternatively driven by density anomalies inferred from seismic tomography -and therefore include both downwellings at subduction zones and upwellings above the superswells- and density anomalies that correspond to subducting slabs only. While the second scenario mostly compresses the active margins of upper plates and leave other areas at rest, the first scenario efficiently compresses passive margins where the geological record reveals their uplift, exhumation, and tectonic inversion.

A volume-filtered formulation to capture particle-shock interactions in multiphase compressible flows

NASA Astrophysics Data System (ADS)

Shallcross, Gregory; Capecelatro, Jesse

2017-11-01

Compressible particle-laden flows are common in engineering systems. Applications include but are not limited to water injection in high-speed jet flows for noise suppression, rocket-plume surface interactions during planetary landing, and explosions during coal mining operations. Numerically, it is challenging to capture these interactions due to the wide range of length and time scales. Additionally, there are many forms of the multiphase compressible flow equations with volume fraction effects, some of which are conflicting in nature. The purpose of this presentation is to develop the capability to accurately capture particle-shock interactions in systems with a large number of particles from dense to dilute regimes. A thorough derivation of the volume filtered equations is presented. The volume filtered equations are then implemented in a high-order, energy-stable Eulerian-Lagrangian framework. We show this framework is capable of decoupling the fluid mesh from the particle size, enabling arbitrary particle size distributions in the presence of shocks. The proposed method is then assessed against particle-laden shock tube data. Quantities of interest include fluid-phase pressure profiles and particle spreading rates. The effect of collisions in 2D and 3D are also evaluated.

Numerical solutions of Navier-Stokes equations for a Butler wing

NASA Technical Reports Server (NTRS)

Abolhassani, J. S.; Tiwari, S. N.

1985-01-01

The flow field is simulated on the surface of a given delta wing (Butler wing) at zero incident in a uniform stream. The simulation is done by integrating a set of flow field equations. This set of equations governs the unsteady, viscous, compressible, heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. These equations are solved by the finite difference method, and results obtained for high-speed freestream conditions are compared with theoretical and experimental results. In this study, the Navier-Stokes equations are solved numerically. These equations are unsteady, compressible, viscous, and three-dimensional without neglecting any terms. The time dependency of the governing equations allows the solution to progress naturally for an arbitrary initial initial guess to an asymptotic steady state, if one exists. The equations are transformed from physical coordinates to the computational coordinates, allowing the solution of the governing equations in a rectangular parallel-piped domain. The equations are solved by the MacCormack time-split technique which is vectorized and programmed to run on the CDC VPS 32 computer.

Application of Compressible Volume of Fluid Model in Simulating the Impact and Solidification of Hollow Spherical ZrO2 Droplet on a Surface

NASA Astrophysics Data System (ADS)

Safaei, Hadi; Emami, Mohsen Davazdah; Jazi, Hamidreza Salimi; Mostaghimi, Javad

2017-12-01

Applications of hollow spherical particles in thermal spraying process have been developed in recent years, accompanied by attempts in the form of experimental and numerical studies to better understand the process of impact of a hollow droplet on a surface. During such process, volume and density of the trapped gas inside droplet change. The numerical models should be able to simulate such changes and their consequent effects. The aim of this study is to numerically simulate the impact of a hollow ZrO2 droplet on a flat surface using the volume of fluid technique for compressible flows. An open-source, finite-volume-based CFD code was used to perform the simulations, where appropriate subprograms were added to handle the studied cases. Simulation results were compared with the available experimental data. Results showed that at high impact velocities ( U 0 > 100 m/s), the compression of trapped gas inside droplet played a significant role in the impact dynamics. In such velocities, the droplet splashed explosively. Compressibility effects result in a more porous splat, compared to the corresponding incompressible model. Moreover, the compressible model predicted a higher spread factor than the incompressible model, due to planetary structure of the splat.

Compressible-Incompressible Two-Phase Flows with Phase Transition: Model Problem

NASA Astrophysics Data System (ADS)

Watanabe, Keiichi

2017-12-01

We study the compressible and incompressible two-phase flows separated by a sharp interface with a phase transition and a surface tension. In particular, we consider the problem in R^N , and the Navier-Stokes-Korteweg equations is used in the upper domain and the Navier-Stokes equations is used in the lower domain. We prove the existence of R -bounded solution operator families for a resolvent problem arising from its model problem. According to Göts and Shibata (Asymptot Anal 90(3-4):207-236, 2014), the regularity of ρ _+ is W^1_q in space, but to solve the kinetic equation: u_Γ \\cdot n_t = [[ρ u

Numerical analysis of exhaust jet secondary combustion in hypersonic flow field

NASA Astrophysics Data System (ADS)

Yang, Tian-Peng; Wang, Jiang-Feng; Zhao, Fa-Ming; Fan, Xiao-Feng; Wang, Yu-Han

2018-05-01

The interaction effect between jet and control surface in supersonic and hypersonic flow is one of the key problems for advanced flight control system. The flow properties of exhaust jet secondary combustion in a hypersonic compression ramp flow field were studied numerically by solving the Navier-Stokes equations with multi-species and combustion reaction effects. The analysis was focused on the flow field structure and the force amplification factor under different jet conditions. Numerical results show that a series of different secondary combustion makes the flow field structure change regularly, and the temperature increases rapidly near the jet exit.

The rheology and composition of cryovolcanic flows on icy satellites

NASA Technical Reports Server (NTRS)

Kargel, Jeffrey S.

1993-01-01

The rheologic properties of terrestrial lavas have been related to morphologic features of their flows, such as levees, banked surfaces, multilobate structures, and compressible folds. These features also have been used to determine rheologies and constrain the compositions of extraterrestrial flows. However, with rare exceptions, such features are not resolvable in Voyager images of the satellites of outer planets. Often only flow length and edge thickness of cryovolcanic flows can be measured reasonably accurately from Voyager images. The semiempirical lava-flow model presented here is a renewed effort to extract useful information from such measurements.

Theoretical prediction of airplane stability derivatives at subcritical speeds

NASA Technical Reports Server (NTRS)

Tulinius, J.; Clever, W.; Nieman, A.; Dunn, K.; Gaither, B.

1973-01-01

The theoretical development and application is described of an analysis for predicting the major static and rotary stability derivatives for a complete airplane. The analysis utilizes potential flow theory to compute the surface flow fields and pressures on any configuration that can be synthesized from arbitrary lifting bodies and nonplanar thick lifting panels. The pressures are integrated to obtain section and total configuration loads and moments due side slip, angle of attack, pitching motion, rolling motion, yawing motion, and control surface deflection. Subcritical compressibility is accounted for by means of the Gothert similarity rule.

Heat transfer, velocity-temperature correlation, and turbulent shear stress from Navier-Stokes computations of shock wave/turbulent boundary layer interaction flows

NASA Technical Reports Server (NTRS)

Wang, C. R.; Hingst, W. R.; Porro, A. R.

1991-01-01

The properties of 2-D shock wave/turbulent boundary layer interaction flows were calculated by using a compressible turbulent Navier-Stokes numerical computational code. Interaction flows caused by oblique shock wave impingement on the turbulent boundary layer flow were considered. The oblique shock waves were induced with shock generators at angles of attack less than 10 degs in supersonic flows. The surface temperatures were kept at near-adiabatic (ratio of wall static temperature to free stream total temperature) and cold wall (ratio of wall static temperature to free stream total temperature) conditions. The computational results were studied for the surface heat transfer, velocity temperature correlation, and turbulent shear stress in the interaction flow fields. Comparisons of the computational results with existing measurements indicated that (1) the surface heat transfer rates and surface pressures could be correlated with Holden's relationship, (2) the mean flow streamwise velocity components and static temperatures could be correlated with Crocco's relationship if flow separation did not occur, and (3) the Baldwin-Lomax turbulence model should be modified for turbulent shear stress computations in the interaction flows.

Progress Towards a Cartesian Cut-Cell Method for Viscous Compressible Flow

NASA Technical Reports Server (NTRS)

Berger, Marsha; Aftosmis, Michael J.

2012-01-01

We present preliminary development of an approach for simulating high Reynolds number steady compressible flow in two space dimensions using a Cartesian cut-cell finite volume method. We consider both laminar and turbulent flow with both low and high cell Reynolds numbers near the wall. The approach solves the full Navier-Stokes equations in all cells, and uses a wall model to address the resolution requirements near boundaries and to mitigate mesh irregularities in cut cells. We present a quadratic wall model for low cell Reynolds numbers. At high cell Reynolds numbers, the quadratic is replaced with a newly developed analytic wall model stemming from solution of a limiting form of the Spalart-Allmaras turbulence model which features a forward evaluation for flow velocity and exactly matches characteristics of the SA turbulence model in the field. We develop multigrid operators which attain convergence rates similar to inviscid multigrid. Investigations focus on preliminary verification and validation of the method. Flows over flat plates and compressible airfoils show good agreement with both theoretical results and experimental data. Mesh convergence studies on sub- and transonic airfoil flows show convergence of surface pressures with wall spacings as large as approx.0.1% chord. With the current analytic wall model, one or two additional refinements near the wall are required to obtain mesh converged values of skin friction.

Squeezing the muscle: compression clothing and muscle metabolism during recovery from high intensity exercise.

PubMed

Sperlich, Billy; Born, Dennis-Peter; Kaskinoro, Kimmo; Kalliokoski, Kari K; Laaksonen, Marko S

2013-01-01

The purpose of this experiment was to investigate skeletal muscle blood flow and glucose uptake in m. biceps (BF) and m. quadriceps femoris (QF) 1) during recovery from high intensity cycle exercise, and 2) while wearing a compression short applying ~37 mmHg to the thigh muscles. Blood flow and glucose uptake were measured in the compressed and non-compressed leg of 6 healthy men by using positron emission tomography. At baseline blood flow in QF (P = 0.79) and BF (P = 0.90) did not differ between the compressed and the non-compressed leg. During recovery muscle blood flow was higher compared to baseline in both compressed (P<0.01) and non-compressed QF (P<0.001) but not in compressed (P = 0.41) and non-compressed BF (P = 0.05; effect size = 2.74). During recovery blood flow was lower in compressed QF (P<0.01) but not in BF (P = 0.26) compared to the non-compressed muscles. During baseline and recovery no differences in blood flow were detected between the superficial and deep parts of QF in both, compressed (baseline P = 0.79; recovery P = 0.68) and non-compressed leg (baseline P = 0.64; recovery P = 0.06). During recovery glucose uptake was higher in QF compared to BF in both conditions (P<0.01) with no difference between the compressed and non-compressed thigh. Glucose uptake was higher in the deep compared to the superficial parts of QF (compression leg P = 0.02). These results demonstrate that wearing compression shorts with ~37 mmHg of external pressure reduces blood flow both in the deep and superficial regions of muscle tissue during recovery from high intensity exercise but does not affect glucose uptake in BF and QF.

Normal- and oblique-shock flow parameters in equilibrium air including attached-shock solutions for surfaces at angles of attack, sweep, and dihedral

NASA Technical Reports Server (NTRS)

Hunt, J. L.; Souders, S. W.

1975-01-01

Normal- and oblique-shock flow parameters for air in thermochemical equilibrium are tabulated as a function of shock angle for altitudes ranging from 15.24 km to 91.44 km in increments of 7.62 km at selected hypersonic speeds. Post-shock parameters tabulated include flow-deflection angle, velocity, Mach number, compressibility factor, isentropic exponent, viscosity, Reynolds number, entropy difference, and static pressure, temperature, density, and enthalpy ratios across the shock. A procedure is presented for obtaining oblique-shock flow properties in equilibrium air on surfaces at various angles of attack, sweep, and dihedral by use of the two-dimensional tabulations. Plots of the flow parameters against flow-deflection angle are presented at altitudes of 30.48, 60.96, and 91.44 km for various stream velocities.

Rheology of the lithosphere and the folding caused by horizontal compression

NASA Astrophysics Data System (ADS)

Birger, B. I.

2015-05-01

The laboratory tests of rock specimens show that transient creep, at which deformations increase with time whereas strain rate decreases occurs when creep strains are sufficiently small. Since plate tectonics only permits small deformations in the lithospheric plates, the creep of the lithosphere is transient (non-steady-state). In this work, we study how the rheology of the lithosphere that possesses elasticity, brittleness (pseudo-plasticity), and creep affects the folding in the Earth's crust. Folding is caused by horizontal compression that results from the collision between the lithospheric plates. The effective viscosity characterizing the transient creep is lower than in the case of a steady-state creep and depends on the characteristic time of the considered process. The allowance for transient creep gives the distribution of the rheological properties of the horizontally compressed lithosphere in which the upper crust is brittle, whereas the lower crust and mantle lithosphere are dominated by transient creep. It is shown that the flows that arise in the lithosphere due to the instability under horizontal compression and cause folding are small-scale. These flows are concentrated in the upper brittle crust, they determine the short-wave Earth's surface topography, penetrate into the lower, creep-dominated crust to a shallow depth, and do not penetrate into the mantle. Therefore, these flows do not deform the Moho.

Finite element flow analysis; Proceedings of the Fourth International Symposium on Finite Element Methods in Flow Problems, Chuo University, Tokyo, Japan, July 26-29, 1982

NASA Astrophysics Data System (ADS)

Kawai, T.

Among the topics discussed are the application of FEM to nonlinear free surface flow, Navier-Stokes shallow water wave equations, incompressible viscous flows and weather prediction, the mathematical analysis and characteristics of FEM, penalty function FEM, convective, viscous, and high Reynolds number FEM analyses, the solution of time-dependent, three-dimensional and incompressible Navier-Stokes equations, turbulent boundary layer flow, FEM modeling of environmental problems over complex terrain, and FEM's application to thermal convection problems and to the flow of polymeric materials in injection molding processes. Also covered are FEMs for compressible flows, including boundary layer flows and transonic flows, hybrid element approaches for wave hydrodynamic loadings, FEM acoustic field analyses, and FEM treatment of free surface flow, shallow water flow, seepage flow, and sediment transport. Boundary element methods and FEM computational technique topics are also discussed. For individual items see A84-25834 to A84-25896

Experimental investigations of the time and flow-direction responses of shear-stress-sensitive liquid crystal coatings

NASA Technical Reports Server (NTRS)

Reda, Daniel C.; Muratore, Joseph J., Jr.; Heineck, James T.

1993-01-01

Time and flow-direction responses of shearstress-sensitive liquid crystal coatings were explored experimentally. For the time-response experiments, coatings were exposed to transient, compressible flows created during the startup and off-design operation of an injector-driven supersonic wind tunnel. Flow transients were visualized with a focusing Schlieren system and recorded with a 1000 frame/sec color video camera. Liquid crystal responses to these changing-shear environments were then recorded with the same video system, documenting color-play response times equal to, or faster than, the time interval between sequential frames (i.e., 1 millisecond). For the flow-direction experiments, a planar test surface was exposed to equal-magnitude and known-direction surface shear stresses generated by both normal and tangential subsonic jet-impingement flows. Under shear, the sense of the angular displacement of the liquid crystal dispersed (reflected) spectrum was found to be a function of the instantaneous direction of the applied shear. This technique thus renders dynamic flow reversals or flow divergences visible over entire test surfaces at image recording rates up to 1 KHz. Extensions of the technique to visualize relatively small changes in surface shear stress direction appear feasible.

Rudder/Fin Seal Investigations for the X-38 Re-Entry Vehicle

NASA Technical Reports Server (NTRS)

Dunlap, Patrick H., Jr.; Steinetz, Bruce M.; Curry, Donald M.

2000-01-01

NASA is currently developing the X-38 vehicle that will be used to demonstrate the technologies required for a crew return vehicle (CRV) for the International Space Station. The X-38 control surfaces require high temperature seals to limit hot gas ingestion and transfer of heat to underlying low-temperature structures to prevent over-temperature of these structures and possible loss of the vehicle. This paper presents results for thermal analyses and flow and compression tests conducted on as-received and thermally exposed seals for the rudder/fin location of the X-38. A thermal analysis of the rudder/fin dual seal assembly based on representative heating rates on the windward surface of the rudder/fin area predicted a peak seal temperature of 1900 F. The temperature-exposed seals were heated in a compressed state at 1900 F corresponding to the predicted peak temperature. Room temperature compression tests were performed to determine load versus linear compression, preload, contact area, stiffness, and resiliency characteristics for the as-received and temperature-exposed seals. Temperature exposure resulted in permanent set and loss of resiliency in these seals. Unit loads and contact pressures for the seals were below the 5 lb/in. and 10 psi limits set to limit the loads on the Shuttle thermal tiles that the seals seal against in the rudder/fin location. Measured seal flow rates for a double seal were about 4.5 times higher than the preliminary seal flow goal. The seal designs examined in this study are expected to be able to endure the high temperatures that they will be exposed to for a single-use life. Tests performed herein combined with future analyses, arc jet tests, and scrubbing tests will be used to select the final seal design for this application.

Rudder/Fin Seal Investigations for the X-38 Re-Entry Vehicle

NASA Technical Reports Server (NTRS)

Dunlap, Patrick H., Jr.; Steinetz, Bruce M.; Curry, Donald M.

2000-01-01

NASA is currently developing the X-38 vehicle that will be used to demonstrate the technologies required for a crew return vehicle (CRV) for the International Space Station. The X-38 control surfaces require high temperature seals to limit hot gas ingestion and transfer of heat to underlying low-temperature structures to prevent over-temperature of these structures and possible loss of the vehicle. This paper presents results for thermal analyses and flow and compression tests conducted on as-received and thermally exposed seals for the rudder/fin location of the X-38. A thermal analysis of the rudder/fin dual seal assembly based on representative heating rates on the windward surface of the rudder/fin area predicted a peak seal temperature of 1900 F. The temperature-exposed seals were heated in a compressed state at 1900 F corresponding to the predicted peak temperature. Room temperature compression tests were performed to determine load versus linear compression, preload, contact area, stiffness, and resiliency characteristics for the as-received and temperature-exposed seals. Temperature exposure resulted in permanent set and loss of resiliency in these seals. Unit loads and contact pressures for the seals were below the five pounds/inch and ten psi limits set to limit the loads on the Shuttle thermal tiles that the seals seal against in the rudder/fin location. Measured seal flow rates for a double seal were about 4.5 times higher than the preliminary seal flow goal. The seal designs examined in this study are expected to be able to endure the high temperatures that they will be exposed to for a single-use life. Tests performed herein combined with future analyses, arc jet tests, and scrubbing tests will be used to select the final seal design for this application.

Transonic Flow Past Cone Cylinders

NASA Technical Reports Server (NTRS)

Solomon, George E

1955-01-01

Experimental results are presented for transonic flow post cone-cylinder, axially symmetric bodies. The drag coefficient and surface Mach number are studied as the free-stream Mach number is varied and, wherever possible, the experimental results are compared with theoretical predictions. Interferometric results for several typical flow configurations are shown and an example of shock-free supersonic-to-subsonic compression is experimentally demonstrated. The theoretical problem of transonic flow past finite cones is discussed briefly and an approximate solution of the axially symmetric transonic equations, valid for a semi-infinite cone, is presented.

Turbulence in Compressible Flows

NASA Technical Reports Server (NTRS)

1997-01-01

Lecture notes for the AGARD Fluid Dynamics Panel (FDP) Special Course on 'Turbulence in Compressible Flows' have been assembled in this report. The following topics were covered: Compressible Turbulent Boundary Layers, Compressible Turbulent Free Shear Layers, Turbulent Combustion, DNS/LES and RANS Simulations of Compressible Turbulent Flows, and Case Studies of Applications of Turbulence Models in Aerospace.

Modeling of Compressible Flow with Friction and Heat Transfer Using the Generalized Fluid System Simulation Program (GFSSP)

NASA Technical Reports Server (NTRS)

Bandyopadhyay, Alak; Majumdar, Alok

2007-01-01

The present paper describes the verification and validation of a quasi one-dimensional pressure based finite volume algorithm, implemented in Generalized Fluid System Simulation Program (GFSSP), for predicting compressible flow with friction, heat transfer and area change. The numerical predictions were compared with two classical solutions of compressible flow, i.e. Fanno and Rayleigh flow. Fanno flow provides an analytical solution of compressible flow in a long slender pipe where incoming subsonic flow can be choked due to friction. On the other hand, Raleigh flow provides analytical solution of frictionless compressible flow with heat transfer where incoming subsonic flow can be choked at the outlet boundary with heat addition to the control volume. Nonuniform grid distribution improves the accuracy of numerical prediction. A benchmark numerical solution of compressible flow in a converging-diverging nozzle with friction and heat transfer has been developed to verify GFSSP's numerical predictions. The numerical predictions compare favorably in all cases.

Modeling of Ice Flow and Internal Layers Along a Flow Line Through Swiss Camp in West Greenland

NASA Technical Reports Server (NTRS)

Wang, W. L.; Zwally, H. Jay; Abdalati, W.; Luo, S.; Koblinsky, Chester J. (Technical Monitor)

2001-01-01

An anisotropic ice flow line model is applied to a flow line through Swiss Camp (69.57 N, 49.28 W) in West Greenland to estimate the dates of internal layers detected by Radio-Echo Sounding measurements. The effect of an anisotropic ice fabric on ice flow is incorporated into the steady state flow line model. The stress-strain rate relationship for anisotropic ice is characterized by an enhancement factor based on the laboratory observations of ice deformation under combined compression and shear stresses. By using present-day data of accumulation rate, surface temperature, surface elevation and ice thickness along the flow line as model inputs, a very close agreement is found between the isochrones generated from the model and the observed internal layers with confirmed dates. The results indicate that this part of Greenland ice sheet is primarily in steady state.

Effects of compressibility on the temperature jump at the interface of layered, spherical-shell convection

NASA Technical Reports Server (NTRS)

Yen, David A.; Zhang, Shuxia; Langenberger, Sherri E.

1988-01-01

Large temperature jumps at the interface of layered convection are important to the argument used against the likelihood of separate circulations in the upper and lower mantles. This problem was studied within the framework of a compressible, constant viscosity spherical-shell model. Both mechanical and thermal coupling configurations are considered. Although the temperature jumps are reduced by compressibility, their magnitudes remain quite large, in the case of mechanical coupling. For thermal coupling, the temperature jumps become smaller but still are substantial, between 500 to 1000 C. In layered spherical-shell convection, flows in the lower mantle are several times greater than the surface velocities.

Numerical Simulation of a High-Lift Configuration Embedded with High Momentum Fluidic Actuators

NASA Technical Reports Server (NTRS)

Vatsa, Veer N.; Duda, Benjamin; Fares, Ehab; Lin, John C.

2016-01-01

Numerical simulations have been performed for a vertical tail configuration with deflected rudder. The suction surface of the main element of this configuration, just upstream of the hinge line, is embedded with an array of 32 fluidic actuators that produce oscillating sweeping jets. Such oscillating jets have been found to be very effective for flow control applications in the past. In the current paper, a high-fidelity computational fluid dynamics (CFD) code known as the PowerFLOW R code is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. A fully compressible version of the PowerFLOW R code valid for high speed flows is used for the present simulations to accurately represent the transonic flow regimes encountered in the flow field due to the actuators operating at higher mass flow (momentum) rates required to mitigate reverse flow regions on a highly-deflected rudder surface. The computed results for the surface pressure and integrated forces compare favorably with measured data. In addition, numerical solutions predict the correct trends in forces with active flow control compared to the no control case. The effect of varying the rudder deflection angle on integrated forces and surface pressures is also presented.

Three-dimensional numerical simulation for plastic injection-compression molding

NASA Astrophysics Data System (ADS)

Zhang, Yun; Yu, Wenjie; Liang, Junjie; Lang, Jianlin; Li, Dequn

2018-03-01

Compared with conventional injection molding, injection-compression molding can mold optical parts with higher precision and lower flow residual stress. However, the melt flow process in a closed cavity becomes more complex because of the moving cavity boundary during compression and the nonlinear problems caused by non-Newtonian polymer melt. In this study, a 3D simulation method was developed for injection-compression molding. In this method, arbitrary Lagrangian- Eulerian was introduced to model the moving-boundary flow problem in the compression stage. The non-Newtonian characteristics and compressibility of the polymer melt were considered. The melt flow and pressure distribution in the cavity were investigated by using the proposed simulation method and compared with those of injection molding. Results reveal that the fountain flow effect becomes significant when the cavity thickness increases during compression. The back flow also plays an important role in the flow pattern and redistribution of cavity pressure. The discrepancy in pressures at different points along the flow path is complicated rather than monotonically decreased in injection molding.

Dynamic interactions between glacier and glacial lake in the Bhutan Himalaya

NASA Astrophysics Data System (ADS)

Tsutaki, S.; Fujita, K.; Yamaguchi, S.; Sakai, A.; Nuimura, T.; Komori, J.; Takenaka, S.; Tshering, P.

2012-04-01

A number of supraglacial lakes formed on the termini of debris-covered glaciers in the Bhutan Himalaya as a result of glacier retreat due to climate change. The terminal part of the lake-terminating glaciers flow faster than that of the land-terminating glaciers because the basal ice motion is enhanced by high subglacial water pressure generated by lake water. Increased ice flux caused by the accelerated glacier flow could be dissipated through the calving process which reduced the glacier thickness. It is important to understand the interaction between lake formation and glacier dynamics. Although glacier flow velocity has been measured by remote-sensing analysis in several regions of the Himalayas, glacier thinning rates have not been observed by neither in-situ nor remote-sensing approaches. The lack of field data raises limitation to interpretations for glacier dynamics. We investigate the influence of the presence/absence of glacial lakes on glacier dynamics and changes in surface elevation. We study two debris-covered glaciers in the Lunana region, the Bhutan Himalaya. Thorthormi Glacier is a land-terminating glacier with some supraglacial lakes while Lugge Glacier is a lake-terminating glaciers. We surveyed the surface elevation of debris-covered areas of the two glaciers in 2004 and 2011 by a differential GPS. Change in surface elevation of the lake-terminating Lugge Glacier (-5.4--2.4 m yr-1) was much more negative than that of the land-terminating Thorthormi Glacier (-3.3-0.6 m yr-1). Surface flow speed of the Thorthormi Glacier measured during 2002-2004 was faster in the upper reaches (~90 m yr-1) and reduced toward the downstream (40 m yr-1). In contrast, the surface flow speed at the Lugge Glacier measured in the same periods was 40-55 m yr-1 and the greatest at the lower most part. Observed spatial distribution of surface flow velocity at both glaciers were evaluated by a two-dimensional numerical flow model. Calculated emergence velocities are 1.9-18.8 m yr-1 at the Thorthormi Glacier while -12.0-2.7 m yr-1 at the Lugge Glacier. This result suggests that decreasing in flow velocity towards the terminus in the Thorthormi Glacier causes compressive flow. It suggests that the compressive flow of the Thorthormi Glacier counterbalanced surface melting, resulting in inhibition of the surface lowering. In contrast, the extensional flow of the Lugge Glacier accelerated the surface lowering. Speed up of glacier terminus induced extensional flow regime causes the thinning of ice and increase in basal motion, which will lead to further flow acceleration. Such positive feedbacks have been found over the ice streams in the polar ice sheets. In this study we showed the observational evidences, in which the similar feedbacks make contrast the terminus behaviors of glaciers in the Bhutan Himalaya. If the supraglacial lake on Thorthormi Glacier expanded, the surface lowering may be accelerated in the future.

Theoretical fluid dynamics

NASA Astrophysics Data System (ADS)

Shivamoggi, B. K.

This book is concerned with a discussion of the dynamical behavior of a fluid, and is addressed primarily to graduate students and researchers in theoretical physics and applied mathematics. A review of basic concepts and equations of fluid dynamics is presented, taking into account a fluid model of systems, the objective of fluid dynamics, the fluid state, description of the flow field, volume forces and surface forces, relative motion near a point, stress-strain relation, equations of fluid flows, surface tension, and a program for analysis of the governing equations. The dynamics of incompressible fluid flows is considered along with the dynamics of compressible fluid flows, the dynamics of viscous fluid flows, hydrodynamic stability, and dynamics of turbulence. Attention is given to the complex-variable method, three-dimensional irrotational flows, vortex flows, rotating flows, water waves, applications to aerodynamics, shock waves, potential flows, the hodograph method, flows at low and high Reynolds numbers, the Jeffrey-Hamel flow, and the capillary instability of a liquid jet.

LOW-VELOCITY COMPRESSIBLE FLOW THEORY

EPA Science Inventory

The widespread application of incompressible flow theory dominates low-velocity fluid dynamics, virtually preventing research into compressible low-velocity flow dynamics. Yet, compressible solutions to simple and well-defined flow problems and a series of contradictions in incom...

Influence of leading edge bluntness on hypersonic flow in a generic internal-compression inlet

NASA Astrophysics Data System (ADS)

Borovoy, V.; Egorov, I.; Mosharov, V.; Radchenko, V.; Skuratov, A.; Struminskaya, I.

2015-06-01

Flow and heat transfer inside a generic inlet are investigated experimentally. The cross section of the inlet is rectangular. The inlet is installed on a flat plat at a significant distance from the leading edge. The experiments are performed in TsAGI wind tunnel UT-1M working in the Ludwieg tube mode at Mach number M∞ = 5 and Reynolds numbers (based on the plate length L = 320 mm) Re∞L = 23 · 106 and 13 · 106. Steady flow duration is 40 ms. Optical panoramic methods are used for investigation of flow outside and inside the inlet as well. For this purpose, the cowl and one of two compressing wedges are made of a transparent material. Heat flux distribution is measured by thin luminescent Temperature Sensitive Paint (TSP). Surface flow and shear stress visualization is performed by viscous oil containing luminophor particles. The investigation shows that at high contraction ratio of the inlet, an increase of plate or cowl bluntness to some critical value leads to sudden change of the flow structure.

A Combined Experimental/Computational Study of Flow in Turbine Blade Cooling Passage

NASA Technical Reports Server (NTRS)

Tse, D. G. N.; Kreskovsky, J. P.; Shamroth, S. J.; Mcgrath, D. B.

1994-01-01

Laser velocimetry was utilized to map the velocity field in a serpentine turbine blade cooling passage at Reynolds and Rotation numbers of up to 25.000 and 0.48. These results were used to assess the combined influence of passage curvature and Coriolis force on the secondary velocity field generated. A Navier-Stokes code (NASTAR) was validated against incompressible test data and then used to simulate the effect of buoyancy. The measurements show a net convection from the low pressure surface to high pressure surface. The interaction of the secondary flows induced by the turns and rotation produces swirl at the turns, which persisted beyond 2 hydraulic diameters downstream of the turns. The incompressible flow field predictions agree well with the measured velocities. With radially outward flow, the buoyancy force causes a further increase in velocity on the high pressure surface and a reduction on the low pressure surface. The results were analyzed in relation to the heat transfer measurements of Wagner et al. (1991). Predicted heat transfer is enhanced on the high pressure surfaces and in turns. The incompressible flow simulation underpredicts heat transfer in these locations. Improvements observed in compressible flow simulation indicate that the buoyancy force may be important.

Unsteady viscous effects in the flow over an oscillating surface. [mathematical model

NASA Technical Reports Server (NTRS)

Lerner, J. I.

1972-01-01

A theoretical model for the interaction of a turbulent boundary layer with an oscillating wavy surface over which a fluid is flowing is developed, with an application to wind-driven water waves and to panel flutter in low supersonic flow. A systematic methodology is developed to obtain the surface pressure distribution by considering separately the effects on the perturbed flow of a mean shear velocity profile, viscous stresses, the turbulent Reynolds stresses, compressibility, and three-dimensionality. The inviscid theory is applied to the wind-water wave problem by specializing to traveling-wave disturbances, and the pressure magnitude and phase shift as a function of the wave phase speed are computed for a logarithmic mean velocity profile and compared with inviscid theory and experiment. The results agree with experimental evidence for the stabilization of the panel motion due to the influence of the unsteady boundary layer.

Effects of compressibility on turbulent relative particle dispersion

NASA Astrophysics Data System (ADS)

Shivamoggi, Bhimsen K.

2016-08-01

In this paper, phenomenological developments are used to explore the effects of compressibility on the relative particle dispersion (RPD) in three-dimensional (3D) fully developed turbulence (FDT). The role played by the compressible FDT cascade physics underlying this process is investigated. Compressibility effects are found to lead to reduction of RPD, development of the ballistic regime and particle clustering, corroborating the laboratory experiment and numerical simulation results (Cressman J. R. et al., New J. Phys., 6 (2004) 53) on the motion of Lagrangian tracers on a surface flow that constitutes a 2D compressible subsystem. These formulations are developed from the scaling relations for compressible FDT and are validated further via an alternative dimensional/scaling development for compressible FDT similar to the one given for incompressible FDT by Batchelor and Townsend (Surveys in Mechanics (Cambridge University Press) 1956, p. 352). The rationale for spatial intermittency effects is legitimized via the nonlinear scaling dependence of RPD on the kinetic-energy dissipation rate.

Interactive computer graphics applications for compressible aerodynamics

NASA Technical Reports Server (NTRS)

Benson, Thomas J.

1994-01-01

Three computer applications have been developed to solve inviscid compressible fluids problems using interactive computer graphics. The first application is a compressible flow calculator which solves for isentropic flow, normal shocks, and oblique shocks or centered expansions produced by two dimensional ramps. The second application couples the solutions generated by the first application to a more graphical presentation of the results to produce a desk top simulator of three compressible flow problems: 1) flow past a single compression ramp; 2) flow past two ramps in series; and 3) flow past two opposed ramps. The third application extends the results of the second to produce a design tool which solves for the flow through supersonic external or mixed compression inlets. The applications were originally developed to run on SGI or IBM workstations running GL graphics. They are currently being extended to solve additional types of flow problems and modified to operate on any X-based workstation.

Development and Application of New Algorithms for the Simulation of Viscous Compressible Flows with Moving Bodies in Three Dimensions.

DTIC Science & Technology

1996-12-01

ranging from academic to industrial demonstrated the utility of the developed procedure for ab initio surface meshing from discrete data, such as...academic to industrial demonstrate the utility of the pro- hypersonic reentry problems, where ray-tracing based on posed procedure for ab initio surface...data input within industrial simulations. The origi- nal CAD dataset had over 500 surface patches, many All of the surface grids shown were obtained

Unsteady three-dimensional flow separation

NASA Technical Reports Server (NTRS)

Hui, W. H.

1988-01-01

A concise mathematical framework is constructed to study the topology of steady 3-D separated flows of an incompressible, or a compressible viscous fluid. Flow separation is defined by the existence of a stream surface which intersects with the body surface. The line of separation is itself a skin-friction line. Flow separation is classified as being either regular or singular, depending respectively on whether the line of separation contains only a finite number of singular points or is a singular line of the skin-friction field. The special cases of 2-D and axisymmetric flow separation are shown to be of singular type. In regular separation it is shown that a line of separation originates from a saddle point of separation of the skin-friction field and ends at nodal points of separation. Unsteady flow separation is defined relative to a coordinate system fixed to the body surface. It is shown that separation of an unsteady 3-D incompressible viscous flow at time t, when viewed from such a frame of reference, is topologically the same as that of the fictitious steady flow obtained by freezing the unsteady flow at the instant t. Examples are given showing effects of various forms of flow unsteadiness on flow separation.

Characterization of bulk and shear properties of basmati and non-basmati rice flour.

PubMed

Jan, Shumaila; Ghoroi, Chinmay; Saxena, Dharmesh Chandra

2018-01-01

Flours are often unstable in relation to their flow performance, which is evident when a free-flowing material ceases to flow and the processing, handling, and production parameters depend on the inherent powder characteristics and their bulk behaviour. The present study was conducted to compare the flowability of basmati and non-basmati rice flour affecting bulk handling, which could be related to its particle size, shape and surface roughness (measured by atomic force microscopy) as well as bulk and shear properties, depending upon the processing conditions. Particle size (171.1-171.9 μm) of both samples was not significantly different. However, the flowability of the non-basmati rice flour was significantly affected by its particle shape (circularity 0.487), surface roughness (124.23 nm) and compressibility (25.32%) in comparison to basmati rice flour (circularity 0.653, surface roughness 113.59 nm and compressibility 21.09%), making it more cohesive than basmati rice flour. Also, basic flow energy was significantly higher in non-basmati flour, thus requiring more energy (147.54 mJ) to flow than basmati rice flour (130.15 mJ). Overall, flowability was analysed by applying three different pressures (3, 6 and 9 kPa), among which non-basmati rice flour was found to be less flowable (flow function coefficient (FFC) 2.33 at 9 kPa) in comparison to basmati (FFC 3.35 at 9 kPa), making bulk handling difficult. This study could be useful in designing processing equipment, hoppers and silos for rice flour handling. © 2017 Society of Chemical Industry. © 2017 Society of Chemical Industry.

Investigation of Injector Slot Geometry on Curved-Diffuser Aerodynamic Performance

NASA Technical Reports Server (NTRS)

Silva, Odlanier

2004-01-01

The Compressor Branch vision is to be recognized as world-class leaders in research for fluid mechanics of compressors. Its mission is to conduct research and develop technology to advance the state of the art of compressors and transfer new technology to U.S. industries. Maintain partnerships with U.S. industries, universities, and other government organizations. Maintain a balance between customers focused and long range research. Flow control comprises enabling technologies to meet compression system performance requirements driven by emissions and fuel reduction goals (e.g., in UEET), missions (e.g., access-to-space), aerodynamically aggressive vehicle configurations (e.g., UAV and future blended wing body configurations with highly distorted inlets), and cost goals (e.g., in VAATE). The compression system requirements include increased efficiency, power-to-weight, and adaptability (i.e., robustness in terms of wide operability, distortion tolerance, and engine system health and reliability). The compressor flow control task comprises efforts to develop, demonstrate, and transfer adaptive flow control technology to industry to increase aerodynamic loading at current blade row loss levels, to enable adaptive1 y wide operability, and to develop plant models for adaptive compression systems. In this context, flow control is the controlled modification of a flow field by a deliberate means beyond the natural (uncontrolled) shaping of the solid surfaces that define the principal flow path. The objective of the compressor flow control task is to develop and apply techniques that control circulation, aerodynamic blockage, and entropy production in order to enhance the performance and operability of compression systems for advanced aero-propulsion applications. This summer I would be working with a curved-diffuser because it simulates what happens with flow in the stator blades in the compressor. With this experiment I will be doing some data analysis and parametric study of the injector slot geometries to get the best aerodynamic performance of it. This includes some data reduction, redesign and fast prototyping of the injector nozzle.

Effects of Medium Temperature and Industrial By-Products on the Key Hardened Properties of High Performance Concrete.

PubMed

Safiuddin, Md; Raman, Sudharshan N; Zain, Muhammad Fauzi Mohd

2015-12-10

The aim of the work reported in this article was to investigate the effects of medium temperature and industrial by-products on the key hardened properties of high performance concrete. Four concrete mixes were prepared based on a water-to-binder ratio of 0.35. Two industrial by-products, silica fume and Class F fly ash, were used separately and together with normal portland cement to produce three concrete mixes in addition to the control mix. The properties of both fresh and hardened concretes were examined in the laboratory. The freshly mixed concrete mixes were tested for slump, slump flow, and V-funnel flow. The hardened concretes were tested for compressive strength and dynamic modulus of elasticity after exposing to 20, 35 and 50 °C. In addition, the initial surface absorption and the rate of moisture movement into the concretes were determined at 20 °C. The performance of the concretes in the fresh state was excellent due to their superior deformability and good segregation resistance. In their hardened state, the highest levels of compressive strength and dynamic modulus of elasticity were produced by silica fume concrete. In addition, silica fume concrete showed the lowest level of initial surface absorption and the lowest rate of moisture movement into the interior of concrete. In comparison, the compressive strength, dynamic modulus of elasticity, initial surface absorption, and moisture movement rate of silica fume-fly ash concrete were close to those of silica fume concrete. Moreover, all concretes provided relatively low compressive strength and dynamic modulus of elasticity when they were exposed to 50 °C. However, the effect of increased temperature was less detrimental for silica fume and silica fume-fly ash concretes in comparison with the control concrete.

Effects of Medium Temperature and Industrial By-Products on the Key Hardened Properties of High Performance Concrete

PubMed Central

Safiuddin, Md.; Raman, Sudharshan N.; Zain, Muhammad Fauzi Mohd.

2015-01-01

The aim of the work reported in this article was to investigate the effects of medium temperature and industrial by-products on the key hardened properties of high performance concrete. Four concrete mixes were prepared based on a water-to-binder ratio of 0.35. Two industrial by-products, silica fume and Class F fly ash, were used separately and together with normal portland cement to produce three concrete mixes in addition to the control mix. The properties of both fresh and hardened concretes were examined in the laboratory. The freshly mixed concrete mixes were tested for slump, slump flow, and V-funnel flow. The hardened concretes were tested for compressive strength and dynamic modulus of elasticity after exposing to 20, 35 and 50 °C. In addition, the initial surface absorption and the rate of moisture movement into the concretes were determined at 20 °C. The performance of the concretes in the fresh state was excellent due to their superior deformability and good segregation resistance. In their hardened state, the highest levels of compressive strength and dynamic modulus of elasticity were produced by silica fume concrete. In addition, silica fume concrete showed the lowest level of initial surface absorption and the lowest rate of moisture movement into the interior of concrete. In comparison, the compressive strength, dynamic modulus of elasticity, initial surface absorption, and moisture movement rate of silica fume-fly ash concrete were close to those of silica fume concrete. Moreover, all concretes provided relatively low compressive strength and dynamic modulus of elasticity when they were exposed to 50 °C. However, the effect of increased temperature was less detrimental for silica fume and silica fume-fly ash concretes in comparison with the control concrete. PMID:28793732

Analytical and experimental evaluation of a 3-D hypersonic fixed-geometry, swept, mixed compression inlet

NASA Technical Reports Server (NTRS)

Agnone, Anthony M.

1987-01-01

The performance of a fixed-geometry, swept, mixed compression hypersonic inlet is presented. The experimental evaluation was conducted for a Mach number of 6.0 and for several angles of attack. The measured surface pressures and pitot pressure surveys at the inlet throat are compared to computations using a three-dimensional Euler code and an integral boundary layer theory. Unique features of the intake design, including the boundary layer control, insure a high inlet performance. The experimental data show the inlet has a high mass averaged total pressure recovery, a high mass capture and nearly uniform flow diffusion. The swept inlet exhibits excellent starting characteristics, and high flow stability at angle of attack.

When mountain belts disrupt mantle flow: from natural evidences to numerical modelling

NASA Astrophysics Data System (ADS)

Yamato, Philippe; Husson, Laurent; Guillaume, Benjamin

2016-04-01

During the Cenozoic, the number of orogens on Earth increased. This observation readily indicates that in the same time, compression in the lithosphere became gradually more and more important. Here, we show that such mountain belts, at plate boundaries, increasingly obstruct plate tectonics, slowing down and reorienting their motions. In turn, it changes the dynamic and kinematic surface conditions of the underlying flowing mantle, which ultimately modifies the pattern of mantle flow. Such forcing could explain many first order features of Cenozoic plate tectonics and mantle flow. Among others, at lithospheric scale, one can cite the compression of passive margins, the important variations in the rates of spreading at oceanic ridges, the initiation of subductions, or the onset of obductions. In the mantle, such changes in boundary conditions redesign the flow pattern and, consequently, disturb the oceanic lithosphere cooling. In order to test this hypothesis we first present thermo-mechanical numerical models of mantle convection above which a lithosphere is resting on top. Our results show that when collision occurs, the mantle flow is strongly modified, which leads to (i) increasing shear stresses below the lithosphere and (ii) a modification of the convection style. In turn, the transition between a "free" convection (mobile lid) and a "disturbed" convection (stagnant - or sluggish - lid) highly impacts the dynamics of the lithosphere at the surface. Thereby, on the basis of these models and a variety of real examples, we show that on the other side of a lithosphere presenting a collision zone, passive margins become squeezed and can undergo compression, which may ultimately evolve into subduction initiation or obduction. We also show that much further, due to the blocking of the lithosphere, spreading rates decrease at the ridge, which may explain a variety of features such as the low magmatism of ultraslow spreading ridges or the departure of slow spreading ridges from the half-space cooling model.

The many impacts of building mountain belts on plate tectonics and mantle flow

NASA Astrophysics Data System (ADS)

Yamato, Philippe; Husson, Laurent

2015-04-01

During the Cenozoic, the number of orogens on Earth increased. This observation readily indicates that in the same time, compression in the lithosphere became gradually more and more important. Such an increase of stresses in the lithosphere can impact on plate tectonics and mantle dynamics. We show that mountain belts at plate boundaries increasingly obstruct plate tectonics, slowing down and reorienting their motions. In turn, this changes the dynamic and kinematic surface conditions of the underlying flowing mantle. Ultimately, this modifies the pattern of mantle flow. This forcing could explain many first order features of Cenozoic plate tectonics and mantle flow. Among these, one can cite the compression of passive margins, the important variations in the rates of spreading at oceanic ridges, or the initiation of subduction, the onset of obduction, for the lithosphere. In the mantle, such change in boundary condition redesigns the pattern of mantle flow and, consequently, the oceanic lithosphere cooling. In order to test this hypothesis we first present thermo-mechanical numerical models of mantle convection above which a lithosphere rests. Our results show that when collision occurs, the mantle flow is highly modified, which leads to (i) increasing shear stresses below the lithosphere and (ii) to a modification of the convection style. In turn, the transition between a 'free' convection (mobile lid) and an 'upset' convection (stagnant -or sluggish- lid) highly impacts the dynamics of the lithosphere at the surface of the Earth. Thereby, on the basis of these models and a variety of real examples, we show that on the other side of a collision zone, passive margins become squeezed and can undergo compression, which may ultimately evolve into subduction or obduction. We also show that much further, due to the blocking of the lithosphere, spreading rates decrease at the ridge, a fact that may explain a variety of features such as the low magmatism of ultraslow spreading ridges or the departure of slow spreading ridges from the half-space cooling model.

A sharp interface method for compressible liquid–vapor flow with phase transition and surface tension

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fechter, Stefan, E-mail: stefan.fechter@iag.uni-stuttgart.de; Munz, Claus-Dieter, E-mail: munz@iag.uni-stuttgart.de; Rohde, Christian, E-mail: Christian.Rohde@mathematik.uni-stuttgart.de

The numerical approximation of non-isothermal liquid–vapor flow within the compressible regime is a difficult task because complex physical effects at the phase interfaces can govern the global flow behavior. We present a sharp interface approach which treats the interface as a shock-wave like discontinuity. Any mixing of fluid phases is avoided by using the flow solver in the bulk regions only, and a ghost-fluid approach close to the interface. The coupling states for the numerical solution in the bulk regions are determined by the solution of local two-phase Riemann problems across the interface. The Riemann solution accounts for the relevantmore » physics by enforcing appropriate jump conditions at the phase boundary. A wide variety of interface effects can be handled in a thermodynamically consistent way. This includes surface tension or mass/energy transfer by phase transition. Moreover, the local normal speed of the interface, which is needed to calculate the time evolution of the interface, is given by the Riemann solution. The interface tracking itself is based on a level-set method. The focus in this paper is the description of the two-phase Riemann solver and its usage within the sharp interface approach. One-dimensional problems are selected to validate the approach. Finally, the three-dimensional simulation of a wobbling droplet and a shock droplet interaction in two dimensions are shown. In both problems phase transition and surface tension determine the global bulk behavior.« less

Calculation of afterbody flows with a composite velocity formulation

NASA Technical Reports Server (NTRS)

Swanson, R. C.; Rubin, S. G.; Khosla, P. K.

1983-01-01

A recently developed technique for numerical solution of the Navier-Stokes equations for subsonic, laminar flows is investigated. It is extended here to allow for the computation of transonic and turbulent flows. The basic approach involves a multiplicative composite of the appropriate velocity representations for the inviscid and viscous flow regions. The resulting equations are structured so that far from the surface of the body the momentum equations lead to the Bernoulli equation for the pressure, while the continuity equation reduces to the familiar potential equation. Close to the body surface, the governing equations and solution techniques are characteristic of those describing interacting boundary layers. The velocity components are computed with a coupled strongly implicity procedure. For transonic flows the artificial compressibility method is used to treat supersonic regions. Calculations are made for both laminar and turbulent flows over axisymmetric afterbody configurations. Present results compare favorably with other numerical solutions and/or experimental data.

Analysis of unsteady compressible viscous layers

NASA Technical Reports Server (NTRS)

Power, G. D.; Verdon, J. M.; Kousen, K. A.

1990-01-01

The development of an analysis to predict the unsteady compressible flows in blade boundary layers and wakes is presented. The equations that govern the flows in these regions are transformed using an unsteady turbulent generalization of the Levy-Lees transformation. The transformed equations are solved using a finite difference technique in which the solution proceeds by marching in time and in the streamwise direction. Both laminar and turbulent flows are studied, the latter using algebraic turbulence and transition models. Laminar solutions for a flat plate are shown to approach classical asymptotic results for both high and low frequency unsteady motions. Turbulent flat-plate results are in qualitative agreement with previous predictions and measurements. Finally, the numerical technique is also applied to the stator and rotor of a low-speed turbine stage to determine unsteady effects on surface heating. The results compare reasonably well with measured heat transfer data and indicate that nonlinear effects have minimal impact on the mean and unsteady components of the flow.

Oscillatory Excitation of Unsteady Compressible Flows over Airfoils at Flight Reynolds Numbers

NASA Technical Reports Server (NTRS)

Seifert, Avi; Pack, LaTunia G.

1999-01-01

An experimental investigation, aimed at delaying flow separation due to the occurrence of a shock-wave-boundary-layer interaction, is reported. The experiment was performed using a NACA 0012 airfoil and a NACA 0015 airfoil at high Reynolds number incompressible and compressible flow conditions. The effects of Mach and Reynolds numbers were identified, using the capabilities of the cryogenic-pressurized facility to maintain one parameter fixed and change the other. Significant Reynolds number effects were identified in the baseline compressible flow conditions even at Reynolds number of 10 and 20 million. The main objectives of the experiment were to study the effects of periodic excitation on airfoil drag-divergence and to alleviate the severe unsteadiness associated with shock-induced separation (known as "buffeting"). Zero-mass-flux oscillatory blowing was introduced through a downstream directed slot located at 10% chord on the upper surface of the NACA 0015 airfoil. The effective frequencies generated 2-4 vortices over the separated region, regardless of the Mach number. Even though the excitation was introduced upstream of the shock-wave, due to experimental limitations, it had pronounced effects downstream of it. Wake deficit (associated with drag) and unsteadiness (associated with buffeting) were significantly reduced. The spectral content of the wake pressure fluctuations indicates of steadier flow throughout the frequency range when excitation was applied. This is especially important at low frequencies which are more likely to interact with the airframe.

Postseismic rebound in fault step-overs caused by pore fluid flow

USGS Publications Warehouse

Peltzer, G.; Rosen, P.; Rogez, F.; Hudnut, K.

1996-01-01

Near-field strain induced by large crustal earthquakes results in changes in pore fluid pressure that dissipate with time and produce surface deformation. Synthetic aperture radar (SAR) interferometry revealed several centimeters of postseismic uplift in pull-apart structures and subsidence in a compressive jog along the Landers, California, 1992 earthquake surface rupture, with a relaxation time of 270 ?? 45 days. Such a postseismic rebound may be explained by the transition of the Poisson's ratio of the deformed volumes of rock from undrained to drained conditions as pore fluid flow allows pore pressure to return to hydrostatic equilibrium.

Numerical modeling of Stokes flows over a superhydrophobic surface containing gas bubbles

NASA Astrophysics Data System (ADS)

Ageev, A. I.; Golubkina, I. V.; Osiptsov, A. N.

2017-10-01

This paper continues the numerical modeling of Stokes flows near cavities of a superhydrophobic surface, occupied by gas bubbles, based on the Boundary Element Method (BEM). The aim of the present study is to estimate the friction reduction (pressure drop) in a microchannel with a bottom superhydrophobic surface, the texture of which is formed by a periodic system of striped rectangular microcavities containing compressible gas bubbles. The model proposed takes into account the streamwise variation of the bubble shift into the cavities, caused by the longitudinal pressure gradient in the channel flow. The solution for the macroscopic (averaged) flow in the microchannel, constructed using an effective slip boundary condition on the superhydrophobic bottom wall, is matched with the solution of the Stokes problem at the microscale of a single cavity containing a gas bubble. The 2D Stokes problems of fluid flow over single cavities containing curved phase interfaces with the condition of zero shear stress are reduced to the boundary integral equations which are solved using the BEM method.

Structure of the Small Amplitude Motion on Transversely Sheared Mean Flows

NASA Technical Reports Server (NTRS)

Goldstein, Marvin E.; Afsar, Mohamed Z.; Leib, Stewart J.

2013-01-01

This paper considers the small amplitude unsteady motion of an inviscid non-heat conducting compressible fluid on a transversely sheared mean flow. It extends a previous result given in Goldstein (1978(b) and 1979(a)) which shows that the hydrodynamic component of the motion is determined by two arbitrary convected quantities in the absence of solid surfaces or other external sources. The result is important because it can be used to specify appropriate boundary conditions for unsteady surface interaction problems on transversely sheared mean flows in the same way that the vortical component of the Kovasznay (1953) decomposition is used to specify these conditions for surface interaction problems on uniform mean flows. But unlike the Kovasznay (1953) case the arbitrary convected quantities no longer bear a simple relation to the physical variables. One purpose of this paper is to derive a formula that relates these quantities to the (physically measurable) vorticity and pressure fluctuations in the flow.

Compressibility and compactibility of granules produced by wet and dry granulation.

PubMed

Bacher, C; Olsen, P M; Bertelsen, P; Sonnergaard, J M

2008-06-24

The bulk properties, compactibility and compressibility of granules produced by wet and dry granulation were compared applying a rotary tablet press, three different morphological forms of calcium carbonate and two particle sizes of sorbitol. Granules from both granulation methods possessed acceptable flow properties; however, the ground (Mikhart) and cubic (Scoralite) calcium carbonate demonstrated better die-filling abilities in the tablet press than the scalenhedral calcium carbonate (Sturcal). The wet processed granules showed in general larger compression properties. This was explained as these granules were mechanical stronger and had a higher initial porosity. In some cases, a large particle surface area of calcium carbonate and sorbitol resulted in a small, insignificant improvement of the consolidation characteristics. A correlation between the compression and compaction characteristics was demonstrated.

A Mass Tracking Formulation for Bubbles in Incompressible Flow

DTIC Science & Technology

2012-10-14

incompressible flow to fully nonlinear compressible flow including the effects of shocks and rarefactions , and then subsequently making a number of...using the ideas from [19] to couple together incompressible flow with fully nonlinear compressible flow including shocks and rarefactions . The results...compressible flow including the effects of shocks and rarefactions , and then subsequently making a number of simplifying assumptions on the air flow

On the finite length-scale of compressible shock-waves formed in free-surface flows of dry granular materials down a slope

NASA Astrophysics Data System (ADS)

Faug, Thierry

2017-04-01

The Rankine-Hugoniot jump conditions traditionally describe the theoretical relationship between the equilibrium state on both sides of a shock-wave. They are based on the crucial assumption that the length-scale needed to adjust the equilibrium state upstream of the shock to downstream of it is too small to be of significance to the problem. They are often used with success to describe the shock-waves in a number of applications found in both fluid and solid mechanics. However, the relations based on jump conditions at singular surfaces may fail to capture some features of the shock-waves formed in complex materials, such as granular matter. This study addresses the particular problem of compressible shock-waves formed in flows of dry granular materials down a slope. This problem is for instance relevant to full-scale geophysical granular flows in interaction with natural obstacles or man-made structures, such as topographical obstacles or mitigation dams respectively. Steady-state jumps formed in granular flows and travelling shock-waves produced at the impact of a granular avalanche-flow with a rigid wall are considered. For both situations, new analytical relations which do not consider that the granular shock-wave shrinks into a singular surface are derived, by using balance equations in their depth-averaged forms for mass and momentum. However, these relations need additional inputs that are closure relations for the size and the shape of the shock-wave, and a relevant constitutive friction law. Small-scale laboratory tests and numerical simulations based on the discrete element method are shortly presented and used to infer crucial information needed for the closure relations. This allows testing some predictive aspects of the simple analytical approach proposed for both steady-state and travelling shock-waves formed in free-surface flows of dry granular materials down a slope.

Experiments in a flighted conveyor comparing shear rates in compressed versus free surface flows

NASA Astrophysics Data System (ADS)

Pohlman, Nicholas; Higgins, Hannah; Krupiarz, Kamila; O'Connor, Ryan

2017-11-01

Uniformity of granular flow rate is critical in industry. Experiments in a flighted conveyor system aim to fill a gap in knowledge of achieving steady mass flow rate by correlating velocity profile data with mass flow rate measurements. High speed images were collected for uniformly-shaped particles in a bottom-driven flow conveyor belt system from which the velocity profiles can be generated. The correlation of mass flow rates from the velocity profiles to the time-dependent mass measurements will determine energy dissipation rates as a function of operating conditions. The velocity profiles as a function of the size of the particles, speed of the belt, and outlet size, will be compared to shear rate relationships found in past experiments that focused on gravity-driven systems. The dimension of the linear shear and type of decaying transition to the stationary bed may appear different due to the compression versus dilation space in open flows. The application of this research can serve to validate simulations in discrete element modeling and physically demonstrate a process that can be further developed and customized for industry applications, such as feeding a biomass conversion reactor. Sponsored by NIU's Office of Student Engagement and Experiential Learning.

CFD Simulations of the IHF Arc-Jet Flow: Compression-Pad Separation Bolt Wedge Tests

NASA Technical Reports Server (NTRS)

Gokcen, Tahir; Skokova, Kristina A.

2017-01-01

This paper reports computational analyses in support of two wedge tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using two different wedge models, each placed in a free jet downstream of a corresponding different conical nozzle in the Ames 60-MW Interaction Heating Facility. Each panel test article included a metallic separation bolt imbedded in Orion compression-pad and heatshield materials, resulting in a circular protuberance over a flat plate. The protuberances produce complex model flowfields, containing shock-shock and shock-boundary layer interactions, and multiple augmented heating regions on the test plate. As part of the test calibration runs, surface pressure and heat flux measurements on water-cooled calibration plates integrated with the wedge models were also obtained. Surface heating distributions on the test articles as well as arc-jet test environment parameters for each test configuration are obtained through computational fluid dynamics simulations, consistent with the facility and calibration measurements. The present analysis comprises simulations of the non-equilibrium flow field in the facility nozzle, test box, and flow field over test articles, and comparisons with the measured calibration data.

A Computational and Experimental Investigation of a Three-Dimensional Hypersonic Scramjet Inlet Flow Field. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Holland, Scott Douglas

1991-01-01

A combined computational and experimental parametric study of the internal aerodynamics of a generic three dimensional sidewall compression scramjet inlet configuration was performed. The study was designed to demonstrate the utility of computational fluid dynamics as a design tool in hypersonic inlet flow fields, to provide a detailed account of the nature and structure of the internal flow interactions, and to provide a comprehensive surface property and flow field database to determine the effects of contraction ratio, cowl position, and Reynolds number on the performance of a hypersonic scramjet inlet configuration.

Three-Dimensional Boundary Layers.

DTIC Science & Technology

1985-02-01

layer edge, We, is seen to increase fast in downstream direction. Near measuring station 9 the wall flow angle exceeds w = 55’, which means that the...leading edge along wing upper and lower surface to the trailing edge. As an excercise , such a boundary layer flow was computed for a simple symmetric...D.I.A. Poll The Development of Intermittent Turbulence on a Swept - Attachment Line Including the Effects of Compressibility. Aero. Qu. (Feb. 1983) 10

Three-dimensional rotational plasma flows near solid surfaces in an axial magnetic field

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gorshunov, N. M., E-mail: gorshunov-nm@nrcki.ru; Potanin, E. P., E-mail: potanin45@yandex.ru

2016-11-15

A rotational flow of a conducting viscous medium near an extended dielectric disk in a uniform axial magnetic field is analyzed in the magnetohydrodynamic (MHD) approach. An analytical solution to the system of nonlinear differential MHD equations of motion in the boundary layer for the general case of different rotation velocities of the disk and medium is obtained using a modified Slezkin–Targ method. A particular case of a medium rotating near a stationary disk imitating the end surface of a laboratory device is considered. The characteristics of a hydrodynamic flow near the disk surface are calculated within the model ofmore » a finite-thickness boundary layer. The influence of the magnetic field on the intensity of the secondary flow is studied. Calculations are performed for a weakly ionized dense plasma flow without allowance for the Hall effect and plasma compressibility. An MHD flow in a rotating cylinder bounded from above by a retarding cap is considered. The results obtained can be used to estimate the influence of the end surfaces on the main azimuthal flow, as well as the intensities of circulating flows in various devices with rotating plasmas, in particular, in plasma centrifuges and laboratory devices designed to study instabilities of rotating plasmas.« less

User Guide for Compressible Flow Toolbox Version 2.1 for Use With MATLAB(Registered Trademark); Version 7

NASA Technical Reports Server (NTRS)

Melcher, Kevin J.

2006-01-01

This report provides a user guide for the Compressible Flow Toolbox, a collection of algorithms that solve almost 300 linear and nonlinear classical compressible flow relations. The algorithms, implemented in the popular MATLAB programming language, are useful for analysis of one-dimensional steady flow with constant entropy, friction, heat transfer, or shock discontinuities. The solutions do not include any gas dissociative effects. The toolbox also contains functions for comparing and validating the equation-solving algorithms against solutions previously published in the open literature. The classical equations solved by the Compressible Flow Toolbox are: isentropic-flow equations, Fanno flow equations (pertaining to flow of an ideal gas in a pipe with friction), Rayleigh flow equations (pertaining to frictionless flow of an ideal gas, with heat transfer, in a pipe of constant cross section.), normal-shock equations, oblique-shock equations, and Prandtl-Meyer expansion equations. At the time this report was published, the Compressible Flow Toolbox was available without cost from the NASA Software Repository.

Numerical Study Comparing RANS and LES Approaches on a Circulation Control Airfoil

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Nishino, Takafumi

2011-01-01

A numerical study over a nominally two-dimensional circulation control airfoil is performed using a large-eddy simulation code and two Reynolds-averaged Navier-Stokes codes. Different Coanda jet blowing conditions are investigated. In addition to investigating the influence of grid density, a comparison is made between incompressible and compressible flow solvers. The incompressible equations are found to yield negligible differences from the compressible equations up to at least a jet exit Mach number of 0.64. The effects of different turbulence models are also studied. Models that do not account for streamline curvature effects tend to predict jet separation from the Coanda surface too late, and can produce non-physical solutions at high blowing rates. Three different turbulence models that account for streamline curvature are compared with each other and with large eddy simulation solutions. All three models are found to predict the Coanda jet separation location reasonably well, but one of the models predicts specific flow field details near the Coanda surface prior to separation much better than the other two. All Reynolds-averaged Navier-Stokes computations produce higher circulation than large eddy simulation computations, with different stagnation point location and greater flow acceleration around the nose onto the upper surface. The precise reasons for the higher circulation are not clear, although it is not solely a function of predicting the jet separation location correctly.

Fuel cell plates with skewed process channels for uniform distribution of stack compression load

DOEpatents

Granata, Jr., Samuel J.; Woodle, Boyd M.

1989-01-01

An electrochemical fuel cell includes an anode electrode, a cathode electrode, an electrolyte matrix sandwiched between electrodes, and a pair of plates above and below the electrodes. The plate above the electrodes has a lower surface with a first group of process gas flow channels formed thereon and the plate below the electrodes has an upper surface with a second group of process gas flow channels formed thereon. The channels of each group extend generally parallel to one another. The improvement comprises the process gas flow channels on the lower surface of the plate above the anode electrode and the process gas flow channels on the upper surface of the plate below the cathode electrode being skewed in opposite directions such that contact areas of the surfaces of the plates through the electrodes are formed in crisscross arrangements. Also, the plates have at least one groove in areas of the surfaces thereof where the channels are absent for holding process gas and increasing electrochemical activity of the fuel cell. The groove in each plate surface intersects with the process channels therein. Also, the opposite surfaces of a bipolar plate for a fuel cell contain first and second arrangements of process gas flow channels in the respective surfaces which are skewed the same amount in opposite directions relative to the longitudinal centerline of the plate.

Effect of compression pressure on inhalation grade lactose as carrier for dry powder inhalations

PubMed Central

Raut, Neha Sureshrao; Jamaiwar, Swapnil; Umekar, Milind Janrao; Kotagale, Nandkishor Ramdas

2016-01-01

Introduction: This study focused on the potential effects of compression forces experienced during lactose (InhaLac 70, 120, and 230) storage and transport on the flowability and aerosol performance in dry powder inhaler formulation. Materials and Methods: Lactose was subjected to typical compression forces 4, 10, and 20 N/cm2. Powder flowability and particle size distribution analysis of un-compressed and compressed lactose was evaluated by Carr's index, Hausner's ratio, the angle of repose and by laser diffraction method. Aerosol performance of un-compressed and compressed lactose was assessed in dispersion studies using glass twin-stage-liquid-impenger at flow rate 40-80 L/min. Results: At compression forces, the flowability of compressed lactose was observed same or slightly improved. Furthermore, compression of lactose caused a decrease in in vitro aerosol dispersion performance. Conclusion: The present study illustrates that, as carrier size increases, a concurrent decrease in drug aerosolization performance was observed. Thus, the compression of the lactose fines onto the surfaces of the larger lactose particles due to compression pressures was hypothesized to be the cause of these observed performance variations. The simulations of storage and transport in an industrial scale can induce significant variations in formulation performance, and it could be a source of batch-to-batch variations. PMID:27014618

Effect of double air injection on performance characteristics of centrifugal compressor

NASA Astrophysics Data System (ADS)

Hirano, Toshiyuki; Ogawa, Tatsuya; Yasui, Ryutaro; Tsujita, Hoshio

2017-02-01

In the operation of a centrifugal compressor of turbocharger, instability phenomena such as rotating stall and surge are induced at a lower flow rate close to the maximum pressure ratio. In this study, the compressed air at the exit of centrifugal compressor was re-circulated and injected to the impeller inlet by using two injection nozzles in order to suppress the surge phenomenon. The most effective circumferential position was examined to reduce the flow rate at the surge inception. Moreover, the influences of the injection on the fluctuating property of the flow field before and after the surge inception were investigated by examining the frequency of static pressure fluctuation on the wall surface and visualizing the compressor wall surface by oil-film visualization technique.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lu, Wei-Yang

Foam materials are used to protect sensitive components from impact loading. In order to predict and simulate the foam performance under various loading conditions, a validated foam model is needed and the mechanical properties of foams need to be characterized. Uniaxial compression and tension tests were conducted for different densities of foams under various temperatures and loading rates. Crush stress, tensile strength, and elastic modulus were obtained. A newly developed confined compression experiment provided data for investigating the foam flow direction. A biaxial tension experiment was also developed to explore the damage surface of a rigid polyurethane foam.

Mathematical modeling of the in-mold coating process for injection-molded thermoplastic parts

NASA Astrophysics Data System (ADS)

Chen, Xu

In-Mold Coating (IMC) has been successfully used for many years for exterior body panels made from compression molded Sheet Molding Compound (SMC). The coating material is a single component reactive fluid, designed to improve the surface quality of SMC moldings in terms of functional and cosmetic properties. When injected onto a cured SMC part, IMC cures and bonds to provide a pain-like surface. Because of its distinct advantages, IMC is being considered for application to injection molded thermoplastic parts. For a successful in mold coating operation, there are two key issues related to the flow of the coating. First, the injection nozzle should be located such that the thermoplastic substrate is totally covered and the potential for air trapping is minimized. The selected location should be cosmetically acceptable since it most likely will leave a mark on the coated surface. The nozzle location also needs to be accessible for easy of maintenance. Secondly, the hydraulic force generated by the coating injection pressure should not exceed the available clamping tonnage. If the clamping force is exceeded, coating leakage will occur. In this study, mathematical models for IMC flow on the compressible thermoplastic substrate have been developed. Finite Difference Method (FDM) is first used to solve the 1 dimensional (1D) IMC flow problem. In order to investigate the application of Control Volume based Finite Element Method (CV/FEM) to more complicated two dimensional IMC flow, that method is first evaluated by solving the 1D IMC flow problem. An analytical solution, which can be obtained when a linear relationship between the coating thickness and coating injection pressure is assumed, is used to verify the numerical results. The mathematical models for the 2 dimensional (2D) IMC flow are based on the generalized Hele-Shaw approximation. It has been found experimentally that the power law viscosity model adequately predicts the rheological behavior of the coating. The compressibility of the substrate is modeled by the 2-domain Tait PVT equation. CV/FEM is used to solve the discretized governing equations. A computer code has been developed to predict the fill pattern of the coating and the injection pressure. A number of experiments have been conducted to verify the numerical predictions of the computer code. It has been found both numerically and experimentally that the substrate thickness plays a significant role on the IMC fill pattern.

Compressible Flow Toolbox

NASA Technical Reports Server (NTRS)

Melcher, Kevin J.

2006-01-01

The Compressible Flow Toolbox is primarily a MATLAB-language implementation of a set of algorithms that solve approximately 280 linear and nonlinear classical equations for compressible flow. The toolbox is useful for analysis of one-dimensional steady flow with either constant entropy, friction, heat transfer, or Mach number greater than 1. The toolbox also contains algorithms for comparing and validating the equation-solving algorithms against solutions previously published in open literature. The classical equations solved by the Compressible Flow Toolbox are as follows: The isentropic-flow equations, The Fanno flow equations (pertaining to flow of an ideal gas in a pipe with friction), The Rayleigh flow equations (pertaining to frictionless flow of an ideal gas, with heat transfer, in a pipe of constant cross section), The normal-shock equations, The oblique-shock equations, and The expansion equations.

Multiparameter Analysis of Gas Transport Phenomena in Shale Gas Reservoirs: Apparent Permeability Characterization.

PubMed

Shen, Yinghao; Pang, Yu; Shen, Ziqi; Tian, Yuanyuan; Ge, Hongkui

2018-02-08

The large amount of nanoscale pores in shale results in the inability to apply Darcy's law. Moreover, the gas adsorption of shale increases the complexity of pore size characterization and thus decreases the accuracy of flow regime estimation. In this study, an apparent permeability model, which describes the adsorptive gas flow behavior in shale by considering the effects of gas adsorption, stress dependence, and non-Darcy flow, is proposed. The pore size distribution, methane adsorption capacity, pore compressibility, and matrix permeability of the Barnett and Eagle Ford shales are measured in the laboratory to determine the critical parameters of gas transport phenomena. The slip coefficients, tortuosity, and surface diffusivity are predicted via the regression analysis of the permeability data. The results indicate that the apparent permeability model, which considers second-order gas slippage, Knudsen diffusion, and surface diffusion, could describe the gas flow behavior in the transition flow regime for nanoporous shale. Second-order gas slippage and surface diffusion play key roles in the gas flow in nanopores for Knudsen numbers ranging from 0.18 to 0.5. Therefore, the gas adsorption and non-Darcy flow effects, which involve gas slippage, Knudsen diffusion, and surface diffusion, are indispensable parameters of the permeability model for shale.

High rate fabrication of compression molded components

DOE Office of Scientific and Technical Information (OSTI.GOV)

Matsen, Marc R.; Negley, Mark A.; Dykstra, William C.

2016-04-19

A method for fabricating a thermoplastic composite component comprises inductively heating a thermoplastic pre-form with a first induction coil by inducing current to flow in susceptor wires disposed throughout the pre-form, inductively heating smart susceptors in a molding tool to a leveling temperature with a second induction coil by applying a high-strength magnetic field having a magnetic flux that passes through surfaces of the smart susceptors, shaping the magnetic flux that passes through surfaces of the smart susceptors to flow substantially parallel to a molding surface of the smart susceptors, placing the heated pre-form between the heated smart susceptors; andmore » applying molding pressure to the pre-form to form the composite component.« less

Responses of three-dimensional flow to variations in the angle of incident wind and profile form of dunes: Greenwich Dunes, Prince Edward Island, Canada

NASA Astrophysics Data System (ADS)

Walker, Ian J.; Hesp, Patrick A.; Davidson-Arnott, Robin G. D.; Bauer, Bernard O.; Namikas, Steven L.; Ollerhead, Jeff

2009-04-01

This study reports the responses of three-dimensional near-surface airflow over a vegetated foredune to variations in the conditions of incident flow during an 8-h experiment. Two parallel measurement transects were established on morphologically different dune profiles: i) a taller, concave-convex West foredune transect with 0.5-m high, densely vegetated (45%), seaward incipient foredune, and ii) a shorter, concave-straight East foredune transect with lower, sparsely vegetated (14%) seaward incipient foredune. Five stations on each transect from the incipient dune to the crest were equipped with ultrasonic anemometers at 0.6 and 1.65 m height and logged at 1 Hz. Incident conditions were recorded from a 4-m tower over a flat beach. Winds increased from 6 m s - 1 to > 20 m s - 1 and were generally obliquely onshore (ENE, 73°). Three sub-events and the population of 10-minute averages of key properties of flow ( U, W, S, CV U) from all sample locations on the East transect ( n = 235) are examined to identify location- and profile-specific responses over 52° of the incident direction of flow (from 11 to 63° onshore). Topographic steering and forcing cause major deviations in the properties and vectors of near-surface flow from the regional wind. Topographic forcing on the concave-straight dune profile increases wind speed and steadiness toward the crest, with speed-up values to 65% in the backshore. Wind speed and steadiness of flow are least responsive to changes in incident angle in the backshore because of stagnation of flow and are most responsive at the lower stoss under pronounced streamline compression. On the steeper concave-convex profile, speed and steadiness decrease toward the crest because of stagnation of flow at the toe and flow expansion at the slope inflection point on the lower stoss. Net downward vertical velocity occurs over both profiles, increases toward the crest, and reflects enhanced turbulent momentum conveyance toward the surface. All of these flow responses are enhanced with faster speeds of incident flow and/or more onshore winds. Significant onshore steering of near-surface vectors of flow (to 37°) occurs and is greatest closer to the surface and during highly oblique winds (~ 15° onshore). Therefore, even subtle effects of streamline compression and amplification of flow under alongshore conditions effectively steer flow and sand transport toward the dune. As topographic forcing and steering cause significant, three-dimensional deviations in near-surface properties of flow, most regional-scale and/or two-dimensional models of dune process-response dynamics are insufficient for characterizing coastal and desert dune sediment budgets and morphodynamics. In particular, deflection of sand transport vectors with greater fetch distances than those derived from regional winds may occur. Coincident flow, transport and morphological response data are required to better quantitatively model these processes.

Recent developments in rotary-wing aerodynamic theory

NASA Technical Reports Server (NTRS)

Johnson, W.

1986-01-01

Current progress in the computational analysis of rotary-wing flowfields is surveyed, and some typical results are presented in graphs. Topics examined include potential theory, rotating coordinate systems, lifting-surface theory (moving singularity, fixed wing, and rotary wing), panel methods (surface singularity representations, integral equations, and compressible flows), transonic theory (the small-disturbance equation), wake analysis (hovering rotor-wake models and transonic blade-vortex interaction), limitations on computational aerodynamics, and viscous-flow methods (dynamic-stall theories and lifting-line theory). It is suggested that the present algorithms and advanced computers make it possible to begin working toward the ultimate goal of turbulent Navier-Stokes calculations for an entire rotorcraft.

A tabulation of pipe length to diameter ratios as a function of Mach number and pressure ratios for compressible flow

NASA Technical Reports Server (NTRS)

Dixon, G. V.; Barringer, S. R.; Gray, C. E.; Leatherman, A. D.

1975-01-01

Computer programs and resulting tabulations are presented of pipeline length-to-diameter ratios as a function of Mach number and pressure ratios for compressible flow. The tabulations are applicable to air, nitrogen, oxygen, and hydrogen for compressible isothermal flow with friction and compressible adiabatic flow with friction. Also included are equations for the determination of weight flow. The tabulations presented cover a wider range of Mach numbers for choked, adiabatic flow than available from commonly used engineering literature. Additional information presented, but which is not available from this literature, is unchoked, adiabatic flow over a wide range of Mach numbers, and choked and unchoked, isothermal flow for a wide range of Mach numbers.

General Equation Set Solver for Compressible and Incompressible Turbomachinery Flows

NASA Technical Reports Server (NTRS)

Sondak, Douglas L.; Dorney, Daniel J.

2002-01-01

Turbomachines for propulsion applications operate with many different working fluids and flow conditions. The flow may be incompressible, such as in the liquid hydrogen pump in a rocket engine, or supersonic, such as in the turbine which may drive the hydrogen pump. Separate codes have traditionally been used for incompressible and compressible flow solvers. The General Equation Set (GES) method can be used to solve both incompressible and compressible flows, and it is not restricted to perfect gases, as are many compressible-flow turbomachinery solvers. An unsteady GES turbomachinery flow solver has been developed and applied to both air and water flows through turbines. It has been shown to be an excellent alternative to maintaining two separate codes.

An immersed boundary formulation for simulating high-speed compressible viscous flows with moving solids

NASA Astrophysics Data System (ADS)

Qu, Yegao; Shi, Ruchao; Batra, Romesh C.

2018-02-01

We present a robust sharp-interface immersed boundary method for numerically studying high speed flows of compressible and viscous fluids interacting with arbitrarily shaped either stationary or moving rigid solids. The Navier-Stokes equations are discretized on a rectangular Cartesian grid based on a low-diffusion flux splitting method for inviscid fluxes and conservative high-order central-difference schemes for the viscous components. Discontinuities such as those introduced by shock waves and contact surfaces are captured by using a high-resolution weighted essentially non-oscillatory (WENO) scheme. Ghost cells in the vicinity of the fluid-solid interface are introduced to satisfy boundary conditions on the interface. Values of variables in the ghost cells are found by using a constrained moving least squares method (CMLS) that eliminates numerical instabilities encountered in the conventional MLS formulation. The solution of the fluid flow and the solid motion equations is advanced in time by using the third-order Runge-Kutta and the implicit Newmark integration schemes, respectively. The performance of the proposed method has been assessed by computing results for the following four problems: shock-boundary layer interaction, supersonic viscous flows past a rigid cylinder, moving piston in a shock tube and lifting off from a flat surface of circular, rectangular and elliptic cylinders triggered by shock waves, and comparing computed results with those available in the literature.

Direct simulation of isothermal-wall supersonic channel flow

NASA Technical Reports Server (NTRS)

Coleman, Gary N.

1993-01-01

The motivation for this work is the fact that in turbulent flows where compressibility effects are important, they are often poorly understood. A few examples of such flows are those associated with astrophysical phenomena and those found in combustion chambers, supersonic diffusers and nozzles, and over high-speed airfoils. For this project, we are primarily interested in compressibility effects near solid surfaces. Our main objective is an improved understanding of the fundamentals of compressible wall-bounded turbulence, which can in turn be used to cast light upon modeling concepts such as the Morkovin hypothesis and the Van Driest transformation. To this end, we have performed a direct numerical simulation (DNS) study of supersonic turbulent flow in a plane channel with constant-temperature walls. All of the relevant spatial and temporal scales are resolved so that no sub grid scale or turbulence model is necessary. The channel geometry was chosen so that finite Mach number effects can be isolated by comparing the present results to well established incompressible channel data. Here the fluid is assumed to be an ideal gas with constant specific heats, constant Prandtl number, and power-law temperature-dependent viscosity. Isothermal-wall boundary conditions are imposed so that a statistically stationary state may be obtained. The flow is driven by a uniform (in space) body force (rather than a mean pressure gradient) to preserve stream wise homogeneity, with the body force defined so that the total mass flux is constant.

Detailed computational procedure for design of cascade blades with prescribed velocity distributions in compressible potential flows

NASA Technical Reports Server (NTRS)

Costello, George R; Cummings, Robert L; Sinnette, John T , Jr

1952-01-01

A detailed step-by-step computational outline is presented for the design of two-dimensional cascade blades having a prescribed velocity distribution on the blade in a potential flow of the usual compressible fluid. The outline is based on the assumption that the magnitude of the velocity in the flow of the usual compressible nonviscous fluid is proportional to the magnitude of the velocity in the flow of a compressible nonviscous fluid with linear pressure-volume relation.

Entrainment of bed material by Earth-surface mass flows: review and reformulation of depth-integrated theory

USGS Publications Warehouse

Iverson, Richard M.; Chaojun Ouyang,

2015-01-01

Earth-surface mass flows such as debris flows, rock avalanches, and dam-break floods can grow greatly in size and destructive potential by entraining bed material they encounter. Increasing use of depth-integrated mass- and momentum-conservation equations to model these erosive flows motivates a review of the underlying theory. Our review indicates that many existing models apply depth-integrated conservation principles incorrectly, leading to spurious inferences about the role of mass and momentum exchanges at flow-bed boundaries. Model discrepancies can be rectified by analyzing conservation of mass and momentum in a two-layer system consisting of a moving upper layer and static lower layer. Our analysis shows that erosion or deposition rates at the interface between layers must in general satisfy three jump conditions. These conditions impose constraints on valid erosion formulas, and they help determine the correct forms of depth-integrated conservation equations. Two of the three jump conditions are closely analogous to Rankine-Hugoniot conditions that describe the behavior of shocks in compressible gasses, and the third jump condition describes shear traction discontinuities that necessarily exist across eroding boundaries. Grain-fluid mixtures commonly behave as compressible materials as they undergo entrainment, because changes in bulk density occur as the mixtures mobilize and merge with an overriding flow. If no bulk density change occurs, then only the shear-traction jump condition applies. Even for this special case, however, accurate formulation of depth-integrated momentum equations requires a clear distinction between boundary shear tractions that exist in the presence or absence of bed erosion.

Application of PDF methods to compressible turbulent flows

NASA Astrophysics Data System (ADS)

Delarue, B. J.; Pope, S. B.

1997-09-01

A particle method applying the probability density function (PDF) approach to turbulent compressible flows is presented. The method is applied to several turbulent flows, including the compressible mixing layer, and good agreement is obtained with experimental data. The PDF equation is solved using a Lagrangian/Monte Carlo method. To accurately account for the effects of compressibility on the flow, the velocity PDF formulation is extended to include thermodynamic variables such as the pressure and the internal energy. The mean pressure, the determination of which has been the object of active research over the last few years, is obtained directly from the particle properties. It is therefore not necessary to link the PDF solver with a finite-volume type solver. The stochastic differential equations (SDE) which model the evolution of particle properties are based on existing second-order closures for compressible turbulence, limited in application to low turbulent Mach number flows. Tests are conducted in decaying isotropic turbulence to compare the performances of the PDF method with the Reynolds-stress closures from which it is derived, and in homogeneous shear flows, at which stage comparison with direct numerical simulation (DNS) data is conducted. The model is then applied to the plane compressible mixing layer, reproducing the well-known decrease in the spreading rate with increasing compressibility. It must be emphasized that the goal of this paper is not as much to assess the performance of models of compressibility effects, as it is to present an innovative and consistent PDF formulation designed for turbulent inhomogeneous compressible flows, with the aim of extending it further to deal with supersonic reacting flows.

Effect of deformation on the thermal conductivity of granular porous media with rough grain surface

NASA Astrophysics Data System (ADS)

Askari, Roohollah; Hejazi, S. Hossein; Sahimi, Muhammad

2017-08-01

Heat transfer in granular porous media is an important phenomenon that is relevant to a wide variety of problems, including geothermal reservoirs and enhanced oil recovery by thermal methods. Resistance to flow of heat in the contact area between the grains strongly influences the effective thermal conductivity of such porous media. Extensive experiments have indicated that the roughness of the grains' surface follows self-affine fractal stochastic functions, and thus, the contact resistance cannot be accounted for by models based on smooth surfaces. Despite the significance of rough contact area, the resistance has been accounted for by a fitting parameter in the models of heat transfer. In this Letter we report on a study of conduction in a packing of particles that contains a fluid of a given conductivity, with each grain having a rough self-affine surface, and is under an external compressive pressure. The deformation of the contact area depends on the fractal dimension that characterizes the grains' rough surface, as well as their Young's modulus. Excellent qualitative agreement is obtained with experimental data. Deformation of granular porous media with grains that have rough self-affine fractal surface is simulated. Thermal contact resistance between grains with rough surfaces is incorporated into the numerical simulation of heat conduction under compressive pressure. By increasing compressive pressure, thermal conductivity is enhanced more in the grains with smoother surfaces and lower Young's modulus. Excellent qualitative agreement is obtained with the experimental data.

A compressibility correction of the pressure strain correlation model in turbulent flow

NASA Astrophysics Data System (ADS)

Klifi, Hechmi; Lili, Taieb

2013-07-01

This paper is devoted to the second-order closure for compressible turbulent flows with special attention paid to modeling the pressure-strain correlation appearing in the Reynolds stress equation. This term appears as the main one responsible for the changes of the turbulence structures that arise from structural compressibility effects. From the analysis and DNS results of Simone et al. and Sarkar, the compressibility effects on the homogeneous turbulence shear flow are parameterized by the gradient Mach number. Several experiment and DNS results suggest that the convective Mach number is appropriate to study the compressibility effects on the mixing layers. The extension of the LRR model recently proposed by Marzougui, Khlifi and Lili for the pressure-strain correlation gives results that are in disagreement with the DNS results of Sarkar for high-speed shear flows. This extension is revised to derive a turbulence model for the pressure-strain correlation in which the compressibility is included in the turbulent Mach number, the gradient Mach number and then the convective Mach number. The behavior of the proposed model is compared to the compressible model of Adumitroiae et al. for the pressure-strain correlation in two turbulent compressible flows: homogeneous shear flow and mixing layers. In compressible homogeneous shear flows, the predicted results are compared with the DNS data of Simone et al. and those of Sarkar. For low compressibility, the two compressible models are similar, but they become substantially different at high compressibilities. The proposed model shows good agreement with all cases of DNS results. Those of Adumitroiae et al. do not reflect any effect of a change in the initial value of the gradient Mach number on the Reynolds stress anisotropy. The models are used to simulate compressible mixing layers. Comparison of our predictions with those of Adumitroiae et al. and with the experimental results of Goebel et al. shows good qualitative agreement.

Flow/Damage Surfaces for Fiber-Reinforced Metals Having Different Periodic Microstructures

NASA Technical Reports Server (NTRS)

Lissenden, Cliff J.; Arnold, Steven M.; Iyer, Saiganesh K.

1998-01-01

Flow/damage surfaces can be defined in terms of stress, inelastic strain rate, and internal variables using a thermodynamics framework. A macroscale definition relevant to thermodynamics and usable in an experimental program is employed to map out surfaces of constant inelastic power in various stress planes. The inelastic flow of a model silicon carbide/ titanium composite system having rectangular, hexagonal, and square diagonal fiber packing arrays subjected to biaxial stresses is quantified by flow/damage surfaces that are determined numerically from micromechanics, using both finite element analysis and the generalized method of cells. Residual stresses from processing are explicitly included and damage in the form of fiber-matrix debonding under transverse tensile and/or shear loading is represented by a simple interface model. The influence of microstructural architecture is largest whenever fiber-matrix debonding is not an issue; for example in the presence of transverse compressive stresses. Additionally, as the fiber volume fraction increases, so does the effect of microstructural architecture. With regard to the micromechanics analysis, the overall inelastic flow predicted by the generalized method of cells is in excellent agreement with that predicted using a large number of displacement-based finite elements.

Flow/Damage Surfaces for Fiber-Reinforced Metals having Different Periodic Microstructures

NASA Technical Reports Server (NTRS)

Lissenden, Cliff J.; Arnold, Steven M.; Iyer, Saiganesh K.

1998-01-01

Flow/damage surfaces can be defined in terms of stress, inelastic strain rate, and internal variables using a thermodynamics framework. A macroscale definition relevant to thermodynamics and usable in an experimental program is employed to map out surfaces of constant inelastic power in various stress planes. The inelastic flow of a model silicon carbide/ titanium composite system having rectangular, hexagonal, and square diagonal fiber packing, arrays subjected to biaxial stresses is quantified by flow/damage surfaces that are determined numerically from micromechanics. using both finite element analysis and the generalized method of cells. Residual stresses from processing are explicitly included and damage in the form of fiber-matrix debonding under transverse tensile and/or shear loading is represented by a simple interface model. The influence of microstructural architecture is largest whenever fiber-matrix debonding is not an issue, for example in the presence of transverse compressive stresses. Additionally, as the fiber volume fraction increases, so does the effect of microstructural architecture. With regard to the micromechanics analysis, the overall inelastic flow predicted by the generalized method of cells is in excellent agreement with that predicted using a large number of displacement-based finite elements.

On the Importance of Spatial Resolution for Flap Side Edge Noise Prediction

NASA Technical Reports Server (NTRS)

Mineck, Raymond E.; Khorrami, Mehdi R.

2017-01-01

A spatial resolution study of flap tip flow and the effects on the farfield noise signature for an 18%-scale, semispan Gulfstream aircraft model are presented. The NASA FUN3D unstructured, compressible Navier-Stokes solver was used to perform the highly resolved, time-dependent, detached eddy simulations of the flow field associated with the flap for this high-fidelity aircraft model. Following our previous work on the same model, the latest computations were undertaken to determine the causes of deficiencies observed in our earlier predictions of the steady and unsteady surface pressures and off-surface flow field at the flap tip regions, in particular the outboard tip area, where the presence of a cavity at the side-edge produces very complex flow features and interactions. The present results show gradual improvement in steady loading at the outboard flap edge region with increasing spatial resolution, yielding more accurate fluctuating surface pressures, off-surface flow field, and farfield noise with improved high-frequency content when compared with wind tunnel measurements. The spatial resolution trends observed in the present study demonstrate that the deficiencies reported in our previous computations are mostly caused by inadequate spatial resolution and are not related to the turbulence model.

FUEL ROD CLUSTERS

DOEpatents

Schultz, A.B.

1959-08-01

A cluster of nuclear fuel rods and a tubular casing therefor through which a coolant flows in heat-exchange contact with the fuel rods is described. The fuel rcds are held in the casing by virtue of the compressive force exerted between longitudinal ribs of the fuel rcds and internal ribs of the casing or the internal surfaces thereof.

SUPIN: A Computational Tool for Supersonic Inlet Design

NASA Technical Reports Server (NTRS)

Slater, John W.

2016-01-01

A computational tool named SUPIN is being developed to design and analyze the aerodynamic performance of supersonic inlets. The inlet types available include the axisymmetric pitot, three-dimensional pitot, axisymmetric outward-turning, two-dimensional single-duct, two-dimensional bifurcated-duct, and streamline-traced inlets. The aerodynamic performance is characterized by the flow rates, total pressure recovery, and drag. The inlet flow-field is divided into parts to provide a framework for the geometry and aerodynamic modeling. Each part of the inlet is defined in terms of geometric factors. The low-fidelity aerodynamic analysis and design methods are based on analytic, empirical, and numerical methods which provide for quick design and analysis. SUPIN provides inlet geometry in the form of coordinates, surface angles, and cross-sectional areas. SUPIN can generate inlet surface grids and three-dimensional, structured volume grids for use with higher-fidelity computational fluid dynamics (CFD) analysis. Capabilities highlighted in this paper include the design and analysis of streamline-traced external-compression inlets, modeling of porous bleed, and the design and analysis of mixed-compression inlets. CFD analyses are used to verify the SUPIN results.

DSMC computations of hypersonic flow separation and re-attachment in the transition to continuum regime

NASA Astrophysics Data System (ADS)

Prakash, Ram; Gai, Sudhir L.; O'Byrne, Sean; Brown, Melrose

2016-11-01

The flow over a `tick' shaped configuration is performed using two Direct Simulation Monte Carlo codes: the DS2V code of Bird and the code from Sandia National Laboratory, called SPARTA. The configuration creates a flow field, where the flow is expanded initially but then is affected by the adverse pressure gradient induced by a compression surface. The flow field is challenging in the sense that the full flow domain is comprised of localized areas spanning continuum and transitional regimes. The present work focuses on the capability of SPARTA to model such flow conditions and also towards a comparative evaluation with results from DS2V. An extensive grid adaptation study is performed using both the codes on a model with a sharp leading edge and the converged results are then compared. The computational predictions are evaluated in terms of surface parameters such as heat flux, shear stress, pressure and velocity slip. SPARTA consistently predicts higher values for these surface properties. The skin friction predictions of both the codes don't give any indication of separation but the velocity slip plots indicate an incipient separation behavior at the corner. The differences in the results are attributed towards the flow resolution at the leading edge that dictates the downstream flow characteristics.

Unified approach for incompressible flows

NASA Astrophysics Data System (ADS)

Chang, Tyne-Hsien

1993-12-01

An unified approach for solving both compressible and incompressible flows was investigated in this study. The difference in CFD code development between incompressible and compressible flows is due to the mathematical characteristics. However, if one can modify the continuity equation for incompressible flows by introducing pseudocompressibility, the governing equations for incompressible flows would have the same mathematical characters as compressible flows. The application of a compressible flow code to solve incompressible flows becomes feasible. Among numerical algorithms developed for compressible flows, the Centered Total Variation Diminishing (CTVD) schemes possess better mathematical properties to damp out the spurious oscillations while providing high-order accuracy for high speed flows. It leads us to believe that CTVD schemes can equally well solve incompressible flows. In this study, the governing equations for incompressible flows include the continuity equation and momentum equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be implemented. In addition, the boundary conditions including physical and numerical boundary conditions must be properly specified to obtain accurate solution. The CFD code for this research is currently in progress. Flow past a circular cylinder will be used for numerical experiments to determine the accuracy and efficiency of the code before applying this code to more specific applications.

Lagrangian particle method for compressible fluid dynamics

NASA Astrophysics Data System (ADS)

Samulyak, Roman; Wang, Xingyu; Chen, Hsin-Chiang

2018-06-01

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface/multiphase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremal points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free interfaces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order. The method is generalizable to coupled hyperbolic-elliptic systems. Numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.

A low diffusive Lagrange-remap scheme for the simulation of violent air-water free-surface flows

NASA Astrophysics Data System (ADS)

Bernard-Champmartin, Aude; De Vuyst, Florian

2014-10-01

In 2002, Després and Lagoutière [17] proposed a low-diffusive advection scheme for pure transport equation problems, which is particularly accurate for step-shaped solutions, and thus suited for interface tracking procedure by a color function. This has been extended by Kokh and Lagoutière [28] in the context of compressible multifluid flows using a five-equation model. In this paper, we explore a simplified variant approach for gas-liquid three-equation models. The Eulerian numerical scheme has two ingredients: a robust remapped Lagrange solver for the solution of the volume-averaged equations, and a low diffusive compressive scheme for the advection of the gas mass fraction. Numerical experiments show the performance of the computational approach on various flow reference problems: dam break, sloshing of a tank filled with water, water-water impact and finally a case of Rayleigh-Taylor instability. One of the advantages of the present interface capturing solver is its natural implementation on parallel processors or computers.

Evaluation of a research circulation control airfoil using Navier-Stokes methods

NASA Technical Reports Server (NTRS)

Shrewsbury, George D.

1987-01-01

The compressible Reynolds time averaged Navier-Stokes equations were used to obtain solutions for flows about a two dimensional circulation control airfoil. The governing equations were written in conservation form for a body-fitted coordinate system and solved using an Alternating Direction Implicit (ADI) procedure. A modified algebraic eddy viscosity model was used to define the turbulent characteristics of the flow, including the wall jet flow over the Coanda surface at the trailing edge. Numerical results are compared to experimental data obtained for a research circulation control airfoil geometry. Excellent agreement with the experimental results was obtained.

Linearized compressible-flow theory for sonic flight speeds

NASA Technical Reports Server (NTRS)

Heaslet, Max A; Lomax, Harvard; Spreiter, John R

1950-01-01

The partial differential equation for the perturbation velocity potential is examined for free-stream Mach numbers close to and equal to one. It is found that, under the assumptions of linearized theory, solutions can be found consistent with the theory for lifting-surface problems both in stationary three-dimensional flow and in unsteady two-dimensional flow. Several examples are solved including a three dimensional swept-back wing and two dimensional harmonically-oscillating wing, both for a free stream Mach number equal to one. Momentum relations for the evaluation of wave and vortex drag are also discussed. (author)

On the Examination of Darcy Permeability a Thin Fibrous Porous Layer

NASA Astrophysics Data System (ADS)

Zhu, Zenghao; Wang, Qiuyun; Wu, Qianhong; Vucbmss Team

2016-11-01

In this paper, we report a novel experimental approach to investigate the Darcy permeability of a soft and thin fibrous porous layer. The project is inspired by recent studies involved compression of very thin porous films and the resultant pore fluid flow inside the confined porous structure. The Darcy permeability plays a critical role during the process, which however, is tricky to measure due to the very thin nature of the porous media. In the current study, a special micro-fluidic device is developed that consists of a rectangular flow channel with adjustable gap height ranging from 20 mm to 0.5 mm. Air is forced through the thin gap filled with testing fibrous materials. By measuring the flow rate and the pressure drop, we have successfully obtained the Darcy permeability of different thin porous sheets at different compression ratios. Furthermore, the surface area of the fibers are evaluated using a Micromeritics® ASAP 2020 (Accelerated Surface Area and Porosimetry) system. We found that, although the functions relating the permeability and porosities are different for different fibrous materials, these functions collapse to a single relationship if one express the permeability as a function of the solid phase surface area per unit volume. This finding provides a useful approach to evaluate the permeability of very thin fibrous porous sheet, which otherwise is difficult to measure directly. This research was supported by the National Science Foundation under Award #1511096.

Crystal and Particle Engineering Strategies for Improving Powder Compression and Flow Properties to Enable Continuous Tablet Manufacturing by Direct Compression.

PubMed

Chattoraj, Sayantan; Sun, Changquan Calvin

2018-04-01

Continuous manufacturing of tablets has many advantages, including batch size flexibility, demand-adaptive scale up or scale down, consistent product quality, small operational foot print, and increased manufacturing efficiency. Simplicity makes direct compression the most suitable process for continuous tablet manufacturing. However, deficiencies in powder flow and compression of active pharmaceutical ingredients (APIs) limit the range of drug loading that can routinely be considered for direct compression. For the widespread adoption of continuous direct compression, effective API engineering strategies to address power flow and compression problems are needed. Appropriate implementation of these strategies would facilitate the design of high-quality robust drug products, as stipulated by the Quality-by-Design framework. Here, several crystal and particle engineering strategies for improving powder flow and compression properties are summarized. The focus is on the underlying materials science, which is the foundation for effective API engineering to enable successful continuous manufacturing by the direct compression process. Copyright © 2018 American Pharmacists Association®. Published by Elsevier Inc. All rights reserved.

A workstation based simulator for teaching compressible aerodynamics

NASA Technical Reports Server (NTRS)

Benson, Thomas J.

1994-01-01

A workstation-based interactive flow simulator has been developed to aid in the teaching of undergraduate compressible aerodynamics. By solving the equations found in NACA 1135, the simulator models three basic fluids problems encountered in supersonic flow: flow past a compression corner, flow past two wedges in series, and flow past two opposed wedges. The study can vary the geometry or flow conditions through a graphical user interface and the new conditions are calculated immediately. Various graphical formats present the results of the flow calculations to the student. The simulator includes interactive questions and answers to aid in both the use of the tool and to develop an understanding of some of the complexities of compressible aerodynamics. A series of help screens make the simulator easy to learn and use.

MR Elastography Can Be Used to Measure Brain Stiffness Changes as a Result of Altered Cranial Venous Drainage During Jugular Compression.

PubMed

Hatt, A; Cheng, S; Tan, K; Sinkus, R; Bilston, L E

2015-10-01

Compressing the internal jugular veins can reverse ventriculomegaly in the syndrome of inappropriately low pressure acute hydrocephalus, and it has been suggested that this works by "stiffening" the brain tissue. Jugular compression may also alter blood and CSF flow in other conditions. We aimed to understand the effect of jugular compression on brain tissue stiffness and CSF flow. The head and neck of 9 healthy volunteers were studied with and without jugular compression. Brain stiffness (shear modulus) was measured by using MR elastography. Phase-contrast MR imaging was used to measure CSF flow in the cerebral aqueduct and blood flow in the neck. The shear moduli of the brain tissue increased with the percentage of blood draining through the internal jugular veins during venous compression. Peak velocity of caudally directed CSF in the aqueduct increased significantly with jugular compression (P < .001). The mean jugular venous flow rate, amplitude, and vessel area were significantly reduced with jugular compression, while cranial arterial flow parameters were unaffected. Jugular compression influences cerebral CSF hydrodynamics in healthy subjects and can increase brain tissue stiffness, but the magnitude of the stiffening depends on the percentage of cranial blood draining through the internal jugular veins during compression—that is, subjects who maintain venous drainage through the internal jugular veins during jugular compression have stiffer brains than those who divert venous blood through alternative pathways. These methods may be useful for studying this phenomenon in patients with the syndrome of inappropriately low-pressure acute hydrocephalus and other conditions. © 2015 by American Journal of Neuroradiology.

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

NASA Technical Reports Server (NTRS)

Chung, T. J. (Editor); Karr, Gerald R. (Editor)

1989-01-01

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Integrated turbomachine oxygen plant

DOEpatents

Anand, Ashok Kumar; DePuy, Richard Anthony; Muthaiah, Veerappan

2014-06-17

An integrated turbomachine oxygen plant includes a turbomachine and an air separation unit. One or more compressor pathways flow compressed air from a compressor through one or more of a combustor and a turbine expander to cool the combustor and/or the turbine expander. An air separation unit is operably connected to the one or more compressor pathways and is configured to separate the compressed air into oxygen and oxygen-depleted air. A method of air separation in an integrated turbomachine oxygen plant includes compressing a flow of air in a compressor of a turbomachine. The compressed flow of air is flowed through one or more of a combustor and a turbine expander of the turbomachine to cool the combustor and/or the turbine expander. The compressed flow of air is directed to an air separation unit and is separated into oxygen and oxygen-depleted air.

CNSFV code development, virtual zone Navier-Stokes computations of oscillating control surfaces and computational support of the laminar flow supersonic wind tunnel

NASA Technical Reports Server (NTRS)

Klopfer, Goetz H.

1993-01-01

The work performed during the past year on this cooperative agreement covered two major areas and two lesser ones. The two major items included further development and validation of the Compressible Navier-Stokes Finite Volume (CNSFV) code and providing computational support for the Laminar Flow Supersonic Wind Tunnel (LFSWT). The two lesser items involve a Navier-Stokes simulation of an oscillating control surface at transonic speeds and improving the basic algorithm used in the CNSFV code for faster convergence rates and more robustness. The work done in all four areas is in support of the High Speed Research Program at NASA Ames Research Center.

An immersed boundary-simplified sphere function-based gas kinetic scheme for simulation of 3D incompressible flows

NASA Astrophysics Data System (ADS)

Yang, L. M.; Shu, C.; Yang, W. M.; Wang, Y.; Wu, J.

2017-08-01

In this work, an immersed boundary-simplified sphere function-based gas kinetic scheme (SGKS) is presented for the simulation of 3D incompressible flows with curved and moving boundaries. At first, the SGKS [Yang et al., "A three-dimensional explicit sphere function-based gas-kinetic flux solver for simulation of inviscid compressible flows," J. Comput. Phys. 295, 322 (2015) and Yang et al., "Development of discrete gas kinetic scheme for simulation of 3D viscous incompressible and compressible flows," J. Comput. Phys. 319, 129 (2016)], which is often applied for the simulation of compressible flows, is simplified to improve the computational efficiency for the simulation of incompressible flows. In the original SGKS, the integral domain along the spherical surface for computing conservative variables and numerical fluxes is usually not symmetric at the cell interface. This leads the expression of numerical fluxes at the cell interface to be relatively complicated. For incompressible flows, the sphere at the cell interface can be approximately considered to be symmetric as shown in this work. Besides that, the energy equation is usually not needed for the simulation of incompressible isothermal flows. With all these simplifications, the simple and explicit formulations for the conservative variables and numerical fluxes at the cell interface can be obtained. Second, to effectively implement the no-slip boundary condition for fluid flow problems with complex geometry as well as moving boundary, the implicit boundary condition-enforced immersed boundary method [Wu and Shu, "Implicit velocity correction-based immersed boundary-lattice Boltzmann method and its applications," J. Comput. Phys. 228, 1963 (2009)] is introduced into the simplified SGKS. That is, the flow field is solved by the simplified SGKS without considering the presence of an immersed body and the no-slip boundary condition is implemented by the immersed boundary method. The accuracy and efficiency of the present scheme are validated by simulating the decaying vortex flow, flow past a stationary and rotating sphere, flow past a stationary torus, and flows over dragonfly flight.

Near-wall modeling of compressible turbulent flow

NASA Technical Reports Server (NTRS)

So, Ronald M. C.

1991-01-01

A near-wall two-equation model for compressible flows is proposed. The model is formulated by relaxing the assumption of dynamic field similarity between compressible and incompressible flows. A postulate is made to justify the extension of incompressible models to ammount for compressibility effects. This requires formulation the turbulent kinetic energy equation in a form similar to its incompressible counterpart. As a result, the compressible dissipation function has to be split into a solenoidal part, which is not sensitive to changes of compressibility indicators, and a dilatational part, which is directly affected by these changes. A model with an explicit dependence on the turbulent Mach number is proposed for the dilatational dissipation rate.

Coherent vertical structures in numerical simulations of buoyant plumes from wildland fires

Treesearch

Philip Cunningham; Scott L. Goodrick; M. Yousuff Hussaini; Rodman R. Linn

2005-01-01

The structure and dynamics of buoyant plumes arising from surface-based heat sources in a vertically sheared ambient atmospheric flow are examined via simulations of a three-dimensional, compressible numerical model. Simple circular heat sources and asymmetric elliptical ring heat sources that are representative of wildland fires of moderate intensity are considered....

Active bypass flow control for a seal in a gas turbine engine

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ebert, Todd A.; Kimmel, Keith D.

An active bypass flow control system for controlling bypass compressed air based upon leakage flow of compressed air flowing past an outer balance seal between a stator and rotor of a first stage of a gas turbine in a gas turbine engine is disclosed. The active bypass flow control system is an adjustable system in which one or more metering devices may be used to control the flow of bypass compressed air as the flow of compressed air past the outer balance seal changes over time as the outer balance seal between the rim cavity and the cooling cavity wears.more » In at least one embodiment, the metering device may include a valve formed from one or more pins movable between open and closed positions in which the one pin at least partially bisects the bypass channel to regulate flow.« less

The least-squares finite element method for low-mach-number compressible viscous flows

NASA Technical Reports Server (NTRS)

Yu, Sheng-Tao

1994-01-01

The present paper reports the development of the Least-Squares Finite Element Method (LSFEM) for simulating compressible viscous flows at low Mach numbers in which the incompressible flows pose as an extreme. Conventional approach requires special treatments for low-speed flows calculations: finite difference and finite volume methods are based on the use of the staggered grid or the preconditioning technique; and, finite element methods rely on the mixed method and the operator-splitting method. In this paper, however, we show that such difficulty does not exist for the LSFEM and no special treatment is needed. The LSFEM always leads to a symmetric, positive-definite matrix through which the compressible flow equations can be effectively solved. Two numerical examples are included to demonstrate the method: first, driven cavity flows at various Reynolds numbers; and, buoyancy-driven flows with significant density variation. Both examples are calculated by using full compressible flow equations.

Dynamic alteration of regional cerebral blood flow during carotid compression and proof of reversibility.

PubMed

Asahi, Kouichi; Hori, M; Hamasaki, N; Sato, S; Nakanishi, H; Kuwatsuru, R; Sasai, K; Aoki, S

2012-01-01

It is difficult to non-invasively visualize changes in regional cerebral blood flow caused by manual compression of the carotid artery. To visualize dynamic changes in regional cerebral blood flow during and after manual compression of the carotid artery. Two healthy volunteers were recruited. Anatomic features and flow directions in the circle of Willis were evaluated with time-of-flight magnetic resonance angiography (MRA) and two-dimensional phase-contrast (2DPC) MRA, respectively. Regional cerebral blood flow was visualized with territorial arterial spin-labeling magnetic resonance imaging (TASL-MRI). TASL-MRI and 2DPC-MRA were performed in three states: at rest, during manual compression of the right carotid artery, and after decompression. In one volunteer, time-space labeling inversion pulse (Time-SLIP) MRA was performed to confirm collateral flow. During manual carotid compression, in one volunteer, the right thalamus changed to be fed only by the vertebrobasilar system, and the right basal ganglia changed to be fed by the left internal carotid artery. In the other volunteer, the right basal ganglia changed to be fed by the vertebrobasilar system. 2DPC-MRA showed that the flow direction changed in the right A1 segment of the anterior cerebral artery and the right posterior communicating artery. Perfusion patterns and flow directions recovered after decompression. Time-SLIP MRA showed pial vessels and dural collateral circulation when the right carotid artery was manually compressed. Use of TASL-MRI and 2DPC-MRA was successful for non-invasive visualization of the dynamic changes in regional cerebral blood flow during and after manual carotid compression.

The unique contribution of manual chest compression-vibrations to airflow during physiotherapy in sedated, fully ventilated children.

PubMed

Gregson, Rachael K; Shannon, Harriet; Stocks, Janet; Cole, Tim J; Peters, Mark J; Main, Eleanor

2012-03-01

This study aimed to quantify the specific effects of manual lung inflations with chest compression-vibrations, commonly used to assist airway clearance in ventilated patients. The hypothesis was that force applied during the compressions made a significant additional contribution to increases in peak expiratory flow and expiratory to inspiratory flow ratio over and above that resulting from accompanying increases in inflation volume. Prospective observational study. Cardiac and general pediatric intensive care. Sedated, fully ventilated children. Customized force-sensing mats and a commercial respiratory monitor recorded force and respiration during physiotherapy. Percentage changes in peak expiratory flow, peak expiratory to inspiratory flow ratios, inflation volume, and peak inflation pressure between baseline and manual inflations with and without compression-vibrations were calculated. Analysis of covariance determined the relative contribution of changes in pressure, volume, and force to influence changes in peak expiratory flow and peak expiratory to inspiratory flow ratio. Data from 105 children were analyzed (median age, 1.3 yrs; range, 1 wk to 15.9 yrs). Force during compressions ranged from 15 to 179 N (median, 46 N). Peak expiratory flow increased on average by 76% during compressions compared with baseline ventilation. Increases in peak expiratory flow were significantly related to increases in inflation volume, peak inflation pressure, and force with peak expiratory flow increasing by, on average, 4% for every 10% increase in inflation volume (p < .001), 5% for every 10% increase in peak inflation pressure (p = .005), and 3% for each 10 N of applied force (p < .001). By contrast, increase in peak expiratory to inspiratory flow ratio was only related to applied force with a 4% increase for each 10 N of force (p < .001). These results provide evidence of the unique contribution of compression forces in increasing peak expiratory flow and peak expiratory to inspiratory flow ratio bias over and above that related to accompanying changes from manual hyperinflations. Force generated during compression-vibrations was the single significant factor in multivariable analysis to explain the increases in expiratory flow bias. Such increases in the expiratory bias provide theoretically optimal physiological conditions for cephalad mucus movement in fully ventilated children.

Aircraft energy efficiency laminar flow control glove flight conceptual design study

NASA Technical Reports Server (NTRS)

Wright, A. S.

1979-01-01

A laminar flow control glove applied to the wing of a short to medium range jet transport with aft mounted engines was designed. A slotted aluminum glove concept and a woven stainless steel mesh porous glove concept suction surfaces were studied. The laminar flow control glove and a dummy glove with a modified supercritical airfoil, ducting, modified wing leading and trailing edges, modified flaps, and an LFC trim tab were applied to the wing after slot spacing suction parameters, and compression power were determined. The results show that a laminar flow control glove can be applied to the wing of a jet transport with an appropriate suction system installed.

Role of Compressibility on Tsunami Propagation

NASA Astrophysics Data System (ADS)

Abdolali, Ali; Kirby, James T.

2017-12-01

In the present paper, we aim to reduce the discrepancies between tsunami arrival times evaluated from tsunami models and real measurements considering the role of ocean compressibility. We perform qualitative studies to reveal the phase speed reduction rate via a modified version of the Mild Slope Equation for Weakly Compressible fluid (MSEWC) proposed by Sammarco et al. (2013). The model is validated against a 3-D computational model. Physical properties of surface gravity waves are studied and compared with those for waves evaluated from an incompressible flow solver over realistic geometry for 2011 Tohoku-oki event, revealing reduction in phase speed.