Lagrangian statistics and flow topology in forced two-dimensional turbulence.

PubMed

Kadoch, B; Del-Castillo-Negrete, D; Bos, W J T; Schneider, K

2011-03-01

A study of the relationship between Lagrangian statistics and flow topology in fluid turbulence is presented. The topology is characterized using the Weiss criterion, which provides a conceptually simple tool to partition the flow into topologically different regions: elliptic (vortex dominated), hyperbolic (deformation dominated), and intermediate (turbulent background). The flow corresponds to forced two-dimensional Navier-Stokes turbulence in doubly periodic and circular bounded domains, the latter with no-slip boundary conditions. In the double periodic domain, the probability density function (pdf) of the Weiss field exhibits a negative skewness consistent with the fact that in periodic domains the flow is dominated by coherent vortex structures. On the other hand, in the circular domain, the elliptic and hyperbolic regions seem to be statistically similar. We follow a Lagrangian approach and obtain the statistics by tracking large ensembles of passively advected tracers. The pdfs of residence time in the topologically different regions are computed introducing the Lagrangian Weiss field, i.e., the Weiss field computed along the particles' trajectories. In elliptic and hyperbolic regions, the pdfs of the residence time have self-similar algebraic decaying tails. In contrast, in the intermediate regions the pdf has exponential decaying tails. The conditional pdfs (with respect to the flow topology) of the Lagrangian velocity exhibit Gaussian-like behavior in the periodic and in the bounded domains. In contrast to the freely decaying turbulence case, the conditional pdfs of the Lagrangian acceleration in forced turbulence show a comparable level of intermittency in both the periodic and the bounded domains. The conditional pdfs of the Lagrangian curvature are characterized, in all cases, by self-similar power-law behavior with a decay exponent of order -2.

Lagrangian chaos in three- dimensional steady buoyancy-driven flows

NASA Astrophysics Data System (ADS)

Contreras, Sebastian; Speetjens, Michel; Clercx, Herman

2016-11-01

Natural convection plays a key role in fluid dynamics owing to its ubiquitous presence in nature and industry. Buoyancy-driven flows are prototypical systems in the study of thermal instabilities and pattern formation. The differentially heated cavity problem has been widely studied for the investigation of buoyancy-induced oscillatory flow. However, far less attention has been devoted to the three-dimensional Lagrangian transport properties in such flows. This study seeks to address this by investigating Lagrangian transport in the steady flow inside a cubic cavity differentially-heated from the side. The theoretical and numerical analysis expands on previously reported similarities between the current flow and lid-driven flows. The Lagrangian dynamics are controlled by the Péclet number (Pe) and the Prandtl number (Pr). Pe controls the behaviour qualitatively in that growing Pe progressively perturbs the integable state (Pe =0), thus paving the way to chaotic dynamics. Pr plays an entirely quantitative role in that Pr<1 and Pr>1 amplifies and diminishes, respectively, the perturbative effect of non-zero Pe. S.C. acknowledges financial support from Consejo Nacional de Ciencia y Tecnología (CONACYT).

Lagrangian transport in a class of three-dimensional buoyancy-driven flows

NASA Astrophysics Data System (ADS)

Contreras, Sebastian; Speetjens, Michel; Clercx, Herman

2017-11-01

The study concerns the Lagrangian dynamics of three-dimensional (3D) buoyancy-driven cavity flows under steady and laminar conditions due to a global temperature gradient imposed via an opposite hot and cold sidewall. This serves as archetypal configuration for natural-convection flows in which gravity is perpendicular to the global temperature gradient. Limited insight into the Lagrangian properties of this class of flows motivates this study. The 3D Lagrangian dynamics are investigated in terms of the generic structure of the Lagrangian flow topology that is described in terms of the Grashof number (Gr) and the Prandtl number (Pr). Gr is the principal control parameter for the flow topology: vanishing Gr yields a state of closed streamlines (integrable state); increasing Gr causes the formation of toroidal coherent structures embedded in chaotic streamlines governed by Hamiltonian mechanisms. Fluid inertia prevails for ``smaller'' Gr. A buoyancy-induced bifurcation of the flow topology occurs for ``larger'' Gr and underlies the emergence of ``secondary rolls'' and secondary tori for ``larger'' Pr. Stagnation points and corresponding manifold interactions are key to the dynamics. S.C. acknowledges financial support from Consejo Nacional de Ciencia y Tecnología (CONACYT).

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

NASA Technical Reports Server (NTRS)

Felici, Helene Marie

1992-01-01

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

A Galerkin least squares approach to viscoelastic flow.

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

Rao, Rekha R.; Schunk, Peter Randall

2015-10-01

A Galerkin/least-squares stabilization technique is applied to a discrete Elastic Viscous Stress Splitting formulation of for viscoelastic flow. From this, a possible viscoelastic stabilization method is proposed. This method is tested with the flow of an Oldroyd-B fluid past a rigid cylinder, where it is found to produce inaccurate drag coefficients. Furthermore, it fails for relatively low Weissenberg number indicating it is not suited for use as a general algorithm. In addition, a decoupled approach is used as a way separating the constitutive equation from the rest of the system. A Pressure Poisson equation is used when the velocity andmore » pressure are sought to be decoupled, but this fails to produce a solution when inflow/outflow boundaries are considered. However, a coupled pressure-velocity equation with a decoupled constitutive equation is successful for the flow past a rigid cylinder and seems to be suitable as a general-use algorithm.« less

Quantitative flow analysis of swimming dynamics with coherent Lagrangian vortices.

PubMed

Huhn, F; van Rees, W M; Gazzola, M; Rossinelli, D; Haller, G; Koumoutsakos, P

2015-08-01

Undulatory swimmers flex their bodies to displace water, and in turn, the flow feeds back into the dynamics of the swimmer. At moderate Reynolds number, the resulting flow structures are characterized by unsteady separation and alternating vortices in the wake. We use the flow field from simulations of a two-dimensional, incompressible viscous flow of an undulatory, self-propelled swimmer and detect the coherent Lagrangian vortices in the wake to dissect the driving momentum transfer mechanisms. The detected material vortex boundary encloses a Lagrangian control volume that serves to track back the vortex fluid and record its circulation and momentum history. We consider two swimming modes: the C-start escape and steady anguilliform swimming. The backward advection of the coherent Lagrangian vortices elucidates the geometry of the vorticity field and allows for monitoring the gain and decay of circulation and momentum transfer in the flow field. For steady swimming, momentum oscillations of the fish can largely be attributed to the momentum exchange with the vortex fluid. For the C-start, an additionally defined jet fluid region turns out to balance the high momentum change of the fish during the rapid start.

Hamiltonian stability for weighted measure and generalized Lagrangian mean curvature flow

NASA Astrophysics Data System (ADS)

Kajigaya, Toru; Kunikawa, Keita

2018-06-01

In this paper, we generalize several results for the Hamiltonian stability and the mean curvature flow of Lagrangian submanifolds in a Kähler-Einstein manifold to more general Kähler manifolds including a Fano manifold equipped with a Kähler form ω ∈ 2 πc1(M) by using the method proposed by Behrndt (2011). Namely, we first consider a weighted measure on a Lagrangian submanifold L in a Kähler manifold M and investigate the variational problem of L for the weighted volume functional. We call a stationary point of the weighted volume functional f-minimal, and define the notion of Hamiltonian f-stability as a local minimizer under Hamiltonian deformations. We show such examples naturally appear in a toric Fano manifold. Moreover, we consider the generalized Lagrangian mean curvature flow in a Fano manifold which is introduced by Behrndt and Smoczyk-Wang. We generalize the result of H. Li, and show that if the initial Lagrangian submanifold is a small Hamiltonian deformation of an f-minimal and Hamiltonian f-stable Lagrangian submanifold, then the generalized MCF converges exponentially fast to an f-minimal Lagrangian submanifold.

A purely Lagrangian method for computing linearly-perturbed flows in spherical geometry

NASA Astrophysics Data System (ADS)

Jaouen, Stéphane

2007-07-01

In many physical applications, one wishes to control the development of multi-dimensional instabilities around a one-dimensional (1D) complex flow. For predicting the growth rates of these perturbations, a general numerical approach is viable which consists in solving simultaneously the one-dimensional equations and their linearized form for three-dimensional perturbations. In Clarisse et al. [J.-M. Clarisse, S. Jaouen, P.-A. Raviart, A Godunov-type method in Lagrangian coordinates for computing linearly-perturbed planar-symmetric flows of gas dynamics, J. Comp. Phys. 198 (2004) 80-105], a class of Godunov-type schemes for planar-symmetric flows of gas dynamics has been proposed. Pursuing this effort, we extend these results to spherically symmetric flows. A new method to derive the Lagrangian perturbation equations, based on the canonical form of systems of conservation laws with zero entropy flux [B. Després, Lagrangian systems of conservation laws. Invariance properties of Lagrangian systems of conservation laws, approximate Riemann solvers and the entropy condition, Numer. Math. 89 (2001) 99-134; B. Després, C. Mazeran, Lagrangian gas dynamics in two dimensions and Lagrangian systems, Arch. Rational Mech. Anal. 178 (2005) 327-372] is also described. It leads to many advantages. First of all, many physical problems we are interested in enter this formalism (gas dynamics, two-temperature plasma equations, ideal magnetohydrodynamics, etc.) whatever is the geometry. Secondly, a class of numerical entropic schemes is available for the basic flow [11]. Last, linearizing and devising numerical schemes for the perturbed flow is straightforward. The numerical capabilities of these methods are illustrated on three test cases of increasing difficulties and we show that - due to its simplicity and its low computational cost - the Linear Perturbations Code (LPC) is a powerful tool to understand and predict the development of hydrodynamic instabilities in the linear regime.

Microfluidic co-flow of Newtonian and viscoelastic fluids for high-resolution separation of microparticles.

PubMed

Tian, Fei; Zhang, Wei; Cai, Lili; Li, Shanshan; Hu, Guoqing; Cong, Yulong; Liu, Chao; Li, Tiejun; Sun, Jiashu

2017-09-12

The microfluidic passive control of microparticles largely relies on the hydrodynamic effects of the carrier media such as Newtonian fluids and viscoelastic fluids. Yet the viscoelastic/Newtonian interfacial effect has been scarcely investigated, especially for high-resolution particle separation. Here we report a microfluidic co-flow of Newtonian (water or PBS) and viscoelastic fluids (PEO) for the size-dependent separation of microparticles. The co-flow condition generates a stable viscoelastic/Newtonian interface, giving rise to the wall-directed elastic lift forces that compete with the center-directed lift forces, and efficiently hinders the migration of microparticles from the Newtonian to the viscoelastic fluid in a size-dependent manner. An almost complete separation of a binary mixture of 1 μm and 2 μm polystyrene particles is achieved by the co-flow of water and a very dilute PEO solution (100 ppm), whereas the sole use of water or PEO could not lead to an efficient separation. This co-flow microfluidic system is also applied for the separation of Staphylococcus aureus (1 μm) from platelets (2-3 μm) with >90% efficiencies and purities.

Advances in the analysis and prediction of turbulent viscoelastic flows

NASA Astrophysics Data System (ADS)

Gatski, T. B.; Thais, L.; Mompean, G.

2014-08-01

It has been well-known for over six decades that the addition of minute amounts of long polymer chains to organic solvents, or water, can lead to significant turbulent drag reduction. This discovery has had many practical applications such as in pipeline fluid transport, oil well operations, vehicle design and submersible vehicle projectiles, and more recently arteriosclerosis treatment. However, it has only been the last twenty-five years that the full utilization of direct numerical simulation of such turbulent viscoelastic flows has been achieved. The unique characteristics of viscoelastic fluid flow are dictated by the nonlinear differential relationship between the flow strain rate field and the extra-stress induced by the additive polymer. A primary motivation for the analysis of these turbulent fluid flows is the understanding of the effect on the dynamic transfer of energy in the turbulent flow due to the presence of the extra-stress field induced by the presence of the viscoelastic polymer chain. Such analyses now utilize direct numerical simulation data of fully developed channel flow for the FENE-P (Finite Extendable Nonlinear Elastic - Peterlin) fluid model. Such multi-scale dynamics suggests an analysis of the transfer of energy between the various component motions that include the turbulent kinetic energy, and the mean polymeric and elastic potential energies. It is shown that the primary effect of the interaction between the turbulent and polymeric fields is to transfer energy from the turbulence to the polymer.

Stochastic Lagrangian dynamics for charged flows in the E-F regions of ionosphere

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

Tang Wenbo; Mahalov, Alex

2013-03-15

We develop a three-dimensional numerical model for the E-F region ionosphere and study the Lagrangian dynamics for plasma flows in this region. Our interest rests on the charge-neutral interactions and the statistics associated with stochastic Lagrangian motion. In particular, we examine the organizing mixing patterns for plasma flows due to polarized gravity wave excitations in the neutral field, using Lagrangian coherent structures (LCS). LCS objectively depict the flow topology-the extracted attractors indicate generation of ionospheric density gradients, due to accumulation of plasma. Using Lagrangian measures such as the finite-time Lyapunov exponents, we locate the Lagrangian skeletons for mixing in plasma,more » hence where charged fronts are expected to appear. With polarized neutral wind, we find that the corresponding plasma velocity is also polarized. Moreover, the polarized velocity alone, coupled with stochastic Lagrangian motion, may give rise to polarized density fronts in plasma. Statistics of these trajectories indicate high level of non-Gaussianity. This includes clear signatures of variance, skewness, and kurtosis of displacements taking polarized structures aligned with the gravity waves, and being anisotropic.« less

Lagrangian flows within reflecting internal waves at a horizontal free-slip surface

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

Zhou, Qi, E-mail: q.zhou@damtp.cam.ac.uk; Diamessis, Peter J.

In this paper sequel to Zhou and Diamessis [“Reflection of an internal gravity wave beam off a horizontal free-slip surface,” Phys. Fluids 25, 036601 (2013)], we consider Lagrangian flows within nonlinear internal waves (IWs) reflecting off a horizontal free-slip rigid lid, the latter being a model of the ocean surface. The problem is approached both analytically using small-amplitude approximations and numerically by tracking Lagrangian fluid particles in direct numerical simulation (DNS) datasets of the Eulerian flow. Inviscid small-amplitude analyses for both plane IWs and IW beams (IWBs) show that Eulerian mean flow due to wave-wave interaction and wave-induced Stokes driftmore » cancels each other out completely at the second order in wave steepness A, i.e., O(A{sup 2}), implying zero Lagrangian mean flow up to that order. However, high-accuracy particle tracking in finite-Reynolds-number fully nonlinear DNS datasets from the work of Zhou and Diamessis suggests that the Euler-Stokes cancelation on O(A{sup 2}) is not complete. This partial cancelation significantly weakens the mean Lagrangian flows but does not entirely eliminate them. As a result, reflecting nonlinear IWBs produce mean Lagrangian drifts on O(A{sup 2}) and thus particle dispersion on O(A{sup 4}). The above findings can be relevant to predicting IW-driven mass transport in the oceanic surface and subsurface region which bears important observational and environmental implications, under circumstances where the effect of Earth rotation can be ignored.« less

Viscoelastic effects on residual oil distribution in flows through pillared microchannels.

PubMed

De, S; Krishnan, P; van der Schaaf, J; Kuipers, J A M; Peters, E A J F; Padding, J T

2018-01-15

Multiphase flow through porous media is important in a number of industrial, natural and biological processes. One application is enhanced oil recovery (EOR), where a resident oil phase is displaced by a Newtonian or polymeric fluid. In EOR, the two-phase immiscible displacement through heterogonous porous media is usually governed by competing viscous and capillary forces, expressed through a Capillary number Ca, and viscosity ratio of the displacing and displaced fluid. However, when viscoelastic displacement fluids are used, elastic forces in the displacement fluid also become significant. It is hypothesized that elastic instabilities are responsible for enhanced oil recovery through an elastic microsweep mechanism. In this work, we use a simplified geometry in the form of a pillared microchannel. We analyze the trapped residual oil size distribution after displacement by a Newtonian fluid, a nearly inelastic shear thinning fluid, and viscoelastic polymers and surfactant solutions. We find that viscoelastic polymers and surfactant solutions can displace more oil compared to Newtonian fluids and nearly inelastic shear thinning polymers at similar Ca numbers. Beyond a critical Ca number, the size of residual oil blobs decreases significantly for viscoelastic fluids. This critical Ca number directly corresponds to flow rates where elastic instabilities occur in single phase flow, suggesting a close link between enhancement of oil recovery and appearance of elastic instabilities. Copyright © 2017 Elsevier Inc. All rights reserved.

Strong polymer-turbulence interactions in viscoelastic turbulent channel flow.

PubMed

Dallas, V; Vassilicos, J C; Hewitt, G F

2010-12-01

This paper is focused on the fundamental mechanism(s) of viscoelastic turbulence that leads to polymer-induced turbulent drag reduction phenomenon. A great challenge in this problem is the computation of viscoelastic turbulent flows, since the understanding of polymer physics is restricted to mechanical models. An effective state-of-the-art numerical method to solve the governing equation for polymers modeled as nonlinear springs, without using any artificial assumptions as usual, was implemented here on a three-dimensional channel flow geometry. The capability of this algorithm to capture the strong polymer-turbulence dynamical interactions is depicted on the results, which are much closer qualitatively to experimental observations. This allowed a more detailed study of the polymer-turbulence interactions, which yields an enhanced picture on a mechanism resulting from the polymer-turbulence energy transfers.

Assimilating Eulerian and Lagrangian data in traffic-flow models

NASA Astrophysics Data System (ADS)

Xia, Chao; Cochrane, Courtney; DeGuire, Joseph; Fan, Gaoyang; Holmes, Emma; McGuirl, Melissa; Murphy, Patrick; Palmer, Jenna; Carter, Paul; Slivinski, Laura; Sandstede, Björn

2017-05-01

Data assimilation of traffic flow remains a challenging problem. One difficulty is that data come from different sources ranging from stationary sensors and camera data to GPS and cell phone data from moving cars. Sensors and cameras give information about traffic density, while GPS data provide information about the positions and velocities of individual cars. Previous methods for assimilating Lagrangian data collected from individual cars relied on specific properties of the underlying computational model or its reformulation in Lagrangian coordinates. These approaches make it hard to assimilate both Eulerian density and Lagrangian positional data simultaneously. In this paper, we propose an alternative approach that allows us to assimilate both Eulerian and Lagrangian data. We show that the proposed algorithm is accurate and works well in different traffic scenarios and regardless of whether ensemble Kalman or particle filters are used. We also show that the algorithm is capable of estimating parameters and assimilating real traffic observations and synthetic observations obtained from microscopic models.

Uncertainty quantification in Eulerian-Lagrangian models for particle-laden flows

NASA Astrophysics Data System (ADS)

Fountoulakis, Vasileios; Jacobs, Gustaaf; Udaykumar, Hs

2017-11-01

A common approach to ameliorate the computational burden in simulations of particle-laden flows is to use a point-particle based Eulerian-Lagrangian model, which traces individual particles in their Lagrangian frame and models particles as mathematical points. The particle motion is determined by Stokes drag law, which is empirically corrected for Reynolds number, Mach number and other parameters. The empirical corrections are subject to uncertainty. Treating them as random variables renders the coupled system of PDEs and ODEs stochastic. An approach to quantify the propagation of this parametric uncertainty to the particle solution variables is proposed. The approach is based on averaging of the governing equations and allows for estimation of the first moments of the quantities of interest. We demonstrate the feasibility of our proposed methodology of uncertainty quantification of particle-laden flows on one-dimensional linear and nonlinear Eulerian-Lagrangian systems. This research is supported by AFOSR under Grant FA9550-16-1-0008.

Lagrangian analysis of multiscale particulate flows with the particle finite element method

NASA Astrophysics Data System (ADS)

Oñate, Eugenio; Celigueta, Miguel Angel; Latorre, Salvador; Casas, Guillermo; Rossi, Riccardo; Rojek, Jerzy

2014-05-01

We present a Lagrangian numerical technique for the analysis of flows incorporating physical particles of different sizes. The numerical approach is based on the particle finite element method (PFEM) which blends concepts from particle-based techniques and the FEM. The basis of the Lagrangian formulation for particulate flows and the procedure for modelling the motion of small and large particles that are submerged in the fluid are described in detail. The numerical technique for analysis of this type of multiscale particulate flows using a stabilized mixed velocity-pressure formulation and the PFEM is also presented. Examples of application of the PFEM to several particulate flows problems are given.

Viscoelastic flow in rotating curved pipes

NASA Astrophysics Data System (ADS)

Chen, Yitung; Chen, Huajun; Zhang, Jinsuo; Zhang, Benzhao

2006-08-01

Fully developed viscoelastic flows in rotating curved pipes with circular cross section are investigated theoretically and numerically employing the Oldroyd-B fluid model. Based on Dean's approximation, a perturbation solution up to the secondary order is obtained. The governing equations are also solved numerically by the finite volume method. The theoretical and numerical solutions agree with each other very well. The results indicate that the rotation, as well as the curvature and elasticity, plays an important role in affecting the friction factor, the secondary flow pattern and intensity. The co-rotation enhances effects of curvature and elasticity on the secondary flow. For the counter-rotation, there is a critical rotational number RΩ', which can make the effect of rotation counteract the effect of curvature and elasticity. Complicated flow behaviors are found at this value. For the relative creeping flow, RΩ' can be estimated according to the expression RΩ'=-4Weδ. Effects of curvature and elasticity at different rotational numbers on both relative creeping flow and inertial flow are also analyzed and discussed.

Investigation of particle lateral migration in sample-sheath flow of viscoelastic fluid and Newtonian fluid.

PubMed

Yuan, Dan; Zhang, Jun; Yan, Sheng; Peng, Gangrou; Zhao, Qianbin; Alici, Gursel; Du, Hejun; Li, Weihua

2016-08-01

In this work, particle lateral migration in sample-sheath flow of viscoelastic fluid and Newtonian fluid was experimentally investigated. The 4.8-μm micro-particles were dispersed in a polyethylene oxide (PEO) viscoelastic solution, and then the solution was injected into a straight rectangular channel with a deionised (DI) water Newtonian sheath flow. Micro-particles suspended in PEO solution migrated laterally to a DI water stream, but migration in the opposite direction from a DI water stream to a PEO solution stream or from one DI water stream to another DI water stream could not be achieved. The lateral migration of particles depends on the viscoelastic properties of the sample fluids. Furthermore, the effects of channel length, flow rate, and PEO concentration were studied. By using viscoelastic sample flow and Newtonian sheath flow, a selective particle lateral migration can be achieved in a simple straight channel, without any external force fields. This particle lateral migration technique could be potentially used in solution exchange fields such as automated cell staining and washing in microfluidic platforms, and holds numerous biomedical applications. © 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Chaos analysis of viscoelastic chaotic flows of polymeric fluids in a micro-channel

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

Lim, C. P.; Lam, Y. C., E-mail: myclam@ntu.edu.sg; BioSystems and Micromechanics

2015-07-15

Many fluids, including biological fluids such as mucus and blood, are viscoelastic. Through the introduction of chaotic flows in a micro-channel and the construction of maps of characteristic chaos parameters, differences in viscoelastic properties of these fluids can be measured. This is demonstrated by creating viscoelastic chaotic flows induced in an H-shaped micro-channel through the steady infusion of a polymeric fluid of polyethylene oxide (PEO) and another immiscible fluid (silicone oil). A protocol for chaos analysis was established and demonstrated for the analysis of the chaotic flows generated by two polymeric fluids of different molecular weight but with similar relaxationmore » times. The flows were shown to be chaotic through the computation of their correlation dimension (D{sub 2}) and the largest Lyapunov exponent (λ{sub 1}), with D{sub 2} being fractional and λ{sub 1} being positive. Contour maps of D{sub 2} and λ{sub 1} of the respective fluids in the operating space, which is defined by the combination of polymeric fluids and silicone oil flow rates, were constructed to represent the characteristic of the chaotic flows generated. It was observed that, albeit being similar, the fluids have generally distinct characteristic maps with some similar trends. The differences in the D{sub 2} and λ{sub 1} maps are indicative of the difference in the molecular weight of the polymers in the fluids because the driving force of the viscoelastic chaotic flows is of molecular origin. This approach in constructing the characteristic maps of chaos parameters can be employed as a diagnostic tool for biological fluids and, more generally, chaotic signals.« less

A new Lagrangian method for three-dimensional steady supersonic flows

NASA Technical Reports Server (NTRS)

Loh, Ching-Yuen; Liou, Meng-Sing

1993-01-01

In this report, the new Lagrangian method introduced by Loh and Hui is extended for three-dimensional, steady supersonic flow computation. The derivation of the conservation form and the solution of the local Riemann solver using the Godunov and the high-resolution TVD (total variation diminished) scheme is presented. This new approach is accurate and robust, capable of handling complicated geometry and interactions between discontinuous waves. Test problems show that the extended Lagrangian method retains all the advantages of the two-dimensional method (e.g., crisp resolution of a slip-surface (contact discontinuity) and automatic grid generation). In this report, we also suggest a novel three dimensional Riemann problem in which interesting and intricate flow features are present.

Reduction of numerical diffusion in three-dimensional vortical flows using a coupled Eulerian/Lagrangian solution procedure

NASA Technical Reports Server (NTRS)

Felici, Helene M.; Drela, Mark

1993-01-01

A new approach based on the coupling of an Eulerian and a Lagrangian solver, aimed at reducing the numerical diffusion errors of standard Eulerian time-marching finite-volume solvers, is presented. The approach is applied to the computation of the secondary flow in two bent pipes and the flow around a 3D wing. Using convective point markers the Lagrangian approach provides a correction of the basic Eulerian solution. The Eulerian flow in turn integrates in time the Lagrangian state-vector. A comparison of coarse and fine grid Eulerian solutions makes it possible to identify numerical diffusion. It is shown that the Eulerian/Lagrangian approach is an effective method for reducing numerical diffusion errors.

Analytical solution of two-fluid electro-osmotic flows of viscoelastic fluids.

PubMed

Afonso, A M; Alves, M A; Pinho, F T

2013-04-01

This paper presents an analytical model that describes a two-fluid electro-osmotic flow of stratified fluids with Newtonian or viscoelastic rheological behavior. This is the principle of operation of an electro-osmotic two-fluid pump as proposed by Brask et al. [Tech. Proc. Nanotech., 1, 190-193, 2003], in which an electrically non-conducting fluid is transported by the interfacial dragging viscous force of a conducting fluid that is driven by electro-osmosis. The electric potential in the conducting fluid and the analytical steady flow solution of the two-fluid electro-osmotic stratified flow in a planar microchannel are presented by assuming a planar interface between the two immiscible fluids with Newtonian or viscoelastic rheological behavior. The effects of fluid rheology, shear viscosity ratio, holdup and interfacial zeta potential are analyzed to show the viability of this technique, where an enhancement of the flow rate is observed as the shear-thinning effects are increased. Copyright © 2012 Elsevier Inc. All rights reserved.

Lagrangian transport properties of pulmonary interfacial flows

PubMed Central

Smith, Bradford J.; Lukens, Sarah; Yamaguchi, Eiichiro; Gaver, Donald P.

2012-01-01

Disease states characterized by airway fluid occlusion and pulmonary surfactant insufficiency, such as respiratory distress syndrome, have a high mortality rate. Understanding the mechanics of airway reopening, particularly involving surfactant transport, may provide an avenue to increase patient survival via optimized mechanical ventilation waveforms. We model the occluded airway as a liquid-filled rigid tube with the fluid phase displaced by a finger of air that propagates with both mean and sinusoidal velocity components. Finite-time Lyapunov exponent (FTLE) fields are employed to analyse the convective transport characteristics, taking note of Lagrangian coherent structures (LCSs) and their effects on transport. The Lagrangian perspective of these techniques reveals flow characteristics that are not readily apparent by observing Eulerian measures. These analysis techniques are applied to surfactant-free velocity fields determined computationally, with the boundary element method, and measured experimentally with micro particle image velocimetry (μ-PIV). We find that the LCS divides the fluid into two regimes, one advected upstream (into the thin residual film) and the other downstream ahead of the advancing bubble. At higher oscillatory frequencies particles originating immediately inside the LCS experience long residence times at the air–liquid interface, which may be conducive to surfactant transport. At high frequencies a well-mixed attractor region is identified; this volume of fluid cyclically travels along the interface and into the bulk fluid. The Lagrangian analysis is applied to velocity data measured with 0.01 mg ml−1 of the clinical pulmonary surfactant Infasurf in the bulk fluid, demonstrating flow field modifications with respect to the surfactant-free system that were not visible in the Eulerian frame. PMID:23049141

A viscoelastic fluid-structure interaction model for carotid arteries under pulsatile flow.

PubMed

Wang, Zhongjie; Wood, Nigel B; Xu, Xiao Yun

2015-05-01

In this study, a fluid-structure interaction model (FSI) incorporating viscoelastic wall behaviour is developed and applied to an idealized model of the carotid artery under pulsatile flow. The shear and bulk moduli of the arterial wall are described by Prony series, where the parameters can be derived from in vivo measurements. The aim is to develop a fully coupled FSI model that can be applied to realistic arterial geometries with normal or pathological viscoelastic wall behaviour. Comparisons between the numerical and analytical solutions for wall displacements demonstrate that the coupled model is capable of predicting the viscoelastic behaviour of carotid arteries. Comparisons are also made between the solid only and FSI viscoelastic models, and the results suggest that the difference in radial displacement between the two models is negligible. Copyright © 2015 John Wiley & Sons, Ltd.

Analytical and numerical study of the electro-osmotic annular flow of viscoelastic fluids.

PubMed

Ferrás, L L; Afonso, A M; Alves, M A; Nóbrega, J M; Pinho, F T

2014-04-15

In this work we present semi-analytical solutions for the electro-osmotic annular flow of viscoelastic fluids modeled by the Linear and Exponential PTT models. The viscoelastic fluid flows in the axial direction between two concentric cylinders under the combined influences of electrokinetic and pressure forcings. The analysis invokes the Debye-Hückel approximation and includes the limit case of pure electro-osmotic flow. The solution is valid for both no slip and slip velocity at the walls and the chosen slip boundary condition is the linear Navier slip velocity model. The combined effects of fluid rheology, electro-osmotic and pressure gradient forcings on the fluid velocity distribution are also discussed. Copyright © 2013 Elsevier Inc. All rights reserved.

Scalable Methods for Eulerian-Lagrangian Simulation Applied to Compressible Multiphase Flows

NASA Astrophysics Data System (ADS)

Zwick, David; Hackl, Jason; Balachandar, S.

2017-11-01

Multiphase flows can be found in countless areas of physics and engineering. Many of these flows can be classified as dispersed two-phase flows, meaning that there are solid particles dispersed in a continuous fluid phase. A common technique for simulating such flow is the Eulerian-Lagrangian method. While useful, this method can suffer from scaling issues on larger problem sizes that are typical of many realistic geometries. Here we present scalable techniques for Eulerian-Lagrangian simulations and apply it to the simulation of a particle bed subjected to expansion waves in a shock tube. The results show that the methods presented here are viable for simulation of larger problems on modern supercomputers. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE-1315138. This work was supported in part by the U.S. Department of Energy under Contract No. DE-NA0002378.

Lagrangian transport near perturbed periodic lines in three-dimensional unsteady flows

NASA Astrophysics Data System (ADS)

Speetjens, Michel

2015-11-01

Periodic lines formed by continuous strings of periodic points are key organizing entities in the Lagrangian flow topology of certain three-dimensional (3D) time-periodic flows. Such lines generically consist of elliptic and/or hyperbolic points and thus give rise to 3D flow topologies made up of families of concentric closed trajectories embedded in chaotic regions. Weak perturbation destroys the periodic lines and causes said trajectories to coalesce into families of concentric tubes. However, emergence of isolated periodic points near the disintegrating periodic lines and/or partitioning of the original lines into elliptic and hyperbolic segments interrupt the tube formation. This yields incomplete tubes that interact with the (chaotic) environment through their open ends, resulting in intricate and essentially 3D flow topologies These phenomena have been observed in various realistic flows yet the underlying mechanisms are to date only partially understood. This study deepens insight into the (perturbed) Lagrangian dynamics of these flows by way of a linearized representation of the equations of motion near the periodic lines. Predictions on the basis of this investigation are in full (qualitative) agreement with observed behavior in the actual flows

Viscoelastic flow modeling in the extrusion of a dough-like fluid

NASA Technical Reports Server (NTRS)

Dhanasekharan, M.; Kokini, J. L.; Janes, H. W. (Principal Investigator)

2000-01-01

This work attempts to investigate the effect of viscoelasticity and three-dimensional geometry in screw channels. The Phan-Thien Tanner (PTT) constitutive equation with simplified model parameters was solved in conjunction with the flow equations. Polyflow, a commercially available finite element code was used to solve the resulting nonlinear partial differential equations. The PTT model predicted one log scale lower pressure buildup compared to the equivalent Newtonian results. However, the velocity profile did not show significant changes for the chosen PTT model parameters. Past Researchers neglected viscoelastic effects and also the three dimensional nature of the flow in extruder channels. The results of this paper provide a starting point for further simulations using more realistic model parameters, which may enable the food engineer to more accurately scale-up and design extrusion processes.

Modeling and Numerical Challenges in Eulerian-Lagrangian Computations of Shock-driven Multiphase Flows

NASA Astrophysics Data System (ADS)

Diggs, Angela; Balachandar, Sivaramakrishnan

2015-06-01

The present work addresses the numerical methods required for particle-gas and particle-particle interactions in Eulerian-Lagrangian simulations of multiphase flow. Local volume fraction as seen by each particle is the quantity of foremost importance in modeling and evaluating such interactions. We consider a general multiphase flow with a distribution of particles inside a fluid flow discretized on an Eulerian grid. Particle volume fraction is needed both as a Lagrangian quantity associated with each particle and also as an Eulerian quantity associated with the flow. In Eulerian Projection (EP) methods, the volume fraction is first obtained within each cell as an Eulerian quantity and then interpolated to each particle. In Lagrangian Projection (LP) methods, the particle volume fraction is obtained at each particle and then projected onto the Eulerian grid. Traditionally, EP methods are used in multiphase flow, but sub-grid resolution can be obtained through use of LP methods. By evaluating the total error and its components we compare the performance of EP and LP methods. The standard von Neumann error analysis technique has been adapted for rigorous evaluation of rate of convergence. The methods presented can be extended to obtain accurate field representations of other Lagrangian quantities. Most importantly, we will show that such careful attention to numerical methodologies is needed in order to capture complex shock interaction with a bed of particles. Supported by U.S. Department of Defense SMART Program and the U.S. Department of Energy PSAAP-II program under Contract No. DE-NA0002378.

Lagrangian transported MDF methods for compressible high speed flows

NASA Astrophysics Data System (ADS)

Gerlinger, Peter

2017-06-01

This paper deals with the application of thermochemical Lagrangian MDF (mass density function) methods for compressible sub- and supersonic RANS (Reynolds Averaged Navier-Stokes) simulations. A new approach to treat molecular transport is presented. This technique on the one hand ensures numerical stability of the particle solver in laminar regions of the flow field (e.g. in the viscous sublayer) and on the other hand takes differential diffusion into account. It is shown in a detailed analysis, that the new method correctly predicts first and second-order moments on the basis of conventional modeling approaches. Moreover, a number of challenges for MDF particle methods in high speed flows is discussed, e.g. high cell aspect ratio grids close to solid walls, wall heat transfer, shock resolution, and problems from statistical noise which may cause artificial shock systems in supersonic flows. A Mach 2 supersonic mixing channel with multiple shock reflection and a model rocket combustor simulation demonstrate the eligibility of this technique to practical applications. Both test cases are simulated successfully for the first time with a hybrid finite-volume (FV)/Lagrangian particle solver (PS).

Growth of viscoelastic wings and the reduction of particle mobility in a viscoelastic shear flow

NASA Astrophysics Data System (ADS)

Murch, William L.; Krishnan, Sreenath; Shaqfeh, Eric S. G.; Iaccarino, Gianluca

2017-10-01

The motion of a rigid spherical particle in a sheared polymeric fluid is studied via experiments and numerical simulations. We study particle mobility in highly elastic fluids, where the deformation due to the sphere's movement and the shear flow both result in significant stretching of the polymer. The shear flow is imposed in a plane perpendicular to the sphere's movement, resulting in regions of high polymer tension in the wake of the sphere that can extend well into the shear flow and gradient directions. We observe that these viscoelastic wake structures, resembling wings, are linked to an increase in the form drag, providing a mechanism for a dramatic decrease in the particle mobility.

Unsteady Boundary-Layer Flow over Jerked Plate Moving in a Free Stream of Viscoelastic Fluid

PubMed Central

Mehmood, Ahmer; Ali, Asif; Saleem, Najma

2014-01-01

This study aims to investigate the unsteady boundary-layer flow of a viscoelastic non-Newtonian fluid over a flat surface. The plate is suddenly jerked to move with uniform velocity in a uniform stream of non-Newtonian fluid. Purely analytic solution to governing nonlinear equation is obtained. The solution is highly accurate and valid for all values of the dimensionless time 0 ≤ τ < ∞. Flow properties of the viscoelastic fluid are discussed through graphs. PMID:24892060

Extreme Lagrangian acceleration in confined turbulent flow.

PubMed

Kadoch, Benjamin; Bos, Wouter J T; Schneider, Kai

2008-05-09

A Lagrangian study of two-dimensional turbulence for two different geometries, a periodic and a confined circular geometry, is presented to investigate the influence of solid boundaries on the Lagrangian dynamics. It is found that the Lagrangian acceleration is even more intermittent in the confined domain than in the periodic domain. The flatness of the Lagrangian acceleration as a function of the radius shows that the influence of the wall on the Lagrangian dynamics becomes negligible in the center of the domain, and it also reveals that the wall is responsible for the increased intermittency. The transition in the Lagrangian statistics between this region, not directly influenced by the walls, and a critical radius which defines a Lagrangian boundary layer is shown to be very sharp with a sudden increase of the acceleration flatness from about 5 to about 20.

Spatio-temporal organization of dynamics in a two-dimensional periodically driven vortex flow: A Lagrangian flow network perspective.

PubMed

Lindner, Michael; Donner, Reik V

2017-03-01

We study the Lagrangian dynamics of passive tracers in a simple model of a driven two-dimensional vortex resembling real-world geophysical flow patterns. Using a discrete approximation of the system's transfer operator, we construct a directed network that describes the exchange of mass between distinct regions of the flow domain. By studying different measures characterizing flow network connectivity at different time-scales, we are able to identify the location of dynamically invariant structures and regions of maximum dispersion. Specifically, our approach allows us to delimit co-existing flow regimes with different dynamics. To validate our findings, we compare several network characteristics to the well-established finite-time Lyapunov exponents and apply a receiver operating characteristic analysis to identify network measures that are particularly useful for unveiling the skeleton of Lagrangian chaos.

V-ONSET: Introducing turbulent multiphase flow facility focusing on Lagrangian interfacial transfer dynamics

NASA Astrophysics Data System (ADS)

Salibindla, Ashwanth; Masuk, Ashik Ullah Mohammad; Ni, Rui

2017-11-01

We have designed and constructed a new vertical water tunnel, V-ONSET, to investigate interfacial mass, momentum and energy transfer between two phases in a Lagrangian frame. This system features an independent control of mean flow and turbulence level. The mean flow opposes the rising/falling velocity of the second phase, ``suspending'' the particles and increasing tracking time in the view area. Strong turbulence is generated by shooting 88 digitally-controlled water jets into the test section. The second phase, either bubbles or oil droplets, can be introduced into the test section through a capillary island. In addition to this flow control system, V-ONSET comes with a 3D two-phase visualization system, consisting of high-speed cameras, two-colored LED system, and in-house Lagrangian particle tracking algorithm. This enables us to acquire the Lagrangian evolution of both phases and the interfacial transfer dynamics in between, paving the way for new closure models for two-phase simulations. Financial support for this project was provided by National Science Foundation under Grant Number: 1653389 and 1705246.

Lagrangian velocity and acceleration correlations of large inertial particles in a closed turbulent flow

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

Machicoane, Nathanaël; Volk, Romain

We investigate the response of large inertial particle to turbulent fluctuations in an inhomogeneous and anisotropic flow. We conduct a Lagrangian study using particles both heavier and lighter than the surrounding fluid, and whose diameters are comparable to the flow integral scale. Both velocity and acceleration correlation functions are analyzed to compute the Lagrangian integral time and the acceleration time scale of such particles. The knowledge of how size and density affect these time scales is crucial in understanding particle dynamics and may permit stochastic process modelization using two-time models (for instance, Sawford’s). As particles are tracked over long timesmore » in the quasi-totality of a closed flow, the mean flow influences their behaviour and also biases the velocity time statistics, in particular the velocity correlation functions. By using a method that allows for the computation of turbulent velocity trajectories, we can obtain unbiased Lagrangian integral time. This is particularly useful in accessing the scale separation for such particles and to comparing it to the case of fluid particles in a similar configuration.« less

Three-Dimensional Mixed Convection Flow of Viscoelastic Fluid with Thermal Radiation and Convective Conditions

PubMed Central

Hayat, Tasawar; Ashraf, Muhammad Bilal; Alsulami, Hamed H.; Alhuthali, Muhammad Shahab

2014-01-01

The objective of present research is to examine the thermal radiation effect in three-dimensional mixed convection flow of viscoelastic fluid. The boundary layer analysis has been discussed for flow by an exponentially stretching surface with convective conditions. The resulting partial differential equations are reduced into a system of nonlinear ordinary differential equations using appropriate transformations. The series solutions are developed through a modern technique known as the homotopy analysis method. The convergent expressions of velocity components and temperature are derived. The solutions obtained are dependent on seven sundry parameters including the viscoelastic parameter, mixed convection parameter, ratio parameter, temperature exponent, Prandtl number, Biot number and radiation parameter. A systematic study is performed to analyze the impacts of these influential parameters on the velocity and temperature, the skin friction coefficients and the local Nusselt number. It is observed that mixed convection parameter in momentum and thermal boundary layers has opposite role. Thermal boundary layer is found to decrease when ratio parameter, Prandtl number and temperature exponent are increased. Local Nusselt number is increasing function of viscoelastic parameter and Biot number. Radiation parameter on the Nusselt number has opposite effects when compared with viscoelastic parameter. PMID:24608594

Three-dimensional mixed convection flow of viscoelastic fluid with thermal radiation and convective conditions.

PubMed

Hayat, Tasawar; Ashraf, Muhammad Bilal; Alsulami, Hamed H; Alhuthali, Muhammad Shahab

2014-01-01

The objective of present research is to examine the thermal radiation effect in three-dimensional mixed convection flow of viscoelastic fluid. The boundary layer analysis has been discussed for flow by an exponentially stretching surface with convective conditions. The resulting partial differential equations are reduced into a system of nonlinear ordinary differential equations using appropriate transformations. The series solutions are developed through a modern technique known as the homotopy analysis method. The convergent expressions of velocity components and temperature are derived. The solutions obtained are dependent on seven sundry parameters including the viscoelastic parameter, mixed convection parameter, ratio parameter, temperature exponent, Prandtl number, Biot number and radiation parameter. A systematic study is performed to analyze the impacts of these influential parameters on the velocity and temperature, the skin friction coefficients and the local Nusselt number. It is observed that mixed convection parameter in momentum and thermal boundary layers has opposite role. Thermal boundary layer is found to decrease when ratio parameter, Prandtl number and temperature exponent are increased. Local Nusselt number is increasing function of viscoelastic parameter and Biot number. Radiation parameter on the Nusselt number has opposite effects when compared with viscoelastic parameter.

Analysis of Lagrangian stretching in turbulent channel flow using a database task-parallel particle tracking approach

NASA Astrophysics Data System (ADS)

Meneveau, Charles; Johnson, Perry; Hamilton, Stephen; Burns, Randal

2016-11-01

An intrinsic property of turbulent flows is the exponential deformation of fluid elements along Lagrangian paths. The production of enstrophy by vorticity stretching follows from a similar mechanism in the Lagrangian view, though the alignment statistics differ and viscosity prevents unbounded growth. In this paper, the stretching properties of fluid elements and vorticity along Lagrangian paths are studied in a channel flow at Reτ = 1000 and compared with prior, known results from isotropic turbulence. To track Lagrangian paths in a public database containing Direct Numerical Simulation (DNS) results, the task-parallel approach previously employed in the isotropic database is extended to the case of flow in a bounded domain. It is shown that above 100 viscous units from the wall, stretching statistics are equal to their isotropic values, in support of the local isotropy hypothesis. Normalized by dissipation rate, the stretching in the buffer layer and below is less efficient due to less favorable alignment statistics. The Cramér function characterizing cumulative Lagrangian stretching statistics shows that overall the channel flow has about half of the stretching per unit dissipation compared with isotropic turbulence. Supported by a National Science Foundation Graduate Research Fellowship Program under Grant No. DGE-1232825, and by National Science Foundation Grants CBET-1507469, ACI-1261715, OCI-1244820 and by JHU IDIES.

Spatial-temporal dynamics of Newtonian and viscoelastic turbulence in channel flow

NASA Astrophysics Data System (ADS)

Wang, Sung-Ning; Shekar, Ashwin; Graham, Michael

2016-11-01

Introducing a trace amount of polymer into liquid turbulent flows can result in substantial reduction of friction drag. This phenomenon has been widely used in fluid transport; however, the mechanism is not well understood. Past studies have found that in minimal domain turbulent simulations, there areoccasional time periods when flow exhibits features such as weaker vortices, lower friction drag and larger log-law slope; these have been denoted as "hibernatingturbulence". Here we address the question of whether similar behavior arises spatio-temporally in extended domains, focusing on turbulence at friction Reynolds numbers near transition and Weissenberg numbers resulting in low-medium drag reduction. By using image analysis and conditional sampling tools, we identify the hibernating states in extended domains and show that they display striking similarity as those in minimal domains. The hibernating states among different Weissenberg numbers exhibit similar flow statistics, suggesting they are unaltered by low to medium viscoelasticity. In addition, the polymer is much less stretched during hibernation. Finally, these hibernating states vanish as Reynolds number increases. However, they reoccur and gradually become dominant with increasing viscoelasticity.

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

NASA Astrophysics Data System (ADS)

Xu, Xiaoyang; Deng, Xiao-Long

2016-04-01

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

a Marker-Based Eulerian-Lagrangian Method for Multiphase Flow with Supersonic Combustion Applications

NASA Astrophysics Data System (ADS)

Fan, Xiaofeng; Wang, Jiangfeng

2016-06-01

The atomization of liquid fuel is a kind of intricate dynamic process from continuous phase to discrete phase. Procedures of fuel spray in supersonic flow are modeled with an Eulerian-Lagrangian computational fluid dynamics methodology. The method combines two distinct techniques and develops an integrated numerical simulation method to simulate the atomization processes. The traditional finite volume method based on stationary (Eulerian) Cartesian grid is used to resolve the flow field, and multi-component Navier-Stokes equations are adopted in present work, with accounting for the mass exchange and heat transfer occupied by vaporization process. The marker-based moving (Lagrangian) grid is utilized to depict the behavior of atomized liquid sprays injected into a gaseous environment, and discrete droplet model 13 is adopted. To verify the current approach, the proposed method is applied to simulate processes of liquid atomization in supersonic cross flow. Three classic breakup models, TAB model, wave model and K-H/R-T hybrid model, are discussed. The numerical results are compared with multiple perspectives quantitatively, including spray penetration height and droplet size distribution. In addition, the complex flow field structures induced by the presence of liquid spray are illustrated and discussed. It is validated that the maker-based Eulerian-Lagrangian method is effective and reliable.

Lagrangian coherent structures separate dynamically distinct regions in fluid flows.

PubMed

Kelley, Douglas H; Allshouse, Michael R; Ouellette, Nicholas T

2013-07-01

Using filter-space techniques, we study the scale-to-scale transport of energy in a quasi-two-dimensional, weakly turbulent fluid flow averaged along the trajectories of fluid elements. We find that although the spatial mean of this Lagrangian-averaged flux is nearly unchanged from its Eulerian counterpart, the spatial structure of the scale-to-scale energy flux changes significantly. In particular, its features appear to correlate with the positions of Lagrangian coherent structures (LCS's). We show that the LCS's tend to lie at zeros of the scale-to-scale flux, and therefore that the LCS's separate regions that have qualitatively different dynamics. Since LCS's are also known to be impenetrable barriers to advection and mixing, we therefore find that the fluid on either side of an LCS is both kinematically and dynamically distinct. Our results extend the utility of LCS's by making clear the role they play in the flow dynamics in addition to the kinematics.

Viscoelastic flow in the lower crust after the 1992 landers, california, earthquake

PubMed

Deng; Gurnis; Kanamori; Hauksson

1998-11-27

Space geodesy showed that broad-scale postseismic deformation occurred after the 1992 Landers earthquake. Three-dimensional modeling shows that afterslip can only explain one horizontal component of the postseismic deformation, whereas viscoelastic flow can explain the horizontal and near-vertical displacements. The viscosity of a weak, about 10-km-thick layer, in the lower crust beneath the rupture zone that controls the rebound is about 10(18) pascal seconds. The viscoelastic behavior of the lower crust may help to explain the extensional structures observed in the Basin and Range province and it may be used for the analysis of earthquake hazard.

Simultaneous Measurements of Geometric and Viscoelastic Properties of Hydrogel Microbeads Using Continuous-Flow Microfluidics with Embedded Electrodes.

PubMed

Niu, Ye; Zhang, Xu; Si, Ting; Zhang, Yuntian; Qi, Lin; Zhao, Gang; Xu, Ronald X; He, Xiaoming; Zhao, Yi

2017-12-01

Geometric and mechanical characterizations of hydrogel materials at the microscale are attracting increasing attention due to their importance in tissue engineering, regenerative medicine, and drug delivery applications. Contemporary approaches for measuring the these properties of hydrogel microbeads suffer from low-throughput, complex system configuration, and measurement inaccuracy. In this work, a continuous-flow device is developed to measure geometric and viscoelastic properties of hydrogel microbeads by flowing the microbeads through a tapered microchannel with an array of interdigitated microelectrodes patterned underneath the channel. The viscoelastic properties are derived from the trajectories of microbeads using a quasi-linear viscoelastic model. The measurement is independent of the applied volumetric flow rate. The results show that the geometric and viscoelastic properties of Ca-alginate hydrogel microbeads can be determined independently and simultaneously. The bulky high-speed optical systems are eliminated, simplifying the system configuration and making it a truly miniaturized device. A throughput of up to 394 microbeads min -1 is achieved. This study may provide a powerful tool for mechanical profiling of hydrogel microbeads to support their wide applications. © 2017 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Lagrangian turbulence: Structures and mixing in admissible model flows

NASA Astrophysics Data System (ADS)

Ottino, Julio M.

1991-12-01

The goal of our research was to bridge the gap between modern ideas from dynamical systems and chaos and more traditional approaches to turbulence. In order to reach this objective we conducted theoretical and computational work on two systems: (1) a perturbed-Kelvin cat eyes flow, and (2) prototype solutions of the Navier-Stokes equations near solid walls. The main results obtained are two-fold: we have been able to produce flows capable of producing complex distributions of vorticity, and we have been able to construct flowfields, based on solutions of the Navier-Stokes equations, which are capable of displaying both Eulerian and Lagrangian turbulence. These results exemplify typical mechanisms of mixing enhancement in transitional flows.

The Effect of Pulsatile Versus Nonpulsatile Blood Flow on Viscoelasticity and Red Blood Cell Aggregation in Extracorporeal Circulation

PubMed Central

Ahn, Chi Bum; Kang, Yang Jun; Kim, Myoung Gon; Yang, Sung; Lim, Choon Hak; Son, Ho Sung; Kim, Ji Sung; Lee, So Young; Son, Kuk Hui; Sun, Kyung

2016-01-01

Background Extracorporeal circulation (ECC) can induce alterations in blood viscoelasticity and cause red blood cell (RBC) aggregation. In this study, the authors evaluated the effects of pump flow pulsatility on blood viscoelasticity and RBC aggregation. Methods Mongrel dogs were randomly assigned to two groups: a nonpulsatile pump group (n=6) or a pulsatile pump group (n=6). After ECC was started at a pump flow rate of 80 mL/kg/min, cardiac fibrillation was induced. Blood sampling was performed before and at 1, 2, and 3 hours after ECC commencement. To eliminate bias induced by hematocrit and plasma, all blood samples were adjusted to a hematocrit of 45% using baseline plasma. Blood viscoelasticity, plasma viscosity, hematocrit, arterial blood gas analysis, central venous O2 saturation, and lactate were measured. Results The blood viscosity and aggregation index decreased abruptly 1 hour after ECC and then remained low during ECC in both groups, but blood elasticity did not change during ECC. Blood viscosity, blood elasticity, plasma viscosity, and the aggregation index were not significantly different in the groups at any time. Hematocrit decreased abruptly 1 hour after ECC in both groups due to dilution by the priming solution used. Conclusion After ECC, blood viscoelasticity and RBC aggregation were not different in the pulsatile and nonpulsatile groups in the adult dog model. Furthermore, pulsatile flow did not have a more harmful effect on blood viscoelasticity or RBC aggregation than nonpulsatile flow. PMID:27298790

The Effect of Pulsatile Versus Nonpulsatile Blood Flow on Viscoelasticity and Red Blood Cell Aggregation in Extracorporeal Circulation.

PubMed

Ahn, Chi Bum; Kang, Yang Jun; Kim, Myoung Gon; Yang, Sung; Lim, Choon Hak; Son, Ho Sung; Kim, Ji Sung; Lee, So Young; Son, Kuk Hui; Sun, Kyung

2016-06-01

Extracorporeal circulation (ECC) can induce alterations in blood viscoelasticity and cause red blood cell (RBC) aggregation. In this study, the authors evaluated the effects of pump flow pulsatility on blood viscoelasticity and RBC aggregation. Mongrel dogs were randomly assigned to two groups: a nonpulsatile pump group (n=6) or a pulsatile pump group (n=6). After ECC was started at a pump flow rate of 80 mL/kg/min, cardiac fibrillation was induced. Blood sampling was performed before and at 1, 2, and 3 hours after ECC commencement. To eliminate bias induced by hematocrit and plasma, all blood samples were adjusted to a hematocrit of 45% using baseline plasma. Blood viscoelasticity, plasma viscosity, hematocrit, arterial blood gas analysis, central venous O2 saturation, and lactate were measured. The blood viscosity and aggregation index decreased abruptly 1 hour after ECC and then remained low during ECC in both groups, but blood elasticity did not change during ECC. Blood viscosity, blood elasticity, plasma viscosity, and the aggregation index were not significantly different in the groups at any time. Hematocrit decreased abruptly 1 hour after ECC in both groups due to dilution by the priming solution used. After ECC, blood viscoelasticity and RBC aggregation were not different in the pulsatile and nonpulsatile groups in the adult dog model. Furthermore, pulsatile flow did not have a more harmful effect on blood viscoelasticity or RBC aggregation than nonpulsatile flow.

Defect-mediated turbulence in ribbons of viscoelastic Taylor-Couette flow.

PubMed

Latrache, Noureddine; Abcha, Nizar; Crumeyrolle, Olivier; Mutabazi, Innocent

2016-04-01

Transition to defect-mediated turbulence in the ribbon patterns observed in a viscoelastic Taylor-Couette flow is investigated when the rotation rate of the inner cylinder is increased while the outer cylinder is fixed. In four polymer solutions with different values of the elasticity number, the defects appear just above the onset of the ribbon pattern and trigger the appearance of disordered oscillations when the rotation rate is increased. The flow structure around the defects is determined and the statistical properties of these defects are analyzed in the framework of the complex Ginzburg-Landau equation.

Grid and subgrid-scale interactions in viscoelastic turbulent flow and implications for modelling

NASA Astrophysics Data System (ADS)

Masoudian, M.; da Silva, C. B.; Pinho, F. T.

2016-06-01

Using direct numerical simulations of turbulent plane channel flow of homogeneous polymer solutions, described by the Finitely Extensible Nonlinear Elastic-Peterlin (FENE-P) rheological constitutive model, a-priori analyses of the filtered momentum and FENE-P constitutive equations are performed. The influence of the polymer additives on the subgrid-scale (SGS) energy is evaluated by comparing the Newtonian and the viscoelastic flows, and a severe suppression of SGS stresses and energy is observed in the viscoelastic flow. All the terms of the transport equation of the SGS kinetic energy for FENE-P fluids are analysed, and an approximated version of this equation for use in future large eddy simulation closures is suggested. The terms responsible for kinetic energy transfer between grid-scale (GS) and SGS energy (split into forward/backward energy transfer) are evaluated in the presence of polymers. It is observed that the probability and intensity of forward scatter events tend to decrease in the presence of polymers.

A new mixed subgrid-scale model for large eddy simulation of turbulent drag-reducing flows of viscoelastic fluids

NASA Astrophysics Data System (ADS)

Li, Feng-Chen; Wang, Lu; Cai, Wei-Hua

2015-07-01

A mixed subgrid-scale (SGS) model based on coherent structures and temporal approximate deconvolution (MCT) is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation (LES) of turbulent drag-reducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence (FHIT) with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation (DNS) results. Compared with the LES results using the temporal approximate deconvolution model (TADM) for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number. For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives. Project supported by the China Postdoctoral Science Foundation (Grant No. 2011M500652), the National Natural Science Foundation of China (Grant Nos. 51276046 and 51206033), and the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20112302110020).

Lagrangian mass-flow investigations of inorganic contaminants in wastewater-impacted streams

USGS Publications Warehouse

Barber, L.B.; Antweiler, Ronald C.; Flynn, J.L.; Keefe, S.H.; Kolpin, D.W.; Roth, D.A.; Schnoebelen, D.J.; Taylor, Howard E.; Verplanck, P.L.

2011-01-01

Understanding the potential effects of increased reliance on wastewater treatment plant (WWTP) effluents to meet municipal, agricultural, and environmental flow requires an understanding of the complex chemical loading characteristics of the WWTPs and the assimilative capacity of receiving waters. Stream ecosystem effects are linked to proportions of WWTP effluent under low-flow conditions as well as the nature of the effluent chemical mixtures. This study quantifies the loading of 58 inorganic constituents (nutrients to rare earth elements) from WWTP discharges relative to upstream landscape-based sources. Stream assimilation capacity was evaluated by Lagrangian sampling, using flow velocities determined from tracer experiments to track the same parcel of water as it moved downstream. Boulder Creek, Colorado and Fourmile Creek, Iowa, representing two different geologic and hydrologic landscapes, were sampled under low-flow conditions in the summer and spring. One-half of the constituents had greater loads from the WWTP effluents than the upstream drainages, and once introduced into the streams, dilution was the predominant assimilation mechanism. Only ammonium and bismuth had significant decreases in mass load downstream from the WWTPs during all samplings. The link between hydrology and water chemistry inherent in Lagrangian sampling allows quantitative assessment of chemical fate across different landscapes. ?? 2011 American Chemical Society.

Crustal deformation, the earthquake cycle, and models of viscoelastic flow in the asthenosphere

NASA Technical Reports Server (NTRS)

Cohen, S. C.; Kramer, M. J.

1983-01-01

The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.

Investigation of transient cavitating flow in viscoelastic pipes

NASA Astrophysics Data System (ADS)

Keramat, A.; Tijsseling, A. S.; Ahmadi, A.

2010-08-01

A study on water hammer in viscoelastic pipes when the fluid pressure drops to liquid vapour pressure is performed. Two important concepts including column separation and the effects of retarded strains in the pipe wall on the fluid response have been investigated separately in recent works, but there is some curiosity as to how the results for pressure and discharge are when column separation occurs in a viscoelastic pipe. For pipes made of plastic such as polyethylene (PE) and polyvinyl chloride (PVC), viscoelasticity is a crucial mechanical property which changes the hydraulic and structural transient responses. Based on previous developments in the analysis of water hammer, a model which is capable of analysing column separation in viscoelastic pipes is presented and used for solving the selected case studies. For the column-separation modelling the Discrete Vapour Cavity Model (DVCM) is utilised and the viscoelasticity property of the pipe wall is modelled by Kelvin-Voigt elements. The effects of viscoelasticity play an important role in the column separation phenomenon because it changes the water hammer fundamental frequency and so affects the time of opening or collapse of the cavities. Verification of the implemented computer code is performed for the effects of viscoelasticity and column separation - separately and simultaneously - using experimental results from the literature. In the provided examples the focus is placed on the simultaneous effect of viscoelasticity and column separation on the hydraulic transient response. The final conclusions drawn are that if rectangular grids are utilised the DVCM gives acceptable predictions of the phenomenon and that the pipe wall material's retarded behaviour strongly dampens the pressure spikes caused by column separation.

Coherent Lagrangian swirls among submesoscale motions.

PubMed

Beron-Vera, F J; Hadjighasem, A; Xia, Q; Olascoaga, M J; Haller, G

2018-03-05

The emergence of coherent Lagrangian swirls (CLSs) among submesoscale motions in the ocean is illustrated. This is done by applying recent nonlinear dynamics tools for Lagrangian coherence detection on a surface flow realization produced by a data-assimilative submesoscale-permitting ocean general circulation model simulation of the Gulf of Mexico. Both mesoscale and submesoscale CLSs are extracted. These extractions prove the relevance of coherent Lagrangian eddies detected in satellite-altimetry-based geostrophic flow data for the arguably more realistic ageostrophic multiscale flow.

Lagrangian averaging with geodesic mean

NASA Astrophysics Data System (ADS)

Oliver, Marcel

2017-11-01

This paper revisits the derivation of the Lagrangian averaged Euler (LAE), or Euler-α equations in the light of an intrinsic definition of the averaged flow map as the geodesic mean on the volume-preserving diffeomorphism group. Under the additional assumption that first-order fluctuations are statistically isotropic and transported by the mean flow as a vector field, averaging of the kinetic energy Lagrangian of an ideal fluid yields the LAE Lagrangian. The derivation presented here assumes a Euclidean spatial domain without boundaries.

Lagrangian averaging with geodesic mean.

PubMed

Oliver, Marcel

2017-11-01

This paper revisits the derivation of the Lagrangian averaged Euler (LAE), or Euler- α equations in the light of an intrinsic definition of the averaged flow map as the geodesic mean on the volume-preserving diffeomorphism group. Under the additional assumption that first-order fluctuations are statistically isotropic and transported by the mean flow as a vector field, averaging of the kinetic energy Lagrangian of an ideal fluid yields the LAE Lagrangian. The derivation presented here assumes a Euclidean spatial domain without boundaries.

A Generalized Eulerian-Lagrangian Analysis, with Application to Liquid Flows with Vapor Bubbles

NASA Technical Reports Server (NTRS)

Dejong, Frederik J.; Meyyappan, Meyya

1993-01-01

Under a NASA MSFC SBIR Phase 2 effort an analysis has been developed for liquid flows with vapor bubbles such as those in liquid rocket engine components. The analysis is based on a combined Eulerian-Lagrangian technique, in which Eulerian conservation equations are solved for the liquid phase, while Lagrangian equations of motion are integrated in computational coordinates for the vapor phase. The novel aspect of the Lagrangian analysis developed under this effort is that it combines features of the so-called particle distribution approach with those of the so-called particle trajectory approach and can, in fact, be considered as a generalization of both of those traditional methods. The result of this generalization is a reduction in CPU time and memory requirements. Particle time step (stability) limitations have been eliminated by semi-implicit integration of the particle equations of motion (and, for certain applications, the particle temperature equation), although practical limitations remain in effect for reasons of accuracy. The analysis has been applied to the simulation of cavitating flow through a single-bladed section of a labyrinth seal. Models for the simulation of bubble formation and growth have been included, as well as models for bubble drag and heat transfer. The results indicate that bubble formation is more or less 'explosive'. for a given flow field, the number density of bubble nucleation sites is very sensitive to the vapor properties and the surface tension. The bubble motion, on the other hand, is much less sensitive to the properties, but is affected strongly by the local pressure gradients in the flow field. In situations where either the material properties or the flow field are not known with sufficient accuracy, parametric studies can be carried out rapidly to assess the effect of the important variables. Future work will include application of the analysis to cavitation in inducer flow fields.

Numerical study on flow over stepped spillway using Lagrangian method

NASA Astrophysics Data System (ADS)

Wang, Junmin; Fu, Lei; Xu, Haibo; Jin, Yeechung

2018-02-01

Flow over stepped spillway has been studied for centuries, due to its unstable and the characteristics of cavity, the simulation of this type of spillway flow is always difficult. Most of the early studies of flow over stepped spillway are based on experiment, while in the recent decades, numerical studies of flow over stepped spillway draw most of the researchers’ attentions due to its simplicity and efficiency. In this study, a new Lagrangian based particle method is introduced to reproduce the phenomenon of flow over stepped spillway, the inherent advantages of this particle based method provide a convincing free surface and velocity profiles compared with previous experimental data. The capacity of this new method is proved and it is anticipated to be an alternative tool of traditional mesh based method in environmental engineering field such as the simulation of flow over stepped spillway.

A Constructive Approach to Regularity of Lagrangian Trajectories for Incompressible Euler Flow in a Bounded Domain

NASA Astrophysics Data System (ADS)

Besse, Nicolas; Frisch, Uriel

2017-04-01

The 3D incompressible Euler equations are an important research topic in the mathematical study of fluid dynamics. Not only is the global regularity for smooth initial data an open issue, but the behaviour may also depend on the presence or absence of boundaries. For a good understanding, it is crucial to carry out, besides mathematical studies, high-accuracy and well-resolved numerical exploration. Such studies can be very demanding in computational resources, but recently it has been shown that very substantial gains can be achieved first, by using Cauchy's Lagrangian formulation of the Euler equations and second, by taking advantage of analyticity results of the Lagrangian trajectories for flows whose initial vorticity is Hölder-continuous. The latter has been known for about 20 years (Serfati in J Math Pures Appl 74:95-104, 1995), but the combination of the two, which makes use of recursion relations among time-Taylor coefficients to obtain constructively the time-Taylor series of the Lagrangian map, has been achieved only recently (Frisch and Zheligovsky in Commun Math Phys 326:499-505, 2014; Podvigina et al. in J Comput Phys 306:320-342, 2016 and references therein). Here we extend this methodology to incompressible Euler flow in an impermeable bounded domain whose boundary may be either analytic or have a regularity between indefinite differentiability and analyticity. Non-constructive regularity results for these cases have already been obtained by Glass et al. (Ann Sci Éc Norm Sup 45:1-51, 2012). Using the invariance of the boundary under the Lagrangian flow, we establish novel recursion relations that include contributions from the boundary. This leads to a constructive proof of time-analyticity of the Lagrangian trajectories with analytic boundaries, which can then be used subsequently for the design of a very high-order Cauchy-Lagrangian method.

An extended Lagrangian method

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

1992-01-01

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method', is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. Meanwhile, it also avoids the inaccuracy incurred due to geometry and variable interpolations used by the previous Lagrangian methods. Unlike the Lagrangian method previously imposed which is valid only for supersonic flows, the present method is general and capable of treating subsonic flows as well as supersonic flows. The method proposed in this paper is robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Shear and shearless Lagrangian structures in compound channels

NASA Astrophysics Data System (ADS)

Enrile, F.; Besio, G.; Stocchino, A.

2018-03-01

Transport processes in a physical model of a natural stream with a composite cross-section (compound channel) are investigated by means of a Lagrangian analysis based on nonlinear dynamical system theory. Two-dimensional free surface Eulerian experimental velocity fields of a uniform flow in a compound channel form the basis for the identification of the so-called Lagrangian Coherent Structures. Lagrangian structures are recognized as the key features that govern particle trajectories. We seek for two particular class of Lagrangian structures: Shear and shearless structures. The former are generated whenever the shear dominates the flow whereas the latter behave as jet-cores. These two type of structures are detected as ridges and trenches of the Finite-Time Lyapunov Exponents fields, respectively. Besides, shearlines computed applying the geodesic theory of transport barriers mark Shear Lagrangian Coherent Structures. So far, the detection of these structures in real experimental flows has not been deeply investigated. Indeed, the present results obtained in a wide range of the controlling parameters clearly show a different behaviour depending on the shallowness of the flow. Shear and Shearless Lagrangian Structures detected from laboratory experiments clearly appear as the flow develops in shallow conditions. The presence of these Lagrangian Structures tends to fade in deep flow conditions.

Lagrangian turbulence near walls: Structures and mixing in admissible model flows

NASA Astrophysics Data System (ADS)

Ottino, J. M.

1989-05-01

The general objective of work during this period was to bridge the gap between modern ideas from dynamical systems and chaos and more traditional approaches to turbulence. In order to reach this objective we conducted theoretical and computational work on two systems: a perturbed Kelvin cat eyes flow, and prototype solutions of the Navier-Stokes equations near solid walls. The main results obtained are two-fold: production flows capable of producing complex distributions of vorticity, and constructed flow fields, based on solutions of the Navier Stokes equations, which are capable of displaying both Eulerian and Lagrangian turbulence.

Transient and steady state viscoelastic rolling contact

NASA Technical Reports Server (NTRS)

Padovan, J.; Paramadilok, O.

1985-01-01

Based on moving total Lagrangian coordinates, a so-called traveling Hughes type contact strategy is developed. Employing the modified contact scheme in conjunction with a traveling finite element strategy, an overall solution methodology is developed to handle transient and steady viscoelastic rolling contact. To verify the scheme, the results of both experimental and analytical benchmarking is presented. The experimental benchmarking includes the handling of rolling tires up to their upper bound behavior, namely the standing wave response.

Lagrangian postprocessing of computational hemodynamics.

PubMed

Shadden, Shawn C; Arzani, Amirhossein

2015-01-01

Recent advances in imaging, modeling, and computing have rapidly expanded our capabilities to model hemodynamics in the large vessels (heart, arteries, and veins). This data encodes a wealth of information that is often under-utilized. Modeling (and measuring) blood flow in the large vessels typically amounts to solving for the time-varying velocity field in a region of interest. Flow in the heart and larger arteries is often complex, and velocity field data provides a starting point for investigating the hemodynamics. This data can be used to perform Lagrangian particle tracking, and other Lagrangian-based postprocessing. As described herein, Lagrangian methods are necessary to understand inherently transient hemodynamic conditions from the fluid mechanics perspective, and to properly understand the biomechanical factors that lead to acute and gradual changes of vascular function and health. The goal of the present paper is to review Lagrangian methods that have been used in post-processing velocity data of cardiovascular flows.

Lagrangian postprocessing of computational hemodynamics

PubMed Central

Shadden, Shawn C.; Arzani, Amirhossein

2014-01-01

Recent advances in imaging, modeling and computing have rapidly expanded our capabilities to model hemodynamics in the large vessels (heart, arteries and veins). This data encodes a wealth of information that is often under-utilized. Modeling (and measuring) blood flow in the large vessels typically amounts to solving for the time-varying velocity field in a region of interest. Flow in the heart and larger arteries is often complex, and velocity field data provides a starting point for investigating the hemodynamics. This data can be used to perform Lagrangian particle tracking, and other Lagrangian-based postprocessing. As described herein, Lagrangian methods are necessary to understand inherently transient hemodynamic conditions from the fluid mechanics perspective, and to properly understand the biomechanical factors that lead to acute and gradual changes of vascular function and health. The goal of the present paper is to review Lagrangian methods that have been used in post-processing velocity data of cardiovascular flows. PMID:25059889

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

NASA Technical Reports Server (NTRS)

Apte, S. V.; Mahesh, K.; Lundgren, T.

2003-01-01

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Finite volume multigrid method of the planar contraction flow of a viscoelastic fluid

NASA Astrophysics Data System (ADS)

Moatssime, H. Al; Esselaoui, D.; Hakim, A.; Raghay, S.

2001-08-01

This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first-order upwind approximation for the viscoelastic stress. A non-uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non-linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss-Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd-B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright

"Lagrangian" for a Non-Lagrangian Field Theory with N=2 Supersymmetry.

PubMed

Gadde, Abhijit; Razamat, Shlomo S; Willett, Brian

2015-10-23

We suggest that at least some of the strongly coupled N=2 quantum field theories in 4D can have a nonconformal N=1 Lagrangian description flowing to them at low energies. In particular, we construct such a description for the N=2 rank one superconformal field theory with E(6) flavor symmetry, for which a Lagrangian description was previously unavailable. We utilize this description to compute several supersymmetric partition functions.

Meshless Lagrangian SPH method applied to isothermal lid-driven cavity flow at low-Re numbers

NASA Astrophysics Data System (ADS)

Fraga Filho, C. A. D.; Chacaltana, J. T. A.; Pinto, W. J. N.

2018-01-01

SPH is a recent particle method applied in the cavities study, without many results available in the literature. The lid-driven cavity flow is a classic problem of the fluid mechanics, extensively explored in the literature and presenting a considerable complexity. The aim of this paper is to present a solution from the Lagrangian viewpoint for this problem. The discretization of the continuum domain is performed using the Lagrangian particles. The physical laws of mass, momentum and energy conservation are presented by the Navier-Stokes equations. A serial numerical code, written in Fortran programming language, has been used to perform the numerical simulations. The application of the SPH and comparison with the literature (mesh methods and a meshless collocation method) have been done. The positions of the primary vortex centre and the non-dimensional velocity profiles passing through the geometric centre of the cavity have been analysed. The numerical Lagrangian results showed a good agreement when compared to the results found in the literature, specifically for { Re} < 100.00 . Suggestions for improvements in the SPH model presented are listed, in the search for better results for flows with higher Reynolds numbers.

Verification assessment of piston boundary conditions for Lagrangian simulation of compressible flow similarity solutions

DOE PAGES

Ramsey, Scott D.; Ivancic, Philip R.; Lilieholm, Jennifer F.

2015-12-10

This work is concerned with the use of similarity solutions of the compressible flow equations as benchmarks or verification test problems for finite-volume compressible flow simulation software. In practice, this effort can be complicated by the infinite spatial/temporal extent of many candidate solutions or “test problems.” Methods can be devised with the intention of ameliorating this inconsistency with the finite nature of computational simulation; the exact strategy will depend on the code and problem archetypes under investigation. For example, self-similar shock wave propagation can be represented in Lagrangian compressible flow simulations as rigid boundary-driven flow, even if no such “piston”more » is present in the counterpart mathematical similarity solution. The purpose of this work is to investigate in detail the methodology of representing self-similar shock wave propagation as a piston-driven flow in the context of various test problems featuring simple closed-form solutions of infinite spatial/temporal extent. The closed-form solutions allow for the derivation of similarly closed-form piston boundary conditions (BCs) for use in Lagrangian compressible flow solvers. Finally, the consequences of utilizing these BCs (as opposed to directly initializing the self-similar solution in a computational spatial grid) are investigated in terms of common code verification analysis metrics (e.g., shock strength/position errors and global convergence rates).« less

Verification assessment of piston boundary conditions for Lagrangian simulation of compressible flow similarity solutions

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

Ramsey, Scott D.; Ivancic, Philip R.; Lilieholm, Jennifer F.

This work is concerned with the use of similarity solutions of the compressible flow equations as benchmarks or verification test problems for finite-volume compressible flow simulation software. In practice, this effort can be complicated by the infinite spatial/temporal extent of many candidate solutions or “test problems.” Methods can be devised with the intention of ameliorating this inconsistency with the finite nature of computational simulation; the exact strategy will depend on the code and problem archetypes under investigation. For example, self-similar shock wave propagation can be represented in Lagrangian compressible flow simulations as rigid boundary-driven flow, even if no such “piston”more » is present in the counterpart mathematical similarity solution. The purpose of this work is to investigate in detail the methodology of representing self-similar shock wave propagation as a piston-driven flow in the context of various test problems featuring simple closed-form solutions of infinite spatial/temporal extent. The closed-form solutions allow for the derivation of similarly closed-form piston boundary conditions (BCs) for use in Lagrangian compressible flow solvers. Finally, the consequences of utilizing these BCs (as opposed to directly initializing the self-similar solution in a computational spatial grid) are investigated in terms of common code verification analysis metrics (e.g., shock strength/position errors and global convergence rates).« less

Investigation of pitchfork bifurcation phenomena effects on heat transfer of viscoelastic flow inside a symmetric sudden expansion

NASA Astrophysics Data System (ADS)

Shahbani-Zahiri, A.; Hassanzadeh, H.; Shahmardan, M. M.; Norouzi, M.

2017-11-01

In this paper, the inertial and non-isothermal flows of the viscoelastic fluid through a planar channel with symmetric sudden expansion are numerically simulated. Effects of pitchfork bifurcation phenomena on the heat transfer rate are examined for the thermally developing and fully developed flow of the viscoelastic fluid inside the expanded part of the planar channel with an expansion ratio of 1:3. The rheological model of exponential Phan Thien-Tanner is used to include both the effects of shear-thinning and elasticity in fluid viscosity. The properties of fluids are temperature-dependent, and the viscous dissipation and heat stored by fluid elasticity are considered in the heat transfer equation. For coupling the governing equations, the PISO algorithm (Pressure Implicit with Splitting of Operator) is applied and the system of equations is linearized using the finite volume method on the collocated grids. The main purpose of this study is to examine the pitchfork bifurcation phenomena and its influences on the temperature distribution, the local and mean Nusselt numbers, and the first and second normal stress differences at different Reynolds, elasticity, and Brinkman numbers. The results show that by increasing the Brinkman number for the heated flow of the viscoelastic fluid inside the expanded part of the channel, the value of the mean Nusselt number is almost linearly decreased. Also, the maximum values of the local Nusselt number for the thermally developing flow and the local Nusselt number of the thermally fully developed flow are decremented by enhancing the Brinkman number.

Tests of dynamic Lagrangian eddy viscosity models in Large Eddy Simulations of flow over three-dimensional bluff bodies

NASA Astrophysics Data System (ADS)

Tseng, Yu-Heng; Meneveau, Charles; Parlange, Marc B.

2004-11-01

Large Eddy Simulations (LES) of atmospheric boundary-layer air movement in urban environments are especially challenging due to complex ground topography. Typically in such applications, fairly coarse grids must be used where the subgrid-scale (SGS) model is expected to play a crucial role. A LES code using pseudo-spectral discretization in horizontal planes and second-order differencing in the vertical is implemented in conjunction with the immersed boundary method to incorporate complex ground topography, with the classic equilibrium log-law boundary condition in the new-wall region, and with several versions of the eddy-viscosity model: (1) the constant-coefficient Smagorinsky model, (2) the dynamic, scale-invariant Lagrangian model, and (3) the dynamic, scale-dependent Lagrangian model. Other planar-averaged type dynamic models are not suitable because spatial averaging is not possible without directions of statistical homogeneity. These SGS models are tested in LES of flow around a square cylinder and of flow over surface-mounted cubes. Effects on the mean flow are documented and found not to be major. Dynamic Lagrangian models give a physically more realistic SGS viscosity field, and in general, the scale-dependent Lagrangian model produces larger Smagorinsky coefficient than the scale-invariant one, leading to reduced distributions of resolved rms velocities especially in the boundary layers near the bluff bodies.

Transient Droplet Behavior and Droplet Breakup during Bulk and Confined Shear Flow in Blends with One Viscoelastic Component: Experiments, Modelling and Simulations

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

Cardinaels, Ruth; Verhulst, Kristof; Moldenaers, Paula

2008-07-07

The transient droplet deformation and droplet orientation after inception of shear, the shape relaxation after cessation of shear and droplet breakup during shear, are microscopically studied, both under bulk and confined conditions. The studied blends contain one viscoelastic Boger fluid phase. A counter rotating setup, based on a Paar Physica MCR300, is used for the droplet visualisation. For bulk shear flow, it is shown that the droplet deformation during startup of shear flow and the shape relaxation after cessation of shear flow are hardly influenced by droplet viscoelasticity, even at moderate to high capillary and Deborah numbers. The effects ofmore » droplet viscoelasticity only become visible close to the critical conditions and a novel break-up mechanism is observed. Matrix viscoelasticity has a more pronounced effect, causing overshoots in the deformation and significantly inhibiting relaxation. However, different applied capillary numbers prior to cessation of shear flow, with the Deborah number fixed, still result in a single master curve for shape retraction, as in fully Newtonian systems. The long tail in the droplet relaxation can be qualitatively described with a phenomenological model for droplet deformation, when using a 5-mode Giesekus model for the fluid rheology. It is found that the shear flow history significantly affects the droplet shape evolution and the breakup process in blends with one viscoelastic component. Confining a droplet between two plates accelerates the droplet deformation kinetics, similar to fully Newtonian systems. However, the increased droplet deformation, due to wall effects, causes the steady state to be reached at a later instant in time. Droplet relaxation is less sensitive to confinement, leading to slower relaxation kinetics only for highly confined droplets. For the blend with a viscoelastic droplet, a non-monotonous trend is found for the critical capillary number as a function of the confinement ratio. Finally

A new Lagrangian random choice method for steady two-dimensional supersonic/hypersonic flow

NASA Technical Reports Server (NTRS)

Loh, C. Y.; Hui, W. H.

1991-01-01

Glimm's (1965) random choice method has been successfully applied to compute steady two-dimensional supersonic/hypersonic flow using a new Lagrangian formulation. The method is easy to program, fast to execute, yet it is very accurate and robust. It requires no grid generation, resolves slipline and shock discontinuities crisply, can handle boundary conditions most easily, and is applicable to hypersonic as well as supersonic flow. It represents an accurate and fast alternative to the existing Eulerian methods. Many computed examples are given.

Parallel computing using a Lagrangian formulation

NASA Technical Reports Server (NTRS)

Liou, May-Fun; Loh, Ching Yuen

1991-01-01

A new Lagrangian formulation of the Euler equation is adopted for the calculation of 2-D supersonic steady flow. The Lagrangian formulation represents the inherent parallelism of the flow field better than the common Eulerian formulation and offers a competitive alternative on parallel computers. The implementation of the Lagrangian formulation on the Thinking Machines Corporation CM-2 Computer is described. The program uses a finite volume, first-order Godunov scheme and exhibits high accuracy in dealing with multidimensional discontinuities (slip-line and shock). By using this formulation, a better than six times speed-up was achieved on a 8192-processor CM-2 over a single processor of a CRAY-2.

Parallel computing using a Lagrangian formulation

NASA Technical Reports Server (NTRS)

Liou, May-Fun; Loh, Ching-Yuen

1992-01-01

This paper adopts a new Lagrangian formulation of the Euler equation for the calculation of two dimensional supersonic steady flow. The Lagrangian formulation represents the inherent parallelism of the flow field better than the common Eulerian formulation and offers a competitive alternative on parallel computers. The implementation of the Lagrangian formulation on the Thinking Machines Corporation CM-2 Computer is described. The program uses a finite volume, first-order Godunov scheme and exhibits high accuracy in dealing with multidimensional discontinuities (slip-line and shock). By using this formulation, we have achieved better than six times speed-up on a 8192-processor CM-2 over a single processor of a CRAY-2.

Uncoupled poroelastic and intrinsic viscoelastic dissipation in cartilage.

PubMed

Han, Guebum; Hess, Cole; Eriten, Melih; Henak, Corinne R

2018-04-26

This paper studies uncoupled poroelastic (flow-dependent) and intrinsic viscoelastic (flow-independent) energy dissipation mechanisms via their dependence on characteristic lengths to understand the root of cartilage's broadband dissipation behavior. Phase shift and dynamic modulus were measured from dynamic microindentation tests conducted on hydrated cartilage at different contact radii, as well as on dehydrated cartilage. Cartilage weight and thickness were recorded during dehydration. Phase shifts revealed poroelastic- and viscoelastic-dominant dissipation regimes in hydrated cartilage. Specifically, phase shift at a relatively small radius was governed by poroviscoelasticity, while phase shift at a relatively large radius was dominantly governed by intrinsic viscoelasticity. The uncoupled dissipation mechanisms demonstrated that intrinsic viscoelastic dissipation provided sustained broadband dissipation for all length scales, and additional poroelastic dissipation increased total dissipation at small length scales. Dehydration decreased intrinsic viscoelastic dissipation of cartilage. The findings demonstrated a possibility to measure poroelastic and intrinsic viscoelastic properties of cartilage at similar microscale lengths. Also they encouraged development of broadband cartilage like-dampers and provided important design parameters to maximize their performance. Copyright © 2018 Elsevier Ltd. All rights reserved.

A mass and momentum conserving unsplit semi-Lagrangian framework for simulating multiphase flows

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

Owkes, Mark, E-mail: mark.owkes@montana.edu; Desjardins, Olivier

In this work, we present a computational methodology for convection and advection that handles discontinuities with second order accuracy and maintains conservation to machine precision. This method can transport a variety of discontinuous quantities and is used in the context of an incompressible gas–liquid flow to transport the phase interface, momentum, and scalars. The proposed method provides a modification to the three-dimensional, unsplit, second-order semi-Lagrangian flux method of Owkes & Desjardins (JCP, 2014). The modification adds a refined grid that provides consistent fluxes of mass and momentum defined on a staggered grid and discrete conservation of mass and momentum, evenmore » for flows with large density ratios. Additionally, the refined grid doubles the resolution of the interface without significantly increasing the computational cost over previous non-conservative schemes. This is possible due to a novel partitioning of the semi-Lagrangian fluxes into a small number of simplices. The proposed scheme is tested using canonical verification tests, rising bubbles, and an atomizing liquid jet.« less

Unsteady Magnetized Flow and Heat Transfer of a Viscoelastic fluid over a Stretching Surface

NASA Astrophysics Data System (ADS)

Ghosh, Sushil Kumar

2017-12-01

This paper is to study the flow of heated ferro-fluid over a stretching sheet under the influence of magnetic field. The fluid considered in the present investigation is a mixture of blood as well as fluid-dispersed magnetic nano particles and under this context blood is found to be the appropriate choice of viscoelastic, Walter's B fluid. The objective of the present work is to study the effect of various parameters found in the mathematical analysis. Taking into account the blood has zero electrical conductivity, magnetization effect has been considered in the governing equation of the present study with the use of ferro-fluid dynamics principle. By introducing appropriate non-dimensional variables into the governing equations of unsteady two-dimensional flow of viscoelastic fluid with heat transfer are converted to a set of ordinary differential equations with appropriate boundary conditions. Newton's linearization technique has been employed for the solution of non-linear ordinary differential equations. Important results found in the present investigation are the substantial influence of ferro-magnetic parameter, Prandlt number and the parameter associated with the thermal conductivity on the flow and heat transfer. It is observed that the presence of magnetic dipole essentially reduces the flow velocity in the vertical direction and that helps to damage the cancer cells in the tumor region.

Bayesian Lagrangian Data Assimilation and Drifter Deployment Strategies

NASA Astrophysics Data System (ADS)

Dutt, A.; Lermusiaux, P. F. J.

2017-12-01

Ocean currents transport a variety of natural (e.g. water masses, phytoplankton, zooplankton, sediments, etc.) and man-made materials and other objects (e.g. pollutants, floating debris, search and rescue, etc.). Lagrangian Coherent Structures (LCSs) or the most influential/persistent material lines in a flow, provide a robust approach to characterize such Lagrangian transports and organize classic trajectories. Using the flow-map stochastic advection and a dynamically-orthogonal decomposition, we develop uncertainty prediction schemes for both Eulerian and Lagrangian variables. We then extend our Bayesian Gaussian Mixture Model (GMM)-DO filter to a joint Eulerian-Lagrangian Bayesian data assimilation scheme. The resulting nonlinear filter allows the simultaneous non-Gaussian estimation of Eulerian variables (e.g. velocity, temperature, salinity, etc.) and Lagrangian variables (e.g. drifter/float positions, trajectories, LCSs, etc.). Its results are showcased using a double-gyre flow with a random frequency, a stochastic flow past a cylinder, and realistic ocean examples. We further show how our Bayesian mutual information and adaptive sampling equations provide a rigorous efficient methodology to plan optimal drifter deployment strategies and predict the optimal times, locations, and types of measurements to be collected.

Modeling of confined turbulent fluid-particle flows using Eulerian and Lagrangian schemes

NASA Technical Reports Server (NTRS)

Adeniji-Fashola, A.; Chen, C. P.

1990-01-01

Two important aspects of fluid-particulate interaction in dilute gas-particle turbulent flows (the turbulent particle dispersion and the turbulence modulation effects) are addressed, using the Eulerian and Lagrangian modeling approaches to describe the particulate phase. Gradient-diffusion approximations are employed in the Eulerian formulation, while a stochastic procedure is utilized to simulate turbulent dispersion in the Lagrangina formulation. The k-epsilon turbulence model is used to characterize the time and length scales of the continuous phase turbulence. Models proposed for both schemes are used to predict turbulent fully-developed gas-solid vertical pipe flow with reasonable accuracy.

Direct numerical simulation of viscoelastic-fluid-based nanofluid turbulent channel flow with heat transfer

NASA Astrophysics Data System (ADS)

Yang, Juan-Cheng; Li, Feng-Chen; Cai, Wei-Hua; Zhang, Hong-Na; Yu, Bo

2015-08-01

Our previous experimental studies have confirmed that viscoelastic-fluid-based nanofluid (VFBN) prepared by suspending nanoparticles in a viscoelastic base fluid (VBF, behaves drag reduction at turbulent flow state) can reduce turbulent flow resistance as compared with water and enhance heat transfer as compared with VBF. Direct numerical simulation (DNS) is performed in this study to explore the mechanisms of heat transfer enhancement (HTE) and flow drag reduction (DR) for the VFBN turbulent flow. The Giesekus model is used as the constitutive equation for VFBN. Our previously proposed thermal dispersion model is adopted to take into account the thermal dispersion effects of nanoparticles in the VFBN turbulent flow. The DNS results show similar behaviors for flow resistance and heat transfer to those obtained in our previous experiments. Detailed analyses are conducted for the turbulent velocity, temperature, and conformation fields obtained by DNSs for different fluid cases, and for the friction factor with viscous, turbulent, and elastic contributions and heat transfer rate with conductive, turbulent and thermal dispersion contributions of nanoparticles, respectively. The mechanisms of HTE and DR of VFBN turbulent flows are then discussed. Based on analogy theory, the ratios of Chilton-Colburn factor to friction factor for different fluid flow cases are investigated, which from another aspect show the significant enhancement in heat transfer performance for some cases of water-based nanofluid and VFBN turbulent flows. Project supported by the National Natural Science Foundation of China (Grant No. 51276046), the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20112302110020), the China Postdoctoral Science Foundation (Grant No. 2014M561037), and the President Fund of University of Chinese Academy of Sciences, China (Grant No. Y3510213N00).

An extended Lagrangian method for subsonic flows

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing; Loh, Ching Y.

1992-01-01

It is well known that fluid motion can be specified by either the Eulerian of Lagrangian description. Most of Computational Fluid Dynamics (CFD) developments over the last three decades have been based on the Eulerian description and considerable progress has been made. In particular, the upwind methods, inspired and guided by the work of Gudonov, have met with many successes in dealing with complex flows, especially where discontinuities exist. However, this shock capturing property has proven to be accurate only when the discontinuity is aligned with one of the grid lines since most upwind methods are strictly formulated in 1-D framework and only formally extended to multi-dimensions. Consequently, the attractive property of crisp resolution of these discontinuities is lost and research on genuine multi-dimensional approach has just been undertaken by several leading researchers. Nevertheless they are still based on the Eulerian description.

On flows of viscoelastic fluids under threshold-slip boundary conditions

NASA Astrophysics Data System (ADS)

Baranovskii, E. S.

2018-03-01

We investigate a boundary-value problem for the steady isothermal flow of an incompressible viscoelastic fluid of Oldroyd type in a 3D bounded domain with impermeable walls. We use the Fujita threshold-slip boundary condition. This condition states that the fluid can slip along a solid surface when the shear stresses reach a certain critical value; otherwise the slipping velocity is zero. Assuming that the flow domain is not rotationally symmetric, we prove an existence theorem for the corresponding slip problem in the framework of weak solutions. The proof uses methods for solving variational inequalities with pseudo-monotone operators and convex functionals, the method of introduction of auxiliary viscosity, as well as a passage-to-limit procedure based on energy estimates of approximate solutions, Korn’s inequality, and compactness arguments. Also, some properties and estimates of weak solutions are established.

Direct experimental visualization of the global Hamiltonian progression of two-dimensional Lagrangian flow topologies from integrable to chaotic state.

PubMed

Baskan, O; Speetjens, M F M; Metcalfe, G; Clercx, H J H

2015-10-01

Countless theoretical/numerical studies on transport and mixing in two-dimensional (2D) unsteady flows lean on the assumption that Hamiltonian mechanisms govern the Lagrangian dynamics of passive tracers. However, experimental studies specifically investigating said mechanisms are rare. Moreover, they typically concern local behavior in specific states (usually far away from the integrable state) and generally expose this indirectly by dye visualization. Laboratory experiments explicitly addressing the global Hamiltonian progression of the Lagrangian flow topology entirely from integrable to chaotic state, i.e., the fundamental route to efficient transport by chaotic advection, appear non-existent. This motivates our study on experimental visualization of this progression by direct measurement of Poincaré sections of passive tracer particles in a representative 2D time-periodic flow. This admits (i) accurate replication of the experimental initial conditions, facilitating true one-to-one comparison of simulated and measured behavior, and (ii) direct experimental investigation of the ensuing Lagrangian dynamics. The analysis reveals a close agreement between computations and observations and thus experimentally validates the full global Hamiltonian progression at a great level of detail.

Lagrangian motion, coherent structures, and lines of persistent material strain.

PubMed

Samelson, R M

2013-01-01

Lagrangian motion in geophysical fluids may be strongly influenced by coherent structures that support distinct regimes in a given flow. The problems of identifying and demarcating Lagrangian regime boundaries associated with dynamical coherent structures in a given velocity field can be studied using approaches originally developed in the context of the abstract geometric theory of ordinary differential equations. An essential insight is that when coherent structures exist in a flow, Lagrangian regime boundaries may often be indicated as material curves on which the Lagrangian-mean principal-axis strain is large. This insight is the foundation of many numerical techniques for identifying such features in complex observed or numerically simulated ocean flows. The basic theoretical ideas are illustrated with a simple, kinematic traveling-wave model. The corresponding numerical algorithms for identifying candidate Lagrangian regime boundaries and lines of principal Lagrangian strain (also called Lagrangian coherent structures) are divided into parcel and bundle schemes; the latter include the finite-time and finite-size Lyapunov exponent/Lagrangian strain (FTLE/FTLS and FSLE/FSLS) metrics. Some aspects and results of oceanographic studies based on these approaches are reviewed, and the results are discussed in the context of oceanographic observations of dynamical coherent structures.

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

PubMed

Holm, Darryl D.

2002-06-01

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

Eulerian and Lagrangian methods for vortex tracking in 2D and 3D flows

NASA Astrophysics Data System (ADS)

Huang, Yangzi; Green, Melissa

2014-11-01

Coherent structures are a key component of unsteady flows in shear layers. Improvement of experimental techniques has led to larger amounts of data and requires of automated procedures for vortex tracking. Many vortex criteria are Eulerian, and identify the structures by an instantaneous local swirling motion in the field, which are indicated by closed or spiral streamlines or pathlines in a reference frame. Alternatively, a Lagrangian Coherent Structures (LCS) analysis is a Lagrangian method based on the quantities calculated along fluid particle trajectories. In the current work, vortex detection is demonstrated on data from the simulation of two cases: a 2D flow with a flat plate undergoing a 45 ° pitch-up maneuver and a 3D wall-bounded turbulence channel flow. Vortices are visualized and tracked by their centers and boundaries using Γ1, the Q criterion, and LCS saddle points. In the cases of 2D flow, saddle points trace showed a rapid acceleration of the structure which indicates the shedding from the plate. For channel flow, saddle points trace shows that average structure convection speed exhibits a similar trend as a function of wall-normal distance as the mean velocity profile, and leads to statistical quantities of vortex dynamics. Dr. Jeff Eldredge and his research group at UCLA are gratefully acknowledged for sharing the database of simulation for the current research. This work was supported by the Air Force Office of Scientific Research under AFOSR Award No. FA9550-14-1-0210.

An Extended Lagrangian Method

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

1995-01-01

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Magnetohydrodynamics effect on convective boundary layer flow and heat transfer of viscoelastic micropolar fluid past a sphere

NASA Astrophysics Data System (ADS)

Amera Aziz, Laila; Kasim, Abdul Rahman Mohd; Zuki Salleh, Mohd; Syahidah Yusoff, Nur; Shafie, Sharidan

2017-09-01

The main interest of this study is to investigate the effect of MHD on the boundary layer flow and heat transfer of viscoelastic micropolar fluid. Governing equations are transformed into dimensionless form in order to reduce their complexity. Then, the stream function is applied to the dimensionless equations to produce partial differential equations which are then solved numerically using the Keller-box method in Fortran programming. The numerical results are compared to published study to ensure the reliability of present results. The effects of selected physical parameters such as the viscoelastic parameter, K, micropolar parameter, K1 and magnetic parameter, M on the flow and heat transfer are discussed and presented in tabular and graphical form. The findings from this study will be of critical importance in the fields of medicine, chemical as well as industrial processes where magnetic field is involved.

Lagrangian solution of supersonic real gas flows

NASA Technical Reports Server (NTRS)

Loh, Ching-Yuen; Liou, Meng-Sing

1993-01-01

The present extention of a Lagrangian approach of the Riemann solution procedure, which was originally proposed for perfect gases, to real gases, is nontrivial and requires the development of an exact real-gas Riemann solver for the Lagrangian form of the conservation laws. Calculations including complex wave interactions of various types were conducted to test the accuracy and robustness of the approach. Attention is given to the case of 2D oblique waves' capture, where a slip line is clearly in evidence; the real gas effect is demonstrated in the case of a generic engine nozzle.

Dynamics and structures of transitional viscoelastic turbulence in channel flow

NASA Astrophysics Data System (ADS)

Shekar, Ashwin; Wang, Sung-Ning; Graham, Michael

2017-11-01

Introducing a trace amount of polymer into turbulent flows can result in a substantial reduction of drag. However, the mechanism is not fully understood at high levels of drag reduction. In this work we perform direct numerical simulations (DNS) of viscoelastic channel flow turbulence using a scheme that guarantees the positive-definiteness of polymer conformation tensor without artificial diffusion. Here we present the results of two parametric studies with the bulk Reynolds number fixed at 2000. First, the Weissenberg number (Wi) is kept at 100 and we vary the viscosity ratio (ratio ratio of the solvent viscosity and the total viscosity). Maximum drag reduction (MDR) is observed with viscosity ratio <0.95. As we decrease the viscosity ratio, i.e. increase polymer concentration, the mean velocity profile is almost invariant. However, this is accompanied by a decrease in velocity fluctuations but the flow stays turbulent. Turbulent kinetic energy budget analysis shows that, in this parameter regime, polymer becomes the major source of velocity fluctuations, replacing the energy transfer from the mean flow. In the second study, we fix the viscosity ratio at 0.95 and trace the Wi up to this regime and present the accompanying changes in flow quantities and structures.

Lagrangian statistics in compressible isotropic homogeneous turbulence

NASA Astrophysics Data System (ADS)

Yang, Yantao; Wang, Jianchun; Shi, Yipeng; Chen, Shiyi

2011-11-01

In this work we conducted the Direct Numerical Simulation (DNS) of a forced compressible isotropic homogeneous turbulence and investigated the flow statistics from the Lagrangian point of view, namely the statistics is computed following the passive tracers trajectories. The numerical method combined the Eulerian field solver which was developed by Wang et al. (2010, J. Comp. Phys., 229, 5257-5279), and a Lagrangian module for tracking the tracers and recording the data. The Lagrangian probability density functions (p.d.f.'s) have then been calculated for both kinetic and thermodynamic quantities. In order to isolate the shearing part from the compressing part of the flow, we employed the Helmholtz decomposition to decompose the flow field (mainly the velocity field) into the solenoidal and compressive parts. The solenoidal part was compared with the incompressible case, while the compressibility effect showed up in the compressive part. The Lagrangian structure functions and cross-correlation between various quantities will also be discussed. This work was supported in part by the China's Turbulence Program under Grant No.2009CB724101.

Correlating Lagrangian structures with forcing in two-dimensional flow

NASA Astrophysics Data System (ADS)

Ouellette, Nicholas; Hogg, Charlie; Liao, Yang

2015-11-01

Lagrangian coherent structures (LCSs) are the dominant transport barriers in unsteady, aperiodic flows, and their role in organizing mixing and transport has been well documented. However, nearly all that is known about LCSs has been gleaned from passive observations: they are computed in a post-processing step after a flow has been observed, and used to understand why the mixing and transport proceeded as it did. Here, we instead take a first step toward controlling the presence or locations of LCSs by studying the relationship between LCSs and external forcing in an experimental quasi-two-dimensional weakly turbulent flow. We find that the likelihood of finding a repelling LCS at a given location is positively correlated with the mean strain rate injected at that point and negatively correlated with the mean speed, and that it is not correlated with the vorticity. We also find that mean time between successive LCSs appearing at a fixed location is related to the structure of the forcing field. Finally, we demonstrate a surprising difference in our results between LCSs computed forward and backwards in time, with forward-time (repelling) LCSs showing much more correlation with the forcing than backwards-time (attracting) LCSs.

Lagrangian Perturbation Approach to the Formation of Large-scale Structure

NASA Astrophysics Data System (ADS)

Buchert, Thomas

The present lecture notes address three columns on which the Lagrangian perturbation approach to cosmological dynamics is based: 1. the formulation of a Lagrangian theory of self-gravitating flows in which the dynamics is described in terms of a single field variable; 2. the procedure, how to obtain the dynamics of Eulerian fields from the Lagrangian picture, and 3. a precise definition of a Newtonian cosmology framework in which Lagrangian perturbation solutions can be studied. While the first is a discussion of the basic equations obtained by transforming the Eulerian evolution and field equations to the Lagrangian picture, the second exemplifies how the Lagrangian theory determines the evolution of Eulerian fields including kinematical variables like expansion, vorticity, as well as the shear and tidal tensors. The third column is based on a specification of initial and boundary conditions, and in particular on the identification of the average flow of an inhomogeneous cosmology with a `Hubble-flow'. Here, we also look at the limits of the Lagrangian perturbation approach as inferred from comparisons with N-body simulations and illustrate some striking properties of the solutions.

Mixed convection flow of viscoelastic fluid by a stretching cylinder with heat transfer.

PubMed

Hayat, Tasawar; Anwar, Muhammad Shoaib; Farooq, Muhammad; Alsaedi, Ahmad

2015-01-01

Flow of viscoelastic fluid due to an impermeable stretching cylinder is discussed. Effects of mixed convection and variable thermal conductivity are present. Thermal conductivity is taken temperature dependent. Nonlinear partial differential system is reduced into the nonlinear ordinary differential system. Resulting nonlinear system is computed for the convergent series solutions. Numerical values of skin friction coefficient and Nusselt number are computed and discussed. The results obtained with the current method are in agreement with previous studies using other methods as well as theoretical ideas. Physical interpretation reflecting the contribution of influential parameters in the present flow is presented. It is hoped that present study serves as a stimulus for modeling further stretching flows especially in polymeric and paper production processes.

Mixed Convection Flow of Viscoelastic Fluid by a Stretching Cylinder with Heat Transfer

PubMed Central

Hayat, Tasawar; Anwar, Muhammad Shoaib; Farooq, Muhammad; Alsaedi, Ahmad

2015-01-01

Flow of viscoelastic fluid due to an impermeable stretching cylinder is discussed. Effects of mixed convection and variable thermal conductivity are present. Thermal conductivity is taken temperature dependent. Nonlinear partial differential system is reduced into the nonlinear ordinary differential system. Resulting nonlinear system is computed for the convergent series solutions. Numerical values of skin friction coefficient and Nusselt number are computed and discussed. The results obtained with the current method are in agreement with previous studies using other methods as well as theoretical ideas. Physical interpretation reflecting the contribution of influential parameters in the present flow is presented. It is hoped that present study serves as a stimulus for modeling further stretching flows especially in polymeric and paper production processes. PMID:25775032

The piecewise-linear predictor-corrector code - A Lagrangian-remap method for astrophysical flows

NASA Technical Reports Server (NTRS)

Lufkin, Eric A.; Hawley, John F.

1993-01-01

We describe a time-explicit finite-difference algorithm for solving the nonlinear fluid equations. The method is similar to existing Eulerian schemes in its use of operator-splitting and artificial viscosity, except that we solve the Lagrangian equations of motion with a predictor-corrector and then remap onto a fixed Eulerian grid. The remap is formulated to eliminate errors associated with coordinate singularities, with a general prescription for remaps of arbitrary order. We perform a comprehensive series of tests on standard problems. Self-convergence tests show that the code has a second-order rate of convergence in smooth, two-dimensional flow, with pressure forces, gravity, and curvilinear geometry included. While not as accurate on idealized problems as high-order Riemann-solving schemes, the predictor-corrector Lagrangian-remap code has great flexibility for application to a variety of astrophysical problems.

Stochastic Simulation of Lagrangian Particle Transport in Turbulent Flows

NASA Astrophysics Data System (ADS)

Sun, Guangyuan

This dissertation presents the development and validation of the One Dimensional Turbulence (ODT) multiphase model in the Lagrangian reference frame. ODT is a stochastic model that captures the full range of length and time scales and provides statistical information on fine-scale turbulent-particle mixing and transport at low computational cost. The flow evolution is governed by a deterministic solution of the viscous processes and a stochastic representation of advection through stochastic domain mapping processes. The three algorithms for Lagrangian particle transport are presented within the context of the ODT approach. The Type-I and -C models consider the particle-eddy interaction as instantaneous and continuous change of the particle position and velocity, respectively. The Type-IC model combines the features of the Type-I and -C models. The models are applied to the multi-phase flows in the homogeneous decaying turbulence and turbulent round jet. Particle dispersion, dispersion coefficients, and velocity statistics are predicted and compared with experimental data. The models accurately reproduces the experimental data sets and capture particle inertial effects and trajectory crossing effect. A new adjustable particle parameter is introduced into the ODT model, and sensitivity analysis is performed to facilitate parameter estimation and selection. A novel algorithm of the two-way momentum coupling between the particle and carrier phases is developed in the ODT multiphase model. Momentum exchange between the phases is accounted for through particle source terms in the viscous diffusion. The source term is implemented in eddy events through a new kernel transformation and an iterative procedure is required for eddy selection. This model is applied to a particle-laden turbulent jet flow, and simulation results are compared with experimental measurements. The effect of particle addition on the velocities of the gas phase is investigated. The development of

Direct experimental visualization of the global Hamiltonian progression of two-dimensional Lagrangian flow topologies from integrable to chaotic state

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

Baskan, O.; Clercx, H. J. H; Speetjens, M. F. M.

Countless theoretical/numerical studies on transport and mixing in two-dimensional (2D) unsteady flows lean on the assumption that Hamiltonian mechanisms govern the Lagrangian dynamics of passive tracers. However, experimental studies specifically investigating said mechanisms are rare. Moreover, they typically concern local behavior in specific states (usually far away from the integrable state) and generally expose this indirectly by dye visualization. Laboratory experiments explicitly addressing the global Hamiltonian progression of the Lagrangian flow topology entirely from integrable to chaotic state, i.e., the fundamental route to efficient transport by chaotic advection, appear non-existent. This motivates our study on experimental visualization of this progressionmore » by direct measurement of Poincaré sections of passive tracer particles in a representative 2D time-periodic flow. This admits (i) accurate replication of the experimental initial conditions, facilitating true one-to-one comparison of simulated and measured behavior, and (ii) direct experimental investigation of the ensuing Lagrangian dynamics. The analysis reveals a close agreement between computations and observations and thus experimentally validates the full global Hamiltonian progression at a great level of detail.« less

An extended Lagrangian method

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

1993-01-01

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method', is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. Meanwhile, it also avoids the inaccuracy incurred due to geometry and variable interpolations used by the previous Lagrangian methods. The present method is general and capable of treating subsonic flows as well as supersonic flows. The method proposed in this paper is robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multidimensional discontinuities with a high level of accuracy, similar to that found in 1D problems.

Computing 3-D steady supersonic flow via a new Lagrangian approach

NASA Technical Reports Server (NTRS)

Loh, C. Y.; Liou, M.-S.

1993-01-01

The new Lagrangian method introduced by Loh and Hui (1990) is extended for 3-D steady supersonic flow computation. Details of the conservation form, the implementation of the local Riemann solver, and the Godunov and the high resolution TVD schemes are presented. The new approach is robust yet accurate, capable of handling complicated geometry and reactions between discontinuous waves. It keeps all the advantages claimed in the 2-D method of Loh and Hui, e.g., crisp resolution for a slip surface (contact discontinuity) and automatic grid generation along the stream.

Nonunitary Lagrangians and Unitary Non-Lagrangian Conformal Field Theories.

PubMed

Buican, Matthew; Laczko, Zoltan

2018-02-23

In various dimensions, we can sometimes compute observables of interacting conformal field theories (CFTs) that are connected to free theories via the renormalization group (RG) flow by computing protected quantities in the free theories. On the other hand, in two dimensions, it is often possible to algebraically construct observables of interacting CFTs using free fields without the need to explicitly construct an underlying RG flow. In this Letter, we begin to extend this idea to higher dimensions by showing that one can compute certain observables of an infinite set of unitary strongly interacting four-dimensional N=2 superconformal field theories (SCFTs) by performing simple calculations involving sets of nonunitary free four-dimensional hypermultiplets. These free fields are distant cousins of the Majorana fermion underlying the two-dimensional Ising model and are not obviously connected to our interacting theories via an RG flow. Rather surprisingly, this construction gives us Lagrangians for particular observables in certain subsectors of many "non-Lagrangian" SCFTs by sacrificing unitarity while preserving the full N=2 superconformal algebra. As a by-product, we find relations between characters in unitary and nonunitary affine Kac-Moody algebras. We conclude by commenting on possible generalizations of our construction.

Nonunitary Lagrangians and Unitary Non-Lagrangian Conformal Field Theories

NASA Astrophysics Data System (ADS)

Buican, Matthew; Laczko, Zoltan

2018-02-01

In various dimensions, we can sometimes compute observables of interacting conformal field theories (CFTs) that are connected to free theories via the renormalization group (RG) flow by computing protected quantities in the free theories. On the other hand, in two dimensions, it is often possible to algebraically construct observables of interacting CFTs using free fields without the need to explicitly construct an underlying RG flow. In this Letter, we begin to extend this idea to higher dimensions by showing that one can compute certain observables of an infinite set of unitary strongly interacting four-dimensional N =2 superconformal field theories (SCFTs) by performing simple calculations involving sets of nonunitary free four-dimensional hypermultiplets. These free fields are distant cousins of the Majorana fermion underlying the two-dimensional Ising model and are not obviously connected to our interacting theories via an RG flow. Rather surprisingly, this construction gives us Lagrangians for particular observables in certain subsectors of many "non-Lagrangian" SCFTs by sacrificing unitarity while preserving the full N =2 superconformal algebra. As a by-product, we find relations between characters in unitary and nonunitary affine Kac-Moody algebras. We conclude by commenting on possible generalizations of our construction.

Particle sedimentation in a sheared viscoelastic fluid

NASA Astrophysics Data System (ADS)

Murch, William L.; Krishnan, Sreenath; Shaqfeh, Eric S. G.; Iaccarino, Gianluca

2017-11-01

Particle suspensions are ubiquitous in engineered processes, biological systems, and natural settings. For an engineering application - whether the intent is to suspend and transport particles (e.g., in hydraulic fracturing fluids) or allow particles to sediment (e.g., in industrial separations processes) - understanding and prediction of the particle mobility is critical. This task is often made challenging by the complex nature of the fluid phase, for example, due to fluid viscoelasticity. In this talk, we focus on a fully 3D flow problem in a viscoelastic fluid: a settling particle with a shear flow applied in the plane perpendicular to gravity (referred to as orthogonal shear). Previously, it has been shown that an orthogonal shear flow can reduce the settling rate of particles in viscoelastic fluids. Using experiments and numerical simulations across a wide range of sedimentation and shear Weissenberg number, this talk will address the underlying physical mechanism responsible for the additional drag experienced by a rigid sphere settling in a confined viscoelastic fluid with orthogonal shear. We will then explore multiple particle effects, and discuss the implications and extensions of this work for particle suspensions. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE-114747 (WLM).

Hyperbolic reformulation of a 1D viscoelastic blood flow model and ADER finite volume schemes

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

Montecinos, Gino I.; Müller, Lucas O.; Toro, Eleuterio F.

2014-06-01

The applicability of ADER finite volume methods to solve hyperbolic balance laws with stiff source terms in the context of well-balanced and non-conservative schemes is extended to solve a one-dimensional blood flow model for viscoelastic vessels, reformulated as a hyperbolic system, via a relaxation time. A criterion for selecting relaxation times is found and an empirical convergence rate assessment is carried out to support this result. The proposed methodology is validated by applying it to a network of viscoelastic vessels for which experimental and numerical results are available. The agreement between the results obtained in the present paper and thosemore » available in the literature is satisfactory. Key features of the present formulation and numerical methodologies, such as accuracy, efficiency and robustness, are fully discussed in the paper.« less

Incomplete Mixing and Reactions - A Lagrangian Approach in a Pure Shear Flow

NASA Astrophysics Data System (ADS)

Paster, A.; Aquino, T.; Bolster, D.

2014-12-01

Incomplete mixing of reactive solutes is well known to slow down reaction rates relative to what would be expected from assuming perfect mixing. As reactions progress in a system and deplete reactant concentrations, initial fluctuations in the concentrations of reactions can be amplified relative to mean background concentrations and lead to spatial segregation of reactants. As the system evolves, in the absence of sufficient mixing, this segregation will increase, leading to a persistence of incomplete mixing that fundamentally changes the effective rate at which overall reactions will progress. On the other hand, non-uniform fluid flows are known to affect mixing between interacting solutes. Thus a natural question arises: Can non-uniform flows sufficiently enhance mixing to suppress incomplete mixing effects, and if so, under what conditions? In this work we address this question by considering one of the simplest possible flows, a laminar pure shear flow, which is known to significantly enhance mixing relative to diffusion alone. To study this system we adapt a novel Lagrangian particle-based random walk method, originally designed to simulate reactions in purely diffusive systems, to the case of advection and diffusion in a shear flow. To interpret the results we develop a semi-analytical solution, by proposing a closure approximation that aims to capture the effect of incomplete mixing. The results obtained via the Lagrangian model and the semi-analytical solutions consistently highlight that if shear effects in the system are not sufficiently strong, incomplete mixing effects initially similar to purely diffusive systems will occur, slowing down the overall reaction rate. Then, at some later time, dependent on the strength of the shear, the system will return to behaving as if it were well-mixed, but represented by a reduced effective reaction rate. If shear effects are sufficiently strong, the incomplete mixing regime never emerges and the system can behave

Incomplete Mixing and Reactions - A Lagrangian Approach in a Pure Shear Flow

NASA Astrophysics Data System (ADS)

Paster, Amir; Bolster, Diogo; Aquino, Tomas

2015-04-01

Incomplete mixing of reactive solutes is well known to slow down reaction rates relative to what would be expected from assuming perfect mixing. As reactions progress in a system and deplete reactant concentrations, initial fluctuations in the concentrations of reactions can be amplified relative to mean background concentrations and lead to spatial segregation of reactants. As the system evolves, in the absence of sufficient mixing, this segregation will increase, leading to a persistence of incomplete mixing that fundamentally changes the effective rate at which overall reactions will progress. On the other hand, nonuniform fluid flows are known to affect mixing between interacting solutes. Thus a natural question arises: Can non-uniform flows sufficiently enhance mixing to suppress incomplete mixing effects, and if so, under what conditions? In this work we address this question by considering one of the simplest possible flows, a laminar pure shear flow, which is known to significantly enhance mixing relative to diffusion alone. To study this system we adapt a novel Lagrangian particle-based random walk method, originally designed to simulate reactions in purely diffusive systems, to the case of advection and diffusion in a shear flow. To interpret the results we develop a semi-analytical solution, by proposing a closure approximation that aims to capture the effect of incomplete mixing. The results obtained via the Lagrangian model and the semi-analytical solutions consistently highlight that if shear effects in the system are not sufficiently strong, incomplete mixing effects initially similar to purely diffusive systems will occur, slowing down the overall reaction rate. Then, at some later time, dependent on the strength of the shear, the system will return to behaving as if it were well-mixed, but represented by a reduced effective reaction rate. If shear effects are sufficiently strong, the incomplete mixing regime never emerges and the system can behave

Correlating Lagrangian structures with forcing in two-dimensional flow

NASA Astrophysics Data System (ADS)

Ouellette, Nicholas T.; Hogg, Charlie A. R.; Liao, Yang

2016-01-01

Lagrangian coherent structures (LCSs) are the dominant transport barriers in unsteady, aperiodic flows, and their role in organizing mixing and transport has been well documented. However, nearly all that is known about LCSs has been gleaned from passive observations: they are computed in a post-processing step after a flow has been observed and used to understand why the mixing and transport proceeded as it did. In many applications, the ability instead to control the presence or location of LCSs via imposed forcing would be valuable. With this goal in mind, we study the relationship between LCSs and external forcing in an experimental quasi-two-dimensional weakly turbulent flow. We find that the likelihood of finding a repelling LCS at a given location is positively correlated with the mean strain rate injected at that point and negatively correlated with the mean speed, and that it is not correlated with the vorticity. We also find that mean time between successive LCSs appearing at a fixed location is related to the structure of the forcing field. Finally, we demonstrate a surprising difference in our results between LCSs computed forward and backward in time, with forward-time (repelling) LCSs showing much more correlation with the forcing than backwards-time (attracting) LCSs.

Swimming & Propulsion in Viscoelastic Media

NASA Astrophysics Data System (ADS)

Arratia, Paulo

2012-02-01

Many microorganisms have evolved within complex fluids, which include soil, intestinal fluid, and mucus. The material properties or rheology of such fluids can strongly affect an organism's swimming behavior. A major challenge is to understand the mechanism of propulsion in media that exhibit both solid- and fluid-like behavior, such as viscoelastic fluids. In this talk, we present experiments that explore the swimming behavior of biological organisms and artificial particles in viscoelastic media. The organism is the nematode Caenorhabditis elegans, a roundworm widely used for biological research that swims by generating traveling waves along its body. Overall, we find that fluid elasticity hinders self-propulsion compared to Newtonian fluids due to the enhanced resistance to flow near hyperbolic points for viscoelastic fluids. As fluid elasticity increases, the nematode's propulsion speed decreases. These results are consistent with recent theoretical models for undulating sheets and cylinders. In order to gain further understanding on propulsion in viscoelastic media, we perform experiments with simple reciprocal artificial `swimmers' (magnetic dumbbell particles) in polymeric and micellar solutions. We find that self-propulsion is possible in viscoelastic media even if the motion is reciprocal.

Using hyperbolic Lagrangian coherent structures to investigate vortices in bioinspired fluid flows

PubMed Central

Green, Melissa A.; Rowley, Clarence W.; Smits, Alexander J.

2010-01-01

We use direct Lyapunov exponents to identify Lagrangian coherent structures (LCSs) in a bioinspired fluid flow: the wakes of rigid pitching panels with a trapezoidal planform geometry chosen to model idealized fish caudal fins. When compared with commonly used Eulerian criteria, the Lagrangian method has previously exhibited the ability to define structure boundaries without relying on a preselected threshold. In addition, qualitative changes in the LCS have previously been shown to correspond to physical changes in the vortex structure. For this paper, digital particle image velocimetry experiments were performed to obtain the time-resolved velocity fields for Strouhal numbers of 0.17 and 0.27. A classic reverse von Kármán vortex street pattern was observed along the midspan of the near wake at low Strouhal number, but at higher Strouhal number the complexity of the wake increased downstream of the trailing edge. The spanwise vortices spread transversely across the wake and lose coherence, and this event was shown to correspond to a qualitative change in the LCS at the same time and location. PMID:20370300

Modeling of Quasi-Four-Phase Flow in Continuous Casting Mold Using Hybrid Eulerian and Lagrangian Approach

NASA Astrophysics Data System (ADS)

Liu, Zhongqiu; Sun, Zhenbang; Li, Baokuan

2017-04-01

Lagrangian tracking model combined with Eulerian multi-phase model is employed to predict the time-dependent argon-steel-slag-air quasi-four-phase flow inside a slab continuous casting mold. The Eulerian approach is used for the description of three phases (molten steel, liquid slag, and air at the top of liquid slag layer). The dispersed argon bubble injected from the SEN is treated in the Lagrangian way. The complex interfacial momentum transfers between various phases are considered. Validation is supported by the measurement data of cold model experiments and industrial practice. Close agreements were achieved for the gas volume fraction, liquid flow pattern, level fluctuation, and exposed slag eye phenomena. Many known phenomena and new predictions were successfully reproduced using this model. The vortex slag entrapment phenomenon at the slag-steel interface was obtained using this model, some small slag drops are sucked deep into the liquid pool of molten steel. Varying gas flow rates have a large effect on the steel flow pattern in the upper recirculation zone. Three typical flow patterns inside the mold with different argon gas flow rates have been obtained: double roll, three roll, and single roll. Effects of argon gas flow rate, casting speed, and slag layer thickness on the exposed slag eye and level fluctuation at the slag-steel interface were studied. A dimensionless value of H ave/ h was proposed to describe the time-averaged level fluctuation of slag-steel interface. The exposed slag eye near the SEN would be formed when the value of H ave/ h is larger than 0.4.

Spectral-clustering approach to Lagrangian vortex detection.

PubMed

Hadjighasem, Alireza; Karrasch, Daniel; Teramoto, Hiroshi; Haller, George

2016-06-01

One of the ubiquitous features of real-life turbulent flows is the existence and persistence of coherent vortices. Here we show that such coherent vortices can be extracted as clusters of Lagrangian trajectories. We carry out the clustering on a weighted graph, with the weights measuring pairwise distances of fluid trajectories in the extended phase space of positions and time. We then extract coherent vortices from the graph using tools from spectral graph theory. Our method locates all coherent vortices in the flow simultaneously, thereby showing high potential for automated vortex tracking. We illustrate the performance of this technique by identifying coherent Lagrangian vortices in several two- and three-dimensional flows.

Ligament Mediated Fragmentation of Viscoelastic Liquids

NASA Astrophysics Data System (ADS)

Keshavarz, Bavand; Houze, Eric C.; Moore, John R.; Koerner, Michael R.; McKinley, Gareth H.

2016-10-01

The breakup and atomization of complex fluids can be markedly different than the analogous processes in a simple Newtonian fluid. Atomization of paint, combustion of fuels containing antimisting agents, as well as physiological processes such as sneezing are common examples in which the atomized liquid contains synthetic or biological macromolecules that result in viscoelastic fluid characteristics. Here, we investigate the ligament-mediated fragmentation dynamics of viscoelastic fluids in three different canonical flows. The size distributions measured in each viscoelastic fragmentation process show a systematic broadening from the Newtonian solvent. In each case, the droplet sizes are well described by Gamma distributions which correspond to a fragmentation-coalescence scenario. We use a prototypical axial step strain experiment together with high-speed video imaging to show that this broadening results from the pronounced change in the corrugated shape of viscoelastic ligaments as they separate from the liquid core. These corrugations saturate in amplitude and the measured distributions for viscoelastic liquids in each process are given by a universal probability density function, corresponding to a Gamma distribution with nmin=4 . The breadth of this size distribution for viscoelastic filaments is shown to be constrained by a geometrical limit which can not be exceeded in ligament-mediated fragmentation phenomena.

Ligament Mediated Fragmentation of Viscoelastic Liquids.

PubMed

Keshavarz, Bavand; Houze, Eric C; Moore, John R; Koerner, Michael R; McKinley, Gareth H

2016-10-07

The breakup and atomization of complex fluids can be markedly different than the analogous processes in a simple Newtonian fluid. Atomization of paint, combustion of fuels containing antimisting agents, as well as physiological processes such as sneezing are common examples in which the atomized liquid contains synthetic or biological macromolecules that result in viscoelastic fluid characteristics. Here, we investigate the ligament-mediated fragmentation dynamics of viscoelastic fluids in three different canonical flows. The size distributions measured in each viscoelastic fragmentation process show a systematic broadening from the Newtonian solvent. In each case, the droplet sizes are well described by Gamma distributions which correspond to a fragmentation-coalescence scenario. We use a prototypical axial step strain experiment together with high-speed video imaging to show that this broadening results from the pronounced change in the corrugated shape of viscoelastic ligaments as they separate from the liquid core. These corrugations saturate in amplitude and the measured distributions for viscoelastic liquids in each process are given by a universal probability density function, corresponding to a Gamma distribution with n_{min}=4. The breadth of this size distribution for viscoelastic filaments is shown to be constrained by a geometrical limit which can not be exceeded in ligament-mediated fragmentation phenomena.

Guidelines for the formulation of Lagrangian stochastic models for particle simulations of single-phase and dispersed two-phase turbulent flows

NASA Astrophysics Data System (ADS)

Minier, Jean-Pierre; Chibbaro, Sergio; Pope, Stephen B.

2014-11-01

In this paper, we establish a set of criteria which are applied to discuss various formulations under which Lagrangian stochastic models can be found. These models are used for the simulation of fluid particles in single-phase turbulence as well as for the fluid seen by discrete particles in dispersed turbulent two-phase flows. The purpose of the present work is to provide guidelines, useful for experts and non-experts alike, which are shown to be helpful to clarify issues related to the form of Lagrangian stochastic models. A central issue is to put forward reliable requirements which must be met by Lagrangian stochastic models and a new element brought by the present analysis is to address the single- and two-phase flow situations from a unified point of view. For that purpose, we consider first the single-phase flow case and check whether models are fully consistent with the structure of the Reynolds-stress models. In the two-phase flow situation, coming up with clear-cut criteria is more difficult and the present choice is to require that the single-phase situation be well-retrieved in the fluid-limit case, elementary predictive abilities be respected and that some simple statistical features of homogeneous fluid turbulence be correctly reproduced. This analysis does not address the question of the relative predictive capacities of different models but concentrates on their formulation since advantages and disadvantages of different formulations are not always clear. Indeed, hidden in the changes from one structure to another are some possible pitfalls which can lead to flaws in the construction of practical models and to physically unsound numerical calculations. A first interest of the present approach is illustrated by considering some models proposed in the literature and by showing that these criteria help to assess whether these Lagrangian stochastic models can be regarded as acceptable descriptions. A second interest is to indicate how future

Scale-by-scale contributions to Lagrangian particle acceleration

NASA Astrophysics Data System (ADS)

Lalescu, Cristian C.; Wilczek, Michael

2017-11-01

Fluctuations on a wide range of scales in both space and time are characteristic of turbulence. Lagrangian particles, advected by the flow, probe these fluctuations along their trajectories. In an effort to isolate the influence of the different scales on Lagrangian statistics, we employ direct numerical simulations (DNS) combined with a filtering approach. Specifically, we study the acceleration statistics of tracers advected in filtered fields to characterize the smallest temporal scales of the flow. Emphasis is put on the acceleration variance as a function of filter scale, along with the scaling properties of the relevant terms of the Navier-Stokes equations. We furthermore discuss scaling ranges for higher-order moments of the tracer acceleration, as well as the influence of the choice of filter on the results. Starting from the Lagrangian tracer acceleration as the short time limit of the Lagrangian velocity increment, we also quantify the influence of filtering on Lagrangian intermittency. Our work complements existing experimental results on intermittency and accelerations of finite-sized, neutrally-buoyant particles: for the passive tracers used in our DNS, feedback effects are neglected such that the spatial averaging effect is cleanly isolated.

Methods to Approach Velocity Data Reduction and Their Effects on Conformation Statistics in Viscoelastic Turbulent Channel Flows

NASA Astrophysics Data System (ADS)

Samanta, Gaurab; Beris, Antony; Handler, Robert; Housiadas, Kostas

2009-03-01

Karhunen-Loeve (KL) analysis of DNS data of viscoelastic turbulent channel flows helps us to reveal more information on the time-dependent dynamics of viscoelastic modification of turbulence [Samanta et. al., J. Turbulence (in press), 2008]. A selected set of KL modes can be used for a data reduction modeling of these flows. However, it is pertinent that verification be done against established DNS results. For this purpose, we did comparisons of velocity and conformations statistics and probability density functions (PDFs) of relevant quantities obtained from DNS and reconstructed fields using selected KL modes and time-dependent coefficients. While the velocity statistics show good agreement between results from DNS and KL reconstructions even with just hundreds of KL modes, tens of thousands of KL modes are required to adequately capture the trace of polymer conformation resulting from DNS. New modifications to KL method have therefore been attempted to account for the differences in conformation statistics. The applicability and impact of these new modified KL methods will be discussed in the perspective of data reduction modeling.

Mode decomposition and Lagrangian structures of the flow dynamics in orbitally shaken bioreactors

NASA Astrophysics Data System (ADS)

Weheliye, Weheliye Hashi; Cagney, Neil; Rodriguez, Gregorio; Micheletti, Martina; Ducci, Andrea

2018-03-01

In this study, two mode decomposition techniques were applied and compared to assess the flow dynamics in an orbital shaken bioreactor (OSB) of cylindrical geometry and flat bottom: proper orthogonal decomposition and dynamic mode decomposition. Particle Image Velocimetry (PIV) experiments were carried out for different operating conditions including fluid height, h, and shaker rotational speed, N. A detailed flow analysis is provided for conditions when the fluid and vessel motions are in-phase (Fr = 0.23) and out-of-phase (Fr = 0.47). PIV measurements in vertical and horizontal planes were combined to reconstruct low order models of the full 3D flow and to determine its Finite-Time Lyapunov Exponent (FTLE) within OSBs. The combined results from the mode decomposition and the FTLE fields provide a useful insight into the flow dynamics and Lagrangian coherent structures in OSBs and offer a valuable tool to optimise bioprocess design in terms of mixing and cell suspension.

Focusing and alignment of erythrocytes in a viscoelastic medium

NASA Astrophysics Data System (ADS)

Go, Taesik; Byeon, Hyeokjun; Lee, Sang Joon

2017-01-01

Viscoelastic fluid flow-induced cross-streamline migration has recently received considerable attention because this process provides simple focusing and alignment over a wide range of flow rates. The lateral migration of particles depends on the channel geometry and physicochemical properties of particles. In this study, digital in-line holographic microscopy (DIHM) is employed to investigate the lateral migration of human erythrocytes induced by viscoelastic fluid flow in a rectangular microchannel. DIHM provides 3D spatial distributions of particles and information on particle orientation in the microchannel. The elastic forces generated in the pressure-driven flows of a viscoelastic fluid push suspended particles away from the walls and enforce erythrocytes to have a fixed orientation. Blood cell deformability influences the lateral focusing and fixed orientation in the microchannel. Different from rigid spheres and hardened erythrocytes, deformable normal erythrocytes disperse from the channel center plane, as the flow rate increases. Furthermore, normal erythrocytes have a higher angle of inclination than hardened erythrocytes in the region near the side-walls of the channel. These results may guide the label-free diagnosis of hematological diseases caused by abnormal erythrocyte deformability.

Getting Things Sorted With Lagrangian Coherent Structures

NASA Astrophysics Data System (ADS)

Atis, Severine; Peacock, Thomas; Environmental Dynamics Laboratory Team

2014-11-01

The dispersion of a tracer in a fluid flow is influenced by the Lagrangian motion of fluid elements. Even in laminar regimes, the irregular chaotic behavior of a fluid flow can lead to effective stirring that rapidly redistributes a tracer throughout the domain. For flows with arbitrary time-dependence, the modern approach of Lagrangian Coherent Structures (LCSs) provide a method for identifying the key material lines that organize flow transport. When the advected tracer particles possess a finite size and nontrivial shape, however, their dynamics can differ markedly from passive tracers, thus affecting the dispersion phenomena. We present details of numerical simulations and laboratory experiments that investigate the behavior of finite size particles in 2-dimensional chaotic flows. We show that the shape and the size of the particles alter the underlying LCSs, facilitating segregation between tracers of different shape in the same flow field.

Flow of viscoelastic fluids around a sharp microfluidic bend: Role of wormlike micellar structure

NASA Astrophysics Data System (ADS)

Hwang, Margaret Y.; Mohammadigoushki, Hadi; Muller, Susan J.

2017-04-01

We examine the flow and instabilities of three viscoelastic fluids—a semidilute aqueous solution of polyethylene oxide (PEO) and two wormlike micellar solutions of cetylpyridinium chloride and sodium salicylate—around a microfluidic 90∘ bend, in which shear deformation and streamline curvature dominate. Similar to results reported by Gulati et al. [S. Gulati et al., Phys. Rev. E 78, 036314 (2008), 10.1103/PhysRevE.78.036314; S. Gulati et al., J. Rheol. 54, 375 (2010), 10.1122/1.3308643] for PEO solutions, we report a critical Weissenberg number (Wi) for the onset of lip vortex formation upstream of the corner. However, the decreased aspect ratio (channel depth to width) results in a slightly higher critical Wi and a vortex that grows more slowly. We consider wormlike micellar solutions of two salt to surfactant concentration ratios R =0.55 and R =0.79 . At R =0.55 , the wormlike micelles are linear and exhibit strong viscoelastic behavior, but at R =0.79 , the wormlike micelles become branched and exhibit shear-banding behavior. Microfluidic experiments on the R =0.55 solution reveal two flow transitions. The first transition, at Wi =6 , is characterized by the formation of a stationary lip vortex upstream of the bend; at the second transition, at Wi =20 , the vortex fluctuates in time and changes size. The R =0.79 solution also exhibits two transitions. The first transition at Wi =4 is characterized by the appearance of two intermittent vortices, one at the lip and one at the far outside corner. Increasing the flow rate to Wi >160 results in a transition to a second unstable regime, where there is only a lip vortex that fluctuates in size. The difference in flow transitions in PEO and wormlike micellar solutions presumably arises from the additional contribution of wormlike micellar breakage and reformation under shear. The flow transitions in wormlike micellar solutions are also significantly affected by chain branching.

3-D Lagrangian-based investigations of the time-dependent cloud cavitating flows around a Clark-Y hydrofoil with special emphasis on shedding process analysis

NASA Astrophysics Data System (ADS)

Cheng, Huai-yu; Long, Xin-ping; Ji, Bin; Liu, Qi; Bai, Xiao-rui

2018-02-01

In the present paper, the unsteady cavitating flow around a 3-D Clark-Y hydrofoil is numerically investigated with the filter-based density correction model (FBDCM), a turbulence model and the Zwart-Gerber-Belamri (ZGB) cavitation model. A reasonable agreement is obtained between the numerical and experimental results. To study the complex flow structures more straightforwardly, a 3-D Lagrangian technology is developed, which can provide the particle tracks and the 3-D Lagrangian coherent structures (LCSs). Combined with the traditional methods based on the Eulerian viewpoint, this technology is used to analyze the attached cavity evolution and the re-entrant jet behavior in detail. At stage I, the collapse of the previous shedding cavity and the growth of a new attached cavity, the significant influence of the collapse both on the suction and pressure sides are captured quite well by the 3-D LCSs, which is underestimated by the traditional methods like the iso-surface of Q-criteria. As a kind of special LCSs, the arching LCSs are observed in the wake, induced by the counter-rotating vortexes. At stage II, with the development of the re-entrant jet, the influence of the cavitation on the pressure side is still not negligible. And with this 3-D Lagrangian technology, the tracks of the re-entrant jet are visualized clearly, moving from the trailing edge to the leading edge. Finally, at stage III, the re-entrant jet collides with the mainstream and finally induces the shedding. The cavitation evolution and the re-entrant jet movement in the whole cycle are well visualized with the 3-D Lagrangian technology. Moreover, the comparison between the LCSs obtained with 2-D and 3-D Lagrangian technologies indicates the advantages of the latter. It is demonstrated that the 3-D Lagrangian technology is a promising tool in the investigation of complex cavitating flows.

A non-conventional discontinuous Lagrangian for viscous flow

PubMed Central

Marner, F.

2017-01-01

Drawing an analogy with quantum mechanics, a new Lagrangian is proposed for a variational formulation of the Navier–Stokes equations which to-date has remained elusive. A key feature is that the resulting Lagrangian is discontinuous in nature, posing additional challenges apropos the mathematical treatment of the related variational problem, all of which are resolvable. In addition to extending Lagrange's formalism to problems involving discontinuous behaviour, it is demonstrated that the associated equations of motion can self-consistently be interpreted within the framework of thermodynamics beyond local equilibrium, with the limiting case recovering the classical Navier–Stokes equations. Perspectives for applying the new formalism to discontinuous physical phenomena such as phase and grain boundaries, shock waves and flame fronts are provided. PMID:28386415

A non-conventional discontinuous Lagrangian for viscous flow.

PubMed

Scholle, M; Marner, F

2017-02-01

Drawing an analogy with quantum mechanics, a new Lagrangian is proposed for a variational formulation of the Navier-Stokes equations which to-date has remained elusive. A key feature is that the resulting Lagrangian is discontinuous in nature, posing additional challenges apropos the mathematical treatment of the related variational problem, all of which are resolvable. In addition to extending Lagrange's formalism to problems involving discontinuous behaviour, it is demonstrated that the associated equations of motion can self-consistently be interpreted within the framework of thermodynamics beyond local equilibrium, with the limiting case recovering the classical Navier-Stokes equations. Perspectives for applying the new formalism to discontinuous physical phenomena such as phase and grain boundaries, shock waves and flame fronts are provided.

Electroosmosis of viscoelastic fluids over charge modulated surfaces in narrow confinements

NASA Astrophysics Data System (ADS)

Ghosh, Uddipta; Chakraborty, Suman

2015-06-01

In the present work, we attempt to analyze the electroosmotic flow of a viscoelastic fluid, following quasi-linear constitutive behavior, over charge modulated surfaces in narrow confinements. We obtain analytical solutions for the flow field for thin electrical double layer (EDL) limit through asymptotic analysis for small Deborah numbers. We show that a combination of matched and regular asymptotic expansion is needed for the thin EDL limit. We subsequently determine the modified Smoluchowski slip velocity for viscoelastic fluids and show that the quasi-linear nature of the constitutive behavior adds to the periodicity of the flow. We also obtain the net throughput in the channel and demonstrate its relative decrement as compared to that of a Newtonian fluid. Our results may have potential implications towards augmenting microfluidic mixing by exploiting electrokinetic transport of viscoelastic fluids over charge modulated surfaces.

A Comparative Study for Flow of Viscoelastic Fluids with Cattaneo-Christov Heat Flux.

PubMed

Hayat, Tasawar; Muhammad, Taseer; Alsaedi, Ahmed; Mustafa, Meraj

2016-01-01

This article examines the impact of Cattaneo-Christov heat flux in flows of viscoelastic fluids. Flow is generated by a linear stretching sheet. Influence of thermal relaxation time in the considered heat flux is seen. Mathematical formulation is presented for the boundary layer approach. Suitable transformations lead to a nonlinear differential system. Convergent series solutions of velocity and temperature are achieved. Impacts of various influential parameters on the velocity and temperature are sketched and discussed. Numerical computations are also performed for the skin friction coefficient and heat transfer rate. Our findings reveal that the temperature profile has an inverse relationship with the thermal relaxation parameter and the Prandtl number. Further the temperature profile and thermal boundary layer thickness are lower for Cattaneo-Christov heat flux model in comparison to the classical Fourier's law of heat conduction.

Two-dimensional Lagrangian simulation of suspended sediment

USGS Publications Warehouse

Schoellhamer, David H.

1988-01-01

A two-dimensional laterally averaged model for suspended sediment transport in steady gradually varied flow that is based on the Lagrangian reference frame is presented. The layered Lagrangian transport model (LLTM) for suspended sediment performs laterally averaged concentration. The elevations of nearly horizontal streamlines and the simulation time step are selected to optimize model stability and efficiency. The computational elements are parcels of water that are moved along the streamlines in the Lagrangian sense and are mixed with neighboring parcels. Three applications show that the LLTM can accurately simulate theoretical and empirical nonequilibrium suspended sediment distributions and slug injections of suspended sediment in a laboratory flume.

Experimental analysis on viscoelasticity-induced migration of RBCs using digital holographic microscopy

NASA Astrophysics Data System (ADS)

Go, Taesik; Byeon, Hyeokjun; Lee, Sang Joon

2016-11-01

Migration of particles in viscoelastic fluids has recently received large attention, because the generated elastic forces in viscoelastic fluids give rise to a simple focusing pattern over a wide range of flow rates. In this study, the vertical focusing and alignment of rigid spherical particles, normal and hardened RBCs in a viscoelastic fluid were experimentally investigated by employing a digital in-line holographic microscopy (DIHM). By the elastic forces, the three different particles are pushed away from the walls and concentrated in the midplane of the rectangular microchannel. Furthermore, most of both RBCs maintain face-on orientation in the microchannel. The effects of deformability of RBC on the viscoelasticity-induced migration and orientation in the channel were also examined. In contrary to non-deformable particles, normal RBCs are dispersed as flow rate increases. In the region near side wall of the microchannel, normal RBCs have edge-on orientation with a large angle of inclination, compared to hardened RBCs. These findings have a strong potential in the design of microfluidic devices for deformability-based separation of cells in viscoelastic fluid flows and label-free diagnoses of certain hematological diseases. This work was supported by the National Research Foundation of Korea (NRF) Grant funded by the Korea government (MSIP) (No. 2008-0061991).

Three dimensional Lagrangian structures in the Antarctic Polar Vortex.

NASA Astrophysics Data System (ADS)

Mancho, Ana M.; Garcia-Garrido, Victor J.; Curbelo, Jezabel; Niang, Coumba; Mechoso, Carlos R.; Wiggins, Stephen

2017-04-01

Dynamical systems theory has supported the description of transport processes in fluid dynamics. For understanding trajectory patterns in chaotic advection the geometrical approach by Poincaré seeks for spatial structures that separate regions corresponding to qualitatively different types of trajectories. These structures have been referred to as Lagrangian Coherent Structures (LCS), which typically in geophysical flows are well described under the approach of incompressible 2D flows. Different tools have been used to visualize LCS. In this presentation we use Lagrangian Descriptors [1,2,3,4] (function M) for visualizing 3D Lagrangian structures in the atmosphere, in particular in the Antarctic Polar Vortex. The function M is computed in a fully 3D incompressible flow obtained from data provided by the European Centre for Medium-Range Weather Forecast and it is represented in 2D surfaces. We discuss the findings during the final warming that took place in the spring of 1979 [5]. This research is supported by MINECO grant MTM2014-56392-R. Support is acknowledged also from CSIC grant COOPB20265, U.S. NSF grant AGS-1245069 and ONR grant No. N00014- 01-1-0769. C. Niang acknowledges Fundacion Mujeres por Africa and ICMAT Severo Ochoa project SEV-2011-0087 for financial support. [1] C. Mendoza, A. M. Mancho. The hidden geometry of ocean flows. Physical Review Letters 105 (2010), 3, 038501-1-038501-4. [2] A. M. Mancho, S. Wiggins, J. Curbelo, C. Mendoza. Lagrangian Descriptors: A Method for Revealing Phase Space Structures of General Time Dependent Dynamical Systems. Communications in Nonlinear Science and Numerical Simulation. 18 (2013) 3530-3557. [3] C. Lopesino, F. Balibrea-Iniesta, S. Wiggins and A. M. Mancho. Lagrangian descriptors for two dimensional, area preserving autonomous and nonautonomous maps. Communications in Nonlinear Science and Numerical Simulations, 27 (2015) (1-3), 40-51. [4] C. Lopesino, F. Balibrea-Iniesta, V. J. García-Garrido, S. Wiggins, and A

A review on the systematic formulation of 3-D multiparameter full waveform inversion in viscoelastic medium

NASA Astrophysics Data System (ADS)

Yang, Pengliang; Brossier, Romain; Métivier, Ludovic; Virieux, Jean

2016-10-01

In this paper, we study 3-D multiparameter full waveform inversion (FWI) in viscoelastic media based on the generalized Maxwell/Zener body including arbitrary number of attenuation mechanisms. We present a frequency-domain energy analysis to establish the stability condition of a full anisotropic viscoelastic system, according to zero-valued boundary condition and the elastic-viscoelastic correspondence principle: the real-valued stiffness matrix becomes a complex-valued one in Fourier domain when seismic attenuation is taken into account. We develop a least-squares optimization approach to linearly relate the quality factor with the anelastic coefficients by estimating a set of constants which are independent of the spatial coordinates, which supplies an explicit incorporation of the parameter Q in the general viscoelastic wave equation. By introducing the Lagrangian multipliers into the matrix expression of the wave equation with implicit time integration, we build a systematic formulation of multiparameter FWI for full anisotropic viscoelastic wave equation, while the equivalent form of the state and adjoint equation with explicit time integration is available to be resolved efficiently. In particular, this formulation lays the foundation for the inversion of the parameter Q in the time domain with full anisotropic viscoelastic properties. In the 3-D isotropic viscoelastic settings, the anelastic coefficients and the quality factors using bulk and shear moduli parametrization can be related to the counterparts using P and S velocity. Gradients with respect to any other parameter of interest can be found by chain rule. Pioneering numerical validations as well as the real applications of this most generic framework will be carried out to disclose the potential of viscoelastic FWI when adequate high-performance computing resources and the field data are available.

Euler-Lagrangian computation for estuarine hydrodynamics

USGS Publications Warehouse

Cheng, Ralph T.

1983-01-01

The transport of conservative and suspended matter in fluid flows is a phenomenon of Lagrangian nature because the process is usually convection dominant. Nearly all numerical investigations of such problems use an Eulerian formulation for the convenience that the computational grids are fixed in space and because the vast majority of field data are collected in an Eulerian reference frame. Several examples are given in this paper to illustrate a modeling approach which combines the advantages of both the Eulerian and Lagrangian computational techniques.

Influence of compressibility on the Lagrangian statistics of vorticity-strain-rate interactions.

PubMed

Danish, Mohammad; Sinha, Sawan Suman; Srinivasan, Balaji

2016-07-01

The objective of this study is to investigate the influence of compressibility on Lagrangian statistics of vorticity and strain-rate interactions. The Lagrangian statistics are extracted from "almost" time-continuous data sets of direct numerical simulations of compressible decaying isotropic turbulence by employing a cubic spline-based Lagrangian particle tracker. We study the influence of compressibility on Lagrangian statistics of alignment in terms of compressibility parameters-turbulent Mach number, normalized dilatation-rate, and flow topology. In comparison to incompressible turbulence, we observe that the presence of compressibility in a flow field weakens the alignment tendency of vorticity toward the largest strain-rate eigenvector. Based on the Lagrangian statistics of alignment conditioned on dilatation and topology, we find that the weakened tendency of alignment observed in compressible turbulence is because of a special group of fluid particles that have an initially negligible dilatation-rate and are associated with stable-focus-stretching topology.

Viscoelastic Taylor-Couette instability as analog of the magnetorotational instability.

PubMed

Bai, Yang; Crumeyrolle, Olivier; Mutabazi, Innocent

2015-09-01

A linear stability analysis and an experimental study of a viscoelastic Taylor-Couette flow corotating in the Keplerian ratio allow us to elucidate the analogy between the viscoelastic instability and the magnetorotational instability (MRI). A generalized Rayleigh criterion allows us to determine the potentially unstable zone to pure-elasticity-driven perturbations. Experiments with a viscoelastic polymer solution yield four modes: one pure-elasticity mode and three elastorotational instability (ERI) modes that represent the MRI-analog modes. The destabilization by the polymer viscosity is evidenced for the ERI modes.

Lagrangians for generalized Argyres-Douglas theories

NASA Astrophysics Data System (ADS)

Benvenuti, Sergio; Giacomelli, Simone

2017-10-01

We continue the study of Lagrangian descriptions of N=2 Argyres-Douglas theories. We use our recent interpretation in terms of sequential confinement to guess the Lagrangians of all the Argyres-Douglas models with Abelian three dimensional mirror. We find classes of four dimensional N=1 quivers that flow in the infrared to generalized Argyres-Douglas theories, such as the ( A k , A kN + N -1) models. We study in detail how the N=1 chiral rings map to the Coulomb and Higgs Branches of the N=2 CFT's. The three dimensional mirror RG flows are shown to land on the N=4 complete graph quivers. We also compactify to three dimensions the gauge theory dual to ( A 1, D 4), and find the expected Abelianization duality with N=4 SQED with 3 flavors.

Acoustic streaming: an arbitrary Lagrangian-Eulerian perspective.

PubMed

Nama, Nitesh; Huang, Tony Jun; Costanzo, Francesco

2017-08-25

We analyse acoustic streaming flows using an arbitrary Lagrangian Eulerian (ALE) perspective. The formulation stems from an explicit separation of time scales resulting in two subproblems: a first-order problem, formulated in terms of the fluid displacement at the fast scale, and a second-order problem, formulated in terms of the Lagrangian flow velocity at the slow time scale. Following a rigorous time-averaging procedure, the second-order problem is shown to be intrinsically steady, and with exact boundary conditions at the oscillating walls. Also, as the second-order problem is solved directly for the Lagrangian velocity, the formulation does not need to employ the notion of Stokes drift, or any associated post-processing, thus facilitating a direct comparison with experiments. Because the first-order problem is formulated in terms of the displacement field, our formulation is directly applicable to more complex fluid-structure interaction problems in microacoustofluidic devices. After the formulation's exposition, we present numerical results that illustrate the advantages of the formulation with respect to current approaches.

Four-way coupling of a three-dimensional debris flow solver to a Lagrangian Particle Simulation: method and first results

NASA Astrophysics Data System (ADS)

von Boetticher, Albrecht; Rickenmann, Dieter; McArdell, Brian; Kirchner, James W.

2017-04-01

Debris flows are dense flowing mixtures of water, clay, silt, sand and coarser particles. They are a common natural hazard in mountain regions and frequently cause severe damage. Modeling debris flows to design protection measures is still challenging due to the complex interactions within the inhomogeneous material mixture, and the sensitivity of the flow process to the channel geometry. The open-source, OpenFOAM-based finite-volume debris flow model debrisInterMixing (von Boetticher et al, 2016) defines rheology parameters based on the material properties of the debris flow mixture to reduce the number of free model parameters. As a simplification in this first model version, gravel was treated as a Coulomb-viscoplastic fluid, neglecting grain-to-grain collisions and the coupling between the coarser gravel grains and the interstitial fluid. Here we present an extension of that solver, accounting for the particle-to-particle and particle-to-boundary contacts with a Lagrangian Particle Simulation composed of spherical grains and a user-defined grain size distribution. The grain collisions of the Lagrangian particles add granular flow behavior to the finite-volume simulation of the continuous phases. The two-way coupling exchanges momentum between the phase-averaged flow in a finite volume cell, and among all individual particles contained in that cell, allowing the user to choose from a number of different drag models. The momentum exchange is implemented in the momentum equation and in the pressure equation (ensuring continuity) of the so-called PISO-loop, resulting in a stable 4-way coupling (particle-to-particle, particle-to-boundary, particle-to-fluid and fluid-to-particle) that represents the granular and viscous flow behavior of debris flow material. We will present simulations that illustrate the relative benefits and drawbacks of explicitly representing grain collisions, compared to the original debrisInterMixing solver.

Network-based study of Lagrangian transport and mixing

NASA Astrophysics Data System (ADS)

Padberg-Gehle, Kathrin; Schneide, Christiane

2017-10-01

Transport and mixing processes in fluid flows are crucially influenced by coherent structures and the characterization of these Lagrangian objects is a topic of intense current research. While established mathematical approaches such as variational methods or transfer-operator-based schemes require full knowledge of the flow field or at least high-resolution trajectory data, this information may not be available in applications. Recently, different computational methods have been proposed to identify coherent behavior in flows directly from Lagrangian trajectory data, that is, numerical or measured time series of particle positions in a fluid flow. In this context, spatio-temporal clustering algorithms have been proven to be very effective for the extraction of coherent sets from sparse and possibly incomplete trajectory data. Inspired by these recent approaches, we consider an unweighted, undirected network, where Lagrangian particle trajectories serve as network nodes. A link is established between two nodes if the respective trajectories come close to each other at least once in the course of time. Classical graph concepts are then employed to analyze the resulting network. In particular, local network measures such as the node degree, the average degree of neighboring nodes, and the clustering coefficient serve as indicators of highly mixing regions, whereas spectral graph partitioning schemes allow us to extract coherent sets. The proposed methodology is very fast to run and we demonstrate its applicability in two geophysical flows - the Bickley jet as well as the Antarctic stratospheric polar vortex.

Investigation of blood flow rheology using second-grade viscoelastic model (Phan-Thien-Tanner) within carotid artery.

PubMed

Ramiar, Abas; Larimi, Morsal Momenti; Ranjbar, Ali Akbar

2017-01-01

Hemodynamic factors, such as Wall Shear Stress (WSS), play a substantial role in arterial diseases. In the larger arteries, such as the carotid artery, interaction between the vessel wall and blood flow affects the distribution of hemodynamic factors. The fluid is considered to be non-Newtonian, whose flow is governed by the equation of a second-grade viscoelastic fluid and the effects of viscoelastic on blood flow in carotid artery is investigated. Pulsatile flow studies were carried out in a 3D model of carotid artery. The governing equations were solved using finite volume C++ based on open source code, OpenFOAM. To describe blood flow, conservation of mass and momentum, a constitutive relation of simplified Phan-Thien-Tanner (sPTT), and appropriate relations were used to explain shear thinning behavior. The first recirculation was observed at t = 0.2 s, in deceleration phase. In the acceleration phase from t = 0.3 s to t = 0.5 s, vortex and recirculation sizes in bulb regions in both ECA and ICA gradually increased. As is observed in the line graphs based on extracted data from ICA, at t = 0.2 s, τyy is the maximum amount of wall shear stress and τxy the minimum one. The maximum shear stress occurred in the inner side of the main branch (inner side of ICA and ECA) because the velocity of blood flow in the inner side of the bulb region was maximum due to the created recirculation zone in the opposite side in this area. The rheology of blood flow and shear stress in various important parts (the area that are in higher rates of WSS such as bifurcation region and the regions after bulb areas in both branches, Line1-4 in Fig. 7) were also analyzed. The investigation of velocity stream line, velocity profile and shear stress in various sections of carotid artery showed that the maximum shear stress occurred in acceleration phase and in the bifurcation region between ECA and ICA which is due to velocity gradients and changes in thinning behavior of blood and

Lattice Boltzmann simulation of viscoelastic flow past a confined free rotating cylinder

NASA Astrophysics Data System (ADS)

Xia, Yi; Zhang, Peijie; Lin, Jianzhong; Ku, Xiaoke; Nie, Deming

2018-05-01

To study the dynamics of rigid body immersed in viscoelastic fluid, an Oldroyd-B fluid flow past an eccentrically situated, free rotating cylinder in a two-dimensional (2D) channel is simulated by a novel lattice Boltzmann method. Two distribution functions are employed, one of which is aimed to solve Navier-Stokes equation and the other to the constitutive equation, respectively. The unified interpolation bounce-back scheme is adopted to treat the moving curved boundary of cylinder, and the novel Galilean invariant momentum exchange method is utilized to obtain the hydrodynamic force and torque exerted on the cylinder. Results show that the center-fixed cylinder rotates inversely in the direction where a cylinder immersed in Newtonian fluid do, which generates a centerline-oriented lift force according to Magnus effect. The cylinder’s eccentricity, flow inertia, fluid elasticity and viscosity would affect the rotation of cylinder in different ways. The cylinder rotates more rapidly when located farther away from the centerline, and slows down when it is too close to the wall. The rotation frequency decreases with increasing Reynolds number, and larger rotation frequency responds to larger Weissenberg number and smaller viscosity ratio, indicating that the fluid elasticity and low solvent viscosity accelerates the flow-induced rotation of cylinder.

A Lagrangian model of Copepod dynamics in turbulent flows

NASA Astrophysics Data System (ADS)

Ardeshiri, Hamidreza; Benkeddad, Ibtissem; Schmitt, Francois G.; Souissi, Sami; Toschi, Federico; Calzavarini, Enrico

2016-04-01

Planktonic copepods are small crustaceans that have the ability to swim by quick powerful jumps. Such an aptness is used to escape from high shear regions, which may be caused either by flow perturbations, produced by a large predator such as fish larave, or by the inherent highly turbulent dynamics of the ocean. Through a combined experimental and numerical study, we investigate the impact of jumping behaviour on the small-scale patchiness of copepods in a turbulent environment. Recorded velocity tracks of copepods displaying escape response jumps in still water are used to define and tune a Lagrangian Copepod (LC) model. The model is further employed to simulate the behaviour of thousands of copepods in a fully developed hydrodynamic turbulent flow obtained by direct numerical simulation of the Navier-Stokes equations. First, we show that the LC velocity statistics is in qualitative agreement with available experimental observations of copepods in turbulence. Second, we quantify the clustering of LC, via the fractal dimension D2. We show that D2 can be as low as 2.3, corresponding to local sheetlike aggregates, and that it critically depends on the shear-rate sensitivity of the proposed LC model. We further investigate the effect of jump intensity, jump orientation and geometrical aspect ratio of the copepods on the small-scale spatial distribution. Possible ecological implications of the observed clustering on encounter rates and mating success are discussed.

Users manual for a one-dimensional Lagrangian transport model

USGS Publications Warehouse

Schoellhamer, D.H.; Jobson, H.E.

1986-01-01

A Users Manual for the Lagrangian Transport Model (LTM) is presented. The LTM uses Lagrangian calculations that are based on a reference frame moving with the river flow. The Lagrangian reference frame eliminates the need to numerically solve the convective term of the convection-diffusion equation and provides significant numerical advantages over the more commonly used Eulerian reference frame. When properly applied, the LTM can simulate riverine transport and decay processes within the accuracy required by most water quality studies. The LTM is applicable to steady or unsteady one-dimensional unidirectional flows in fixed channels with tributary and lateral inflows. Application of the LTM is relatively simple and optional capabilities improve the model 's convenience. Appendices give file formats and three example LTM applications that include the incorporation of the QUAL II water quality model 's reaction kinetics into the LTM. (Author 's abstract)

Estimates of Lagrangian particle transport by wave groups: forward transport by Stokes drift and backward transport by the return flow

NASA Astrophysics Data System (ADS)

van den Bremer, Ton S.; Taylor, Paul H.

2014-11-01

Although the literature has examined Stokes drift, the net Lagrangian transport by particles due to of surface gravity waves, in great detail, the motion of fluid particles transported by surface gravity wave groups has received considerably less attention. In practice nevertheless, the wave field on the open sea often has a group-like structure. The motion of particles is different, as particles at sufficient depth are transported backwards by the Eulerian return current that was first described by Longuet-Higgins & Stewart (1962) and forms an inseparable counterpart of Stokes drift for wave groups ensuring the (irrotational) mass balance holds. We use WKB theory to study the variation of the Lagrangian transport by the return current with depth distinguishing two-dimensional seas, three-dimensional seas, infinite depth and finite depth. We then provide dimensional estimates of the net horizontal Lagrangian transport by the Stokes drift on the one hand and the return flow on the other hand for realistic sea states in all four cases. Finally we propose a simple scaling relationship for the transition depth: the depth above which Lagrangian particles are transported forwards by the Stokes drift and below which such particles are transported backwards by the return current.

The poro-viscoelastic properties of trabecular bone: a micro computed tomography-based finite element study.

PubMed

Sandino, Clara; McErlain, David D; Schipilow, John; Boyd, Steven K

2015-04-01

Bone is a porous structure with a solid phase that contains hydroxyapatite and collagen. Due to its composition, bone is often represented either as a poroelastic or as a viscoelastic material; however, the poro-viscoelastic formulation that allows integrating the effect of both the fluid flow and the collagen on the mechanical response of the tissue, has not been applied yet. The objective of this study was to develop a micro computed tomography (µCT)-based finite element (FE) model of trabecular bone that includes both the poroelastic and the viscoelastic nature of the tissue. Cubes of trabecular bone (N=25) from human distal tibia were scanned with µCT and stress relaxation experiments were conducted. The µCT images were the basis for sample specific FE models, and the stress relaxation experiments were simulated applying a poro-viscoelastic formulation. The model considers two scales of the tissue: the intertrabecular pore and the lacunar-canalicular pore scales. Independent viscoelastic and poroelastic models were also developed to determine their contribution to the poro-viscoelastic model. All the experiments exhibited a similar relaxation trend. The average reaction force before relaxation was 9.28 × 10(2)N (SD ± 5.11 × 10(2)N), and after relaxation was 4.69 × 10(2)N (SD ± 2.88 × 10(2)N). The slope of the regression line between the force before and after relaxation was 1.92 (R(2)=0.96). The poro-viscoelastic models captured 49% of the variability of the experimental data before relaxation and 33% after relaxation. The relaxation predicted with viscoelastic models was similar to the poro-viscoelastic ones; however, the poroelastic formulation underestimated the reaction force before relaxation. These data suggest that the contribution of viscoelasticity (fluid flow-independent mechanism) to the mechanical response of the tissue is significantly greater than the contribution of the poroelasticity (fluid flow-dependent mechanism). Copyright © 2015

Model and Comparative Study for Flow of Viscoelastic Nanofluids with Cattaneo-Christov Double Diffusion

PubMed Central

Hayat, Tasawar; Aziz, Arsalan; Muhammad, Taseer; Alsaedi, Ahmed

2017-01-01

Here two classes of viscoelastic fluids have been analyzed in the presence of Cattaneo-Christov double diffusion expressions of heat and mass transfer. A linearly stretched sheet has been used to create the flow. Thermal and concentration diffusions are characterized firstly by introducing Cattaneo-Christov fluxes. Novel features regarding Brownian motion and thermophoresis are retained. The conversion of nonlinear partial differential system to nonlinear ordinary differential system has been taken into place by using suitable transformations. The resulting nonlinear systems have been solved via convergent approach. Graphs have been sketched in order to investigate how the velocity, temperature and concentration profiles are affected by distinct physical flow parameters. Numerical values of skin friction coefficient and heat and mass transfer rates at the wall are also computed and discussed. Our observations demonstrate that the temperature and concentration fields are decreasing functions of thermal and concentration relaxation parameters. PMID:28046011

Actuation of flexoelectric membranes in viscoelastic fluids with applications to outer hair cells

PubMed Central

Herrera-Valencia, E. E.; Rey, Alejandro D.

2014-01-01

Liquid crystal flexoelectric actuation uses an imposed electric field to create membrane bending, and it is used by the outer hair cells (OHCs) located in the inner ear, whose role is to amplify sound through generation of mechanical power. Oscillations in the OHC membranes create periodic viscoelastic flows in the contacting fluid media. A key objective of this work on flexoelectric actuation relevant to OHCs is to find the relations and impact of the electromechanical properties of the membrane, the rheological properties of the viscoelastic media, and the frequency response of the generated mechanical power output. The model developed and used in this work is based on the integration of: (i) the flexoelectric membrane shape equation applied to a circular membrane attached to the inner surface of a circular capillary and (ii) the coupled capillary flow of contacting viscoelastic phases, such that the membrane flexoelectric oscillations drive periodic viscoelastic capillary flows, as in OHCs. By applying the Fourier transform formalism to the governing equation, analytical expressions for the transfer function associated with the curvature and electrical field and for the power dissipation of elastic storage energy were found. PMID:25332388

A Lagrangian dynamic subgrid-scale model turbulence

NASA Technical Reports Server (NTRS)

Meneveau, C.; Lund, T. S.; Cabot, W.

1994-01-01

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Variational Lagrangian data assimilation in open channel networks

NASA Astrophysics Data System (ADS)

Wu, Qingfang; Tinka, Andrew; Weekly, Kevin; Beard, Jonathan; Bayen, Alexandre M.

2015-04-01

This article presents a data assimilation method in a tidal system, where data from both Lagrangian drifters and Eulerian flow sensors were fused to estimate water velocity. The system is modeled by first-order, hyperbolic partial differential equations subject to periodic forcing. The estimation problem can then be formulated as the minimization of the difference between the observed variables and model outputs, and eventually provide the velocity and water stage of the hydrodynamic system. The governing equations are linearized and discretized using an implicit discretization scheme, resulting in linear equality constraints in the optimization program. Thus, the flow estimation can be formed as an optimization problem and efficiently solved. The effectiveness of the proposed method was substantiated by a large-scale field experiment in the Sacramento-San Joaquin River Delta in California. A fleet of 100 sensors developed at the University of California, Berkeley, were deployed in Walnut Grove, CA, to collect a set of Lagrangian data, a time series of positions as the sensors moved through the water. Measurements were also taken from Eulerian sensors in the region, provided by the United States Geological Survey. It is shown that the proposed method can effectively integrate Lagrangian and Eulerian measurement data, resulting in a suited estimation of the flow variables within the hydraulic system.

Flow stagnation volume and abdominal aortic aneurysm growth: Insights from patient-specific computational flow dynamics of Lagrangian-coherent structures.

PubMed

Joly, Florian; Soulez, Gilles; Garcia, Damien; Lessard, Simon; Kauffmann, Claude

2018-01-01

Abdominal aortic aneurysms (AAA) are localized, commonly-occurring dilations of the aorta. When equilibrium between blood pressure (loading) and wall mechanical resistance is lost, rupture ensues, and patient death follows, if not treated immediately. Experimental and numerical analyses of flow patterns in arteries show direct correlations between wall shear stress and wall mechano-adaptation with the development of zones prone to thrombus formation. For further insights into AAA flow topology/growth interaction, a workout of patient-specific computational flow dynamics (CFD) is proposed to compute finite-time Lyapunov exponents and extract Lagrangian-coherent structures (LCS). This computational model was first compared with 4-D phase-contrast magnetic resonance imaging (MRI) in 5 patients. To better understand the impact of flow topology and transport on AAA growth, hyperbolic, repelling LCS were computed in 1 patient during 8-year follow-up, including 9 volumetric morphologic AAA measures by computed tomography-angiography (CTA). LCS defined barriers to Lagrangian jet cores entering AAA. Domains enclosed between LCS and the aortic wall were considered to be stagnation zones. Their evolution was studied during AAA growth. Good correlation - 2-D cross-correlation coefficients of 0.65, 0.86 and 0.082 (min, max, SD) - was obtained between numerical simulations and 4-D MRI acquisitions in 6 specific cross-sections from 4 patients. In follow-up study, LCS divided AAA lumens into 3 dynamically-isolated zones: 2 stagnation volumes lying in dilated portions of the AAA, and circulating volume connecting the inlet to the outlet. The volume of each zone was tracked over time. Although circulating volume remained unchanged during 8-year follow-up, the AAA lumen and main stagnation zones grew significantly (8 cm 3 /year and 6 cm 3 /year, respectively). This study reveals that transient transport topology can be quantified in patient-specific AAA during disease progression

Micro-Macro Simulation of Viscoelastic Fluids in Three Dimensions

NASA Astrophysics Data System (ADS)

Rüttgers, Alexander; Griebel, Michael

2012-11-01

The development of the chemical industry resulted in various complex fluids that cannot be correctly described by classical fluid mechanics. For instance, this includes paint, engine oils with polymeric additives and toothpaste. We currently perform multiscale viscoelastic flow simulations for which we have coupled our three-dimensional Navier-Stokes solver NaSt3dGPF with the stochastic Brownian configuration field method on the micro-scale. In this method, we represent a viscoelastic fluid as a dumbbell system immersed in a three-dimensional Newtonian liquid which leads to a six-dimensional problem in space. The approach requires large computational resources and therefore depends on an efficient parallelisation strategy. Our flow solver is parallelised with a domain decomposition approach using MPI. It shows excellent scale-up results for up to 128 processors. In this talk, we present simulation results for viscoelastic fluids in square-square contractions due to their relevance for many engineering applications such as extrusion. Another aspect of the talk is the parallel implementation in NaSt3dGPF and the parallel scale-up and speed-up behaviour.

A Lagrangian Analysis of a Developing and Non-Developing Disturbance Observed During the PREDICT Experiment

DTIC Science & Technology

2012-12-03

paper provides an introduction of Lagrangian techniques for locating flow boundaries that encompass regions of recirculation in time- dependent flows...the low- to mid- level embryonic vortex from adverse conditions, while the 1The glossary on NOAA’s Hurricane Research Division’s web - site uses...wave or disturbance. This paper provides an introduction of Lagrangian techniques for locating flow boundaries that encompass regions of recirculation

Chaotic Lagrangian models for turbulent relative dispersion.

PubMed

Lacorata, Guglielmo; Vulpiani, Angelo

2017-04-01

A deterministic multiscale dynamical system is introduced and discussed as a prototype model for relative dispersion in stationary, homogeneous, and isotropic turbulence. Unlike stochastic diffusion models, here trajectory transport and mixing properties are entirely controlled by Lagrangian chaos. The anomalous "sweeping effect," a known drawback common to kinematic simulations, is removed through the use of quasi-Lagrangian coordinates. Lagrangian dispersion statistics of the model are accurately analyzed by computing the finite-scale Lyapunov exponent (FSLE), which is the optimal measure of the scaling properties of dispersion. FSLE scaling exponents provide a severe test to decide whether model simulations are in agreement with theoretical expectations and/or observation. The results of our numerical experiments cover a wide range of "Reynolds numbers" and show that chaotic deterministic flows can be very efficient, and numerically low-cost, models of turbulent trajectories in stationary, homogeneous, and isotropic conditions. The mathematics of the model is relatively simple, and, in a geophysical context, potential applications may regard small-scale parametrization issues in general circulation models, mixed layer, and/or boundary layer turbulence models as well as Lagrangian predictability studies.

Chaotic Lagrangian models for turbulent relative dispersion

NASA Astrophysics Data System (ADS)

Lacorata, Guglielmo; Vulpiani, Angelo

2017-04-01

A deterministic multiscale dynamical system is introduced and discussed as a prototype model for relative dispersion in stationary, homogeneous, and isotropic turbulence. Unlike stochastic diffusion models, here trajectory transport and mixing properties are entirely controlled by Lagrangian chaos. The anomalous "sweeping effect," a known drawback common to kinematic simulations, is removed through the use of quasi-Lagrangian coordinates. Lagrangian dispersion statistics of the model are accurately analyzed by computing the finite-scale Lyapunov exponent (FSLE), which is the optimal measure of the scaling properties of dispersion. FSLE scaling exponents provide a severe test to decide whether model simulations are in agreement with theoretical expectations and/or observation. The results of our numerical experiments cover a wide range of "Reynolds numbers" and show that chaotic deterministic flows can be very efficient, and numerically low-cost, models of turbulent trajectories in stationary, homogeneous, and isotropic conditions. The mathematics of the model is relatively simple, and, in a geophysical context, potential applications may regard small-scale parametrization issues in general circulation models, mixed layer, and/or boundary layer turbulence models as well as Lagrangian predictability studies.

Lagrangian methods of cosmic web classification

NASA Astrophysics Data System (ADS)

Fisher, J. D.; Faltenbacher, A.; Johnson, M. S. T.

2016-05-01

The cosmic web defines the large-scale distribution of matter we see in the Universe today. Classifying the cosmic web into voids, sheets, filaments and nodes allows one to explore structure formation and the role environmental factors have on halo and galaxy properties. While existing studies of cosmic web classification concentrate on grid-based methods, this work explores a Lagrangian approach where the V-web algorithm proposed by Hoffman et al. is implemented with techniques borrowed from smoothed particle hydrodynamics. The Lagrangian approach allows one to classify individual objects (e.g. particles or haloes) based on properties of their nearest neighbours in an adaptive manner. It can be applied directly to a halo sample which dramatically reduces computational cost and potentially allows an application of this classification scheme to observed galaxy samples. Finally, the Lagrangian nature admits a straightforward inclusion of the Hubble flow negating the necessity of a visually defined threshold value which is commonly employed by grid-based classification methods.

Viscoelastic tides: models for use in Celestial Mechanics

NASA Astrophysics Data System (ADS)

Ragazzo, C.; Ruiz, L. S.

2017-05-01

This paper contains equations for the motion of linear viscoelastic bodies interacting under gravity. The equations are fully three dimensional and allow for the integration of the spin, the orbit, and the deformation of each body. The goal is to present good models for the tidal forces that take into account the possibly different rheology of each body. The equations are obtained within a finite dimension Lagrangian framework with dissipation function. The main contribution is a procedure to associate to each spring-dashpot model, which defines the rheology of a body, a potential and a dissipation function for the body deformation variables. The theory is applied to the Earth (solid part plus oceans) and a comparison between model and observation of the following quantities is made: norm of the Love numbers, rate of tidal energy dissipation, Chandler period, and Earth-Moon distance increase.

An adaptive reconstruction for Lagrangian, direct-forcing, immersed-boundary methods

NASA Astrophysics Data System (ADS)

Posa, Antonio; Vanella, Marcos; Balaras, Elias

2017-12-01

Lagrangian, direct-forcing, immersed boundary (IB) methods have been receiving increased attention due to their robustness in complex fluid-structure interaction problems. They are very sensitive, however, on the selection of the Lagrangian grid, which is typically used to define a solid or flexible body immersed in a fluid flow. In the present work we propose a cost-efficient solution to this problem without compromising accuracy. Central to our approach is the use of isoparametric mapping to bridge the relative resolution requirements of Lagrangian IB, and Eulerian grids. With this approach, the density of surface Lagrangian markers, which is essential to properly enforce boundary conditions, is adapted dynamically based on the characteristics of the underlying Eulerian grid. The markers are not stored and the Lagrangian data-structure is not modified. The proposed scheme is implemented in the framework of a moving least squares reconstruction formulation, but it can be adapted to any Lagrangian, direct-forcing formulation. The accuracy and robustness of the approach is demonstrated in a variety of test cases of increasing complexity.

Pure axial flow of viscoelastic fluids in rectangular microchannels under combined effects of electro-osmosis and hydrodynamics

NASA Astrophysics Data System (ADS)

Reshadi, Milad; Saidi, Mohammad Hassan; Ebrahimi, Abbas

2018-02-01

This paper presents an analysis of the combined electro-osmotic and pressure-driven axial flows of viscoelastic fluids in a rectangular microchannel with arbitrary aspect ratios. The rheological behavior of the fluid is described by the complete form of Phan-Thien-Tanner (PTT) model with the Gordon-Schowalter convected derivative which covers the upper convected Maxwell, Johnson-Segalman and FENE-P models. Our numerical simulation is based on the computation of 2D Poisson-Boltzmann, Cauchy momentum and PTT constitutive equations. The solution of these governing nonlinear coupled set of equations is obtained by using the second-order central finite difference method in a non-uniform grid system and is verified against 1D analytical solution of the velocity profile with less than 0.06% relative error. Also, a parametric study is carried out to investigate the effect of channel aspect ratio (width to height), wall zeta potential and the Debye-Hückel parameter on 2D velocity profile, volumetric flow rate and the Poiseuille number in the mixed EO/PD flows of viscoelastic fluids with different Weissenberg numbers. Our results show that, for low channel aspect ratios, the previous 1D analytical models underestimate the velocity profile at the channel half-width centerline in the case of favorable pressure gradients and overestimate it in the case of adverse pressure gradients. The results reveal that the inapplicability of the Debye-Hückel approximation at high zeta potentials is more significant for higher Weissenberg number fluids. Also, it is found that, under the specified values of electrokinetic parameters, there is a threshold for velocity scale ratio in which the Poiseuille number is approximately independent of channel aspect ratio.

Guidelines for the formulation of Lagrangian stochastic models for particle simulations of single-phase and dispersed two-phase turbulent flows

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

Minier, Jean-Pierre, E-mail: Jean-Pierre.Minier@edf.fr; Chibbaro, Sergio; Pope, Stephen B.

In this paper, we establish a set of criteria which are applied to discuss various formulations under which Lagrangian stochastic models can be found. These models are used for the simulation of fluid particles in single-phase turbulence as well as for the fluid seen by discrete particles in dispersed turbulent two-phase flows. The purpose of the present work is to provide guidelines, useful for experts and non-experts alike, which are shown to be helpful to clarify issues related to the form of Lagrangian stochastic models. A central issue is to put forward reliable requirements which must be met by Lagrangianmore » stochastic models and a new element brought by the present analysis is to address the single- and two-phase flow situations from a unified point of view. For that purpose, we consider first the single-phase flow case and check whether models are fully consistent with the structure of the Reynolds-stress models. In the two-phase flow situation, coming up with clear-cut criteria is more difficult and the present choice is to require that the single-phase situation be well-retrieved in the fluid-limit case, elementary predictive abilities be respected and that some simple statistical features of homogeneous fluid turbulence be correctly reproduced. This analysis does not address the question of the relative predictive capacities of different models but concentrates on their formulation since advantages and disadvantages of different formulations are not always clear. Indeed, hidden in the changes from one structure to another are some possible pitfalls which can lead to flaws in the construction of practical models and to physically unsound numerical calculations. A first interest of the present approach is illustrated by considering some models proposed in the literature and by showing that these criteria help to assess whether these Lagrangian stochastic models can be regarded as acceptable descriptions. A second interest is to indicate how future

Postseismic viscoelastic surface deformation and stress. Part 1: Theoretical considerations, displacement and strain calculations

NASA Technical Reports Server (NTRS)

Cohen, S. C.

1979-01-01

A model of viscoelastic deformations associated with earthquakes is presented. A strike-slip fault is represented by a rectangular dislocation in a viscoelastic layer (lithosphere) lying over a viscoelastic half-space (asthenosphere). Deformations occur on three time scales. The initial response is governed by the instantaneous elastic properties of the earth. A slower response is associated with viscoelastic relaxation of the lithosphere and a yet slower response is due to viscoelastic relaxation of the asthenosphere. The major conceptual contribution is the inclusion of lithospheric viscoelastic properties into a dislocation model of earthquake related deformations and stresses. Numerical calculations using typical fault parameters reveal that the postseismic displacements and strains are small compared to the coseismic ones near the fault, but become significant further away. Moreover, the directional sense of the deformations attributable to the elastic response, the lithospheric viscoelastic softening, and the asthenospheric viscoelastic flow may differ and depend on location and model details. The results and theoretical arguments suggest that the stress changes accompanying lithospheric relaxation may also be in a different sense than and be larger than the strain changes.

An updated Lagrangian particle hydrodynamics (ULPH) for Newtonian fluids

NASA Astrophysics Data System (ADS)

Tu, Qingsong; Li, Shaofan

2017-11-01

In this work, we have developed an updated Lagrangian particle hydrodynamics (ULPH) for Newtonian fluid. Unlike the smoothed particle hydrodynamics, the non-local particle hydrodynamics formulation proposed here is consistent and convergence. Unlike the state-based peridynamics, the discrete particle dynamics proposed here has no internal material bond between particles, and it is not formulated with respect to initial or a fixed referential configuration. In specific, we have shown that (1) the non-local update Lagrangian particle hydrodynamics formulation converges to the conventional local fluid mechanics formulation; (2) the non-local updated Lagrangian particle hydrodynamics can capture arbitrary flow discontinuities without any changes in the formulation, and (3) the proposed non-local particle hydrodynamics is computationally efficient and robust.

Lagrangian particle method for compressible fluid dynamics

NASA Astrophysics Data System (ADS)

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

2018-06-01

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

A Godunov-like point-centered essentially Lagrangian hydrodynamic approach

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

Morgan, Nathaniel R.; Waltz, Jacob I.; Burton, Donald E.

We present an essentially Lagrangian hydrodynamic scheme suitable for modeling complex compressible flows on tetrahedron meshes. The scheme reduces to a purely Lagrangian approach when the flow is linear or if the mesh size is equal to zero; as a result, we use the term essentially Lagrangian for the proposed approach. The motivation for developing a hydrodynamic method for tetrahedron meshes is because tetrahedron meshes have some advantages over other mesh topologies. Notable advantages include reduced complexity in generating conformal meshes, reduced complexity in mesh reconnection, and preserving tetrahedron cells with automatic mesh refinement. A challenge, however, is tetrahedron meshesmore » do not correctly deform with a lower order (i.e. piecewise constant) staggered-grid hydrodynamic scheme (SGH) or with a cell-centered hydrodynamic (CCH) scheme. The SGH and CCH approaches calculate the strain via the tetrahedron, which can cause artificial stiffness on large deformation problems. To resolve the stiffness problem, we adopt the point-centered hydrodynamic approach (PCH) and calculate the evolution of the flow via an integration path around the node. The PCH approach stores the conserved variables (mass, momentum, and total energy) at the node. The evolution equations for momentum and total energy are discretized using an edge-based finite element (FE) approach with linear basis functions. A multidirectional Riemann-like problem is introduced at the center of the tetrahedron to account for discontinuities in the flow such as a shock. Conservation is enforced at each tetrahedron center. The multidimensional Riemann-like problem used here is based on Lagrangian CCH work [8, 19, 37, 38, 44] and recent Lagrangian SGH work [33-35, 39, 45]. In addition, an approximate 1D Riemann problem is solved on each face of the nodal control volume to advect mass, momentum, and total energy. The 1D Riemann problem produces fluxes [18] that remove a volume error in the PCH

A Godunov-like point-centered essentially Lagrangian hydrodynamic approach

DOE PAGES

Morgan, Nathaniel R.; Waltz, Jacob I.; Burton, Donald E.; ...

2014-10-28

We present an essentially Lagrangian hydrodynamic scheme suitable for modeling complex compressible flows on tetrahedron meshes. The scheme reduces to a purely Lagrangian approach when the flow is linear or if the mesh size is equal to zero; as a result, we use the term essentially Lagrangian for the proposed approach. The motivation for developing a hydrodynamic method for tetrahedron meshes is because tetrahedron meshes have some advantages over other mesh topologies. Notable advantages include reduced complexity in generating conformal meshes, reduced complexity in mesh reconnection, and preserving tetrahedron cells with automatic mesh refinement. A challenge, however, is tetrahedron meshesmore » do not correctly deform with a lower order (i.e. piecewise constant) staggered-grid hydrodynamic scheme (SGH) or with a cell-centered hydrodynamic (CCH) scheme. The SGH and CCH approaches calculate the strain via the tetrahedron, which can cause artificial stiffness on large deformation problems. To resolve the stiffness problem, we adopt the point-centered hydrodynamic approach (PCH) and calculate the evolution of the flow via an integration path around the node. The PCH approach stores the conserved variables (mass, momentum, and total energy) at the node. The evolution equations for momentum and total energy are discretized using an edge-based finite element (FE) approach with linear basis functions. A multidirectional Riemann-like problem is introduced at the center of the tetrahedron to account for discontinuities in the flow such as a shock. Conservation is enforced at each tetrahedron center. The multidimensional Riemann-like problem used here is based on Lagrangian CCH work [8, 19, 37, 38, 44] and recent Lagrangian SGH work [33-35, 39, 45]. In addition, an approximate 1D Riemann problem is solved on each face of the nodal control volume to advect mass, momentum, and total energy. The 1D Riemann problem produces fluxes [18] that remove a volume error in the PCH

Approximation of Viscoelastic Stresses from Newtonian Turbulent Kinematics

DTIC Science & Technology

1988-09-01

birefringence of polyethylene oxide solutions in a four roll mill. J.Poly.Sci.:Poly.Phys.Ed. 14, 1111-1119. Dandridge, A., Meeten , G.H., Layec-Raphalen, M.N...flows. Poly. Comm. 25, 144-146. Metzner, A.B., & Astarita, G . 1967 External flow of viscoelastic materials: fluid property restrictions on the use of...dumbbell model for dilute solutions. Rheol.Acta 23, 151-162. Philippoff, W. 1956 Flow-birefringence and stress. Nature 178 , 811-812. Ryskin, G . 1987a

Lateral migration of a viscoelastic drop in a Newtonian fluid in a shear flow near a wall

PubMed Central

Mukherjee, Swarnajay; Sarkar, Kausik

2014-01-01

Wall induced lateral migration of a viscoelastic (FENE-MCR) drop in a Newtonian fluid is investigated. Just like a Newtonian drop, a viscoelastic drop reaches a quasi-steady state where the lateral velocity only depends on the instantaneous distance from the wall. The drop migration velocity and the deformation scale inversely with the square and the cube of the distance from the wall, respectively. The migration velocity varies non-monotonically with increasing viscoelasticity (increasing Deborah number); initially increasing and then decreasing. An analytical explanation has been given of the effects by computing the migration velocity as arising from an image stresslet field due to the drop. The semi-analytical expression matches well with the simulated migration velocity away from the wall. It contains a viscoelastic stresslet component apart from those arising from interfacial tension and viscosity ratio. The migration dynamics is a result of the competition between the viscous (interfacial tension and viscosity ratio) and the viscoelastic effects. The viscoelastic stresslet contribution towards the migration velocity steadily increases. But the interfacial stresslet—arising purely from the drop shape—first increases and then decreases with rising Deborah number causing the migration velocity to be non-monotonic. The geometric effect of the interfacial stresslet is caused by a corresponding nonmonotonic variation of the drop inclination. High viscosity ratio is briefly considered to show that the drop viscoelasticity could stabilize a drop against breakup, and the increase in migration velocity due to viscoelasticity is larger compared to the viscosity-matched case. PMID:25378894

Pearling Instabilities of a Viscoelastic Thread

NASA Astrophysics Data System (ADS)

Deblais, A.; Velikov, K. P.; Bonn, D.

2018-05-01

Pearling instabilities of slender viscoelastic threads have received much attention, but remain incompletely understood. We study the instabilities in polymer solutions subject to uniaxial elongational flow. Two distinctly different instabilites are observed: beads on a string and blistering. The beads-on-a-string structure arises from a capillary instability whereas the blistering instability has a different origin: it is due to a coupling between stress and polymer concentration. By varying the temperature to change the solution properties we elucidate the interplay between flow and phase separation.

A comparative study on the breakup of Newtonian and viscoelastic liquid films

NASA Astrophysics Data System (ADS)

Qian, Lijuan; Song, Shaobo; Jiang, Lisha; Li, Xiaolu; Lin, Jianzhong

2018-05-01

The breakup of viscoelastic liquid films are investigated experimentally and analytically. The breakup phenomena of viscoelastic liquid film are recorded by the time resolved high speed camera. Video images reveal the difference behavior of liquid bubble breakup for Newtonian and viscoelastic liquid. For the Newtonian liquid, cylindrical ligaments are stretched into droplets with large distributions of drop size. For the viscoelastic liquid, the pinch-off point is located on the liquid connections to the nozzle and finally the main part of the ligament no longer elongates. Furthermore, a dispersion relation based on the stability analysis is involved to predict the ligament length and drop mean size after breakup for liquid film. The calculated ligament length is validated by the measured drop mean size at higher air-to-liquid mass flow ratio.

Comparing High-latitude Ionospheric and Thermospheric Lagrangian Coherent Structures

NASA Astrophysics Data System (ADS)

Wang, N.; Ramirez, U.; Flores, F.; Okic, D.; Datta-Barua, S.

2015-12-01

Lagrangian Coherent Structures (LCSs) are invisible boundaries in time varying flow fields that may be subject to mixing and turbulence. The LCS is defined by the local maxima of the finite time Lyapunov exponent (FTLE), a scalar field quantifying the degree of stretching of fluid elements over the flow domain. Although the thermosphere is dominated by neutral wind processes and the ionosphere is governed by plasma electrodynamics, we can compare the LCS in the two modeled flow fields to yield insight into transport and interaction processes in the high-latitude IT system. For obtaining thermospheric LCS, we use the Horizontal Wind Model 2014 (HWM14) [1] at a single altitude to generate the two-dimensional velocity field. The FTLE computation is applied to study the flow field of the neutral wind, and to visualize the forward-time Lagrangian Coherent Structures in the flow domain. The time-varying structures indicate a possible thermospheric LCS ridge in the auroral oval area. The results of a two-day run during a geomagnetically quiet period show that the structures are diurnally quasi-periodic, thus that solar radiation influences the neutral wind flow field. To find the LCS in the high-latitude ionospheric drifts, the Weimer 2001 [2] polar electric potential model and the International Geomagnetic Reference Field 11 [3] are used to compute the ExB drift flow field in ionosphere. As with the neutral winds, the Lagrangian Coherent Structures are obtained by applying the FTLE computation. The relationship between the thermospheric and ionospheric LCS is analyzed by comparing overlapping FTLE maps. Both a publicly available FTLE solver [4] and a custom-built FTLE computation are used and compared for validation [5]. Comparing the modeled IT LCSs on a quiet day with the modeled IT LCSs on a storm day indicates important factors on the structure and time evolution of the LCS.

An online-coupled NWP/ACT model with conserved Lagrangian levels

NASA Astrophysics Data System (ADS)

Sørensen, B.; Kaas, E.; Lauritzen, P. H.

2012-04-01

Numerical weather and climate modelling is under constant development. Semi-implicit semi-Lagrangian (SISL) models have proven to be numerically efficient in both short-range weather forecasts and climate models, due to the ability to use long time steps. Chemical/aerosol feedback mechanism are becoming more and more relevant in NWP as well as climate models, since the biogenic and anthropogenic emissions can have a direct effect on the dynamics and radiative properties of the atmosphere. To include chemical feedback mechanisms in the NWP models, on-line coupling is crucial. In 3D semi-Lagrangian schemes with quasi-Lagrangian vertical coordinates the Lagrangian levels are remapped to Eulerian model levels each time step. This remapping introduces an undesirable tendency to smooth sharp gradients and creates unphysical numerical diffusion in the vertical distribution. A semi-Lagrangian advection method is introduced, it combines an inherently mass conserving 2D semi-Lagrangian scheme, with a SISL scheme employing both hybrid vertical coordinates and a fully Lagrangian vertical coordinate. This minimizes the vertical diffusion and thus potentially improves the simulation of the vertical profiles of moisture, clouds, and chemical constituents. Since the Lagrangian levels suffer from traditional Lagrangian limitations caused by the convergence and divergence of the flow, remappings to the Eulerian model levels are generally still required - but this need only be applied after a number of time steps - unless dynamic remapping methods are used. For this several different remapping methods has been implemented. The combined scheme is mass conserving, consistent, and multi-tracer efficient.

Lagrangian and Eulerian statistics obtained from direct numerical simulations of homogeneous turbulence

NASA Technical Reports Server (NTRS)

Squires, Kyle D.; Eaton, John K.

1991-01-01

Direct numerical simulation is used to study dispersion in decaying isotropic turbulence and homogeneous shear flow. Both Lagrangian and Eulerian data are presented allowing direct comparison, but at fairly low Reynolds number. The quantities presented include properties of the dispersion tensor, isoprobability contours of particle displacement, Lagrangian and Eulerian velocity autocorrelations and time scale ratios, and the eddy diffusivity tensor. The Lagrangian time microscale is found to be consistently larger than the Eulerian microscale, presumably due to the advection of the small scales by the large scales in the Eulerian reference frame.

Deformation and relaxation of an incompressible viscoelastic body with surface viscoelasticity

NASA Astrophysics Data System (ADS)

Liu, Liping; Yu, Miao; Lin, Hao; Foty, Ramsey

2017-01-01

Measuring mechanical properties of cells or cell aggregates has proven to be an involved process due to their geometrical and structural complexity. Past measurements are based on material models that completely neglect the elasticity of either the surface membrane or the interior bulk. In this work, we consider general material models to account for both surface and bulk viscoelasticity. The boundary value problems are formulated for deformations and relaxations of a closed viscoelastic surface coupled with viscoelastic media inside and outside of the surface. The linearized surface elasticity models are derived for the constant surface tension model and the Helfrich-Canham bending model for coupling with the bulk viscoelasticity. For quasi-spherical surfaces, explicit solutions are obtained for the deformation, stress-strain and relaxation behaviors under a variety of loading conditions. These solutions can be applied to extract the intrinsic surface and bulk viscoelastic properties of biological cells or cell aggregates in the indentation, electro-deformation and relaxation experiments.

Transport induced by mean-eddy interaction: I. Theory, and relation to Lagrangian lobe dynamics

NASA Astrophysics Data System (ADS)

Ide, Kayo; Wiggins, Stephen

2015-02-01

In this paper we develop a method for the estimation of Transport Induced by the Mean-Eddy interaction (TIME) in two-dimensional unsteady flows. The method is based on the dynamical systems approach to fluid transport and can be viewed as a hybrid combination of Lagrangian and Eulerian methods. The (Eulerian) boundaries across which we consider (Lagrangian) transport are kinematically defined by appropriately chosen streamlines of the mean flow. By evaluating the impact of the mean-eddy interaction on transport, the TIME method can be used as a diagnostic tool for transport processes that occur during a specified time interval along a specified boundary segment. We introduce two types of TIME functions: one that quantifies the accumulation of flow properties and another that measures the displacement of the transport geometry. The spatial geometry of transport is described by the so-called pseudo-lobes, and temporal evolution of transport by their dynamics. In the case where the TIME functions are evaluated along a separatrix, the pseudo-lobes have a relationship to the lobes of Lagrangian transport theory. In fact, one of the TIME functions is identical to the Melnikov function that is used to measure the distance, at leading order in a small parameter, between the two invariant manifolds that define the Lagrangian lobes. We contrast the similarities and differences between the TIME and Lagrangian lobe dynamics in detail. An application of the TIME method is carried out for inter-gyre transport in the wind-driven oceanic circulation model and a comparison with the Lagrangian transport theory is made.

Some Lagrangians for systems without a Lagrangian

NASA Astrophysics Data System (ADS)

Nucci, M. C.; Leach, P. G. L.

2011-03-01

We demonstrate how to construct many different Lagrangians for two famous examples that were deemed by Douglas (1941 Trans. Am. Math. Soc. 50 71-128) not to have a Lagrangian. Following Bateman's dictum (1931 Phys. Rev. 38 815-9), we determine different sets of equations that are compatible with those of Douglas and derivable from a variational principle.

Predictability of the Lagrangian Motion in the Upper Ocean

NASA Astrophysics Data System (ADS)

Piterbarg, L. I.; Griffa, A.; Griffa, A.; Mariano, A. J.; Ozgokmen, T. M.; Ryan, E. H.

2001-12-01

better than the Kalman-based algorithm in strong shear flows. An important component of our research is the optimal predictor location problem; Where should floats be launched in order to minimize the Lagrangian prediction error? Preliminary Lagrangian sampling results for different flow scenarios will be presented.

GPU acceleration of Eulerian-Lagrangian particle-laden turbulent flow simulations

NASA Astrophysics Data System (ADS)

Richter, David; Sweet, James; Thain, Douglas

2017-11-01

The Lagrangian point-particle approximation is a popular numerical technique for representing dispersed phases whose properties can substantially deviate from the local fluid. In many cases, particularly in the limit of one-way coupled systems, large numbers of particles are desired; this may be either because many physical particles are present (e.g. LES of an entire cloud), or because the use of many particles increases statistical convergence (e.g. high-order statistics). Solving the trajectories of very large numbers of particles can be problematic in traditional MPI implementations, however, and this study reports the benefits of using graphical processing units (GPUs) to integrate the particle equations of motion while preserving the original MPI version of the Eulerian flow solver. It is found that GPU acceleration becomes cost effective around one million particles, and performance enhancements of up to 15x can be achieved when O(108) particles are computed on the GPU rather than the CPU cluster. Optimizations and limitations will be discussed, as will prospects for expanding to two- and four-way coupled systems. ONR Grant No. N00014-16-1-2472.

Viscoelastic Vibration Dampers for Turbomachine Blades

NASA Technical Reports Server (NTRS)

Nguyen, Nhan

2003-01-01

Simple viscoelastic dampers have been invented for use on the root attachments of turbomachine blades. These dampers suppress bending- and torsion-mode blade vibrations, which are excited by unsteady aerodynamic forces during operation. In suppressing vibrations, these dampers reduce fatigue (thereby prolonging blade lifetimes) while reducing noise. These dampers can be installed in new turbomachines or in previously constructed turbomachines, without need for structural modifications. Moreover, because these dampers are not exposed to flows, they do not affect the aerodynamic performances of turbomachines. Figure 1 depicts a basic turbomachine rotor, which includes multiple blades affixed to a hub by means of dovetail root attachments. Prior to mounting of the blades, thin layers of a viscoelastic material are applied to selected areas of the blade roots. Once the blades have been installed in the hub and the rotor is set into rotation, centrifugal force compresses these layers between the mating load-bearing surfaces of the hub and the blade root. The layers of viscoelastic material provide load paths through which the vibration energy of the blade can be dissipated. The viscoelasticity of the material converts mechanical vibration energy into shear strain energy and then from shear strain energy to heat. Of the viscoelastic materials that have been considered thus far for this application, the one of choice is a commercial polyurethane that is available in tape form, coated on one side with an adhesive that facilitates bonding to blade roots. The thickness of the tape can be chosen to suit the specific application. The typical thickness of 0.012 in. (.0.3 mm) is small enough that the tape can fit in the clearance between the mating blade-root and hub surfaces in a typical turbomachine. In an experiment, a blade was mounted in a test fixture designed to simulate the blade-end conditions that prevail in a turbocompressor. Vibrations were excited in the blade by

Dynamic viscoelasticity measurement under alternative torque using electromagnetically spinning method with quadruple electromagnets.

PubMed

Matsuura, Yusuke; Hirano, Taichi; Sakai, Keiji

2017-07-01

In this study, we developed a novel type of rheological measurement system. Here, a spherical probe is driven to rotate periodically by applying torques using quadruple electromagnets in a noncontact manner. Moreover, this system is an enhancement of our electromagnetically spinning (EMS) viscometer, which is widely used for measuring rheological flow curves in various industrial fields. The quadruple EMS method provides the frequency spectrum of viscoelasticity, in addition to shear viscosity, in a steady flow by switching the operation modes of the driving torque. We show the results obtained for Newtonian fluids and viscoelastic materials and demonstrate the validity of the system.

Viscoelastic behavior of lung and chest wall in dogs determined by flow interruption.

PubMed

Similowski, T; Levy, P; Corbeil, C; Albala, M; Pariente, R; Derenne, J P; Bates, J H; Jonson, B; Milic-Emili, J

1989-12-01

Pulmonary and chest wall mechanics were studied in six anesthetized paralyzed dogs, by use of the technique of rapid airway occlusion during constant flow inflation. Analysis of the pressure changes after flow interruption allowed us to partition the overall resistance of the lung (Rl) and chest wall (Rw) and total respiratory system (Rrs) into two components, one (Rinit) reflecting in the lung airway resistance (Raw), the other (delta R) reflecting primarily the viscoelastic properties of the pulmonary and chest wall tissues. The effects of varying inspiratory flow and inflation volume were interpreted in terms of frequency dependence of resistance, by using a spring-and-dashpot model previously proposed and substantiated by Bates et al. (Proc. 9th Annu. Conf. IEEE Med. Biol. Soc., 1987, vol. 3, p. 1802-1803). We observed that 1) Raw and Rw,init were nearly equal and small relative to Rl and Rw (both were unaffected by flow); 2) Rrs,init decreased slightly with increasing volume; 3) both delta Rl and delta Rw decreased with increasing flow and increased with increasing lung volume. These changes were manifestations of frequency dependence of delta R, as it is predicted by the model; 4) Rrs, Rl, and Rw followed the same trends as delta R. These results corroborate data previously reported in the literature with the use of different techniques to measure airways and pulmonary tissue resistances and confirm that the use of Rl to assess bronchial reactivity is problematic. The interrupter techniques provides a convenient way to obtain Raw values, as well as analogs of lung and chest wall tissue resistances in intact dogs.

Mechanochemical pattern formation in simple models of active viscoelastic fluids and solids

NASA Astrophysics Data System (ADS)

Alonso, Sergio; Radszuweit, Markus; Engel, Harald; Bär, Markus

2017-11-01

The cytoskeleton of the organism Physarum polycephalum is a prominent example of a complex active viscoelastic material wherein stresses induce flows along the organism as a result of the action of molecular motors and their regulation by calcium ions. Experiments in Physarum polycephalum have revealed a rich variety of mechanochemical patterns including standing, traveling and rotating waves that arise from instabilities of spatially homogeneous states without gradients in stresses and resulting flows. Herein, we investigate simple models where an active stress induced by molecular motors is coupled to a model describing the passive viscoelastic properties of the cellular material. Specifically, two models for viscoelastic fluids (Maxwell and Jeffrey model) and two models for viscoelastic solids (Kelvin-Voigt and Standard model) are investigated. Our focus is on the analysis of the conditions that cause destabilization of spatially homogeneous states and the related onset of mechano-chemical waves and patterns. We carry out linear stability analyses and numerical simulations in one spatial dimension for different models. In general, sufficiently strong activity leads to waves and patterns. The primary instability is stationary for all active fluids considered, whereas all active solids have an oscillatory primary instability. All instabilities found are of long-wavelength nature reflecting the conservation of the total calcium concentration in the models studied.

Unsteady force estimation using a Lagrangian drift-volume approach

NASA Astrophysics Data System (ADS)

McPhaden, Cameron J.; Rival, David E.

2018-04-01

A novel Lagrangian force estimation technique for unsteady fluid flows has been developed, using the concept of a Darwinian drift volume to measure unsteady forces on accelerating bodies. The construct of added mass in viscous flows, calculated from a series of drift volumes, is used to calculate the reaction force on an accelerating circular flat plate, containing highly-separated, vortical flow. The net displacement of fluid contained within the drift volumes is, through Darwin's drift-volume added-mass proposition, equal to the added mass of the plate and provides the reaction force of the fluid on the body. The resultant unsteady force estimates from the proposed technique are shown to align with the measured drag force associated with a rapid acceleration. The critical aspects of understanding unsteady flows, relating to peak and time-resolved forces, often lie within the acceleration phase of the motions, which are well-captured by the drift-volume approach. Therefore, this Lagrangian added-mass estimation technique opens the door to fluid-dynamic analyses in areas that, until now, were inaccessible by conventional means.

Stanley Corrsin Award Lecture: Lagrangian Measurements in Turbulence: From Fundamentals to Applications

NASA Astrophysics Data System (ADS)

Bodenschatz, Eberhard

2014-11-01

In my talk I shall present results from particle tracking experiments in turbulence. After a short review of the history of the field, I shall summarize the most recent technological advances that range form low and high-density particle tracking to direct measurements of the Lagrangian evolution of vorticity. I shall embark on a journey that describes the discoveries made possible by this new technology in the last 15 years. I present results that challenge our understanding of turbulence and show how Lagrangian particle tracking can help us ask questions on turbulent flows that so far were hidden. I shall show how Lagrangian particle tracking may provide important insights into the reversibility of turbulent flows, on vorticity generation, the energy cascade and turbulent mixing. I shall describe the consequences of inertial particle transport on rain formation and end with an outlook on how Lagrangian particle tracking experiments on non-stationary flows in real-world situations may provide high quality data that can support real world engineering problems. I am very thankful for the support by Cornell University, the National Science Foundation, the Research Corporation, the Alfred P. Sloan Foundation, the Kavli Institute for Theoretical Physics, the German Research Foundation, the European Union and the Max Planck Society. I very gratefully acknowledge the excellent partnership with many colleagues in the field of fluid mechanics and turbulence.

Lagrangian particle method for compressible fluid dynamics

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

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

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

Lagrangian particle method for compressible fluid dynamics

DOE PAGES

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

2018-02-09

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

Designing with non-linear viscoelastic fluids

NASA Astrophysics Data System (ADS)

Schuh, Jonathon; Lee, Yong Hoon; Allison, James; Ewoldt, Randy

2017-11-01

Material design is typically limited to hard materials or simple fluids; however, design with more complex materials can provide ways to enhance performance. Using the Criminale-Ericksen-Filbey (CEF) constitutive model in the thin film lubrication limit, we derive a modified Reynolds Equation (based on asymptotic analysis) that includes shear thinning, first normal stress, and terminal regime viscoelastic effects. This allows for designing non-linear viscoelastic fluids in thin-film creeping flow scenarios, i.e. optimizing the shape of rheological material properties to achieve different design objectives. We solve the modified Reynolds equation using the pseudo-spectral method, and describe a case study in full-film lubricated sliding where optimal fluid properties are identified. These material-agnostic property targets can then guide formulation of complex fluids which may use polymeric, colloidal, or other creative approaches to achieve the desired non-Newtonian properties.

Viscoelastic and elastomeric active matter: linear instability and nonlinear dynamics

NASA Astrophysics Data System (ADS)

Hemingway, Ewan J.; Cates, M. E.; Marchetti, M. C.; Fielding, S. M.

We consider a continuum model of active viscoelastic matter, whereby a model of an active nematic liquid-crystal is coupled to a minimal model of polymer dynamics with a viscoelastic relaxation time τc. To explore the resulting interplay between active and polymeric dynamics, we first generalise a linear stability analysis (from earlier studies without polymer) to derive criteria for the onset of spontaneous flow. Perhaps surprisingly, our results show that the spontaneous flow instability persists even for divergent polymer relaxation times. We explore the novel dynamical states to which these instabilities lead by means of nonlinear numerical simulations. This reveals oscillatory shear-banded states in 1D, and activity-driven turbulence in 2D, even in the limit τc --> ∞ . Adding polymer can also have calming effects, increasing the net throughput of spontaneous flow along a channel in a new type of ''drag-reduction'', an effect that may have implications for cytoplasmic streaming processes within the cell.

Second order upwind Lagrangian particle method for Euler equations

DOE PAGES

Samulyak, Roman; Chen, Hsin -Chiang; Yu, Kwangmin

2016-06-01

A new second order upwind Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multiphase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) an upwind second-order particle-based algorithm with limiter, providing accuracy and longmore » term stability, and (c) accurate resolution of states at free interfaces. In conclusion, numerical verification tests demonstrating the convergence order for fixed domain and free surface problems are presented.« less

Second order upwind Lagrangian particle method for Euler equations

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

Samulyak, Roman; Chen, Hsin -Chiang; Yu, Kwangmin

A new second order upwind Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multiphase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) an upwind second-order particle-based algorithm with limiter, providing accuracy and longmore » term stability, and (c) accurate resolution of states at free interfaces. In conclusion, numerical verification tests demonstrating the convergence order for fixed domain and free surface problems are presented.« less

Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M = 9.2 Sumatra earthquake

USGS Publications Warehouse

Pollitz, F.; Banerjee, P.; Grijalva, K.; Nagarajan, B.; Burgmann, R.

2008-01-01

The 2004 M=9.2 Sumatra-Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ???1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress-strain relation; an aspherical perturbation in viscoelastic structure; a 'static'mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra-Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation. ?? Journal compilation ?? 2008 RAS.

Dynamics and Control of Newtonian and Viscoelastic Fluids

NASA Astrophysics Data System (ADS)

Lieu, Binh K.

Transition to turbulence represents one of the most intriguing natural phenomena. Flows that are smooth and ordered may become complex and disordered as the flow strength increases. This process is known as transition to turbulence. In this dissertation, we develop theoretical and computational tools for analysis and control of transition and turbulence in shear flows of Newtonian, such as air and water, and complex viscoelastic fluids, such as polymers and molten plastics. Part I of the dissertation is devoted to the design and verification of sensor-free and feedback-based strategies for controlling the onset of turbulence in channel flows of Newtonian fluids. We use high fidelity simulations of the nonlinear flow dynamics to demonstrate the effectiveness of our model-based approach to flow control design. In Part II, we utilize systems theoretic tools to study transition and turbulence in channel flows of viscoelastic fluids. For flows with strong elastic forces, we demonstrate that flow fluctuations can experience significant amplification even in the absence of inertia. We use our theoretical developments to uncover the underlying physical mechanism that leads to this high amplification. For turbulent flows with polymer additives, we develop a model-based method for analyzing the influence of polymers on drag reduction. We demonstrate that our approach predicts drag reducing trends observed in full-scale numerical simulations. In Part III, we develop mathematical framework and computational tools for calculating frequency responses of spatially distributed systems. Using state-of-the-art automatic spectral collocation techniques and new integral formulation, we show that our approach yields more reliable and accurate solutions than currently available methods.

Dissipative inertial transport patterns near coherent Lagrangian eddies in the ocean.

PubMed

Beron-Vera, Francisco J; Olascoaga, María J; Haller, George; Farazmand, Mohammad; Triñanes, Joaquín; Wang, Yan

2015-08-01

Recent developments in dynamical systems theory have revealed long-lived and coherent Lagrangian (i.e., material) eddies in incompressible, satellite-derived surface ocean velocity fields. Paradoxically, observed drifting buoys and floating matter tend to create dissipative-looking patterns near oceanic eddies, which appear to be inconsistent with the conservative fluid particle patterns created by coherent Lagrangian eddies. Here, we show that inclusion of inertial effects (i.e., those produced by the buoyancy and size finiteness of an object) in a rotating two-dimensional incompressible flow context resolves this paradox. Specifically, we obtain that anticyclonic coherent Lagrangian eddies attract (repel) negatively (positively) buoyant finite-size particles, while cyclonic coherent Lagrangian eddies attract (repel) positively (negatively) buoyant finite-size particles. We show how these results explain dissipative-looking satellite-tracked surface drifter and subsurface float trajectories, as well as satellite-derived Sargassum distributions.

Using Lagrangian Coherent Structures to understand coastal water quality

NASA Astrophysics Data System (ADS)

Fiorentino, L. A.; Olascoaga, M. J.; Reniers, A.; Feng, Z.; Beron-Vera, F. J.; MacMahan, J. H.

2012-09-01

The accumulation of pollutants near the shoreline can result in low quality coastal water with negative effects on human health. To understand the role of mixing by tidal flows in coastal water quality we study the nearshore Lagrangian circulation. Specifically, we reveal Lagrangian Coherent Structures (LCSs), i.e., distinguished material curves which shape global mixing patterns and thus act as skeletons of the Lagrangian circulation. This is done using the recently developed geodesic theory of transport barriers. Particular focus is placed on Hobie Beach, a recreational subtropical marine beach located in Virginia Key, Miami, Florida. According to studies of water quality, Hobie Beach is characterized by high microbial levels. Possible sources of pollution in Hobie Beach include human bather shedding, dog fecal matter, runoff, and sand efflux at high tides. Consistent with the patterns formed by satellite-tracked drifter trajectories, the LCSs extracted from simulated currents reveal a Lagrangian circulation favoring the retention near the shoreline of pollutants released along the shoreline, which can help explain the low quality water registered at Hobie Beach.

Theoretical study of the flow in a fluid damper containing high viscosity silicone oil: Effects of shear-thinning and viscoelasticity

NASA Astrophysics Data System (ADS)

Syrakos, Alexandros; Dimakopoulos, Yannis; Tsamopoulos, John

2018-03-01

The flow inside a fluid damper where a piston reciprocates sinusoidally inside an outer casing containing high-viscosity silicone oil is simulated using a finite volume method, at various excitation frequencies. The oil is modeled by the Carreau-Yasuda (CY) and Phan-Thien and Tanner (PTT) constitutive equations. Both models account for shear-thinning, but only the PTT model accounts for elasticity. The CY and other generalised Newtonian models have been previously used in theoretical studies of fluid dampers, but the present study is the first to perform full two-dimensional (axisymmetric) simulations employing a viscoelastic constitutive equation. It is found that the CY and PTT predictions are similar when the excitation frequency is low, but at medium and higher frequencies, the CY model fails to describe important phenomena that are predicted by the PTT model and observed in experimental studies found in the literature, such as the hysteresis of the force-displacement and force-velocity loops. Elastic effects are quantified by applying a decomposition of the damper force into elastic and viscous components, inspired from large amplitude oscillatory shear theory. The CY model also overestimates the damper force relative to the PTT model because it underpredicts the flow development length inside the piston-cylinder gap. It is thus concluded that (a) fluid elasticity must be accounted for and (b) theoretical approaches that rely on the assumption of one-dimensional flow in the piston-cylinder gap are of limited accuracy, even if they account for fluid viscoelasticity. The consequences of using lower-viscosity silicone oil are also briefly examined.

Gradient-based Optimization for Poroelastic and Viscoelastic MR Elastography

PubMed Central

Tan, Likun; McGarry, Matthew D.J.; Van Houten, Elijah E.W.; Ji, Ming; Solamen, Ligin; Weaver, John B.

2017-01-01

We describe an efficient gradient computation for solving inverse problems arising in magnetic resonance elastography (MRE). The algorithm can be considered as a generalized ‘adjoint method’ based on a Lagrangian formulation. One requirement for the classic adjoint method is assurance of the self-adjoint property of the stiffness matrix in the elasticity problem. In this paper, we show this property is no longer a necessary condition in our algorithm, but the computational performance can be as efficient as the classic method, which involves only two forward solutions and is independent of the number of parameters to be estimated. The algorithm is developed and implemented in material property reconstructions using poroelastic and viscoelastic modeling. Various gradient- and Hessian-based optimization techniques have been tested on simulation, phantom and in vivo brain data. The numerical results show the feasibility and the efficiency of the proposed scheme for gradient calculation. PMID:27608454

Lagrangian Particle Tracking in a Discontinuous Galerkin Method for Hypersonic Reentry Flows in Dusty Environments

NASA Astrophysics Data System (ADS)

Ching, Eric; Lv, Yu; Ihme, Matthias

2017-11-01

Recent interest in human-scale missions to Mars has sparked active research into high-fidelity simulations of reentry flows. A key feature of the Mars atmosphere is the high levels of suspended dust particles, which can not only enhance erosion of thermal protection systems but also transfer energy and momentum to the shock layer, increasing surface heat fluxes. Second-order finite-volume schemes are typically employed for hypersonic flow simulations, but such schemes suffer from a number of limitations. An attractive alternative is discontinuous Galerkin methods, which benefit from arbitrarily high spatial order of accuracy, geometric flexibility, and other advantages. As such, a Lagrangian particle method is developed in a discontinuous Galerkin framework to enable the computation of particle-laden hypersonic flows. Two-way coupling between the carrier and disperse phases is considered, and an efficient particle search algorithm compatible with unstructured curved meshes is proposed. In addition, variable thermodynamic properties are considered to accommodate high-temperature gases. The performance of the particle method is demonstrated in several test cases, with focus on the accurate prediction of particle trajectories and heating augmentation. Financial support from a Stanford Graduate Fellowship and the NASA Early Career Faculty program are gratefully acknowledged.

A Lagrangian particle method with remeshing for tracer transport on the sphere

DOE PAGES

Bosler, Peter Andrew; Kent, James; Krasny, Robert; ...

2017-03-30

A Lagrangian particle method (called LPM) based on the flow map is presented for tracer transport on the sphere. The particles carry tracer values and are located at the centers and vertices of triangular Lagrangian panels. Remeshing is applied to control particle disorder and two schemes are compared, one using direct tracer interpolation and another using inverse flow map interpolation with sampling of the initial tracer density. Test cases include a moving-vortices flow and reversing-deformational flow with both zero and nonzero divergence, as well as smooth and discontinuous tracers. We examine the accuracy of the computed tracer density and tracermore » integral, and preservation of nonlinear correlation in a pair of tracers. Here, we compare results obtained using LPM and the Lin–Rood finite-volume scheme. An adaptive particle/panel refinement scheme is demonstrated.« less

A Lagrangian particle method with remeshing for tracer transport on the sphere

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

Bosler, Peter Andrew; Kent, James; Krasny, Robert

A Lagrangian particle method (called LPM) based on the flow map is presented for tracer transport on the sphere. The particles carry tracer values and are located at the centers and vertices of triangular Lagrangian panels. Remeshing is applied to control particle disorder and two schemes are compared, one using direct tracer interpolation and another using inverse flow map interpolation with sampling of the initial tracer density. Test cases include a moving-vortices flow and reversing-deformational flow with both zero and nonzero divergence, as well as smooth and discontinuous tracers. We examine the accuracy of the computed tracer density and tracermore » integral, and preservation of nonlinear correlation in a pair of tracers. Here, we compare results obtained using LPM and the Lin–Rood finite-volume scheme. An adaptive particle/panel refinement scheme is demonstrated.« less

Viscoelasticity of blood and viscoelastic blood analogues for use in polydymethylsiloxane in vitro models of the circulatory system.

PubMed

Campo-Deaño, Laura; Dullens, Roel P A; Aarts, Dirk G A L; Pinho, Fernando T; Oliveira, Mónica S N

2013-01-01

The non-Newtonian properties of blood are of great importance since they are closely related with incident cardiovascular diseases. A good understanding of the hemodynamics through the main vessels of the human circulatory system is thus fundamental in the detection and especially in the treatment of these diseases. Very often such studies take place in vitro for convenience and better flow control and these generally require blood analogue solutions that not only adequately mimic the viscoelastic properties of blood but also minimize undesirable optical distortions arising from vessel curvature that could interfere in flow visualizations or particle image velocimetry measurements. In this work, we present the viscoelastic moduli of whole human blood obtained by means of passive microrheology experiments. These results and existing shear and extensional rheological data for whole human blood in the literature enabled us to develop solutions with rheological behavior analogous to real whole blood and with a refractive index suited for PDMS (polydymethylsiloxane) micro- and milli-channels. In addition, these blood analogues can be modified in order to obtain a larger range of refractive indices from 1.38 to 1.43 to match the refractive index of several materials other than PDMS.

Viscoelasticity of blood and viscoelastic blood analogues for use in polydymethylsiloxane in vitro models of the circulatory system

PubMed Central

Campo-Deaño, Laura; Dullens, Roel P. A.; Aarts, Dirk G. A. L.; Pinho, Fernando T.; Oliveira, Mónica S. N.

2013-01-01

The non-Newtonian properties of blood are of great importance since they are closely related with incident cardiovascular diseases. A good understanding of the hemodynamics through the main vessels of the human circulatory system is thus fundamental in the detection and especially in the treatment of these diseases. Very often such studies take place in vitro for convenience and better flow control and these generally require blood analogue solutions that not only adequately mimic the viscoelastic properties of blood but also minimize undesirable optical distortions arising from vessel curvature that could interfere in flow visualizations or particle image velocimetry measurements. In this work, we present the viscoelastic moduli of whole human blood obtained by means of passive microrheology experiments. These results and existing shear and extensional rheological data for whole human blood in the literature enabled us to develop solutions with rheological behavior analogous to real whole blood and with a refractive index suited for PDMS (polydymethylsiloxane) micro- and milli-channels. In addition, these blood analogues can be modified in order to obtain a larger range of refractive indices from 1.38 to 1.43 to match the refractive index of several materials other than PDMS. PMID:24404022

Generalized Lagrangian coherent structures

NASA Astrophysics Data System (ADS)

Balasuriya, Sanjeeva; Ouellette, Nicholas T.; Rypina, Irina I.

2018-06-01

The notion of a Lagrangian Coherent Structure (LCS) is by now well established as a way to capture transient coherent transport dynamics in unsteady and aperiodic fluid flows that are known over finite time. We show that the concept of an LCS can be generalized to capture coherence in other quantities of interest that are transported by, but not fully locked to, the fluid. Such quantities include those with dynamic, biological, chemical, or thermodynamic relevance, such as temperature, pollutant concentration, vorticity, kinetic energy, plankton density, and so on. We provide a conceptual framework for identifying the Generalized Lagrangian Coherent Structures (GLCSs) associated with such evolving quantities. We show how LCSs can be seen as a special case within this framework, and provide an overarching discussion of various methods for identifying LCSs. The utility of this more general viewpoint is highlighted through a variety of examples. We also show that although LCSs approximate GLCSs in certain limiting situations under restrictive assumptions on how the velocity field affects the additional quantities of interest, LCSs are not in general sufficient to describe their coherent transport.

Influence of gas injection on viscous and viscoelastic properties of Xanthan gum.

PubMed

Bobade, Veena; Cheetham, Madalyn; Hashim, Jamal; Eshtiaghi, Nicky

2018-05-01

Xanthan gum is widely used as a model fluid for sludge to mimic the rheological behaviour under various conditions including impact of gas injection in sludge. However, there is no study to show the influence of gas injection on rheological properties of xanthan gum specifically at the concentrations at which it is used as a model fluid for sludge with solids concentration above 2%. In this paper, the rheological properties of aqueous xanthan gum solutions at different concentrations were measured over a range of gas injection flow rates. The effect of gas injection on both the flow and viscoelastic behaviour of Xanthan gum (using two different methods - a creep test and a time sweep test) was evaluated. The viscosity curve of different solid concentrations of digested sludge and waste activated sludge were compared with different solid concentrations of Xanthan gum and the results showed that Xanthan gum can mimic the flow behaviour of sludge in flow regime. The results in linear viscoelastic regime showed that increasing gas flow rate increases storage modulus (G'), indicating an increase in the intermolecular associations within the material structure leading to an increase in material strength and solid behaviour. Similarly, in creep test an increase in the gas flow rate decreased strain%, signifying that the material has become more resistant to flow. Both observed behaviour is opposite to what occurs in sludge under similar conditions. The results of both the creep test and the time sweep test indicated that choosing Xanthan gum aqueous solution as a transparent model fluid for sludge in viscoelastic regime under similar conditions involving gas injection in a concentration range studied is not feasible. However Xanthan gum can be used as a model material for sludge in flow regime; because it shows a similar behaviour to sludge. Copyright © 2018 Elsevier Ltd. All rights reserved.

Numerical simulation of viscoelastic layer rearrangement in polymer melts using OpenFOAM®

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

Köpplmayr, Thomas, E-mail: tkoepplmayr@gmail.com; Mayrhofer, Elias

In addition to their shear-thinning behavior, polymer melts are characterized by first and second normal stress differences, which cause secondary motions. Polymer coextrusion processes involve viscoelastic two-phase flows that influence layer formation. Using polymer melts with different pigmentation makes visible the layers deformed by second normal stress differences. We used a new solver for the OpenFOAM CFD toolbox which handles viscoelastic two-phase flows. A derivative of the volume-of-fluid (VoF) methodology was employed to describe the interface. Different types of polymer melt, such as polyethylene (PE), polypropylene (PP) and polyethylene terephthalate (PET) were investigated. In a coextrusion process, the less viscousmore » phase usually tends to encapsulate the more viscous one. However, the different viscoelastic properties of the melts also influence interface deformation. The materials were characterized by small-amplitude oscillatory-shear rheometry, and a multimode Giesekus model was used to fit shear viscosity, storage and loss modulus. Our simulations also took interfacial tension into account. Experimental observations and corresponding numerical simulations were found to be in good accordance.« less

A combined Eulerian-Lagrangian two-phase flow analysis of SSME HPOTP nozzle plug trajectories. II - Results

NASA Technical Reports Server (NTRS)

Mcconnaughey, P. K.; Garcia, R.; Dejong, F. J.; Sabnis, J. S.; Pribik, D. A.

1989-01-01

An analysis of Space Shuttle Main Engine high-pressure oxygen turbopump nozzle plug trajectories has been performed, using a Lagrangian method to track nozzle plug particles expelled from a turbine through a high Reynolds number flow in a turnaround duct with turning vanes. Axisymmetric and parametric analyses reveal that if nozzle plugs exited the turbine they would probably impact the LOX heat exchanger with impact velocities which are significantly less than the penetration velocity. The finding that only slight to moderate damage will result from nozzle plug failure in flight is supported by the results of a hot-fire engine test with induced nozzle plug failures.

Relativistic thermodynamics, a Lagrangian field theory for general flows including rotation

NASA Astrophysics Data System (ADS)

Frønsdal, Christian

Any theory that is based on an action principle has a much greater predictive power than one that does not have such a formulation. The formulation of a dynamical theory of General Relativity, including matter, is here viewed as a problem of coupling Einstein’s theory of pure gravity to an independently chosen and well-defined field theory of matter. It is well known that this is accomplished in a most natural way when both theories are formulated as relativistic, Lagrangian field theories, as is the case with Einstein-Maxwell theory. Special matter models of this type have been available; here a more general thermodynamical model that allows for vortex flows is presented. In a wider context, the problem of subjecting hydrodynamics and thermodynamics to an action principle is one that has been pursued for at least 150 years. A solution to this problem has been known for some time, but only under the strong restriction to potential flows. A variational principle for general flows has become available. It represents a development of the Navier-Stokes-Fourier approach to fluid dynamics. The principal innovation is the recognition that two kinds of flow velocity fields are needed, one the gradient of a scalar field and the other the time derivative of a vector field, the latter closely associated with vorticity. In the relativistic theory that is presented here, the latter is the Hodge dual of an exact 3-form, well known as the notoph field of Ogievetskij and Palubarinov, the B-field of Kalb and Ramond and the vorticity field of Lund and Regge. The total number of degrees of freedom of a unary system, including the density and the two velocity fields is 4, as expected — as in classical hydrodynamics. In this paper, we do not reduce Einstein’s dynamical equation for the metric to phenomenology, which would have denied the relevance of any intrinsic dynamics for the matter sector, nor do we abandon the equation of continuity - the very soul of hydrodynamics.

Comparison of formation of visco-elastic masses and their properties between zeins and kafirins.

PubMed

Taylor, Janet; Anyango, Joseph O; Muhiwa, Peter J; Oguntoyinbo, Segun I; Taylor, John R N

2018-04-15

Zeins of differing sub-class composition much more readily formed visco-elastic masses in water or acetic acid solutions than equivalent kafirin preparations. Visco-elastic masses could be formed from both zein and kafirin preparations by coacervation from glacial acetic acid. Dissolving the prolamins in glacial acetic acid apparently enabled protonation and complete solvation. Stress-relaxation analysis of coacervated zein and kafirin visco-elastic masses showed they were initially soft. With storage, they became much firmer. Zein masses exhibited predominantly viscous flow properties, whereas kafirin masses were more elastic. The γ-sub-class is apparently necessary for the retention of visco-elastic mass softness with kafirin and zein, and for elastic recovery of kafirin. Generally, regardless of water or acetic acid treatment, all the zein preparations had similar FTIR spectra, with greater α-helical conformation, than the kafirin preparations which were also similar to each other. Kafirin visco-elastic masses have a much higher elastic character than zein masses. Copyright © 2017 Elsevier Ltd. All rights reserved.

Experimental design for drifting buoy Lagrangian test

NASA Technical Reports Server (NTRS)

Saunders, P. M.

1975-01-01

A test of instrumentation fabricated to measure the performance of a free drifting buoy as a (Lagrangian) current meter is described. Specifically it is proposed to distinguish between the trajectory of a drogued buoy and the trajectory of the water at the level of the drogue by measuring the flow relative to the drogue.

Viscoelasticity measurement of gel formed at the liquid-liquid reactive interfaces

NASA Astrophysics Data System (ADS)

Ujiie, Tomohiro

2012-11-01

We have experimentally studied a reacting liquid flow with gel formation by using viscous fingering (VF) as a flow field. Here, two systems were employed. In one system, sodium polyacrylate (SPA) solution and ferric ion solution were used as the more and less viscous liquids, respectively. In another system, xthantan gum (XG) solution and the ferric ion solution were used as the more and less viscous liquids, respectively. We showed that influence of gel formation on VF were qualitatively different in these two systems. We consider that the difference in the two systems will be caused by the difference in the properties of the gels. Therefore, we have measured the rheological properties of the gels by means of a rheometer. In the present study, viscoelasticity measurement was performed by two methods. One is the method which uses Double Wall Ring sensor (TA instrument) and another is the method using parallel plate. In both viscoelasticity measurements, the behavior of the formed gel was qualitatively consistent. We have found that the gel in the SPA system shows viscoelastic fluid like behavior. Moreover, we have found that the gel in the XG system shows solid like behavior.

Viscoelastic fluid-structure interactions between a flexible cylinder and wormlike micelle solution

NASA Astrophysics Data System (ADS)

Dey, Anita A.; Modarres-Sadeghi, Yahya; Rothstein, Jonathan P.

2018-06-01

It is well known that when a flexible or flexibly mounted structure is placed perpendicular to the flow of a Newtonian fluid, it can oscillate due to the shedding of separated vortices at high Reynolds numbers. Unlike Newtonian fluids, the flow of viscoelastic fluids can become unstable even at infinitesimal Reynolds numbers due to a purely elastic flow instability that can occur at large Weissenberg numbers. Recent work has shown that these elastic flow instabilities can drive the motion of flexible sheets. The fluctuating fluid forces exerted on the structure from the elastic flow instabilities can lead to a coupling between an oscillatory structural motion and the state of stress in the fluid flow. In this paper, we present the results of an investigation into the flow of a viscoelastic wormlike micelle solution past a flexible circular cylinder. The time variation of the flow field and the state of stress in the fluid are shown using a combination of particle image tracking and flow-induced birefringence images. The static and dynamic responses of the flexible cylinder are presented for a range of flow velocities. The nonlinear dynamics of the structural motion is studied to better understand an observed transition from a symmetric to an asymmetric structural deformation and oscillation behavior.

Effects of Helicity on Lagrangian and Eulerian Time Correlations in Turbulence

NASA Technical Reports Server (NTRS)

Rubinstein, Robert; Zhou, Ye

1998-01-01

Taylor series expansions of turbulent time correlation functions are applied to show that helicity influences Eulerian time correlations more strongly than Lagrangian time correlations: to second order in time, the helicity effect on Lagrangian time correlations vanishes, but the helicity effect on Eulerian time correlations is nonzero. Fourier analysis shows that the helicity effect on Eulerian time correlations is confined to the largest inertial range scales. Some implications for sound radiation by swirling flows are discussed.

Learn the Lagrangian: A Vector-Valued RKHS Approach to Identifying Lagrangian Systems.

PubMed

Cheng, Ching-An; Huang, Han-Pang

2016-12-01

We study the modeling of Lagrangian systems with multiple degrees of freedom. Based on system dynamics, canonical parametric models require ad hoc derivations and sometimes simplification for a computable solution; on the other hand, due to the lack of prior knowledge in the system's structure, modern nonparametric models in machine learning face the curse of dimensionality, especially in learning large systems. In this paper, we bridge this gap by unifying the theories of Lagrangian systems and vector-valued reproducing kernel Hilbert space. We reformulate Lagrangian systems with kernels that embed the governing Euler-Lagrange equation-the Lagrangian kernels-and show that these kernels span a subspace capturing the Lagrangian's projection as inverse dynamics. By such property, our model uses only inputs and outputs as in machine learning and inherits the structured form as in system dynamics, thereby removing the need for the mundane derivations for new systems as well as the generalization problem in learning from scratches. In effect, it learns the system's Lagrangian, a simpler task than directly learning the dynamics. To demonstrate, we applied the proposed kernel to identify the robot inverse dynamics in simulations and experiments. Our results present a competitive novel approach to identifying Lagrangian systems, despite using only inputs and outputs.

Differential geometry based solvation model II: Lagrangian formulation.

PubMed

Chen, Zhan; Baker, Nathan A; Wei, G W

2011-12-01

Solvation is an elementary process in nature and is of paramount importance to more sophisticated chemical, biological and biomolecular processes. The understanding of solvation is an essential prerequisite for the quantitative description and analysis of biomolecular systems. This work presents a Lagrangian formulation of our differential geometry based solvation models. The Lagrangian representation of biomolecular surfaces has a few utilities/advantages. First, it provides an essential basis for biomolecular visualization, surface electrostatic potential map and visual perception of biomolecules. Additionally, it is consistent with the conventional setting of implicit solvent theories and thus, many existing theoretical algorithms and computational software packages can be directly employed. Finally, the Lagrangian representation does not need to resort to artificially enlarged van der Waals radii as often required by the Eulerian representation in solvation analysis. The main goal of the present work is to analyze the connection, similarity and difference between the Eulerian and Lagrangian formalisms of the solvation model. Such analysis is important to the understanding of the differential geometry based solvation model. The present model extends the scaled particle theory of nonpolar solvation model with a solvent-solute interaction potential. The nonpolar solvation model is completed with a Poisson-Boltzmann (PB) theory based polar solvation model. The differential geometry theory of surfaces is employed to provide a natural description of solvent-solute interfaces. The optimization of the total free energy functional, which encompasses the polar and nonpolar contributions, leads to coupled potential driven geometric flow and PB equations. Due to the development of singularities and nonsmooth manifolds in the Lagrangian representation, the resulting potential-driven geometric flow equation is embedded into the Eulerian representation for the purpose of

Differential geometry based solvation model II: Lagrangian formulation

PubMed Central

Chen, Zhan; Baker, Nathan A.; Wei, G. W.

2010-01-01

Solvation is an elementary process in nature and is of paramount importance to more sophisticated chemical, biological and biomolecular processes. The understanding of solvation is an essential prerequisite for the quantitative description and analysis of biomolecular systems. This work presents a Lagrangian formulation of our differential geometry based solvation model. The Lagrangian representation of biomolecular surfaces has a few utilities/advantages. First, it provides an essential basis for biomolecular visualization, surface electrostatic potential map and visual perception of biomolecules. Additionally, it is consistent with the conventional setting of implicit solvent theories and thus, many existing theoretical algorithms and computational software packages can be directly employed. Finally, the Lagrangian representation does not need to resort to artificially enlarged van der Waals radii as often required by the Eulerian representation in solvation analysis. The main goal of the present work is to analyze the connection, similarity and difference between the Eulerian and Lagrangian formalisms of the solvation model. Such analysis is important to the understanding of the differential geometry based solvation model. The present model extends the scaled particle theory (SPT) of nonpolar solvation model with a solvent-solute interaction potential. The nonpolar solvation model is completed with a Poisson-Boltzmann (PB) theory based polar solvation model. The differential geometry theory of surfaces is employed to provide a natural description of solvent-solute interfaces. The minimization of the total free energy functional, which encompasses the polar and nonpolar contributions, leads to coupled potential driven geometric flow and Poisson-Boltzmann equations. Due to the development of singularities and nonsmooth manifolds in the Lagrangian representation, the resulting potential-driven geometric flow equation is embedded into the Eulerian representation for

Geometric decomposition of the conformation tensor in viscoelastic turbulence

NASA Astrophysics Data System (ADS)

Hameduddin, Ismail; Meneveau, Charles; Zaki, Tamer A.; Gayme, Dennice F.

2018-05-01

This work introduces a mathematical approach to analysing the polymer dynamics in turbulent viscoelastic flows that uses a new geometric decomposition of the conformation tensor, along with associated scalar measures of the polymer fluctuations. The approach circumvents an inherent difficulty in traditional Reynolds decompositions of the conformation tensor: the fluctuating tensor fields are not positive-definite and so do not retain the physical meaning of the tensor. The geometric decomposition of the conformation tensor yields both mean and fluctuating tensor fields that are positive-definite. The fluctuating tensor in the present decomposition has a clear physical interpretation as a polymer deformation relative to the mean configuration. Scalar measures of this fluctuating conformation tensor are developed based on the non-Euclidean geometry of the set of positive-definite tensors. Drag-reduced viscoelastic turbulent channel flow is then used an example case study. The conformation tensor field, obtained using direct numerical simulations, is analysed using the proposed framework.

Lagrangian coherent structures along atmospheric rivers.

PubMed

Garaboa-Paz, Daniel; Eiras-Barca, Jorge; Huhn, Florian; Pérez-Muñuzuri, Vicente

2015-06-01

We show that filamentous Atmospheric Rivers (ARs) over the Northern Atlantic Ocean are closely linked to attracting Lagrangian Coherent Structures (LCSs) in the large scale wind field. The detected LCSs represent lines of attraction in the evolving flow with a significant impact on all passive tracers. Using Finite-Time Lyapunov Exponents, we extract LCSs from a two-dimensional flow derived from water vapor flux of atmospheric reanalysis data and compare them to the three-dimensional LCS obtained from the wind flow. We correlate the typical filamentous water vapor patterns of ARs with LCSs and find that LCSs bound the filaments on the back side. Passive advective transport of water vapor in the AR from tropical latitudes is potentially possible.

A new Lagrangian method for real gases at supersonic speed

NASA Technical Reports Server (NTRS)

Loh, C. Y.; Liou, Meng-Sing

1992-01-01

With the renewed interest in high speed flights, the real gas effect is of theoretical as well as practical importance. In the past decade, upwind splittings or Godunov-type Riemann solutions have received tremendous attention and as a result significant progress has been made both in the ideal and non-ideal gas. In this paper, we propose a new approach that is formulated using the Lagrangian description, for the calculation of supersonic/hypersonic real gas inviscid flows. This new formulation avoids the grid generation step which is automatically obtained as the solution procedure marches in the 'time-like' direction. As a result, no remapping is required and the accuracy is faithfully maintained in the Lagrangian level. In this paper, we give numerical results for a variety of real gas problems consisting of essential elements in high speed flows, such as shock waves, expansion waves, slip surfaces and their interactions. Finally, calculations for flows in a generic inlet and nozzle are presented.

Droplet breakup dynamics of weakly viscoelastic fluids

NASA Astrophysics Data System (ADS)

Marshall, Kristin; Walker, Travis

2016-11-01

The addition of macromolecules to solvent, even in dilute quantities, can alter a fluid's response in an extensional flow. For low-viscosity fluids, the presence of elasticity may not be apparent when measured using a standard rotational rheometer, yet it may still alter the response of a fluid when undergoing an extensional deformation, especially at small length scales where elastic effects are enhanced. Applications such as microfluidics necessitate investigating the dynamics of fluids with elastic properties that are not pronounced at large length scales. In the present work, a microfluidic cross-slot configuration is used to study the effects of elasticity on droplet breakup. Droplet breakup and the subsequent iterated-stretching - where beads form along a filament connecting two primary droplets - were observed for a variety of material and flow conditions. We present a relationship on the modes of bead formation and how and when these modes will form based on key parameters such as the properties of the outer continuous-phase fluid. The results are vital not only for simulating the droplet breakup of weakly viscoelastic fluids but also for understanding how the droplet breakup event can be used for characterizing the extensional properties of weakly-viscoelastic fluids.

Lagrangian description of warm plasmas

NASA Technical Reports Server (NTRS)

Kim, H.

1970-01-01

Efforts are described to extend the averaged Lagrangian method of describing small signal wave propagation and nonlinear wave interaction, developed by earlier workers for cold plasmas, to the more general conditions of warm collisionless plasmas, and to demonstrate particularly the effectiveness of the method in analyzing wave-wave interactions. The theory is developed for both the microscopic description and the hydrodynamic approximation to plasma behavior. First, a microscopic Lagrangian is formulated rigorously, and expanded in terms of perturbations about equilibrium. Two methods are then described for deriving a hydrodynamic Lagrangian. In the first of these, the Lagrangian is obtained by velocity integration of the exact microscopic Lagrangian. In the second, the expanded hydrodynamic Lagrangian is obtained directly from the expanded microscopic Lagrangian. As applications of the microscopic Lagrangian, the small-signal dispersion relations and the coupled mode equations are derived for all possible waves in a warm infinite, weakly inhomogeneous magnetoplasma, and their interactions are examined.

N =1 Lagrangians for generalized Argyres-Douglas theories

NASA Astrophysics Data System (ADS)

Agarwal, Prarit; Sciarappa, Antonio; Song, Jaewon

2017-10-01

We find N = 1 Lagrangian gauge theories that flow to generalized ArgyresDouglas theories with N = 2 supersymmetry. We find that certain SU quiver gauge theories flow to generalized Argyres-Douglas theories of type ( A k-1 , A mk-1) and ( I m,km , S). We also find quiver gauge theories of SO/Sp gauge groups flowing to the ( A 2 m-1 , D 2 mk+1), ( A 2 m , D 2 m( k-1)+ k ) and D m(2 k + 2) m(2 k + 2) [ m] theories.

{lambda} elements for singular problems in CFD: Viscoelastic fluids

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

Wong, K.K.; Surana, K.S.

1996-10-01

This paper presents two dimensional {lambda} element formulation for viscoelastic fluid flow containing point singularities in the flow field. The flow of viscoelastic fluid even without singularities are a difficult class of problems for increasing Deborah number or Weissenburg number due to increased dominance of convective terms and thus increased hyperbolicity. In the present work the equations of fluid motion and the constitutive laws are recast in the form of a first order system of coupled equations with the use of auxiliary variables. The velocity, pressure and stresses are interpolated using equal order C{sup 0} {lambda} element approximations. The Leastmore » Squares Finite Element Method (LSFEM) is used to construct the integral form (error functional I) corresponding to these equations. The error functional is constructed by taking the integrated sum of the squares of the errors or residuals (over the whole discretization) resulting when the element approximation is substituted into these equations. The conditions resulting from the minimization of the error functional are satisfied by using Newton`s method with line search. LSFEM has much superior performance when dealing with non-linear and convection dominated problems.« less

On steady motion of viscoelastic fluid of Oldroyd type

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

Baranovskii, E. S., E-mail: esbaranovskii@gmail.com

2014-06-01

We consider a mathematical model describing the steady motion of a viscoelastic medium of Oldroyd type under the Navier slip condition at the boundary. In the rheological relation, we use the objective regularized Jaumann derivative. We prove the solubility of the corresponding boundary-value problem in the weak setting. We obtain an estimate for the norm of a solution in terms of the data of the problem. We show that the solution set is sequentially weakly closed. Furthermore, we give an analytic solution of the boundary-value problem describing the flow of a viscoelastic fluid in a flat channel under a slipmore » condition at the walls. Bibliography: 13 titles. (paper)« less

Viscoelastic representation of surface waves in patchy saturated poroelastic media

NASA Astrophysics Data System (ADS)

Zhang, Yu; Xu, Yixian; Xia, Jianghai; Ping, Ping; Zhang, Shuangxi

2014-08-01

Wave-induced flow is observed as the dominated factor for P wave propagation at seismic frequencies. This mechanism has a mesoscopic scale nature. The inhomogeneous unsaturated patches are regarded larger than the pore size, but smaller than the wavelength. Surface wave, e.g., Rayleigh wave, which propagates along the free surface, generated by the interfering of body waves is also affected by the mesoscopic loss mechanisms. Recent studies have reported that the effect of the wave-induced flow in wave propagation shows a relaxation behavior. Viscoelastic equivalent relaxation function associated with the wave mode can describe the kinetic nature of the attenuation. In this paper, the equivalent viscoelastic relaxation functions are extended to take into account the free surface for the Rayleigh surface wave propagation in patchy saturated poroelastic media. Numerical results for the frequency-dependent velocity and attenuation and the time-dependent dynamical responses for the equivalent Rayleigh surface wave propagation along an interface between vacuum and patchy saturated porous media are reported in the low-frequency range (0.1-1,000 Hz). The results show that the dispersion and attenuation and kinetic characteristics of the mesoscopic loss effect for the surface wave can be effectively represented in the equivalent viscoelastic media. The simulation of surface wave propagation within mesoscopic patches requires solving Biot's differential equations in very small grid spaces, involving the conversion of the fast P wave energy diffusion into the Biot slow wave. This procedure requires a very large amount of computer consumption. An efficient equivalent approach for this patchy saturated poroelastic media shows a more convenient way to solve the single phase viscoelastic differential equations.

Linear viscoelasticity and thermorheological simplicity of n-hexadecane fluids under oscillatory shear via non-equilibrium molecular dynamics simulations.

PubMed

Tseng, Huan-Chang; Wu, Jiann-Shing; Chang, Rong-Yeu

2010-04-28

A small amplitude oscillatory shear flows with the classic characteristic of a phase shift when using non-equilibrium molecular dynamics simulations for n-hexadecane fluids. In a suitable range of strain amplitude, the fluid possesses significant linear viscoelastic behavior. Non-linear viscoelastic behavior of strain thinning, which means the dynamic modulus monotonously decreased with increasing strain amplitudes, was found at extreme strain amplitudes. Under isobaric conditions, different temperatures strongly affected the range of linear viscoelasticity and the slope of strain thinning. The fluid's phase states, containing solid-, liquid-, and gel-like states, can be distinguished through a criterion of the viscoelastic spectrum. As a result, a particular condition for the viscoelastic behavior of n-hexadecane molecules approaching that of the Rouse chain was obtained. Besides, more importantly, evidence of thermorheologically simple materials was presented in which the relaxation modulus obeys the time-temperature superposition principle. Therefore, using shift factors from the time-temperature superposition principle, the estimated Arrhenius flow activation energy was in good agreement with related experimental values. Furthermore, one relaxation modulus master curve well exhibited both transition and terminal zones. Especially regarding non-equilibrium thermodynamic states, variations in the density, with respect to frequencies, were revealed.

Stochastic partial differential fluid equations as a diffusive limit of deterministic Lagrangian multi-time dynamics.

PubMed

Cotter, C J; Gottwald, G A; Holm, D D

2017-09-01

In Holm (Holm 2015 Proc. R. Soc. A 471 , 20140963. (doi:10.1098/rspa.2014.0963)), stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small-scale map. We employ homogenization theory to derive effective slow stochastic particle dynamics for the resolved mean part, thereby obtaining stochastic fluid partial equations in the Eulerian formulation. To justify the application of rigorous homogenization theory, we assume mildly chaotic fast small-scale dynamics, as well as a centring condition. The latter requires that the mean of the fluctuating deviations is small, when pulled back to the mean flow.

A Thermodynamic Theory of Solid Viscoelasticity. Part II:; Nonlinear Thermo-viscoelasticity

NASA Technical Reports Server (NTRS)

Freed, Alan D.; Leonov, Arkady I.; Gray, Hugh R. (Technical Monitor)

2002-01-01

This paper, second in the series of three papers, develops a general, nonlinear, non-isothermal, compressible theory for finite rubber viscoelasticity and specifies it in a form convenient for solving problems important to the rubber, tire, automobile, and air-space industries, among others. Based on the quasi-linear approach of non-equilibrium thermodynamics, a general nonlinear theory of differential type has been developed for arbitrary non-isothermal deformations of viscoelastic solids. In this theory, the constitutive equations were presented as the sum of a rubber elastic (equilibrium) and a liquid type viscoelastic (non-equilibrium) terms. These equations have then been simplified using several modeling and simplicity arguments.

Hemodynamics in a Pediatric Ascending Aorta Using a Viscoelastic Pediatric Blood Model

PubMed Central

Good, Bryan C.; Deutsch, Steven; Manning, Keefe B.

2015-01-01

Congenital heart disease is the leading cause of infant death in the United States with over 36,000 newborns affected each year. Despite this growing problem there are few mechanical circulatory support devices designed specifically for pediatric and neonate patients. Previous research has been done investigating pediatric ventricular assist devices (PVADs) assuming blood to be a Newtonian fluid in computational fluid dynamics (CFD) simulations, ignoring its viscoelastic and shear-thinning properties. In contrast to adult VADs, PVADs may be more susceptible to hemolysis and thrombosis due to altered flow into the aorta, and therefore, a more accurate blood model should be used. A CFD solver that incorporates a modified Oldroyd-B model designed specifically for pediatric blood is used to investigate important hemodynamic parameters in a pediatric aortic model under pulsatile flow conditions. These results are compared to Newtonian blood simulations at three physiological pediatric hematocrits. Minor differences are seen in both velocity and WSS during early stages of the cardiac systole between the Newtonian and viscoelastic models. During diastole, significant differences are seen in the velocities in the descending aorta (up to 12%) and in the aortic branches (up to 30%) between the two models. Additionally, peak wall shear stress (WSS) differences are seen between the models throughout the cardiac cycle. At the onset of diastole, peak WSS differences of 43% are seen between the Newtonian and viscoelastic model and between the 20 and 60% hematocrit viscoelastic models at peak systole of 41%. PMID:26159560

An investigation on the motion and deformation of viscoelastic drops descending in another viscoelastic media

NASA Astrophysics Data System (ADS)

Davoodi, M.; Norouzi, M.

2016-10-01

In the present study, an investigation of the motion and shape deformation of drops is carried out in creeping flow to highlight the effect of viscoelastic properties on the problem. A perturbation method is employed to derive an analytical solution for the general case that both interior and exterior fluids are viscoelastic, both fluids obeying the Giesekus model. An experiment is also performed for the limiting case of an immiscible drop of a 0.03% (w/w) polyacrylamide in an 80:20 glycerol/water solution falling through a viscous Newtonian silicon oil (410 cP polydimethylsiloxane oil) in order to check the accuracy of the analytical solution. It is shown that the addition of elastic properties to the interior fluid may cause a decrease in the terminal velocity of the droplet while an increase in the elastic properties of the exterior fluid results in the opposite behavior and increases the terminal velocity. The well-known spherical shape of creeping drops for Newtonian fluids is modified by elasticity into either prolate or oblate shapes. Using the analytical solution, it is shown that normal stresses play a key role on the final steady-state shape of the drops. To keep the drops spherical in viscoelastic phases, it is shown that the effect of normal stresses on the interior and exterior media can cancel out under certain conditions. The results presented here may be of interest to industries dealing with petroleum and medicine processing, paint and power-plant related fields where knowledge of the shape and terminal velocity of descending droplets is of great importance.

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

NASA Astrophysics Data System (ADS)

Schamel, Hans

2004-03-01

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

Measurement of the viscoelastic compliance of the eustachian tube using a modified forced-response test.

PubMed

Ghadiali, Samir N; Federspiel, William J; Swarts, J Douglas; Doyle, William J

2002-01-01

Eustachian tube compliance (ETC) was suggested to be an important determinate of function. Previous attempts to quantify ETC used summary measures that are not clearly related to the physical properties of the system. Here, we present a new method for measuring ETC that conforms more closely to the engineering definition of compliance. The forced response test was modified to include oscillations in applied flow after the forced tubal opening. Pressure and flow were recorded during the standard and modified test in 12 anesthetized cynomolgus monkeys. The resulting pressure-flow, hysteresis loops were compared with those predicted by a simple fluid-structure model of the Eustachian tube with linear-elastic or viscoelastic properties. The tubal compliance index (TCI) and a viscoelastic compliance (C(v)) were calculated from these data for each monkey. The behavior of a viscoelastic, but not a linear elastic model accurately reproduced the experimental data for the monkey. The TCI and C(v) were linearly related, but the shared variance in these measures was only 63%. This new method for measuring ETC captures all information contained in the traditional TCI, but also provides information regarding the contribution of wall viscosity to Eustachian tube mechanics.

On tide-induced Lagrangian residual current and residual transport: 1. Lagrangian residual current

USGS Publications Warehouse

Feng, Shizuo; Cheng, Ralph T.; Pangen, Xi

1986-01-01

Residual currents in tidal estuaries and coastal embayments have been recognized as fundamental factors which affect the long-term transport processes. It has been pointed out by previous studies that it is more relevant to use a Lagrangian mean velocity than an Eulerian mean velocity to determine the movements of water masses. Under weakly nonlinear approximation, the parameter k, which is the ratio of the net displacement of a labeled water mass in one tidal cycle to the tidal excursion, is assumed to be small. Solutions for tides, tidal current, and residual current have been considered for two-dimensional, barotropic estuaries and coastal seas. Particular attention has been paid to the distinction between the Lagrangian and Eulerian residual currents. When k is small, the first-order Lagrangian residual is shown to be the sum of the Eulerian residual current and the Stokes drift. The Lagrangian residual drift velocity or the second-order Lagrangian residual current has been shown to be dependent on the phase of tidal current. The Lagrangian drift velocity is induced by nonlinear interactions between tides, tidal currents, and the first-order residual currents, and it takes the form of an ellipse on a hodograph plane. Several examples are given to further demonstrate the unique properties of the Lagrangian residual current.

About non standard Lagrangians in cosmology

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

Dimitrijevic, Dragoljub D.; Milosevic, Milan

A review of non standard Lagrangians present in modern cosmological models will be considered. Well known example of non standard Lagrangian is Dirac-Born-Infeld (DBI) type Lagrangian for tachyon field. Another type of non standard Lagrangian under consideration contains scalar field which describes open p-adic string tachyon and is called p-adic string theory Lagrangian. We will investigate homogenous cases of both DBI and p-adic fields and obtain Lagrangians of the standard type which have the same equations of motions as aforementioned non standard one.

Lagrangian coherent structures at the onset of hyperchaos in the two-dimensional Navier-Stokes equations.

PubMed

Miranda, Rodrigo A; Rempel, Erico L; Chian, Abraham C-L; Seehafer, Norbert; Toledo, Benjamin A; Muñoz, Pablo R

2013-09-01

We study a transition to hyperchaos in the two-dimensional incompressible Navier-Stokes equations with periodic boundary conditions and an external forcing term. Bifurcation diagrams are constructed by varying the Reynolds number, and a transition to hyperchaos (HC) is identified. Before the onset of HC, there is coexistence of two chaotic attractors and a hyperchaotic saddle. After the transition to HC, the two chaotic attractors merge with the hyperchaotic saddle, generating random switching between chaos and hyperchaos, which is responsible for intermittent bursts in the time series of energy and enstrophy. The chaotic mixing properties of the flow are characterized by detecting Lagrangian coherent structures. After the transition to HC, the flow displays complex Lagrangian patterns and an increase in the level of Lagrangian chaoticity during the bursty periods that can be predicted statistically by the hyperchaotic saddle prior to HC transition.

Stochastic partial differential fluid equations as a diffusive limit of deterministic Lagrangian multi-time dynamics

PubMed Central

Cotter, C. J.

2017-01-01

In Holm (Holm 2015 Proc. R. Soc. A 471, 20140963. (doi:10.1098/rspa.2014.0963)), stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small-scale map. We employ homogenization theory to derive effective slow stochastic particle dynamics for the resolved mean part, thereby obtaining stochastic fluid partial equations in the Eulerian formulation. To justify the application of rigorous homogenization theory, we assume mildly chaotic fast small-scale dynamics, as well as a centring condition. The latter requires that the mean of the fluctuating deviations is small, when pulled back to the mean flow. PMID:28989316

Spatial-Temporal dynamics of Newtonian and viscoelastic turbulence

NASA Astrophysics Data System (ADS)

Wang, Sung-Ning; Graham, Michael

2015-11-01

Introducing a trace amount of polymer into liquid turbulent flows can result in substantial reduction of friction drag. This phenomenon has been widely used in fluid transport, such as the Alaska crude oil pipeline. However, the mechanism is not well understood. We conduct direct numerical simulations of Newtonian and viscoelastic turbulence in large domains, in which the flow shows different characteristics in different regions. In some areas the drag is low and vortex motions are quiescent, while in other areas the drag is higher and the motions are more active. To identify these regions, we apply a statistical method, k-means clustering, which partitions the observations into k clusters by assigning each observation to its nearest centroid. The resulting partition maximizes the between-cluster variance. In the simulations, the observations are the instantaneous wall shear rate. Regions with different levels of drag are automatically identified by the partitioning algorithm. We find that the velocity profiles of the centroids exhibit characteristics similar to the individual coherent structures observed in minimal domain simulations. In addition, as viscoelasticity increases, polymer stretch becomes strongly correlated with wall shear stress. This work was supported by NSF grant CBET-1510291.

Lagrangian analysis of the laminar flat plate boundary layer

NASA Astrophysics Data System (ADS)

Gabr, Mohammad

2016-10-01

The flow properties at the leading edge of a flat plate represent a singularity to the Blasius laminar boundary layer equations; by applying the Lagrangian approach, the leading edge velocity profiles of the laminar boundary layer over a flat plate are studied. Experimental observations as well as the theoretical analysis show an exact Gaussian distribution curve as the original starting profile of the laminar flow. Comparisons between the Blasius solution and the Gaussian curve solution are carried out providing a new insight into the physics of the laminar flow.

Effect of long-time immersion of soft denture liners in water on viscoelastic properties.

PubMed

Iwasaki, Naohiko; Yamaki, Chisato; Takahashi, Hidekazu; Oki, Meiko; Suzuki, Tetsuya

2017-09-26

Aim of this study was to investigate the effect of long-time immersion of soft denture liners in 37°C water on viscoelastic properties. Six silicone-based and two acrylic resin-based soft denture liners were selected. Cylindrical specimens were stored in distilled water at 37°C for 6 months. Viscoelastic properties, which were instantaneous and delayed elastic displacements, viscous flow, and residual displacement, were determined using a creep meter, and analyzed with 2-way analysis of variance and Tukey's comparison (α=0.05). Viscoelastic properties and their time-dependent changes were varied among materials examined. The observed viscoelastic properties of three from six silicone-based liners did not significantly change after 6-month immersion, but those of two acrylic resin-based liners significantly changed with the increase of immersion time. However, the sum of initial instantaneous elastic displacement and delayed elastic displacement of two acrylic resin-based liners during 6-month immersion changed less than 10%, which might indicate clinically sufficient elastic performance.

Size-Based Separation of Particles and Cells Utilizing Viscoelastic Effects in Straight Microchannels.

PubMed

Liu, Chao; Xue, Chundong; Chen, Xiaodong; Shan, Lei; Tian, Yu; Hu, Guoqing

2015-06-16

Viscoelasticity-induced particle migration has recently received increasing attention due to its ability to obtain high-quality focusing over a wide range of flow rates. However, its application is limited to low throughput regime since the particles can defocus as flow rate increases. Using an engineered carrier medium with constant and low viscosity and strong elasticity, the sample flow rates are improved to be 1 order of magnitude higher than those in existing studies. Utilizing differential focusing of particles of different sizes, here, we present sheathless particle/cell separation in simple straight microchannels that possess excellent parallelizability for further throughput enhancement. The present method can be implemented over a wide range of particle/cell sizes and flow rates. We successfully separate small particles from larger particles, MCF-7 cells from red blood cells (RBCs), and Escherichia coli (E. coli) bacteria from RBCs in different straight microchannels. The proposed method could broaden the applications of viscoelastic microfluidic devices to particle/cell separation due to the enhanced sample throughput and simple channel design.

SPH modeling and simulation of spherical particles interacting in a viscoelastic matrix

NASA Astrophysics Data System (ADS)

Vázquez-Quesada, A.; Ellero, M.

2017-12-01

In this work, we extend the three-dimensional Smoothed Particle Hydrodynamics (SPH) non-colloidal particulate model previously developed for Newtonian suspending media in Vázquez-Quesada and Ellero ["Rheology and microstructure of non-colloidal suspensions under shear studied with smoothed particle hydrodynamics," J. Non-Newtonian Fluid Mech. 233, 37-47 (2016)] to viscoelastic matrices. For the solvent medium, the coarse-grained SPH viscoelastic formulation proposed in Vázquez-Quesada, Ellero, and Español ["Smoothed particle hydrodynamic model for viscoelastic fluids with thermal fluctuations," Phys. Rev. E 79, 056707 (2009)] is adopted. The property of this particular set of equations is that they are entirely derived within the general equation for non-equilibrium reversible-irreversible coupling formalism and therefore enjoy automatically thermodynamic consistency. The viscoelastic model is derived through a physical specification of a conformation-tensor-dependent entropy function for the fluid particles. In the simple case of suspended Hookean dumbbells, this delivers a specific SPH discretization of the Oldroyd-B constitutive equation. We validate the suspended particle model by studying the dynamics of single and mutually interacting "noncolloidal" rigid spheres under shear flow and in the presence of confinement. Numerical results agree well with available numerical and experimental data. It is straightforward to extend the particulate model to Brownian conditions and to more complex viscoelastic solvents.

Guided waves propagating in sandwich structures made of anisotropic, viscoelastic, composite materials

NASA Astrophysics Data System (ADS)

Castaings, Michel; Hosten, Bernard

2003-05-01

The propagation of Lamb-like waves in sandwich plates made of anisotropic and viscoelastic material layers is studied. A semi-analytical model is described and used for predicting the dispersion curves (phase velocity, energy velocity, and complex wave-number) and the through-thickness distribution fields (displacement, stress, and energy flow). Guided modes propagating along a test-sandwich plate are shown to be quite different than classical Lamb modes, because this structure does not have the mirror symmetry, contrary to most of composite material plates. Moreover, the viscoelastic material properties imply complex roots of the dispersion equation to be found that lead to connections between some of the dispersion curves, meaning that some of the modes get coupled together. Gradual variation from zero to nominal values of the imaginary parts of the viscoelastic moduli shows that the mode coupling depends on the level of material viscoelasticity, except for one particular case where this phenomenon exists whether the medium is viscoelastic or not. The model is used to quantify the sensitivity of both the dispersion curves and the through-thickness mode shapes to the level of material viscoelasticity, and to physically explain the mode-coupling phenomenon. Finite element software is also used to confirm results obtained for the purely elastic structure. Finally, experiments are made using ultrasonic, air-coupled transducers for generating and detecting guided modes in the test-sandwich structure. The mode-coupling phenomenon is then confirmed, and the potential of the air-coupled system for developing single-sided, contactless, NDT applications of such structures is discussed.

Cattaneo-Christov double-diffusion theory for three-dimensional flow of viscoelastic nanofluid with the effect of heat generation/absorption

NASA Astrophysics Data System (ADS)

Hayat, Tasawar; Qayyum, Sajid; Shehzad, Sabir Ali; Alsaedi, Ahmed

2018-03-01

The present research article focuses on three-dimensional flow of viscoelastic(second grade) nanofluid in the presence of Cattaneo-Christov double-diffusion theory. Flow caused is due to stretching sheet. Characteristics of heat transfer are interpreted by considering the heat generation/absorption. Nanofluid theory comprises of Brownian motion and thermophoresis. Cattaneo-Christov double-diffusion theory is introduced in the energy and concentration expressions. Such diffusions are developed as a part of formulating the thermal and solutal relaxation times framework. Suitable variables are implemented for the conversion of partial differential systems into a sets of ordinary differential equations. The transformed expressions have been explored through homotopic algorithm. Behavior of sundry variables on the velocities, temperature and concentration are scrutinized graphically. Numerical values of skin friction coefficients are also calculated and examined. Here thermal field enhances for heat generation parameter while reverse situation is noticed for heat absorption parameter.

The effect of viscoelasticity on the stability of a pulmonary airway liquid layer

NASA Astrophysics Data System (ADS)

Halpern, David; Fujioka, Hideki; Grotberg, James B.

2010-01-01

The lungs consist of a network of bifurcating airways that are lined with a thin liquid film. This film is a bilayer consisting of a mucus layer on top of a periciliary fluid layer. Mucus is a non-Newtonian fluid possessing viscoelastic characteristics. Surface tension induces flows within the layer, which may cause the lung's airways to close due to liquid plug formation if the liquid film is sufficiently thick. The stability of the liquid layer is also influenced by the viscoelastic nature of the liquid, which is modeled using the Oldroyd-B constitutive equation or as a Jeffreys fluid. To examine the role of mucus alone, a single layer of a viscoelastic fluid is considered. A system of nonlinear evolution equations is derived using lubrication theory for the film thickness and the film flow rate. A uniform film is initially perturbed and a normal mode analysis is carried out that shows that the growth rate g for a viscoelastic layer is larger than for a Newtonian fluid with the same viscosity. Closure occurs if the minimum core radius, Rmin(t), reaches zero within one breath. Solutions of the nonlinear evolution equations reveal that Rmin normally decreases to zero faster with increasing relaxation time parameter, the Weissenberg number We. For small values of the dimensionless film thickness parameter ɛ, the closure time, tc, increases slightly with We, while for moderate values of ɛ, ranging from 14% to 18% of the tube radius, tc decreases rapidly with We provided the solvent viscosity is sufficiently small. Viscoelasticity was found to have little effect for ɛ >0.18, indicating the strong influence of surface tension. The film thickness parameter ɛ and the Weissenberg number We also have a significant effect on the maximum shear stress on tube wall, max(τw), and thus, potentially, an impact on cell damage. Max(τw) increases with ɛ for fixed We, and it decreases with increasing We for small We provided the solvent viscosity parameter is sufficiently

Chemical control of the viscoelastic properties of vinylogous urethane vitrimers

PubMed Central

Denissen, Wim; Droesbeke, Martijn; Nicolaÿ, Renaud; Leibler, Ludwik; Winne, Johan M.; Du Prez, Filip E.

2017-01-01

Vinylogous urethane based vitrimers are polymer networks that have the intrinsic property to undergo network rearrangements, stress relaxation and viscoelastic flow, mediated by rapid addition/elimination reactions of free chain end amines. Here we show that the covalent exchange kinetics significantly can be influenced by combination with various simple additives. As anticipated, the exchange reactions on network level can be further accelerated using either Brønsted or Lewis acid additives. Remarkably, however, a strong inhibitory effect is observed when a base is added to the polymer matrix. These effects have been mechanistically rationalized, guided by low-molecular weight kinetic model experiments. Thus, vitrimer elastomer materials can be rationally designed to display a wide range of viscoelastic properties. PMID:28317893

Ultrasonic measurements of breast viscoelasticity.

PubMed

Sridhar, Mallika; Insana, Michael F

2007-12-01

In vivo measurements of the viscoelastic properties of breast tissue are described. Ultrasonic echo frames were recorded from volunteers at 5 fps while applying a uniaxial compressive force (1-20 N) within a 1 s ramp time and holding the force constant for up to 200 s. A time series of strain images was formed from the echo data, spatially averaged viscous creep curves were computed, and viscoelastic strain parameters were estimated by fitting creep curves to a second-order Voigt model. The useful strain bandwidth from this quasi-static ramp stimulus was 10(-2) < or = omega < or = 10(0) rad/s (0.0016-0.16 Hz). The stress-strain curves for normal glandular tissues are linear when the surface force applied is between 2 and 5 N. In this range, the creep response was characteristic of biphasic viscoelastic polymers, settling to a constant strain (arrheodictic) after 100 s. The average model-based retardance time constants for the viscoelastic response were 3.2 +/- 0.8 and 42.0 +/- 28 s. Also, the viscoelastic strain amplitude was approximately equal to that of the elastic strain. Above 5 N of applied force, however, the response of glandular tissue became increasingly nonlinear and rheodictic, i.e., tissue creep never reached a plateau. Contrasting in vivo breast measurements with those in gelatin hydrogels, preliminary ideas regarding the mechanisms for viscoelastic contrast are emerging.

Laboratory experiment on the 3D tide-induced Lagrangian residual current using the PIV technique

NASA Astrophysics Data System (ADS)

Chen, Yang; Jiang, Wensheng; Chen, Xu; Wang, Tao; Bian, Changwei

2017-12-01

The 3D structure of the tide-induced Lagrangian residual current was studied using the particle image velocimetry (PIV) technique in a long shallow narrow tank in the laboratory. At the mouth of the tank, a wave generator was used to make periodic wave which represents the tide movement, and at the head of the tank, a laterally sloping topography with the length of one fifth of the water tank was installed, above which the tide-induced Lagrangian residual current was studied. Under the weakly nonlinear condition in the present experiment setup, the results show that the Lagrangian residual velocity (LRV) field has a three-layer structure. The residual current flows inwards (towards the head) in the bottom layer and flows outwards in the middle layer, while in the surface layer, it flows inwards along the shallow side of the sloping topography and outwards along the deep side. The depth-averaged and breadth-averaged LRV are also analyzed based on the 3D LRV observations. Our results are in good agreement with the previous experiment studies, the analytical solutions with similar conditions and the observational results in real bays. Moreover, the volume flux comparison between the Lagrangian and Eulerian residual currents shows that the Eulerian residual velocity violates the mass conservation law while the LRV truly represents the inter-tidal water transport. This work enriches the laboratory studies of the LRV and offers valuable references for the LRV studies in real bays.

Evaluation of the HF-Radar network system around Taiwan using normalized cumulative Lagrangian separation.

NASA Astrophysics Data System (ADS)

Fredj, Erick; Kohut, Josh; Roarty, Hugh; Lai, Jian-Wu

2017-04-01

The Lagrangian separation distance between the endpoints of simulated and observed drifter trajectories is often used to assess the performance of numerical particle trajectory models. However, the separation distance fails to indicate relative model performance in weak and strong current regions, such as over continental shelves and the adjacent deep ocean. A skill score described in detail by (Lui et.al. 2011) was applied to estimate the cumulative Lagrangian separation distances normalized by the associated cumulative trajectory lengths. In contrast, the Lagrangian separation distance alone gives a misleading result. The proposed dimensionless skill score is particularly useful when the number of drifter trajectories is limited and neither a conventional Eulerian-based velocity nor a Lagrangian based probability density function may be estimated. The skill score assesses The Taiwan Ocean Radar Observing System (TOROS) performance. TOROS consists of 17 SeaSonde type radars around the Taiwan Island. The currents off Taiwan are significantly influenced by the nearby Kuroshio current. The main stream of the Kuroshio flows along the east coast of Taiwan to the north throughout the year. Sometimes its branch current also bypasses the south end of Taiwan and goes north along the west coast of Taiwan. The Kuroshio is also prone to seasonal change in its speed of flow, current capacity, distribution width, and depth. The evaluations of HF-Radar National Taiwanese network performance using Lagrangian drifter records demonstrated the high quality and robustness of TOROS HF-Radar data using a purely trajectory-based non-dimensional index. Yonggang Liu and Robert H. Weisberg, "Evaluation of trajectory modeling in different dynamic regions using normalized cumulative Lagrangian separation", Journal of Geophysical Research, Vol. 116, C09013, doi:10.1029/2010JC006837, 2011

Negative vortices: The formation of vortex rings with reversed rotation in viscoelastic liquids

NASA Astrophysics Data System (ADS)

Palacios-Morales, Carlos; Barbosa, Christophe; Solorio, Francisco; Zenit, Roberto

2015-05-01

The formation process of vortex rings in a viscoelastic liquid is studied experimentally considering a piston-cylinder arrangement. Initially, a vortex ring begins to form as fluid is injected from the cylinder into the tank in a manner similar to that observed for Newtonian liquids. For later times, when the piston ceases its motion, the flow changes dramatically. A secondary vortex with reversed spinning direction appears and grows to be as large in size as the original one. The formation process is studied by contrasting the evolution with that obtained for Newtonian liquids with equivalent Reynolds numbers and stroke ratios. We argue that the reversing flow, or negative vortex, results from the combined action of shear and extension rates produced during the vortex formation, in a process similar to that observed behind ascending bubbles and falling spheres in viscoelastic media.

A Thermodynamic Theory Of Solid Viscoelasticity. Part 1: Linear Viscoelasticity.

NASA Technical Reports Server (NTRS)

Freed, Alan D.; Leonov, Arkady I.

2002-01-01

The present series of three consecutive papers develops a general theory for linear and finite solid viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the general approach is specified in a form convenient for solving problems important for many industries that involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear viscoelasticity. The most important is the thermodynamic modeling of thermo-rheological complexity , i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A special structure of interaction matrices is established for different physical mechanisms contributed to the normal relaxation modes. This structure seems to be in accord with observations, and creates a simple mathematical framework for both continuum and molecular theories of the thermo-rheological complex relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers in the rubbery and glassy regions.

A semi-Lagrangian transport method for kinetic problems with application to dense-to-dilute polydisperse reacting spray flows

NASA Astrophysics Data System (ADS)

Doisneau, François; Arienti, Marco; Oefelein, Joseph C.

2017-01-01

For sprays, as described by a kinetic disperse phase model strongly coupled to the Navier-Stokes equations, the resolution strategy is constrained by accuracy objectives, robustness needs, and the computing architecture. In order to leverage the good properties of the Eulerian formalism, we introduce a deterministic particle-based numerical method to solve transport in physical space, which is simple to adapt to the many types of closures and moment systems. The method is inspired by the semi-Lagrangian schemes, developed for Gas Dynamics. We show how semi-Lagrangian formulations are relevant for a disperse phase far from equilibrium and where the particle-particle coupling barely influences the transport; i.e., when particle pressure is negligible. The particle behavior is indeed close to free streaming. The new method uses the assumption of parcel transport and avoids to compute fluxes and their limiters, which makes it robust. It is a deterministic resolution method so that it does not require efforts on statistical convergence, noise control, or post-processing. All couplings are done among data under the form of Eulerian fields, which allows one to use efficient algorithms and to anticipate the computational load. This makes the method both accurate and efficient in the context of parallel computing. After a complete verification of the new transport method on various academic test cases, we demonstrate the overall strategy's ability to solve a strongly-coupled liquid jet with fine spatial resolution and we apply it to the case of high-fidelity Large Eddy Simulation of a dense spray flow. A fuel spray is simulated after atomization at Diesel engine combustion chamber conditions. The large, parallel, strongly coupled computation proves the efficiency of the method for dense, polydisperse, reacting spray flows.

A semi-Lagrangian transport method for kinetic problems with application to dense-to-dilute polydisperse reacting spray flows

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

Doisneau, François, E-mail: fdoisne@sandia.gov; Arienti, Marco, E-mail: marient@sandia.gov; Oefelein, Joseph C., E-mail: oefelei@sandia.gov

For sprays, as described by a kinetic disperse phase model strongly coupled to the Navier–Stokes equations, the resolution strategy is constrained by accuracy objectives, robustness needs, and the computing architecture. In order to leverage the good properties of the Eulerian formalism, we introduce a deterministic particle-based numerical method to solve transport in physical space, which is simple to adapt to the many types of closures and moment systems. The method is inspired by the semi-Lagrangian schemes, developed for Gas Dynamics. We show how semi-Lagrangian formulations are relevant for a disperse phase far from equilibrium and where the particle–particle coupling barelymore » influences the transport; i.e., when particle pressure is negligible. The particle behavior is indeed close to free streaming. The new method uses the assumption of parcel transport and avoids to compute fluxes and their limiters, which makes it robust. It is a deterministic resolution method so that it does not require efforts on statistical convergence, noise control, or post-processing. All couplings are done among data under the form of Eulerian fields, which allows one to use efficient algorithms and to anticipate the computational load. This makes the method both accurate and efficient in the context of parallel computing. After a complete verification of the new transport method on various academic test cases, we demonstrate the overall strategy's ability to solve a strongly-coupled liquid jet with fine spatial resolution and we apply it to the case of high-fidelity Large Eddy Simulation of a dense spray flow. A fuel spray is simulated after atomization at Diesel engine combustion chamber conditions. The large, parallel, strongly coupled computation proves the efficiency of the method for dense, polydisperse, reacting spray flows.« less

Lagrangian formulation and symmetrical description of liquid dynamics.

PubMed

Trachenko, K

2017-12-01

Theoretical description of liquids has been primarily based on the hydrodynamic approach and its generalization to the solid-like regime. We show that the same liquid properties can be derived starting from solid-like equations and generalizing them to account for the hydrodynamic flow. Both approaches predict propagating shear waves with the notable gap in k-space. This gives an important symmetry of liquids regarding their description. We subsequently construct a two-field Lagrangian of liquid dynamics where the dissipative hydrodynamic and solid-like terms are treated on equal footing. The Lagrangian predicts two gapped waves propagating in opposite space-time directions. The dissipative and mass terms compete by promoting gaps in k-space and energy, respectively. When bare mass is close to the field hopping frequency, both gaps close and the dissipative term annihilates the bare mass.

Lagrangian formulation and symmetrical description of liquid dynamics

NASA Astrophysics Data System (ADS)

Trachenko, K.

2017-12-01

Theoretical description of liquids has been primarily based on the hydrodynamic approach and its generalization to the solid-like regime. We show that the same liquid properties can be derived starting from solid-like equations and generalizing them to account for the hydrodynamic flow. Both approaches predict propagating shear waves with the notable gap in k -space. This gives an important symmetry of liquids regarding their description. We subsequently construct a two-field Lagrangian of liquid dynamics where the dissipative hydrodynamic and solid-like terms are treated on equal footing. The Lagrangian predicts two gapped waves propagating in opposite space-time directions. The dissipative and mass terms compete by promoting gaps in k -space and energy, respectively. When bare mass is close to the field hopping frequency, both gaps close and the dissipative term annihilates the bare mass.

Viscoelastic Properties of Human Tracheal Tissues.

PubMed

Safshekan, Farzaneh; Tafazzoli-Shadpour, Mohammad; Abdouss, Majid; Shadmehr, Mohammad B

2017-01-01

The physiological performance of trachea is highly dependent on its mechanical behavior, and therefore, the mechanical properties of its components. Mechanical characterization of trachea is key to succeed in new treatments such as tissue engineering, which requires the utilization of scaffolds which are mechanically compatible with the native human trachea. In this study, after isolating human trachea samples from brain-dead cases and proper storage, we assessed the viscoelastic properties of tracheal cartilage, smooth muscle, and connective tissue based on stress relaxation tests (at 5% and 10% strains for cartilage and 20%, 30%, and 40% for smooth muscle and connective tissue). After investigation of viscoelastic linearity, constitutive models including Prony series for linear viscoelasticity and quasi-linear viscoelastic, modified superposition, and Schapery models for nonlinear viscoelasticity were fitted to the experimental data to find the best model for each tissue. We also investigated the effect of age on the viscoelastic behavior of tracheal tissues. Based on the results, all three tissues exhibited a (nonsignificant) decrease in relaxation rate with increasing the strain, indicating viscoelastic nonlinearity which was most evident for cartilage and with the least effect for connective tissue. The three-term Prony model was selected for describing the linear viscoelasticity. Among different models, the modified superposition model was best able to capture the relaxation behavior of the three tracheal components. We observed a general (but not significant) stiffening of tracheal cartilage and connective tissue with aging. No change in the stress relaxation percentage with aging was observed. The results of this study may be useful in the design and fabrication of tracheal tissue engineering scaffolds.

Ultrasonic measurements of breast viscoelasticity

PubMed Central

Sridhar, Mallika; Insana, Michael F.

2009-01-01

In vivo measurements of the viscoelastic properties of breast tissue are described. Ultrasonic echo frames were recorded from volunteers at 5 fps while applying a uniaxial compressive force (1–20 N) within a 1 s ramp time and holding the force constant for up to 200 s. A time series of strain images was formed from the echo data, spatially averaged viscous creep curves were computed, and viscoelastic strain parameters were estimated by fitting creep curves to a second-order Voigt model. The useful strain bandwidth from this quasi-static ramp stimulus was 10−2 ≤ ω ≤ 100 rad/s (0.0016–0.16 Hz). The stress-strain curves for normal glandular tissues are linear when the surface force applied is between 2 and 5 N. In this range, the creep response was characteristic of biphasic viscoelastic polymers, settling to a constant strain (arrheodictic) after 100 s. The average model-based retardance time constants for the viscoelastic response were 3.2±0.8 and 42.0±28 s. Also, the viscoelastic strain amplitude was approximately equal to that of the elastic strain. Above 5 N of applied force, however, the response of glandular tissue became increasingly nonlinear and rheodictic, i.e., tissue creep never reached a plateau. Contrasting in vivo breast measurements with those in gelatin hydrogels, preliminary ideas regarding the mechanisms for viscoelastic contrast are emerging. PMID:18196803

Form of the manifestly covariant Lagrangian

NASA Astrophysics Data System (ADS)

Johns, Oliver Davis

1985-10-01

The preferred form for the manifestly covariant Lagrangian function of a single, charged particle in a given electromagnetic field is the subject of some disagreement in the textbooks. Some authors use a ``homogeneous'' Lagrangian and others use a ``modified'' form in which the covariant Hamiltonian function is made to be nonzero. We argue in favor of the ``homogeneous'' form. We show that the covariant Lagrangian theories can be understood only if one is careful to distinguish quantities evaluated on the varied (in the sense of the calculus of variations) world lines from quantities evaluated on the unvaried world lines. By making this distinction, we are able to derive the Hamilton-Jacobi and Klein-Gordon equations from the ``homogeneous'' Lagrangian, even though the covariant Hamiltonian function is identically zero on all world lines. The derivation of the Klein-Gordon equation in particular gives Lagrangian theoretical support to the derivations found in standard quantum texts, and is also shown to be consistent with the Feynman path-integral method. We conclude that the ``homogeneous'' Lagrangian is a completely adequate basis for covariant Lagrangian theory both in classical and quantum mechanics. The article also explores the analogy with the Fermat theorem of optics, and illustrates a simple invariant notation for the Lagrangian and other four-vector equations.

LSPRAY-V: A Lagrangian Spray Module

NASA Technical Reports Server (NTRS)

Raju, M. S.

2015-01-01

LSPRAY-V is a Lagrangian spray solver developed for application with unstructured grids and massively parallel computers. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray encountered over a wide range of operating conditions in modern aircraft engine development. It could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. With the development of LSPRAY-V, we have advanced the state-of-the-art in spray computations in several important ways.

Floating shock fitting via Lagrangian adaptive meshes

NASA Technical Reports Server (NTRS)

Vanrosendale, John

1995-01-01

In recent work we have formulated a new approach to compressible flow simulation, combining the advantages of shock-fitting and shock-capturing. Using a cell-centered on Roe scheme discretization on unstructured meshes, we warp the mesh while marching to steady state, so that mesh edges align with shocks and other discontinuities. This new algorithm, the Shock-fitting Lagrangian Adaptive Method (SLAM), is, in effect, a reliable shock-capturing algorithm which yields shock-fitted accuracy at convergence.

LSPRAY-III: A Lagrangian Spray Module

NASA Technical Reports Server (NTRS)

Raju, M. S.

2008-01-01

LSPRAY-III is a Lagrangian spray solver developed for application with parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type for the gas flow grid representation. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray because of its importance in aerospace application. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. With the development of LSPRAY-III, we have advanced the state-of-the-art in spray computations in several important ways.

LSPRAY-II: A Lagrangian Spray Module

NASA Technical Reports Server (NTRS)

Raju, M. S.

2004-01-01

LSPRAY-II is a Lagrangian spray solver developed for application with parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo Probability Density Function (PDF) solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type for the gas flow grid representation. It is mainly designed to predict the flow, thermal and transport properties of a rapidly vaporizing spray because of its importance in aerospace application. The manual provides the user with an understanding of various models involved in the spray formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. With the development of LSPRAY-II, we have advanced the state-of-the-art in spray computations in several important ways.

Lagrangian analysis of premixed turbulent combustion in hydrogen-air flames

NASA Astrophysics Data System (ADS)

Darragh, Ryan; Poludnenko, Alexei; Hamlington, Peter

2016-11-01

Lagrangian analysis has long been a tool used to analyze non-reacting turbulent flows, and has recently gained attention in the reacting flow and combustion communities. The approach itself allows one to separate local molecular effects, such as those due to reactions or diffusion, from turbulent advective effects along fluid pathlines, or trajectories. Accurate calculation of these trajectories can, however, be rather difficult due to the chaotic nature of turbulent flows and the added complexity of reactions. In order to determine resolution requirements and verify the numerical algorithm, extensive tests are described in this talk for prescribed steady, unsteady, and chaotic flows, as well as for direct numerical simulations (DNS) of non-reacting homogeneous isotropic turbulence. The Lagrangian analysis is then applied to DNS of premixed hydrogen-air flames at two different turbulence intensities for both single- and multi-step chemical mechanisms. Non-monotonic temperature and fuel-mass fraction evolutions are found to exist along trajectories passing through the flame brush. Such non-monotonicity is shown to be due to molecular diffusion resulting from large spatial gradients created by turbulent advection. This work was supported by the Air Force Office of Scientific Research (AFOSR) under Award No. FA9550-14-1-0273, and the Department of Defense (DoD) High Performance Computing Modernization Program (HPCMP) under a Frontier project award.

Continuous and Pulsatile Pediatric Ventricular Assist Device Hemodynamics with a Viscoelastic Blood Model

PubMed Central

Good, Bryan C.; Deutsch, Steven; Manning, Keefe B.

2015-01-01

Purpose To investigate the effects of pulsatile and continuous pediatric ventricular assist (PVAD) flow and pediatric blood viscoelasticity on hemodynamics in a pediatric aortic graft model. Methods Hemodynamic parameters of pulsatility, along with velocity and wall shear stress (WSS), are analyzed and compared between Newtonian and viscoelastic blood models at a range of physiological pediatric hematocrits using computational fluid dynamics. Results Both pulsatile and continuous PVAD flow lead to a decrease in pulsatility (surplus hemodynamic energy (SHE), ergs/cm3) compared to healthy aortic flow but with continuous PVAD pulsatility up to 2.4 times lower than pulsatile PVAD pulsatility at each aortic outlet. Significant differences are also seen between the two flow modes in velocity and WSS. The higher velocity jet during systole with pulsatile flow leads to higher WSSs at the anastomotic toe and at the aortic branch bifurcations. The lower velocity but continuous flow jet leads to a much different flow field and higher WSSs into diastole. Under a range of physiological pediatric hematocrit (20-60%), both velocity and WSS can vary significantly with the higher hematocrit blood model generally leading to higher peak WSSs but also lower WSSs in regions of flow separation. Conclusions The large decrease in pulsatility seen from continuous PVAD flow could lead to complications in pediatric vascular development while the high WSSs during peak systole from pulsatile PVAD flow could lead to blood damage. Both flow modes lead to similar regions prone to intimal hyperplasia (IH) resulting from low time-averaged WSS (TAWSS) and high oscillatory shear index (OSI). PMID:26643646

A Combined Eulerian-Lagrangian Data Representation for Large-Scale Applications.

PubMed

Sauer, Franz; Xie, Jinrong; Ma, Kwan-Liu

2017-10-01

The Eulerian and Lagrangian reference frames each provide a unique perspective when studying and visualizing results from scientific systems. As a result, many large-scale simulations produce data in both formats, and analysis tasks that simultaneously utilize information from both representations are becoming increasingly popular. However, due to their fundamentally different nature, drawing correlations between these data formats is a computationally difficult task, especially in a large-scale setting. In this work, we present a new data representation which combines both reference frames into a joint Eulerian-Lagrangian format. By reorganizing Lagrangian information according to the Eulerian simulation grid into a "unit cell" based approach, we can provide an efficient out-of-core means of sampling, querying, and operating with both representations simultaneously. We also extend this design to generate multi-resolution subsets of the full data to suit the viewer's needs and provide a fast flow-aware trajectory construction scheme. We demonstrate the effectiveness of our method using three large-scale real world scientific datasets and provide insight into the types of performance gains that can be achieved.

A Thermodynamic Theory of Solid Viscoelasticity. Part 3: Nonlinear Glassy Viscoelasticity, Stability Constraints, Specifications

NASA Technical Reports Server (NTRS)

Freed, Alan; Leonov, Arkady I.

2002-01-01

This paper, the last in the series, continues developing the nonlinear constitutive relations for non-isothermal, compressible, solid viscoelasticity. We initially discuss a single integral approach, more suitable for the glassy state of rubber-like materials, with basic functionals involved in the thermodynamic description for this type of viscoelasticity. Then we switch our attention to analyzing stability constraints, imposed on the general formulation of the nonlinear theory of solid viscoelasticity. Finally, we discuss specific (known from the literature or new) expressions for material functions that are involved in the constitutive formulations of both the rubber-like and glassy-like, complementary parts of the theory.

Flow and fracturing of viscoelastic media under diffusion-driven bubble growth: An analogue experiment for eruptive volcanic conduits

NASA Astrophysics Data System (ADS)

Taddeucci, J.; Spieler, O.; Ichihara, M.; Dingwell, D. B.; Scarlato, P.

2006-03-01

To visualize the behavior of erupting magma in volcanic conduits, we performed shock tube experiments on the ductile-brittle response of a viscoelastic medium to diffusion-driven bubble expansion. A sample of shear-thinning magma analogue is saturated by gas Ar under high pressure. On rapid decompression, Ar supersaturation causes bubbles to nucleate, grow, and coalesce in the sample, forcing it to expand, flow, and fracture. Experimental variables include saturation pressure and duration, and shape and lubrication of the flow path. Bubble growth in the experiments controls both flow and fracturing, and is consistent with physical models of magma vesiculation. Two types of fractures are observed: i) sharp fractures along the uppermost rim of the sample, and ii) fractures pervasively diffused throughout the sample. Rim fractures open when shear stress accumulates and strain rate is highest at the margin of the flow (a process already inferred from observations and models to occur in magma). Pervasive fractures originate when wall-friction retards expansion of the sample, causing pressure to build-up in the bubbles. When bubble pressure overcomes wall-friction and the tensile strength of the porous sample, fractures open with a range of morphologies. Both types of fracture open normally to flow direction, and both may heal as the flow proceeds. These experiments also illustrate how the development of pervasive fractures allows exsolving gas to escape from the sample before the generation of a permeable network via other processes, e.g., bubble coalescence. This is an observation that potentially impact the degassing of magma and the transition between explosive and effusive eruptions.

Droplet impact on soft viscoelastic surfaces.

PubMed

Chen, Longquan; Bonaccurso, Elmar; Deng, Peigang; Zhang, Haibo

2016-12-01

In this work, we experimentally investigate the impact of water droplets onto soft viscoelastic surfaces with a wide range of impact velocities. Several impact phenomena, which depend on the dynamic interaction between the droplets and viscoelastic surfaces, have been identified and analyzed. At low We, complete rebound is observed when the impact velocity is between a lower and an upper threshold, beyond which droplets are deposited on the surface after impact. At intermediate We, entrapment of an air bubble inside the impinging droplets is found on soft surfaces, while a bubble entrapment on the surface is observed on rigid surfaces. At high We, partial rebound is only identified on the most rigid surface at We≳92. Rebounding droplets behave similarly to elastic drops rebounding on superhydrophobic surfaces and the impact process is independent of surface viscoelasticity. Further, surface viscoelasticity does not influence drop spreading after impact-as the surfaces behave like rigid surfaces-but it does affect drop recoiling. Also, the postimpact drop oscillation on soft viscoelastic surfaces is influenced by dynamic wettability of these surfaces. Comparing sessile drop oscillation with a damped harmonic oscillator allows us to conclude that surface viscoelasticity affects the damping coefficient and liquid surface tension sets the spring constant of the system.

Viscoelastic dynamic arterial response.

PubMed

Charalambous, Haralambia P; Roussis, Panayiotis C; Giannakopoulos, Antonios E

2017-10-01

Arteries undergo large deformations under applied intraluminal pressure and may exhibit small hysteresis due to creep or relaxation process. The mechanical response of arteries depends, among others, on their topology along the arterial tree. Viscoelasticity of arterial tissues, which is the topic investigated in this study, is mainly a characteristic mechanical response of arteries that are located away from the heart and have increased smooth muscle cells content. The arterial wall viscosity is simulated by adopting a generalized Maxwell model and the method of internal variables, as proposed by Bonet and Holzapfel et al. The total stresses consist of elastic long-term stresses and viscoelastic stresses, requiring an iterative procedure for their calculation. The cross-section of the artery is modeled as a circular ring, consisting of a single homogenized layer, under a time-varying blood pressure. Two different loading approximations for the aortic pressure vs time are considered. A novel numerical method is developed in order to solve the controlling integro-differential equation. A large number of numerical investigations are performed and typical response time-profiles are presented in pictorial form. Results suggest that the viscoelastic arterial response is mainly affected by the ratio of the relaxation time to the characteristic time of the response and by the pressure-time approximation. Numerical examples, based on data available in the literature, are conducted. The investigation presented in this study reveals the effect of each material parameter on the viscoelastic arterial response. Thus, a better understanding of the behavior of viscoelastic arteries is achieved. Copyright © 2017 Elsevier Ltd. All rights reserved.

Lagrangian numerical methods for ocean biogeochemical simulations

NASA Astrophysics Data System (ADS)

Paparella, Francesco; Popolizio, Marina

2018-05-01

We propose two closely-related Lagrangian numerical methods for the simulation of physical processes involving advection, reaction and diffusion. The methods are intended to be used in settings where the flow is nearly incompressible and the Péclet numbers are so high that resolving all the scales of motion is unfeasible. This is commonplace in ocean flows. Our methods consist in augmenting the method of characteristics, which is suitable for advection-reaction problems, with couplings among nearby particles, producing fluxes that mimic diffusion, or unresolved small-scale transport. The methods conserve mass, obey the maximum principle, and allow to tune the strength of the diffusive terms down to zero, while avoiding unwanted numerical dissipation effects.

Direct Lagrangian tracking simulations of particles in vertically-developing atmospheric clouds

NASA Astrophysics Data System (ADS)

Onishi, Ryo; Kunishima, Yuichi

2017-11-01

We have been developing the Lagrangian Cloud Simulator (LCS), which follows the so-called Euler-Lagrangian framework, where flow motion and scalar transportations (i.e., temperature and humidity) are computed with the Euler method and particle motion with the Lagrangian method. The LCS simulation considers the hydrodynamic interaction between approaching particles for robust collision detection. This leads to reliable simulations of collision growth of cloud droplets. Recently the activation process, in which aerosol particles become tiny liquid droplets, has been implemented in the LCS. The present LCS can therefore consider the whole warm-rain precipitation processes -activation, condensation, collision and drop precipitation. In this talk, after briefly introducing the LCS, we will show kinematic simulations using the LCS for quasi-one dimensional domain, i.e., vertically elongated 3D domain. They are compared with one-dimensional kinematic simulations using a spectral-bin cloud microphysics scheme, which is based on the Euler method. The comparisons show fairly good agreement with small discrepancies, the source of which will be presented. The Lagrangian statistics, obtained for the first time for the vertical domain, will be the center of discussion. This research was supported by MEXT as ``Exploratory Challenge on Post-K computer'' (Frontiers of Basic Science: Challenging the Limits).

Evaluation of wastewater contaminant transport in surface waters using verified Lagrangian sampling.

PubMed

Antweiler, Ronald C; Writer, Jeffrey H; Murphy, Sheila F

2014-02-01

Contaminants released from wastewater treatment plants can persist in surface waters for substantial distances. Much research has gone into evaluating the fate and transport of these contaminants, but this work has often assumed constant flow from wastewater treatment plants. However, effluent discharge commonly varies widely over a 24-hour period, and this variation controls contaminant loading and can profoundly influence interpretations of environmental data. We show that methodologies relying on the normalization of downstream data to conservative elements can give spurious results, and should not be used unless it can be verified that the same parcel of water was sampled. Lagrangian sampling, which in theory samples the same water parcel as it moves downstream (the Lagrangian parcel), links hydrologic and chemical transformation processes so that the in-stream fate of wastewater contaminants can be quantitatively evaluated. However, precise Lagrangian sampling is difficult, and small deviations - such as missing the Lagrangian parcel by less than 1h - can cause large differences in measured concentrations of all dissolved compounds at downstream sites, leading to erroneous conclusions regarding in-stream processes controlling the fate and transport of wastewater contaminants. Therefore, we have developed a method termed "verified Lagrangian" sampling, which can be used to determine if the Lagrangian parcel was actually sampled, and if it was not, a means for correcting the data to reflect the concentrations which would have been obtained had the Lagrangian parcel been sampled. To apply the method, it is necessary to have concentration data for a number of conservative constituents from the upstream, effluent, and downstream sites, along with upstream and effluent concentrations that are constant over the short-term (typically 2-4h). These corrections can subsequently be applied to all data, including non-conservative constituents. Finally, we show how data

Viscoelastic Postseismic Rebound to Strike-Slip Earthquakes in Regions of Oblique Plate Convergence

NASA Technical Reports Server (NTRS)

Cohen, Steven C.

1999-01-01

According to the slip partitioning concept, the trench parallel component of relative plate motion in regions of oblique convergence is accommodated by strike-slip faulting in the overriding continental lithosphere. The pattern of postseismic surface deformation due to viscoelastic flow in the lower crust and asthenosphere following a major earthquake on such a fault is modified from that predicted from the conventual elastic layer over viscoelastic halfspace model by the presence of the subducting slab. The predicted effects, such as a partial suppression of the postseismic velocities by 1 cm/yr or more immediately following a moderate to great earthquake, are potentially detectable using contemporary geodetic techniques.

Viscoelastic and elastomeric active matter: Linear instability and nonlinear dynamics.

PubMed

Hemingway, E J; Cates, M E; Fielding, S M

2016-03-01

We consider a continuum model of active viscoelastic matter, whereby an active nematic liquid crystal is coupled to a minimal model of polymer dynamics with a viscoelastic relaxation time τ(C). To explore the resulting interplay between active and polymeric dynamics, we first generalize a linear stability analysis (from earlier studies without polymer) to derive criteria for the onset of spontaneous heterogeneous flows (strain rate) and/or deformations (strain). We find two modes of instability. The first is a viscous mode, associated with strain rate perturbations. It dominates for relatively small values of τ(C) and is a simple generalization of the instability known previously without polymer. The second is an elastomeric mode, associated with strain perturbations, which dominates at large τ(C) and persists even as τ(C)→∞. We explore the dynamical states to which these instabilities lead by means of direct numerical simulations. These reveal oscillatory shear-banded states in one dimension and activity-driven turbulence in two dimensions even in the elastomeric limit τ(C)→∞. Adding polymer can also have calming effects, increasing the net throughput of spontaneous flow along a channel in a type of drag reduction. The effect of including strong antagonistic coupling between the nematic and polymer is examined numerically, revealing a rich array of spontaneously flowing states.

Viscoelastic and elastomeric active matter: Linear instability and nonlinear dynamics

NASA Astrophysics Data System (ADS)

Hemingway, E. J.; Cates, M. E.; Fielding, S. M.

2016-03-01

We consider a continuum model of active viscoelastic matter, whereby an active nematic liquid crystal is coupled to a minimal model of polymer dynamics with a viscoelastic relaxation time τC. To explore the resulting interplay between active and polymeric dynamics, we first generalize a linear stability analysis (from earlier studies without polymer) to derive criteria for the onset of spontaneous heterogeneous flows (strain rate) and/or deformations (strain). We find two modes of instability. The first is a viscous mode, associated with strain rate perturbations. It dominates for relatively small values of τC and is a simple generalization of the instability known previously without polymer. The second is an elastomeric mode, associated with strain perturbations, which dominates at large τC and persists even as τC→∞ . We explore the dynamical states to which these instabilities lead by means of direct numerical simulations. These reveal oscillatory shear-banded states in one dimension and activity-driven turbulence in two dimensions even in the elastomeric limit τC→∞ . Adding polymer can also have calming effects, increasing the net throughput of spontaneous flow along a channel in a type of drag reduction. The effect of including strong antagonistic coupling between the nematic and polymer is examined numerically, revealing a rich array of spontaneously flowing states.

Direct numerical simulation of particle alignment in viscoelastic fluids

NASA Astrophysics Data System (ADS)

Hulsen, Martien; Jaensson, Nick; Anderson, Patrick

2016-11-01

Rigid particles suspended in viscoelastic fluids under shear can align in string-like structures in flow direction. To unravel this phenomenon, we present 3D direct numerical simulations of the alignment of two and three rigid, non-Brownian particles in a shear flow of a viscoelastic fluid. The equations are solved on moving, boundary-fitted meshes, which are locally refined to accurately describe the polymer stresses around and in between the particles. A small minimal gap size between the particles is introduced. The Giesekus model is used and the effect of the Weissenberg number, shear thinning and solvent viscosity is investigated. Alignment of two and three particles is observed. Morphology plots have been created for various combinations of fluid parameters. Alignment is mainly governed by the value of the elasticity parameter S, defined as half of the ratio between the first normal stress difference and shear stress of the suspending fluid. Alignment appears to occur above a critical value of S, which decreases with increasing shear thinning. This result, together with simulations of a shear-thinning Carreau fluid, leads us to the conclusion that normal stress differences are essential for particle alignment to occur, but it is also strongly promoted by shear thinning.

Optical method of caustics applied in viscoelastic fracture analysis

NASA Astrophysics Data System (ADS)

Gao, Guiyun; Li, Zheng; Xu, Jie

2011-05-01

The optical method of caustics is developed here to study the fracture of viscoelastic materials. By adopting a distribution of viscoelastic stress fields near the crack tip, the method of caustics is used to determine the viscoelastic fracture parameters from the caustic patterns near the crack tip. Two viscoelastic materials are studied. These are PMMA and ternary composites of HDPE/POE-g-MA/CaCO 3. The transmitted and reflective methods of caustics are performed separately to investigate viscoelastic fracture behaviors. The stress intensity factors (SIFs) versus time is determined by a series of shadow spot patterns combined with viscoelastic parameters evaluated by creep tests. In order to understand the viscoelastic fracture mechanisms of HDPE/POE-g-MA/CaCO 3 composites, their fracture surfaces are observed by a Scanning Electron Microscope (SEM). The results indicate that the method of caustics can be used to characterize the fracture behaviors of viscoelastic materials and further to optimize the design of polymer composites.

Numerical Simulations of Homogeneous Turbulence Using Lagrangian-Averaged Navier-Stokes Equations

NASA Technical Reports Server (NTRS)

Mohseni, Kamran; Shkoller, Steve; Kosovic, Branko; Marsden, Jerrold E.; Carati, Daniele; Wray, Alan; Rogallo, Robert

2000-01-01

The Lagrangian-averaged Navier-Stokes (LANS) equations are numerically evaluated as a turbulence closure. They are derived from a novel Lagrangian averaging procedure on the space of all volume-preserving maps and can be viewed as a numerical algorithm which removes the energy content from the small scales (smaller than some a priori fixed spatial scale alpha) using a dispersive rather than dissipative mechanism, thus maintaining the crucial features of the large scale flow. We examine the modeling capabilities of the LANS equations for decaying homogeneous turbulence, ascertain their ability to track the energy spectrum of fully resolved direct numerical simulations (DNS), compare the relative energy decay rates, and compare LANS with well-accepted large eddy simulation (LES) models.

Recent progress of particle migration in viscoelastic fluids.

PubMed

Yuan, Dan; Zhao, Qianbin; Yan, Sheng; Tang, Shi-Yang; Alici, Gursel; Zhang, Jun; Li, Weihua

2018-02-13

Recently, research on particle migration in non-Newtonian viscoelastic fluids has gained considerable attention. In a viscoelastic fluid, three dimensional (3D) particle focusing can be easily realized in simple channels without the need for any external force fields or complex microchannel structures compared with that in a Newtonian fluid. Due to its promising properties for particle precise focusing and manipulation, this field has been developed rapidly, and research on the field has been shifted from fundamentals to applications. This review will elaborate the recent progress of particle migration in viscoelastic fluids, especially on the aspect of applications. The hydrodynamic forces on the micro/nano particles in viscoelastic fluids are discussed. Next, we elaborate the basic particle migration in viscoelasticity-dominant fluids and elasto-inertial fluids in straight channels. After that, a comprehensive review on the applications of viscoelasticity-induced particle migration (particle separation, cell deformability measurement and alignment, particle solution exchange, rheometry-on-a-chip and others) is presented; finally, we thrash out some perspectives on the future directions of particle migration in viscoelastic fluids.

Eulerian and Lagrangian approaches to multidimensional condensation and collection

NASA Astrophysics Data System (ADS)

Li, Xiang-Yu; Brandenburg, A.; Haugen, N. E. L.; Svensson, G.

2017-06-01

Turbulence is argued to play a crucial role in cloud droplet growth. The combined problem of turbulence and cloud droplet growth is numerically challenging. Here an Eulerian scheme based on the Smoluchowski equation is compared with two Lagrangian superparticle (or superdroplet) schemes in the presence of condensation and collection. The growth processes are studied either separately or in combination using either two-dimensional turbulence, a steady flow or just gravitational acceleration without gas flow. Good agreement between the different schemes for the time evolution of the size spectra is observed in the presence of gravity or turbulence. The Lagrangian superparticle schemes are found to be superior over the Eulerian one in terms of computational performance. However, it is shown that the use of interpolation schemes such as the cloud-in-cell algorithm is detrimental in connection with superparticle or superdroplet approaches. Furthermore, the use of symmetric over asymmetric collection schemes is shown to reduce the amount of scatter in the results. For the Eulerian scheme, gravitational collection is rather sensitive to the mass bin resolution, but not so in the case with turbulence.