Cascades and Dissipative Anomalies in Compressible Fluid Turbulence

NASA Astrophysics Data System (ADS)

Eyink, Gregory L.; Drivas, Theodore D.

2018-02-01

We investigate dissipative anomalies in a turbulent fluid governed by the compressible Navier-Stokes equation. We follow an exact approach pioneered by Onsager, which we explain as a nonperturbative application of the principle of renormalization-group invariance. In the limit of high Reynolds and Péclet numbers, the flow realizations are found to be described as distributional or "coarse-grained" solutions of the compressible Euler equations, with standard conservation laws broken by turbulent anomalies. The anomalous dissipation of kinetic energy is shown to be due not only to local cascade but also to a distinct mechanism called pressure-work defect. Irreversible heating in stationary, planar shocks with an ideal-gas equation of state exemplifies the second mechanism. Entropy conservation anomalies are also found to occur via two mechanisms: an anomalous input of negative entropy (negentropy) by pressure work and a cascade of negentropy to small scales. We derive "4 /5 th-law"-type expressions for the anomalies, which allow us to characterize the singularities (structure-function scaling exponents) required to sustain the cascades. We compare our approach with alternative theories and empirical evidence. It is argued that the "Big Power Law in the Sky" observed in electron density scintillations in the interstellar medium is a manifestation of a forward negentropy cascade or an inverse cascade of usual thermodynamic entropy.

Compression in Working Memory and Its Relationship with Fluid Intelligence

ERIC Educational Resources Information Center

Chekaf, Mustapha; Gauvrit, Nicolas; Guida, Alessandro; Mathy, Fabien

2018-01-01

Working memory has been shown to be strongly related to fluid intelligence; however, our goal is to shed further light on the process of information compression in working memory as a determining factor of fluid intelligence. Our main hypothesis was that compression in working memory is an excellent indicator for studying the relationship between…

Relation Between Pore Size and the Compressibility of a Confined Fluid

PubMed Central

Gor, Gennady Y.; Siderius, Daniel W.; Rasmussen, Christopher J.; Krekelberg, William P.; Shen, Vincent K.; Bernstein, Noam

2015-01-01

When a fluid is confined to a nanopore, its thermodynamic properties differ from the properties of a bulk fluid, so measuring such properties of the confined fluid can provide information about the pore sizes. Here we report a simple relation between the pore size and isothermal compressibility of argon confined in these pores. Compressibility is calculated from the fluctuations of the number of particles in the grand canonical ensemble using two different simulation techniques: conventional grand-canonical Monte Carlo and grand-canonical ensemble transition-matrix Monte Carlo. Our results provide a theoretical framework for extracting the information on the pore sizes of fluid-saturated samples by measuring the compressibility from ultrasonic experiments. PMID:26590541

Shock formation and the ideal shape of ramp compression waves

NASA Astrophysics Data System (ADS)

Swift, Damian C.; Kraus, Richard G.; Loomis, Eric N.; Hicks, Damien G.; McNaney, James M.; Johnson, Randall P.

2008-12-01

We derive expressions for shock formation based on the local curvature of the flow characteristics during dynamic compression. Given a specific ramp adiabat, calculated for instance from the equation of state for a substance, the ideal nonlinear shape for an applied ramp loading history can be determined. We discuss the region affected by lateral release, which can be presented in compact form for the ideal loading history. Example calculations are given for representative metals and plastic ablators. Continuum dynamics (hydrocode) simulations were in good agreement with the algebraic forms. Example applications are presented for several classes of laser-loading experiment, identifying conditions where shocks are desired but not formed, and where long-duration ramps are desired.

Hypersonic Magneto-Fluid-Dynamic Compression in Cylindrical Inlet

NASA Technical Reports Server (NTRS)

Shang, Joseph S.; Chang, Chau-Lyan

2007-01-01

Hypersonic magneto-fluid-dynamic interaction has been successfully performed as a virtual leading-edge strake and a virtual cowl of a cylindrical inlet. In a side-by-side experimental and computational study, the magnitude of the induced compression was found to be depended on configuration and electrode placement. To better understand the interacting phenomenon the present investigation is focused on a direct current discharge at the leading edge of a cylindrical inlet for which validating experimental data is available. The present computational result is obtained by solving the magneto-fluid-dynamics equations at the low magnetic Reynolds number limit and using a nonequilibrium weakly ionized gas model based on the drift-diffusion theory. The numerical simulation provides a detailed description of the intriguing physics. After validation with experimental measurements, the computed results further quantify the effectiveness of a magnet-fluid-dynamic compression for a hypersonic cylindrical inlet. At a minuscule power input to a direct current surface discharge of 8.14 watts per square centimeter of electrode area produces an additional compression of 6.7 percent for a constant cross-section cylindrical inlet.

An immersed boundary method for fluid-structure interaction with compressible multiphase flows

NASA Astrophysics Data System (ADS)

Wang, Li; Currao, Gaetano M. D.; Han, Feng; Neely, Andrew J.; Young, John; Tian, Fang-Bao

2017-10-01

This paper presents a two-dimensional immersed boundary method for fluid-structure interaction with compressible multiphase flows involving large structure deformations. This method involves three important parts: flow solver, structure solver and fluid-structure interaction coupling. In the flow solver, the compressible multiphase Navier-Stokes equations for ideal gases are solved by a finite difference method based on a staggered Cartesian mesh, where a fifth-order accuracy Weighted Essentially Non-Oscillation (WENO) scheme is used to handle spatial discretization of the convective term, a fourth-order central difference scheme is employed to discretize the viscous term, the third-order TVD Runge-Kutta scheme is used to discretize the temporal term, and the level-set method is adopted to capture the multi-material interface. In this work, the structure considered is a geometrically non-linear beam which is solved by using a finite element method based on the absolute nodal coordinate formulation (ANCF). The fluid dynamics and the structure motion are coupled in a partitioned iterative manner with a feedback penalty immersed boundary method where the flow dynamics is defined on a fixed Lagrangian grid and the structure dynamics is described on a global coordinate. We perform several validation cases (including fluid over a cylinder, structure dynamics, flow induced vibration of a flexible plate, deformation of a flexible panel induced by shock waves in a shock tube, an inclined flexible plate in a hypersonic flow, and shock-induced collapse of a cylindrical helium cavity in the air), and compare the results with experimental and other numerical data. The present results agree well with the published data and the current experiment. Finally, we further demonstrate the versatility of the present method by applying it to a flexible plate interacting with multiphase flows.

On the Singular Incompressible Limit of Inviscid Compressible Fluids

NASA Astrophysics Data System (ADS)

Secchi, P.

We consider the Euler equations of barotropic inviscid compressible fluids in a bounded domain. It is well known that, as the Mach number goes to zero, the compressible flows approximate the solution of the equations of motion of inviscid, incompressible fluids. In this paper we discuss, for the boundary case, the different kinds of convergence under various assumptions on the data, in particular the weak convergence in the case of uniformly bounded initial data and the strong convergence in the norm of the data space.

A Finite-Volume approach for compressible single- and two-phase flows in flexible pipelines with fluid-structure interaction

NASA Astrophysics Data System (ADS)

Daude, F.; Galon, P.

2018-06-01

A Finite-Volume scheme for the numerical computations of compressible single- and two-phase flows in flexible pipelines is proposed based on an approximate Godunov-type approach. The spatial discretization is here obtained using the HLLC scheme. In addition, the numerical treatment of abrupt changes in area and network including several pipelines connected at junctions is also considered. The proposed approach is based on the integral form of the governing equations making it possible to tackle general equations of state. A coupled approach for the resolution of fluid-structure interaction of compressible fluid flowing in flexible pipes is considered. The structural problem is solved using Euler-Bernoulli beam finite elements. The present Finite-Volume method is applied to ideal gas and two-phase steam-water based on the Homogeneous Equilibrium Model (HEM) in conjunction with a tabulated equation of state in order to demonstrate its ability to tackle general equations of state. The extensive application of the scheme for both shock tube and other transient flow problems demonstrates its capability to resolve such problems accurately and robustly. Finally, the proposed 1-D fluid-structure interaction model appears to be computationally efficient.

Fluid Compressibility Effects on the Dynamic Response of Hydrostatic Journal Bearings

NASA Technical Reports Server (NTRS)

Sanandres, Luis A.

1991-01-01

A theoretical analysis for the dynamic performance characteristics of laminar flow, capillar/orifice compensated hydrostatic journal bearings is presented. The analysis considers in detail the effect of fluid compressibility in the bearing recesses. At high frequency excitations beyond a break frequency, the bearing hydrostatic stiffness increases sharply and it is accompanied by a rapid decrease in direct damping. Also, the potential of pneumatic hammer instability (negative damping) at low frequencies is likely to occur in hydrostatic bearing applications handling highly compressible fluids. Useful design criteria to avoid undesirable dynamic operating conditions at low and high frequencies are determined. The effect of fluid recess compressibility is brought into perspective, and found to be of utmost importance on the entire frequency spectrum response and stability characteristics of hydrostatic/hybrid journal bearings.

Exact density functional theory for ideal polymer fluids with nearest neighbor bonding constraints.

PubMed

Woodward, Clifford E; Forsman, Jan

2008-08-07

We present a new density functional theory of ideal polymer fluids, assuming nearest-neighbor bonding constraints. The free energy functional is expressed in terms of end site densities of chain segments and thus has a simpler mathematical structure than previously used expressions using multipoint distributions. This work is based on a formalism proposed by Tripathi and Chapman [Phys. Rev. Lett. 94, 087801 (2005)]. Those authors obtain an approximate free energy functional for ideal polymers in terms of monomer site densities. Calculations on both repulsive and attractive surfaces show that their theory is reasonably accurate in some cases, but does differ significantly from the exact result for longer polymers with attractive surfaces. We suggest that segment end site densities, rather than monomer site densities, are the preferred choice of "site functions" for expressing the free energy functional of polymer fluids. We illustrate the application of our theory to derive an expression for the free energy of an ideal fluid of infinitely long polymers.

Thermodynamical effects and high resolution methods for compressible fluid flows

NASA Astrophysics Data System (ADS)

Li, Jiequan; Wang, Yue

2017-08-01

One of the fundamental differences of compressible fluid flows from incompressible fluid flows is the involvement of thermodynamics. This difference should be manifested in the design of numerical schemes. Unfortunately, the role of entropy, expressing irreversibility, is often neglected even though the entropy inequality, as a conceptual derivative, is verified for some first order schemes. In this paper, we refine the GRP solver to illustrate how the thermodynamical variation is integrated into the design of high resolution methods for compressible fluid flows and demonstrate numerically the importance of thermodynamic effects in the resolution of strong waves. As a by-product, we show that the GRP solver works for generic equations of state, and is independent of technical arguments.

Metering System for Compressible Fluids.

DTIC Science & Technology

1995-04-10

pressure switch and a low pass pressure switch are included in 5 line with the compressible fluid cylinder; consequently, the density of the...Once the pressure in first container 30 reaches the preset pressure for pressure switch 58, inlet valves 20 and 24 are closed and outlet valves 36...is allowed to drop to the preset pressure for pressure switch 60, at which time outlet valves 36 and 40 are closed, inlet valves 20 and 24 are

Simulating coupled dynamics of a rigid-flexible multibody system and compressible fluid

NASA Astrophysics Data System (ADS)

Hu, Wei; Tian, Qiang; Hu, HaiYan

2018-04-01

As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the smoothed particle hydrodynamics (SPH) method is used to model the compressible fluid, the natural coordinate formulation (NCF) and absolute nodal coordinate formulation (ANCF) are used to model the rigid and flexible bodies, respectively. In order to model the compressible fluid properly and efficiently via SPH method, three measures are taken as follows. The first is to use the Riemann solver to cope with the fluid compressibility, the second is to define virtual particles of SPH to model the dynamic interaction between the fluid and the multibody system, and the third is to impose the boundary conditions of periodical inflow and outflow to reduce the number of SPH particles involved in the computation process. Afterwards, a parallel computation strategy is proposed based on the graphics processing unit (GPU) to detect the neighboring SPH particles and to solve the dynamic equations of SPH particles in order to improve the computation efficiency. Meanwhile, the generalized-alpha algorithm is used to solve the dynamic equations of the multibody system. Finally, four case studies are given to validate the proposed parallel computation approach.

Transition-metal alloying of γ'-Ni3Al : Effects on the ideal uniaxial compressive strength from first-principles calculations

NASA Astrophysics Data System (ADS)

Wen, Minru; Wang, Chong-Yu

2018-01-01

The addition of transition-metal (TM) elements into the γ' precipitate phase of a Ni-based single-crystal superalloy can significantly affect its mechanical properties, including the intrinsic mechanical property of compressive strength. Using first-principles density functional calculations, the effects of 3 d (Sc-Zn), 4 d (Y-Cd), and 5 d (Hf-Au) TM alloying elements on the ideal uniaxial compressive strength of γ'-Ni3Al were investigated. The stress-strain relationships of pure Ni3Al under [100], [110], and [111] compressive loads and the site occupancy behavior of TM elements in Ni3Al were previously studied using a total-energy method based on density functional theory. Our results showed that the capacity of TM elements for strengthening the ideal compressive strength was associated with the d -electron number. The alloying elements with half-filled d bands (i.e., Cr, Mo, W, Tc, and Re) manifested the greatest efficacy for improving the ideal strength of Ni3Al under a deformation along the weakest compressive direction. Furthermore, the charge redistribution of Ni3Al doped with 5 d elements were also analyzed to understand the strengthening mechanisms of TM elements in the γ'-Ni3Al phase.

A Geophysical Flow Experiment in a Compressible Critical Fluid

NASA Technical Reports Server (NTRS)

Hegseth, John; Garcia, Laudelino

1996-01-01

The first objective of this experiment is to build an experimental system in which, in analogy to a geophysical system, a compressible fluid in a spherical annulus becomes radially stratified in density through an A.C. electric field. When this density gradient is demonstrated, the system will be augmented so that the fluid can be driven by heating and rotation and tested in preparation for a microgravity experiment. This apparatus consists of a spherical capacitor filled with critical fluid in a temperature controlled environment. To make the fluid critical, the apparatus will be operated near the critical pressure, critical density, and critical temperature of the fluid. This will result in a highly compressible fluid because of the properties of the fluid near its critical point. A high voltage A.C. source applied across the capacitor will create a spherically symmetric central force because of the dielectric properties of the fluid in an electric field gradient. This central force will induce a spherically symmetric density gradient that is analogous to a geophysical fluid system. To generate such a density gradient the system must be small (approx. 1 inch diameter). This small cell will also be capable of driving the critical fluid by heating and rotation. Since a spherically symmetric density gradient can only be made in microgravity, another small cell, of the same geometry, will be built that uses incompressible fluid. The driving of the fluid by rotation and heating in these small cells will be developed. The resulting instabilities from the driving in these two systems will then be studied. The second objective is to study the pattern forming instabilities (bifurcations) resulting from the well controlled experimental conditions in the critical fluid cell. This experiment will come close to producing conditions that are geophysically similar and will be studied as the driving parameters are changed.

The energy density distribution of an ideal gas and Bernoulli’s equations

NASA Astrophysics Data System (ADS)

Santos, Leonardo S. F.

2018-05-01

This work discusses the energy density distribution in an ideal gas and the consequences of Bernoulli’s equation and the corresponding relation for compressible fluids. The aim of this work is to study how Bernoulli’s equation determines the energy flow in a fluid, although Bernoulli’s equation does not describe the energy density itself. The model from molecular dynamic considerations that describes an ideal gas at rest with uniform density is modified to explore the gas in motion with non-uniform density and gravitational effects. The difference between the component of the speed of a particle that is parallel to the gas speed and the gas speed itself is called ‘parallel random speed’. The pressure from the ‘parallel random speed’ is denominated as parallel pressure. The modified model predicts that the energy density is the sum of kinetic and potential gravitational energy densities plus two terms with static and parallel pressures. The application of Bernoulli’s equation and the corresponding relation for compressible fluids in the energy density expression has resulted in two new formulations. For incompressible and compressible gas, the energy density expressions are written as a function of stagnation, static and parallel pressures, without any dependence on kinetic or gravitational potential energy densities. These expressions of the energy density are the main contributions of this work. When the parallel pressure was uniform, the energy density distribution for incompressible approximation and compressible gas did not converge to zero for the limit of null static pressure. This result is rather unusual because the temperature tends to zero for null pressure. When the gas was considered incompressible and the parallel pressure was equal to static pressure, the energy density maintained this unusual behaviour with small pressures. If the parallel pressure was equal to static pressure, the energy density converged to zero for the limit of the

Fluid displacement from intraluminal thrombus of abdominal aortic aneurysm as a result of uniform compression.

PubMed

van Noort, Kim; Schuurmann, Richte Cl; Wermelink, Bryan; Slump, Cornelis H; Kuijpers, Karel C; de Vries, Jean-Paul Pm

2017-10-01

Objectives The results after aneurysm repair with an endovascular aneurysm sealing (EVAS) system are dependent on the stability of the aneurysm sac and particularly the intraluminal abdominal aortic thrombus (ILT). The postprocedural ILT volume is decreased compared with preprocedural ILT volume in aortic aneurysm patients treated with EVAS. We hypothesize that ILT is not stable in all patients and pressurization of the ILT may result in displacement of fluids from the ILT, no differently than serum is displaced from whole blood when it settles. To date, the mechanism and quantification of fluid displacement from ILT are unknown. Methods The study included 21 patients who underwent elective open abdominal aortic aneurysm repair. The ILT was harvested as a routine procedure during the operation. After excision of a histologic sample of the ILT specimen in four patients, ILT volume was measured and the ILT was compressed in a dedicated compression setup designed to apply uniform compression of 200 mmHg for 5 min. After compression, the volumes of the remaining thrombus and the displaced fluid were measured. Results The median (interquartile-range) of ILT volume before compression was 60 (66) mL, and a median of 5.7 (8.4) mL of fluid was displaced from the ILT after compression, resulting in a median thrombus volume decrease of 11% (10%). Fluid components can be up to 31% of the entire ILT volume. Histologic examination of four ILT specimens showed a reduction of the medial layer of the ILT after compression, which was the result of compression of fluid-containing canaliculi. Conclusions Applying pressure of 200 mmHg to abdominal aortic aneurysm ILT resulted in the displacement of fluid, with a large variation among patients. Fluid displacement may result in decrease of ILT volume during and after EVAS, which might have implications on pre-EVAS volume planning and on stability of the endobags during follow-up which may lead to migration, endoleak or both.

Decay of the compressible magneto-micropolar fluids

NASA Astrophysics Data System (ADS)

Zhang, Peixin

2018-02-01

This paper considers the large-time behavior of solutions to the Cauchy problem on the compressible magneto-micropolar fluid system under small perturbation in regular Sobolev space. Based on the time-weighted energy estimate, the asymptotic stability of the steady state with the strictly positive constant density, vanishing velocity, micro-rotational velocity, and magnetic field is established.

The pointwise estimates of diffusion wave of the compressible micropolar fluids

NASA Astrophysics Data System (ADS)

Wu, Zhigang; Wang, Weike

2018-09-01

The pointwise estimates for the compressible micropolar fluids in dimension three are given, which exhibit generalized Huygens' principle for the fluid density and fluid momentum as the compressible Navier-Stokes equation, while the micro-rational momentum behaves like the fluid momentum of the Euler equation with damping. To circumvent the complexity from 7 × 7 Green's matrix, we use the decomposition of fluid part and electromagnetic part for the momentums to study three smaller Green's matrices. The following from this decomposition is that we have to deal with the new problem that the nonlinear terms contain nonlocal operators. We solve it by using the natural match of these new Green's functions and the nonlinear terms. Moreover, to derive the different pointwise estimates for different unknown variables such that the estimate of each unknown variable is in agreement with its Green's function, we develop some new estimates on the nonlinear interplay between different waves.

Conjugate Compressible Fluid Flow and Heat Transfer in Ducts

NASA Technical Reports Server (NTRS)

Cross, M. F.

2011-01-01

A computational approach to modeling transient, compressible fluid flow with heat transfer in long, narrow ducts is presented. The primary application of the model is for analyzing fluid flow and heat transfer in solid propellant rocket motor nozzle joints during motor start-up, but the approach is relevant to a wide range of analyses involving rapid pressurization and filling of ducts. Fluid flow is modeled through solution of the spatially one-dimensional, transient Euler equations. Source terms are included in the governing equations to account for the effects of wall friction and heat transfer. The equation solver is fully-implicit, thus providing greater flexibility than an explicit solver. This approach allows for resolution of pressure wave effects on the flow as well as for fast calculation of the steady-state solution when a quasi-steady approach is sufficient. Solution of the one-dimensional Euler equations with source terms significantly reduces computational run times compared to general purpose computational fluid dynamics packages solving the Navier-Stokes equations with resolved boundary layers. In addition, conjugate heat transfer is more readily implemented using the approach described in this paper than with most general purpose computational fluid dynamics packages. The compressible flow code has been integrated with a transient heat transfer solver to analyze heat transfer between the fluid and surrounding structure. Conjugate fluid flow and heat transfer solutions are presented. The author is unaware of any previous work available in the open literature which uses the same approach described in this paper.

The Role of Fluid Compression in Particle Energization during Magnetic Reconnection

NASA Astrophysics Data System (ADS)

Li, X.; Guo, F.; Li, H.; Li, S.

2017-12-01

Theories of particle transport and acceleration have shown that fluid compression is the leading mechanism for particle energization. However, the role of compression in particle energization during magnetic reconnection is unclear. We present a cluster of studies to clarify and show the effect of fluid compression in accelerating particles to high energies during magnetic reconnection. Using fully kinetic reconnection simulations, we show that fluid compression is the leading mechanism for high-energy particle energization. We find that the compressional energization is more important in a low-beta plasma or in a reconnection layer with a weak guide field (the magnetic field component perpendicular to the reconnecting magnetic field), which are relevant to solar flares. Our analysis on 3D kinetic simulations shows that the self-generated turbulence scatters particles and enhances the particle diffusion processes in the acceleration regions. Based on these results, we then study large-scale reconnection acceleration by solving the particle transport equation in a large-scale reconnection layer evolved with MHD simulations. Due to the compressional effect, particles are accelerated to high energies and develop power-law energy distributions. This study clarifies the nature of particle acceleration in reconnection layer and is important to understand particle energization during large-scale acceleration such as solar flares.

Velocity relaxation of a particle in a confined compressible fluid

NASA Astrophysics Data System (ADS)

Tatsumi, Rei; Yamamoto, Ryoichi

2013-05-01

The velocity relaxation of an impulsively forced spherical particle in a fluid confined by two parallel plane walls is studied using a direct numerical simulation approach. During the relaxation process, the momentum of the particle is transmitted in the ambient fluid by viscous diffusion and sound wave propagation, and the fluid flow accompanied by each mechanism has a different character and affects the particle motion differently. Because of the bounding walls, viscous diffusion is hampered, and the accompanying shear flow is gradually diminished. However, the sound wave is repeatedly reflected and spreads diffusely. As a result, the particle motion is governed by the sound wave and backtracks differently in a bulk fluid. The time when the backtracking of the particle occurs changes non-monotonically with respect to the compressibility factor ɛ = ν/ac and is minimized at the characteristic compressibility factor. This factor depends on the wall spacing, and the dependence is different at small and large wall spacing regions based on the different mechanisms causing the backtracking.

Fully-coupled aeroelastic simulation with fluid compressibility — For application to vocal fold vibration

PubMed Central

Yang, Jubiao; Wang, Xingshi; Krane, Michael; Zhang, Lucy T.

2017-01-01

In this study, a fully-coupled fluid–structure interaction model is developed for studying dynamic interactions between compressible fluid and aeroelastic structures. The technique is built based on the modified Immersed Finite Element Method (mIFEM), a robust numerical technique to simulate fluid–structure interactions that has capabilities to simulate high Reynolds number flows and handles large density disparities between the fluid and the solid. For accurate assessment of this intricate dynamic process between compressible fluid, such as air and aeroelastic structures, we included in the model the fluid compressibility in an isentropic process and a solid contact model. The accuracy of the compressible fluid solver is verified by examining acoustic wave propagations in a closed and an open duct, respectively. The fully-coupled fluid–structure interaction model is then used to simulate and analyze vocal folds vibrations using compressible air interacting with vocal folds that are represented as layered viscoelastic structures. Using physiological geometric and parametric setup, we are able to obtain a self-sustained vocal fold vibration with a constant inflow pressure. Parametric studies are also performed to study the effects of lung pressure and vocal fold tissue stiffness in vocal folds vibrations. All the case studies produce expected airflow behavior and a sustained vibration, which provide verification and confidence in our future studies of realistic acoustical studies of the phonation process. PMID:29527067

Formation of Hydro-acoustic Waves in Dissipative Coupled Weakly Compressible Fluids

NASA Astrophysics Data System (ADS)

Abdolali, A.; Kirby, J. T., Jr.; Bellotti, G.

2014-12-01

Recent advances in deep sea measurement technology provide an increasing opportunity to detect and interpret hydro-acoustic waves as a component in improved Tsunami Early Warning Systems (TEWS). For the idealized case of a homogeneous water column above a moving but otherwise rigid bottom (in terms of assessing acoustic wave interaction), the description of the infinite family of acoustic modes is characterized by local water depth at source area; i.e. the period of the first acoustic mode is given by four times the required time for sound to travel from the seabed to the surface. Spreading off from earthquake zone, the dominant spectrum is filtered and enriched by seamounts and barriers. This study focuses on the characteristics of hydro-acoustic waves generated by sudden sea bottom motion in a weakly compressible fluid coupled with an underlying sedimentary layer, where the added complexity of the sediment layer rheology leads to both the lowering of dominant spectral peaks and wave attenuation across the full spectrum. To overcome the computational difficulties of three-dimensional models, we derive a depth integrated equation valid for varying water depth and sediment thickness. Damping behavior of the two layered system is initially taken into account by introducing the viscosity of fluid-like sedimentary layer. We show that low frequency pressure waves which are precursor components of tsunamis contain information of seafloor motion.

Pressure and compressibility factor of bidisperse magnetic fluids

NASA Astrophysics Data System (ADS)

Minina, Elena S.; Blaak, Ronald; Kantorovich, Sofia S.

2018-04-01

In this work, we investigate the pressure and compressibility factors of bidisperse magnetic fluids with relatively weak dipolar interactions and different granulometric compositions. In order to study these properties, we employ the method of diagram expansion, taking into account two possible scenarios: (1) dipolar particles repel each other as hard spheres; (2) the polymer shell on the surface of the particles is modelled through a soft-sphere approximation. The theoretical predictions of the pressure and compressibility factors of bidisperse ferrofluids at different granulometric compositions are supported by data obtained by means of molecular dynamics computer simulations, which we also carried out for these systems. Both theory and simulations reveal that the pressure and compressibility factors decrease with growing dipolar correlations in the system, namely with an increasing fraction of large particles. We also demonstrate that even if dipolar interactions are too weak for any self-assembly to take place, the interparticle correlations lead to a qualitative change in the behaviour of the compressibility factors when compared to that of non-dipolar spheres, making the dependence monotonic.

Exact Theory of Compressible Fluid Turbulence

NASA Astrophysics Data System (ADS)

Drivas, Theodore; Eyink, Gregory

2017-11-01

We obtain exact results for compressible turbulence with any equation of state, using coarse-graining/filtering. We find two mechanisms of turbulent kinetic energy dissipation: scale-local energy cascade and ``pressure-work defect'', or pressure-work at viscous scales exceeding that in the inertial-range. Planar shocks in an ideal gas dissipate all kinetic energy by pressure-work defect, but the effect is omitted by standard LES modeling of pressure-dilatation. We also obtain a novel inverse cascade of thermodynamic entropy, injected by microscopic entropy production, cascaded upscale, and removed by large-scale cooling. This nonlinear process is missed by the Kovasznay linear mode decomposition, treating entropy as a passive scalar. For small Mach number we recover the incompressible ``negentropy cascade'' predicted by Obukhov. We derive exact Kolmogorov 4/5th-type laws for energy and entropy cascades, constraining scaling exponents of velocity, density, and internal energy to sub-Kolmogorov values. Although precise exponents and detailed physics are Mach-dependent, our exact results hold at all Mach numbers. Flow realizations at infinite Reynolds are ``dissipative weak solutions'' of compressible Euler equations, similarly as Onsager proposed for incompressible turbulence.

Causal implications of viscous damping in compressible fluid flows

PubMed

Jordan; Meyer; Puri

2000-12-01

Classically, a compressible, isothermal, viscous fluid is regarded as a mathematical continuum and its motion is governed by the linearized continuity, Navier-Stokes, and state equations. Unfortunately, solutions of this system are of a diffusive nature and hence do not satisfy causality. However, in the case of a half-space of fluid set to motion by a harmonically vibrating plate the classical equation of motion can, under suitable conditions, be approximated by the damped wave equation. Since this equation is hyperbolic, the resulting solutions satisfy causal requirements. In this work the Laplace transform and other analytical and numerical tools are used to investigate this apparent contradiction. To this end the exact solutions, as well as their special and limiting cases, are found and compared for the two models. The effects of the physical parameters on the solutions and associated quantities are also studied. It is shown that propagating wave fronts are only possible under the hyperbolic model and that the concept of phase speed has different meanings in the two formulations. In addition, discontinuities and shock waves are noted and a physical system is modeled under both formulations. Overall, it is shown that the hyperbolic form gives a more realistic description of the physical problem than does the classical theory. Lastly, a simple mechanical analog is given and connections to viscoelastic fluids are noted. In particular, the research presented here supports the notion that linear compressible, isothermal, viscous fluids can, at least in terms of causality, be better characterized as a type of viscoelastic fluid.

Compressible Fluids Interacting with a Linear-Elastic Shell

NASA Astrophysics Data System (ADS)

Breit, Dominic; Schwarzacher, Sebastian

2018-05-01

We study the Navier-Stokes equations governing the motion of an isentropic compressible fluid in three dimensions interacting with a flexible shell of Koiter type. The latter one constitutes a moving part of the boundary of the physical domain. Its deformation is modeled by a linearized version of Koiter's elastic energy. We show the existence of weak solutions to the corresponding system of PDEs provided the adiabatic exponent satisfies {γ > 12/7} ({γ >1 } in two dimensions). The solution exists until the moving boundary approaches a self-intersection. This provides a compressible counterpart of the results in Lengeler and Růžičkaka (Arch Ration Mech Anal 211(1):205-255, 2014) on incompressible Navier-Stokes equations.

University Students Explaining Adiabatic Compression of an Ideal Gas--A New Phenomenon in Introductory Thermal Physics

ERIC Educational Resources Information Center

Leinonen, Risto; Asikainen, Mervi A.; Hirvonen, Pekka E.

2012-01-01

This study focuses on second-year university students' explanations and reasoning related to adiabatic compression of an ideal gas. The phenomenon was new to the students, but it was one which they should have been capable of explaining using their previous upper secondary school knowledge. The students' explanations and reasoning were…

Advanced computational techniques for incompressible/compressible fluid-structure interactions

NASA Astrophysics Data System (ADS)

Kumar, Vinod

2005-07-01

Fluid-Structure Interaction (FSI) problems are of great importance to many fields of engineering and pose tremendous challenges to numerical analyst. This thesis addresses some of the hurdles faced for both 2D and 3D real life time-dependent FSI problems with particular emphasis on parachute systems. The techniques developed here would help improve the design of parachutes and are of direct relevance to several other FSI problems. The fluid system is solved using the Deforming-Spatial-Domain/Stabilized Space-Time (DSD/SST) finite element formulation for the Navier-Stokes equations of incompressible and compressible flows. The structural dynamics solver is based on a total Lagrangian finite element formulation. Newton-Raphson method is employed to linearize the otherwise nonlinear system resulting from the fluid and structure formulations. The fluid and structural systems are solved in decoupled fashion at each nonlinear iteration. While rigorous coupling methods are desirable for FSI simulations, the decoupled solution techniques provide sufficient convergence in the time-dependent problems considered here. In this thesis, common problems in the FSI simulations of parachutes are discussed and possible remedies for a few of them are presented. Further, the effects of the porosity model on the aerodynamic forces of round parachutes are analyzed. Techniques for solving compressible FSI problems are also discussed. Subsequently, a better stabilization technique is proposed to efficiently capture and accurately predict the shocks in supersonic flows. The numerical examples simulated here require high performance computing. Therefore, numerical tools using distributed memory supercomputers with message passing interface (MPI) libraries were developed.

Local invariants in non-ideal flows of neutral fluids and two-fluid plasmas

NASA Astrophysics Data System (ADS)

Zhu, Jian-Zhou

2018-03-01

The main objective is the locally invariant geometric object of any (magneto-)fluid dynamics with forcing and damping (nonideal), while more attention is paid to the untouched dynamical properties of two-fluid fashion. Specifically, local structures, beyond the well-known "frozen-in" to the barotropic flows of the generalized vorticities, of the two-fluid model of plasma flows are presented. More general non-barotropic situations are also considered. A modified Euler equation [T. Tao, "Finite time blowup for Lagrangian modifications of the three-dimensional Euler equation," Ann. PDE 2, 9 (2016)] is also accordingly analyzed and remarked from the angle of view of the two-fluid model, with emphasis on the local structures. The local constraints of high-order differential forms such as helicity, among others, find simple formulation for possible practices in modeling the dynamics. Thus, the Cauchy invariants equation [N. Besse and U. Frisch, "Geometric formulation of the Cauchy invariants for incompressible Euler flow in flat and curved spaces," J. Fluid Mech. 825, 412 (2017)] may be enabled to find applications in non-ideal flows. Some formal examples are offered to demonstrate the calculations, and particularly interestingly the two-dimensional-three-component (2D3C) or the 2D passive scalar problem presents that a locally invariant Θ = 2θζ, with θ and ζ being, respectively, the scalar value of the "vertical velocity" (or the passive scalar) and the "vertical vorticity," may be used as if it were the spatial density of the globally invariant helicity, providing a Lagrangian prescription to control the latter in some situations of studying its physical effects in rapidly rotating flows (ubiquitous in atmosphere of astrophysical objects) with marked 2D3C vortical modes or in purely 2D passive scalars.

A fast numerical method for ideal fluid flow in domains with multiple stirrers

NASA Astrophysics Data System (ADS)

Nasser, Mohamed M. S.; Green, Christopher C.

2018-03-01

A collection of arbitrarily-shaped solid objects, each moving at a constant speed, can be used to mix or stir ideal fluid, and can give rise to interesting flow patterns. Assuming these systems of fluid stirrers are two-dimensional, the mathematical problem of resolving the flow field—given a particular distribution of any finite number of stirrers of specified shape and speed—can be formulated as a Riemann-Hilbert (R-H) problem. We show that this R-H problem can be solved numerically using a fast and accurate algorithm for any finite number of stirrers based around a boundary integral equation with the generalized Neumann kernel. Various systems of fluid stirrers are considered, and our numerical scheme is shown to handle highly multiply connected domains (i.e. systems of many fluid stirrers) with minimal computational expense.

CFD analysis of thermally induced thermodynamic losses in the reciprocating compression and expansion of real gases

NASA Astrophysics Data System (ADS)

Taleb, Aly I.; Sapin, Paul; Barfuß, Christoph; Fabris, Drazen; Markides, Christos N.

2017-03-01

compared to the real-gas model for heavier gases. This discrepancy is most pronounced at rotational speeds where the losses are highest. The real-gas model predicts a peak loss of 8.9% of the compression work, while the ideal-gas model predicts a peak loss of 5.7%. These differences in the work loss are due to the fact that the gas behaves less ideally during expansion than during compression, with the compressibility factor being lower during compression. This behaviour cannot be captured with the ideal-gas law. It is concluded that real-gas effects must be taken into account in order to predict accurately the thermally induced loss mechanism when using heavy fluid molecules in such devices.

Methods for compressible fluid simulation on GPUs using high-order finite differences

NASA Astrophysics Data System (ADS)

Pekkilä, Johannes; Väisälä, Miikka S.; Käpylä, Maarit J.; Käpylä, Petri J.; Anjum, Omer

2017-08-01

We focus on implementing and optimizing a sixth-order finite-difference solver for simulating compressible fluids on a GPU using third-order Runge-Kutta integration. Since graphics processing units perform well in data-parallel tasks, this makes them an attractive platform for fluid simulation. However, high-order stencil computation is memory-intensive with respect to both main memory and the caches of the GPU. We present two approaches for simulating compressible fluids using 55-point and 19-point stencils. We seek to reduce the requirements for memory bandwidth and cache size in our methods by using cache blocking and decomposing a latency-bound kernel into several bandwidth-bound kernels. Our fastest implementation is bandwidth-bound and integrates 343 million grid points per second on a Tesla K40t GPU, achieving a 3 . 6 × speedup over a comparable hydrodynamics solver benchmarked on two Intel Xeon E5-2690v3 processors. Our alternative GPU implementation is latency-bound and achieves the rate of 168 million updates per second.

Efficient Simulation of Compressible, Viscous Fluids using Multi-rate Time Integration

NASA Astrophysics Data System (ADS)

Mikida, Cory; Kloeckner, Andreas; Bodony, Daniel

2017-11-01

In the numerical simulation of problems of compressible, viscous fluids with single-rate time integrators, the global timestep used is limited to that of the finest mesh point or fastest physical process. This talk discusses the application of multi-rate Adams-Bashforth (MRAB) integrators to an overset mesh framework to solve compressible viscous fluid problems of varying scale with improved efficiency, with emphasis on the strategy of timescale separation and the application of the resulting numerical method to two sample problems: subsonic viscous flow over a cylinder and a viscous jet in crossflow. The results presented indicate the numerical efficacy of MRAB integrators, outline a number of outstanding code challenges, demonstrate the expected reduction in time enabled by MRAB, and emphasize the need for proper load balancing through spatial decomposition in order for parallel runs to achieve the predicted time-saving benefit. This material is based in part upon work supported by the Department of Energy, National Nuclear Security Administration, under Award Number DE-NA0002374.

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

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.

Multidomain approach for calculating compressible flows

NASA Technical Reports Server (NTRS)

Cambier, L.; Chazzi, W.; Veuillot, J. P.; Viviand, H.

1982-01-01

A multidomain approach for calculating compressible flows by using unsteady or pseudo-unsteady methods is presented. This approach is based on a general technique of connecting together two domains in which hyperbolic systems (that may differ) are solved with the aid of compatibility relations associated with these systems. Some examples of this approach's application to calculating transonic flows in ideal fluids are shown, particularly the adjustment of shock waves. The approach is then applied to treating a shock/boundary layer interaction problem in a transonic channel.

Lagrangian particle method for compressible fluid dynamics

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

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

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 / multi-phase 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 extremalmore » points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. 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. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.« less

Lagrangian particle method for compressible fluid dynamics

DOE PAGES

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

2018-02-09

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 / multi-phase 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 extremalmore » points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. 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. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.« less

Ideal magnetohydrodynamic theory for localized interchange modes in toroidal anisotropic plasmas

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

Shi, Tonghui, E-mail: thshi@ipp.ac.cn; Wan, B. N.; Sun, Y.

2016-08-15

Ideal magnetohydrodynamic theory for localized interchange modes is developed for toroidal plasmas with anisotropic pressure. The work extends the existing theories of Johnson and Hastie [Phys. Fluids 31, 1609 (1988)], etc., to the low n mode case, where n is the toroidal mode number. Also, the plasma compressibility is included, so that the coupling of the parallel motion to perpendicular one, i.e., the so-called apparent mass effect, is investigated in the anisotropic pressure case. The singular layer equation is obtained, and the generalized Mercier's criterion is derived.

Solution of weakly compressible isothermal flow in landfill gas collection networks

NASA Astrophysics Data System (ADS)

Nec, Y.; Huculak, G.

2017-12-01

Pipe networks collecting gas in sanitary landfills operate under the regime of a weakly compressible isothermal flow of ideal gas. The effect of compressibility has been traditionally neglected in this application in favour of simplicity, thereby creating a conceptual incongruity between the flow equations and thermodynamic equation of state. Here the flow is solved by generalisation of the classic Darcy-Weisbach equation for an incompressible steady flow in a pipe to an ordinary differential equation, permitting continuous variation of density, viscosity and related fluid parameters, as well as head loss or gain due to gravity, in isothermal flow. The differential equation is solved analytically in the case of ideal gas for a single edge in the network. Thereafter the solution is used in an algorithm developed to construct the flow equations automatically for a network characterised by an incidence matrix, and determine pressure distribution, flow rates and all associated parameters therein.

Comment on "Proposal of a critical test of the Navier-Stokes-Fourier paradigm for compressible fluid continua".

PubMed

Felderhof, B U

2013-08-01

Recently, a critical test of the Navier-Stokes-Fourier equations for compressible fluid continua was proposed [H. Brenner, Phys. Rev. E 87, 013014 (2013)]. It was shown that the equations of bivelocity hydrodynamics imply that a compressible fluid in an isolated rotating circular cylinder attains a nonequilibrium steady state with a nonuniform temperature increasing radially with distance from the axis. We demonstrate that statistical mechanical arguments, involving Hamiltonian dynamics and ergodicity due to irregularity of the wall, lead instead to a thermal equilibrium state with uniform temperature. This is the situation to be expected in experiment.

Electrical resistivity of fluid methane multiply shock compressed to 147 GPa

NASA Astrophysics Data System (ADS)

Wang, Yi-Gao; Liu, Fu-Sheng; Liu, Qi-Jun; Wang, Wen-Peng

2018-01-01

Shock wave experiments were carried out to measure the electrical resistivity of fluid methane. The pressure range of 89-147 GPa and the temperature range from 1800 to 2600 K were achieved with a two-stage light-gas gun. We obtained a minimum electrical resistivity value of 4.5 × 10-2 Ω cm at pressure and temperature of 147 GPa and 2600 K, which is two orders of magnitude higher than that of hydrogen under similar conditions. The data are interpreted in terms of a continuous transition from insulator to semiconductor state. One possibility reason is chemical decomposition of methane in the shock compression process. Along density and temperature increase with Hugoniot pressure, dissociation of fluid methane increases continuously to form a H2-rich fluid.

Hydraulic permeability of multilayered collagen gel scaffolds under plastic compression-induced unidirectional fluid flow.

PubMed

Serpooshan, Vahid; Quinn, Thomas M; Muja, Naser; Nazhat, Showan N

2013-01-01

Under conditions of free fluid flow, highly hydrated fibrillar collagen gels expel fluid and undergo gravity driven consolidation (self-compression; SC). This process can be accelerated by the application of a compressive stress (plastic compression; PC) in order to generate dense collagen scaffolds for tissue engineering. To define the microstructural evolution of collagen gels under PC, this study applied a two-layer micromechanical model that was previously developed to measure hydraulic permeability (k) under SC. Radially confined PC resulted in unidirectional fluid flow through the gel and the formation of a dense lamella at the fluid expulsion boundary which was confirmed by confocal microscopy of collagen immunoreactivity. Gel mass loss due to PC and subsequent SC were measured and applied to Darcy's law to calculate the thickness of the lamella and hydrated layer, as well as their relative permeabilities. Increasing PC level resulted in a significant increase in mass loss fraction and lamellar thickness, while the thickness of the hydrated layer dramatically decreased. Permeability of lamella also decreased from 1.8×10(-15) to 1.0×10(-15) m(2) in response to an increase in PC level. Ongoing SC, following PC, resulted in a uniform decrease in mass loss and k with increasing PC level and as a function SC time. Experimental k data were in close agreement with those estimated by the Happel model. Calculation of average k values for various two-layer microstructures indicated that they each approached 10(-15)-10(-14) m(2) at equilibrium. In summary, the two-layer micromechanical model can be used to define the microstructure and permeability of multi-layered biomimetic scaffolds generated by PC. Copyright © 2012 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Fluid mechanics in fluids at rest.

PubMed

Brenner, Howard

2012-07-01

Using readily available experimental thermophoretic particle-velocity data it is shown, contrary to current teachings, that for the case of compressible flows independent dye- and particle-tracer velocity measurements of the local fluid velocity at a point in a flowing fluid do not generally result in the same fluid velocity measure. Rather, tracer-velocity equality holds only for incompressible flows. For compressible fluids, each type of tracer is shown to monitor a fundamentally different fluid velocity, with (i) a dye (or any other such molecular-tagging scheme) measuring the fluid's mass velocity v appearing in the continuity equation and (ii) a small, physicochemically and thermally inert, macroscopic (i.e., non-Brownian), solid particle measuring the fluid's volume velocity v(v). The term "compressibility" as used here includes not only pressure effects on density, but also temperature effects thereon. (For example, owing to a liquid's generally nonzero isobaric coefficient of thermal expansion, nonisothermal liquid flows are to be regarded as compressible despite the general perception of liquids as being incompressible.) Recognition of the fact that two independent fluid velocities, mass- and volume-based, are formally required to model continuum fluid behavior impacts on the foundations of contemporary (monovelocity) fluid mechanics. Included therein are the Navier-Stokes-Fourier equations, which are now seen to apply only to incompressible fluids (a fact well-known, empirically, to experimental gas kineticists). The findings of a difference in tracer velocities heralds the introduction into fluid mechanics of a general bipartite theory of fluid mechanics, bivelocity hydrodynamics [Brenner, Int. J. Eng. Sci. 54, 67 (2012)], differing from conventional hydrodynamics in situations entailing compressible flows and reducing to conventional hydrodynamics when the flow is incompressible, while being applicable to both liquids and gases.

Astrophysical fluid dynamics

NASA Astrophysics Data System (ADS)

Ogilvie, Gordon I.

2016-06-01

> These lecture notes and example problems are based on a course given at the University of Cambridge in Part III of the Mathematical Tripos. Fluid dynamics is involved in a very wide range of astrophysical phenomena, such as the formation and internal dynamics of stars and giant planets, the workings of jets and accretion discs around stars and black holes and the dynamics of the expanding Universe. Effects that can be important in astrophysical fluids include compressibility, self-gravitation and the dynamical influence of the magnetic field that is `frozen in' to a highly conducting plasma. The basic models introduced and applied in this course are Newtonian gas dynamics and magnetohydrodynamics (MHD) for an ideal compressible fluid. The mathematical structure of the governing equations and the associated conservation laws are explored in some detail because of their importance for both analytical and numerical methods of solution, as well as for physical interpretation. Linear and nonlinear waves, including shocks and other discontinuities, are discussed. The spherical blast wave resulting from a supernova, and involving a strong shock, is a classic problem that can be solved analytically. Steady solutions with spherical or axial symmetry reveal the physics of winds and jets from stars and discs. The linearized equations determine the oscillation modes of astrophysical bodies, as well as their stability and their response to tidal forcing.

The Investigation of Ghost Fluid Method for Simulating the Compressible Two-Medium Flow

NASA Astrophysics Data System (ADS)

Lu, Hai Tian; Zhao, Ning; Wang, Donghong

2016-06-01

In this paper, we investigate the conservation error of the two-dimensional compressible two-medium flow simulated by the front tracking method. As the improved versions of the original ghost fluid method, the modified ghost fluid method and the real ghost fluid method are selected to define the interface boundary conditions, respectively, to show different effects on the conservation error. A Riemann problem is constructed along the normal direction of the interface in the front tracking method, with the goal of obtaining an efficient procedure to track the explicit sharp interface precisely. The corresponding Riemann solutions are also used directly in these improved ghost fluid methods. Extensive numerical examples including the sod tube and the shock-bubble interaction are tested to calculate the conservation error. It is found that these two ghost fluid methods have distinctive performances for different initial conditions of the flow field, and the related conclusions are made to suggest the best choice for the combination.

Effect of pore geometry on the compressibility of a confined simple fluid

NASA Astrophysics Data System (ADS)

Dobrzanski, Christopher D.; Maximov, Max A.; Gor, Gennady Y.

2018-02-01

Fluids confined in nanopores exhibit properties different from the properties of the same fluids in bulk; among these properties is the isothermal compressibility or elastic modulus. The modulus of a fluid in nanopores can be extracted from ultrasonic experiments or calculated from molecular simulations. Using Monte Carlo simulations in the grand canonical ensemble, we calculated the modulus for liquid argon at its normal boiling point (87.3 K) adsorbed in model silica pores of two different morphologies and various sizes. For spherical pores, for all the pore sizes (diameters) exceeding 2 nm, we obtained a logarithmic dependence of fluid modulus on the vapor pressure. Calculation of the modulus at saturation showed that the modulus of the fluid in spherical pores is a linear function of the reciprocal pore size. The calculation of the modulus of the fluid in cylindrical pores appeared too scattered to make quantitative conclusions. We performed additional simulations at higher temperature (119.6 K), at which Monte Carlo insertions and removals become more efficient. The results of the simulations at higher temperature confirmed both regularities for cylindrical pores and showed quantitative difference between the fluid moduli in pores of different geometries. Both of the observed regularities for the modulus stem from the Tait-Murnaghan equation applied to the confined fluid. Our results, along with the development of the effective medium theories for nanoporous media, set the groundwork for analysis of the experimentally measured elastic properties of fluid-saturated nanoporous materials.

The role of fluid compression in energy conversion and particle energization during magnetic reconnection

NASA Astrophysics Data System (ADS)

Li, X.; Guo, F.; Li, G.; Li, H.

2016-12-01

Theories of particle transport and acceleration have shown that fluid compression is the leading mechanism for particle acceleration and plasma energization. However, the role of compression in particle acceleration during magnetic reconnection is unclear. We use two approaches to study this issue. First, using fully kinetic simulations, we quantitatively calculate the effect of compression in energy conversion and particle energization during magnetic reconnection for a range of plasma beta and guide field. We show that compression has an important contribution for the energy conversion between the bulk kinetic energy and the internal energy when the guide field is smaller than the reconnecting component. Based on this result, we then study the large-scale reconnection acceleration by solving the Parker's transport equation in a background reconnecting flow provided by MHD simulations. Due to the compression effect, the simulations suggest fast particle acceleration to high energies in the reconnection layer. This study clarifies the nature of particle acceleration in reconnection layer, and may be important to understand particle acceleration and plasma energization during solar flares.

Propagation of exponential shock wave in an axisymmetric rotating non-ideal dusty gas

NASA Astrophysics Data System (ADS)

Nath, G.

2016-09-01

One-dimensional unsteady isothermal and adiabatic flow behind a strong exponential shock wave propagating in a rotational axisymmetric mixture of non-ideal gas and small solid particles, which has variable azimuthal and axial fluid velocities, is analyzed. The shock wave is driven out by a piston moving with time according to exponential law. The azimuthal and axial components of the fluid velocity in the ambient medium are assumed to be varying and obeying exponential laws. In the present work, small solid particles are considered as pseudo-fluid with the assumption that the equilibrium flow-conditions are maintained in the flow-field, and the viscous-stress and heat conduction of the mixture are negligible. Solutions are obtained in both the cases, when the flow between the shock and the piston is isothermal or adiabatic by taking into account the components of vorticity vector and compressibility. It is found that the assumption of zero temperature gradient brings a profound change in the density, axial component of vorticity vector and compressibility distributions as compared to that of the adiabatic case. To investigate the behavior of the flow variables and the influence on the shock wave propagation by the parameter of non-idealness of the gas overline{b} in the mixture as well as by the mass concentration of solid particles in the mixture Kp and by the ratio of the density of solid particles to the initial density of the gas G1 are worked out in detail. It is interesting to note that the shock strength increases with an increase in G1 ; whereas it decreases with an increase in overline{b} . Also, a comparison between the solutions in the cases of isothermal and adiabatic flows is made.

Self-propulsion of free solid bodies with internal rotors via localized singular vortex shedding in planar ideal fluids

NASA Astrophysics Data System (ADS)

Tallapragada, P.; Kelly, S. D.

2015-11-01

Diverse mechanisms for animal locomotion in fluids rely on vortex shedding to generate propulsive forces. This is a complex phenomenon that depends essentially on fluid viscosity, but its influence can be modeled in an inviscid setting by introducing localized velocity constraints to systems comprising solid bodies interacting with ideal fluids. In the present paper, we invoke an unsteady version of the Kutta condition from inviscid airfoil theory and a more primitive stagnation condition to model vortex shedding from a geometrically contrasting pair of free planar bodies representing idealizations of swimming animals or robotic vehicles. We demonstrate with simulations that these constraints are sufficient to enable both bodies to propel themselves with very limited actuation. The solitary actuator in each case is a momentum wheel internal to the body, underscoring the symmetry-breaking role played by vortex shedding in converting periodic variations in a generic swimmer's angular momentum to forward locomotion. The velocity constraints are imposed discretely in time, resulting in the shedding of discrete vortices; we observe the roll-up of these vortices into distinctive wake structures observed in viscous models and physical experiments.

Nonlinear theory of magnetohydrodynamic flows of a compressible fluid in the shallow water approximation

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

Klimachkov, D. A., E-mail: klimchakovdmitry@gmail.com; Petrosyan, A. S., E-mail: apetrosy@iki.rssi.ru

2016-09-15

Shallow water magnetohydrodynamic (MHD) theory describing incompressible flows of plasma is generalized to the case of compressible flows. A system of MHD equations is obtained that describes the flow of a thin layer of compressible rotating plasma in a gravitational field in the shallow water approximation. The system of quasilinear hyperbolic equations obtained admits a complete simple wave analysis and a solution to the initial discontinuity decay problem in the simplest version of nonrotating flows. In the new equations, sound waves are filtered out, and the dependence of density on pressure on large scales is taken into account that describesmore » static compressibility phenomena. In the equations obtained, the mass conservation law is formulated for a variable that nontrivially depends on the shape of the lower boundary, the characteristic vertical scale of the flow, and the scale of heights at which the variation of density becomes significant. A simple wave theory is developed for the system of equations obtained. All self-similar discontinuous solutions and all continuous centered self-similar solutions of the system are obtained. The initial discontinuity decay problem is solved explicitly for compressible MHD equations in the shallow water approximation. It is shown that there exist five different configurations that provide a solution to the initial discontinuity decay problem. For each configuration, conditions are found that are necessary and sufficient for its implementation. Differences between incompressible and compressible cases are analyzed. In spite of the formal similarity between the solutions in the classical case of MHD flows of an incompressible and compressible fluids, the nonlinear dynamics described by the solutions are essentially different due to the difference in the expressions for the squared propagation velocity of weak perturbations. In addition, the solutions obtained describe new physical phenomena related to the dependence

Fluid Mechanics

NASA Astrophysics Data System (ADS)

Pnueli, David; Gutfinger, Chaim

1997-01-01

This text is intended for the study of fluid mechanics at an intermediate level. The presentation starts with basic concepts, in order to form a sound conceptual structure that can support engineering applications and encourage further learning. The presentation is exact, incorporating both the mathematics involved and the physics needed to understand the various phenomena in fluid mechanics. Where a didactical choice must be made between the two, the physics prevails. Throughout the book the authors have tried to reach a balance between exact presentation, intuitive grasp of new ideas, and creative applications of concepts. This approach is reflected in the examples presented in the text and in the exercises given at the end of each chapter. Subjects treated are hydrostatics, viscous flow, similitude and order of magnitude, creeping flow, potential flow, boundary layer flow, turbulent flow, compressible flow, and non-Newtonian flows. This book is ideal for advanced undergraduate students in mechanical, chemical, aerospace, and civil engineering. Solutions manual available.

A variational principle for compressible fluid mechanics: Discussion of the multi-dimensional theory

NASA Technical Reports Server (NTRS)

Prozan, R. J.

1982-01-01

The variational principle for compressible fluid mechanics previously introduced is extended to two dimensional flow. The analysis is stable, exactly conservative, adaptable to coarse or fine grids, and very fast. Solutions for two dimensional problems are included. The excellent behavior and results lend further credence to the variational concept and its applicability to the numerical analysis of complex flow fields.

Quasiconservation laws for compressible three-dimensional Navier-Stokes flow.

PubMed

Gibbon, J D; Holm, D D

2012-10-01

We formulate the quasi-Lagrangian fluid transport dynamics of mass density ρ and the projection q=ω·∇ρ of the vorticity ω onto the density gradient, as determined by the three-dimensional compressible Navier-Stokes equations for an ideal gas, although the results apply for an arbitrary equation of state. It turns out that the quasi-Lagrangian transport of q cannot cross a level set of ρ. That is, in this formulation, level sets of ρ (isopycnals) are impermeable to the transport of the projection q.

Energy transfer in compressible magnetohydrodynamic turbulence for isothermal self-gravitating fluids

NASA Astrophysics Data System (ADS)

Banerjee, Supratik; Kritsuk, Alexei G.

2018-02-01

Three-dimensional, compressible, magnetohydrodynamic turbulence of an isothermal, self-gravitating fluid is analyzed using two-point statistics in the asymptotic limit of large Reynolds numbers (both kinetic and magnetic). Following an alternative formulation proposed by Banerjee and Galtier [Phys. Rev. E 93, 033120 (2016), 10.1103/PhysRevE.93.033120; J. Phys. A: Math. Theor. 50, 015501 (2017), 10.1088/1751-8113/50/1/015501], an exact relation has been derived for the total energy transfer. This approach results in a simpler relation expressed entirely in terms of mixed second-order structure functions. The kinetic, thermodynamic, magnetic, and gravitational contributions to the energy transfer rate can be easily separated in the present form. By construction, the new formalism includes such additional effects as global rotation, the Hall term in the induction equation, etc. The analysis shows that solid-body rotation cannot alter the energy flux rate of compressible turbulence. However, the contribution of a uniform background magnetic field to the flux is shown to be nontrivial unlike in the incompressible case. Finally, the compressible, turbulent energy flux rate does not vanish completely due to simple alignments, which leads to a zero turbulent energy flux rate in the incompressible case.

Energy transfer in compressible magnetohydrodynamic turbulence for isothermal self-gravitating fluids.

PubMed

Banerjee, Supratik; Kritsuk, Alexei G

2018-02-01

Three-dimensional, compressible, magnetohydrodynamic turbulence of an isothermal, self-gravitating fluid is analyzed using two-point statistics in the asymptotic limit of large Reynolds numbers (both kinetic and magnetic). Following an alternative formulation proposed by Banerjee and Galtier [Phys. Rev. E 93, 033120 (2016)2470-004510.1103/PhysRevE.93.033120; J. Phys. A: Math. Theor. 50, 015501 (2017)1751-811310.1088/1751-8113/50/1/015501], an exact relation has been derived for the total energy transfer. This approach results in a simpler relation expressed entirely in terms of mixed second-order structure functions. The kinetic, thermodynamic, magnetic, and gravitational contributions to the energy transfer rate can be easily separated in the present form. By construction, the new formalism includes such additional effects as global rotation, the Hall term in the induction equation, etc. The analysis shows that solid-body rotation cannot alter the energy flux rate of compressible turbulence. However, the contribution of a uniform background magnetic field to the flux is shown to be nontrivial unlike in the incompressible case. Finally, the compressible, turbulent energy flux rate does not vanish completely due to simple alignments, which leads to a zero turbulent energy flux rate in the incompressible case.

Astrophysical fluid simulations of thermally ideal gases with non-constant adiabatic index: numerical implementation

NASA Astrophysics Data System (ADS)

Vaidya, B.; Mignone, A.; Bodo, G.; Massaglia, S.

2015-08-01

Context. An equation of state (EoS) is a relation between thermodynamic state variables and it is essential for closing the set of equations describing a fluid system. Although an ideal EoS with a constant adiabatic index Γ is the preferred choice owing to its simplistic implementation, many astrophysical fluid simulations may benefit from a more sophisticated treatment that can account for diverse chemical processes. Aims: In the present work we first review the basic thermodynamic principles of a gas mixture in terms of its thermal and caloric EoS by including effects like ionization, dissociation, and temperature dependent degrees of freedom such as molecular vibrations and rotations. The formulation is revisited in the context of plasmas that are either in equilibrium conditions (local thermodynamic- or collisional excitation-equilibria) or described by non-equilibrium chemistry coupled to optically thin radiative cooling. We then present a numerical implementation of thermally ideal gases obeying a more general caloric EoS with non-constant adiabatic index in Godunov-type numerical schemes. Methods: We discuss the necessary modifications to the Riemann solver and to the conversion between total energy and pressure (or vice versa) routinely invoked in Godunov-type schemes. We then present two different approaches for computing the EoS. The first employs root-finder methods and it is best suited for EoS in analytical form. The second is based on lookup tables and interpolation and results in a more computationally efficient approach, although care must be taken to ensure thermodynamic consistency. Results: A number of selected benchmarks demonstrate that the employment of a non-ideal EoS can lead to important differences in the solution when the temperature range is 500-104 K where dissociation and ionization occur. The implementation of selected EoS introduces additional computational costs although the employment of lookup table methods (when possible) can

Canonical fluid thermodynamics. [variational principles of stability for compressible adiabatic flow

NASA Technical Reports Server (NTRS)

Schmid, L. A.

1974-01-01

The space-time integral of the thermodynamic pressure plays in a certain sense the role of the thermodynamic potential for compressible adiabatic flow. The stability criterion can be converted into a variational minimum principle by requiring the molar free-enthalpy and temperature to be generalized velocities. In the fluid context, the definition of proper-time differentiation involves the fluid velocity expressed in terms of three particle identity parameters. The pressure function is then converted into a functional which is the Lagrangian density of the variational principle. Being also a minimum principle, the variational principle provides a means for comparing the relative stability of different flows. For boundary conditions with a high degree of symmetry, as in the case of a uniformly expanding spherical gas box, the most stable flow is a rectilinear flow for which the world-trajectory of each particle is a straight line. Since the behavior of the interior of a freely expanding cosmic cloud may be expected to be similar to that of the fluid in the spherical box of gas, this suggests that the cosmic principle is a consequence of the laws of thermodynamics, rather than just an ad hoc postulate.

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.

Compressible, multiphase semi-implicit method with moment of fluid interface representation

DOE PAGES

Jemison, Matthew; Sussman, Mark; Arienti, Marco

2014-09-16

A unified method for simulating multiphase flows using an exactly mass, momentum, and energy conserving Cell-Integrated Semi-Lagrangian advection algorithm is presented. The deforming material boundaries are represented using the moment-of-fluid method. Our new algorithm uses a semi-implicit pressure update scheme that asymptotically preserves the standard incompressible pressure projection method in the limit of infinite sound speed. The asymptotically preserving attribute makes the new method applicable to compressible and incompressible flows including stiff materials; enabling large time steps characteristic of incompressible flow algorithms rather than the small time steps required by explicit methods. Moreover, shocks are captured and material discontinuities aremore » tracked, without the aid of any approximate or exact Riemann solvers. As a result, wimulations of underwater explosions and fluid jetting in one, two, and three dimensions are presented which illustrate the effectiveness of the new algorithm at efficiently computing multiphase flows containing shock waves and material discontinuities with large “impedance mismatch.”« less

A volume-of-fluid method for simulation of compressible axisymmetric multi-material flow

NASA Astrophysics Data System (ADS)

de Niem, D.; Kührt, E.; Motschmann, U.

2007-02-01

A two-dimensional Eulerian hydrodynamic method for the numerical simulation of inviscid compressible axisymmetric multi-material flow in external force fields for the situation of pure fluids separated by macroscopic interfaces is presented. The method combines an implicit Lagrangian step with an explicit Eulerian advection step. Individual materials obey separate energy equations, fulfill general equations of state, and may possess different temperatures. Material volume is tracked using a piecewise linear volume-of-fluid method. An overshoot-free logically simple and economic material advection algorithm for cylinder coordinates is derived, in an algebraic formulation. New aspects arising in the case of more than two materials such as the material ordering strategy during transport are presented. One- and two-dimensional numerical examples are given.

An object-oriented and quadrilateral-mesh based solution adaptive algorithm for compressible multi-fluid flows

NASA Astrophysics Data System (ADS)

Zheng, H. W.; Shu, C.; Chew, Y. T.

2008-07-01

In this paper, an object-oriented and quadrilateral-mesh based solution adaptive algorithm for the simulation of compressible multi-fluid flows is presented. The HLLC scheme (Harten, Lax and van Leer approximate Riemann solver with the Contact wave restored) is extended to adaptively solve the compressible multi-fluid flows under complex geometry on unstructured mesh. It is also extended to the second-order of accuracy by using MUSCL extrapolation. The node, edge and cell are arranged in such an object-oriented manner that each of them inherits from a basic object. A home-made double link list is designed to manage these objects so that the inserting of new objects and removing of the existing objects (nodes, edges and cells) are independent of the number of objects and only of the complexity of O( 1). In addition, the cells with different levels are further stored in different lists. This avoids the recursive calculation of solution of mother (non-leaf) cells. Thus, high efficiency is obtained due to these features. Besides, as compared to other cell-edge adaptive methods, the separation of nodes would reduce the memory requirement of redundant nodes, especially in the cases where the level number is large or the space dimension is three. Five two-dimensional examples are used to examine its performance. These examples include vortex evolution problem, interface only problem under structured mesh and unstructured mesh, bubble explosion under the water, bubble-shock interaction, and shock-interface interaction inside the cylindrical vessel. Numerical results indicate that there is no oscillation of pressure and velocity across the interface and it is feasible to apply it to solve compressible multi-fluid flows with large density ratio (1000) and strong shock wave (the pressure ratio is 10,000) interaction with the interface.

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.

Symmetric flows for compressible heat-conducting fluids with temperature dependent viscosity coefficients

NASA Astrophysics Data System (ADS)

Wan, Ling; Wang, Tao

2017-06-01

We consider the Navier-Stokes equations for compressible heat-conducting ideal polytropic gases in a bounded annular domain when the viscosity and thermal conductivity coefficients are general smooth functions of temperature. A global-in-time, spherically or cylindrically symmetric, classical solution to the initial boundary value problem is shown to exist uniquely and converge exponentially to the constant state as the time tends to infinity under certain assumptions on the initial data and the adiabatic exponent γ. The initial data can be large if γ is sufficiently close to 1. These results are of Nishida-Smoller type and extend the work (Liu et al. (2014) [16]) restricted to the one-dimensional flows.

Steady Secondary Flows Generated by Periodic Compression and Expansion of an Ideal Gas in a Pulse Tube

NASA Technical Reports Server (NTRS)

Lee, Jeffrey M.

1999-01-01

This study establishes a consistent set of differential equations for use in describing the steady secondary flows generated by periodic compression and expansion of an ideal gas in pulse tubes. Also considered is heat transfer between the gas and the tube wall of finite thickness. A small-amplitude series expansion solution in the inverse Strouhal number is proposed for the two-dimensional axisymmetric mass, momentum and energy equations. The anelastic approach applies when shock and acoustic energies are small compared with the energy needed to compress and expand the gas. An analytic solution to the ordered series is obtained in the strong temperature limit where the zeroth-order temperature is constant. The solution shows steady velocities increase linearly for small Valensi number and can be of order I for large Valensi number. A conversion of steady work flow to heat flow occurs whenever temperature, velocity or phase angle gradients are present. Steady enthalpy flow is reduced by heat transfer and is scaled by the Prandtl times Valensi numbers. Particle velocities from a smoke-wire experiment were compared with predictions for the basic and orifice pulse tube configurations. The theory accurately predicted the observed steady streaming.

Fluid Absorption and Release of Nonwovens and their Response to Compression

NASA Astrophysics Data System (ADS)

Bateny, Fatemeh

Fluid handling is a key property in one of the major nonwoven applications in absorbent product such as wipes, hygiene products, and baby diapers. These products are subjected to various levels of compression in real-use. The aim of this study was to investigate the liquid absorption and release properties of nonwovens to establish the absorption structure-property relationship at various compression levels. A comprehensive methodology, considering various flow directions, was employed to establish the relationship by decoupling the effect of structural parameters and material properties in two phases of this study respectively. In the first phase, the mechanism of absorption by pore structure was investigated through considering various fiber cross-sectional size and shape, as well as heterogeneous layered structures having a pore size reduction and expansion. In the second phase, the mechanism of absorption by fiber and consequent swelling was evaluated in view of fluid diffusion into the rayon fibers in samples having different percentages of PET fiber (non-absorbent) and rayon fiber (absorbent). The analysis of absorption and release properties through the entire dissertation was based on the pore characteristics of the nonwovens by measuring the average pore sizes, pore size distribution, and solidity. The investigation revealed that the absorption and release properties of nonwovens are governed by their pore characteristics. In homogeneous non-layered nonwoven fabrics, maximum absorption is mainly governed by the available pore volume. Absorbency rate is determined according to pore size and the maximum rate of absorption is achieved at a specific range of pore sizes. This indicates that an in-depth understanding of the absorption and release properties brings about valuable information for the absorbent product engineering.

Process of establishing a plane-wave system on ice cover over a dipole moving uniformly in an ideal fluid column

NASA Astrophysics Data System (ADS)

Il'ichev, A. T.; Savin, A. S.

2017-12-01

We consider a planar evolution problem for perturbations of the ice cover by a dipole starting its uniform rectilinear horizontal motion in a column of an initially stationary fluid. Using asymptotic Fourier analysis, we show that at supercritical velocities, waves of two types form on the water-ice interface. We describe the process of establishing these waves during the dipole motion. We assume that the fluid is ideal and incompressible and its motion is potential. The ice cover is modeled by the Kirchhoff-Love plate.

Isoscalar compression modes within fluid dynamic approach

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

Kolomietz, V. M.; Cyclotron Institute, Texas A and M University, College Station, Texas 77843-3366; Shlomo, S.

2000-06-01

We study the nuclear isoscalar monopole and dipole compression modes in nuclei within the fluid dynamic approach (FDA) with and without the effect of relaxation. For a wide region of the medium and heavy nuclei, the FDA predicts that the isoscalar giant monopole resonance (ISGMR) and the isoscalar giant dipole resonance (ISGDR) exhaust about 90% of the corresponding model-independent sum rules. In the case of neglecting the effect of relaxation, the FDA, when adjusted to reproduce the centroid energy E0 of the ISGMR, results with centroid energy E1 of the ISGDR which is in agreement with the predictions of themore » self-consistent Hartree-Fock random-phase approximation calculations and the scaling model but significantly larger than the experimental value. We also show that the FDA leads to the correct hydrodynamic limit for the ratio (E1/E0){sub FDA}. We find that the ratio (E1/E0){sub FDA} depends on the relaxation time and approaches the preliminary experimental value (E1/E0){sub exp}=1.5{+-}0.1 in a short relaxation time limit. (c) 2000 The American Physical Society.« less

Design Space Approach in Optimization of Fluid Bed Granulation and Tablets Compression Process

PubMed Central

Djuriš, Jelena; Medarević, Djordje; Krstić, Marko; Vasiljević, Ivana; Mašić, Ivana; Ibrić, Svetlana

2012-01-01

The aim of this study was to optimize fluid bed granulation and tablets compression processes using design space approach. Type of diluent, binder concentration, temperature during mixing, granulation and drying, spray rate, and atomization pressure were recognized as critical formulation and process parameters. They were varied in the first set of experiments in order to estimate their influences on critical quality attributes, that is, granules characteristics (size distribution, flowability, bulk density, tapped density, Carr's index, Hausner's ratio, and moisture content) using Plackett-Burman experimental design. Type of diluent and atomization pressure were selected as the most important parameters. In the second set of experiments, design space for process parameters (atomization pressure and compression force) and its influence on tablets characteristics was developed. Percent of paracetamol released and tablets hardness were determined as critical quality attributes. Artificial neural networks (ANNs) were applied in order to determine design space. ANNs models showed that atomization pressure influences mostly on the dissolution profile, whereas compression force affects mainly the tablets hardness. Based on the obtained ANNs models, it is possible to predict tablet hardness and paracetamol release profile for any combination of analyzed factors. PMID:22919295

Self-Similar Apical Sharpening of an Ideal Perfecting Conducting Fluid Subject to Maxwell Stresses

NASA Astrophysics Data System (ADS)

Zhou, Chengzhe; Troian, Sandra M.

2016-11-01

We examine the apical behavior of an ideal, perfectly conducting incompressible fluid surrounded by vacuum in circumstances where the capillary, Maxwell and inertial forces contribute to formation of a liquid cone. A previous model based on potential flow describes a family of self-similar solutions with conic cusps whose interior angles approach the Taylor cone angle. These solutions were obtained by matching powers of the leading order terms in the velocity and electric field potential to the asymptotic form dictated by a stationary cone shape. In re-examining this earlier work, we have found a more important, neglected leading order term in the velocity and field potentials, which satisfies the governing, interfacial and far-field conditions as well. This term allows for the development of additional self-similar, sharpening apical shapes, including time reversed solutions for conic tip recoil after fluid ejection. We outline the boundary-element technique for solving the exact similarity solutions, which have parametric dependence on the far-field conditions, and discuss consequences of our findings.

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.

University Students Explaining Adiabatic Compression of an Ideal Gas—A New Phenomenon in Introductory Thermal Physics

NASA Astrophysics Data System (ADS)

Leinonen, Risto; Asikainen, Mervi A.; Hirvonen, Pekka E.

2012-12-01

This study focuses on second-year university students' explanations and reasoning related to adiabatic compression of an ideal gas. The phenomenon was new to the students, but it was one which they should have been capable of explaining using their previous upper secondary school knowledge. The students' explanations and reasoning were investigated with the aid of paper and pencil tests ( n = 86) and semi-structured interviews ( n = 5) at the start of a thermal physics course at the University of Eastern Finland. The paper and pencil test revealed that the students had difficulties in applying content taught during earlier education in a new context: only a few of them were able to produce a correct explanation for the phenomenon. A majority of the students used either explanations with invalid but physically correct models, such as the ideal gas law or a microscopic model, or erroneous dependencies between quantities. The results also indicated that students had problems in seeing deficiencies or inconsistencies in their reasoning, in both test and interview situations. We suggest in our conclusion that the contents of upper secondary school thermal physics courses should be carefully examined to locate the best emphases for different laws, principles, concepts, and models. In particular, the limitations of models should be made explicit in teaching and students should be guided towards critical scientific thinking, including metaconceptual awareness.

Entropy stable high order discontinuous Galerkin methods for ideal compressible MHD on structured meshes

NASA Astrophysics Data System (ADS)

Liu, Yong; Shu, Chi-Wang; Zhang, Mengping

2018-02-01

We present a discontinuous Galerkin (DG) scheme with suitable quadrature rules [15] for ideal compressible magnetohydrodynamic (MHD) equations on structural meshes. The semi-discrete scheme is analyzed to be entropy stable by using the symmetrizable version of the equations as introduced by Godunov [32], the entropy stable DG framework with suitable quadrature rules [15], the entropy conservative flux in [14] inside each cell and the entropy dissipative approximate Godunov type numerical flux at cell interfaces to make the scheme entropy stable. The main difficulty in the generalization of the results in [15] is the appearance of the non-conservative "source terms" added in the modified MHD model introduced by Godunov [32], which do not exist in the general hyperbolic system studied in [15]. Special care must be taken to discretize these "source terms" adequately so that the resulting DG scheme satisfies entropy stability. Total variation diminishing / bounded (TVD/TVB) limiters and bound-preserving limiters are applied to control spurious oscillations. We demonstrate the accuracy and robustness of this new scheme on standard MHD examples.

Global Well-Posedness and Decay Rates of Strong Solutions to a Non-Conservative Compressible Two-Fluid Model

NASA Astrophysics Data System (ADS)

Evje, Steinar; Wang, Wenjun; Wen, Huanyao

2016-09-01

In this paper, we consider a compressible two-fluid model with constant viscosity coefficients and unequal pressure functions {P^+neq P^-}. As mentioned in the seminal work by Bresch, Desjardins, et al. (Arch Rational Mech Anal 196:599-629, 2010) for the compressible two-fluid model, where {P^+=P^-} (common pressure) is used and capillarity effects are accounted for in terms of a third-order derivative of density, the case of constant viscosity coefficients cannot be handled in their settings. Besides, their analysis relies on a special choice for the density-dependent viscosity [refer also to another reference (Commun Math Phys 309:737-755, 2012) by Bresch, Huang and Li for a study of the same model in one dimension but without capillarity effects]. In this work, we obtain the global solution and its optimal decay rate (in time) with constant viscosity coefficients and some smallness assumptions. In particular, capillary pressure is taken into account in the sense that {Δ P=P^+ - P^-=fneq 0} where the difference function {f} is assumed to be a strictly decreasing function near the equilibrium relative to the fluid corresponding to {P^-}. This assumption plays an key role in the analysis and appears to have an essential stabilization effect on the model in question.

Compressible Analysis of Bénard Convection of Magneto Rotatory Couple-Stress Fluid

NASA Astrophysics Data System (ADS)

Mehta, C. B.; Singh, M.

2018-02-01

Thermal Instability (Benard's Convection) in the presence of uniform rotation and uniform magnetic field (separately) is studied. Using the linearized stability theory and normal mode analyses the dispersion relation is obtained in each case. In the case of rotatory Benard's stationary convection compressibility and rotation postpone the onset of convection whereas the couple-stress have duel character onset of convection depending on rotation parameter. While in the absence of rotation couple-stress always postpones the onset of convection. On the other hand, magnetic field on thermal instability problem on couple-stress fluid for stationary convection couple-stress parameter and magnetic field postpones the onset of convection. The effect of compressibility also postpones the onset of convection in both cases as rotation and magnetic field. Graphs have been plotted by giving numerical values to the parameters to depict the stationary characteristics. Further, the magnetic field and rotation are found to introduce oscillatory modes which were non-existent in their absence and then the principle of exchange of stability is valid. The sufficient conditions for non-existence of overstability are also obtained.

An Immersed Boundary Method for Solving the Compressible Navier-Stokes Equations with Fluid Structure Interaction

NASA Technical Reports Server (NTRS)

Brehm, Christoph; Barad, Michael F.; Kiris, Cetin C.

2016-01-01

An immersed boundary method for the compressible Navier-Stokes equation and the additional infrastructure that is needed to solve moving boundary problems and fully coupled fluid-structure interaction is described. All the methods described in this paper were implemented in NASA's LAVA solver framework. The underlying immersed boundary method is based on the locally stabilized immersed boundary method that was previously introduced by the authors. In the present paper this method is extended to account for all aspects that are involved for fluid structure interaction simulations, such as fast geometry queries and stencil computations, the treatment of freshly cleared cells, and the coupling of the computational fluid dynamics solver with a linear structural finite element method. The current approach is validated for moving boundary problems with prescribed body motion and fully coupled fluid structure interaction problems in 2D and 3D. As part of the validation procedure, results from the second AIAA aeroelastic prediction workshop are also presented. The current paper is regarded as a proof of concept study, while more advanced methods for fluid structure interaction are currently being investigated, such as geometric and material nonlinearities, and advanced coupling approaches.

Ideal Particle Sizes for Inhaled Steroids Targeting Vocal Granulomas: Preliminary Study Using Computational Fluid Dynamics.

PubMed

Perkins, Elizabeth L; Basu, Saikat; Garcia, Guilherme J M; Buckmire, Robert A; Shah, Rupali N; Kimbell, Julia S

2018-03-01

Objectives Vocal fold granulomas are benign lesions of the larynx commonly caused by gastroesophageal reflux, intubation, and phonotrauma. Current medical therapy includes inhaled corticosteroids to target inflammation that leads to granuloma formation. Particle sizes of commonly prescribed inhalers range over 1 to 4 µm. The study objective was to use computational fluid dynamics to investigate deposition patterns over a range of particle sizes of inhaled corticosteroids targeting the larynx and vocal fold granulomas. Study Design Retrospective, case-specific computational study. Setting Tertiary academic center. Subjects/Methods A 3-dimensional anatomically realistic computational model of a normal adult airway from mouth to trachea was constructed from 3 computed tomography scans. Virtual granulomas of varying sizes and positions along the vocal fold were incorporated into the base model. Assuming steady-state, inspiratory, turbulent airflow at 30 L/min, computational fluid dynamics was used to simulate respiratory transport and deposition of inhaled corticosteroid particles ranging over 1 to 20 µm. Results Laryngeal deposition in the base model peaked for particle sizes 8 to 10 µm (2.8%-3.5%). Ideal sizes ranged over 6 to 10, 7 to 13, and 7 to 14 µm for small, medium, and large granuloma sizes, respectively. Glottic deposition was maximal at 10.8% for 9-µm-sized particles for the large posterior granuloma, 3 times the normal model (3.5%). Conclusion As the virtual granuloma size increased and the location became more posterior, glottic deposition and ideal particle size generally increased. This preliminary study suggests that inhalers with larger particle sizes, such as fluticasone propionate dry-powder inhaler, may improve laryngeal drug deposition. Most commercially available inhalers have smaller particles than suggested here.

Lagrangian fluid description with simple applications in compressible plasma and gas dynamics

NASA Astrophysics Data System (ADS)

Schamel, Hans

2004-03-01

The Lagrangian fluid description, in which the dynamics of fluids is formulated in terms of trajectories of fluid elements, not only presents an alternative to the more common Eulerian description but has its own merits and advantages. This aspect, which seems to be not fully explored yet, is getting increasing attention in fluid dynamics and related areas as Lagrangian codes and experimental techniques are developed utilizing the Lagrangian point of view with the ultimate goal of a deeper understanding of flow dynamics. In this tutorial review we report on recent progress made in the analysis of compressible, more or less perfect flows such as plasmas and dilute gases. The equations of motion are exploited to get further insight into the formation and evolution of coherent structures, which often exhibit a singular or collapse type behavior occurring in finite time. It is argued that this technique of solution has a broad applicability due to the simplicity and generality of equations used. The focus is on four different topics, the physics of which being governed by simple fluid equations subject to initial and/or boundary conditions. Whenever possible also experimental results are mentioned. In the expansion of a semi-infinite plasma into a vacuum the energetic ion peak propagating supersonically towards the vacuum-as seen in laboratory experiments-is interpreted by means of the Lagrangian fluid description as a relic of a wave breaking scenario of the corresponding inviscid ion dynamics. The inclusion of viscosity is shown numerically to stabilize the associated density collapse giving rise to a well defined fast ion peak reminiscent of adhesive matter. In purely convection driven flows the Lagrangian flow velocity is given by its initial value and hence the Lagrangian velocity gradient tensor can be evaluated accurately to find out the appearance of singularities in density and vorticity and the emergence of new structures such as wavelets in one-dimension (1D

Using a tracer technique to identify the extent of non-ideal flows in the continuous mixing of non-Newtonian fluids

NASA Astrophysics Data System (ADS)

Patel, D.; Ein-Mozaffari, F.; Mehrvar, M.

2013-05-01

The identification of non-ideal flows in a continuous-flow mixing of non-Newtonian fluids is a challenging task for various chemical industries: plastic manufacturing, water and wastewater treatment, and pulp and paper manufacturing. Non-ideal flows such as channelling, recirculation, and dead zones significantly affect the performance of continuous-flow mixing systems. Therefore, the main objective of this paper was to develop an identification protocol to measure non-ideal flows in the continuous-flow mixing system. The extent of non-ideal flows was quantified using a dynamic model that incorporated channelling, recirculation, and dead volume in the mixing vessel. To estimate the dynamic model parameters, the system was excited using a frequency-modulated random binary input by injecting the saline solution (as a tracer) into the fresh feed stream prior to being pumped into the mixing vessel. The injection of the tracer was controlled by a computer-controlled on-off solenoid valve. Using the trace technique, the extent of channelling and the effective mixed volume were successfully determined and used as mixing quality criteria. Such identification procedures can be applied at various areas of chemical engineering in order to improve the mixing quality.

Magnetic power piston fluid compressor

NASA Technical Reports Server (NTRS)

Gasser, Max G. (Inventor)

1994-01-01

A compressor with no moving parts in the traditional sense having a housing having an inlet end allowing a low pressure fluid to enter and an outlet end allowing a high pressure fluid to exit is described. Within the compressor housing is at least one compression stage to increase the pressure of the fluid within the housing. The compression stage has a quantity of magnetic powder within the housing, is supported by a screen that allows passage of the fluid, and a coil for selectively providing a magnetic field across the magnetic powder such that when the magnetic field is not present the individual particles of the powder are separated allowing the fluid to flow through the powder and when the magnetic field is present the individual particles of the powder pack together causing the powder mass to expand preventing the fluid from flowing through the powder and causing a pressure pulse to compress the fluid.

On controlling nonlinear dissipation in high order filter methods for ideal and non-ideal MHD

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sjogreen, B.

2004-01-01

The newly developed adaptive numerical dissipation control in spatially high order filter schemes for the compressible Euler and Navier-Stokes equations has been recently extended to the ideal and non-ideal magnetohydrodynamics (MHD) equations. These filter schemes are applicable to complex unsteady MHD high-speed shock/shear/turbulence problems. They also provide a natural and efficient way for the minimization of Div(B) numerical error. The adaptive numerical dissipation mechanism consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. The numerical dissipation considered consists of high order linear dissipation for the suppression of high frequency oscillation and the nonlinear dissipative portion of high-resolution shock-capturing methods for discontinuity capturing. The applicable nonlinear dissipative portion of high-resolution shock-capturing methods is very general. The objective of this paper is to investigate the performance of three commonly used types of nonlinear numerical dissipation for both the ideal and non-ideal MHD.

An oscillation free shock-capturing method for compressible van der Waals supercritical fluid flows

DOE PAGES

Pantano, C.; Saurel, R.; Schmitt, T.

2017-02-01

Numerical solutions of the Euler equations using real gas equations of state (EOS) often exhibit serious inaccuracies. The focus here is the van der Waals EOS and its variants (often used in supercritical fluid computations). The problems are not related to a lack of convexity of the EOS since the EOS are considered in their domain of convexity at any mesh point and at any time. The difficulties appear as soon as a density discontinuity is present with the rest of the fluid in mechanical equilibrium and typically result in spurious pressure and velocity oscillations. This is reminiscent of well-knownmore » pressure oscillations occurring with ideal gas mixtures when a mass fraction discontinuity is present, which can be interpreted as a discontinuity in the EOS parameters. We are concerned with pressure oscillations that appear just for a single fluid each time a density discontinuity is present. As a result, the combination of density in a nonlinear fashion in the EOS with diffusion by the numerical method results in violation of mechanical equilibrium conditions which are not easy to eliminate, even under grid refinement.« less

An oscillation free shock-capturing method for compressible van der Waals supercritical fluid flows

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

Pantano, C.; Saurel, R.; Schmitt, T.

Numerical solutions of the Euler equations using real gas equations of state (EOS) often exhibit serious inaccuracies. The focus here is the van der Waals EOS and its variants (often used in supercritical fluid computations). The problems are not related to a lack of convexity of the EOS since the EOS are considered in their domain of convexity at any mesh point and at any time. The difficulties appear as soon as a density discontinuity is present with the rest of the fluid in mechanical equilibrium and typically result in spurious pressure and velocity oscillations. This is reminiscent of well-knownmore » pressure oscillations occurring with ideal gas mixtures when a mass fraction discontinuity is present, which can be interpreted as a discontinuity in the EOS parameters. We are concerned with pressure oscillations that appear just for a single fluid each time a density discontinuity is present. As a result, the combination of density in a nonlinear fashion in the EOS with diffusion by the numerical method results in violation of mechanical equilibrium conditions which are not easy to eliminate, even under grid refinement.« less

Biot-Savart helicity versus physical helicity: A topological description of ideal flows

NASA Astrophysics Data System (ADS)

Sahihi, Taliya; Eshraghi, Homayoon

2014-08-01

For an isentropic (thus compressible) flow, fluid trajectories are considered as orbits of a family of one parameter, smooth, orientation-preserving, and nonsingular diffeomorphisms on a compact and smooth-boundary domain in the Euclidian 3-space which necessarily preserve a finite measure, later interpreted as the fluid mass. Under such diffeomorphisms the Biot-Savart helicity of the pushforward of a divergence-free and tangent to the boundary vector field is proved to be conserved and since these circumstances present an isentropic flow, the conservation of the "Biot-Savart helicity" is established for such flows. On the other hand, the well known helicity conservation in ideal flows which here we call it "physical helicity" is found to be an independent constant with respect to the Biot-Savart helicity. The difference between these two helicities reflects some topological features of the domain as well as the velocity and vorticity fields which is discussed and is shown for simply connected domains the two helicities coincide. The energy variation of the vorticity field is shown to be formally the same as for the incompressible flow obtained before. For fluid domains consisting of several disjoint solid tori, at each time, the harmonic knot subspace of smooth vector fields on the fluid domain is found to have two independent base sets with a special type of orthogonality between these two bases by which a topological description of the vortex and velocity fields depending on the helicity difference is achieved since this difference is shown to depend only on the harmonic knot parts of velocity, vorticity, and its Biot-Savart vector field. For an ideal magnetohydrodynamics (MHD) flow three independent constant helicities are reviewed while the helicity of magnetic potential is generalized for non-simply connected domains by inserting a special harmonic knot field in the dynamics of the magnetic potential. It is proved that the harmonic knot part of the vorticity

Universal penetration test apparatus with fluid penetration sensor

DOEpatents

Johnson, Phillip W.; Stampfer, Joseph F.; Bradley, Orvil D.

1999-01-01

A universal penetration test apparatus for measuring resistance of a material to a challenge fluid. The apparatus includes a pad saturated with the challenge fluid. The apparatus includes a compression assembly for compressing the material between the pad and a compression member. The apparatus also includes a sensor mechanism for automatically detecting when the challenge fluid penetrates the material.

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...

Flow-through compression cell for small-angle and ultra-small-angle neutron scattering measurements

NASA Astrophysics Data System (ADS)

Hjelm, Rex P.; Taylor, Mark A.; Frash, Luke P.; Hawley, Marilyn E.; Ding, Mei; Xu, Hongwu; Barker, John; Olds, Daniel; Heath, Jason; Dewers, Thomas

2018-05-01

In situ measurements of geological materials under compression and with hydrostatic fluid pressure are important in understanding their behavior under field conditions, which in turn provides critical information for application-driven research. In particular, understanding the role of nano- to micro-scale porosity in the subsurface liquid and gas flow is critical for the high-fidelity characterization of the transport and more efficient extraction of the associated energy resources. In other applications, where parts are produced by the consolidation of powders by compression, the resulting porosity and crystallite orientation (texture) may affect its in-use characteristics. Small-angle neutron scattering (SANS) and ultra SANS are ideal probes for characterization of these porous structures over the nano to micro length scales. Here we show the design, realization, and performance of a novel neutron scattering sample environment, a specially designed compression cell, which provides compressive stress and hydrostatic pressures with effective stress up to 60 MPa, using the neutron beam to probe the effects of stress vectors parallel to the neutron beam. We demonstrate that the neutron optics is suitable for the experimental objectives and that the system is highly stable to the stress and pressure conditions of the measurements.

Universal penetration test apparatus with fluid penetration sensor

DOEpatents

Johnson, P.W.; Stampfer, J.F.; Bradley, O.D.

1999-02-02

A universal penetration test apparatus is described for measuring resistance of a material to a challenge fluid. The apparatus includes a pad saturated with the challenge fluid. The apparatus includes a compression assembly for compressing the material between the pad and a compression member. The apparatus also includes a sensor mechanism for automatically detecting when the challenge fluid penetrates the material. 23 figs.

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.

Energy recovery during expansion of compressed gas using power plant low-quality heat sources

DOEpatents

Ochs, Thomas L [Albany, OR; O'Connor, William K [Lebanon, OR

2006-03-07

A method of recovering energy from a cool compressed gas, compressed liquid, vapor, or supercritical fluid is disclosed which includes incrementally expanding the compressed gas, compressed liquid, vapor, or supercritical fluid through a plurality of expansion engines and heating the gas, vapor, compressed liquid, or supercritical fluid entering at least one of the expansion engines with a low quality heat source. Expansion engines such as turbines and multiple expansions with heating are disclosed.

Lagrangian averages, averaged Lagrangians, and the mean effects of fluctuations in fluid dynamics.

PubMed

Holm, Darryl D.

2002-06-01

We begin by placing the generalized Lagrangian mean (GLM) equations for a compressible adiabatic fluid into the Euler-Poincare (EP) variational framework of fluid dynamics, for an averaged Lagrangian. This is the Lagrangian averaged Euler-Poincare (LAEP) theorem. Next, we derive a set of approximate small amplitude GLM equations (glm equations) at second order in the fluctuating displacement of a Lagrangian trajectory from its mean position. These equations express the linear and nonlinear back-reaction effects on the Eulerian mean fluid quantities by the fluctuating displacements of the Lagrangian trajectories in terms of their Eulerian second moments. The derivation of the glm equations uses the linearized relations between Eulerian and Lagrangian fluctuations, in the tradition of Lagrangian stability analysis for fluids. The glm derivation also uses the method of averaged Lagrangians, in the tradition of wave, mean flow interaction. Next, the new glm EP motion equations for incompressible ideal fluids are compared with the Euler-alpha turbulence closure equations. An alpha model is a GLM (or glm) fluid theory with a Taylor hypothesis closure. Such closures are based on the linearized fluctuation relations that determine the dynamics of the Lagrangian statistical quantities in the Euler-alpha equations. Thus, by using the LAEP theorem, we bridge between the GLM equations and the Euler-alpha closure equations, through the small-amplitude glm approximation in the EP variational framework. We conclude by highlighting a new application of the GLM, glm, and alpha-model results for Lagrangian averaged ideal magnetohydrodynamics. (c) 2002 American Institute of Physics.

Compressible viscous flows generated by oscillating flexible cylinders

NASA Astrophysics Data System (ADS)

Van Eysden, Cornelis A.; Sader, John E.

2009-01-01

The fluid dynamics of oscillating elastic beams underpin the operation of many modern technological devices ranging from micromechanical sensors to the atomic force microscope. While viscous effects are widely acknowledged to have a strong influence on these dynamics, fluid compressibility is commonly neglected. Here, we theoretically study the three-dimensional flow fields that are generated by the motion of flexible cylinders immersed in viscous compressible fluids and discuss the implications of compressibility in practice. We consider cylinders of circular cross section and flat blades of zero thickness that are executing flexural and torsional oscillations of arbitrary wave number. Exact analytical solutions are derived for these flow fields and their resulting hydrodynamic loads.

Potential Flow Model for Compressible Stratified Rayleigh-Taylor Instability

NASA Astrophysics Data System (ADS)

Rydquist, Grant; Reckinger, Scott; Owkes, Mark; Wieland, Scott

2017-11-01

The Rayleigh-Taylor Instability (RTI) is an instability that occurs when a heavy fluid lies on top of a lighter fluid in a gravitational field, or a gravity-like acceleration. It occurs in many fluid flows of a highly compressive nature. In this study potential flow analysis (PFA) is used to model the early stages of RTI growth for compressible fluids. In the localized region near the bubble tip, the effects of vorticity are negligible, so PFA is applicable, as opposed to later stages where the induced velocity due to vortices generated from the growth of the instability dominate the flow. The incompressible PFA is extended for compressibility effects by applying the growth rate and the associated perturbation spatial decay from compressible linear stability theory. The PFA model predicts theoretical values for a bubble terminal velocity for single-mode compressible RTI, dependent upon the Atwood (A) and Mach (M) numbers, which is a parameter that measures both the strength of the stratification and intrinsic compressibility. The theoretical bubble terminal velocities are compared against numerical simulations. The PFA model correctly predicts the M dependence at high A, but the model must be further extended to include additional physics to capture the behavior at low A. Undergraduate Scholars Program - Montana State University.

A FLUID SORBENT RECYCLING DEVICE FOR INDUSTRIAL FLUID USERS

EPA Science Inventory

A roller compression Extractor® that extracts fluids from reusable sorbent pads was evaluated as a method of waste reduction. The extraction device, evaluated for industrial fluid users in New Jersey, was found to be effective in recycling unpleated sorbent pads, especially ...

Amniotic fluid stem cells: an ideal resource for therapeutic application in bone tissue engineering.

PubMed

Pantalone, A; Antonucci, I; Guelfi, M; Pantalone, P; Usuelli, F G; Stuppia, L; Salini, V

2016-07-01

Skeletal diseases, both degenerative and secondary to trauma, infections or tumors, represent an ideal target for regenerative medicine and in the last years, stem cells have been considered as good candidates for in vitro and in vivo bone regeneration. To date, several stem cell sources, such as adult mesenchymal stem cells, embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs), have shown significant osteogenic potential. In this narrative review, we analyze the possible advantages of the use of AFSCs in the treatment of skeletal diseases, especially through the application of tissue engineering and biomaterials. Among the different sources of stem cells, great attention has been recently devoted to amniotic fluid-derived stem cells (AFSC) characterized by high renewal capacity and ability to differentiate along several different lineages. Due to these features, AFSCs represent an interesting model for regenerative medicine, also considering their low immunogenicity and the absence of tumor formation after transplantation in nude mice.

A computational fluid dynamics (CFD) study of WEB-treated aneurysms: Can CFD predict WEB "compression" during follow-up?

PubMed

Caroff, Jildaz; Mihalea, Cristian; Da Ros, Valerio; Yagi, Takanobu; Iacobucci, Marta; Ikka, Léon; Moret, Jacques; Spelle, Laurent

2017-07-01

Recent reports have revealed a worsening of aneurysm occlusion between WEB treatment baseline and angiographic follow-up due to "compression" of the device. We utilized computational fluid dynamics (CFD) in order to determine whether the underlying mechanism of this worsening is flow related. We included data from all consecutive patients treated in our institution with a WEB for unruptured aneurysms located either at the middle cerebral artery or basilar tip. The CFD study was performed using pre-operative 3D rotational angiography. From digital subtraction follow-up angiographies patients were dichotomized into two groups: one with WEB "compression" and one without. We performed statistical analyses to determine a potential correlation between WEB compression and CFD inflow ratio. Between July 2012 and June 2015, a total of 22 unruptured middle cerebral artery or basilar tip aneurysms were treated with a WEB device in our department. Three patients were excluded from the analysis and the mean follow-up period was 17months. Eleven WEBs presented "compression" during follow-up. Interestingly, device "compression" was statistically correlated to the CFD inflow ratio (P=0.018), although not to aneurysm volume, aspect ratio or neck size. The mechanisms underlying the worsening of aneurysm occlusion in WEB-treated patients due to device compression are most likely complex as well as multifactorial. However, it is apparent from our pilot study that a high arterial inflow is, at least, partially involved. Further theoretical and animal research studies are needed to increase our understanding of this phenomenon. Copyright © 2017 Elsevier Masson SAS. All rights reserved.

Compressible convection in geophysical fluids: comparison of anelastic, anelastic liquid and full numerical simulations

NASA Astrophysics Data System (ADS)

Curbelo, Jezabel; Alboussiere, Thierry; Labrosse, Stephane; Dubuffet, Fabien; Ricard, Yanick

2015-11-01

In this talk we describe the numerical method implemented to study convection in a fully compressible two-dimensional model, which may be reduced to the different simplifications such as the anelastic approximation and the anelastic liquid approximation. Various equations of state are considered, from the ideal gas equation to equations related to liquid or solid condensed matter. We are particularly interested in the total value and spatial distribution of viscous dissipation. We analyze the solutions obtained with each approximation in a wide range of dimensionless parameters and compare the domain of validity of each of them. The authors are grateful to the LABEX Lyon Institute of Origins (ANR-10-LABX-0066) of the Universite de Lyon for its financial support ``Investissements d'Avenir'' (ANR-11-IDEX-0007) of the French government operated by the National Research Agency (ANR).

The importance of robust error control in data compression applications

NASA Technical Reports Server (NTRS)

Woolley, S. I.

1993-01-01

Data compression has become an increasingly popular option as advances in information technology have placed further demands on data storage capabilities. With compression ratios as high as 100:1 the benefits are clear; however, the inherent intolerance of many compression formats to error events should be given careful consideration. If we consider that efficiently compressed data will ideally contain no redundancy, then the introduction of a channel error must result in a change of understanding from that of the original source. While the prefix property of codes such as Huffman enables resynchronisation, this is not sufficient to arrest propagating errors in an adaptive environment. Arithmetic, Lempel-Ziv, discrete cosine transform (DCT) and fractal methods are similarly prone to error propagating behaviors. It is, therefore, essential that compression implementations provide sufficient combatant error control in order to maintain data integrity. Ideally, this control should be derived from a full understanding of the prevailing error mechanisms and their interaction with both the system configuration and the compression schemes in use.

A Generalized Fluid System Simulation Program to Model Flow Distribution in Fluid Networks

NASA Technical Reports Server (NTRS)

Majumdar, Alok; Bailey, John W.; Schallhorn, Paul; Steadman, Todd

1998-01-01

This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.

Self-charging metering and dispensing device for fluids

NASA Technical Reports Server (NTRS)

Hooper, S. L.; Setzer, D. (Inventor)

1984-01-01

A self-metering and dispensing device for fluids obtained from a pressurized fluid supply is discussed. Tubing and valving means permit the introduction of fluid into and discharge from a closed cylindrical reservoir. The reservoir contains a slideably disposed piston co-acting with a coil compression spring, with piston travel determining the amount of fluid in the reservoir. Once the determined amount of fluid is introduced into the reservoir, the fluid is discharged by the force of the coil compression spring acting upon the piston.

Hybrid fluid-particle simulation of whistler-mode waves in a compressed dipole magnetic field: Implications for dayside high-latitude chorus

NASA Astrophysics Data System (ADS)

da Silva, C. L.; Wu, S.; Denton, R. E.; Hudson, M. K.; Millan, R. M.

2017-01-01

In this work we present a methodology for simulating whistler-mode waves self-consistently generated by electron temperature anisotropy in the inner magnetosphere. We present simulation results using a hybrid fluid/particle-in-cell code that treats the hot, anisotropic (i.e., ring current) electron population as particles and the background (i.e., the cold and inertialess) electrons as fluid. Since the hot electrons are only a small fraction of the total population, warm (and isotropic) particle electrons are added to the simulation to increase the fraction of particles with mass, providing a more accurate characterization of the wave dispersion relation. Ions are treated as a fixed background of positive charge density. The plasma transport equations are coupled to Maxwell's equations and solved in a meridional plane (a 2-D simulation with 3-D fields). We use a curvilinear coordinate system that follows the topological curvature of Earth's geomagnetic field lines, based on an analytic expression for a compressed dipole magnetic field. Hence, we are able to simulate whistler wave generation at dawn (pure dipole field lines) and dayside (compressed dipole) by simply adjusting one scalar quantity. We demonstrate how, on the dayside, whistler-mode waves can be locally generated at a range of high latitudes, within pockets of minimum magnetic field, and propagate equatorward. The obtained dayside waves (in a compressed dipole field) have similar amplitude and frequency content to their dawn sector counterparts (in a pure dipole field) but tend to propagate more field aligned.

Convection in a Very Compressible Fluid: Comparison of Simulations With Experiments

NASA Technical Reports Server (NTRS)

Meyer, H.; Furukawa, A.; Onuki, A.; Kogan, A. B.

2003-01-01

The time profile (Delta)T(t) of the temperature difference, measured across a very compressible fluid layer of supercritical He-3 after the start of a heat flow, shows a damped oscillatory behavior before steady state convection is reached. The results for (Delta)T(t) obtained from numerical simulations and from laboratory experiments are compared over a temperature range where the compressibility varies by a factor of approx. = 40. First the steady-state convective heat current j(sup conv) as a function of the Rayleigh number R(alpha) is presented, and the agreement is found to be good. Second, the shape of the time profile and two characteristic times in the transient part of (Delta)T(t) from simulations and experiments are compared, namely: 1) t(sub osc), the oscillatory period and 2) t(sub p), the time of the first peak after starting the heat flow. These times, scaled by the diffusive time tau(sub D) versus R(alpha), are presented. The agreement is good for t(sup osc)/tau(sub D), where the results collapse on a single curve showing a powerlaw behavior. The simulation hence confirms the universal scaling behavior found experimentally. However for t(sub p)/tau(sub D), where the experimental data also collapse on a single curve, the simulation results show systematic departures from such a behavior. A possible reason for some of the disagreements, both in the time profile and in t(sub p) is discussed. In the Appendix a third characteristic time, t(sub m), between the first peak and the first oscillation minimum is plotted and a comparison between the results of experiments and simulations is made.

Fluid Film Bearing Code Development

NASA Technical Reports Server (NTRS)

1995-01-01

The next generation of rocket engine turbopumps is being developed by industry through Government-directed contracts. These turbopumps will use fluid film bearings because they eliminate the life and shaft-speed limitations of rolling-element bearings, increase turbopump design flexibility, and reduce the need for turbopump overhauls and maintenance. The design of the fluid film bearings for these turbopumps, however, requires sophisticated analysis tools to model the complex physical behavior characteristic of fluid film bearings operating at high speeds with low viscosity fluids. State-of-the-art analysis and design tools are being developed at the Texas A&M University under a grant guided by the NASA Lewis Research Center. The latest version of the code, HYDROFLEXT, is a thermohydrodynamic bulk flow analysis with fluid compressibility, full inertia, and fully developed turbulence models. It can predict the static and dynamic force response of rigid and flexible pad hydrodynamic bearings and of rigid and tilting pad hydrostatic bearings. The Texas A&M code is a comprehensive analysis tool, incorporating key fluid phenomenon pertinent to bearings that operate at high speeds with low-viscosity fluids typical of those used in rocket engine turbopumps. Specifically, the energy equation was implemented into the code to enable fluid properties to vary with temperature and pressure. This is particularly important for cryogenic fluids because their properties are sensitive to temperature as well as pressure. As shown in the figure, predicted bearing mass flow rates vary significantly depending on the fluid model used. Because cryogens are semicompressible fluids and the bearing dynamic characteristics are highly sensitive to fluid compressibility, fluid compressibility effects are also modeled. The code contains fluid properties for liquid hydrogen, liquid oxygen, and liquid nitrogen as well as for water and air. Other fluids can be handled by the code provided that the

Effect of viscosity on the wave propagation: Experimental determination of compression and expansion pulse wave velocity in fluid-fill elastic tube.

PubMed

Stojadinović, Bojana; Tenne, Tamar; Zikich, Dragoslav; Rajković, Nemanja; Milošević, Nebojša; Lazović, Biljana; Žikić, Dejan

2015-11-26

The velocity by which the disturbance travels through the medium is the wave velocity. Pulse wave velocity is one of the main parameters in hemodynamics. The study of wave propagation through the fluid-fill elastic tube is of great importance for the proper biophysical understanding of the nature of blood flow through of cardiovascular system. The effect of viscosity on the pulse wave velocity is generally ignored. In this paper we present the results of experimental measurements of pulse wave velocity (PWV) of compression and expansion waves in elastic tube. The solutions with different density and viscosity were used in the experiment. Biophysical model of the circulatory flow is designed to perform measurements. Experimental results show that the PWV of the expansion waves is higher than the compression waves during the same experimental conditions. It was found that the change in viscosity causes a change of PWV for both waves. We found a relationship between PWV, fluid density and viscosity. Copyright © 2015 Elsevier Ltd. All rights reserved.

Fluid-structure finite-element vibrational analysis

NASA Technical Reports Server (NTRS)

Feng, G. C.; Kiefling, L.

1974-01-01

A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.

Ideal heat transfer conditions for tubular solar receivers with different design constraints

NASA Astrophysics Data System (ADS)

Kim, Jin-Soo; Potter, Daniel; Gardner, Wilson; Too, Yen Chean Soo; Padilla, Ricardo Vasquez

2017-06-01

The optimum heat transfer condition for a tubular type solar receiver was investigated for various receiver pipe size, heat transfer fluid, and design requirement and constraint(s). Heat transfer of a single plain receiver pipe exposed to concentrated solar energy was modelled along the flow path of the heat transfer fluid. Three different working fluids, molten salt, sodium, and supercritical carbon dioxide (sCO2) were considered in the case studies with different design conditions. The optimized ideal heat transfer condition was identified through fast iterative heat transfer calculations solving for all relevant radiation, conduction and convection heat transfers throughout the entire discretized tubular receiver. The ideal condition giving the best performance was obtained by finding the highest acceptable solar energy flux optimally distributed to meet different constraint(s), such as maximum allowable material temperature of receiver, maximum allowable film temperature of heat transfer fluid, and maximum allowable stress of receiver pipe material. The level of fluid side turbulence (represented by pressure drop in this study) was also optimized to give the highest net power production. As the outcome of the study gives information on the most ideal heat transfer condition, it can be used as a useful guideline for optimal design of a real receiver and solar field in a combined manner. The ideal heat transfer condition is especially important for high temperature tubular receivers (e.g. for supplying heat to high efficiency Brayton cycle turbines) where the system design and performance is tightly constrained by the receiver pipe material strength.

Ideal glass transitions in thin films: An energy landscape perspective

NASA Astrophysics Data System (ADS)

Truskett, Thomas M.; Ganesan, Venkat

2003-07-01

We introduce a mean-field model for the potential energy landscape of a thin fluid film confined between parallel substrates. The model predicts how the number of accessible basins on the energy landscape and, consequently, the film's ideal glass transition temperature depend on bulk pressure, film thickness, and the strength of the fluid-fluid and fluid-substrate interactions. The predictions are in qualitative agreement with the experimental trends for the kinetic glass transition temperature of thin films, suggesting the utility of landscape-based approaches for studying the behavior of confined fluids.

Bifurcation and response analysis of a nonlinear flexible rotating disc immersed in bounded compressible fluid

NASA Astrophysics Data System (ADS)

Remigius, W. Dheelibun; Sarkar, Sunetra; Gupta, Sayan

2017-03-01

Use of heavy gases in centrifugal compressors for enhanced oil extraction have made the impellers susceptible to failures through acousto-elastic instabilities. This study focusses on understanding the dynamical behavior of such systems by considering the effects of the bounded fluid housed in a casing on a rotating disc. First, a mathematical model is developed that incorporates the interaction between the rotating impeller - modelled as a flexible disc - and the bounded compressible fluid medium in which it is immersed. The nonlinear effects arising due to large deformations of the disc have been included in the formulation so as to capture the post flutter behavior. A bifurcation analysis is carried out with the disc rotational speed as the bifurcation parameter to investigate the dynamical behavior of the coupled system and estimate the stability boundaries. Parametric studies reveal that the relative strengths of the various dissipation mechanisms in the coupled system play a significant role that affect the bifurcation route and the post flutter behavior in the acousto-elastic system.

Steady fall of isothermal, resistive-viscous, compressible fluid across magnetic field

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

Low, B. C., E-mail: low@ucar.edu; Egan, A. K., E-mail: andrea.egan@colorado.edu

This is a basic MHD study of the steady fall of an infinite, vertical slab of isothermal, resistive-viscous, compressible fluid across a dipped magnetic field in uniform gravity. This double-diffusion steady flow in unbounded space poses a nonlinear but numerically tractable, one-dimensional (1D) free-boundary problem, assuming constant coefficients of resistivity and viscosity. The steady flow is determined by a dimensionless number μ{sub 1} proportional to the triple product of the two diffusion coefficients and the square of the linear total mass. For a sufficiently large μ{sub 1}, the Lorentz, viscous, fluid-pressure, and gravitational forces pack and collimate the fluid intomore » a steady flow of a finite width defined by the two zero-pressure free-boundaries of the slab with vacuum. The viscous force is essential in this collimation effect. The study conjectures that in the regime μ{sub 1}→0, the 1D steady state exists only for μ{sub 1}∈Ω, a spectrum of an infinite number of discrete values, including μ{sub 1} = 0 that corresponds to two steady states, the classical zero-resistivity static slab of Kippenhahn and Schlüter [R. Kippenhahn and A. Schlüter, Z. Astrophys. 43, 36 (1957)] and its recent generalization [B. C. Low et al., Astrophys. J. 755, 34 (2012)] to admit an inviscid resistive flow. The pair of zero-pressure boundaries of each of the μ{sub 1}→0 steady-state slabs are located at infinity. Computational evidence suggests that the Ω steady-states are densely distributed around μ{sub 1} = 0, as an accumulation point, but are sparsely separated by open intervals of μ{sub 1}-values for which the slab must be either time-dependent or spatially multi-dimensional. The widths of these intervals are vanishingly small as μ{sub 1}→0. This topological structure of physical states is similar to that described by Landau and Liftshitz [L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Addison-Wesley, Reading, MA, 1959)] to explain

Creation of an idealized nasopharynx geometry for accurate computational fluid dynamics simulations of nasal airflow in patient-specific models lacking the nasopharynx anatomy

PubMed Central

Borojeni, Azadeh A.T.; Frank-Ito, Dennis O.; Kimbell, Julia S.; Rhee, John S.; Garcia, Guilherme J. M.

2016-01-01

Virtual surgery planning based on computational fluid dynamics (CFD) simulations has the potential to improve surgical outcomes for nasal airway obstruction (NAO) patients, but the benefits of virtual surgery planning must outweigh the risks of radiation exposure. Cone beam computed tomography (CBCT) scans represent an attractive imaging modality for virtual surgery planning due to lower costs and lower radiation exposures compared with conventional CT scans. However, to minimize the radiation exposure, the CBCT sinusitis protocol sometimes images only the nasal cavity, excluding the nasopharynx. The goal of this study was to develop an idealized nasopharynx geometry for accurate representation of outlet boundary conditions when the nasopharynx geometry is unavailable. Anatomically-accurate models of the nasopharynx created from thirty CT scans were intersected with planes rotated at different angles to obtain an average geometry. Cross sections of the idealized nasopharynx were approximated as ellipses with cross-sectional areas and aspect ratios equal to the average in the actual patient-specific models. CFD simulations were performed to investigate whether nasal airflow patterns were affected when the CT-based nasopharynx was replaced by the idealized nasopharynx in 10 NAO patients. Despite the simple form of the idealized geometry, all biophysical variables (nasal resistance, airflow rate, and heat fluxes) were very similar in the idealized vs. patient-specific models. The results confirmed the expectation that the nasopharynx geometry has a minimal effect in the nasal airflow patterns during inspiration. The idealized nasopharynx geometry will be useful in future CFD studies of nasal airflow based on medical images that exclude the nasopharynx. PMID:27525807

Dynamics of Vortex and Magnetic Lines in Ideal Hydrodynamics and MHD

NASA Astrophysics Data System (ADS)

Kuznetsov, E. A.; Ruban, V. P.

Vortex line and magnetic line representations are introduced for description of flows in ideal hydrodynamics and MHD, respectively. For incompressible fluids it is shown that the equations of motion for vorticity φ and magnetic field with the help of this transformation follow from the variational principle. By means of this representation it is possible to integrate the system of hydrodynamic type with the Hamiltonian H=|φ|dr. It is also demonstrated that these representations allow to remove from the noncanonical Poisson brackets, defined on the space of divergence-free vector fields, degeneracy connected with the vorticity frozenness for the Euler equation and with magnetic field frozenness for ideal MHD. For MHD a new Weber type transformation is found. It is shown how this transformation can be obtained from the two-fluid model when electrons and ions can be considered as two independent fluids. The Weber type transformation for ideal MHD gives the whole Lagrangian vector invariant. When this invariant is absent this transformation coincides with the Clebsch representation analog introduced in [1].

On the flow of a compressible fluid by the hodograph method I : unification and extension of present-day results

NASA Technical Reports Server (NTRS)

Garrick, I E; Kaplan, Carl

1944-01-01

Elementary basic solutions of the equations of motion of a compressible fluid in the hodograph variables are developed and used to provide a basis for comparison, in the form of velocity correction formulas, of corresponding compressible and incompressible flows. The known approximate results of Chaplygin, Von Karman and Tsien, Temple and Yarwood, and Prandtl and Glauert are unified by means of the analysis of the present paper. Two new types of approximations, obtained from the basic solutions, are introduced; they possess certain desirable features of the other approximations and appear preferable as a basis for extrapolation into the range of high stream Mach numbers and large disturbances to the main stream. Tables and figures giving velocity and pressure-coefficient correction factors are included in order to facilitate the practical application of the results.

Effects of real fluid properties on axial turbine meanline design and off-design analysis

NASA Astrophysics Data System (ADS)

MacLean, Cameron

The effects of real fluid properties on axial turbine meanline analysis have been investigated employing two meanline analysis codes, namely Turbine Meanline Design (TMLD) and Turbine Meanline Off-Design (TMLO). The previously developed TMLD code assumed the working fluid was an ideal gas. Therefore it was modified to use real fluid properties. TMLO was then developed from TMLD Both codes can be run using either the ideal gas assumption or real fluid properties. TMLD was employed for the meanline design of several axial turbines for a range of inlet conditions, using both the ideal gas assumption and real fluid properties. The resulting designs were compared to see the effects of real fluid properties. Meanline designs, generated using the ideal gas assumption, were then analysed with TMLO using real fluid properties. This was done over a range of inlet conditions that correspond to varying degrees of departure from ideal gas conditions. The goal was to show how machines designed with the ideal gas assumption would perform with the real working fluid. The working fluid used in both investigations was supercritical carbon dioxide. Results from the investigation show that real fluid properties had a strong effect on the gas path areas of the turbine designs as well as the performance of turbines designed using the ideal gas assumption. Specifically, power output and the velocities of the working fluid were affected. It was found that accounting for losses tended to lessen the effects of the real fluid properties.

Current Results and Proposed Activities in Microgravity Fluid Dynamics

NASA Technical Reports Server (NTRS)

Polezhaev, V. I.

1996-01-01

The Institute for Problems in Mechanics' Laboratory work in mathematical and physical modelling of fluid mechanics develops models, methods, and software for analysis of fluid flow, instability analysis, direct numerical modelling and semi-empirical models of turbulence, as well as experimental research and verification of these models and their applications in technological fluid dynamics, microgravity fluid mechanics, geophysics, and a number of engineering problems. This paper presents an overview of the results in microgravity fluid dynamics research during the last two years. Nonlinear problems of weakly compressible and compressible fluid flows are discussed.

On the Representation of Aquifer Compressibility in General Subsurface Flow Codes: How an Alternate Definition of Aquifer Compressibility Matches Results from the Groundwater Flow Equation

NASA Astrophysics Data System (ADS)

Birdsell, D.; Karra, S.; Rajaram, H.

2017-12-01

The governing equations for subsurface flow codes in deformable porous media are derived from the fluid mass balance equation. One class of these codes, which we call general subsurface flow (GSF) codes, does not explicitly track the motion of the solid porous media but does accept general constitutive relations for porosity, density, and fluid flux. Examples of GSF codes include PFLOTRAN, FEHM, STOMP, and TOUGH2. Meanwhile, analytical and numerical solutions based on the groundwater flow equation have assumed forms for porosity, density, and fluid flux. We review the derivation of the groundwater flow equation, which uses the form of Darcy's equation that accounts for the velocity of fluids with respect to solids and defines the soil matrix compressibility accordingly. We then show how GSF codes have a different governing equation if they use the form of Darcy's equation that is written only in terms of fluid velocity. The difference is seen in the porosity change, which is part of the specific storage term in the groundwater flow equation. We propose an alternative definition of soil matrix compressibility to correct for the untracked solid velocity. Simulation results show significantly less error for our new compressibility definition than the traditional compressibility when compared to analytical solutions from the groundwater literature. For example, the error in one calculation for a pumped sandstone aquifer goes from 940 to <70 Pa when the new compressibility is used. Code users and developers need to be aware of assumptions in the governing equations and constitutive relations in subsurface flow codes, and our newly-proposed compressibility function should be incorporated into GSF codes.

On the Representation of Aquifer Compressibility in General Subsurface Flow Codes: How an Alternate Definition of Aquifer Compressibility Matches Results from the Groundwater Flow Equation

NASA Astrophysics Data System (ADS)

Birdsell, D.; Karra, S.; Rajaram, H.

2016-12-01

The governing equations for subsurface flow codes in deformable porous media are derived from the fluid mass balance equation. One class of these codes, which we call general subsurface flow (GSF) codes, does not explicitly track the motion of the solid porous media but does accept general constitutive relations for porosity, density, and fluid flux. Examples of GSF codes include PFLOTRAN, FEHM, STOMP, and TOUGH2. Meanwhile, analytical and numerical solutions based on the groundwater flow equation have assumed forms for porosity, density, and fluid flux. We review the derivation of the groundwater flow equation, which uses the form of Darcy's equation that accounts for the velocity of fluids with respect to solids and defines the soil matrix compressibility accordingly. We then show how GSF codes have a different governing equation if they use the form of Darcy's equation that is written only in terms of fluid velocity. The difference is seen in the porosity change, which is part of the specific storage term in the groundwater flow equation. We propose an alternative definition of soil matrix compressibility to correct for the untracked solid velocity. Simulation results show significantly less error for our new compressibility definition than the traditional compressibility when compared to analytical solutions from the groundwater literature. For example, the error in one calculation for a pumped sandstone aquifer goes from 940 to <70 Pa when the new compressibility is used. Code users and developers need to be aware of assumptions in the governing equations and constitutive relations in subsurface flow codes, and our newly-proposed compressibility function should be incorporated into GSF codes.

When an Adiabatic Irreversible Expansion or Compression Becomes Reversible

ERIC Educational Resources Information Center

Anacleto, Joaquim; Ferreira, J. M.; Soares, A. A.

2009-01-01

This paper aims to contribute to a better understanding of the concepts of a "reversible process" and "entropy". For this purpose, an adiabatic irreversible expansion or compression is analysed, by considering that an ideal gas is expanded (compressed), from an initial pressure P[subscript i] to a final pressure P[subscript f], by being placed in…

Sphere based fluid systems

NASA Technical Reports Server (NTRS)

Elleman, Daniel D. (Inventor); Wang, Taylor G. (Inventor)

1989-01-01

Systems are described for using multiple closely-packed spheres. In one system for passing fluid, a multiplicity of spheres lie within a container, with all of the spheres having the same outside diameter and with the spheres being closely nested in one another to create multiple interstitial passages of a known size and configuration and smooth walls. The container has an inlet and outlet for passing fluid through the interstitial passages formed between the nested spheres. The small interstitial passages can be used to filter out material, especially biological material such as cells in a fluid, where the cells can be easily destroyed if passed across sharp edges. The outer surface of the spheres can contain a material that absorbs a constitutent in the flowing fluid, such as a particular contamination gas, or can contain a catalyst to chemically react the fluid passing therethrough, the use of multiple small spheres assuring a large area of contact of these surfaces of the spheres with the fluid. In a system for storing and releasing a fluid such as hydrogen as a fuel, the spheres can include a hollow shell containing the fluid to be stored, and located within a compressable container that can be compressed to break the shells and release the stored fluid.

Simple cubic equation of state applied to hard-sphere, Lennard-Jones fluids, simple fluids and solids

NASA Astrophysics Data System (ADS)

Sun, Jiu-Xun; Cai, Ling-Cang; Wu, Qiang; Jin, Ke

2013-09-01

Based on the expansion and extension of the virial equation of state (EOS) of hard-sphere fluids solved by the Percus-Yevick integration equation, a universal cubic (UC) EOS is developed. The UC EOS is applied to model hard-sphere and Lennard-Jones (LJ) fluids, simple Ar and N2 liquids at low temperatures, and supercritical Ar and N2 fluids at high temperatures, as well as ten solids, respectively. The three parameters are determined for the hard-sphere fluid by fitting molecular dynamics (MD) simulation data of the third to eighth virial coefficients in the literature; for other fluids by fitting isothermal compression data; and for solids by using the Einstein model. The results show that the UC EOS gives better results than the Carnahan-Starling EOS for compressibility of hard-sphere fluids. The Helmholtz free energy and internal energy for LJ fluids are predicted and compared with MD simulation data. The calculated pressures for simple Ar and N2 liquids are compared with experimental data. The agreement is fairly good. Eight three-parameter EOSs are applied to describe isothermals of ten typical solids. It is shown that the UC EOS gives the best precision with correct behavior at high-pressure limitation. The UC EOS considering thermal effects is used to analytically evaluate the isobaric thermal expansivity and isothermal compressibility coefficients. The results are in good agreement with experimental data.

The deformation of gum metal under nanoindentation and sub-micron pillar compression

NASA Astrophysics Data System (ADS)

Withey, Elizabeth Ann

Reaching ideal strength has proven to be difficult in most materials. Dislocation slip, phase transformations, twinning, and fracture all tend to occur at stresses well below the ideal strength of a material. Only on very small scales has it been possible to approach ideal strength. Thus, it was of great interest when a set of beta-Ti alloys, Gum Metal, were found to have a bulk yield strength close to half of its ideal strength. However, some recent studies have questioned the reliability of this claim. Several studies have suggested Gum Metal deforms by dislocation slip. Others have suggested the possibility of transformation-induced plasticity. The present study was undertaken in order to help clarify if and how Gum Metal can reach ideal strength. Two different experiments, ex situ nanoindentation and quantitative in situ nanopillar compression in a transmission electron microscope to correlate real-time deformation behavior, were performed on a single composition of Gum Metal, Ti-23Nb-0.7Ta-2Zr-1.20 at. %, obtained from Toyota Central R&D Laboratories. Nanoindented specimens were thinned from the bottom surface until the pits of multiple indentations became electron-transparent allowing for qualitative analysis of the deformation microstructure in both fully cold-worked and solution-treated specimens. Real-time load-displacement behavior from the nanopillar compression tests was correlated with real-time video recorded during each compression to determine both the compressive strength of each pillar and the timing and strengths of different deformation behaviors observed. Combining the results from both experiments provided several important conclusions. First, Gum Metal approaches and can attain ideal strength in nanopillars regardless of processing condition. While dislocations exist in Gum Metal, they can be tightly pinned by obstacles with spacing less than ˜20 nm, which should inhibit their motion at strengths below the ideal shear strength. The plastic

The Effect of Viscosity of a Fluid on the Frequency Response of a Viscoelastic Plate Loaded by This Fluid

NASA Astrophysics Data System (ADS)

Zamanov, A. D.; Ismailov, M. I.; Akbarov, S. D.

2018-03-01

A hydroviscoelastic system consisting of a viscoelastic plate and a half-plane filled with a viscous fluid is considered. The effect of viscosity of the fluid on the frequency response of the system and its dependence on the rheological parameters of plate material are estimated. The problem on forced vibrations of the system in the plane strain state is investigated using the exact equations of viscoelastodynamics for describing the motion of the plate and linearized Navier-Stokes equations for describing the flow of the fluid. The results found in the cases of nonviscous compressible and Newtonian compressible viscous fluids are compared.

Advanced Fluid Reduced Order Models for Compressible Flow.

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

Tezaur, Irina Kalashnikova; Fike, Jeffrey A.; Carlberg, Kevin Thomas

This report summarizes fiscal year (FY) 2017 progress towards developing and implementing within the SPARC in-house finite volume flow solver advanced fluid reduced order models (ROMs) for compressible captive-carriage flow problems of interest to Sandia National Laboratories for the design and qualification of nuclear weapons components. The proposed projection-based model order reduction (MOR) approach, known as the Proper Orthogonal Decomposition (POD)/Least- Squares Petrov-Galerkin (LSPG) method, can substantially reduce the CPU-time requirement for these simulations, thereby enabling advanced analyses such as uncertainty quantification and de- sign optimization. Following a description of the project objectives and FY17 targets, we overview briefly themore » POD/LSPG approach to model reduction implemented within SPARC . We then study the viability of these ROMs for long-time predictive simulations in the context of a two-dimensional viscous laminar cavity problem, and describe some FY17 enhancements to the proposed model reduction methodology that led to ROMs with improved predictive capabilities. Also described in this report are some FY17 efforts pursued in parallel to the primary objective of determining whether the ROMs in SPARC are viable for the targeted application. These include the implemen- tation and verification of some higher-order finite volume discretization methods within SPARC (towards using the code to study the viability of ROMs on three-dimensional cavity problems) and a novel structure-preserving constrained POD/LSPG formulation that can improve the accuracy of projection-based reduced order models. We conclude the report by summarizing the key takeaways from our FY17 findings, and providing some perspectives for future work.« less

A parametric analysis of waves propagating in a porous solid saturated by a three-phase fluid.

PubMed

Santos, Juan E; Savioli, Gabriela B

2015-11-01

This paper presents an analysis of a model for the propagation of waves in a poroelastic solid saturated by a three-phase viscous, compressible fluid. The constitutive relations and the equations of motion are stated first. Then a plane wave analysis determines the phase velocities and attenuation coefficients of the four compressional waves and one shear wave that propagate in this type of medium. A procedure to compute the elastic constants in the constitutive relations is defined next. Assuming the knowledge of the shear modulus of the dry matrix, the other elastic constants in the stress-strain relations are determined by employing ideal gedanken experiments generalizing those of Biot's theory for single-phase fluids. These experiments yield expressions for the elastic constants in terms of the properties of the individual solid and fluids phases. Finally the phase velocities and attenuation coefficients of all waves are computed for a sample of Berea sandstone saturated by oil, gas, and water.

Dynamic compression of dense oxide (Gd 3Ga 5O 12) from 0.4 to 2.6 TPa: Universal Hugoniot of fluid metals

DOE PAGES

Ozaki, N.; Nellis, W. J.; Mashimo, T.; ...

2016-05-19

Materials at high pressures and temperatures are of great current interest for warm dense matter physics, planetary sciences, and inertial fusion energy research. Shock-compression equation-of-state data and optical reflectivities of the fluid dense oxide, Gd 3Ga 5O 12 (GGG), were measured at extremely high pressures up to 2.6 TPa (26 Mbar) generated by high-power laser irradiation and magnetically-driven hypervelocity impacts. Above 0.75 TPa, the GGG Hugoniot data approach/reach a universal linear line of fluid metals, and the optical reflectivity most likely reaches a constant value indicating that GGG undergoes a crossover from fluid semiconductor to poor metal with minimum metallicmore » conductivity (MMC). These results suggest that most fluid compounds, e.g., strong planetary oxides, reach a common state on the universal Hugoniot of fluid metals (UHFM) with MMC at sufficiently extreme pressures and temperatures. Lastly, the systematic behaviors of warm dense fluid would be useful benchmarks for developing theoretical equation-of-state and transport models in the warm dense matter regime in determining computational predictions.« less

Review on Applications of NanoFluids used in Vapour Compression Refrigeration System for Cop Enhancement

NASA Astrophysics Data System (ADS)

Veera Raghavalu, K.; Govindha Rasu, N.

2018-03-01

The present research paper focuses on the use of Nano additive refrigerants in vapor compression refrigeration system (VCRS) because of their amazing development during Thermo Physical along with heat transfer potential to improve the coefficient of performance (COP) and reliability of refrigeration system. Furthermore, challenges and future instructions of performance enhancement of VCRS using Nano additive refrigerants were presented. Lubricant oil is essential in the entire vapour compression refrigeration systems, mostly for the efficient function of the compressor. But, some assign of the oil is entire the cycle oil circulates with the refrigerant. Presently, an assortment of investigation is going on in the field of the Nano-particles like metals, oxides, carbon Nano-tubes or carbides. Nano-lubricants are unique type of Nano-fluids which are varieties of Nano-particles, lubricants and have a wide variety in the fields of refrigeration systems. This paper, has been done on the application of Nano-particles balanced in lubricating oils of refrigerating systems are reviewed. The aim of this investigation is to study and find which type of lubricant oil works better with Nano-particles in the area of refrigeration. From the review of literature, it has been observed that Nano-particles mixed with mineral oil gives enhanced results than polyolester (POE) oil.

Fluid driven torsional dipole seismic source

DOEpatents

Hardee, Harry C.

1991-01-01

A compressible fluid powered oscillating downhole seismic source device capable of periodically generating uncontaminated horizontally-propagated, shear waves is provided. A compressible fluid generated oscillation is created within the device which imparts an oscillation to a housing when the device is installed in a housing such as the cylinder off an existing downhole tool, thereby a torsional seismic source is established. Horizontal waves are transferred to the surrounding bore hole medium through downhole clamping.

High Order Approximations for Compressible Fluid Dynamics on Unstructured and Cartesian Meshes

NASA Technical Reports Server (NTRS)

Barth, Timothy (Editor); Deconinck, Herman (Editor)

1999-01-01

The development of high-order accurate numerical discretization techniques for irregular domains and meshes is often cited as one of the remaining challenges facing the field of computational fluid dynamics. In structural mechanics, the advantages of high-order finite element approximation are widely recognized. This is especially true when high-order element approximation is combined with element refinement (h-p refinement). In computational fluid dynamics, high-order discretization methods are infrequently used in the computation of compressible fluid flow. The hyperbolic nature of the governing equations and the presence of solution discontinuities makes high-order accuracy difficult to achieve. Consequently, second-order accurate methods are still predominately used in industrial applications even though evidence suggests that high-order methods may offer a way to significantly improve the resolution and accuracy for these calculations. To address this important topic, a special course was jointly organized by the Applied Vehicle Technology Panel of NATO's Research and Technology Organization (RTO), the von Karman Institute for Fluid Dynamics, and the Numerical Aerospace Simulation Division at the NASA Ames Research Center. The NATO RTO sponsored course entitled "Higher Order Discretization Methods in Computational Fluid Dynamics" was held September 14-18, 1998 at the von Karman Institute for Fluid Dynamics in Belgium and September 21-25, 1998 at the NASA Ames Research Center in the United States. During this special course, lecturers from Europe and the United States gave a series of comprehensive lectures on advanced topics related to the high-order numerical discretization of partial differential equations with primary emphasis given to computational fluid dynamics (CFD). Additional consideration was given to topics in computational physics such as the high-order discretization of the Hamilton-Jacobi, Helmholtz, and elasticity equations. This volume consists

Comparison of cardiac magnetic resonance imaging and bio-impedance spectroscopy for the assessment of fluid displacement induced by external leg compression.

PubMed

Saporito, Salvatore; Dovancescu, Silviu; Herold, Ingeborg H F; van den Bosch, Harrie C M; van Assen, Hans C; Aarts, Ronald M; Korsten, Hendrikus H M; Mischi, Massimo

2017-01-01

Heart failure is marked by frequent hospital admissions, often as a consequence of pulmonary congestion. Current gold standard techniques for thoracic fluid measurement require invasive heamodynamic access and therefore they are not suitable for continuous monitoring. Changes in thoracic impedance (TI) may enable non-invasive early detection of congestion and prevention of unplanned hospitalizations. However, the usefulness of TI to assess thoracic fluid status is limited by inter-subject variability and by the lack of reliable normalization methods. Indicator dilution methods allow absolute fluid volume estimation; cardiac magnetic resonance (CMR) has been recently proposed to apply indicator dilution methods in a minimally-invasive manner. In this study, we aim to compare bio-impedance spectroscopy (BIS) and CMR for the assessment of thoracic fluid status, and to determine their ability to detect fluid displacement induced by a leg compression procedure in healthy volunteers. A pressure gradient was applied across each subject's legs for 5 min (100-60 mmHg, distal to proximal). Each subject underwent a continuous TI-BIS measurement during the procedure, and repeated CMR-based indicator dilution measurements on a 1.5 T scanner at baseline, during compression, and after pressure release. The Cole-Cole and the local density random walk models were used for parameter extraction from TI-BIS and indicator dilution measurements, respectively. Intra-thoracic blood volume index (ITBI) derived from CMR, and extracellular fluid resistance (R E ) from TI-BIS, were considered as thoracic fluid status measures. Eight healthy volunteers were included in this study. An increase in ITBI of 45.2  ±  47.2 ml m -2 was observed after the leg inflation (13.1  ±  15.1% w.r.t. baseline, p  <  0.05), while a decrease of  -0.84  ±  0.39 Ω in R E (-1.7  ±  0.9% w.r.t. baseline, p  <  0.05) was observed. ITBV and R E normalized by

Effect of rock rheology on fluid leak- off during hydraulic fracturing

NASA Astrophysics Data System (ADS)

Yarushina, V. M.; Bercovici, D.; Oristaglio, M. L.

2012-04-01

In this communication, we evaluate the effect of rock rheology on fluid leak­off during hydraulic fracturing of reservoirs. Fluid leak-off in hydraulic fracturing is often nonlinear. The simple linear model developed by Carter (1957) for flow of fracturing fluid into a reservoir has three different regions in the fractured zone: a filter cake on the fracture face, formed by solid additives from the fracturing fluid; a filtrate zone affected by invasion of the fracturing fluid; and a reservoir zone with the original formation fluid. The width of each zone, as well as its permeability and pressure drop, is assumed to remain constant. Physical intuition suggests some straightforward corrections to this classical theory to take into account the pressure dependence of permeability, the compressibility or non-Newtonian rheology of fracturing fluid, and the radial (versus linear) geometry of fluid leak­off from the borehole. All of these refinements, however, still assume that the reservoir rock adjacent to the fracture face is non­deformable. Although the effect of poroelastic stress changes on leak-off is usually thought to be negligible, at the very high fluid pressures used in hydraulic fracturing, where the stresses exceed the rock strength, elastic rheology may not be the best choice. For example, calculations show that perfectly elastic rock formations do not undergo the degree of compaction typically seen in sedimentary basins. Therefore, pseudo-elastic or elastoplastic models are used to fit observed porosity profiles with depth. Starting from balance equations for mass and momentum for fluid and rock, we derive a hydraulic flow equation coupled with a porosity equation describing rock compaction. The result resembles a pressure diffusion equation with the total compressibility being a sum of fluid, rock and pore-space compressibilities. With linear elastic rheology, the bulk formation compressibility is dominated by fluid compressibility. But the possibility

The Statistical Mechanics of Ideal MHD Turbulence

NASA Technical Reports Server (NTRS)

Shebalin, John V.

2003-01-01

Turbulence is a universal, nonlinear phenomenon found in all energetic fluid and plasma motion. In particular. understanding magneto hydrodynamic (MHD) turbulence and incorporating its effects in the computation and prediction of the flow of ionized gases in space, for example, are great challenges that must be met if such computations and predictions are to be meaningful. Although a general solution to the "problem of turbulence" does not exist in closed form, numerical integrations allow us to explore the phase space of solutions for both ideal and dissipative flows. For homogeneous, incompressible turbulence, Fourier methods are appropriate, and phase space is defined by the Fourier coefficients of the physical fields. In the case of ideal MHD flows, a fairly robust statistical mechanics has been developed, in which the symmetry and ergodic properties of phase space is understood. A discussion of these properties will illuminate our principal discovery: Coherent structure and randomness co-exist in ideal MHD turbulence. For dissipative flows, as opposed to ideal flows, progress beyond the dimensional analysis of Kolmogorov has been difficult. Here, some possible future directions that draw on the ideal results will also be discussed. Our conclusion will be that while ideal turbulence is now well understood, real turbulence still presents great challenges.

Effect of pairwise additivity on finite-temperature behavior of classical ideal gas

NASA Astrophysics Data System (ADS)

Shekaari, Ashkan; Jafari, Mahmoud

2018-05-01

Finite-temperature molecular dynamics simulations have been applied to inquire into the effect of pairwise additivity on the behavior of classical ideal gas within the temperature range of T = 250-4000 K via applying a variety of pair potentials and then examining the temperature dependence of a number of thermodynamical properties. Examining the compressibility factor reveals the most deviation from ideal-gas behavior for the Lennard-Jones system mainly due to the presence of both the attractive and repulsive terms. The systems with either attractive or repulsive intermolecular potentials are found to present no resemblance to real gases, but the most similarity to the ideal one as temperature rises.

Vapor Compression Cycle Design Program (CYCLE_D)

National Institute of Standards and Technology Data Gateway

SRD 49 NIST Vapor Compression Cycle Design Program (CYCLE_D) (PC database for purchase)   The CYCLE_D database package simulates the vapor compression refrigeration cycles. It is fully compatible with REFPROP 9.0 and covers the 62 single-compound refrigerants . Fluids can be used in mixtures comprising up to five components.

Condensation and critical exponents of an ideal non-Abelian gas

NASA Astrophysics Data System (ADS)

Talaei, Zahra; Mirza, Behrouz; Mohammadzadeh, Hosein

2017-11-01

We investigate an ideal gas obeying non-Abelian statistics and derive the expressions for some thermodynamic quantities. It is found that thermodynamic quantities are finite at the condensation point where their derivatives diverge and, near this point, they behave as \\vert T-Tc\\vert^{-ρ} in which Tc denotes the condensation temperature and ρ is a critical exponent. The critical exponents related to the heat capacity and compressibility are obtained by fitting numerical results and others are obtained using the scaling law hypothesis for a three-dimensional non-Abelian ideal gas. This set of critical exponents introduces a new universality class.

The global existence and large time behavior of smooth compressible fluid in an infinitely expanding ball, III: The 3-D Boltzmann equation

NASA Astrophysics Data System (ADS)

Yin, Huicheng; Zhao, Wenbin

2018-01-01

This paper is a continuation of the works in [35] and [37], where the authors have established the global existence of smooth compressible flows in infinitely expanding balls for inviscid gases and viscid gases, respectively. In this paper, we are concerned with the global existence and large time behavior of compressible Boltzmann gases in an infinitely expanding ball. Such a problem is one of the interesting models in studying the theory of global smooth solutions to multidimensional compressible gases with time dependent boundaries and vacuum states at infinite time. Due to the conservation of mass, the fluid in the expanding ball becomes rarefied and eventually tends to a vacuum state meanwhile there are no appearances of vacuum domains in any part of the expansive ball, which is easily observed in finite time. In the present paper, we will confirm this physical phenomenon for the Boltzmann equation by obtaining the exact lower and upper bound on the macroscopic density function.

Pancreatic fluid collections: What is the ideal imaging technique?

PubMed

Dhaka, Narendra; Samanta, Jayanta; Kochhar, Suman; Kalra, Navin; Appasani, Sreekanth; Manrai, Manish; Kochhar, Rakesh

2015-12-28

Pancreatic fluid collections (PFCs) are seen in up to 50% of cases of acute pancreatitis. The Revised Atlanta classification categorized these collections on the basis of duration of disease and contents, whether liquid alone or a mixture of fluid and necrotic debris. Management of these different types of collections differs because of the variable quantity of debris; while patients with pseudocysts can be drained by straight-forward stent placement, walled-off necrosis requires multi-disciplinary approach. Differentiating these collections on the basis of clinical severity alone is not reliable, so imaging is primarily performed. Contrast-enhanced computed tomography is the commonly used modality for the diagnosis and assessment of proportion of solid contents in PFCs; however with certain limitations such as use of iodinated contrast material especially in renal failure patients and radiation exposure. Magnetic resonance imaging (MRI) performs better than computed tomography (CT) in characterization of pancreatic/peripancreatic fluid collections especially for quantification of solid debris and fat necrosis (seen as fat density globules), and is an alternative in those situations where CT is contraindicated. Also magnetic resonance cholangiopancreatography is highly sensitive for detecting pancreatic duct disruption and choledocholithiasis. Endoscopic ultrasound is an evolving technique with higher reproducibility for fluid-to-debris component estimation with the added advantage of being a single stage procedure for both diagnosis (solid debris delineation) and management (drainage of collection) in the same sitting. Recently role of diffusion weighted MRI and positron emission tomography/CT with (18)F-FDG labeled autologous leukocytes is also emerging for detection of infection noninvasively. Comparative studies between these imaging modalities are still limited. However we look forward to a time when this gap in literature will be fulfilled.

Generalized Fluid System Simulation Program, Version 6.0

NASA Technical Reports Server (NTRS)

Majumdar, A. K.; LeClair, A. C.; Moore, R.; Schallhorn, P. A.

2016-01-01

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct

Fluid channeling system

NASA Technical Reports Server (NTRS)

Davis, Donald Y. (Inventor); Hitch, Bradley D. (Inventor)

1994-01-01

A fluid channeling system includes a fluid ejector, a heat exchanger, and a fluid pump disposed in series flow communication The ejector includes a primary inlet for receiving a primary fluid, and a secondary inlet for receiving a secondary fluid which is mixed with the primary fluid and discharged therefrom as ejector discharge. Heat is removed from the ejector discharge in the heat exchanger, and the heat exchanger discharge is compressed in the fluid pump and channeled to the ejector secondary inlet as the secondary fluid In an exemplary embodiment, the temperature of the primary fluid is greater than the maximum operating temperature of a fluid motor powering the fluid pump using a portion of the ejector discharge, with the secondary fluid being mixed with the primary fluid so that the ejector discharge temperature is equal to about the maximum operating temperature of the fluid motor.

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.

Dynamic compression of dense oxide (Gd3Ga5O12) from 0.4 to 2.6 TPa: Universal Hugoniot of fluid metals

PubMed Central

Ozaki, N.; Nellis, W. J.; Mashimo, T.; Ramzan, M.; Ahuja, R.; Kaewmaraya, T.; Kimura, T.; Knudson, M.; Miyanishi, K.; Sakawa, Y.; Sano, T.; Kodama, R.

2016-01-01

Materials at high pressures and temperatures are of great current interest for warm dense matter physics, planetary sciences, and inertial fusion energy research. Shock-compression equation-of-state data and optical reflectivities of the fluid dense oxide, Gd3Ga5O12 (GGG), were measured at extremely high pressures up to 2.6 TPa (26 Mbar) generated by high-power laser irradiation and magnetically-driven hypervelocity impacts. Above 0.75 TPa, the GGG Hugoniot data approach/reach a universal linear line of fluid metals, and the optical reflectivity most likely reaches a constant value indicating that GGG undergoes a crossover from fluid semiconductor to poor metal with minimum metallic conductivity (MMC). These results suggest that most fluid compounds, e.g., strong planetary oxides, reach a common state on the universal Hugoniot of fluid metals (UHFM) with MMC at sufficiently extreme pressures and temperatures. The systematic behaviors of warm dense fluid would be useful benchmarks for developing theoretical equation-of-state and transport models in the warm dense matter regime in determining computational predictions. PMID:27193942

Dynamic compression of dense oxide (Gd3Ga5O12) from 0.4 to 2.6 TPa: Universal Hugoniot of fluid metals.

PubMed

Ozaki, N; Nellis, W J; Mashimo, T; Ramzan, M; Ahuja, R; Kaewmaraya, T; Kimura, T; Knudson, M; Miyanishi, K; Sakawa, Y; Sano, T; Kodama, R

2016-05-19

Materials at high pressures and temperatures are of great current interest for warm dense matter physics, planetary sciences, and inertial fusion energy research. Shock-compression equation-of-state data and optical reflectivities of the fluid dense oxide, Gd3Ga5O12 (GGG), were measured at extremely high pressures up to 2.6 TPa (26 Mbar) generated by high-power laser irradiation and magnetically-driven hypervelocity impacts. Above 0.75 TPa, the GGG Hugoniot data approach/reach a universal linear line of fluid metals, and the optical reflectivity most likely reaches a constant value indicating that GGG undergoes a crossover from fluid semiconductor to poor metal with minimum metallic conductivity (MMC). These results suggest that most fluid compounds, e.g., strong planetary oxides, reach a common state on the universal Hugoniot of fluid metals (UHFM) with MMC at sufficiently extreme pressures and temperatures. The systematic behaviors of warm dense fluid would be useful benchmarks for developing theoretical equation-of-state and transport models in the warm dense matter regime in determining computational predictions.

Rock deformation models and fluid leak-off in hydraulic fracturing

NASA Astrophysics Data System (ADS)

Yarushina, Viktoriya M.; Bercovici, David; Oristaglio, Michael L.

2013-09-01

Fluid loss into reservoir rocks during hydraulic fracturing is modelled via a poro-elastoplastic pressure diffusion equation in which the total compressibility is a sum of fluid, rock and pore space compressibilities. Inclusion of pore compressibility and porosity-dependent permeability in the model leads to a strong pressure dependence of leak-off (i.e. drainage rate). Dilation of the matrix due to fluid invasion causes higher rates of fluid leak-off. The present model is appropriate for naturally fractured and tight gas reservoirs as well as for soft and poorly consolidated formations whose mechanical behaviour departs from simple elastic laws. Enhancement of the leak-off coefficient by dilation, predicted by the new model, may help explain the low percentage recovery of fracturing fluid (usually between 5 and 50 per cent) in shale gas stimulation by hydraulic fracturing.

Variational principles for stochastic fluid dynamics

PubMed Central

Holm, Darryl D.

2015-01-01

This paper derives stochastic partial differential equations (SPDEs) for fluid dynamics from a stochastic variational principle (SVP). The paper proceeds by taking variations in the SVP to derive stochastic Stratonovich fluid equations; writing their Itô representation; and then investigating the properties of these stochastic fluid models in comparison with each other, and with the corresponding deterministic fluid models. The circulation properties of the stochastic Stratonovich fluid equations are found to closely mimic those of the deterministic ideal fluid models. As with deterministic ideal flows, motion along the stochastic Stratonovich paths also preserves the helicity of the vortex field lines in incompressible stochastic flows. However, these Stratonovich properties are not apparent in the equivalent Itô representation, because they are disguised by the quadratic covariation drift term arising in the Stratonovich to Itô transformation. This term is a geometric generalization of the quadratic covariation drift term already found for scalar densities in Stratonovich's famous 1966 paper. The paper also derives motion equations for two examples of stochastic geophysical fluid dynamics; namely, the Euler–Boussinesq and quasi-geostropic approximations. PMID:27547083

Thermodynamic Behavior in Decaying Compressible Turbulence with Initially Dominant Temperature Fluctuations

NASA Astrophysics Data System (ADS)

Cai, X. D.; O'Brien, Edward E.; Ladeinde, Foluso

1996-11-01

Direct numerical simulation of decaying, isotropic, compressible turbulence in three dimensions is used to examine the behavior of fluctuations in density, temperature, and pressure when the initial conditions include temperature fluctuations larger than pressure fluctuations. The numerical procedure is described elsewhere (Ladeinde, F. et al.,) Phys. Fluids 7(11), pp. 2848 (1995), the initial turbulence Mach number range is subsonic, 0.3 to 0.7, and, following Ghosh and Matthaeus(Ghosh, S. and Matthaeus, W. H. Phys. Fluids A, pp. 148 (1991)), the initial compressible turbulence is characterized as a: mostly solenoidal, b: random, or c: longitudinal. These cases represent, respectively, ratios of initial kinetic energy in the compressible modes to total initial kinetic energy, say \\chi_0, which are either a: very small, b: about 0.6, or c: near unity. Thermodynamic scalings at the lowest values of initial Mach number and \\chi0 follow the predictions of Zank and Matthaeus (Zank, G. P. and Matthaeus, W. H. Phys. Fluids A(3), pp. 69 (1991)), but not otherwise. The relationship between \\chi, Mach number, and compressible pressure predicted by Sarkar et al.(Sarkar, S. et al.,) J. Fluid Mech. 227, pp. 473 (1991) applies, on average, to all cases computed.

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.

Conformal mapping for the Helmholtz equation: acoustic wave scattering by a two dimensional inclusion with irregular shape in an ideal fluid.

PubMed

Liu, Gang; Jayathilake, Pahala G; Khoo, Boo Cheong; Han, Feng; Liu, Dian Kui

2012-02-01

The complex variables method with mapping function was extended to solve the linear acoustic wave scattering by an inclusion with sharp/smooth corners in an infinite ideal fluid domain. The improved solutions of Helmholtz equation, shown as Bessel function with mapping function as the argument and fractional order Bessel function, were analytically obtained. Based on the mapping function, the initial geometry as well as the original physical vector can be transformed into the corresponding expressions inside the mapping plane. As all the physical vectors are calculated in the mapping plane (η,η), this method can lead to potential vast savings of computational resources and memory. In this work, the results are validated against several published works in the literature. The different geometries of the inclusion with sharp corners based on the proposed mapping functions for irregular polygons are studied and discussed. The findings show that the variation of angles and frequencies of the incident waves have significant influence on the bistatic scattering pattern and the far-field form factor for the pressure in the fluid. © 2012 Acoustical Society of America

MHD simulation of plasma compression experiments

NASA Astrophysics Data System (ADS)

Reynolds, Meritt; Barsky, Sandra; de Vietien, Peter

2017-10-01

General Fusion (GF) is working to build a magnetized target fusion (MTF) power plant based on compression of magnetically-confined plasma by liquid metal. GF is testing this compression concept by collapsing solid aluminum liners onto plasmas formed by coaxial helicity injection in a series of experiments called PCS (Plasma Compression, Small). We simulate the PCS experiments using the finite-volume MHD code VAC. The single-fluid plasma model includes temperature-dependent resistivity and anisotropic heat transport. The time-dependent curvilinear mesh for MHD simulation is derived from LS-DYNA simulations of actual field tests of liner implosion. We will discuss how 3D simulations reproduced instability observed in the PCS13 experiment and correctly predicted stabilization of PCS14 by ramping the shaft current during compression. We will also present a comparison of simulated Mirnov and x-ray diagnostics with experimental measurements indicating that PCS14 compressed well to a linear compression ratio of 2.5:1.

Reduced viscosity interpreted for fluid/gas mixtures

NASA Technical Reports Server (NTRS)

Lewis, D. H.

1981-01-01

Analysis predicts decrease in fluid viscosity by comparing pressure profile of fluid/gas mixture with that of power-law fluid. Fluid is taken to be viscous, non-Newtonian, and incompressible; the gas to be ideal; the flow to be inertia-free, isothermal, and one dimensional. Analysis assists in design of flow systems for petroleum, coal, polymers, and other materials.

The Statistical Mechanics of Ideal Homogeneous Turbulence

NASA Technical Reports Server (NTRS)

Shebalin, John V.

2002-01-01

Plasmas, such as those found in the space environment or in plasma confinement devices, are often modeled as electrically conducting fluids. When fluids and plasmas are energetically stirred, regions of highly nonlinear, chaotic behavior known as turbulence arise. Understanding the fundamental nature of turbulence is a long-standing theoretical challenge. The present work describes a statistical theory concerning a certain class of nonlinear, finite dimensional, dynamical models of turbulence. These models arise when the partial differential equations describing incompressible, ideal (i.e., nondissipative) homogeneous fluid and magnetofluid (i.e., plasma) turbulence are Fourier transformed into a very large set of ordinary differential equations. These equations define a divergenceless flow in a high-dimensional phase space, which allows for the existence of a Liouville theorem, guaranteeing a distribution function based on constants of the motion (integral invariants). The novelty of these particular dynamical systems is that there are integral invariants other than the energy, and that some of these invariants behave like pseudoscalars under two of the discrete symmetry transformations of physics, parity, and charge conjugation. In this work the 'rugged invariants' of ideal homogeneous turbulence are shown to be the only significant scalar and pseudoscalar invariants. The discovery that pseudoscalar invariants cause symmetries of the original equations to be dynamically broken and induce a nonergodic structure on the associated phase space is the primary result presented here. Applicability of this result to dissipative turbulence is also discussed.

Quasi-isentropic compression using compressed water flow generated by underwater electrical explosion of a wire array

NASA Astrophysics Data System (ADS)

Gurovich, V.; Virozub, A.; Rososhek, A.; Bland, S.; Spielman, R. B.; Krasik, Ya. E.

2018-05-01

A major experimental research area in material equation-of-state today involves the use of off-Hugoniot measurements rather than shock experiments that give only Hugoniot data. There is a wide range of applications using quasi-isentropic compression of matter including the direct measurement of the complete isentrope of materials in a single experiment and minimizing the heating of flyer plates for high-velocity shock measurements. We propose a novel approach to generating quasi-isentropic compression of matter. Using analytical modeling and hydrodynamic simulations, we show that a working fluid composed of compressed water, generated by an underwater electrical explosion of a planar wire array, might be used to efficiently drive the quasi-isentropic compression of a copper target to pressures ˜2 × 1011 Pa without any complex target designs.

Broken Ergodicity in Ideal, Homogeneous, Incompressible Turbulence

NASA Technical Reports Server (NTRS)

Morin, Lee; Shebalin, John; Fu, Terry; Nguyen, Phu; Shum, Victor

2010-01-01

We discuss the statistical mechanics of numerical models of ideal homogeneous, incompressible turbulence and their relevance for dissipative fluids and magnetofluids. These numerical models are based on Fourier series and the relevant statistical theory predicts that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. However, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We explain this phenomena in terms of broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We review the theoretical basis of broken ergodicity, apply it to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence, and show new results from simulations using GPU (graphical processing unit) computers.

Noise Production of an Idealized Two-Dimensional Fish School

NASA Astrophysics Data System (ADS)

Wagenhoffer, Nathan; Moored, Keith; Jaworski, Justin

2017-11-01

The analysis of quiet bio-inspired propulsive concepts requires a rapid, unified computational framework that integrates the coupled fluid-solid dynamics of swimmers and their wakes with the resulting noise generation. Such a framework is presented for two-dimensional flows, where the fluid motion is modeled by an unsteady boundary element method with a vortex-particle wake. The unsteady surface forces from the potential flow solver are then passed to an acoustic boundary element solver to predict the radiated sound in low-Mach-number flows. The coupled flow-acoustic solver is validated against canonical vortex-sound problems. A diamond arrangement of four airfoils are subjected to traveling wave kinematics representing a known idealized pattern for a school of fish, and the airfoil motion and inflow values are derived from the range of Strouhal values common to many natural swimmers. The coupled flow-acoustic solver estimates and analyzes the hydrodynamic performance and noise production of the idealized school of swimmers.

Polarization and Compressibility of Oblique Kinetic Alfven Waves

NASA Technical Reports Server (NTRS)

Hunana, Peter; Goldstein, M. L.; Passot, T.; Sulem, P. L.; Laveder, D.; Zank, G. P.

2012-01-01

Even though solar wind, as a collisionless plasma, is properly described by the kineticMaxwell-Vlasov description, it can be argued that much of our understanding of solar wind observational data comes from an interpretation and numerical modeling which is based on a fluid description of magnetohydrodynamics. In recent years, there has been a significant interest in better understanding the importance of kinetic effects, i.e. the differences between the kinetic and usual fluid descriptions. Here we concentrate on physical properties of oblique kinetic Alfvn waves (KAWs), which are often recognized as one of the key ingredients in the solar wind turbulence cascade. We use three different fluid models with various degrees of complexity and calculate polarization and magnetic compressibility of oblique KAWs (propagation angle q = 88), which we compare to solutions derived from linear kinetic theory. We explore a wide range of possible proton plasma b = [0.1,10.0] and a wide range of length scales krL = [0.001,10.0]. It is shown that the classical isotropic two-fluid model is very compressible in comparison with kinetic theory and that the largest discrepancy occurs at scales larger than the proton gyroscale. We also show that the two-fluid model contains a large error in the polarization of electric field, even at scales krL 1. Furthermore, to understand these discrepancies between the two-fluid model and the kinetic theory, we employ two versions of the Landau fluid model that incorporate linear low-frequency kinetic effects such as Landau damping and finite Larmor radius (FLR) corrections into the fluid description. It is shown that Landau damping significantly reduces the magnetic compressibility and that FLR corrections (i.e. nongyrotropic contributions) are required to correctly capture the polarization.We also show that, in addition to Landau damping, FLR corrections are necessary to accurately describe the damping rate of KAWs. We conclude that kinetic effects

Bi-directionally draining pore fluid extraction vessel

DOEpatents

Prizio, Joseph; Ritt, Alexander; Mower, Timothy E.; Rodine, Lonn

1991-01-01

The invention is used to extract pore fluid from porous solids through a combination of mechanical compression and inert-gas injection and comprises a piston for axially compressing samples to force water out, and top and bottom drainage plates for capturing the exuded water and using inert gas to force water to exit when the limits of mechanical compression have been reached.

Density distribution function of a self-gravitating isothermal compressible turbulent fluid in the context of molecular clouds ensembles

NASA Astrophysics Data System (ADS)

Donkov, Sava; Stefanov, Ivan Z.

2018-03-01

We have set ourselves the task of obtaining the probability distribution function of the mass density of a self-gravitating isothermal compressible turbulent fluid from its physics. We have done this in the context of a new notion: the molecular clouds ensemble. We have applied a new approach that takes into account the fractal nature of the fluid. Using the medium equations, under the assumption of steady state, we show that the total energy per unit mass is an invariant with respect to the fractal scales. As a next step we obtain a non-linear integral equation for the dimensionless scale Q which is the third root of the integral of the probability distribution function. It is solved approximately up to the leading-order term in the series expansion. We obtain two solutions. They are power-law distributions with different slopes: the first one is -1.5 at low densities, corresponding to an equilibrium between all energies at a given scale, and the second one is -2 at high densities, corresponding to a free fall at small scales.

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.

Micromechanics of compression failures in open hole composite laminates

NASA Technical Reports Server (NTRS)

Guynn, E. Gail; Bradley, Walter L.

1987-01-01

The high strength-to-weight ratio of composite materials is ideally suited for aerospace applications where they already are used in commercial and military aircraft secondary structures and will soon be used for heavily loaded primary structures. One area impeding the widespread application of composites is their inherent weakness in compressive strength when compared to the tensile properties of the same material. Furthermore, these airframe designs typically contain many bolted or riveted joints, as well as electrical and hydraulic control lines. These applications produce areas of stress concentration, and thus, further complicate the compression failure problem. Open hole compression failures which represent a typical failure mode for composite materials are addressed.

Diffusion-driven fluid dynamics in ideal gases and plasmas

NASA Astrophysics Data System (ADS)

Vold, E. L.; Yin, L.; Taitano, W.; Molvig, K.; Albright, B. J.

2018-06-01

The classical transport theory based on Chapman-Enskog methods provides self-consistent approximations for the kinetic flux of mass, heat, and momentum in a fluid limit characterized with a small Knudsen number. The species mass fluxes relative to the center of mass, or "diffusive fluxes," are expressed as functions of known gradient quantities with kinetic coefficients evaluated using similar analyses for mixtures of gases or plasma components. The sum over species of the diffusive mass fluxes is constrained to be zero in the Lagrange frame, and thus results in a non-zero molar flux leading to a pressure perturbation. At an interface between two species initially in pressure equilibrium, the pressure perturbation driven by the diffusive molar flux induces a center of mass velocity directed from the species of greater atomic mass towards the lighter atomic mass species. As the ratio of the species particle masses increases, this center of mass velocity carries an increasingly greater portion of the mass across the interface and for a particle mass ratio greater than about two, the center of mass velocity carries more mass than the gradient driven diffusion flux. Early time transients across an interface between two species in a 1D plasma regime and initially in equilibrium are compared using three methods; a fluid code with closure in a classical transport approximation, a particle in cell simulation, and an implicit Fokker-Planck solver for the particle distribution functions. The early time transient phenomenology is shown to be similar in each of the computational simulation methods, including a pressure perturbation associated with the stationary "induced" component of the center of mass velocity which decays to pressure equilibrium during diffusion. At early times, the diffusive process generates pressure and velocity waves which propagate outward from the interface and are required to maintain momentum conservation. The energy in the outgoing waves dissipates

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.

COMPRESSIBLE FLOW, ENTRAINMENT, AND MEGAPLUME

EPA Science Inventory

It is generally believed that low Mach number, i.e., low-velocity, flow may be assumed to be incompressible flow. Under steady-state conditions, an exact equation of continuity may then be used to show that such flow is non-divergent. However, a rigorous, compressible fluid-dynam...

Proposal of a critical test of the Navier-Stokes-Fourier paradigm for compressible fluid continua.

PubMed

Brenner, Howard

2013-01-01

A critical, albeit simple experimental and/or molecular-dynamic (MD) simulation test is proposed whose outcome would, in principle, establish the viability of the Navier-Stokes-Fourier (NSF) equations for compressible fluid continua. The latter equation set, despite its longevity as constituting the fundamental paradigm of continuum fluid mechanics, has recently been criticized on the basis of its failure to properly incorporate volume transport phenomena-as embodied in the proposed bivelocity paradigm [H. Brenner, Int. J. Eng. Sci. 54, 67 (2012)]-into its formulation. Were the experimental or simulation results found to accord, even only qualitatively, with bivelocity predictions, the temperature distribution in a gas-filled, thermodynamically and mechanically isolated circular cylinder undergoing steady rigid-body rotation in an inertial reference frame would not be uniform; rather, the temperature would be higher at the cylinder wall than along the axis of rotation. This radial temperature nonuniformity contrasts with the uniformity of the temperature predicted by the NSF paradigm for these same circumstances. Easily attainable rates of rotation in centrifuges and readily available tools for measuring the expected temperature differences render experimental execution of the proposed scheme straightforward in principle. As such, measurement-via experiment or MD simulation-of, say, the temperature difference ΔT between the gas at the wall and along the axis of rotation would provide quantitative tests of both the NSF and bivelocity hydrodynamic models, whose respective solutions for the stated set of circumstances are derived in this paper. Independently of the correctness of the bivelocity model, any temperature difference observed during the proposed experiment or simulation, irrespective of magnitude, would preclude the possibility of the NSF paradigm being correct for fluid continua, except for incompressible flows.

An added-mass partition algorithm for fluid–structure interactions of compressible fluids and nonlinear solids

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

Banks, J.W., E-mail: banksj3@rpi.edu; Henshaw, W.D., E-mail: henshw@rpi.edu; Kapila, A.K., E-mail: kapila@rpi.edu

We describe an added-mass partitioned (AMP) algorithm for solving fluid–structure interaction (FSI) problems involving inviscid compressible fluids interacting with nonlinear solids that undergo large rotations and displacements. The computational approach is a mixed Eulerian–Lagrangian scheme that makes use of deforming composite grids (DCG) to treat large changes in the geometry in an accurate, flexible, and robust manner. The current work extends the AMP algorithm developed in Banks et al. [1] for linearly elasticity to the case of nonlinear solids. To ensure stability for the case of light solids, the new AMP algorithm embeds an approximate solution of a nonlinear fluid–solidmore » Riemann (FSR) problem into the interface treatment. The solution to the FSR problem is derived and shown to be of a similar form to that derived for linear solids: the state on the interface being fundamentally an impedance-weighted average of the fluid and solid states. Numerical simulations demonstrate that the AMP algorithm is stable even for light solids when added-mass effects are large. The accuracy and stability of the AMP scheme is verified by comparison to an exact solution using the method of analytical solutions and to a semi-analytical solution that is obtained for a rotating solid disk immersed in a fluid. The scheme is applied to the simulation of a planar shock impacting a light elliptical-shaped solid, and comparisons are made between solutions of the FSI problem for a neo-Hookean solid, a linearly elastic solid, and a rigid solid. The ability of the approach to handle large deformations is demonstrated for a problem of a high-speed flow past a light, thin, and flexible solid beam.« less

Negative Compressibility and Inverse Problem for Spinning Gas

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

Vasily Geyko and Nathaniel J. Fisch

2013-01-11

A spinning ideal gas in a cylinder with a smooth surface is shown to have unusual properties. First, under compression parallel to the axis of rotation, the spinning gas exhibits negative compressibility because energy can be stored in the rotation. Second, the spinning breaks the symmetry under which partial pressures of a mixture of gases simply add proportional to the constituent number densities. Thus, remarkably, in a mixture of spinning gases, an inverse problem can be formulated such that the gas constituents can be determined through external measurements only.

Adaptive mesh fluid simulations on GPU

NASA Astrophysics Data System (ADS)

Wang, Peng; Abel, Tom; Kaehler, Ralf

2010-10-01

We describe an implementation of compressible inviscid fluid solvers with block-structured adaptive mesh refinement on Graphics Processing Units using NVIDIA's CUDA. We show that a class of high resolution shock capturing schemes can be mapped naturally on this architecture. Using the method of lines approach with the second order total variation diminishing Runge-Kutta time integration scheme, piecewise linear reconstruction, and a Harten-Lax-van Leer Riemann solver, we achieve an overall speedup of approximately 10 times faster execution on one graphics card as compared to a single core on the host computer. We attain this speedup in uniform grid runs as well as in problems with deep AMR hierarchies. Our framework can readily be applied to more general systems of conservation laws and extended to higher order shock capturing schemes. This is shown directly by an implementation of a magneto-hydrodynamic solver and comparing its performance to the pure hydrodynamic case. Finally, we also combined our CUDA parallel scheme with MPI to make the code run on GPU clusters. Close to ideal speedup is observed on up to four GPUs.

The J-S model versus a non-ideal MHD theory

NASA Astrophysics Data System (ADS)

Franchi, Franca; Lazzari, Barbara; Nibbi, Roberta

2015-07-01

A new non-ideal electromagnetic interpretation of the J-S type viscoelastic model for polymeric fluids is given and a generalized resisto-elastic magnetohydrodynamic scenario for collisionless plasmas is proposed. The influence of the new theory on the incompressible transverse Alfvén waves is thoroughly investigated.

Methods for producing films using supercritical fluid

DOEpatents

Yonker, Clement R.; Fulton, John L.

2004-06-15

A method for forming a continuous film on a substrate surface that involves depositing particles onto a substrate surface and contacting the particle-deposited substrate surface with a supercritical fluid under conditions sufficient for forming a continuous film from the deposited particles. The particles may have a mean particle size of less 1 micron. The method may be performed by providing a pressure vessel that can contain a compressible fluid. A particle-deposited substrate is provided in the pressure vessel and the compressible fluid is maintained at a supercritical or sub-critical state sufficient for forming a film from the deposited particles. The T.sub.g of particles may be reduced by subjecting the particles to the methods detailed in the present disclosure.

Fluid Dynamics for Physicists

NASA Astrophysics Data System (ADS)

Faber, T. E.

1995-08-01

This textbook provides an accessible and comprehensive account of fluid dynamics that emphasizes fundamental physical principles and stresses connections with other branches of physics. Beginning with a basic introduction, the book goes on to cover many topics not typically treated in texts, such as compressible flow and shock waves, sound attenuation and bulk viscosity, solitary waves and ship waves, thermal convection, instabilities, turbulence, and the behavior of anisotropic, non-Newtonian and quantum fluids. Undergraduate or graduate students in physics or engineering who are taking courses in fluid dynamics will find this book invaluable.

Canonical fluid thermodynamics

NASA Technical Reports Server (NTRS)

Schmid, L. A.

1972-01-01

The space-time integral of the thermodynamic pressure plays the role of the thermodynamic potential for compressible, adiabatic flow in the sense that the pressure integral for stable flow is less than for all slightly different flows. This stability criterion can be converted into a variational minimum principle by requiring the molar free-enthalpy and the temperature, which are the arguments of the pressure function, to be generalized velocities, that is, the proper-time derivatives of scalar spare-time functions which are generalized coordinates in the canonical formalism. In a fluid context, proper-time differentiation must be expressed in terms of three independent quantities that specify the fluid velocity. This can be done in several ways, all of which lead to different variants (canonical transformations) of the same constraint-free action integral whose Euler-Lagrange equations are just the well-known equations of motion for adiabatic compressible flow.

Edge compression manifold apparatus

DOEpatents

Renzi, Ronald F.

2004-12-21

A manifold for connecting external capillaries to the inlet and/or outlet ports of a microfluidic device for high pressure applications is provided. The fluid connector for coupling at least one fluid conduit to a corresponding port of a substrate that includes: (i) a manifold comprising one or more channels extending therethrough wherein each channel is at least partially threaded, (ii) one or more threaded ferrules each defining a bore extending therethrough with each ferrule supporting a fluid conduit wherein each ferrule is threaded into a channel of the manifold, (iii) a substrate having one or more ports on its upper surface wherein the substrate is positioned below the manifold so that the one or more ports is aligned with the one or more channels of the manifold, and (iv) device to apply an axial compressive force to the substrate to couple the one or more ports of the substrate to a corresponding proximal end of a fluid conduit.

Edge compression manifold apparatus

DOEpatents

Renzi, Ronald F [Tracy, CA

2007-02-27

A manifold for connecting external capillaries to the inlet and/or outlet ports of a microfluidic device for high pressure applications is provided. The fluid connector for coupling at least one fluid conduit to a corresponding port of a substrate that includes: (i) a manifold comprising one or more channels extending therethrough wherein each channel is at least partially threaded, (ii) one or more threaded ferrules each defining a bore extending therethrough with each ferrule supporting a fluid conduit wherein each ferrule is threaded into a channel of the manifold, (iii) a substrate having one or more ports on its upper surface wherein the substrate is positioned below the manifold so that the one or more ports is aligned with the one or more channels of the manifold, and (iv) device to apply an axial compressive force to the substrate to couple the one or more ports of the substrate to a corresponding proximal end of a fluid conduit.

Electromotive force in strongly compressible magnetohydrodynamic turbulence

NASA Astrophysics Data System (ADS)

Yokoi, N.

2017-12-01

Variable density fluid turbulence is ubiquitous in geo-fluids, not to mention in astrophysics. Depending on the source of density variation, variable density fluid turbulence may be divided into two categories: the weak compressible (entropy mode) turbulence for slow flow and the strong compressible (acoustic mode) turbulence for fast flow. In the strong compressible turbulence, the pressure fluctuation induces a strong density fluctuation ρ ', which is represented by the density variance <ρ'2> (<·> denotes the ensemble average). The turbulent effect on the large-scale magnetic-field B induction is represented by the turbulent electromotive force (EMF) (u': velocity fluctuation, b': magnetic-field fluctuation). In the usual treatment in the dynamo theory, the expression for the EMF has been obtained in the framework of incompressible or weak compressible turbulence, where only the variation of the mean density <ρ>, if any, is taken into account. We see from the equation of the density fluctuation ρ', the density variance <ρ'2> is generated by the large mean density variation ∂<ρ> coupled with the turbulent mass flux <ρ'u'>. This means that in the region where the mean density steeply changes, the density variance effect becomes relevant for the magnetic field evolution. This situation is typically the case for phenomena associated with shocks and compositional discontinuities. With the aid of the analytical theory of inhomogeneous compressible magnetohydrodynamic (MHD) turbulence, the expression for the turbulent electromotive force is investigated. It is shown that, among others, an obliqueness (misalignment) between the mean density gradient ∂<ρ> and the mean magnetic field B may contribute to the EMF as ≈χ B×∂<ρ> with the turbulent transport coefficient χ proportional to the density variance (χ <ρ'2>). This density variance effect is expected to strongly affect the EMF near the interface, and changes the transport

Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations

NASA Astrophysics Data System (ADS)

Derigs, Dominik; Winters, Andrew R.; Gassner, Gregor J.; Walch, Stefanie; Bohm, Marvin

2018-07-01

The paper presents two contributions in the context of the numerical simulation of magnetized fluid dynamics. First, we show how to extend the ideal magnetohydrodynamics (MHD) equations with an inbuilt magnetic field divergence cleaning mechanism in such a way that the resulting model is consistent with the second law of thermodynamics. As a byproduct of these derivations, we show that not all of the commonly used divergence cleaning extensions of the ideal MHD equations are thermodynamically consistent. Secondly, we present a numerical scheme obtained by constructing a specific finite volume discretization that is consistent with the discrete thermodynamic entropy. It includes a mechanism to control the discrete divergence error of the magnetic field by construction and is Galilean invariant. We implement the new high-order MHD solver in the adaptive mesh refinement code FLASH where we compare the divergence cleaning efficiency to the constrained transport solver available in FLASH (unsplit staggered mesh scheme).

Geometry-driven distributed compression of the plenoptic function: performance bounds and constructive algorithms.

PubMed

Gehrig, Nicolas; Dragotti, Pier Luigi

2009-03-01

In this paper, we study the sampling and the distributed compression of the data acquired by a camera sensor network. The effective design of these sampling and compression schemes requires, however, the understanding of the structure of the acquired data. To this end, we show that the a priori knowledge of the configuration of the camera sensor network can lead to an effective estimation of such structure and to the design of effective distributed compression algorithms. For idealized scenarios, we derive the fundamental performance bounds of a camera sensor network and clarify the connection between sampling and distributed compression. We then present a distributed compression algorithm that takes advantage of the structure of the data and that outperforms independent compression algorithms on real multiview images.

On compressible and piezo-viscous flow in thin porous media.

PubMed

Pérez-Ràfols, F; Wall, P; Almqvist, A

2018-01-01

In this paper, we study flow through thin porous media as in, e.g. seals or fractures. It is often useful to know the permeability of such systems. In the context of incompressible and iso-viscous fluids, the permeability is the constant of proportionality relating the total flow through the media to the pressure drop. In this work, we show that it is also relevant to define a constant permeability when compressible and/or piezo-viscous fluids are considered. More precisely, we show that the corresponding nonlinear equation describing the flow of any compressible and piezo-viscous fluid can be transformed into a single linear equation. Indeed, this linear equation is the same as the one describing the flow of an incompressible and iso-viscous fluid. By this transformation, the total flow can be expressed as the product of the permeability and a nonlinear function of pressure, which represents a generalized pressure drop.

Physics-Based Computational Algorithm for the Multi-Fluid Plasma Model

DTIC Science & Technology

2014-06-30

applying it to study laser - 20 Physics-Based Multi-Fluid Plasma Algorithm Shumlak Figure 6: Blended finite element method applied to the species...separation problem in capsule implosions. Number densities and electric field are shown after the laser drive has compressed the multi-fluid plasma and...6 after the laser drive has started the compression. A separation clearly develops. The solution is found using an explicit advance (CFL=1) for the

Heterogeneous Compression of Large Collections of Evolutionary Trees.

PubMed

Matthews, Suzanne J

2015-01-01

Compressing heterogeneous collections of trees is an open problem in computational phylogenetics. In a heterogeneous tree collection, each tree can contain a unique set of taxa. An ideal compression method would allow for the efficient archival of large tree collections and enable scientists to identify common evolutionary relationships over disparate analyses. In this paper, we extend TreeZip to compress heterogeneous collections of trees. TreeZip is the most efficient algorithm for compressing homogeneous tree collections. To the best of our knowledge, no other domain-based compression algorithm exists for large heterogeneous tree collections or enable their rapid analysis. Our experimental results indicate that TreeZip averages 89.03 percent (72.69 percent) space savings on unweighted (weighted) collections of trees when the level of heterogeneity in a collection is moderate. The organization of the TRZ file allows for efficient computations over heterogeneous data. For example, consensus trees can be computed in mere seconds. Lastly, combining the TreeZip compressed (TRZ) file with general-purpose compression yields average space savings of 97.34 percent (81.43 percent) on unweighted (weighted) collections of trees. Our results lead us to believe that TreeZip will prove invaluable in the efficient archival of tree collections, and enables scientists to develop novel methods for relating heterogeneous collections of trees.

Simulating compressible-incompressible two-phase flows

NASA Astrophysics Data System (ADS)

Denner, Fabian; van Wachem, Berend

2017-11-01

Simulating compressible gas-liquid flows, e.g. air-water flows, presents considerable numerical issues and requires substantial computational resources, particularly because of the stiff equation of state for the liquid and the different Mach number regimes. Treating the liquid phase (low Mach number) as incompressible, yet concurrently considering the gas phase (high Mach number) as compressible, can improve the computational performance of such simulations significantly without sacrificing important physical mechanisms. A pressure-based algorithm for the simulation of two-phase flows is presented, in which a compressible and an incompressible fluid are separated by a sharp interface. The algorithm is based on a coupled finite-volume framework, discretised in conservative form, with a compressive VOF method to represent the interface. The bulk phases are coupled via a novel acoustically-conservative interface discretisation method that retains the acoustic properties of the compressible phase and does not require a Riemann solver. Representative test cases are presented to scrutinize the proposed algorithm, including the reflection of acoustic waves at the compressible-incompressible interface, shock-drop interaction and gas-liquid flows with surface tension. Financial support from the EPSRC (Grant EP/M021556/1) is gratefully acknowledged.

Dynamics and stability of a 2D ideal vortex under external strain

NASA Astrophysics Data System (ADS)

Hurst, N. C.; Danielson, J. R.; Dubin, D. H. E.; Surko, C. M.

2017-11-01

The behavior of an initially axisymmetric 2D ideal vortex under an externally imposed strain flow is studied experimentally. The experiments are carried out using electron plasmas confined in a Penning-Malmberg trap; here, the dynamics of the plasma density transverse to the field are directly analogous to the dynamics of vorticity in a 2D ideal fluid. An external strain flow is applied using boundary conditions in a way that is consistent with 2D fluid dynamics. Data are compared to predictions from a theory assuming a piecewise constant elliptical vorticity distribution. Excellent agreement is found for quasi-flat profiles, whereas the dynamics of smooth profiles feature modified stability limits and inviscid damping of periodic elliptical distortions. This work supported by U.S. DOE Grants DE-SC0002451 and DE-SC0016532, and NSF Grant PHY-1414570.

Initial eccentricity and constituent quark number scaling of elliptic flow in ideal and viscous dynamics

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

Chaudhuri, A. K.

2010-04-15

In the Israel-Stewart theory of dissipative hydrodynamics, the scaling properties of elliptic flow in Au+Au collisions are studied. The initial energy density of the fluid was fixed to reproduce STAR data on phi-meson multiplicity in 0-5% Au+Au collisions such that, irrespective of fluid viscosity, entropy at the freeze-out is similar in ideal or in viscous evolution. The initial eccentricity or constituent quark number scaling is only approximate in ideal or minimally viscous (eta/s=1/4pi) fluid. Eccentricity scaling becomes nearly exact in more viscous fluid (eta/s>=0.12). However, in more viscous fluid, constituent quark number scaled elliptic flow for mesons and baryons splitsmore » into separate scaling functions. Simulated flows also do not exhibit 'universal scaling'; that is, elliptic flow scaled by the constituent quark number and charged particles v{sub 2} is not a single function of transverse kinetic energy scaled by the quark number. From a study of the violation of universal scaling, we obtain an estimate of quark-gluon plasma viscosity, eta/s=0.12+-0.03. The error is statistical only. The systematic error in eta/s could be as large.« less

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.

Anomalous sorption of supercritical fluids on polymer thin films.

PubMed

Wang, Xiaochu; Sanchez, Isaac C

2006-10-24

Unusual sorption has been reported in thin polymer films exposed to near-critical CO2. When the supercritical fluid approaches the critical point, the film appears to thicken, but it is not clear whether the film swells or there is an adsorption layer on the film surface. A combination of the gradient theory of inhomogeneous systems and the Sanchez-Lacombe equation of state has been used to investigate this phenomenon. It is shown analytically that surface adsorption on an attractive surface is proportional to the compressibility of the fluid. We have also investigated numerically the sorption of supercritical CO2 on poly(dimethylsiloxane) and polyisobutylene, and supercritical 1,1-difluoroethane on polystyrene. By calculating the Gibbs adsorption and adsorption layer thickness of the supercritical fluids, we found in all cases (different substrates, different supercritical fluids) that maximum adsorption occurs when the supercritical fluid is near its compressibility maximum.

Effect of initial densities in the lattice Boltzmann model for non-ideal fluid with curved interface

NASA Astrophysics Data System (ADS)

Gong, Jiaming; Oshima, Nobuyuki

2017-06-01

The effect of initial densities in a free energy based two-phase-flow lattice Boltzmann method for non-ideal fluids with a curved interface was investigated in the present work. To investigate this effect, the initial densities in the liquid and gas phases coming from the saturation points and the equilibrium state were adopted in the simulation of a static droplet in an open and a closed system. For the purpose of simplicity and easier comparison, the closed system is fabricated by the implementation of the periodic boundary condition at the inlet and outlet of a gas channel, and the open system is fabricated by the implementation of a constant flux boundary condition at the inlet and a free-out boundary condition at the outlet of the same gas channel. By comparing the simulation results from the two types of initial densities in the open and closed systems, it is proven that the commonly used saturation initial densities setting is the reason for droplet mass and volume variation which occurred in the simulation, particularly in the open system with a constant flux boundary condition. Such problems are believed to come from the curvature effect of the surface tension and can be greatly reduced by adopting the initial densities in the two phases from equilibrium state.

Microfluidic pressure sensing using trapped air compression.

PubMed

Srivastava, Nimisha; Burns, Mark A

2007-05-01

We have developed a microfluidic method for measuring the fluid pressure head experienced at any location inside a microchannel. The principal component is a microfabricated sealed chamber with a single inlet and no exit; the entrance to the single inlet is positioned at the location where pressure is to be measured. The pressure measurement is then based on monitoring the movement of a liquid-air interface as it compresses air trapped inside the microfabricated sealed chamber and calculating the pressure using the ideal gas law. The method has been used to measure the pressure of the air stream and continuous liquid flow inside microfluidic channels (d approximately 50 microm). Further, a pressure drop has also been measured using multiple microfabricated sealed chambers. For air pressure, a resolution of 700 Pa within a full-scale range of 700-100 kPa was obtained. For liquids, pressure drops as low as 70 Pa were obtained in an operating range from 70 Pa to 10 kPa. Since the method primarily uses a microfluidic sealed chamber, it does not require additional fabrication steps and may easily be incorporated in several lab-on-a-chip fluidic applications for laminar as well as turbulent flow conditions.

Microfluidic pressure sensing using trapped air compression

PubMed Central

Srivastava, Nimisha; Burns, Mark A.

2010-01-01

We have developed a microfluidic method for measuring the fluid pressure head experienced at any location inside a microchannel. The principal component is a microfabricated sealed chamber with a single inlet and no exit; the entrance to the single inlet is positioned at the location where pressure is to be measured. The pressure measurement is then based on monitoring the movement of a liquid–air interface as it compresses air trapped inside the microfabricated sealed chamber and calculating the pressure using the ideal gas law. The method has been used to measure the pressure of the air stream and continuous liquid flow inside microfluidic channels (d ~ 50 μm). Further, a pressure drop has also been measured using multiple microfabricated sealed chambers. For air pressure, a resolution of 700 Pa within a full-scale range of 700–100 kPa was obtained. For liquids, pressure drops as low as 70 Pa were obtained in an operating range from 70 Pa to 10 kPa. Since the method primarily uses a microfluidic sealed chamber, it does not require additional fabrication steps and may easily be incorporated in several lab-on-a-chip fluidic applications for laminar as well as turbulent flow conditions. PMID:17476384

On the Flow of a Compressible Fluid by the Hodograph Method. II - Fundamental Set of Particular Flow Solutions of the Chaplygin Differential Equation

NASA Technical Reports Server (NTRS)

Garrick, I. E.; Kaplan, Carl

1944-01-01

The differential equation of Chaplygin's jet problem is utilized to give a systematic development of particular solutions of the hodograph flow equations, which extends the treatment of Chaplygin into the supersonic range and completes the set of particular solutions. The particular solutions serve to place on a reasonable basis the use of velocity correction formulas for the comparison of incompressible and compressible flows. It is shown that the geometric-mean type of velocity correction formula introduced in part I has significance as an over-all type of approximation in the subsonic range. A brief review of general conditions limiting the potential flow of an adiabatic compressible fluid is given and application is made to the particular solutions, yielding conditions for the existence of singular loci in the supersonic range. The combining of particular solutions in accordance with prescribed boundary flow conditions is not treated in the present paper.

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.