Fluid absorption solar energy receiver

NASA Technical Reports Server (NTRS)

Bair, Edward J.

1993-01-01

A conventional solar dynamic system transmits solar energy to the flowing fluid of a thermodynamic cycle through structures which contain the gas and thermal energy storage material. Such a heat transfer mechanism dictates that the structure operate at a higher temperature than the fluid. This investigation reports on a fluid absorption receiver where only a part of the solar energy is transmitted to the structure. The other part is absorbed directly by the fluid. By proportioning these two heat transfer paths the energy to the structure can preheat the fluid, while the energy absorbed directly by the fluid raises the fluid to its final working temperature. The surface temperatures need not exceed the output temperature of the fluid. This makes the output temperature of the gas the maximum temperature in the system. The gas can have local maximum temperatures higher than the output working temperature. However local high temperatures are quickly equilibrated, and since the gas does not emit radiation, local high temperatures do not result in a radiative heat loss. Thermal radiation, thermal conductivity, and heat exchange with the gas all help equilibrate the surface temperature.

A Spalart-Allmaras local correlation-based transition model for Thermo-fuid dynamics

NASA Astrophysics Data System (ADS)

D'Alessandro, V.; Garbuglia, F.; Montelpare, S.; Zoppi, A.

2017-11-01

The study of innovative energy systems often involves complex fluid flows problems and the Computational Fluid-Dynamics (CFD) is one of the main tools of analysis. It is important to put in evidence that in several energy systems the flow field experiences the laminar-to-turbulent transition. Direct Numerical Simulations (DNS) or Large Eddy Simulation (LES) are able to predict the flow transition but they are still inapplicable to the study of real problems due to the significant computational resources requirements. Differently standard Reynolds Averaged Navier Stokes (RANS) approaches are not always reliable since they assume a fully turbulent regime. In order to overcome this drawback in the recent years some locally formulated transition RANS models have been developed. In this work, we present a local correlation-based transition approach adding two equations that control the laminar-toturbulent transition process -γ and \\[\\overset{}{\\mathop{{{\\operatorname{Re}}θ, \\text{t}}}} \\] - to the well-known Spalart-Allmaras (SA) turbulence model. The new model was implemented within OpenFOAM code. The energy equation is also implemented in order to evaluate the model performance in thermal-fluid dynamics applications. In all the considered cases a very good agreement between numerical and experimental data was observed.

Statistical Decoupling of a Lagrangian Fluid Parcel in Newtonian Cosmology

NASA Astrophysics Data System (ADS)

Wang, Xin; Szalay, Alex

2016-03-01

The Lagrangian dynamics of a single fluid element within a self-gravitational matter field is intrinsically non-local due to the presence of the tidal force. This complicates the theoretical investigation of the nonlinear evolution of various cosmic objects, e.g., dark matter halos, in the context of Lagrangian fluid dynamics, since fluid parcels with given initial density and shape may evolve differently depending on their environments. In this paper, we provide a statistical solution that could decouple this environmental dependence. After deriving the evolution equation for the probability distribution of the matter field, our method produces a set of closed ordinary differential equations whose solution is uniquely determined by the initial condition of the fluid element. Mathematically, it corresponds to the projected characteristic curve of the transport equation of the density-weighted probability density function (ρPDF). Consequently it is guaranteed that the one-point ρPDF would be preserved by evolving these local, yet nonlinear, curves with the same set of initial data as the real system. Physically, these trajectories describe the mean evolution averaged over all environments by substituting the tidal tensor with its conditional average. For Gaussian distributed dynamical variables, this mean tidal tensor is simply proportional to the velocity shear tensor, and the dynamical system would recover the prediction of the Zel’dovich approximation (ZA) with the further assumption of the linearized continuity equation. For a weakly non-Gaussian field, the averaged tidal tensor could be expanded perturbatively as a function of all relevant dynamical variables whose coefficients are determined by the statistics of the field.

STATISTICAL DECOUPLING OF A LAGRANGIAN FLUID PARCEL IN NEWTONIAN COSMOLOGY

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

Wang, Xin; Szalay, Alex, E-mail: xwang@cita.utoronto.ca

The Lagrangian dynamics of a single fluid element within a self-gravitational matter field is intrinsically non-local due to the presence of the tidal force. This complicates the theoretical investigation of the nonlinear evolution of various cosmic objects, e.g., dark matter halos, in the context of Lagrangian fluid dynamics, since fluid parcels with given initial density and shape may evolve differently depending on their environments. In this paper, we provide a statistical solution that could decouple this environmental dependence. After deriving the evolution equation for the probability distribution of the matter field, our method produces a set of closed ordinary differentialmore » equations whose solution is uniquely determined by the initial condition of the fluid element. Mathematically, it corresponds to the projected characteristic curve of the transport equation of the density-weighted probability density function (ρPDF). Consequently it is guaranteed that the one-point ρPDF would be preserved by evolving these local, yet nonlinear, curves with the same set of initial data as the real system. Physically, these trajectories describe the mean evolution averaged over all environments by substituting the tidal tensor with its conditional average. For Gaussian distributed dynamical variables, this mean tidal tensor is simply proportional to the velocity shear tensor, and the dynamical system would recover the prediction of the Zel’dovich approximation (ZA) with the further assumption of the linearized continuity equation. For a weakly non-Gaussian field, the averaged tidal tensor could be expanded perturbatively as a function of all relevant dynamical variables whose coefficients are determined by the statistics of the field.« less

Real-Time Maps of Fluid Flow Fields in Porous Biomaterials

PubMed Central

Mack, Julia J.; Youssef, Khalid; Noel, Onika D.V.; Lake, Michael P.; Wu, Ashley; Iruela-Arispe, M. Luisa; Bouchard, Louis-S.

2013-01-01

Mechanical forces such as fluid shear have been shown to enhance cell growth and differentiation, but knowledge of their mechanistic effect on cells is limited because the local flow patterns and associated metrics are not precisely known. Here we present real-time, noninvasive measures of local hydrodynamics in 3D biomaterials based on nuclear magnetic resonance. Microflow maps were further used to derive pressure, shear and fluid permeability fields. Finally, remodeling of collagen gels in response to precise fluid flow parameters was correlated with structural changes. It is anticipated that accurate flow maps within 3D matrices will be a critical step towards understanding cell behavior in response to controlled flow dynamics. PMID:23245922

Influence of mechanical rock properties and fracture healing rate on crustal fluid flow dynamics

NASA Astrophysics Data System (ADS)

Sachau, Till; Bons, Paul; Gomez-Rivas, Enrique; Koehn, Daniel; de Riese, Tamara

2016-04-01

Fluid flow in the Earth's crust is very slow over extended periods of time, during which it occurs within the connected pore space of rocks. If the fluid production rate exceeds a certain threshold, matrix permeability alone is insufficient to drain the fluid volume and fluid pressure builds up, thereby reducing the effective stress supported by the rock matrix. Hydraulic fractures form once the effective pressure exceeds the tensile strength of the rock matrix and act subsequently as highly effective fluid conduits. Once local fluid pressure is sufficiently low again, flow ceases and fractures begin to heal. Since fluid flow is controlled by the alternation of fracture permeability and matrix permeability, the flow rate in the system is strongly discontinuous and occurs in intermittent pulses. Resulting hydraulic fracture networks are largely self-organized: opening and subsequent healing of hydraulic fractures depends on the local fluid pressure and on the time-span between fluid pulses. We simulate this process with a computer model and describe the resulting dynamics statistically. Special interest is given to a) the spatially and temporally discontinuous formation and closure of fractures and fracture networks and b) the total flow rate over time. The computer model consists of a crustal-scale dual-porosity setup. Control parameters are the pressure- and time-dependent fracture healing rate, and the strength and the permeability of the intact rock. Statistical analysis involves determination of the multifractal properties and of the power spectral density of the temporal development of the total drainage rate and hydraulic fractures. References Bons, P. D. (2001). The formation of large quartz veins by rapid ascent of fluids in mobile hydrofractures. Tectonophysics, 336, 1-17. Miller, S. a., & Nur, A. (2000). Permeability as a toggle switch in fluid-controlled crustal processes. Earth and Planetary Science Letters, 183(1-2), 133-146. Sachau, T., Bons, P. D., & Gomez-Rivas, E. (2015). Transport efficiency and dynamics of hydraulic fracture networks. Frontiers in Physics, 3.

Stalactite growth as a free-boundary problem: a geometric law and its platonic ideal.

PubMed

Short, Martin B; Baygents, James C; Beck, J Warren; Stone, David A; Toomey, Rickard S; Goldstein, Raymond E

2005-01-14

The chemical mechanisms underlying the growth of cave formations such as stalactites are well known, yet no theory has yet been proposed which successfully accounts for the dynamic evolution of their shapes. Here we consider the interplay of thin-film fluid dynamics, calcium carbonate chemistry, and CO2 transport in the cave to show that stalactites evolve according to a novel local geometric growth law which exhibits extreme amplification at the tip as a consequence of the locally-varying fluid layer thickness. Studies of this model show that a broad class of initial conditions is attracted to an ideal shape which is strikingly close to a statistical average of natural stalactites.

Domain decomposition methods in computational fluid dynamics

NASA Technical Reports Server (NTRS)

Gropp, William D.; Keyes, David E.

1991-01-01

The divide-and-conquer paradigm of iterative domain decomposition, or substructuring, has become a practical tool in computational fluid dynamic applications because of its flexibility in accommodating adaptive refinement through locally uniform (or quasi-uniform) grids, its ability to exploit multiple discretizations of the operator equations, and the modular pathway it provides towards parallelism. These features are illustrated on the classic model problem of flow over a backstep using Newton's method as the nonlinear iteration. Multiple discretizations (second-order in the operator and first-order in the preconditioner) and locally uniform mesh refinement pay dividends separately, and they can be combined synergistically. Sample performance results are included from an Intel iPSC/860 hypercube implementation.

Water pumping in mantle shear zones

PubMed Central

Précigout, Jacques; Prigent, Cécile; Palasse, Laurie; Pochon, Anthony

2017-01-01

Water plays an important role in geological processes. Providing constraints on what may influence the distribution of aqueous fluids is thus crucial to understanding how water impacts Earth's geodynamics. Here we demonstrate that ductile flow exerts a dynamic control on water-rich fluid circulation in mantle shear zones. Based on amphibole distribution and using dislocation slip-systems as a proxy for syn-tectonic water content in olivine, we highlight fluid accumulation around fine-grained layers dominated by grain-size-sensitive creep. This fluid aggregation correlates with dislocation creep-accommodated strain that localizes in water-rich layers. We also give evidence of cracking induced by fluid pressure where the highest amount of water is expected. These results emphasize long-term fluid pumping attributed to creep cavitation and associated phase nucleation during grain size reduction. Considering the ubiquitous process of grain size reduction during strain localization, our findings shed light on multiple fluid reservoirs in the crust and mantle. PMID:28593947

Experimental Observations of Multiscale Dynamics of Viscous Fluid Behavior: Implications in Volcanic Systems

NASA Astrophysics Data System (ADS)

Arciniega-Ceballos, A.; Spina, L.; Scheu, B.; Dingwell, D. B.

2015-12-01

We have investigated the dynamics of Newtonian fluids with viscosities (10-1000 Pa s; corresponding to mafic to intermediate silicate melts) during slow decompression, in a Plexiglas shock tube. As an analogue fluid we used silicon oil saturated with Argon gas for 72 hours. Slow decompression, dropping from 10 MPa to ambient pressure, acts as the excitation mechanism, initiating several processes with their own distinct timescales. The evolution of this multi-timescale phenomenon generates complex non-stationary microseismic signals, which have been recorded with 7 high-dynamic piezoelectric sensors located along the conduit. Correlation analysis of these time series with the associated high-speed imaging enables characterization of distinct phases of the dynamics of these viscous fluids and the extraction of the time and the frequency characteristics of the individual processes. We have identified fluid-solid elastic interaction, degassing, fluid mass expansion and flow, bubble nucleation, growth, coalescence and collapse, foam building and vertical wagging. All these processes (in fine and coarse scales) are sequentially coupled in time, occur within specific pressure intervals, and exhibit a localized distribution in space. Their coexistence and interactions constitute the stress field and driving forces that determine the dynamics of the system. Our observations point to the great potential of this experimental approach in the understanding of volcanic processes and volcanic seismicity.

Episodic Tremor and Slip Explained by Fluid-Enhanced Microfracturing and Sealing

NASA Astrophysics Data System (ADS)

Bernaudin, M.; Gueydan, F.

2018-04-01

Episodic tremor and slow-slip events at the deep extension of plate boundary faults illuminate seismic to aseismic processes around the brittle-ductile transition. These events occur in volumes characterized by overpressurized fluids and by near failure shear stress conditions. We present a new modeling approach based on a ductile grain size-sensitive rheology with microfracturing and sealing, which provides a mechanical and field-based explanation of such phenomena. We also model pore fluid pressure variation as a function of changes in porosity/permeability and strain rate-dependent fluid pumping. The fluid-enhanced dynamic evolution of microstructures defines cycles of ductile strain localization and implies increase in pore fluid pressure. We propose that slow-slip events are ductile processes related to transient strain localization, while nonvolcanic tremor corresponds to fracturing of the whole rock at the peak of pore fluid pressure. Our model shows that the availability of fluids and the efficiency of fluid pumping control the occurrence and the P-T conditions of episodic tremor and slip.

Effective field theory of dissipative fluids (II): classical limit, dynamical KMS symmetry and entropy current

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

Glorioso, Paolo; Crossley, Michael; Liu, Hong

2017-09-20

Here in this paper we further develop the fluctuating hydrodynamics proposed in a number of ways. We first work out in detail the classical limit of the hydrodynamical action, which exhibits many simplifications. In particular, this enables a transparent formulation of the action in physical spacetime in the presence of arbitrary external fields. It also helps to clarify issues related to field redefinitions and frame choices. We then propose that the action is invariant under a Z2 symmetry to which we refer as the dynamical KMS symmetry. The dynamical KMS symmetry is physically equivalent to the previously proposed local KMSmore » condition in the classical limit, but is more convenient to implement and more general. It is applicable to any states in local equilibrium rather than just thermal density matrix perturbed by external background fields. Finally we elaborate the formulation for a conformal fluid, which contains some new features, and work out the explicit form of the entropy current to second order in derivatives for a neutral conformal fluid.« less

Energy dissipation in flows through curved spaces.

PubMed

Debus, J-D; Mendoza, M; Succi, S; Herrmann, H J

2017-02-14

Fluid dynamics in intrinsically curved geometries is encountered in many physical systems in nature, ranging from microscopic bio-membranes all the way up to general relativity at cosmological scales. Despite the diversity of applications, all of these systems share a common feature: the free motion of particles is affected by inertial forces originating from the curvature of the embedding space. Here we reveal a fundamental process underlying fluid dynamics in curved spaces: the free motion of fluids, in the complete absence of solid walls or obstacles, exhibits loss of energy due exclusively to the intrinsic curvature of space. We find that local sources of curvature generate viscous stresses as a result of the inertial forces. The curvature- induced viscous forces are shown to cause hitherto unnoticed and yet appreciable energy dissipation, which might play a significant role for a variety of physical systems involving fluid dynamics in curved spaces.

Self-regulation in self-propelled nematic fluids.

PubMed

Baskaran, A; Marchetti, M C

2012-09-01

We consider the hydrodynamic theory of an active fluid of self-propelled particles with nematic aligning interactions. This class of materials has polar symmetry at the microscopic level, but forms macrostates of nematic symmetry. We highlight three key features of the dynamics. First, as in polar active fluids, the control parameter for the order-disorder transition, namely the density, is dynamically convected by the order parameter via active currents. The resulting dynamical self-regulation of the order parameter is a generic property of active fluids and destabilizes the uniform nematic state near the mean-field transition. Secondly, curvature-driven currents render the system unstable deep in the nematic state, as found previously. Finally, and unique to self-propelled nematics, nematic order induces local polar order that in turn leads to the growth of density fluctuations. We propose this as a possible mechanism for the smectic order of polar clusters seen in numerical simulations.

Beyond the standard two-film theory: Computational fluid dynamics simulations for carbon dioxide capture in a wetted wall column

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

Wang, Chao; Xu, Zhijie; Lai, Canhai

The standard two-film theory (STFT) is a diffusion-based mechanism that can be used to describe gas mass transfer across liquid film. Fundamental assumptions of the STFT impose serious limitations on its ability to predict mass transfer coefficients. To better understand gas absorption across liquid film in practical situations, a multiphase computational fluid dynamics (CFD) model fully equipped with mass transport and chemistry capabilities has been developed for solvent-based carbon dioxide (CO 2) capture to predict the CO 2 mass transfer coefficient in a wetted wall column. The hydrodynamics is modeled using a volume of fluid method, and the diffusive andmore » reactive mass transfer between the two phases is modeled by adopting a one-fluid formulation. We demonstrate that the proposed CFD model can naturally account for the influence of many important factors on the overall mass transfer that cannot be quantitatively explained by the STFT, such as the local variation in fluid velocities and properties, flow instabilities, and complex geometries. The CFD model also can predict the local mass transfer coefficient variation along the column height, which the STFT typically does not consider.« less

Beyond the standard two-film theory: Computational fluid dynamics simulations for carbon dioxide capture in a wetted wall column

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

Wang, Chao; Xu, Zhijie; Lai, Canhai

The standard two-film theory (STFT) is a diffusion-based mechanism that can be used to describe gas mass transfer across liquid film. Fundamental assumptions of the STFT impose serious limitations on its ability to predict mass transfer coefficients. To better understand gas absorption across liquid film in practical situations, a multiphase computational fluid dynamics (CFD) model fully equipped with mass transport and chemistry capabilities has been developed for solvent-based carbon dioxide (CO2) capture to predict the CO2 mass transfer coefficient in a wetted wall column. The hydrodynamics is modeled using a volume of fluid method, and the diffusive and reactive massmore » transfer between the two phases is modeled by adopting a one-fluid formulation. We demonstrate that the proposed CFD model can naturally account for the influence of many important factors on the overall mass transfer that cannot be quantitatively explained by the STFT, such as the local variation in fluid velocities and properties, flow instabilities, and complex geometries. The CFD model also can predict the local mass transfer coefficient variation along the column height, which the STFT typically does not consider.« less

Beyond the standard two-film theory: Computational fluid dynamics simulations for carbon dioxide capture in a wetted wall column

DOE PAGES

Wang, Chao; Xu, Zhijie; Lai, Canhai; ...

2018-03-27

The standard two-film theory (STFT) is a diffusion-based mechanism that can be used to describe gas mass transfer across liquid film. Fundamental assumptions of the STFT impose serious limitations on its ability to predict mass transfer coefficients. To better understand gas absorption across liquid film in practical situations, a multiphase computational fluid dynamics (CFD) model fully equipped with mass transport and chemistry capabilities has been developed for solvent-based carbon dioxide (CO 2) capture to predict the CO 2 mass transfer coefficient in a wetted wall column. The hydrodynamics is modeled using a volume of fluid method, and the diffusive andmore » reactive mass transfer between the two phases is modeled by adopting a one-fluid formulation. We demonstrate that the proposed CFD model can naturally account for the influence of many important factors on the overall mass transfer that cannot be quantitatively explained by the STFT, such as the local variation in fluid velocities and properties, flow instabilities, and complex geometries. The CFD model also can predict the local mass transfer coefficient variation along the column height, which the STFT typically does not consider.« less

Dynamic, Hot Surface Ignition of Aircraft Fuels and Hydraulic Fluids

DTIC Science & Technology

1980-10-01

fuels on a heated stainless steel surface. Higher local surface air speeds necessitated higher surface temperatures for ignition of an applied fluid._-7...Aircraft Fuels ( stainless steel surface) 8. Air Speed and Surface Material Effects on Hot Surface 21 Ignition Temperature of Aircraft Fuels (Titanium...Material Effects on Hot Surface 26 Ignition Temperature of Aircraft Hydraulic Fluids ( Stainless steel surface) 11. Air Speed and Surface Material

Localization and Instability in Sheared Granular Materials: Role of Pore Fluids and Non-monotonic Rate Dependent Rheology

NASA Astrophysics Data System (ADS)

Ma, X.; Elbanna, A. E.; Kothari, K.

2017-12-01

Fault zone dynamics hold the key to resolving many outstanding geophysical problems including the heat flow paradox, discrepancy between fault static and dynamic strength, and energy partitioning. Most fault zones that generate tectonic events are gouge filled and fluid saturated posing the need for formulating gouge-specific constitutive models that capture spatially heterogeneous compaction and dilation, non-monotonic rate dependence, and transition between localized and distributed deformation. In this presentation, we focus primarily on elucidating microscopic underpinnings for shear banding and stick-slip instabilities in sheared saturated granular materials and explore their implications for earthquake dynamics. We use a non-equilibrium thermodynamics model, the Shear Transformation Zone theory, to investigate the dynamics of strain localization and its connection to stability of sliding in the presence and absence of pore fluids. We also consider the possible influence of self-induced mechanical vibrations as well as the role of external acoustic vibrations as analogue for triggering by a distant event. For the dry case, our results suggest that at low and intermediate strain rates, persistent shear bands develop only in the absence of vibrations. Vibrations tend to fluidize the granular network and de-localize slip at these rates. Stick-slip is only observed for rough grains and it is confined to the shear band. At high strain rates, stick-slip disappears and the different systems exhibit similar stress-slip response. Changing the vibration intensity, duration or time of application alters the system response and may cause long-lasting rheological changes. The presence of pore fluids modifies the stick slip pattern and may lead to both loss and development of slip instability depending on the value of the confining pressure, imposed strain rate and hydraulic parameters. We analyze these observations in terms of possible transitions between rate strengthening and rate weakening response facilitated by a competition between shear induced dilation and acoustic compaction. We discuss the implications of our results on dynamic triggering, quiescence and strength evolution in gouge filled fault zones.

Relativistic fluid dynamics with spin

NASA Astrophysics Data System (ADS)

Florkowski, Wojciech; Friman, Bengt; Jaiswal, Amaresh; Speranza, Enrico

2018-04-01

Using the conservation laws for charge, energy, momentum, and angular momentum, we derive hydrodynamic equations for the charge density, local temperature, and fluid velocity, as well as for the polarization tensor, starting from local equilibrium distribution functions for particles and antiparticles with spin 1/2. The resulting set of differential equations extends the standard picture of perfect-fluid hydrodynamics with a conserved entropy current in a minimal way. This framework can be used in space-time analyses of the evolution of spin and polarization in various physical systems including high-energy nuclear collisions. We demonstrate that a stationary vortex, which exhibits vorticity-spin alignment, corresponds to a special solution of the spin-hydrodynamical equations.

Microscopic theory of topologically entangled fluids of rigid macromolecules

NASA Astrophysics Data System (ADS)

Sussman, Daniel M.; Schweizer, Kenneth S.

2011-06-01

We present a first-principles theory for the slow dynamics of a fluid of entangling rigid crosses of zero excluded volume based on a generalization of the dynamic mean-field approach of Szamel for infinitely thin nonrotating rods. The latter theory exactly includes topological constraints at the two-body collision level and self-consistently renormalizes an effective diffusion tensor to account for many-body effects. Remarkably, it predicts scaling laws consistent with the phenomenological reptation-tube predictions of Doi and Edwards for the long-time diffusion and the localization length in the heavily entangled limit. We generalize this approach to a different macromolecular architecture, infinitely thin three-dimensional crosses, and also extend the range of densities over which a dynamic localization length can be calculated for rods. Ideal gases of nonrotating crosses have recently received attention in computer simulations and are relevant as a simple model of both a strong-glass former and entangling star-branched polymers. Comparisons of our theory with these simulations reveal reasonable agreement for the magnitude and reduced density dependence of the localization length and also the self-diffusion constant if the consequences of local density fluctuations are taken into account.

Single DNA molecules on freestanding and supported cationic lipid bilayers: diverse conformational dynamics controlled by the local bilayer properties

NASA Astrophysics Data System (ADS)

Herold, Christoph; Schwille, Petra; Petrov, Eugene P.

2016-02-01

We present experimental results on the interaction of DNA macromolecules with cationic lipid membranes with different properties, including freestanding membranes in the fluid and gel state, and supported lipid membranes in the fluid state and under conditions of fluid-gel phase coexistence. We observe diverse conformational dynamics of membrane-bound DNA molecules controlled by the local properties of the lipid bilayer. In case of fluid-state freestanding lipid membranes, the behaviour of DNA on the membrane is controlled by the membrane charge density: whereas DNA bound to weakly charged membranes predominantly behaves as a 2D random coil, an increase in the membrane charge density leads to membrane-driven irreversible DNA collapse and formation of subresolution-sized DNA globules. On the other hand, electrostatic binding of DNA macromolecules to gel-state freestanding membranes leads to completely arrested diffusion and conformational dynamics of membrane-adsorbed DNA. A drastically different picture is observed in case of DNA interaction with supported cationic lipid bilayers: When the supported bilayer is in the fluid state, membrane-bound DNA molecules undergo 2D translational Brownian motion and conformational fluctuations, irrespectively of the charge density of the supported bilayer. At the same time, when the supported cationic membrane shows fluid-gel phase coexistence, membrane-bound DNA molecules are strongly attracted to micrometre-sized gel-phase domains enriched with the cationic lipid, which results in 2D compaction of the membrane-bound macromolecules. This DNA compaction, however, is fully reversible, and disappears as soon as the membrane is heated above the fluid-gel coexistence. We also discuss possible biological implications of our experimental findings.

Predicting Flow Reversals in a Computational Fluid Dynamics Simulated Thermosyphon Using Data Assimilation.

PubMed

Reagan, Andrew J; Dubief, Yves; Dodds, Peter Sheridan; Danforth, Christopher M

2016-01-01

A thermal convection loop is a annular chamber filled with water, heated on the bottom half and cooled on the top half. With sufficiently large forcing of heat, the direction of fluid flow in the loop oscillates chaotically, dynamics analogous to the Earth's weather. As is the case for state-of-the-art weather models, we only observe the statistics over a small region of state space, making prediction difficult. To overcome this challenge, data assimilation (DA) methods, and specifically ensemble methods, use the computational model itself to estimate the uncertainty of the model to optimally combine these observations into an initial condition for predicting the future state. Here, we build and verify four distinct DA methods, and then, we perform a twin model experiment with the computational fluid dynamics simulation of the loop using the Ensemble Transform Kalman Filter (ETKF) to assimilate observations and predict flow reversals. We show that using adaptively shaped localized covariance outperforms static localized covariance with the ETKF, and allows for the use of less observations in predicting flow reversals. We also show that a Dynamic Mode Decomposition (DMD) of the temperature and velocity fields recovers the low dimensional system underlying reversals, finding specific modes which together are predictive of reversal direction.

Predicting Flow Reversals in a Computational Fluid Dynamics Simulated Thermosyphon Using Data Assimilation

PubMed Central

Reagan, Andrew J.; Dubief, Yves; Dodds, Peter Sheridan; Danforth, Christopher M.

2016-01-01

A thermal convection loop is a annular chamber filled with water, heated on the bottom half and cooled on the top half. With sufficiently large forcing of heat, the direction of fluid flow in the loop oscillates chaotically, dynamics analogous to the Earth’s weather. As is the case for state-of-the-art weather models, we only observe the statistics over a small region of state space, making prediction difficult. To overcome this challenge, data assimilation (DA) methods, and specifically ensemble methods, use the computational model itself to estimate the uncertainty of the model to optimally combine these observations into an initial condition for predicting the future state. Here, we build and verify four distinct DA methods, and then, we perform a twin model experiment with the computational fluid dynamics simulation of the loop using the Ensemble Transform Kalman Filter (ETKF) to assimilate observations and predict flow reversals. We show that using adaptively shaped localized covariance outperforms static localized covariance with the ETKF, and allows for the use of less observations in predicting flow reversals. We also show that a Dynamic Mode Decomposition (DMD) of the temperature and velocity fields recovers the low dimensional system underlying reversals, finding specific modes which together are predictive of reversal direction. PMID:26849061

Enforcing the Courant-Friedrichs-Lewy condition in explicitly conservative local time stepping schemes

NASA Astrophysics Data System (ADS)

Gnedin, Nickolay Y.; Semenov, Vadim A.; Kravtsov, Andrey V.

2018-04-01

An optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubic "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a constraint on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.

Nanoscopic dynamics of phospholipid in unilamellar vesicles: Effect of gel to fluid phase transition

DOE PAGES

Sharma, V. K.; Mamontov, E.; Anunciado, D. B.; ...

2015-03-04

Dynamics of phospholipids in unilamellar vesicles (ULV) is of interest in biology, medical, and food sciences since these molecules are widely used as biocompatible agents and a mimic of cell membrane systems. We have investigated the nanoscopic dynamics of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) phospholipid in ULV as a function of temperature using elastic and quasielastic neutron scattering (QENS). The dependence of the signal on the scattering momentum transfer, which is a critical advantage of neutron scattering techniques, allows the detailed analysis of the lipid motions that cannot be carried out by other means. In agreement with a differential scanning calorimetry measurement, amore » sharp rise in the elastic scattering intensity below ca. 296 K indicates a phase transition from the high-temperature fluid phase to the low-temperature solid gel phase. The microscopic lipid dynamics exhibits qualitative differences between the solid gel phase (in a measurement at 280 K) and the fluid phase (in a measurement at a physiological temperature of 310 K). The data analysis invariably shows the presence of two distinct motions: the whole lipid molecule motion within a monolayer, or lateral diffusion, and the relatively faster internal motion of the DMPC molecule. The lateral diffusion of the whole lipid molecule is found to be Fickian in character, whereas the internal lipid motions are of localized character, consistent with the structure of the vesicles. The lateral motion slows down by an order of magnitude in the solid gel phase, whereas for the internal motion not only the time scale, but also the character of the motion changes upon the phase transition. In the solid gel phase, the lipids are more ordered and undergo uniaxial rotational motion. However, in the fluid phase, the hydrogen atoms of the lipid tails undergo confined translation diffusion rather than uniaxial rotational diffusion. The localized translational diffusion of the hydrogen atoms of the lipid tails is a manifestation of the flexibility of the chains acquired in the fluid phase. Because of this flexibility, both the local diffusivity and the confinement volume for the hydrogen atoms increase linearly from near the lipid s polar head group to the end of its hydrophobic tail. Our results present a quantitative and detailed picture of the effect of the gel-fluid phase transition on the nanoscopic lipid dynamics in ULV. Lastly, the data analysis approach developed here has a potential for probing the dynamic response of lipids to the presence of additional cell membrane components.« less

Light-cone reduction vs. TsT transformations: a fluid dynamics perspective

NASA Astrophysics Data System (ADS)

Dutta, Suvankar; Krishna, Hare

2018-05-01

We compute constitutive relations for a charged (2+1) dimensional Schrödinger fluid up to first order in derivative expansion, using holographic techniques. Starting with a locally boosted, asymptotically AdS, 4 + 1 dimensional charged black brane geometry, we uplift that to ten dimensions and perform TsT transformations to obtain an effective five dimensional local black brane solution with asymptotically Schrödinger isometries. By suitably implementing the holographic techniques, we compute the constitutive relations for the effective fluid living on the boundary of this space-time and extract first order transport coefficients from these relations. Schrödinger fluid can also be obtained by reducing a charged relativistic conformal fluid over light-cone. It turns out that both the approaches result the same system at the end. Fluid obtained by light-cone reduction satisfies a restricted class of thermodynamics. Here, we see that the charged fluid obtained holographically also belongs to the same restricted class.

RELATIONSHIPS BETWEEN FLUID VORTICITY, KINETIC HELICITY, AND MAGNETIC FIELD ON SMALL-SCALES (QUIET-NETWORK) ON THE SUN

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

Sangeetha, C. R.; Rajaguru, S. P., E-mail: crsangeetha@iiap.res.in

We derive horizontal fluid motions on the solar surface over large areas covering the quiet-Sun magnetic network from local correlation tracking of convective granules imaged in continuum intensity and Doppler velocity by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory . From these we calculate the horizontal divergence, the vertical component of vorticity, and the kinetic helicity of fluid motions. We study the correlations between fluid divergence and vorticity, and between vorticity (kinetic helicity) and the magnetic field. We find that the vorticity (kinetic helicity) around small-scale fields exhibits a hemispherical pattern (in sign) similar tomore » that followed by the magnetic helicity of large-scale active regions (containing sunspots). We identify this pattern to be a result of the Coriolis force acting on supergranular-scale flows (both the outflows and inflows), consistent with earlier studies using local helioseismology. Furthermore, we show that the magnetic fields cause transfer of vorticity from supergranular inflow regions to outflow regions, and that they tend to suppress the vortical motions around them when magnetic flux densities exceed about 300 G (from HMI). We also show that such an action of the magnetic fields leads to marked changes in the correlations between fluid divergence and vorticity. These results are speculated to be of importance to local dynamo action (if present) and to the dynamical evolution of magnetic helicity at the small-scale.« less

Local mass and energy transports in evaporation processes from a vapor-liquid interface in a slit pore based on molecular dynamics

NASA Astrophysics Data System (ADS)

Fujiwara, K.; Shibahara, M.

2018-02-01

Molecular evaporation processes from a vapor-liquid interface formed in a slit-like pore were examined based on the classical molecular dynamics method, in order to elucidate a molecular mechanism of local mass and energy transports in a slit. The calculation system consisted of monatomic molecules and atoms which interact through the 12-6 Lennard-Jones potential. At first, a liquid was situated in a slit with a vapor-liquid interface, and instantaneous amounts of the mass and energy fluxes defined locally in the slit were obtained in two dimensions to reveal local fluctuation properties of the fluid in equilibrium states. Then, imposing a temperature gradient in the calculation system, non-equilibrium evaporation processes in the slit were investigated in details based on the local mass and energy fluxes. In this study, we focused on the fluid which is in the vicinity of the solid surface and in contact with the vapor phase. In the non-equilibrium evaporation processes, the results revealed that the local energy transport mechanism in the vicinity of the solid surface is different from that of the vapor phase, especially in the case of the relatively strong fluid-solid interaction. The results also revealed that the local mass transport in the vicinity of the solid surface can be interpreted based on the mechanism of the local energy transport, and the mechanism provides valuable information about pictures of the evaporation phenomena especially in the vicinity of the hydrophilic surfaces. It suggests that evaluating and changing this mechanism of the local energy transport are necessary to control the local mass flux more precisely in the vicinity of the solid surface.

Dynamic nuclear polarization enhanced nuclear magnetic resonance and electron spin resonance studies of hydration and local water dynamics in micelle and vesicle assemblies.

PubMed

McCarney, Evan R; Armstrong, Brandon D; Kausik, Ravinath; Han, Songi

2008-09-16

We present a unique analysis tool for the selective detection of local water inside soft molecular assemblies (hydrophobic cores, vesicular bilayers, and micellar structures) suspended in bulk water. Through the use of dynamic nuclear polarization (DNP), the (1)H NMR signal of water is amplified, as it interacts with stable radicals that possess approximately 658 times higher spin polarization. We utilized stable nitroxide radicals covalently attached along the hydrophobic tail of stearic acid molecules that incorporate themselves into surfactant-based micelle or vesicle structures. Here, we present a study of local water content and fluid viscosity inside oleate micelles and vesicles and Triton X-100 micelles to serve as model systems for soft molecular assemblies. This approach is unique because the amplification of the NMR signal is performed in bulk solution and under ambient conditions with site-specific spin labels that only detect the water that is directly interacting with the localized spin labels. Continuous wave (cw) electron spin resonance (ESR) analysis provides rotational dynamics of the spin-labeled molecular chain segments and local polarity parameters that can be related to hydration properties, whereas we show that DNP-enhanced (1)H NMR analysis of fluid samples directly provides translational water dynamics and permeability of the local environment probed by the spin label. Our technique therefore has the potential to become a powerful analysis tool, complementary to cw ESR, to study hydration characteristics of surfactant assemblies, lipid bilayers, or protein aggregates, where water dynamics is a key parameter of their structure and function. In this study, we find that there is significant penetration of water inside the oleate micelles with a higher average local water viscosity (approximately 1.8 cP) than in bulk water, and Triton X-100 micelles and oleate vesicle bilayers mostly exclude water while allowing for considerable surfactant chain motion and measurable water permeation through the soft structure.

Locomotion at Low Reynolds Number: Dynamics in Newtonian and Non-Newtonian Systems with Biomedical Applications

NASA Astrophysics Data System (ADS)

Gagnon, David A.

Swimming microorganisms such as bacteria, spermatozoa, algae, and nematodes are critical to ubiquitous biological phenomena such as disease and infection, ecosystem dynamics, and mammalian fertilization. While there has been much scientific and practical interest in studying these swimmers in Newtonian (water-like) fluids, there are fewer systematic experimental studies on swimming through non-Newtonian (non-water-like) fluids with biologically-relevant mechanical properties. These organisms commonly swim through viscoelastic, structured, or shear-rate-dependent fluids, such as blood, mucus, and living tissues. Furthermore, the small length scales of these organisms dictate that their motion is dominated by viscous forces and inertia is negligible. Using rheology, microscopy, particle tracking, and image processing techniques, we examine the interaction of low Reynolds number swimmers and non-Newtonian fluids including viscoelastic, locally-anisotropic, and shear-thinning fluids. We then apply our understanding of locomotion to the study of the genetic disease Spinal Muscular Atrophy.

LES of a Jet Excited by the Localized Arc Filament Plasma Actuators

NASA Technical Reports Server (NTRS)

Brown, Clifford A.

2011-01-01

The fluid dynamics of a high-speed jet are governed by the instability waves that form in the free-shear boundary layer of the jet. Jet excitation manipulates the growth and saturation of particular instability waves to control the unsteady flow structures that characterize the energy cascade in the jet.The results may include jet noise mitigation or a reduction in the infrared signature of the jet. The Localized Arc Filament Plasma Actuators (LAFPA) have demonstrated the ability to excite a high-speed jets in laboratory experiments. Extending and optimizing this excitation technology, however, is a complex process that will require many tests and trials. Computational simulations can play an important role in understanding and optimizing this actuator technology for real-world applications. Previous research has focused on developing a suitable actuator model and coupling it with the appropriate computational fluid dynamics (CFD) methods using two-dimensional spatial flow approximations. This work is now extended to three-dimensions (3-D) in space. The actuator model is adapted to a series of discrete actuators and a 3-D LES simulation of an excited jet is run. The results are used to study the fluid dynamics near the actuator and in the jet plume.

Simulation of windblown dust transport from a mine tailings impoundment using a computational fluid dynamics model

NASA Astrophysics Data System (ADS)

Stovern, Michael; Felix, Omar; Csavina, Janae; Rine, Kyle P.; Russell, MacKenzie R.; Jones, Robert M.; King, Matt; Betterton, Eric A.; Sáez, A. Eduardo

2014-09-01

Mining operations are potential sources of airborne particulate metal and metalloid contaminants through both direct smelter emissions and wind erosion of mine tailings. The warmer, drier conditions predicted for the Southwestern US by climate models may make contaminated atmospheric dust and aerosols increasingly important, due to potential deleterious effects on human health and ecology. Dust emissions and dispersion of dust and aerosol from the Iron King Mine tailings in Dewey-Humboldt, Arizona, a Superfund site, are currently being investigated through in situ field measurements and computational fluid dynamics modeling. These tailings are heavily contaminated with lead and arsenic. Using a computational fluid dynamics model, we model dust transport from the mine tailings to the surrounding region. The model includes gaseous plume dispersion to simulate the transport of the fine aerosols, while individual particle transport is used to track the trajectories of larger particles and to monitor their deposition locations. In order to improve the accuracy of the dust transport simulations, both regional topographical features and local weather patterns have been incorporated into the model simulations. Results show that local topography and wind velocity profiles are the major factors that control deposition.

Simulation of windblown dust transport from a mine tailings impoundment using a computational fluid dynamics model.

PubMed

Stovern, Michael; Felix, Omar; Csavina, Janae; Rine, Kyle P; Russell, MacKenzie R; Jones, Robert M; King, Matt; Betterton, Eric A; Sáez, A Eduardo

2014-09-01

Mining operations are potential sources of airborne particulate metal and metalloid contaminants through both direct smelter emissions and wind erosion of mine tailings. The warmer, drier conditions predicted for the Southwestern US by climate models may make contaminated atmospheric dust and aerosols increasingly important, due to potential deleterious effects on human health and ecology. Dust emissions and dispersion of dust and aerosol from the Iron King Mine tailings in Dewey-Humboldt, Arizona, a Superfund site, are currently being investigated through in situ field measurements and computational fluid dynamics modeling. These tailings are heavily contaminated with lead and arsenic. Using a computational fluid dynamics model, we model dust transport from the mine tailings to the surrounding region. The model includes gaseous plume dispersion to simulate the transport of the fine aerosols, while individual particle transport is used to track the trajectories of larger particles and to monitor their deposition locations. In order to improve the accuracy of the dust transport simulations, both regional topographical features and local weather patterns have been incorporated into the model simulations. Results show that local topography and wind velocity profiles are the major factors that control deposition.

Simulation of windblown dust transport from a mine tailings impoundment using a computational fluid dynamics model

PubMed Central

Stovern, Michael; Felix, Omar; Csavina, Janae; Rine, Kyle P.; Russell, MacKenzie R.; Jones, Robert M.; King, Matt; Betterton, Eric A.; Sáez, A. Eduardo

2014-01-01

Mining operations are potential sources of airborne particulate metal and metalloid contaminants through both direct smelter emissions and wind erosion of mine tailings. The warmer, drier conditions predicted for the Southwestern US by climate models may make contaminated atmospheric dust and aerosols increasingly important, due to potential deleterious effects on human health and ecology. Dust emissions and dispersion of dust and aerosol from the Iron King Mine tailings in Dewey-Humboldt, Arizona, a Superfund site, are currently being investigated through in situ field measurements and computational fluid dynamics modeling. These tailings are heavily contaminated with lead and arsenic. Using a computational fluid dynamics model, we model dust transport from the mine tailings to the surrounding region. The model includes gaseous plume dispersion to simulate the transport of the fine aerosols, while individual particle transport is used to track the trajectories of larger particles and to monitor their deposition locations. In order to improve the accuracy of the dust transport simulations, both regional topographical features and local weather patterns have been incorporated into the model simulations. Results show that local topography and wind velocity profiles are the major factors that control deposition. PMID:25621085

Tectonic and hydrological controls on multiscale deformations in the Levant: numerical modeling and theoretical analysis

NASA Astrophysics Data System (ADS)

Belferman, Mariana; Katsman, Regina; Agnon, Amotz; Ben Avraham, Zvi

2016-04-01

Understanding the role of the dynamics of water bodies in triggering deformations in the upper crust and subsequently leading to earthquakes has been attracting considerable attention. We suggest that dynamic changes in the levels of the water bodies occupying tectonic depressions along the Dead Sea Transform (DST) cause significant variations in the shallow crustal stress field and affect local fault systems in a way that eventually leads to earthquakes. This mechanism and its spatial and temporal scales differ from those in tectonically-driven deformations. In this study we present a new thermo-mechanical model, constructed using the finite element method, and extended by including a fluid flow component in the upper crust. The latter is modeled on a basis of two-way poroelastic coupling with the momentum equation. This coupling is essential for capturing fluid flow evolution induced by dynamic water loading in the DST depressions and to resolve porosity changes. All the components of the model, namely elasticity, creep, plasticity, heat transfer, and fluid flow, have been extensively verified and presented in the study. The two-way coupling between localized plastic volumetric deformations and enhanced fluid flow is addressed, as well as the role of variability of the rheological and the hydrological parameters in inducing deformations in specific faulting environments. Correlations with historical and contemporary earthquakes in the region are discussed.

The effect of fluids on the frictional behavior of calcite gouge

NASA Astrophysics Data System (ADS)

Rempe, M.; Di Toro, G.; Mitchell, T. M.; Hirose, T.; Smith, S. A. F.; Renner, J.

2016-12-01

The presence of fluids in fault zones affects the faults' strength and the nucleation and propagation of earthquakes due to mechanical or physico-chemical weakening effects. To better understand the effect of pore fluids on the frictional behavior of gouge-bearing faults, a series of intermediate- to high-velocity experiments was conducted using the Phv rotary-shear apparatus (Kochi Core Center, Japan) equipped with a servo-controlled pore-fluid pressure system. Calcite gouge was sheared up to several meters displacement at room-humidity (dry) and water-saturated conditions. The pore-fluid factor, λ=pf/σn, ranged from 0.15 to 0.7 and the effective normal stress, σn,eff=σn-pf, from 1 to 12 MPa. Sheared samples were analyzed using scanning electron microscopy and Raman spectroscopy. The steady-state shear stress is lower for saturated than for dry gouges sliding at V=1 mm/s, possibly due to higher intergranular lubrication and/or accelerated subcritical crack growth, as evidenced also by the observed higher degree of compaction. At V=1 m/s, dry gouges show a pronounced strengthening phase preceding the onset of dynamic weakening; saturated gouges weaken abruptly. The higher λ, the lower the peak and steady-state shear stress, but -counterintuitively- the less localized deformation. Degree of weakening and localization might be influenced by insufficient drainage at high λ. In undrained experiments, the shear stress is slightly decreased likely due to thermal pressurization of the pore fluid, but the onset of dynamic weakening is not accelerated, indicating that dynamic weakening is due to more efficient mechanisms. For example, amorphous carbon may lubricate the slip surfaces of dry and saturated calcite gouges and cause dynamic weakening, yet Raman spectra only show the presence of disordered carbon on the principal slip surface. Furthermore, the presence of small recrystallized grains suggests that strain accommodation during steady-state slip might occur by non-frictional processes, such as grain-boundary sliding aided by diffusion creep.

Fluid-Solid Interaction and Multiscale Dynamic Processes: Experimental Approach

NASA Astrophysics Data System (ADS)

Arciniega-Ceballos, Alejandra; Spina, Laura; Mendo-Pérez, Gerardo M.; Guzmán-Vázquez, Enrique; Scheu, Bettina; Sánchez-Sesma, Francisco J.; Dingwell, Donald B.

2017-04-01

The speed and the style of a pressure drop in fluid-filled conduits determines the dynamics of multiscale processes and the elastic interaction between the fluid and the confining solid. To observe this dynamics we performed experiments using fluid-filled transparent tubes (15-50 cm long, 2-4 cm diameter and 0.3-1 cm thickness) instrumented with high-dynamic piezoelectric sensors and filmed the evolution of these processes with a high speed camera. We analyzed the response of Newtonian fluids to slow and sudden pressure drops from 3 bar-10 MPa to ambient pressure. We used fluids with viscosities of mafic to intermediate silicate melts of 1 to 1000 Pa s and water. The processes observed are fluid mass expansion, fluid flow, jets, bubbles nucleation, growth, coalescence and collapse, degassing, foam building at the surface and vertical wagging. All these processes (in fine and coarse scales) are triggered by the pressure drop and are sequentially coupled in time while interacting with the solid. During slow decompression, the multiscale processes are recognized occurring within specific pressure intervals, and exhibit a localized distribution along the conduit. In this, degassing predominates near the surface and may present piston-like oscillations. In contrast, during sudden decompression the fluid-flow reaches higher velocities, the dynamics is dominated by a sequence of gas-packet pulses driving jets of the gas-fluid mixture. The evolution of this multiscale phenomenon generates complex non-stationary microseismic signals recorded along the conduit. We discuss distinctive characteristics of these signals depending on the decompression style and compare them with synthetics. These synthetics are obtained numerically under an averaging modeling scheme, that accounted for the stress-strain of the cyclic dynamic interaction between the fluid and the solid wall, assuming an incompressible and viscous fluid that flows while the elastic solid responds oscillating. Analysis of time series, both experimental and synthetics, synchronized with high-speed imaging enables the explanation and interpretation of distinct phases of the dynamics of these fluids and the extraction of time and frequency characteristics of the individual processes. We observed that the effects of both, pressure drop triggering function and viscosity, control the characteristics of the micro-signals in time and frequency. This suggests the great potential that experimental and numerical approaches provide to untangle from field volcanic seismograms the multiscale processes of the stress field, driving forces and fluid-rock interaction that determine the volcanic conduit dynamics.

Multiscale Behavior of Viscous Fluids Dynamics: Experimental Observations

NASA Astrophysics Data System (ADS)

Arciniega-Ceballos, Alejandra; Spina, Laura; Scheu, Bettina; Dingwell, Donald B.

2016-04-01

The dynamics of Newtonian fluids with viscosities of mafic to intermediate silicate melts (10-1000 Pa s) during slow decompression present multi-time scale processes. To observe these processes we have performed several experiments on silicon oil saturated with Argon gas for 72 hours, in a Plexiglas autoclave. The slow decompression, dropping from 10 MPa to ambient pressure, acting as the excitation mechanism, triggered several processes with their own distinct timescales. These processes generate complex non-stationary microseismic signals, which have been recorded with 7 high-dynamic piezoelectric sensors located along the conduit flanked by high-speed video recordings. The analysis in time and frequency of these time series and their correlation with the associated high-speed imaging enables the characterization of distinct phases and the extraction of the individual processes during the evolution of decompression of these viscous fluids. We have observed fluid-solid elastic interaction, degassing, fluid mass expansion and flow, bubble nucleation, growth, coalescence and collapse, foam building and vertical wagging. All these processes (in fine and coarse scales) are sequentially coupled in time, occur within specific pressure intervals, and exhibit a localized distribution along the conduit. Their coexistence and interactions constitute the stress field and driving forces that determine the dynamics of the conduit system. Our observations point to the great potential of this experimental approach in the understanding of volcanic conduit dynamics and volcanic seismicity.

Integration of Extended MHD and Kinetic Effects in Global Magnetosphere Models

NASA Astrophysics Data System (ADS)

Germaschewski, K.; Wang, L.; Maynard, K. R. M.; Raeder, J.; Bhattacharjee, A.

2015-12-01

Computational models of Earth's geospace environment are an important tool to investigate the science of the coupled solar-wind -- magnetosphere -- ionosphere system, complementing satellite and ground observations with a global perspective. They are also crucial in understanding and predicting space weather, in particular under extreme conditions. Traditionally, global models have employed the one-fluid MHD approximation, which captures large-scale dynamics quite well. However, in Earth's nearly collisionless plasma environment it breaks down on small scales, where ion and electron dynamics and kinetic effects become important, and greatly change the reconnection dynamics. A number of approaches have recently been taken to advance global modeling, e.g., including multiple ion species, adding Hall physics in a Generalized Ohm's Law, embedding local PIC simulations into a larger fluid domain and also some work on simulating the entire system with hybrid or fully kinetic models, the latter however being to computationally expensive to be run at realistic parameters. We will present an alternate approach, ie., a multi-fluid moment model that is derived rigorously from the Vlasov-Maxwell system. The advantage is that the computational cost remains managable, as we are still solving fluid equations. While the evolution equation for each moment is exact, it depends on the next higher-order moment, so that truncating the hiearchy and closing the system to capture the essential kinetic physics is crucial. We implement 5-moment (density, momentum, scalar pressure) and 10-moment (includes pressure tensor) versions of the model, and use local approximations for the heat flux to close the system. We test these closures by local simulations where we can compare directly to PIC / hybrid codes, and employ them in global simulations using the next-generation OpenGGCM to contrast them to MHD / Hall-MHD results and compare with observations.

u-w formulation for dynamic problems in large deformation regime solved through an implicit meshfree scheme

NASA Astrophysics Data System (ADS)

Navas, Pedro; Sanavia, Lorenzo; López-Querol, Susana; Yu, Rena C.

2017-12-01

Solving dynamic problems for fluid saturated porous media at large deformation regime is an interesting but complex issue. An implicit time integration scheme is herein developed within the framework of the u-w (solid displacement-relative fluid displacement) formulation for the Biot's equations. In particular, liquid water saturated porous media is considered and the linearization of the linear momentum equations taking into account all the inertia terms for both solid and fluid phases is for the first time presented. The spatial discretization is carried out through a meshfree method, in which the shape functions are based on the principle of local maximum entropy LME. The current methodology is firstly validated with the dynamic consolidation of a soil column and the plastic shear band formulation of a square domain loaded by a rigid footing. The feasibility of this new numerical approach for solving large deformation dynamic problems is finally demonstrated through the application to an embankment problem subjected to an earthquake.

Amplification without instability: applying fluid dynamical insights in chemistry and biology

NASA Astrophysics Data System (ADS)

McCoy, Jonathan H.

2013-11-01

While amplification of small perturbations often arises from instability, transient amplification is possible locally even in asymptotically stable systems. That is, knowledge of a system's stability properties can mislead one's intuition for its transient behaviors. This insight, which has an interesting history in fluid dynamics, has more recently been rediscovered in ecology. Surprisingly, many nonlinear fluid dynamical and ecological systems share linear features associated with transient amplification of noise. This paper aims to establish that these features are widespread in many other disciplines concerned with noisy systems, especially chemistry, cell biology and molecular biology. Here, using classic nonlinear systems and the graphical language of network science, we explore how the noise amplification problem can be reframed in terms of activatory and inhibitory interactions between dynamical variables. The interaction patterns considered here are found in a great variety of systems, ranging from autocatalytic reactions and activator-inhibitor systems to influential models of nerve conduction, glycolysis, cell signaling and circadian rhythms.

Coarse-grained forms for equations describing the microscopic motion of particles in a fluid.

PubMed

Das, Shankar P; Yoshimori, Akira

2013-10-01

Exact equations of motion for the microscopically defined collective density ρ(x,t) and the momentum density ĝ(x,t) of a fluid have been obtained in the past starting from the corresponding Langevin equations representing the dynamics of the fluid particles. In the present work we average these exact equations of microscopic dynamics over the local equilibrium distribution to obtain stochastic partial differential equations for the coarse-grained densities with smooth spatial and temporal dependence. In particular, we consider Dean's exact balance equation for the microscopic density of a system of interacting Brownian particles to obtain the basic equation of the dynamic density functional theory with noise. Our analysis demonstrates that on thermal averaging the dependence of the exact equations on the bare interaction potential is converted to dependence on the corresponding thermodynamic direct correlation functions in the coarse-grained equations.

A generalized transport-velocity formulation for smoothed particle hydrodynamics

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

Zhang, Chi; Hu, Xiangyu Y., E-mail: xiangyu.hu@tum.de; Adams, Nikolaus A.

The standard smoothed particle hydrodynamics (SPH) method suffers from tensile instability. In fluid-dynamics simulations this instability leads to particle clumping and void regions when negative pressure occurs. In solid-dynamics simulations, it results in unphysical structure fragmentation. In this work the transport-velocity formulation of Adami et al. (2013) is generalized for providing a solution of this long-standing problem. Other than imposing a global background pressure, a variable background pressure is used to modify the particle transport velocity and eliminate the tensile instability completely. Furthermore, such a modification is localized by defining a shortened smoothing length. The generalized formulation is suitable formore » fluid and solid materials with and without free surfaces. The results of extensive numerical tests on both fluid and solid dynamics problems indicate that the new method provides a unified approach for multi-physics SPH simulations.« less

Analysis of the autonomous problem about coupled active non-Newtonian multi-seepage in sparse medium

NASA Astrophysics Data System (ADS)

Deng, Shuxian; Li, Hongen

2017-10-01

The flow field of non-Newtonian fluid in sparse medium was analyzed by computational fluid dynamics (CFD) method. The results show that the axial velocity and radial velocity of the non-Newtonian fluid are larger than those of the Newtonian fluid due to the coupling of the viscosity of the non-Newtonian fluid and the shear rate, and the tangential velocity is less than that of the Newtonian fluid. These differences lead to the difference in the sparse medium Non-Newtonian fluids are of a special nature. The influence of the weight function on the global existence and blasting of the problem is discussed by analyzing the non-Newtonian percolation equation with nonlocal and weighted non-local Dirichlet boundary conditions. According to the non-Newtonian percolation equation, we define the weak solution of the problem and expound the local existence of the weak solution. Then we construct the test function and prove the weak comparison principle by using the Grown well inequality. The overall existence and blasting are analyzed by constructing the upper and lower solutions.

Enforcing the Courant–Friedrichs–Lewy condition in explicitly conservative local time stepping schemes

DOE PAGES

Gnedin, Nickolay Y.; Semenov, Vadim A.; Kravtsov, Andrey V.

2018-01-30

In this study, an optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubicmore » "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a condition on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.« less

Droplets and the three-phase contact line at the nano-scale. Statics and dynamics

NASA Astrophysics Data System (ADS)

Yatsyshin, Petr; Sibley, David; Savva, Nikos; Kalliadasis, Serafim

2014-11-01

Understanding the behaviour of the solid-liquid-vapour contact line at the scale of several tens of molecular diameters is important in wetting hydrodynamics with applications in micro- and nano-fluidics, including the design of lab-on-a-chip devices and surfaces with specific wetting properties. Due to the fluid inhomogeneity at the nano-scale, the application of continuum-mechanical approaches is limited, and a natural way to remedy this is to seek descriptions accounting for the non-local molecular-level interactions. Density Functional Theory (DFT) for fluids offers a statistical-mechanical framework based on expressing the free energy of the fluid-solid pair as a functional of the spatially varying fluid density. DFT allows us to investigate small drops deposited on planar substrates whilst keeping track of the microscopic structural details of the fluid. Starting from a model of intermolecular forces, we systematically obtain interfaces, surface tensions, and the microscopic contact angle. Using a dynamic extension of equilibrium DFT, we investigate the diffusion-driven evolution of the three-phase contact line to gain insight into the dynamic behaviour of the microscopic contact angle, which is still under debate.

Alfvén ionization in an MHD-gas interactions code

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

Wilson, A. D.; Diver, D. A.

A numerical model of partially ionized plasmas is developed in order to capture their evolving ionization fractions as a result of Alfvén ionization (AI). The mechanism of, and the parameter regime necessary for, AI is discussed and an expression for the AI rate based on fluid parameters, from a gas-MHD model, is derived. This AI term is added to an existing MHD-gas interactions' code, and the result is a linear, 2D, two-fluid model that includes momentum transfer between charged and neutral species as well as an ionization rate that depends on the velocity fields of both fluids. The dynamics ofmore » waves propagating through such a partially ionized plasma are investigated, and it is found that AI has a significant influence on the fluid dynamics as well as both the local and global ionization fraction.« less

Detecting Rotational Superradiance in Fluid Laboratories

NASA Astrophysics Data System (ADS)

Cardoso, Vitor; Coutant, Antonin; Richartz, Mauricio; Weinfurtner, Silke

2016-12-01

Rotational superradiance was predicted theoretically decades ago, and is chiefly responsible for a number of important effects and phenomenology in black-hole physics. However, rotational superradiance has never been observed experimentally. Here, with the aim of probing superradiance in the lab, we investigate the behavior of sound and surface waves in fluids resting in a circular basin at the center of which a rotating cylinder is placed. We show that with a suitable choice for the material of the cylinder, surface and sound waves are amplified. Two types of instabilities are studied: one sets in whenever superradiant modes are confined near the rotating cylinder and the other, which does not rely on confinement, corresponds to a local excitation of the cylinder. Our findings are experimentally testable in existing fluid laboratories and, hence, offer experimental exploration and comparison of dynamical instabilities arising from rapidly rotating boundary layers in astrophysical as well as in fluid dynamical systems.

Effective field theory of dissipative fluids

DOE PAGES

Crossley, Michael; Glorioso, Paolo; Liu, Hong

2017-09-20

We develop an effctive fi eld theory for dissipative fluids which governs the dynamics of long-lived gapless modes associated with conserved quantities. The resulting theory gives a path integral formulation of fluctuating hydrodynamics which systematically incorporates nonlinear interactions of noises. The dynamical variables are mappings between a "fluid spacetime" and the physical spacetime and an essential aspect of our formulation is to identify the appropriate symmetries in the fluid spacetime. The theory applies to nonlinear disturbances around a general density matrix. For a thermal density matrix, we require an additional Z2 symmetry, to which we refer as the local KMSmore » condition. This leads to the standard constraints of hydrodynamics, as well as a nonlinear generalization of the Onsager relations. It also leads to an emergent supersymmetry in the classical statistical regime, and a higher derivative deformation of supersymmetry in the full quantum regime.« less

Effective field theory of dissipative fluids

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

Crossley, Michael; Glorioso, Paolo; Liu, Hong

We develop an effctive fi eld theory for dissipative fluids which governs the dynamics of long-lived gapless modes associated with conserved quantities. The resulting theory gives a path integral formulation of fluctuating hydrodynamics which systematically incorporates nonlinear interactions of noises. The dynamical variables are mappings between a "fluid spacetime" and the physical spacetime and an essential aspect of our formulation is to identify the appropriate symmetries in the fluid spacetime. The theory applies to nonlinear disturbances around a general density matrix. For a thermal density matrix, we require an additional Z2 symmetry, to which we refer as the local KMSmore » condition. This leads to the standard constraints of hydrodynamics, as well as a nonlinear generalization of the Onsager relations. It also leads to an emergent supersymmetry in the classical statistical regime, and a higher derivative deformation of supersymmetry in the full quantum regime.« less

Active learning of constitutive relation from mesoscopic dynamics for macroscopic modeling of non-Newtonian flows

NASA Astrophysics Data System (ADS)

Zhao, Lifei; Li, Zhen; Caswell, Bruce; Ouyang, Jie; Karniadakis, George Em

2018-06-01

We simulate complex fluids by means of an on-the-fly coupling of the bulk rheology to the underlying microstructure dynamics. In particular, a continuum model of polymeric fluids is constructed without a pre-specified constitutive relation, but instead it is actively learned from mesoscopic simulations where the dynamics of polymer chains is explicitly computed. To couple the bulk rheology of polymeric fluids and the microscale dynamics of polymer chains, the continuum approach (based on the finite volume method) provides the transient flow field as inputs for the (mesoscopic) dissipative particle dynamics (DPD), and in turn DPD returns an effective constitutive relation to close the continuum equations. In this multiscale modeling procedure, we employ an active learning strategy based on Gaussian process regression (GPR) to minimize the number of expensive DPD simulations, where adaptively selected DPD simulations are performed only as necessary. Numerical experiments are carried out for flow past a circular cylinder of a non-Newtonian fluid, modeled at the mesoscopic level by bead-spring chains. The results show that only five DPD simulations are required to achieve an effective closure of the continuum equations at Reynolds number Re = 10. Furthermore, when Re is increased to 100, only one additional DPD simulation is required for constructing an extended GPR-informed model closure. Compared to traditional message-passing multiscale approaches, applying an active learning scheme to multiscale modeling of non-Newtonian fluids can significantly increase the computational efficiency. Although the method demonstrated here obtains only a local viscosity from the polymer dynamics, it can be extended to other multiscale models of complex fluids whose macro-rheology is unknown.

Computational Fluid Dynamics of Developing Avian Outflow Tract Heart Valves

PubMed Central

Bharadwaj, Koonal N.; Spitz, Cassie; Shekhar, Akshay; Yalcin, Huseyin C.; Butcher, Jonathan T.

2012-01-01

Hemodynamic forces play an important role in sculpting the embryonic heart and its valves. Alteration of blood flow patterns through the hearts of embryonic animal models lead to malformations that resemble some clinical congenital heart defects, but the precise mechanisms are poorly understood. Quantitative understanding of the local fluid forces acting in the heart has been elusive because of the extremely small and rapidly changing anatomy. In this study, we combine multiple imaging modalities with computational simulation to rigorously quantify the hemodynamic environment within the developing outflow tract (OFT) and its eventual aortic and pulmonary valves. In vivo Doppler ultrasound generated velocity profiles were applied to Micro-Computed Tomography generated 3D OFT lumen geometries from Hamburger-Hamilton (HH) stage 16 to 30 chick embryos. Computational fluid dynamics simulation initial conditions were iterated until local flow profiles converged with in vivo Doppler flow measurements. Results suggested that flow in the early tubular OFT (HH16 and HH23) was best approximated by Poiseuille flow, while later embryonic OFT septation (HH27, HH30) was mimicked by plug flow conditions. Peak wall shear stress (WSS) values increased from 18.16 dynes/cm2 at HH16 to 671.24 dynes/cm2 at HH30. Spatiotemporally averaged WSS values also showed a monotonic increase from 3.03 dynes/cm2 at HH16 to 136.50 dynes/cm2 at HH30. Simulated velocity streamlines in the early heart suggest a lack of mixing, which differed from classical ink injections. Changes in local flow patterns preceded and correlated with key morphogenetic events such as OFT septation and valve formation. This novel method to quantify local dynamic hemodynamics parameters affords insight into sculpting role of blood flow in the embryonic heart and provides a quantitative baseline dataset for future research. PMID:22535311

Reference Computational Meshing Strategy for Computational Fluid Dynamics Simulation of Departure from Nucleate BoilingReference Computational Meshing Strategy for Computational Fluid Dynamics Simulation of Departure from Nucleate Boiling

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

Pointer, William David

The objective of this effort is to establish a strategy and process for generation of suitable computational mesh for computational fluid dynamics simulations of departure from nucleate boiling in a 5 by 5 fuel rod assembly held in place by PWR mixing vane spacer grids. This mesh generation process will support ongoing efforts to develop, demonstrate and validate advanced multi-phase computational fluid dynamics methods that enable more robust identification of dryout conditions and DNB occurrence.Building upon prior efforts and experience, multiple computational meshes were developed using the native mesh generation capabilities of the commercial CFD code STAR-CCM+. These meshes weremore » used to simulate two test cases from the Westinghouse 5 by 5 rod bundle facility. The sensitivity of predicted quantities of interest to the mesh resolution was then established using two evaluation methods, the Grid Convergence Index method and the Least Squares method. This evaluation suggests that the Least Squares method can reliably establish the uncertainty associated with local parameters such as vector velocity components at a point in the domain or surface averaged quantities such as outlet velocity magnitude. However, neither method is suitable for characterization of uncertainty in global extrema such as peak fuel surface temperature, primarily because such parameters are not necessarily associated with a fixed point in space. This shortcoming is significant because the current generation algorithm for identification of DNB event conditions relies on identification of such global extrema. Ongoing efforts to identify DNB based on local surface conditions will address this challenge« less

Geophysical Fluid Dynamics Laboratory Open Days at the Woods Hole Oceanographic Institution

NASA Astrophysics Data System (ADS)

Hyatt, Jason; Cenedese, Claudia; Jensen, Anders

2015-11-01

This event was hosted for one week for two consecutive years in 2013 and 2014. It targeted postdocs, graduate students, K-12 students and local community participation. The Geophysical Fluid Dynamics Laboratory at the Woods Hole Oceanographic Institution hosted 10 hands-on demonstrations and displays, with something for all ages, to share the excitement of fluid mechanics and oceanography. The demonstrations/experiments spanned as many fluid mechanics problems as possible in all fields of oceanography and gave insight into using fluids laboratory experiments as a research tool. The chosen experiments were `simple' yet exciting for a 6 year old child, a high school student, a graduate student, and a postdoctoral fellow from different disciplines within oceanography. The laboratory is a perfect environment in which to create excitement and stimulate curiosity. Even what we consider `simple' experiments can fascinate and generate interesting questions from both a 6 year old child and a physics professor. How does an avalanche happen? How does a bath tub vortex form? What happens to waves when they break? How does a hurricane move? Hands-on activities in the fluid dynamics laboratory helped students of all ages in answering these and other intriguing questions. The laboratory experiments/demonstrations were accompanied by `live' videos to assist in the interpretation of the demonstrations. Posters illustrated the oceanographic/scientific applicability and the location on Earth where the dynamics in the experiments occur. Support was given by the WHOI Doherty Chair in Education.

Spontaneous flow in polar active fluids: the effect of a phenomenological self propulsion-like term.

PubMed

Bonelli, Francesco; Gonnella, Giuseppe; Tiribocchi, Adriano; Marenduzzo, Davide

2016-01-01

We present hybrid lattice Boltzmann simulations of extensile and contractile active fluids where we incorporate phenomenologically the tendency of active particles such as cell and bacteria, to move, or swim, along the local orientation. Quite surprisingly, we show that the interplay between alignment and activity can lead to completely different results, according to geometry (periodic boundary conditions or confinement between flat walls) and nature of the activity (extensile or contractile). An interesting generic outcome is that the alignment interaction can transform stationary active patterns into continuously moving ones: the dynamics of these evolving patterns can be oscillatory or chaotic according to the strength of the alignment term. Our results suggest that flow-polarisation alignment can have important consequences on the collective dynamics of active fluids and active gel.

New views of granular mass flows

USGS Publications Warehouse

Iverson, R.M.; Vallance, J.W.

2001-01-01

Concentrated grain-fluid mixtures in rock avalanches, debris flows, and pyroclastic flows do not behave as simple materials with fixed rheologies. Instead, rheology evolves as mixture agitation, grain concentration, and fluid-pressure change during flow initiation, transit, and deposition. Throughout a flow, however, normal forces on planes parallel to the free upper surface approximately balance the weight of the superincumbent mixture, and the Coulomb friction rule describes bulk intergranular shear stresses on such planes. Pore-fluid pressure can temporarily or locally enhance mixture mobility by reducing Coulomb friction and transferring shear stress to the fluid phase. Initial conditions, boundary conditions, and grain comminution and sorting can influence pore-fluid pressures and cause variations in flow dynamics and deposits.

Response to “Comment on ‘General rotating quantum vortex filaments in the low-temperature Svistunov model of the local induction approximation’” [Phys. Fluids 26, 119101 (2014)

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

Van Gorder, Robert A., E-mail: rav@knights.ucf.edu

2014-11-15

In R. A. Van Gorder, “General rotating quantum vortex filaments in the low-temperature Svistunov model of the local induction approximation,” Phys. Fluids 26, 065105 (2014) I discussed properties of generalized vortex filaments exhibiting purely rotational motion under the low-temperature Svistunov model of the local induction approximation. Such solutions are stationary in terms of translational motion. In the Comment [N. Hietala, “Comment on ‘General rotating quantum vortex filaments in the low-temperature Svistunov model of the local induction approximation’ [Phys. Fluids 26, 065105 (2014)],” Phys. Fluids 26, 119101 (2014)], the author criticizes my paper for not including translational motion (although it wasmore » clearly stated that the filament motion was assumed rotational). As it turns out, if one is interested in studying the geometric structure of solutions (which was the point of my paper), one obtains the needed qualitative results on the structure of such solutions by studying the purely rotational case. Nevertheless, in this Response I shall discuss the vortex filaments that have both rotational and translational motions. I then briefly discuss why one might want to study such generalized rotating filament solutions, in contrast to simple the standard helical or planar examples (which are really special cases). I also discuss how one can study the time evolution of filaments which exhibit more complicated dynamics than pure translation and rotation. Doing this, one can study non-stationary solutions which initially appear purely rotational and gradually display other dynamics as the filaments evolve.« less

Fluid Dynamics of Magnetic Nanoparticles in Simulated Blood Vessels

NASA Astrophysics Data System (ADS)

Blue, Lauren; Sewell, Mary Kathryn; Brazel, Christopher S.

2008-11-01

Magnetic nanoparticles (MNPs) can be used to locally target therapies and offer the benefit of using an AC magnetic field to combine hyperthermia treatment with the triggered release of therapeutic agents. Here, we investigate localization of MNPs in a simulated environment to understand the relationship between magnetic field intensity and bulk fluid dynamics to determine MNP retention in a simulated blood vessel. As MNPs travel through blood vessels, they can be slowed or trapped in a specific area by applying a magnetic field. Magnetic cobalt ferrite nanoparticles were synthesized and labeled with a fluorescent rhodamine tag to visualize patterns in a flow cell, as monitored by a fluorescence microscope. Particle retention was determined as a function of flow rate, concentration, and magnetic field strength. Understanding the relationship between magnetic field intensity, flow behavior and nanoparticle characteristics will aid in the development of therapeutic systems specifically targeted to diseased tissue.

A generalised porous medium approach to study thermo-fluid dynamics in human eyes.

PubMed

Mauro, Alessandro; Massarotti, Nicola; Salahudeen, Mohamed; Romano, Mario R; Romano, Vito; Nithiarasu, Perumal

2018-03-22

The present work describes the application of the generalised porous medium model to study heat and fluid flow in healthy and glaucomatous eyes of different subject specimens, considering the presence of ocular cavities and porous tissues. The 2D computational model, implemented into the open-source software OpenFOAM, has been verified against benchmark data for mixed convection in domains partially filled with a porous medium. The verified model has been employed to simulate the thermo-fluid dynamic phenomena occurring in the anterior section of four patient-specific human eyes, considering the presence of anterior chamber (AC), trabecular meshwork (TM), Schlemm's canal (SC), and collector channels (CC). The computational domains of the eye are extracted from tomographic images. The dependence of TM porosity and permeability on intraocular pressure (IOP) has been analysed in detail, and the differences between healthy and glaucomatous eye conditions have been highlighted, proving that the different physiological conditions of patients have a significant influence on the thermo-fluid dynamic phenomena. The influence of different eye positions (supine and standing) on thermo-fluid dynamic variables has been also investigated: results are presented in terms of velocity, pressure, temperature, friction coefficient and local Nusselt number. The results clearly indicate that porosity and permeability of TM are two important parameters that affect eye pressure distribution. Graphical abstract Velocity contours and vectors for healthy eyes (top) and glaucomatous eyes (bottom) for standing position.

Advances in Quantum Trajectory Approaches to Dynamics

NASA Astrophysics Data System (ADS)

Askar, Attila

2001-03-01

The quantum fluid dynamics (QFD) formulation is based on the separation of the amplitude and phase of the complex wave function in Schrodinger's equation. The approach leads to conservation laws for an equivalent "gas continuum". The Lagrangian [1] representation corresponds to following the particles of the fluid continuum, i. e. calculating "quantum trajectories". The Eulerian [2] representation on the other hand, amounts to observing the dynamics of the gas continuum at the points of a fixed coordinate frame. The combination of several factors leads to a most encouraging computational efficiency. QFD enables the numerical analysis to deal with near monotonic amplitude and phase functions. The Lagrangian description concentrates the computation effort to regions of highest probability as an optimal adaptive grid. The Eulerian representation allows the study of multi-coordinate problems as a set of one-dimensional problems within an alternating direction methodology. An explicit time integrator limits the increase in computational effort with the number of discrete points to linear. Discretization of the space via local finite elements [1,2] and global radial functions [3] will be discussed. Applications include wave packets in four-dimensional quadratic potentials and two coordinate photo-dissociation problems for NOCl and NO2. [1] "Quantum fluid dynamics (QFD) in the Lagrangian representation with applications to photo-dissociation problems", F. Sales, A. Askar and H. A. Rabitz, J. Chem. Phys. 11, 2423 (1999) [2] "Multidimensional wave-packet dynamics within the fluid dynamical formulation of the Schrodinger equation", B. Dey, A. Askar and H. A. Rabitz, J. Chem. Phys. 109, 8770 (1998) [3] "Solution of the quantum fluid dynamics equations with radial basis function interpolation", Xu-Guang Hu, Tak-San Ho, H. A. Rabitz and A. Askar, Phys. Rev. E. 61, 5967 (2000)

General connected and reconnected fields in plasmas

NASA Astrophysics Data System (ADS)

Mahajan, Swadesh M.; Asenjo, Felipe A.

2018-02-01

For plasma dynamics, more encompassing than the magnetohydrodynamical (MHD) approximation, the foundational concepts of "magnetic reconnection" may require deep revisions because, in the larger dynamics, magnetic field is no longer connected to the fluid lines; it is replaced by more general fields (one for each plasma specie) that are weighted combination of the electromagnetic and the thermal-vortical fields. We study the two-fluid plasma dynamics plasma expressed in two different sets of variables: the two-fluid (2F) description in terms of individual fluid velocities, and the one-fluid (1F) variables comprising the plasma bulk motion and plasma current. In the 2F description, a Connection Theorem is readily established; we show that, for each specie, there exists a Generalized (Magnetofluid/Electro-Vortic) field that is frozen-in the fluid and consequently remains, forever, connected to the flow. This field is an expression of the unification of the electromagnetic, and fluid forces (kinematic and thermal) for each specie. Since the magnetic field, by itself, is not connected in the first place, its reconnection is never forbidden and does not require any external agency (like resistivity). In fact, a magnetic field reconnection (local destruction) must be interpreted simply as a consequence of the preservation of the dynamical structure of the unified field. In the 1F plasma description, however, it is shown that there is no exact physically meaningful Connection Theorem; a general and exact field does not exist, which remains connected to the bulk plasma flow. It is also shown that the helicity conservation and the existence of a Connected field follow from the same dynamical structure; the dynamics must be expressible as an ideal Ohm's law with a physical velocity. This new perspective, emerging from the analysis of the post MHD physics, must force us to reexamine the meaning as well as our understanding of magnetic reconnection.

Unraveling earthquake stresses: Insights from dynamically triggered and induced earthquakes

NASA Astrophysics Data System (ADS)

Velasco, A. A.; Alfaro-Diaz, R. A.

2017-12-01

Induced seismicity, earthquakes caused by anthropogenic activity, has more than doubled in the last several years resulting from practices related to oil and gas production. Furthermore, large earthquakes have been shown to promote the triggering of other events within two fault lengths (static triggering), due to static stresses caused by physical movement along the fault, and also remotely from the passage of seismic waves (dynamic triggering). Thus, in order to understand the mechanisms for earthquake failure, we investigate regions where natural, induced, and dynamically triggered events occur, and specifically target Oklahoma. We first analyze data from EarthScope's USArray Transportable Array (TA) and local seismic networks implementing an optimized (STA/LTA) detector in order to develop local detection and earthquake catalogs. After we identify triggered events through statistical analysis, and perform a stress analysis to gain insight on the stress-states leading to triggered earthquake failure. We use our observations to determine the role of different transient stresses in contributing to natural and induced seismicity by comparing these stresses to regional stress orientation. We also delineate critically stressed regions of triggered seismicity that may indicate areas susceptible to earthquake hazards associated with sustained fluid injection in provinces of induced seismicity. Anthropogenic injection and extraction activity can alter the stress state and fluid flow within production basins. By analyzing the stress release of these ancient faults caused by dynamic stresses, we may be able to determine if fluids are solely responsible for increased seismic activity in induced regions.

Mount St. Helens (Washington, USA) and World Trade Center (New York, USA) collapse: a fluid dynamic analogy

NASA Astrophysics Data System (ADS)

Doronzo, Domenico; de Tullio, Marco; Pascazio, Giuseppe; Dellino, Pierfrancesco

2013-04-01

When a skyscraper collapses, the non-fragmented material is rapidly deposited close to the source, whereas the fragmented counterpart is loaded turbulently in the associated currents. Indeed, on impact with the ground, collapses of volcanic columns, domes, or sectors of volcanoes generate thick deposits of coarser material, and from there on the finer material is suspended over the landscape, to be re-deposited far away in thin deposits. Here, we explore the multiphase fluid dynamic behavior of the World Trade Center (New York, USA) collapse, which on 11 September 2001 followed the fragmentation of the Twin Towers, and generated shear dusty currents. These currents had a multiphase and turbulent behavior, and resemble the volcanic flow generated during the 18 May 1980 explosive eruption of Mount St. Helens (Washington, USA), in which a sector of the volcano collapsed, then a highly mobile, multiphase turbulent current followed and heavily interacted with the surrounding landscape. This analogy allows to focus on the comparison between volcanic and skyscraper collapse. A computational fluid dynamic investigation, along with a locally refined Cartesian grid, are adopted to simulate numerically the propagation of the 11 September dusty currents in Manhattan. Results of flow dynamic pressure, the parameter of volcanic hazard, and particle deposition reveal that the pressure can locally increase up to a factor 10 because of flow-building interaction. Also, the surrounding buildings make the urban setting as of a high turbulence and exponential decay of deposit thickness.

Analog geometry in an expanding fluid from AdS/CFT perspective

NASA Astrophysics Data System (ADS)

Bilić, Neven; Domazet, Silvije; Tolić, Dijana

2015-04-01

The dynamics of an expanding hadron fluid at temperatures below the chiral transition is studied in the framework of AdS/CFT correspondence. We establish a correspondence between the asymptotic AdS geometry in the 4 + 1 dimensional bulk with the analog spacetime geometry on its 3 + 1 dimensional boundary with the background fluid undergoing a spherical Bjorken type expansion. The analog metric tensor on the boundary depends locally on the soft pion dispersion relation and the four-velocity of the fluid. The AdS/CFT correspondence provides a relation between the pion velocity and the critical temperature of the chiral phase transition.

Computational Fluid Dynamics Simulations of Hemodynamics in Plaque Erosion

PubMed Central

Campbell, Ian C.; Timmins, Lucas H.; Giddens, Don P.; Virmani, Renu; Veneziani, Alessandro; Rab, S. Tanveer; Samady, Habib; McDaniel, Michael C.; Finn, Aloke V.; Taylor, W. Robert; Oshinski, John N.

2013-01-01

Purpose We investigated whether local hemodynamics were associated with sites of plaque erosion and hypothesized that patients with plaque erosion have locally elevated WSS magnitude in regions where erosion has occurred. Methods We generated 3D, patient-specific models of coronary arteries from biplane angiographic images in 3 human patients with plaque erosion diagnosed by optical coherence tomography (OCT). Using computational fluid dynamics, we simulated pulsatile blood flow and calculated both wall shear stress (WSS) and oscillatory shear index (OSI). We also investigated anatomic features of plaque erosion sites by examining branching and local curvature in x-ray angiograms of barium-perfused autopsy hearts. Results Neither high nor low magnitudes of mean WSS were associated with sites of plaque erosion. OSI and local curvature were also not associated with erosion. Anatomically, 8 of 13 hearts had a nearby bifurcation upstream of the site of plaque erosion. Conclusions This study provides preliminary evidence that neither hemodynamics nor anatomy are predictors of plaque erosion, based upon a very unique dataset. Our sample sizes are small, but this dataset suggests that high magnitudes of wall shear stress, one potential mechanism for inducing plaque erosion, are not necessary for erosion to occur. PMID:24223678

Fluid dynamic propagation of initial baryon number perturbations on a Bjorken flow background

DOE PAGES

Floerchinger, Stefan; Martinez, Mauricio

2015-12-11

Baryon number density perturbations offer a possible route to experimentally measure baryon number susceptibilities and heat conductivity of the quark gluon plasma. We study the fluid dynamical evolution of local and event-by-event fluctuations of baryon number density, flow velocity, and energy density on top of a (generalized) Bjorken expansion. To that end we use a background-fluctuation splitting and a Bessel-Fourier decomposition for the fluctuating part of the fluid dynamical fields with respect to the azimuthal angle, the radius in the transverse plane, and rapidity. Here, we examine how the time evolution of linear perturbations depends on the equation of statemore » as well as on shear viscosity, bulk viscosity, and heat conductivity for modes with different azimuthal, radial, and rapidity wave numbers. Finally we discuss how this information is accessible to experiments in terms of the transverse and rapidity dependence of correlation functions for baryonic particles in high energy nuclear collisions.« less

Cloud fluid models of gas dynamics and star formation in galaxies

NASA Technical Reports Server (NTRS)

Struck-Marcell, Curtis; Scalo, John M.; Appleton, P. N.

1987-01-01

The large dynamic range of star formation in galaxies, and the apparently complex environmental influences involved in triggering or suppressing star formation, challenges the understanding. The key to this understanding may be the detailed study of simple physical models for the dominant nonlinear interactions in interstellar cloud systems. One such model is described, a generalized Oort model cloud fluid, and two simple applications of it are explored. The first of these is the relaxation of an isolated volume of cloud fluid following a disturbance. Though very idealized, this closed box study suggests a physical mechanism for starbursts, which is based on the approximate commensurability of massive cloud lifetimes and cloud collisional growth times. The second application is to the modeling of colliding ring galaxies. In this case, the driving processes operating on a dynamical timescale interact with the local cloud processes operating on the above timescale. The results is a variety of interesting nonequilibrium behaviors, including spatial variations of star formation that do not depend monotonically on gas density.

Synchronization crossover of polariton condensates in weakly disordered lattices

NASA Astrophysics Data System (ADS)

Ohadi, H.; del Valle-Inclan Redondo, Y.; Ramsay, A. J.; Hatzopoulos, Z.; Liew, T. C. H.; Eastham, P. R.; Savvidis, P. G.; Baumberg, J. J.

2018-05-01

We demonstrate that the synchronization of a lattice of solid-state condensates when intersite tunneling is switched on depends strongly on the weak local disorder. This finding is vital for implementation of condensate arrays as computation devices. The condensates here are nonlinear bosonic fluids of exciton-polaritons trapped in a weakly disordered Bose-Hubbard potential, where the nearest-neighboring tunneling rate (Josephson coupling) can be dynamically tuned. The system can thus be tuned from a localized to a delocalized fluid as the number density or the Josephson coupling between nearest neighbors increases. The localized fluid is observed as a lattice of unsynchronized condensates emitting at different energies set by the disorder potential. In the delocalized phase, the condensates synchronize and long-range order appears, evidenced by narrowing of momentum and energy distributions, new diffraction peaks in momentum space, and spatial coherence between condensates. Our paper identifies similarities and differences of this nonequilibrium crossover to the traditional Bose-glass to superfluid transition in atomic condensates.

Thermally driven mass flows in the convection zone of the sun

NASA Technical Reports Server (NTRS)

Dijkhuis, G. C.

1973-01-01

A formulation of the fluid dynamics of convective regions is developed which leads to an analytical description of the solar rotation, the Evershed flow, and the supergranulation. The starting point of the present formulation is the mixing length picture of convective equilibrium, but the earlier point mass model for convective molecules is replaced here by a model with both inertia and intrinsic moment of inertia. This extension introduces three rotational degrees of freedom into the dynamics of individual convective molecules, which enter into the dynamical equations for a mixing length fluid in the form of a separate vector field which we term the spin field. It is shown that for convective molecules having a spherically symmetric mass distribution, the spin field is proportional to the local vorticity.

On locally and nonlocally related potential systems

NASA Astrophysics Data System (ADS)

Cheviakov, Alexei F.; Bluman, George W.

2010-07-01

For any partial differential equation (PDE) system, a local conservation law yields potential equations in terms of some potential variable, which normally is a nonlocal variable. The current paper examines situations when such a potential variable is a local variable, i.e., is a function of the independent and dependent variables of a given PDE system, and their derivatives. In the case of two independent variables, a simple necessary and sufficient condition is presented for the locality of such a potential variable, and this is illustrated by several examples. As a particular example, two-dimensional reductions of equilibrium equations for fluid and plasma dynamics are considered. It is shown that such reductions with respect to helical, axial, and translational symmetries have conservation laws which yield local potential variables. This leads to showing that the well-known Johnson-Frieman-Kruskal-Oberman (JFKO) and Bragg-Hawthorne (Grad-Shafranov) equations are locally related to the corresponding helically and axially symmetric PDE systems of fluid/plasma dynamics. For the axially symmetric case, local symmetry classifications and arising invariant solutions are compared for the original PDE system and the Bragg-Hawthorne (potential) equation. The potential equation is shown to have additional symmetries, denoted as restricted symmetries. Restricted symmetries leave invariant a family of solutions of a given PDE system but not the whole solution manifold, and hence are not symmetries of the given PDE system. Corresponding reductions are shown to yield solutions, which are not obtained as invariant solutions from local symmetry reduction.

3D MR flow analysis in realistic rapid-prototyping model systems of the thoracic aorta: comparison with in vivo data and computational fluid dynamics in identical vessel geometries.

PubMed

Canstein, C; Cachot, P; Faust, A; Stalder, A F; Bock, J; Frydrychowicz, A; Küffer, J; Hennig, J; Markl, M

2008-03-01

The knowledge of local vascular anatomy and function in the human body is of high interest for the diagnosis and treatment of cardiovascular disease. A comprehensive analysis of the hemodynamics in the thoracic aorta is presented based on the integration of flow-sensitive 4D MRI with state-of-the-art rapid prototyping technology and computational fluid dynamics (CFD). Rapid prototyping was used to transform aortic geometries as measured by contrast-enhanced MR angiography into realistic vascular models with large anatomical coverage. Integration into a flow circuit with patient-specific pulsatile in-flow conditions and application of flow-sensitive 4D MRI permitted detailed analysis of local and global 3D flow dynamics in a realistic vascular geometry. Visualization of characteristic 3D flow patterns and quantitative comparisons of the in vitro experiments with in vivo data and CFD simulations in identical vascular geometries were performed to evaluate the accuracy of vascular model systems. The results indicate the potential of such patient-specific model systems for detailed experimental simulation of realistic vascular hemodynamics. Further studies are warranted to examine the influence of refined boundary conditions of the human circulatory system such as fluid-wall interaction and their effect on normal and pathological blood flow characteristics associated with vascular geometry. (c) 2008 Wiley-Liss, Inc.

Quantification and Control of Wall Effects in Porous Media Experiments

NASA Astrophysics Data System (ADS)

Roth, E. J.; Mays, D. C.; Neupauer, R.; Crimaldi, J. P.

2017-12-01

Fluid flow dynamics in porous media are dominated by media heterogeneity. This heterogeneity can create preferential pathways in which local seepage velocities dwarf system seepage velocities, further complicating an already incomplete understanding of dispersive processes. In physical models of porous media flows, apparatus walls introduce preferential flow paths (i.e., wall effects) that may overwhelm other naturally occurring preferential pathways within the apparatus, leading to deceptive results. We used planar laser-induced fluorescence (PLIF) in conjunction with refractive index matched (RIM) porous media and pore fluid to observe fluid dynamics in the porous media, with particular attention to the region near the apparatus walls in a 17 cm x 8 cm x 7 cm uniform flow cell. Hexagonal close packed spheres were used to create an isotropic, homogenous porous media field in the interior of the apparatus. Visualization of the movement of a fluorescent dye revealed the influence of the wall in creating higher permeability preferential flow paths in an otherwise homogenous media packing. These preferential flow paths extended approximately one half of one sphere diameter from the wall for homogenously packed regions, with a quickly diminishing effect on flow dynamics for homogenous media adjacent to the preferential pathway, but with major influence on flow dynamics for adjoining heterogeneous regions. Multiple approaches to mitigate wall effects were investigated, and a modified wall was created such that the fluid dynamics near the wall mimics the fluid dynamics within the homogenous porous media. This research supports the design of a two-dimensional experimental apparatus that will simulate engineered pumping schemes for use in contaminant remediation. However, this research could benefit the design of fixed bed reactors or other engineering challenges in which vessel walls contribute to unwanted preferential flow.

Conserving the linear momentum in stochastic dynamics: Dissipative particle dynamics as a general strategy to achieve local thermostatization in molecular dynamics simulations.

PubMed

Passler, Peter P; Hofer, Thomas S

2017-02-15

Stochastic dynamics is a widely employed strategy to achieve local thermostatization in molecular dynamics simulation studies; however, it suffers from an inherent violation of momentum conservation. Although this short-coming has little impact on structural and short-time dynamic properties, it can be shown that dynamics in the long-time limit such as diffusion is strongly dependent on the respective thermostat setting. Application of the methodically similar dissipative particle dynamics (DPD) provides a simple, effective strategy to ensure the advantages of local, stochastic thermostatization while at the same time the linear momentum of the system remains conserved. In this work, the key parameters to employ the DPD thermostats in the framework of periodic boundary conditions are investigated, in particular the dependence of the system properties on the size of the DPD-region as well as the treatment of forces near the cutoff. Structural and dynamical data for light and heavy water as well as a Lennard-Jones fluid have been compared to simulations executed via stochastic dynamics as well as via use of the widely employed Nose-Hoover chain and Berendsen thermostats. It is demonstrated that a small size of the DPD region is sufficient to achieve local thermalization, while at the same time artifacts in the self-diffusion characteristic for stochastic dynamics are eliminated. © 2016 Wiley Periodicals, Inc. © 2016 Wiley Periodicals, Inc.

Methodology Development of a Gas-Liquid Dynamic Flow Regime Transition Model

NASA Astrophysics Data System (ADS)

Doup, Benjamin Casey

Current reactor safety analysis codes, such as RELAP5, TRACE, and CATHARE, use flow regime maps or flow regime transition criteria that were developed for static fully-developed two-phase flows to choose interfacial transfer models that are necessary to solve the two-fluid model. The flow regime is therefore difficult to identify near the flow regime transitions, in developing two-phase flows, and in transient two-phase flows. Interfacial area transport equations were developed to more accurately predict the dynamic nature of two-phase flows. However, other model coefficients are still flow regime dependent. Therefore, an accurate prediction of the flow regime is still important. In the current work, the methodology for the development of a dynamic flow regime transition model that uses the void fraction and interfacial area concentration obtained by solving three-field the two-fluid model and two-group interfacial area transport equation is investigated. To develop this model, detailed local experimental data are obtained, the two-group interfacial area transport equations are revised, and a dynamic flow regime transition model is evaluated using a computational fluid dynamics model. Local experimental data is acquired for 63 different flow conditions in bubbly, cap-bubbly, slug, and churn-turbulent flow regimes. The measured parameters are the group-1 and group-2 bubble number frequency, void fraction, interfacial area concentration, and interfacial bubble velocities. The measurements are benchmarked by comparing the prediction of the superficial gas velocities, determined using the local measurements with those determined from volumetric flow rate measurements and the agreement is generally within +/-20%. The repeatability four-sensor probe construction process is within +/-10%. The repeatability of the measurement process is within +/-7%. The symmetry of the test section is examined and the average agreement is within +/-5.3% at z/D = 10 and +/-3.4% at z/D = 32. Revised source/sink terms for the two-group interfacial area transport equations are derived and fit to area-averaged experimental data to determine new model coefficients. The average agreement between this model and the experiment data for the void fraction and interfacial area concentration is 10.6% and 15.7%, respectively. This revised two-group interfacial area transport equation and the three-field two-fluid model are used to solve for the group-1 and group-2 interfacial area concentration and void fraction. These values and a dynamic flow regime transition model are used to classify the flow regimes. The flow regimes determined using this model are compared with the flow regimes based on the experimental data and on a flow regime map using Mishima and Ishii's (1984) transition criteria. The dynamic flow regime transition model is shown to predict the flow regimes dynamically and has improved the prediction of the flow regime over that using a flow regime map. Safety codes often employ the one-dimensional two-fluid model to model two-phase flows. The area-averaged relative velocity correlation necessary to close this model is derived from the drift flux model. The effects of the necessary assumptions used to derive this correlation are investigated using local measurements and these effects are found to have a limited impact on the prediction of the area-averaged relative velocity.

Smoothed dissipative particle dynamics with angular momentum conservation

NASA Astrophysics Data System (ADS)

Müller, Kathrin; Fedosov, Dmitry A.; Gompper, Gerhard

2015-01-01

Smoothed dissipative particle dynamics (SDPD) combines two popular mesoscopic techniques, the smoothed particle hydrodynamics and dissipative particle dynamics (DPD) methods, and can be considered as an improved dissipative particle dynamics approach. Despite several advantages of the SDPD method over the conventional DPD model, the original formulation of SDPD by Español and Revenga (2003) [9], lacks angular momentum conservation, leading to unphysical results for problems where the conservation of angular momentum is essential. To overcome this limitation, we extend the SDPD method by introducing a particle spin variable such that local and global angular momentum conservation is restored. The new SDPD formulation (SDPD+a) is directly derived from the Navier-Stokes equation for fluids with spin, while thermal fluctuations are incorporated similarly to the DPD method. We test the new SDPD method and demonstrate that it properly reproduces fluid transport coefficients. Also, SDPD with angular momentum conservation is validated using two problems: (i) the Taylor-Couette flow with two immiscible fluids and (ii) a tank-treading vesicle in shear flow with a viscosity contrast between inner and outer fluids. For both problems, the new SDPD method leads to simulation predictions in agreement with the corresponding analytical theories, while the original SDPD method fails to capture properly physical characteristics of the systems due to violation of angular momentum conservation. In conclusion, the extended SDPD method with angular momentum conservation provides a new approach to tackle fluid problems such as multiphase flows and vesicle/cell suspensions, where the conservation of angular momentum is essential.

Unified description of the slip phenomena in sheared polymer films: A molecular dynamics study

NASA Astrophysics Data System (ADS)

Priezjev, Nikolai

2010-03-01

The dynamic behavior of the slip length in shear flow of polymer melts past atomically smooth surfaces is investigated using MD simulations. The polymer melt was modeled as a collection of FENE-LJ bead-spring chains. We consider shear flow conditions at low pressures and weak wall-fluid interaction energy so that fluid velocity profiles are linear throughout the channel at all shear rates examined. In agreement with earlier studies we confirm that for shear- thinning fluids the slip length passes through a local minimum at low shear rates and then increases rapidly at higher shear rates. We found that the rate dependence of the slip length depends on the lattice orientation at high shear rates. The MD results show that the ratio of slip length to viscosity follows a master curve when plotted as a function of a single variable that depends on the structure factor, contact density and temperature of the first fluid layer near the solid wall. The universal dependence of the slip length holds for a number of parameters of the interface: fluid density and structure (chain length), wall-fluid interaction energy, wall density, lattice orientation, thermal or solid walls.

Effect of surface roughness and size of beam on squeeze-film damping—Molecular dynamics simulation study

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

Kim, Hojin; Strachan, Alejandro

2015-11-28

We use large-scale molecular dynamics (MD) to characterize fluid damping between a substrate and an approaching beam. We focus on the near contact regime where squeeze film (where fluid gap is comparable to the mean free path of the gas molecules) and many-body effects in the fluid become dominant. The MD simulations provide explicit description of many-body and non-equilibrium processes in the fluid as well as the surface topography. We study how surface roughness and beam width increases the damping coefficient due to their effect on fluid mobility. We find that the explicit simulations are in good agreement with priormore » direct simulation Monte Carlo results except at near-contact conditions where many-body effects in the compressed fluid lead the increased damping and weaker dependence on beam width. We also show that velocity distributions near the beam edges and for short gaps deviate from the Boltzmann distribution indicating a degree of local non-equilibrium. These results will be useful to parameterize compact models used for microsystem device-level simulations and provide insight into mesoscale simulations of near-contact damping.« less

Vortex formation through inertial wave focusing

NASA Astrophysics Data System (ADS)

Duran-Matute, Matias; Flor, Jan-Bert; Godeferd, Fabien

2011-11-01

We present a novel experimental and numerical study on the formation of columnar vortical structures by inertial waves in a rotating fluid. Two inertial-wave cones are generated by a vertically oscillating torus in a fluid in solid body rotation At the tip of the cones, there is a singular point towards which the energy of the waves gets focused. The particularity of this configuration, as compared to those of previous experiments (e.g. oscillating sphere or disc), is that the singular point's position within the fluid leads to complex non-linear wave interaction, which may lead to the formation of a localized vortex that expands in the vertical in the form of a Taylor column. Using detailed PIV measurements we consider the flow evolution from the localized wave overturning motion to the Taylor column formation as well as the inertial wave dynamics during this process, The results are discussed in the context of turbulence in rotating fluids. We acknowledge financial support from projects ANR ANISO and CIBLE.

Role of string-like collective atomic motion on diffusion and structural relaxation in glass forming Cu-Zr alloys

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

Zhang, Hao; Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta T6G 2V4; Zhong, Cheng

2015-04-28

We investigate Cu-Zr liquid alloys using molecular dynamics simulation and well-accepted embedded atom method potentials over a wide range of chemical composition and temperature as model metallic glass-forming (GF) liquids. As with other types of GF materials, the dynamics of these complex liquids are characterized by “dynamic heterogeneity” in the form of transient polymeric clusters of highly mobile atoms that are composed in turn of atomic clusters exhibiting string-like cooperative motion. In accordance with the string model of relaxation, an extension of the Adam-Gibbs (AG) model, changes in the activation free energy ΔG{sub a} with temperature of both the Cumore » and Zr diffusion coefficients D, and the alpha structural relaxation time τ{sub α} can be described to a good approximation by changes in the average string length, L. In particular, we confirm that the strings are a concrete realization of the abstract “cooperatively rearranging regions” of AG. We also find coexisting clusters of relatively “immobile” atoms that exhibit predominantly icosahedral local packing rather than the low symmetry packing of “mobile” atoms. These two distinct types of dynamic heterogeneity are then associated with different fluid structural states. Glass-forming liquids are thus analogous to polycrystalline materials where the icosahedrally packed regions correspond to crystal grains, and the strings reside in the relatively disordered grain boundary-like regions exterior to these locally well-ordered regions. A dynamic equilibrium between localized (“immobile”) and wandering (“mobile”) particles exists in the liquid so that the dynamic heterogeneity can be considered to be type of self-assembly process. We also characterize changes in the local atomic free volume in the course of string-like atomic motion to better understand the initiation and propagation of these fluid excitations.« less

Creep cavitation can establish a dynamic granular fluid pump in ductile shear zones.

PubMed

Fusseis, F; Regenauer-Lieb, K; Liu, J; Hough, R M; De Carlo, F

2009-06-18

The feedback between fluid migration and rock deformation in mid-crustal shear zones is acknowledged as being critical for earthquake nucleation, the initiation of subduction zones and the formation of mineral deposits. The importance of this poorly understood feedback is further highlighted by evidence for shear-zone-controlled advective flow of fluids in the ductile lower crust and the recognition that deformation-induced grain-scale porosity is a key to large-scale geodynamics. Fluid migration in the middle crust cannot be explained in terms of classical concepts. The environment is considered too hot for a dynamic fracture-sustained permeability as in the upper crust, and fluid pathways are generally too deformed to be controlled by equilibrium wetting angles that apply to hotter, deeper environments. Here we present evidence that mechanical and chemical potentials control a syndeformational porosity generation in mid-crustal shear zones. High-resolution synchrotron X-ray tomography and scanning electron microscopy observations allow us to formulate a model for fluid migration in shear zones where a permeable porosity is dynamically created by viscous grain-boundary sliding, creep cavitation, dissolution and precipitation. We propose that syndeformational fluid migration in our 'granular fluid pump' model is a self-sustained process controlled by the explicit role of the rate of entropy production of the underlying irreversible mechanical and chemical microprocesses. The model explains fluid transfer through the middle crust, where strain localization in the creep regime is required for plate tectonics, the formation of giant ore deposits, mantle degassing and earthquake nucleation. Our findings provide a key component for the understanding of creep instabilities in the middle crust.

Modeling the formation of porphyry-copper ores

USGS Publications Warehouse

Ingebritsen, Steven E.

2012-01-01

Porphyry-copper ore systems, the source of much of the world's copper and molybdenum, form when metal-bearing fluids are expelled from shallow, degassing magmas. On page 1613 of this issue, Weis et al. (1) demonstrate that self-organizing processes focus metal deposition. Specifically, their simulation studies indicate that ores develop as consequences of dynamic variations in rock permeability driven by injection of volatile species from rising magmas. Scenarios with a static permeability structure could not reproduce key field observations, whereas dynamic permeability responses to magmatic-fluid injection localized a metal-precipitation front where enrichment by a factor of 103 could be achieved [for an overview of their numerical-simulation model CSMP++, see (2)].

Dynamics of Active Microfilaments

NASA Astrophysics Data System (ADS)

Ling, Feng; Guo, Hanliang; Kanso, Eva

2017-11-01

Soft elastic filaments are ubiquitous in natural and artificial systems at various length scales, and their interactions within and between filaments and their environments provide a persistent source of curiosity due to both the complexity of their behaviors and the relative mathematical simplicity of their structures. Specifically, a deeper understanding of the dynamic characteristics of microscopic filaments in viscous fluids is relevant to many biophysical and physiological processes. Here we start with the Cosserat model that allows all six possible modes of deformation for an elastic rod, and focus on the case of inextensible filaments submerged in viscous fluids by ignoring inertial effects and using local resistive force theory for fluid-filament interactions. We verify our simulations against special analytic solutions and present some results on the active internal control of cilia and flagella motion. We conclude by commenting on the utility of this general framework for studying other cellular and sub-cellular physical processes such as systems involving protein filaments.

Numerical modeling of local scour around hydraulic structure in sandy beds by dynamic mesh method

NASA Astrophysics Data System (ADS)

Fan, Fei; Liang, Bingchen; Bai, Yuchuan; Zhu, Zhixia; Zhu, Yanjun

2017-10-01

Local scour, a non-negligible factor in hydraulic engineering, endangers the safety of hydraulic structures. In this work, a numerical model for simulating local scour was constructed, based on the open source code computational fluid dynamics model OpenFOAM. We consider both the bedload and suspended load sediment transport in the scour model and adopt the dynamic mesh method to simulate the evolution of the bed elevation. We use the finite area method to project data between the three-dimensional flow model and the two-dimensional (2D) scour model. We also improved the 2D sand slide method and added it to the scour model to correct the bed bathymetry when the bed slope angle exceeds the angle of repose. Moreover, to validate our scour model, we conducted and compared the results of three experiments with those of the developed model. The validation results show that our developed model can reliably simulate local scour.

Local stability of galactic discs in modified dynamics

NASA Astrophysics Data System (ADS)

Shenavar, Hossein; Ghafourian, Neda

2018-04-01

The local stability of stellar and fluid discs, under a new modified dynamical model, is surveyed by using WKB approximation. The exact form of the modified Toomre criterion is derived for both types of systems and it is shown that the new model is, in all situations, more locally stable than Newtonian model. In addition, it has been proved that the central surface density of the galaxies plays an important role in the local stability in the sense that low surface brightness (LSB) galaxies are more stable than high surface brightness (HSBs). Furthermore, the growth rate in the new model is found to be lower than the Newtonian one. We found that, according to this model, the local instability is related to the ratio of surface density of the disc to a critical surface density Σcrit. We provide observational evidence to support this result based on star formation rate in HSBs and LSBs.

STARDUST-U experiments on fluid-dynamic conditions affecting dust mobilization during LOVAs

NASA Astrophysics Data System (ADS)

Poggi, L. A.; Malizia, A.; Ciparisse, J. F.; Tieri, F.; Gelfusa, M.; Murari, A.; Del Papa, C.; Giovannangeli, I.; Gaudio, P.

2016-07-01

Since 2006 the Quantum Electronics and Plasma Physics (QEP) Research Group together with ENEA FusTech of Frascati have been working on dust re-suspension inside tokamaks and its potential capability to jeopardize the integrity of future fusion nuclear plants (i.e. ITER or DEMO) and to be a risk for the health of the operators. Actually, this team is working with the improved version of the "STARDUST" facility, i.e. "STARDUST-Upgrade". STARDUST-U facility has four new air inlet ports that allow the experimental replication of Loss of Vacuum Accidents (LOVAs). The experimental campaign to detect the different pressurization rates, local air velocity, temperature, have been carried out from all the ports in different accident conditions and the principal results will be analyzed and compared with the numerical simulations obtained through a CFD (Computational Fluid Dynamic) code. This preliminary thermo fluid-dynamic analysis of the accident is crucial for numerical model development and validation, and for the incoming experimental campaign of dust resuspension inside STARDUST-U due to well-defined accidents presented in this paper.

Characterization of local fluid flow in 3D porous construct characterized by Fourier domain Doppler optical coherence tomography

NASA Astrophysics Data System (ADS)

Bagnaninchi, P. O.; Yang, Y.; El Haj, A.; Hinds, M. T.; Wang, R. K.

2007-02-01

In order to achieve functional tissue with the correct biomechanical properties it is critical to stimulate mechanically the cells. Perfusion bioreactor induces fluid shear stress that has been well characterized for two-dimensional culture where both simulation and experimental data are available. However these results can't be directly translated to tissue engineering that makes use of complex three-dimensional porous scaffold. Moreover, stimulated cells produce extensive extra-cellular matrix (ECM) that alter dramatically the micro-architecture of the constructs, changing the local flow dynamic. In this study a Fourier domain Doppler optical coherent tomography (FD-DOCT) system working at 1300nm with a bandwidth of 50nm has been used to determine the local flow rate inside different types of porous scaffolds used in tissue engineering. Local flow rates can then be linearly related, for Newtonian fluid, to the fluid shear stress occurring on the pores wall. Porous chitosan scaffolds (\\fgr 1.5mm x 3mm) with and without a central 250 μm microchannel have been produced by a freeze-drying technique. This techniques allow us to determine the actual shear stress applied to the cells and to optimise the input flow rate consequently, but also to relate the change of the flow distribution to the amount of ECM production allowing the monitoring of tissue formation.

A Smoluchowski model of crystallization dynamics of small colloidal clusters

NASA Astrophysics Data System (ADS)

Beltran-Villegas, Daniel J.; Sehgal, Ray M.; Maroudas, Dimitrios; Ford, David M.; Bevan, Michael A.

2011-10-01

We investigate the dynamics of colloidal crystallization in a 32-particle system at a fixed value of interparticle depletion attraction that produces coexisting fluid and solid phases. Free energy landscapes (FELs) and diffusivity landscapes (DLs) are obtained as coefficients of 1D Smoluchowski equations using as order parameters either the radius of gyration or the average crystallinity. FELs and DLs are estimated by fitting the Smoluchowski equations to Brownian dynamics (BD) simulations using either linear fits to locally initiated trajectories or global fits to unbiased trajectories using Bayesian inference. The resulting FELs are compared to Monte Carlo Umbrella Sampling results. The accuracy of the FELs and DLs for modeling colloidal crystallization dynamics is evaluated by comparing mean first-passage times from BD simulations with analytical predictions using the FEL and DL models. While the 1D models accurately capture dynamics near the free energy minimum fluid and crystal configurations, predictions near the transition region are not quantitatively accurate. A preliminary investigation of ensemble averaged 2D order parameter trajectories suggests that 2D models are required to capture crystallization dynamics in the transition region.

THREE-DIMENSIONAL COMPUTATIONAL FLUID DYNAMICS SIMULATIONS OF LOCAL PARTICLE DEPOSITION PATTERNS IN LUNG AIRWAYS

EPA Science Inventory

EPA has identified respirable particulate matter (PM) as a significant threat to human health, particularly in the elderly, in children, and in persons with respiratory disease. However, deposition of PM in the respiratory system is highly variable, depending upon particle chara...

Finite-time barriers to reaction front propagation

NASA Astrophysics Data System (ADS)

Locke, Rory; Mahoney, John; Mitchell, Kevin

2015-11-01

Front propagation in advection-reaction-diffusion systems gives rise to rich geometric patterns. It has been shown for time-independent and time-periodic fluid flows that invariant manifolds, termed burning invariant manifolds (BIMs), serve as one-sided dynamical barriers to the propagation of reaction front. More recently, theoretical work has suggested that one-sided barriers, termed burning Lagrangian Coherent structures (bLCSs), exist for fluid velocity data prescribed over a finite time interval, with no assumption on the time-dependence of the flow. In this presentation, we use a time-varying fluid ``wind'' in a double-vortex channel flow to demonstrate that bLCSs form the (locally) most attracting or repelling fronts.

Non-local features of a hydrodynamic pilot-wave system

NASA Astrophysics Data System (ADS)

Nachbin, Andre; Couchman, Miles; Bush, John

2016-11-01

A droplet walking on the surface of a vibrating fluid bath constitutes a pilot-wave system of the form envisaged for quantum dynamics by Louis de Broglie: a particle moves in resonance with its guiding wave field. We here present an examination of pilot-wave hydrodynamics in a confined domain. Specifically, we present a one-dimensional water wave model that describes droplets walking in single and multiple cavities. The cavities are separated by a submerged barrier, and so allow for the study of tunneling. They also highlight the non-local dynamical features arising due to the spatially-extended wave field. Results from computational simulations are complemented by laboratory experiments.

Original analytical model of the hydrodynamic loads applied on the half-bridge of a circular settling tank

NASA Astrophysics Data System (ADS)

Oanta, Emil M.; Dascalescu, Anca-Elena; Sabau, Adrian

2016-12-01

The paper presents an original analytical model of the hydrodynamic loads applied on the half-bridge of a circular settling tank. The calculus domain is defined using analytical geometry and the calculus of the local dynamic pressure is based on the radius from the center of the settling tank to the current area, i.e. the relative velocity of the fluid and the depth where the current area is located, i.e. the density of the fluid. Calculus of the local drag forces uses the discrete frontal cross sectional areas of the submerged structure in contact with the fluid. In the last stage is performed the reduction of the local drag forces in the appropriate points belonging to the main beam. This class of loads is producing the flexure of the main beam in a horizontal plane and additional twisting moments along this structure. Taking into account the hydrodynamic loads, the results of the theoretical models, i.e. the analytical model and the finite element model, may have an increased accuracy.

SSME main combustion chamber and nozzle flowfield analysis

NASA Technical Reports Server (NTRS)

Farmer, R. C.; Wang, T. S.; Smith, S. D.; Prozan, R. J.

1986-01-01

An investigation is presented of the computational fluid dynamics (CFD) tools which would accurately analyze main combustion chamber and nozzle flow. The importance of combustion phenomena and local variations in mixture ratio are fully appreciated; however, the computational aspects of the gas dynamics involved were the sole issues addressed. The CFD analyses made are first compared with conventional nozzle analyses to determine the accuracy for steady flows, and then transient analyses are discussed.

Mesh quality control for multiply-refined tetrahedral grids

NASA Technical Reports Server (NTRS)

Biswas, Rupak; Strawn, Roger

1994-01-01

A new algorithm for controlling the quality of multiply-refined tetrahedral meshes is presented in this paper. The basic dynamic mesh adaption procedure allows localized grid refinement and coarsening to efficiently capture aerodynamic flow features in computational fluid dynamics problems; however, repeated application of the procedure may significantly deteriorate the quality of the mesh. Results presented show the effectiveness of this mesh quality algorithm and its potential in the area of helicopter aerodynamics and acoustics.

High strain rate behavior of saturated and non-saturated sandstone: implications for earthquake mechanisms.

NASA Astrophysics Data System (ADS)

Aben, F. M.; Doan, M. L.; Gratier, J. P.; Renard, F.

2015-12-01

Damage zones of active faults control their resistance to rupture and transport properties. Hence, knowing the damage's origin is crucial to shed light on the (paleo)seismic behavior of the fault. Coseismic damage in the damage zone occurs by stress-wave loading of a passing earthquake rupture tip, resulting in dynamic (high strain rate) loading and subsequent dynamic fracturing or pulverization. Recently, interest in this type of damage has increased and several experimental studies were performed on dry rock specimens to search for pulverization-controlling parameters. However, the influence of fluids in during dynamic loading needs to be constrained. Hence, we have performed compressional dynamic loading experiments on water saturated and oven dried Vosges sandstone samples using a Split Hopkinson Pressure Bar apparatus. Due to the high porosity in these rocks, close to 20%, the effect of fluids should be clear. Afterwards, microstructural analyses have been applied on thin sections. Water saturated samples reveal dynamic mechanical behavior that follows linear poro-elasticity for undrained conditions: the peak strength of the sample decreases by 30-50% and the accumulated strain increases relative to the dry samples that were tested under similar conditions. The mechanical behavior of partially saturated samples falls in between. Microstructural studies on thin section show that fractures are restricted to some quartz grains while other quartz grains remain intact, similar to co-seismically damaged sandstones observed in the field. Most deformation is accommodated by inter-granular processes, thereby appointing an important role to the cement matrix in between grains. Intra-granular fracture damage is highest for the saturated samples. The presence of pore fluids in the rocks lower the dynamic peak strength, especially since fast dynamic loading does not allow for time-dependent fluid dissipation. Thus, fluid-saturated rocks would show undrained mechanical behavior, creating local overpressure in the pore that breaks the inter-granular cement. This strength-decreasing effect provides an explanation for the presence of pulverized and coseismically damaged rocks at depth and extends the range of dynamic stress where dynamic damage can occur in fault zones.

Multiphase fluid-solid coupled analysis of shock-bubble-stone interaction in shockwave lithotripsy.

PubMed

Wang, Kevin G

2017-10-01

A novel multiphase fluid-solid-coupled computational framework is applied to investigate the interaction of a kidney stone immersed in liquid with a lithotripsy shock wave (LSW) and a gas bubble near the stone. The main objective is to elucidate the effects of a bubble in the shock path to the elastic and fracture behaviors of the stone. The computational framework couples a finite volume 2-phase computational fluid dynamics solver with a finite element computational solid dynamics solver. The surface of the stone is represented as a dynamic embedded boundary in the computational fluid dynamics solver. The evolution of the bubble surface is captured by solving the level set equation. The interface conditions at the surfaces of the stone and the bubble are enforced through the construction and solution of local fluid-solid and 2-fluid Riemann problems. This computational framework is first verified for 3 example problems including a 1D multimaterial Riemann problem, a 3D shock-stone interaction problem, and a 3D shock-bubble interaction problem. Next, a series of shock-bubble-stone-coupled simulations are presented. This study suggests that the dynamic response of a bubble to LSW varies dramatically depending on its initial size. Bubbles with an initial radius smaller than a threshold collapse within 1 μs after the passage of LSW, whereas larger bubbles do not. For a typical LSW generated by an electrohydraulic lithotripter (p max  = 35.0MPa, p min  =- 10.1MPa), this threshold is approximately 0.12mm. Moreover, this study suggests that a noncollapsing bubble imposes a negative effect on stone fracture as it shields part of the LSW from the stone. On the other hand, a collapsing bubble may promote fracture on the proximal surface of the stone, yet hinder fracture from stone interior. Copyright © 2016 John Wiley & Sons, Ltd.

Optimization of a hydrodynamic separator using a multiscale computational fluid dynamics approach.

PubMed

Schmitt, Vivien; Dufresne, Matthieu; Vazquez, Jose; Fischer, Martin; Morin, Antoine

2013-01-01

This article deals with the optimization of a hydrodynamic separator working on the tangential separation mechanism along a screen. The aim of this study is to optimize the shape of the device to avoid clogging. A multiscale approach is used. This methodology combines measurements and computational fluid dynamics (CFD). A local model enables us to observe the different phenomena occurring at the orifice scale, which shows the potential of expanded metal screens. A global model is used to simulate the flow within the device using a conceptual model of the screen (porous wall). After validation against the experimental measurements, the global model was used to investigate the influence of deflectors and disk plates in the structure.

Nature of local magma storage zones and geometry of conduit systems below balsatic eruption sites - Pu'u 'O'o, Kilauea East Rift, Hawaii, example

NASA Technical Reports Server (NTRS)

Wilson, Lionel; Head, James W., III

1988-01-01

The fluid dynamics of the well-documented eruptive episodes at Pu'u 'O'o, Kilauea are used to investigate quantitatively the size and shape of the shallow conduit system beneath the vent. The possible geometry of this region is considered. The dynamics of the eruptive episodes is used to place restrictions on the size and shape of the region and thermal calculations are used to show that the geometry is consistent with the region being the fluid residue of the partially cooled, major preepisode 1 dike. The Pu'u 'O'o example is used to illustrate some general properties of shallow magma storage zones.

Frozen reaction fronts in steady flows: A burning-invariant-manifold perspective

NASA Astrophysics Data System (ADS)

Mahoney, John R.; Li, John; Boyer, Carleen; Solomon, Tom; Mitchell, Kevin A.

2015-12-01

The dynamics of fronts, such as chemical reaction fronts, propagating in two-dimensional fluid flows can be remarkably rich and varied. For time-invariant flows, the front dynamics may simplify, settling in to a steady state in which the reacted domain is static, and the front appears "frozen." Our central result is that these frozen fronts in the two-dimensional fluid are composed of segments of burning invariant manifolds, invariant manifolds of front-element dynamics in x y θ space, where θ is the front orientation. Burning invariant manifolds (BIMs) have been identified previously as important local barriers to front propagation in fluid flows. The relevance of BIMs for frozen fronts rests in their ability, under appropriate conditions, to form global barriers, separating reacted domains from nonreacted domains for all time. The second main result of this paper is an understanding of bifurcations that lead from a nonfrozen state to a frozen state, as well as bifurcations that change the topological structure of the frozen front. Although the primary results of this study apply to general fluid flows, our analysis focuses on a chain of vortices in a channel flow with an imposed wind. For this system, we present both experimental and numerical studies that support the theoretical analysis developed here.

Global two-fluid turbulence simulations of L-H transitions and edge localized mode dynamics in the COMPASS-D tokamak

NASA Astrophysics Data System (ADS)

Thyagaraja, A.; Valovič, M.; Knight, P. J.

2010-04-01

It is shown that the transition from L-mode to H-mode regimes in tokamaks can be reproduced using a two-fluid, fully electromagnetic, plasma model when a suitable particle sink is added at the edge. Such a model is implemented in the CUTIE code [A. Thyagaraja et al., Eur. J. Mech. B/Fluids 23, 475 (2004)] and is illustrated on plasma parameters that mimic those in the COMPASS-D tokamak with electron cyclotron resonance heating [Fielding et al., Plasma Phys. Contr. Fusion 42, A191 (2000)]. In particular, it is shown that holding the heating power, current, and magnetic field constant and increasing the fuelling rate to raise the plasma density leads spontaneously to the formation of an edge transport barrier (ETB) which occurs going from low to higher density experimentally. In the following quiescent period in which the stored energy of the plasma rises linearly with time, a dynamical transition occurs in the simulation with the appearance of features resembling strong edge localized modes. The simulation qualitatively reproduces many features observed in the experiment. Its relative robustness suggests that some, at least of the observed characteristics of ETBs and L-H transitions, can be captured in the global electromagnetic turbulence model.

Quantum droplets of light in the presence of synthetic magnetic fields

NASA Astrophysics Data System (ADS)

Wilson, Kali; Westerberg, Niclas; Valiente, Manuel; Duncan, Callum; Wright, Ewan; Ohberg, Patrik; Faccio, Daniele

2017-04-01

Recently, quantum droplets have been demonstrated in dipolar Bose-Einstein condensates, where the long range (nonlocal) attractive interaction is counterbalanced by a local repulsive interaction. In this work, we investigate the formation of quantum droplets in a two-dimensional nonlocal fluid of light. Fluids of light allow us to control the geometry of the system, and thus introduce vorticity which in turn creates an artificial magnetic field for the quantum droplet. In a quantum fluid of light, the photons comprising the fluid are treated as a gas of interacting Bose-particles, where the nonlocal interaction comes from the nonlinearity inherent in the material, in our case an attractive third-order thermo-optical nonlinearity. In contrast to matter-wave droplets, photon fluid droplets are not stabilised by local particle-particle scattering, but from the quantum pressure itself, i.e., a balance between diffraction and the nonlocal nonlinearity. We will present a numerical and analytical investigation of the ground state of these droplets and of their subsequent dynamics under the influence of a self-induced artificial magnetic field, and discuss experimental work with the possibility to include artificial gauge interactions between droplets.

Dynamics of poroelastic foams

NASA Astrophysics Data System (ADS)

Forterre, Yoel; Sobac, Benjamin

2010-11-01

Soft poroelastic structures are widespread in biological tissues such as cartilaginous joints in bones, blood-filled placentae or plant organs. Here we investigate the dynamics of open elastic foams immersed in viscous fluids, as model soft poroelastic materials. The experiment consists in slowly compacting blocs of polyurethane solid foam embedded in silicon oil-tanks and studying their relaxation to equilibrium when the confining stress is suddenly released. Measurements of the local fluid pressure and foam velocity field are compared with a simple two-phase flow approach. For small initial compactions, the results show quantitative agreement with the classical diffusion theory of soil consolidation (Terzaghi, Biot). On the other hand, for large initial compactions, the dynamics exhibits long relaxation times and decompaction fronts, which are mainly controlled by the highly non-linear mechanical response of the foam. The analogy between this process and the evaporation of a polymer melt close to the glass transition will be briefly discussed.

Localized stress fluctuations drive shear thickening in dense suspensions

NASA Astrophysics Data System (ADS)

Rathee, Vikram; Blair, Daniel L.; Urbach, Jeffrey S.

2017-08-01

Dense particulate suspensions exhibit a dramatic increase in average viscosity above a critical, material-dependent shear stress. This thickening changes from continuous to discontinuous as the concentration is increased. Using direct measurements of spatially resolved surface stresses in the continuous thickening regime, we report the existence of clearly defined dynamic localized regions of substantially increased stress that appear intermittently at stresses above the critical stress. With increasing applied stress, these regions occupy an increasing fraction of the system, and the increase accounts quantitatively for the observed shear thickening. The regions represent high-viscosity fluid phases, with a size determined by the distance between the shearing surfaces and a viscosity that is nearly independent of shear rate but that increases rapidly with concentration. Thus, we find that continuous shear thickening arises from increasingly frequent localized discontinuous transitions between distinct fluid phases with widely differing viscosities.

Local-Scale Simulations of Nucleate Boiling on Micrometer Featured Surfaces: Preprint

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

Sitaraman, Hariswaran; Moreno, Gilberto; Narumanchi, Sreekant V

2017-08-03

A high-fidelity computational fluid dynamics (CFD)-based model for bubble nucleation of the refrigerant HFE7100 on micrometer-featured surfaces is presented in this work. The single-fluid incompressible Navier-Stokes equations, along with energy transport and natural convection effects are solved on a featured surface resolved grid. An a priori cavity detection method is employed to convert raw profilometer data of a surface into well-defined cavities. The cavity information and surface morphology are represented in the CFD model by geometric mesh deformations. Surface morphology is observed to initiate buoyancy-driven convection in the liquid phase, which in turn results in faster nucleation of cavities. Simulationsmore » pertaining to a generic rough surface show a trend where smaller size cavities nucleate with higher wall superheat. This local-scale model will serve as a self-consistent connection to larger device scale continuum models where local feature representation is not possible.« less

Local-Scale Simulations of Nucleate Boiling on Micrometer-Featured Surfaces

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

Sitaraman, Hariswaran; Moreno, Gilberto; Narumanchi, Sreekant V

2017-07-12

A high-fidelity computational fluid dynamics (CFD)-based model for bubble nucleation of the refrigerant HFE7100 on micrometer-featured surfaces is presented in this work. The single-fluid incompressible Navier-Stokes equations, along with energy transport and natural convection effects are solved on a featured surface resolved grid. An a priori cavity detection method is employed to convert raw profilometer data of a surface into well-defined cavities. The cavity information and surface morphology are represented in the CFD model by geometric mesh deformations. Surface morphology is observed to initiate buoyancy-driven convection in the liquid phase, which in turn results in faster nucleation of cavities. Simulationsmore » pertaining to a generic rough surface show a trend where smaller size cavities nucleate with higher wall superheat. This local-scale model will serve as a self-consistent connection to larger device scale continuum models where local feature representation is not possible.« less

Dynamics of Nanoparticles in Entangled Polymer Solutions

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

Nath, Pooja; Mangal, Rahul; Kohle, Ferdinand

The mean square displacement < r 2 > of nanoparticle probes dispersed in simple isotropic liquids and in polymer solutions is interrogated using fluorescence correlation spectroscopy and single-particle tracking (SPT) experiments. Probe dynamics in different regimes of particle diameter (d), relative to characteristic polymer length scales, including the correlation length (ξ), the entanglement mesh size (a), and the radius of gyration (R g), are investigated. In simple fluids and for polymer solutions in which d >> R g, long-time particle dynamics obey random-walk statistics < r 2 >:t, with the bulk zero-shear viscosity of the polymer solution determining the frictionalmore » resistance to particle motion. In contrast, in polymer solutions with d < R g, polymer molecules in solution exert noncontinuum resistances to particle motion and nanoparticle probes appear to interact hydrodynamically only with a local fluid medium with effective drag comparable to that of a solution of polymer chain segments with sizes similar to those of the nanoparticle probes. Under these conditions, the nanoparticles exhibit orders of magnitude faster dynamics than those expected from continuum predictions based on the Stokes–Einstein relation. SPT measurements further show that when d > a, nanoparticle dynamics transition from diffusive to subdiffusive on long timescales, reminiscent of particle transport in a field with obstructions. This last finding is in stark contrast to the nanoparticle dynamics observed in entangled polymer melts, where X-ray photon correlation spectroscopy measurements reveal faster but hyperdiffusive dynamics. As a result, we analyze these results with the help of the hopping model for particle dynamics in polymers proposed by Cai et al. and, on that basis, discuss the physical origins of the local drag experienced by the nanoparticles in entangled polymer solutions.« less

Dynamics of Nanoparticles in Entangled Polymer Solutions

DOE PAGES

Nath, Pooja; Mangal, Rahul; Kohle, Ferdinand; ...

2017-12-01

The mean square displacement < r 2 > of nanoparticle probes dispersed in simple isotropic liquids and in polymer solutions is interrogated using fluorescence correlation spectroscopy and single-particle tracking (SPT) experiments. Probe dynamics in different regimes of particle diameter (d), relative to characteristic polymer length scales, including the correlation length (ξ), the entanglement mesh size (a), and the radius of gyration (R g), are investigated. In simple fluids and for polymer solutions in which d >> R g, long-time particle dynamics obey random-walk statistics < r 2 >:t, with the bulk zero-shear viscosity of the polymer solution determining the frictionalmore » resistance to particle motion. In contrast, in polymer solutions with d < R g, polymer molecules in solution exert noncontinuum resistances to particle motion and nanoparticle probes appear to interact hydrodynamically only with a local fluid medium with effective drag comparable to that of a solution of polymer chain segments with sizes similar to those of the nanoparticle probes. Under these conditions, the nanoparticles exhibit orders of magnitude faster dynamics than those expected from continuum predictions based on the Stokes–Einstein relation. SPT measurements further show that when d > a, nanoparticle dynamics transition from diffusive to subdiffusive on long timescales, reminiscent of particle transport in a field with obstructions. This last finding is in stark contrast to the nanoparticle dynamics observed in entangled polymer melts, where X-ray photon correlation spectroscopy measurements reveal faster but hyperdiffusive dynamics. As a result, we analyze these results with the help of the hopping model for particle dynamics in polymers proposed by Cai et al. and, on that basis, discuss the physical origins of the local drag experienced by the nanoparticles in entangled polymer solutions.« less

A dissipative particle dynamics method for arbitrarily complex geometries

NASA Astrophysics Data System (ADS)

Li, Zhen; Bian, Xin; Tang, Yu-Hang; Karniadakis, George Em

2018-02-01

Dissipative particle dynamics (DPD) is an effective Lagrangian method for modeling complex fluids in the mesoscale regime but so far it has been limited to relatively simple geometries. Here, we formulate a local detection method for DPD involving arbitrarily shaped geometric three-dimensional domains. By introducing an indicator variable of boundary volume fraction (BVF) for each fluid particle, the boundary of arbitrary-shape objects is detected on-the-fly for the moving fluid particles using only the local particle configuration. Therefore, this approach eliminates the need of an analytical description of the boundary and geometry of objects in DPD simulations and makes it possible to load the geometry of a system directly from experimental images or computer-aided designs/drawings. More specifically, the BVF of a fluid particle is defined by the weighted summation over its neighboring particles within a cutoff distance. Wall penetration is inferred from the value of the BVF and prevented by a predictor-corrector algorithm. The no-slip boundary condition is achieved by employing effective dissipative coefficients for liquid-solid interactions. Quantitative evaluations of the new method are performed for the plane Poiseuille flow, the plane Couette flow and the Wannier flow in a cylindrical domain and compared with their corresponding analytical solutions and (high-order) spectral element solution of the Navier-Stokes equations. We verify that the proposed method yields correct no-slip boundary conditions for velocity and generates negligible fluctuations of density and temperature in the vicinity of the wall surface. Moreover, we construct a very complex 3D geometry - the "Brown Pacman" microfluidic device - to explicitly demonstrate how to construct a DPD system with complex geometry directly from loading a graphical image. Subsequently, we simulate the flow of a surfactant solution through this complex microfluidic device using the new method. Its effectiveness is demonstrated by examining the rich dynamics of surfactant micelles, which are flowing around multiple small cylinders and stenotic regions in the microfluidic device without wall penetration. In addition to stationary arbitrary-shape objects, the new method is particularly useful for problems involving moving and deformable boundaries, because it only uses local information of neighboring particles and satisfies the desired boundary conditions on-the-fly.

Hydrodynamically Coupled Brownian Dynamics: A coarse-grain particle-based Brownian dynamics technique with hydrodynamic interactions for modeling self-developing flow of polymer solutions

NASA Astrophysics Data System (ADS)

Ahuja, V. R.; van der Gucht, J.; Briels, W. J.

2018-01-01

We present a novel coarse-grain particle-based simulation technique for modeling self-developing flow of dilute and semi-dilute polymer solutions. The central idea in this paper is the two-way coupling between a mesoscopic polymer model and a phenomenological fluid model. As our polymer model, we choose Responsive Particle Dynamics (RaPiD), a Brownian dynamics method, which formulates the so-called "conservative" and "transient" pair-potentials through which the polymers interact besides experiencing random forces in accordance with the fluctuation dissipation theorem. In addition to these interactions, our polymer blobs are also influenced by the background solvent velocity field, which we calculate by solving the Navier-Stokes equation discretized on a moving grid of fluid blobs using the Smoothed Particle Hydrodynamics (SPH) technique. While the polymers experience this frictional force opposing their motion relative to the background flow field, our fluid blobs also in turn are influenced by the motion of the polymers through an interaction term. This makes our technique a two-way coupling algorithm. We have constructed this interaction term in such a way that momentum is conserved locally, thereby preserving long range hydrodynamics. Furthermore, we have derived pairwise fluctuation terms for the velocities of the fluid blobs using the Fokker-Planck equation, which have been alternatively derived using the General Equation for the Non-Equilibrium Reversible-Irreversible Coupling (GENERIC) approach in Smoothed Dissipative Particle Dynamics (SDPD) literature. These velocity fluctuations for the fluid may be incorporated into the velocity updates for our fluid blobs to obtain a thermodynamically consistent distribution of velocities. In cases where these fluctuations are insignificant, however, these additional terms may well be dropped out as they are in a standard SPH simulation. We have applied our technique to study the rheology of two different concentrations of our model linear polymer solutions. The results show that the polymers and the fluid are coupled very well with each other, showing no lag between their velocities. Furthermore, our results show non-Newtonian shear thinning and the characteristic flattening of the Poiseuille flow profile typically observed for polymer solutions.

Hydrodynamically Coupled Brownian Dynamics: A coarse-grain particle-based Brownian dynamics technique with hydrodynamic interactions for modeling self-developing flow of polymer solutions.

PubMed

Ahuja, V R; van der Gucht, J; Briels, W J

2018-01-21

We present a novel coarse-grain particle-based simulation technique for modeling self-developing flow of dilute and semi-dilute polymer solutions. The central idea in this paper is the two-way coupling between a mesoscopic polymer model and a phenomenological fluid model. As our polymer model, we choose Responsive Particle Dynamics (RaPiD), a Brownian dynamics method, which formulates the so-called "conservative" and "transient" pair-potentials through which the polymers interact besides experiencing random forces in accordance with the fluctuation dissipation theorem. In addition to these interactions, our polymer blobs are also influenced by the background solvent velocity field, which we calculate by solving the Navier-Stokes equation discretized on a moving grid of fluid blobs using the Smoothed Particle Hydrodynamics (SPH) technique. While the polymers experience this frictional force opposing their motion relative to the background flow field, our fluid blobs also in turn are influenced by the motion of the polymers through an interaction term. This makes our technique a two-way coupling algorithm. We have constructed this interaction term in such a way that momentum is conserved locally, thereby preserving long range hydrodynamics. Furthermore, we have derived pairwise fluctuation terms for the velocities of the fluid blobs using the Fokker-Planck equation, which have been alternatively derived using the General Equation for the Non-Equilibrium Reversible-Irreversible Coupling (GENERIC) approach in Smoothed Dissipative Particle Dynamics (SDPD) literature. These velocity fluctuations for the fluid may be incorporated into the velocity updates for our fluid blobs to obtain a thermodynamically consistent distribution of velocities. In cases where these fluctuations are insignificant, however, these additional terms may well be dropped out as they are in a standard SPH simulation. We have applied our technique to study the rheology of two different concentrations of our model linear polymer solutions. The results show that the polymers and the fluid are coupled very well with each other, showing no lag between their velocities. Furthermore, our results show non-Newtonian shear thinning and the characteristic flattening of the Poiseuille flow profile typically observed for polymer solutions.

Investigation of wind behaviour around high-rise buildings

NASA Astrophysics Data System (ADS)

Mat Isa, Norasikin; Fitriah Nasir, Nurul; Sadikin, Azmahani; Ariff Hairul Bahara, Jamil

2017-09-01

A study on the investigation of wind behaviour around the high-rise buildings is done through an experiment using a wind tunnel and computational fluid dynamics. High-rise buildings refer to buildings or structures that have more than 12 floors. Wind is invisible to the naked eye; thus, it is hard to see and analyse its flow around and over buildings without the use of proper methods, such as the use of wind tunnel and computational fluid dynamics software.The study was conducted on buildings located in Presint 4, Putrajaya, Malaysia which is the Ministry of Rural and Regional Development, Ministry of Information Communications and Culture, Ministry of Urban Wellbeing, Housing and Local Government and the Ministry of Women, Family, and Community by making scaled models of the buildings. The parameters in which this study is conducted on are, four different wind velocities used based on the seasonal monsoons, and wind direction. ANSYS Fluent workbench software is used to compute the simulations in order to achieve the objectives of this study. The data from the computational fluid dynamics are validated with the experiment done through the wind tunnel. From the results obtained through the use of the computation fluid dynamics, this study can identify the characteristics of wind around buildings, including boundary layer of the buildings, separation flow, wake region and etc. Then analyses is conducted on the occurance resulting from the wind that passes the buildings based on the velocity difference between before and after the wind passes the buildings.

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

Müller, Kathrin, E-mail: k.mueller@fz-juelich.de; Fedosov, Dmitry A., E-mail: d.fedosov@fz-juelich.de; Gompper, Gerhard, E-mail: g.gompper@fz-juelich.de

Smoothed dissipative particle dynamics (SDPD) combines two popular mesoscopic techniques, the smoothed particle hydrodynamics and dissipative particle dynamics (DPD) methods, and can be considered as an improved dissipative particle dynamics approach. Despite several advantages of the SDPD method over the conventional DPD model, the original formulation of SDPD by Español and Revenga (2003) [9], lacks angular momentum conservation, leading to unphysical results for problems where the conservation of angular momentum is essential. To overcome this limitation, we extend the SDPD method by introducing a particle spin variable such that local and global angular momentum conservation is restored. The new SDPDmore » formulation (SDPD+a) is directly derived from the Navier–Stokes equation for fluids with spin, while thermal fluctuations are incorporated similarly to the DPD method. We test the new SDPD method and demonstrate that it properly reproduces fluid transport coefficients. Also, SDPD with angular momentum conservation is validated using two problems: (i) the Taylor–Couette flow with two immiscible fluids and (ii) a tank-treading vesicle in shear flow with a viscosity contrast between inner and outer fluids. For both problems, the new SDPD method leads to simulation predictions in agreement with the corresponding analytical theories, while the original SDPD method fails to capture properly physical characteristics of the systems due to violation of angular momentum conservation. In conclusion, the extended SDPD method with angular momentum conservation provides a new approach to tackle fluid problems such as multiphase flows and vesicle/cell suspensions, where the conservation of angular momentum is essential.« less

Spectral analysis of the turbulent mixing of two fluids

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

Steinkamp, M.J.

1996-02-01

The authors describe a spectral approach to the investigation of fluid instability, generalized turbulence, and the interpenetration of fluids across an interface. The technique also applies to a single fluid with large variations in density. Departures of fluctuating velocity components from the local mean are far subsonic, but the mean Mach number can be large. Validity of the description is demonstrated by comparisons with experiments on turbulent mixing due to the late stages of Rayleigh-Taylor instability, when the dynamics become approximately self-similar in response to a constant body force. Generic forms for anisotropic spectral structure are described and used asmore » a basis for deriving spectrally integrated moment equations that can be incorporated into computer codes for scientific and engineering analyses.« less

Implementing a Loosely Coupled Fluid Structure Interaction Finite Element Model in PHASTA

NASA Astrophysics Data System (ADS)

Pope, David

Fluid Structure Interaction problems are an important multi-physics phenomenon in the design of aerospace vehicles and other engineering applications. A variety of computational fluid dynamics solvers capable of resolving the fluid dynamics exist. PHASTA is one such computational fluid dynamics solver. Enhancing the capability of PHASTA to resolve Fluid-Structure Interaction first requires implementing a structural dynamics solver. The implementation also requires a correction of the mesh used to solve the fluid equations to account for the deformation of the structure. This results in mesh motion and causes the need for an Arbitrary Lagrangian-Eulerian modification to the fluid dynamics equations currently implemented in PHASTA. With the implementation of both structural dynamics physics, mesh correction, and the Arbitrary Lagrangian-Eulerian modification of the fluid dynamics equations, PHASTA is made capable of solving Fluid-Structure Interaction problems.

General framework for fluctuating dynamic density functional theory

NASA Astrophysics Data System (ADS)

Durán-Olivencia, Miguel A.; Yatsyshin, Peter; Goddard, Benjamin D.; Kalliadasis, Serafim

2017-12-01

We introduce a versatile bottom-up derivation of a formal theoretical framework to describe (passive) soft-matter systems out of equilibrium subject to fluctuations. We provide a unique connection between the constituent-particle dynamics of real systems and the time evolution equation of their measurable (coarse-grained) quantities, such as local density and velocity. The starting point is the full Hamiltonian description of a system of colloidal particles immersed in a fluid of identical bath particles. Then, we average out the bath via Zwanzig’s projection-operator techniques and obtain the stochastic Langevin equations governing the colloidal-particle dynamics. Introducing the appropriate definition of the local number and momentum density fields yields a generalisation of the Dean-Kawasaki (DK) model, which resembles the stochastic Navier-Stokes description of a fluid. Nevertheless, the DK equation still contains all the microscopic information and, for that reason, does not represent the dynamical law of observable quantities. We address this controversial feature of the DK description by carrying out a nonequilibrium ensemble average. Adopting a natural decomposition into local-equilibrium and nonequilibrium contribution, where the former is related to a generalised version of the canonical distribution, we finally obtain the fluctuating-hydrodynamic equation governing the time-evolution of the mesoscopic density and momentum fields. Along the way, we outline the connection between the ad hoc energy functional introduced in previous DK derivations and the free-energy functional from classical density-functional theory. The resultant equation has the structure of a dynamical density-functional theory (DDFT) with an additional fluctuating force coming from the random interactions with the bath. We show that our fluctuating DDFT formalism corresponds to a particular version of the fluctuating Navier-Stokes equations, originally derived by Landau and Lifshitz. Our framework thus provides the formal apparatus for ab initio derivations of fluctuating DDFT equations capable of describing the dynamics of soft-matter systems in and out of equilibrium.

Experimental quantification of the fluid dynamics in blood-processing devices through 4D-flow imaging: A pilot study on a real oxygenator/heat-exchanger module.

PubMed

Piatti, Filippo; Palumbo, Maria Chiara; Consolo, Filippo; Pluchinotta, Francesca; Greiser, Andreas; Sturla, Francesco; Votta, Emiliano; Siryk, Sergii V; Vismara, Riccardo; Fiore, Gianfranco Beniamino; Lombardi, Massimo; Redaelli, Alberto

2018-02-08

The performance of blood-processing devices largely depends on the associated fluid dynamics, which hence represents a key aspect in their design and optimization. To this aim, two approaches are currently adopted: computational fluid-dynamics, which yields highly resolved three-dimensional data but relies on simplifying assumptions, and in vitro experiments, which typically involve the direct video-acquisition of the flow field and provide 2D data only. We propose a novel method that exploits space- and time-resolved magnetic resonance imaging (4D-flow) to quantify the complex 3D flow field in blood-processing devices and to overcome these limitations. We tested our method on a real device that integrates an oxygenator and a heat exchanger. A dedicated mock loop was implemented, and novel 4D-flow sequences with sub-millimetric spatial resolution and region-dependent velocity encodings were defined. Automated in house software was developed to quantify the complex 3D flow field within the different regions of the device: region-dependent flow rates, pressure drops, paths of the working fluid and wall shear stresses were computed. Our analysis highlighted the effects of fine geometrical features of the device on the local fluid-dynamics, which would be unlikely observed by current in vitro approaches. Also, the effects of non-idealities on the flow field distribution were captured, thanks to the absence of the simplifying assumptions that typically characterize numerical models. To the best of our knowledge, our approach is the first of its kind and could be extended to the analysis of a broad range of clinically relevant devices. Copyright © 2017 Elsevier Ltd. All rights reserved.

Marangoni-induced symmetry-breaking pattern selection on viscous fluids

NASA Astrophysics Data System (ADS)

Shen, Li; Denner, Fabian; Morgan, Neal; van Wachem, Berend; Dini, Daniele

2016-11-01

Symmetry breaking transitions on curved surfaces are found in a wide range of dissipative systems, ranging from asymmetric cell divisions to structure formation in thin films. Inherent within the nonlinearities are the associated curvilinear geometry, the elastic stretching, bending and the various fluid dynamical processes. We present a generalised Swift-Hohenberg pattern selection theory on a thin, curved and viscous films in the presence of non-trivial Marangoni effect. Testing the theory with experiments on soap bubbles, we observe the film pattern selection to mimic that of the elastic wrinkling morphology on a curved elastic bilayer in regions of slow viscous flow. By examining the local state of damping of surface capillary waves we attempt to establish an equivalence between the Marangoni fluid dynamics and the nonlinear elastic shell theory above the critical wavenumber of the instabilities and propose a possible explanation for the perceived elastic-fluidic duality. The authors acknowledge the financial support of the Shell University Technology Centre for fuels and lubricants.

High-temperature thermocline TES combining sensible and latent heat - CFD modeling and experimental validation

NASA Astrophysics Data System (ADS)

Zavattoni, Simone A.; Geissbühler, Lukas; Barbato, Maurizio C.; Zanganeh, Giw; Haselbacher, Andreas; Steinfeld, Aldo

2017-06-01

The concept of combined sensible/latent heat thermal energy storage (TES) has been exploited to mitigate an intrinsic thermocline TES systems drawback of heat transfer fluid outflow temperature reduction during discharging. In this study, the combined sensible/latent TES prototype under investigation is constituted by a packed bed of rocks and a small amount of encapsulated phase change material (AlSi12) as sensible heat and latent heat sections respectively. The thermo-fluid dynamics behavior of the combined TES prototype was analyzed by means of a computational fluid dynamics approach. Due to the small value of the characteristic vessel-to-particles diameter ratio, the effect of radial void-fraction variation, also known as channeling, was accounted for. Both the sensible and the latent heat sections of the storage were modeled as porous media under the assumption of local thermal non-equilibrium (LTNE). The commercial code ANSYS Fluent 15.0 was used to solve the model's constitutive conservation and transport equations obtaining a fairly good agreement with reference experimental measurements.

Fluid Dynamics of Jets with Applications to V/STOL.

DTIC Science & Technology

1982-01-01

Velocity coefficient Rr Local Reynolds number v Cinematic viscosity R Particle’s Reynolds number C1,C2 Non dimensional coordinates U0 Jet exit... realisme de la solution tant sur le plan fonctionnel que technologique, le recollement 6tant spontan# et stable et l’architecture extr~mement rustique. Les

Thermal Hydraulic Computational Fluid Dynamics Simulations and Experimental Investigation of Deformed Fuel Assemblies

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

Mays, Brian; Jackson, R. Brian

2017-03-08

The project, Toward a Longer Life Core: Thermal Hydraulic CFD Simulations and Experimental Investigation of Deformed Fuel Assemblies, DOE Project code DE-NE0008321, was a verification and validation project for flow and heat transfer through wire wrapped simulated liquid metal fuel assemblies that included both experiments and computational fluid dynamics simulations of those experiments. This project was a two year collaboration between AREVA, TerraPower, Argonne National Laboratory and Texas A&M University. Experiments were performed by AREVA and Texas A&M University. Numerical simulations of these experiments were performed by TerraPower and Argonne National Lab. Project management was performed by AREVA Federal Services.more » The first of a kind project resulted in the production of both local point temperature measurements and local flow mixing experiment data paired with numerical simulation benchmarking of the experiments. The project experiments included the largest wire-wrapped pin assembly Mass Index of Refraction (MIR) experiment in the world, the first known wire-wrapped assembly experiment with deformed duct geometries and the largest numerical simulations ever produced for wire-wrapped bundles.« less

Modelling deformation of partially melted rock using a poroviscoelastic rheology with dynamic power law viscosity

NASA Astrophysics Data System (ADS)

Simakin, A.; Ghassemi, A.

2005-03-01

A poroviscoelastic constitutive model is developed and used to study coupled rock deformation and fluid flow. The model allows the relaxation of both shear and symmetric components of the effective stress. Experimental results are usually interpreted in terms of the power law viscous material. However, in this work the effect of strain damage on viscosity is considered by treating the viscosity as a dynamic time-dependent parameter that varies proportionally to the second invariant of the strain rate. Healing is also taken into account so that the dynamic power law viscosity has a constant asymptotic at a given strain rate. The theoretical model is implemented in a finite element (FE) formulation that couples fluid flow and mechanical equilibrium equations. The FE method is applied to numerically study the triaxial compression of partially melted rocks at elevated PT conditions. It is found that the numerically calculated stress-strain curves demonstrate maxima similar to those observed in laboratory experiments. Also, the computed pattern of melt redistribution and strain localization at the contact between the rock sample and a stiff spacer is qualitatively similar to the experimental observations. The results also indicate that the matrix sensitivity to damage affects the scale of strain localization and melt redistribution.

Stability and stabilization of the lattice Boltzmann method

NASA Astrophysics Data System (ADS)

Brownlee, R. A.; Gorban, A. N.; Levesley, J.

2007-03-01

We revisit the classical stability versus accuracy dilemma for the lattice Boltzmann methods (LBM). Our goal is a stable method of second-order accuracy for fluid dynamics based on the lattice Bhatnager-Gross-Krook method (LBGK). The LBGK scheme can be recognized as a discrete dynamical system generated by free flight and entropic involution. In this framework the stability and accuracy analysis are more natural. We find the necessary and sufficient conditions for second-order accurate fluid dynamics modeling. In particular, it is proven that in order to guarantee second-order accuracy the distribution should belong to a distinguished surface—the invariant film (up to second order in the time step). This surface is the trajectory of the (quasi)equilibrium distribution surface under free flight. The main instability mechanisms are identified. The simplest recipes for stabilization add no artificial dissipation (up to second order) and provide second-order accuracy of the method. Two other prescriptions add some artificial dissipation locally and prevent the system from loss of positivity and local blowup. Demonstration of the proposed stable LBGK schemes are provided by the numerical simulation of a one-dimensional (1D) shock tube and the unsteady 2D flow around a square cylinder up to Reynolds number Rẽ20000 .

Working research codes into fluid dynamics education: a science gateway approach

NASA Astrophysics Data System (ADS)

Mason, Lachlan; Hetherington, James; O'Reilly, Martin; Yong, May; Jersakova, Radka; Grieve, Stuart; Perez-Suarez, David; Klapaukh, Roman; Craster, Richard V.; Matar, Omar K.

2017-11-01

Research codes are effective for illustrating complex concepts in educational fluid dynamics courses, compared to textbook examples, an interactive three-dimensional visualisation can bring a problem to life! Various barriers, however, prevent the adoption of research codes in teaching: codes are typically created for highly-specific `once-off' calculations and, as such, have no user interface and a steep learning curve. Moreover, a code may require access to high-performance computing resources that are not readily available in the classroom. This project allows academics to rapidly work research codes into their teaching via a minimalist `science gateway' framework. The gateway is a simple, yet flexible, web interface allowing students to construct and run simulations, as well as view and share their output. Behind the scenes, the common operations of job configuration, submission, monitoring and post-processing are customisable at the level of shell scripting. In this talk, we demonstrate the creation of an example teaching gateway connected to the Code BLUE fluid dynamics software. Student simulations can be run via a third-party cloud computing provider or a local high-performance cluster. EPSRC, UK, MEMPHIS program Grant (EP/K003976/1), RAEng Research Chair (OKM).

Entropy density of spacetime and the Navier-Stokes fluid dynamics of null surfaces

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

Padmanabhan, T.

2011-02-15

It has been known for several decades that Einstein's field equations, when projected onto a null surface, exhibit a structure very similar to the nonrelativistic Navier-Stokes equation. I show that this result arises quite naturally when gravitational dynamics is viewed as an emergent phenomenon. Extremizing the spacetime entropy density associated with the null surfaces leads to a set of equations which, when viewed in the local inertial frame, becomes identical to the Navier-Stokes equation. This is in contrast to the usual description of the Damour-Navier-Stokes equation in a general coordinate system, in which there appears a Lie derivative rather thanmore » a convective derivative. I discuss this difference, its importance, and why it is more appropriate to view the equation in a local inertial frame. The viscous force on fluid, arising from the gradient of the viscous stress-tensor, involves the second derivatives of the metric and does not vanish in the local inertial frame, while the viscous stress-tensor itself vanishes so that inertial observers detect no dissipation. We thus provide an entropy extremization principle that leads to the Damour-Navier-Stokes equation, which makes the hydrodynamical analogy with gravity completely natural and obvious. Several implications of these results are discussed.« less

Discovering Coherent Structures Using Local Causal States

NASA Astrophysics Data System (ADS)

Rupe, Adam; Crutchfield, James P.; Kashinath, Karthik; Prabhat, Mr.

2017-11-01

Coherent structures were introduced in the study of fluid dynamics and were initially defined as regions characterized by high levels of coherent vorticity, i.e. regions where instantaneously space and phase correlated vorticity are high. In a more general spatiotemporal setting, coherent structures can be seen as localized broken symmetries which persist in time. Building off the computational mechanics framework, which integrates tools from computation and information theory to capture pattern and structure in nonlinear dynamical systems, we introduce a theory of coherent structures, in the more general sense. Central to computational mechanics is the causal equivalence relation, and a local spatiotemporal generalization of it is used to construct the local causal states, which are utilized to uncover a system's spatiotemporal symmetries. Coherent structures are then identified as persistent, localized deviations from these symmetries. We illustrate how novel patterns and structures can be discovered in cellular automata and outline the path from them to laminar, transitional and turbulent flows. Funded by Intel through the Big Data Center at LBNL and the IPCC at UC Davis.

Dynamic hydraulic fluid stimulation regulated intramedullary pressure.

PubMed

Hu, Minyi; Serra-Hsu, Frederick; Bethel, Neville; Lin, Liangjun; Ferreri, Suzanne; Cheng, Jiqi; Qin, Yi-Xian

2013-11-01

Physical signals within the bone, i.e. generated from mechanical loading, have the potential to initiate skeletal adaptation. Strong evidence has pointed to bone fluid flow (BFF) as a media between an external load and the bone cells, in which altered velocity and pressure can ultimately initiate the mechanotransduction and the remodeling process within the bone. Load-induced BFF can be altered by factors such as intramedullary pressure (ImP) and/or bone matrix strain, mediating bone adaptation. Previous studies have shown that BFF induced by ImP alone, with minimum bone strain, can initiate bone remodeling. However, identifying induced ImP dynamics and bone strain factor in vivo using a non-invasive method still remains challenging. To apply ImP as a means for alteration of BFF, it was hypothesized that non-invasive dynamic hydraulic stimulation (DHS) can induce local ImP with minimal bone strain to potentially elicit osteogenic adaptive responses via bone-muscle coupling. The goal of this study was to evaluate the immediate effects on local and distant ImP and strain in response to a range of loading frequencies using DHS. Simultaneous femoral and tibial ImP and bone strain values were measured in three 15-month-old female Sprague Dawley rats during DHS loading on the tibia with frequencies of 1Hz to 10Hz. DHS showed noticeable effects on ImP induction in the stimulated tibia in a nonlinear fashion in response to DHS over the range of loading frequencies, where they peaked at 2Hz. DHS at various loading frequencies generated minimal bone strain in the tibiae. Maximal bone strain measured at all loading frequencies was less than 8με. No detectable induction of ImP or bone strain was observed in the femur. This study suggested that oscillatory DHS may regulate the local fluid dynamics with minimal mechanical strain in the bone, which serves critically in bone adaptation. These results clearly implied DHS's potential as an effective, non-invasive intervention for osteopenia and osteoporosis treatments. © 2013. Published by Elsevier Inc. All rights reserved.

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.

Fluid dynamics of out of equilibrium boost invariant plasmas

NASA Astrophysics Data System (ADS)

Blaizot, Jean-Paul; Yan, Li

2018-05-01

By solving a simple kinetic equation, in the relaxation time approximation, and for a particular set of moments of the distribution function, we establish a set of equations which, on the one hand, capture exactly the dynamics of the kinetic equation, and, on the other hand, coincide with the hierarchy of equations of viscous hydrodynamics, to arbitrary order in the viscous corrections. This correspondence sheds light on the underlying mechanism responsible for the apparent success of hydrodynamics in regimes that are far from local equilibrium.

Dynamic subfilter-scale stress model for large-eddy simulations

NASA Astrophysics Data System (ADS)

Rouhi, A.; Piomelli, U.; Geurts, B. J.

2016-08-01

We present a modification of the integral length-scale approximation (ILSA) model originally proposed by Piomelli et al. [Piomelli et al., J. Fluid Mech. 766, 499 (2015), 10.1017/jfm.2015.29] and apply it to plane channel flow and a backward-facing step. In the ILSA models the length scale is expressed in terms of the integral length scale of turbulence and is determined by the flow characteristics, decoupled from the simulation grid. In the original formulation the model coefficient was constant, determined by requiring a desired global contribution of the unresolved subfilter scales (SFSs) to the dissipation rate, known as SFS activity; its value was found by a set of coarse-grid calculations. Here we develop two modifications. We de-fine a measure of SFS activity (based on turbulent stresses), which adds to the robustness of the model, particularly at high Reynolds numbers, and removes the need for the prior coarse-grid calculations: The model coefficient can be computed dynamically and adapt to large-scale unsteadiness. Furthermore, the desired level of SFS activity is now enforced locally (and not integrated over the entire volume, as in the original model), providing better control over model activity and also improving the near-wall behavior of the model. Application of the local ILSA to channel flow and a backward-facing step and comparison with the original ILSA and with the dynamic model of Germano et al. [Germano et al., Phys. Fluids A 3, 1760 (1991), 10.1063/1.857955] show better control over the model contribution in the local ILSA, while the positive properties of the original formulation (including its higher accuracy compared to the dynamic model on coarse grids) are maintained. The backward-facing step also highlights the advantage of the decoupling of the model length scale from the mesh.

Computational fluid dynamics (CFD) insights into agitation stress methods in biopharmaceutical development.

PubMed

Bai, Ge; Bee, Jared S; Biddlecombe, James G; Chen, Quanmin; Leach, W Thomas

2012-02-28

Agitation of small amounts of liquid is performed routinely in biopharmaceutical process, formulation, and packaging development. Protein degradation commonly results from agitation, but the specific stress responsible or degradation mechanism is usually not well understood. Characterization of the agitation stress methods is critical to identifying protein degradation mechanisms or specific sensitivities. In this study, computational fluid dynamics (CFD) was used to model agitation of 1 mL of fluid by four types of common laboratory agitation instruments, including a rotator, orbital shaker, magnetic stirrer and vortex mixer. Fluid stresses in the bulk liquid and near interfaces were identified, quantified and compared. The vortex mixer provides the most intense stresses overall, while the stir bar system presented locally intense shear proximal to the hydrophobic stir bar surface. The rotator provides gentler fluid stresses, but the air-water interfacial area and surface stresses are relatively high given its low rotational frequency. The orbital shaker provides intermediate-level stresses but with the advantage of a large stable platform for consistent vial-to-vial homogeneity. Selection of experimental agitation methods with targeted types and intensities of stresses can facilitate better understanding of protein degradation mechanisms and predictability for "real world" applications. Copyright © 2011 Elsevier B.V. All rights reserved.

Mingus Discontinuous Multiphysics

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

Pat Notz, Dan Turner

Mingus provides hybrid coupled local/non-local mechanics analysis capabilities that extend several traditional methods to applications with inherent discontinuities. Its primary features include adaptations of solid mechanics, fluid dynamics and digital image correlation that naturally accommodate dijointed data or irregular solution fields by assimilating a variety of discretizations (such as control volume finite elements, peridynamics and meshless control point clouds). The goal of this software is to provide an analysis framework form multiphysics engineering problems with an integrated image correlation capability that can be used for experimental validation and model

Local ice melting by an antifreeze protein.

PubMed

Calvaresi, Matteo; Höfinger, Siegfried; Zerbetto, Francesco

2012-07-09

Antifreeze proteins, AFP, impede freezing of bodily fluids and damaging of cellular tissues by low temperatures. Adsorption-inhibition mechanisms have been developed to explain their functioning. Using in silico Molecular Dynamics, we show that type I AFP can also induce melting of the local ice surface. Simulations of antifreeze-positive and antifreeze-negative mutants show a clear correlation between melting induction and antifreeze activity. The presence of local melting adds a function to type I AFPs that is unique to these proteins. It may also explain some apparently conflicting experimental results where binding to ice appears both quasipermanent and reversible.

An EQT-cDFT approach to determine thermodynamic properties of confined fluids.

PubMed

Mashayak, S Y; Motevaselian, M H; Aluru, N R

2015-06-28

We present a continuum-based approach to predict the structure and thermodynamic properties of confined fluids at multiple length-scales, ranging from a few angstroms to macro-meters. The continuum approach is based on the empirical potential-based quasi-continuum theory (EQT) and classical density functional theory (cDFT). EQT is a simple and fast approach to predict inhomogeneous density and potential profiles of confined fluids. We use EQT potentials to construct a grand potential functional for cDFT. The EQT-cDFT-based grand potential can be used to predict various thermodynamic properties of confined fluids. In this work, we demonstrate the EQT-cDFT approach by simulating Lennard-Jones fluids, namely, methane and argon, confined inside slit-like channels of graphene. We show that the EQT-cDFT can accurately predict the structure and thermodynamic properties, such as density profiles, adsorption, local pressure tensor, surface tension, and solvation force, of confined fluids as compared to the molecular dynamics simulation results.

Criteria for Analyzing the Quality Estimate shake using two different impellers, “Rushton Turbine” and “PBT 45”

NASA Astrophysics Data System (ADS)

Oliveira, E. C. M.; Silva, M. O. M.; Nunhez, J. R.; Ribeiro, e. L. C. L. J.

2016-08-01

The Computational Fluid Dynamics (CFD) is used to numerically simulate the behavior of fluid flow, and all the laws that rule the study of fluids, both the energy and mass transfer, chemical reactions, hydraulic behavior, and other applications. Among these innumerable applications, the local behavior of a mixing tank can be estimated using the CFD techniques. In order to follow-up studies in the literature giving them increased credibility, this project aimed to determine whether both the scheduling Potency/Volume and Torque/Volume maintain the necessary level of agitation in the mix tank to extrapolate the volume 10 times and 100 times, using an axial impeller blades inclined at 45 degrees type (PBT 45), and a radial impeller (Rushton turbine) with the use of the CFD. With the results, was analyzed if a criterion is superior to the other in terms of maintaining the level of agitation in scheduling. The case study made it possible to achieve the main objective, deepen studies in Computational Fluid Dynamics. Therefore, managed to the largest domain in tools offered by CFX-ANSYS 14.0 software, allowing future work to be conducted.

Inertial effects during irreversible meniscus reconfiguration in angular pores

NASA Astrophysics Data System (ADS)

Ferrari, Andrea; Lunati, Ivan

2014-12-01

In porous media, the dynamics of the invading front between two immiscible fluids is often characterized by abrupt reconfigurations caused by local instabilities of the interface. As a prototype of these phenomena we consider the dynamics of a meniscus in a corner as it can be encountered in angular pores. We investigate this process in detail by means of direct numerical simulations that solve the Navier-Stokes equations in the pore space and employ the Volume of Fluid method (VOF) to track the evolution of the interface. We show that for a quasi-static displacement, the numerically calculated surface energy agrees well with the analytical solutions that we have derived for pores with circular and square cross sections. However, the spontaneous reconfigurations are irreversible and cannot be controlled by the injection rate: they are characterized by the amount of surface energy that is spontaneously released and transformed into kinetic energy. The resulting local velocities can be orders of magnitude larger than the injection velocity and they induce damped oscillations of the interface that possess their own time scales and depend only on fluid properties and pore geometry. In complex media (we consider a network of cubic pores) reconfigurations are so frequent and oscillations last long enough that increasing inertial effects leads to a different fluid distribution by influencing the selection of the next pore to be invaded. This calls into question simple pore-filling rules based only on capillary forces. Also, we demonstrate that inertial effects during irreversible reconfigurations can influence the work done by the external forces that is related to the pressure drop in Darcy's law. This suggests that these phenomena have to be considered when upscaling multiphase flow because local oscillations of the menisci affect macroscopic quantities and modify the constitutive relationships to be used in macro-scale models. These results can be extrapolated to other interface instabilities that are at the origin of fast pore-scale events, such as Haines jumps, snap-off and coalescence.

Studying dynamic processes in liquids by TEM/STEM/DTEM

NASA Astrophysics Data System (ADS)

Abellan, Patricia; Evans, James; Woehl, Taylor; Jungjohann, Katherine; Parent, Lucas; Arslan, Ilke; Ristenpart, William; Browning, Nigel; Mater. Sci. Group Team; Microsc. Group Team; Catal. Sci. Group Collaboration; Ristenpart Res. Group Collaboration

2013-03-01

In order to study dynamic phenomena such as corrosion or catalysis, extreme environmental conditions must be reproduced around the specimen - these include high-temperatures, high-pressures, specific oxidizing/reducing atmospheres or a liquid environment. The use of environmental stages specifically designed to fit in any transmission electron microscope (TEM) allows us to apply the distinct capabilities of each instrument to study dynamic processes. Localized gas/fluid conditions are created around the sample and separated from the high vacuum inside the microscope using hermetically sealed windowed-cells. Advanced capabilities of these techniques include spatial resolutions of ~1 Angstrom or better in aberration corrected instruments or temporal resolutions in the microsecond-nanosecond range in a dynamic TEM (DTEM). Here, unique qualities of the DTEM that benefit the in-situ experiments with gas/fluid environmental cells will be discussed. We also present our results with a liquid stage allowing atomic resolution imaging of nanomaterials in a colloidal suspension, core EEL spectra acquisition, continuous flow, controlled growth of nanocrystals and systematic calibration of the effect of the electron dose on silver nuclei formation.

Absorption and folding of melittin onto lipid bilayer membranes via unbiased atomic detail microsecond molecular dynamics simulation.

PubMed

Chen, Charles H; Wiedman, Gregory; Khan, Ayesha; Ulmschneider, Martin B

2014-09-01

Unbiased molecular simulation is a powerful tool to study the atomic details driving functional structural changes or folding pathways of highly fluid systems, which present great challenges experimentally. Here we apply unbiased long-timescale molecular dynamics simulation to study the ab initio folding and partitioning of melittin, a template amphiphilic membrane active peptide. The simulations reveal that the peptide binds strongly to the lipid bilayer in an unstructured configuration. Interfacial folding results in a localized bilayer deformation. Akin to purely hydrophobic transmembrane segments the surface bound native helical conformer is highly resistant against thermal denaturation. Circular dichroism spectroscopy experiments confirm the strong binding and thermostability of the peptide. The study highlights the utility of molecular dynamics simulations for studying transient mechanisms in fluid lipid bilayer systems. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova. Copyright © 2014. Published by Elsevier B.V.

Local support against gravity in magnetoturbulent fluids

NASA Astrophysics Data System (ADS)

Schmidt, W.; Collins, D. C.; Kritsuk, A. G.

2013-06-01

Comparisons of the integrated thermal pressure support of gas against its gravitational potential energy lead to critical mass scales for gravitational instability such as the Jeans and the Bonnor-Ebert masses, which play an important role in the analysis of many physical systems, including the heuristics of numerical simulations. In a strict theoretical sense, however, neither the Jeans nor the Bonnor-Ebert mass is meaningful when applied locally to substructure in a self-gravitating turbulent medium. For this reason, we investigate the local support by thermal pressure, turbulence and magnetic fields against gravitational compression through an approach that is independent of these concepts. At the centre of our approach is the dynamical equation for the divergence of the velocity field. We carry out a statistical analysis of the source terms of the local compression rate (the negative time derivative of the divergence) for simulations of forced self-gravitating turbulence in periodic boxes with zero, weak and moderately strong mean magnetic fields (measured by the averages of the magnetic and thermal pressures). We also consider the amplification of the magnetic field energy by shear and by compression. Thereby, we are able to demonstrate that the support against gravity is dominated by thermal pressure fluctuations, although magnetic pressure also yields a significant contribution. The net effect of turbulence in the highly supersonic regime, however, is to enhance compression rather than supporting overdense gas even if the vorticity is very high. This is incommensurate with the support of the highly dynamical substructures in magnetoturbulent fluids being determined by local virial equilibria of volume energies without surface stresses.

Modelling cavitation erosion using fluid–material interaction simulations

PubMed Central

Chahine, Georges L.; Hsiao, Chao-Tsung

2015-01-01

Material deformation and pitting from cavitation bubble collapse is investigated using fluid and material dynamics and their interaction. In the fluid, a novel hybrid approach, which links a boundary element method and a compressible finite difference method, is used to capture non-spherical bubble dynamics and resulting liquid pressures efficiently and accurately. The bubble dynamics is intimately coupled with a finite-element structure model to enable fluid/structure interaction simulations. Bubble collapse loads the material with high impulsive pressures, which result from shock waves and bubble re-entrant jet direct impact on the material surface. The shock wave loading can be from the re-entrant jet impact on the opposite side of the bubble, the fast primary collapse of the bubble, and/or the collapse of the remaining bubble ring. This produces high stress waves, which propagate inside the material, cause deformation, and eventually failure. A permanent deformation or pit is formed when the local equivalent stresses exceed the material yield stress. The pressure loading depends on bubble dynamics parameters such as the size of the bubble at its maximum volume, the bubble standoff distance from the material wall and the pressure driving the bubble collapse. The effects of standoff and material type on the pressure loading and resulting pit formation are highlighted and the effects of bubble interaction on pressure loading and material deformation are preliminarily discussed. PMID:26442140

B-spline Method in Fluid Dynamics

NASA Technical Reports Server (NTRS)

Botella, Olivier; Shariff, Karim; Mansour, Nagi N. (Technical Monitor)

2001-01-01

B-spline functions are bases for piecewise polynomials that possess attractive properties for complex flow simulations : they have compact support, provide a straightforward handling of boundary conditions and grid nonuniformities, and yield numerical schemes with high resolving power, where the order of accuracy is a mere input parameter. This paper reviews the progress made on the development and application of B-spline numerical methods to computational fluid dynamics problems. Basic B-spline approximation properties is investigated, and their relationship with conventional numerical methods is reviewed. Some fundamental developments towards efficient complex geometry spline methods are covered, such as local interpolation methods, fast solution algorithms on cartesian grid, non-conformal block-structured discretization, formulation of spline bases of higher continuity over triangulation, and treatment of pressure oscillations in Navier-Stokes equations. Application of some of these techniques to the computation of viscous incompressible flows is presented.

Flow Mapping in a Gas-Solid Riser via Computer Automated Radioactive Particle Tracking (CARPT)

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

Muthanna Al-Dahhan; Milorad P. Dudukovic; Satish Bhusarapu

2005-06-04

Statement of the Problem: Developing and disseminating a general and experimentally validated model for turbulent multiphase fluid dynamics suitable for engineering design purposes in industrial scale applications of riser reactors and pneumatic conveying, require collecting reliable data on solids trajectories, velocities ? averaged and instantaneous, solids holdup distribution and solids fluxes in the riser as a function of operating conditions. Such data are currently not available on the same system. Multiphase Fluid Dynamics Research Consortium (MFDRC) was established to address these issues on a chosen example of circulating fluidized bed (CFB) reactor, which is widely used in petroleum and chemicalmore » industry including coal combustion. This project addresses the problem of lacking reliable data to advance CFB technology. Project Objectives: The objective of this project is to advance the understanding of the solids flow pattern and mixing in a well-developed flow region of a gas-solid riser, operated at different gas flow rates and solids loading using the state-of-the-art non-intrusive measurements. This work creates an insight and reliable database for local solids fluid-dynamic quantities in a pilot-plant scale CFB, which can then be used to validate/develop phenomenological models for the riser. This study also attempts to provide benchmark data for validation of Computational Fluid Dynamic (CFD) codes and their current closures. Technical Approach: Non-Invasive Computer Automated Radioactive Particle Tracking (CARPT) technique provides complete Eulerian solids flow field (time average velocity map and various turbulence parameters such as the Reynolds stresses, turbulent kinetic energy, and eddy diffusivities). It also gives directly the Lagrangian information of solids flow and yields the true solids residence time distribution (RTD). Another radiation based technique, Computed Tomography (CT) yields detailed time averaged local holdup profiles at various planes. Together, these two techniques can provide the needed local solids flow dynamic information for the same setup under identical operating conditions, and the data obtained can be used as a benchmark for development, and refinement of the appropriate riser models. For the above reasons these two techniques were implemented in this study on a fully developed section of the riser. To derive the global mixing information in the riser, accurate solids RTD is needed and was obtained by monitoring the entry and exit of a single radioactive tracer. Other global parameters such as Cycle Time Distribution (CTD), overall solids holdup in the riser, solids recycle percentage at the bottom section of the riser were evaluated from different solids travel time distributions. Besides, to measure accurately and in-situ the overall solids mass flux, a novel method was applied.« less

Comparison of the effects of antihypertensive treatment with angiotensin II blockade and beta-blockade on carotid wall structure and haemodynamics: protocol and baseline demographics.

PubMed

Ariff, Ben; Stanton, Alice; Barratt, Dean; Augst, Alex; Glor, Fadi; Poulter, Neil; Sever, Peter; Xu, Yun; Hughes, Alun; Thom, Simon A Mc G

2002-06-01

Several systemic factors have been shown to contribute to the acceleration of large vessel atheroma. Correction of these factors leads to a reduction in the progression of plaque formation and associated arterial wall thickness. Atheroma remains, however, a focal disease, developing at characteristic sites within the arterial tree. These sites are typically at areas of vessel branching or marked curvature, and correspond to regions of high tensile stress and low sheer stress, leading to the hypothesis that local haemodynamic factors and vessel wall mechanics potentiate the focal development of atheroma. Current assessment of vascular haemodynamics suffers from an inability to handle complex flow, and does not allow accurate determination of locally varying flow, and shear stress patterns. The application of computational fluid dynamic (CFD) flow simulation techniques to ultrasound and local pressure data, however, allows a comprehensive, non-invasive appraisal of haemodynamic flow parameters to be performed. The Candesartan cilexetil and Atenolol Carotid Haemodynamic Endpoint Trial (CACHET) study compares the effects of two antihypertensive regimens, one b-blocker-based, the other angiotensin receptor blocker based, on carotid intima-media thickness. The collection of ultrasound and pressure data on each subject provides a unique opportunity to apply these data to the CFD model to study the effects of these antihypertensive regimens on local fluid dynamics. This will lead to a greater understanding of the relationship of these factors to atheroma formation and regression.

Bifurcations and degenerate periodic points in a three dimensional chaotic fluid flow

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

Smith, L. D., E-mail: lachlan.smith@monash.edu; CSIRO Mineral Resources, Clayton, Victoria 3800; Rudman, M.

2016-05-15

Analysis of the periodic points of a conservative periodic dynamical system uncovers the basic kinematic structure of the transport dynamics and identifies regions of local stability or chaos. While elliptic and hyperbolic points typically govern such behaviour in 3D systems, degenerate (parabolic) points also play an important role. These points represent a bifurcation in local stability and Lagrangian topology. In this study, we consider the ramifications of the two types of degenerate periodic points that occur in a model 3D fluid flow. (1) Period-tripling bifurcations occur when the local rotation angle associated with elliptic points is reversed, creating a reversalmore » in the orientation of associated Lagrangian structures. Even though a single unstable point is created, the bifurcation in local stability has a large influence on local transport and the global arrangement of manifolds as the unstable degenerate point has three stable and three unstable directions, similar to hyperbolic points, and occurs at the intersection of three hyperbolic periodic lines. The presence of period-tripling bifurcation points indicates regions of both chaos and confinement, with the extent of each depending on the nature of the associated manifold intersections. (2) The second type of bifurcation occurs when periodic lines become tangent to local or global invariant surfaces. This bifurcation creates both saddle–centre bifurcations which can create both chaotic and stable regions, and period-doubling bifurcations which are a common route to chaos in 2D systems. We provide conditions for the occurrence of these tangent bifurcations in 3D conservative systems, as well as constraints on the possible types of tangent bifurcation that can occur based on topological considerations.« less

Coalescence and Interaction of Solitons in the Coupled Korteweg-de Vries System

NASA Astrophysics Data System (ADS)

Chung, Wai Choi; Chow, Kwok Wing

2017-11-01

There are many physical systems which are governed by the classical Korteweg-de Vries equation. One of the prominent examples is the shallow water wave in fluid dynamics. In recent years, a coupled Korteweg-de Vries system has been proposed to describe fluids in a two-layer flow, and coherent structures in terms of solitons are found. We studied the coupled Korteweg-de Vries system by means of the Hirota bilinear method. Soliton and breather solutions are constructed. Localized pulses which result from the coupling of waves can be formed. The structure of the localized pulses becomes asymmetric as the control parameter varies. The coalescence and interaction of solitons in the coupled Korteweg-de Vries system will be discussed. Partial financial support has been provided by the Research Grants Council contract HKU 17200815.

Learning to classify wakes from local sensory information

NASA Astrophysics Data System (ADS)

Alsalman, Mohamad; Colvert, Brendan; Kanso, Eva; Kanso Team

2017-11-01

Aquatic organisms exhibit remarkable abilities to sense local flow signals contained in their fluid environment and to surmise the origins of these flows. For example, fish can discern the information contained in various flow structures and utilize this information for obstacle avoidance and prey tracking. Flow structures created by flapping and swimming bodies are well characterized in the fluid dynamics literature; however, such characterization relies on classical methods that use an external observer to reconstruct global flow fields. The reconstructed flows, or wakes, are then classified according to the unsteady vortex patterns. Here, we propose a new approach for wake identification: we classify the wakes resulting from a flapping airfoil by applying machine learning algorithms to local flow information. In particular, we simulate the wakes of an oscillating airfoil in an incoming flow, extract the downstream vorticity information, and train a classifier to learn the different flow structures and classify new ones. This data-driven approach provides a promising framework for underwater navigation and detection in application to autonomous bio-inspired vehicles.

Electroosmotic flow and Joule heating in preparative continuous annular electrochromatography.

PubMed

Laskowski, René; Bart, Hans-Jörg

2015-09-01

An openFOAM "computational fluid dynamic" simulation model was developed for the description of local interaction of hydrodynamics and Joule heating in annular electrochromatography. A local decline of electrical conductivity of the background eluent is caused by an electrokinetic migration of ions resulting in higher Joule heat generation. The model equations consider the Navier-Stokes equation for incompressible fluids, the energy equation for stationary temperature fields, and the mass transfer equation for the electrokinetic flow. The simulations were embedded in commercial ANSYS Fluent software and in open-source environment openFOAM. The annular gap (1 mm width) contained an inorganic C8 reverse-phase monolith as stationary phase prepared by an in situ sol-gel process. The process temperature generated by Joule heating was determined by thermal camera system. The local hydrodynamics in the prototype was detected by a gravimetric contact-free measurement method and experimental and simulated values matched quite well. © 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Apparatus for characterizing the temporo-spatial properties of a dynamic fluid front and method thereof

DOEpatents

Battiste, Richard L.

2007-12-25

Methods and apparatus are described for characterizing the temporal-spatial properties of a dynamic fluid front within a mold space while the mold space is being filled with fluid. A method includes providing a mold defining a mold space and having one or more openings into the mold space; heating a plurality of temperature sensors that extend into the mold space; injecting a fluid into the mold space through the openings, the fluid experiencing a dynamic fluid front while filling the mold space with the fluid; and characterizing temporal-spatial properties of the dynamic fluid front by monitoring a temperature of each of the plurality of heated temperature sensors while the mold space is being filled with the fluid. An apparatus includes a mold defining a mold space; one or more openings for introducing a fluid into the mold space and filling the mold space with the fluid, the fluid experiencing a dynamic fluid front while filling the mold space; a plurality of heated temperature sensors extending into the mold space; and a computer coupled to the plurality of heated temperature sensors for characterizing the temporal-spatial properties of the dynamic fluid front.

Apparatus for characterizing the temporo-spatial properties of a dynamic fluid front and method thereof

DOEpatents

Battiste, Richard L

2013-12-31

Methods and apparatus are described for characterizing the temporal-spatial properties of a dynamic fluid front within a mold space while the mold space is being filled with fluid. A method includes providing a mold defining a mold space and having one or more openings into the mold space; heating a plurality of temperature sensors that extend into the mold space; injecting a fluid into th emold space through the openings, the fluid experiencing a dynamic fluid front while filling the mold space with a fluid; and characterizing temporal-spatial properties of the dynamic fluid front by monitoring a termperature of each of the plurality of heated temperature sensors while the mold space is being filled with the fluid. An apparatus includes a mold defining a mold space; one or more openings for introducing a fluid into th emold space and filling the mold space with the fluid, the fluid experiencing a dynamic fluid front while filling the mold space; a plurality of heated temperature sensors extending into the mold space; and a computer coupled to the plurality of heated temperature sensors for characterizing the temporal-spatial properties of the dynamic fluid front.

Endocytic reawakening of motility in jammed epithelia

NASA Astrophysics Data System (ADS)

Malinverno, Chiara; Corallino, Salvatore; Giavazzi, Fabio; Bergert, Martin; Li, Qingsen; Leoni, Marco; Disanza, Andrea; Frittoli, Emanuela; Oldani, Amanda; Martini, Emanuele; Lendenmann, Tobias; Deflorian, Gianluca; Beznoussenko, Galina V.; Poulikakos, Dimos; Ong, Kok Haur; Uroz, Marina; Trepat, Xavier; Parazzoli, Dario; Maiuri, Paolo; Yu, Weimiao; Ferrari, Aldo; Cerbino, Roberto; Scita, Giorgio

2017-05-01

Dynamics of epithelial monolayers has recently been interpreted in terms of a jamming or rigidity transition. How cells control such phase transitions is, however, unknown. Here we show that RAB5A, a key endocytic protein, is sufficient to induce large-scale, coordinated motility over tens of cells, and ballistic motion in otherwise kinetically arrested monolayers. This is linked to increased traction forces and to the extension of cell protrusions, which align with local velocity. Molecularly, impairing endocytosis, macropinocytosis or increasing fluid efflux abrogates RAB5A-induced collective motility. A simple model based on mechanical junctional tension and an active cell reorientation mechanism for the velocity of self-propelled cells identifies regimes of monolayer dynamics that explain endocytic reawakening of locomotion in terms of a combination of large-scale directed migration and local unjamming. These changes in multicellular dynamics enable collectives to migrate under physical constraints and may be exploited by tumours for interstitial dissemination.

Complex dynamics induced by strong confinement - From tracer diffusion in strongly heterogeneous media to glassy relaxation of dense fluids in narrow slits

NASA Astrophysics Data System (ADS)

Mandal, Suvendu; Spanner-Denzer, Markus; Leitmann, Sebastian; Franosch, Thomas

2017-08-01

We provide an overview of recent advances of the complex dynamics of particles in strong confinements. The first paradigm is the Lorentz model where tracers explore a quenched disordered host structure. Such systems naturally occur as limiting cases of binary glass-forming systems if the dynamics of one component is much faster than the other. For a certain critical density of the host structure the tracers undergo a localization transition which constitutes a critical phenomenon. A series of predictions in the vicinity of the transition have been elaborated and tested versus computer simulations. Analytical progress is achieved for small obstacle densities. The second paradigm is a dense strongly interacting liquid confined to a narrow slab. Then the glass transition depends nonmonotonically on the separation of the plates due to an interplay of local packing and layering. Very small slab widths allow to address certain features of the statics and dynamics analytically.

High-throughput microfluidics to control and measure signaling dynamics in single yeast cells

PubMed Central

Hansen, Anders S.; Hao, Nan; O'Shea, Erin K.

2015-01-01

Microfluidics coupled to quantitative time-lapse fluorescence microscopy is transforming our ability to control, measure, and understand signaling dynamics in single living cells. Here we describe a pipeline that incorporates multiplexed microfluidic cell culture, automated programmable fluid handling for cell perturbation, quantitative time-lapse microscopy, and computational analysis of time-lapse movies. We illustrate how this setup can be used to control the nuclear localization of the budding yeast transcription factor Msn2. Using this protocol, we generate oscillations of Msn2 localization and measure the dynamic gene expression response of individual genes in single cells. The protocol allows a single researcher to perform up to 20 different experiments in a single day, whilst collecting data for thousands of single cells. Compared to other protocols, the present protocol is relatively easy to adopt and higher-throughput. The protocol can be widely used to control and monitor single-cell signaling dynamics in other signal transduction systems in microorganisms. PMID:26158443

Role of deep crustal fluids in the genesis of intraplate earthquakes in the Kachchh region, northwestern India

NASA Astrophysics Data System (ADS)

Pavan Kumar, G.; Mahesh, P.; Nagar, Mehul; Mahender, E.; Kumar, Virendhar; Mohan, Kapil; Ravi Kumar, M.

2017-05-01

Fluids play a prominent role in the genesis of earthquakes, particularly in intraplate settings. In this study, we present evidence for a highly heterogeneous nature of electrical conductivity in the crust and uppermost mantle beneath the Kachchh rift basin of northwestern India, which is host to large, deadly intraplate earthquakes. We interpret our results of high conductive zones inferred from magnetotelluric and 3-D local earthquake tomography investigations in terms of a fluid reservoir in the upper mantle. The South Wagad Fault (SWF) imaged as a near-vertical north dipping low resistivity zone traversing the entire crust and an elongated south dipping conductor demarcating the North Wagad Fault (NWF) serve as conduits for fluid flow from the reservoir to the middle to lower crustal depths. Importantly, the epicentral zone of the 2001 main shock is characterized as a fluid saturated zone at the rooting of NWF onto the SWF.