MHD stagnation point flow and heat transfer of a nanofluid over a permeable nonlinear stretching/shrinking sheet with viscous dissipation effect

NASA Astrophysics Data System (ADS)

Jusoh, Rahimah; Nazar, Roslinda

2018-04-01

The magnetohydrodynamic (MHD) stagnation point flow and heat transfer of an electrically conducting nanofluid over a nonlinear stretching/shrinking sheet is studied numerically. Mathematical modelling and analysis are attended in the presence of viscous dissipation. Appropriate similarity transformations are used to reduce the boundary layer equations for momentum, energy and concentration into a set of ordinary differential equations. The reduced equations are solved numerically using the built in bvp4c function in Matlab. The numerical and graphical results on the effects of various parameters on the velocity and temperature profiles as well as the skin friction coefficient and the local Nusselt number are analyzed and discussed in this paper. The study discovers the existence of dual solutions for a certain range of the suction parameter. The conducted stability analysis reveals that the first solution is stable and feasible, while the second solution is unstable.

Thermally Radiative Rotating Magneto-Nanofluid Flow over an Exponential Sheet with Heat Generation and Viscous Dissipation: A Comparative Study

NASA Astrophysics Data System (ADS)

Sagheer, M.; Bilal, M.; Hussain, S.; Ahmed, R. N.

2018-03-01

This article examines a mathematical model to analyze the rotating flow of three-dimensional water based nanofluid over a convectively heated exponentially stretching sheet in the presence of transverse magnetic field with additional effects of thermal radiation, Joule heating and viscous dissipation. Silver (Ag), copper (Cu), copper oxide (CuO), aluminum oxide (Al 2 O 3 ) and titanium dioxide (TiO 2 ) have been taken under consideration as the nanoparticles and water (H 2 O) as the base fluid. Using suitable similarity transformations, the governing partial differential equations (PDEs) of the modeled problem are transformed to the ordinary differential equations (ODEs). These ODEs are then solved numerically by applying the shooting method. For the particular situation, the results are compared with the available literature. The effects of different nanoparticles on the temperature distribution are also discussed graphically and numerically. It is witnessed that the skin friction coefficient is maximum for silver based nanofluid. Also, the velocity profile is found to diminish for the increasing values of the magnetic parameter.

Wave-induced drift of large floating sheets

NASA Astrophysics Data System (ADS)

Christensen, K. H.; Weber, J. E.

In this article we study the wave-induced drift of large, flexible shallow floating objects, referred to as sheets. When surface waves propagate through a sheet, they provide a mean stress on the sheet, resulting in a mean drift. In response, the sheet generates an Ekman current. The drift velocity of the sheet is determined by (i) the wave-induced stress, (ii) the viscous stress due to the Ekman current, and (iii) the Coriolis force. The sheet velocity and the current beneath the sheet are determined for constant and depth-varying eddy viscosities.

A climate model with cryodynamics and geodynamics

NASA Technical Reports Server (NTRS)

Ghil, M.; Le Treut, H.

1981-01-01

A simplified, zero-dimensional model of the climatic system is presented which attempts to incorporate mechanisms important on the time scale of glaciation cycles: 10,000 to 100,000 years. The ocean-atmosphere radiation balance, continental ice sheet plastic flow, and upper mantle viscous flow are taken into account, with stress on the interaction between the ice sheets and the upper mantle. The model exhibits free, self-sustained oscillations of an amplitude and period comparable to those found in the paleoclimatic record of glaciations, offering mild support for the idea that unforced oscillations can actually exist in the real climatic system itself. The careful study of the interplay between internal mechanisms and external forcing is held to represent an interesting challenge to the theory of ice ages.

Influence of nonlinear thermal radiation and viscous dissipation on three-dimensional flow of Jeffrey nano fluid over a stretching sheet in the presence of Joule heating

NASA Astrophysics Data System (ADS)

Ganesh Kumar, K.; Rudraswamy, N. G.; Gireesha, B. J.; Krishnamurthy, M. R.

2017-09-01

Present exploration discusses the combined effect of viscous dissipation and Joule heating on three dimensional flow and heat transfer of a Jeffrey nanofluid in the presence of nonlinear thermal radiation. Here the flow is generated over bidirectional stretching sheet in the presence of applied magnetic field by accounting thermophoresis and Brownian motion of nanoparticles. Suitable similarity transformations are employed to reduce the governing partial differential equations into coupled nonlinear ordinary differential equations. These nonlinear ordinary differential equations are solved numerically by using the Runge-Kutta-Fehlberg fourth-fifth order method with shooting technique. Graphically results are presented and discussed for various parameters. Validation of the current method is proved by comparing our results with the existing results under limiting situations. It can be concluded that combined effect of Joule and viscous heating increases the temperature profile and thermal boundary layer thickness.

An orthotropic viscoelastic model for the passive myocardium: continuum basis and numerical treatment.

PubMed

Gültekin, Osman; Sommer, Gerhard; Holzapfel, Gerhard A

2016-11-01

This study deals with the viscoelastic constitutive modeling and the respective computational analysis of the human passive myocardium. We start by recapitulating the locally orthotropic inner structure of the human myocardial tissue and model the mechanical response through invariants and structure tensors associated with three orthonormal basis vectors. In accordance with recent experimental findings the ventricular myocardial tissue is assumed to be incompressible, thick-walled, orthotropic and viscoelastic. In particular, one spring element coupled with Maxwell elements in parallel endows the model with viscoelastic features such that four dashpots describe the viscous response due to matrix, fiber, sheet and fiber-sheet fragments. In order to alleviate the numerical obstacles, the strictly incompressible model is altered by decomposing the free-energy function into volumetric-isochoric elastic and isochoric-viscoelastic parts along with the multiplicative split of the deformation gradient which enables the three-field mixed finite element method. The crucial aspect of the viscoelastic formulation is linked to the rate equations of the viscous overstresses resulting from a 3-D analogy of a generalized 1-D Maxwell model. We provide algorithmic updates for second Piola-Kirchhoff stress and elasticity tensors. In the sequel, we address some numerical aspects of the constitutive model by applying it to elastic, cyclic and relaxation test data obtained from biaxial extension and triaxial shear tests whereby we assess the fitting capacity of the model. With the tissue parameters identified, we conduct (elastic and viscoelastic) finite element simulations for an ellipsoidal geometry retrieved from a human specimen.

A fiber-reinforced-fluid model of anisotropic plant root cell growth

NASA Astrophysics Data System (ADS)

Jensen, Oliver E.; Dyson, Rosemary J.

2009-11-01

We present a theoretical model of a single cell in the expansion zone of the primary root of the plant Arabidopsis thaliana. The cell undergoes rapid elongation with approximately constant radius. Growth is driven by high internal turgor pressure causing viscous stretching of the cell wall, with embedded cellulose microfibrils providing the wall with strongly anisotropic properties. We represent the cell as a thin cylindrical fiber-reinforced viscous sheet between rigid end plates. Asymptotic reduction of the governing equations, under simple sets of assumptions about fiber and wall properties, yields variants of the traditional Lockhart equation that relates the axial cell growth rate to the internal pressure. The model provides insights into the geometric and biomechanical parameters underlying bulk quantities such as wall extensibility and shows how either dynamical changes in wall material properties or passive fibre reorientation may suppress cell elongation.

Four large-scale field-aligned current systmes in the dayside high-latitude region

NASA Technical Reports Server (NTRS)

Ohtani, S.; Potemra, T. A.; Newell, P.T.; Zanetti, L. J.; Iijima, T.; Watanabe, M.; Blomberg, L. G.; Elphinstone, R. D.; Murphree, J. S.; Yamauchi, M.

1995-01-01

A system of four current sheets of large-scale field-aligned currents (FACs) was discovered in the data set of simultaneous Viking and Defense Meteorological Satellire Program-F7 (DMSP-F7) crossing of the dayside high-latitude region. This paper reports four examples of this system that were observed in the prenoon sector. The flow polarities of FACs are upward, downward, upward, and downward, from equatorward to poleward. The lowest-latitude upward current is flowing mostly in the central plasma sheet (CPS) precipitation region, often overlapping with the boundary plasma sheet (BPS) at its poleward edge, andis interpreted as a region 2 current. The pair of downward and upward FACs in the middle of te structure are collocated with structured electron precipitation. The precipitation of high-energy (greater than 1 keV) electrons is more intense in the lower-latitude downward current sheet. The highest-latitude downward flowing current sheet is located in a weak, low-energy particle precipitation region, suggesting that this current is flowing on open field lines. Simulaneous observations in the postnoon local time sector reveal the standard three-sheet structure of FACs, sometimes described as region 2, region 1, and mantle (referred to the midday region O) currents. A high correlation was found between the occurrence of the four FAC sheet structure and negative interplanetary magnetic field (IMF) B(sub Y). We discuss the FAC structurein terms of three types of convection cells: the merging, viscous, andlobe cells. During strongly negative IMF B(sub Y), two convection reversals exist in the prenoon sector; one is inside the viscous cell, and the other is between the viscous cell and the lobe cell. This structure of convection flow is supported by the Viking electric field and auroral UV image data. Based on the convection pattern, the four FAC sheet structure is interpreted as the latitude overlap of midday and morning FAC systems. We suggest that the for-current sheet structure is common in a certain prenoon localtime sector during strongly negative IMF B(sub Y).

Oblique impact of dense granular sheets

NASA Astrophysics Data System (ADS)

Ellowitz, Jake; Guttenberg, Nicholas; Jaeger, Heinrich M.; Nagel, Sidney R.; Zhang, Wendy W.

2013-11-01

Motivated by experiments showing impacts of granular jets with non-circular cross sections produce thin ejecta sheets with anisotropic shapes, we study what happens when two sheets containing densely packed, rigid grains traveling at the same speed collide asymmetrically. Discrete particle simulations and a continuum frictional fluid model yield the same steady-state solution of two exit streams emerging from incident streams. When the incident angle Δθ is less than Δθc =120° +/-10° , the exit streams' angles differ from that measured in water sheet experiments. Below Δθc , the exit angles from granular and water sheet impacts agree. This correspondence is surprising because 2D Euler jet impact, the idealization relevant for both situations, is ill posed: a generic Δθ value permits a continuous family of solutions. Our finding that granular and water sheet impacts evolve into the same member of the solution family suggests previous proposals that perturbations such as viscous drag, surface tension or air entrapment select the actual outcome are not correct. Currently at Department of Physics, University of Oregon, Eugene, OR 97403.

Long discontinuous fiber composite structure: Forming and structural mechanics

NASA Technical Reports Server (NTRS)

Pipes, R. B.; Santare, M. H.; Otoole, B. J.; Beaussart, A. J.; Deheer, D. C.; Okine, R. K.

1991-01-01

Cost effective composite structure has motivated the investigation of several new approaches to develop composite structure from innovative material forms. Among the promising new approaches is the conversion of planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. In the present study, the authors have established a framework which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. The initial study focuses upon the establishment of micromechanics models for prediction of the effective anisotropic viscosities of the oriented fiber assembly in a viscous matrix. Next, the developed constitutive relation is employed through an analogy with incompressible elasticity to exercise the finite element technique for determination of local fiber orientation and laminate thickness after forming. Results are presented for the stretch bending of a curved beam from an arbitrary composite laminate and the bulging of a clamped sheet. Structural analyses are conducted to determine the effect of microstructure on the performance of curved beams manufactured from long discontinuous fiber composites. For the purposes of this study, several curved beams with ideal and non-ideal microstructures are compared for response under pure bending. Material parameters are determined from a separate microstructural analysis.

Combined effects of slip and convective boundary condition on MHD 3D stretched flow of nanofluid through porous media inspired by non-linear thermal radiation

NASA Astrophysics Data System (ADS)

Nayak, M. K.; Shaw, Sachin; Pandey, V. S.; Chamkha, Ali J.

2018-02-01

In the present study, the main concern is to investigate the magnetohydrodynamic nanofluid flow subject to porous matrix and convective heating past a permeable linear stretching sheet. In addition, the influence of velocity slip, viscous dissipation, Joule heating and non-linear thermal radiation are considered. A new micro-convection model known as the Patel model is implemented for considerable enhancement of the thermal conductivity and hence, the heat transfer capability of nanofluids. Moreover, a convective heat transfer model is introduced where the bottom surface of the sheet gets heated due to a convection mechanism from a hot fluid of particular temperature. The numerical results of the transformed governing differential equations have been obtained by using fourth-order Runge-Kutta method along with shooting approach and secant method is used for better approximation. In the present analysis, base fluids such as water and Ethylene glycol and Copper, Silver and Aluminum oxide nanoparticles are considered. Results of the present investigation show that inclusion of porous matrix contributes to slow down the fluid velocity and diminution of wall shear stress (axial as well as transverse). Drag force due to magnetic field strength, velocity slip and imposed fluid suction impede the fluid motion and upsurge the heat transfer rate from the surface. In addition, rise in viscous dissipation widens the thermal boundary layer.

MHD Jeffrey nanofluid past a stretching sheet with viscous dissipation effect

NASA Astrophysics Data System (ADS)

Zokri, S. M.; Arifin, N. S.; Salleh, M. Z.; Kasim, A. R. M.; Mohammad, N. F.; Yusoff, W. N. S. W.

2017-09-01

This study investigates the influence of viscous dissipation on magnetohydrodynamic (MHD) flow of Jeffrey nanofluid over a stretching sheet with convective boundary conditions. The nonlinear partial differential equations are reduced into the nonlinear ordinary differential equations by utilizing the similarity transformation variables. The Runge-Kutta Fehlberg method is used to solve the problem numerically. The numerical solutions obtained are presented graphically for several dimensionless parameters such as Brownian motion, Lewis number and Eckert number on the specified temperature and concentration profiles. It is noted that the temperature profile is accelerated due to increasing values of Brownian motion parameter and Eckert number. In contrast, both the Brownian motion parameter and Lewis number have caused the deceleration in the concentration profiles.

Large-scale deformation associated with ridge subduction

USGS Publications Warehouse

Geist, E.L.; Fisher, M.A.; Scholl, D. W.

1993-01-01

Continuum models are used to investigate the large-scale deformation associated with the subduction of aseismic ridges. Formulated in the horizontal plane using thin viscous sheet theory, these models measure the horizontal transmission of stress through the arc lithosphere accompanying ridge subduction. Modelling was used to compare the Tonga arc and Louisville ridge collision with the New Hebrides arc and d'Entrecasteaux ridge collision, which have disparate arc-ridge intersection speeds but otherwise similar characteristics. Models of both systems indicate that diffuse deformation (low values of the effective stress-strain exponent n) are required to explain the observed deformation. -from Authors

Grain-scale alignment of melt in sheared partially molten rocks: implications for viscous anisotropy

NASA Astrophysics Data System (ADS)

Pec, Matej; Quintanilla-Terminel, Alejandra; Holtzman, Benjamin; Zimmerman, Mark; Kohlstedt, David

2016-04-01

Presence of melt significantly influences rheological properties of partially molten rocks by providing fast diffusional pathways. Under stress, melt aligns at the grain scale and this alignment induces viscous anisotropy in the deforming aggregate. One of the consequences of viscous anisotropy is melt segregation into melt-rich sheets oriented at low angle to the shear plane on much larger scales than the grain scale. The magnitude and orientation of viscous anisotropy with respect to the applied stress are important parameters for constitutive models (Takei and Holtzman 2009) that must be constrained by experimental studies. In this contribution, we analyze the shape preferred orientation (SPO) of individual grain-scale melt pockets in deformed partially molten mantle rocks. The starting materials were obtained by isostatically hot-pressing olivine + basalt and olivine + chromite + basalt powders. These partially molten rocks were deformed in general shear or torsion at a confining pressure, Pc = 300 MPa, temperature, T = 1200° - 1250° C, and strain rates of 10-3 - 10-5 s-1to finite shear strains, γ, of 0.5 - 5. After the experiment, high resolution backscattered electron images were obtained using a SEM equipped with a field emission gun. Individual melt pockets were segmented and their SPO analyzed using the paror and surfor methods and Fourier transforms (Heilbronner and Barret 2014). Melt segregation into melt-rich sheets inclined at 15° -20° antithetic with respect to the shear plane occurs in three-phase system (olivine + chromite + basalt) and in two-phase systems (olivine + basalt) twisted to high strain. The SPO of individual melt pockets within the melt-rich bands is moderately strong (b/a ≈ 0.8) and is always steeper (20° -40°) than the average melt-rich band orientation. In the two-phase system (olivine + basalt) sheared to lower strains, no distinct melt-rich sheets are observed. Individual grain-scale melt pockets are oriented at 45° -55° antithetic with respect to the shear plane (i.e., sub-perpendicular to σ3) with a strong SPO (b/a ≈ 0.7) that decreases with increasing finite strain. Our observations of melt alignment at low strains are in agreement with observations performed on analogue materials (borneol, Takei 2010) and provide further constraints for the orientation of viscous anisotropy in the Earth's mantle. The systematic difference in grain-scale melt alignment between samples in which melt segregation did and did not occur - irrespective of the deformation geometry and mineralogy - suggests that melt segregation into bands leads to local stress rotation within the samples.

Wrinkling instability of an inhomogeneously stretched viscous sheet

NASA Astrophysics Data System (ADS)

Srinivasan, Siddarth; Wei, Zhiyan; Mahadevan, L.

2017-07-01

Motivated by the redrawing of hot glass into thin sheets, we investigate the shape and stability of a thin viscous sheet that is inhomogeneously stretched in an imposed nonuniform temperature field. We first determine the associated base flow by solving the long-time-scale stretching flow of a flat sheet as a function of two dimensionless parameters: the normalized stretching velocity α and a dimensionless width of the heating zone β . This allows us to determine the conditions for the onset of an out-of-plane wrinkling instability stated in terms of an eigenvalue problem for a linear partial differential equation governing the displacement of the midsurface of the sheet. We show that the sheet can become unstable in two regions that are upstream and downstream of the heating zone where the minimum in-plane stress is negative. This yields the shape and growth rates of the most unstable buckling mode in both regions for various values of the stretching velocity and heating zone width. A transition from stationary to oscillatory unstable modes is found in the upstream region with increasing β , while the downstream region is always stationary. We show that the wrinkling instability can be entirely suppressed when the surface tension is large enough relative to the magnitude of the in-plane stress. Finally, we present an operating diagram that indicates regions of the parameter space that result in a required outlet sheet thickness upon stretching while simultaneously minimizing or suppressing the out-of-plane buckling, a result that is relevant for the glass redraw method used to create ultrathin glass sheets.

Sea-level and solid-Earth deformation feedbacks in ice sheet modelling

NASA Astrophysics Data System (ADS)

Konrad, Hannes; Sasgen, Ingo; Klemann, Volker; Thoma, Malte; Grosfeld, Klaus; Martinec, Zdeněk

2014-05-01

The interactions of ice sheets with the sea level and the solid Earth are important factors for the stability of the ice shelves and the tributary inland ice (e.g. Thomas and Bentley, 1978; Gomez et al, 2012). First, changes in ice extent and ice thickness induce viscoelastic deformation of the Earth surface and Earth's gravity field. In turn, global and local changes in sea level and bathymetry affect the grounding line and, subsequently, alter the ice dynamic behaviour. Here, we investigate these feedbacks for a synthetic ice sheet configuration as well as for the Antarctic ice sheet using a three-dimensional thermomechanical ice sheet and shelf model, coupled to a viscoelastic solid-Earth and gravitationally self-consistent sea-level model. The respective ice sheet undergoes a forcing from rising sea level, warming ocean, and/or changing surface mass balance. The coupling is realized by exchanging ice thickness, Earth surface deformation and sea level periodically. We apply several sets of viscoelastic Earth parameters to our coupled model, e.g. simulating a low-viscous upper mantle present at the Antarctic Peninsula (Ivins et al., 2011). Special focus of our study lies on the evolution of Earth surface deformation and local sea level changes, as well as on the accompanying grounding line evolution. N. Gomez, D. Pollard, J. X. Mitrovica, P. Huybers, and P. U. Clark 2012. Evolution of a coupled marine ice sheet-sea level model, J. Geophys. Res., 117, F01013, doi:10.1029/2011JF002128. E. R. Ivins, M. M. Watkins, D.-N. Yuan, R. Dietrich, G. Casassa, and A. Rülke 2011. On-land ice loss and glacial isostatic adjustment at the Drake Passage: 2003-2009, J. Geophys. Res. 116, B02403, doi: 10.1029/2010JB007607 R. H. Thomas and C. R. Bentley 1978. A model for Holocene retreat of the West Antarctic Ice Sheet, Quaternary Research, 10 (2), pages 150-170, doi: 10.1016/0033-5894(78)90098-4.

Analogue modeling of 3-D structural segmentation in fold-and-thrust belts: interactions between frictional and viscous provinces in foreland basins

NASA Astrophysics Data System (ADS)

Borderie, Sandra; Graveleau, Fabien; Witt, César; Vendeville, Bruno C.

2016-04-01

Accretionary wedges are generally segmented both across and along strike because of diverse factors including tectonic and stratigraphic inheritance. In fold-and-thrust belts, along-strike stratigraphic changes in the foreland sequence are classically observed and cause a curvature of the deformation front. Although the parameters controlling this curvature are well documented, the structural interactions and mutual influences between adjacent provinces are much less analyzed. To investigate this question, we deformed analogue models in a compressional box equipped with digital cameras and a topographic measurement apparatus. Models where shortened above a basal frictional detachment (glass microbeads) and segmentation was tested by having a region in which we added an interbedded viscous level (silicone polymer) within the sedimentary cover (dry sand). By changing the number (2 or 3) and the relative width of the purely frictional and viscous provinces, our goal was to characterize geometrically and kinematically the interactions between the viscous and the purely frictional provinces. We used a commercial geomodeller to generate 3-D geometrical models. The results indicate that regardless of the relative width of the purely frictional vs. viscous provinces, the deformation style in the frictional province is not influenced by the presence of the adjacent viscous province. On the contrary, the structural style and the deformation kinematics in the viscous province is significantly impacted by the presence or absence of an adjacent purely frictional province. At first order, the deformation style in the viscous province depends on its width, and three structural styles can be defined along strike. Far from the frictional area, structures are primarily of salt-massif type, and they do not seem to be influenced by the frictional wedge province. Towards the frictional province, deformation changes gradually to a zone of purely forethrusts (foreland verging), and finally to a highly faulted zone with both fore- and backthrusts (hinterland verging). In addition, a kinematic analysis indicates that narrow viscous provinces are strongly influenced by the presence of an adjacent frictional province. Indeed, propagation of shallow thrusts occurs in sequence and the deformation front reaches lately the external décollement pinchout. On the contrary, the deformation front of the wide viscous provinces propagates rapidly to the external décollement pinchout, then younger thrusts form out of sequence. Along-strike segmentation also affects the deep structures (thrusts detaching on the basal frictional décollement). In the viscous province, the presence of an upper viscous décollement opposes the advance of the basal deformation front. There, the rear of the wedge is characterized by imbrications of thrusts sheets (antiformal stacks), and the deep deformation front is convex towards the hinterland. Our experiments allow to better understand the dynamics of salt-controlled fold-and-thrust belts such as in the Huallaga (Peru) and Kuqa (China) basins or the Franklin Mountains (NW Canada).

The development of a mixing layer under the action of weak streamwise vortices

NASA Technical Reports Server (NTRS)

Goldstein, Marvin E.; Mathew, Joseph

1993-01-01

The action of weak, streamwise vortices on a plane, incompressible, steady mixing layer is examined in the large Reynolds-number limit. The outer, inviscid region is bounded by a vortex sheet to which the viscous region is confined. It is shown that the local linear analysis becomes invalid at streamwise distances O(epsilon(sup -1)), where epsilon is much less than 1 is the cross flow amplitude, and a new nonlinear analysis is constructed for this region. Numerical solutions of the nonlinear problem show that the vortex sheet undergoes an O(1) change in position and that the solution is ultimately terminated by the appearance of a singularity. The corresponding viscous layer shows downstream thickening, but appears to remain well behaved up to the singular location.

The development of a mixing layer under the action of weak streamwise vortices

NASA Technical Reports Server (NTRS)

Goldstein, M. E.; Mathew, Joseph

1993-01-01

The action of weak, streamwise vortices on a plane, incompressible, steady mixing layer is examined in the large Reynolds number limit. The outer, inviscid region is bounded by a vortex sheet to which the viscous region is confined. It is shown that the local linear analysis becomes invalid at streamwise distances O(epsilon sup -1), where (epsilon much less than 1) is the crossflow amplitude, and a new nonlinear analysis is constructed for this region. Numerical solutions of the nonlinear problem show that the vortex sheet undergoes an O(1) change in position and that the solution is ultimately terminated by a breakdown in the numerical procedure. The corresponding viscous layer shows downstream thickening, but appears to remain well behaved up to the terminal location.

Stretching a Curved Surface in a Viscous Fluid

NASA Astrophysics Data System (ADS)

Sajid, M.; N., Ali; T., Javed; Z., Abbas

2010-02-01

This work is concerned with the viscous flow due to a curved stretching sheet. The similarity solution of the problem is obtained numerically by a shooting method using the Runge-Kutta algorithm. The physical quantities of interest like the fluid velocity and skin friction coefficient are obtained and discussed under the influence of dimensionless curvature. It is evident from the results that dimensionless curvature causes an increase in boundary layer thickness and a decrease in the skin friction coefficient.

A model for spiral flows in basal ice and the formation of subglacial flutes based on a Reiner-Rivlin rheology for glacial ice

NASA Astrophysics Data System (ADS)

Schoof, Christian G.; Clarke, Garry K. C.

2008-05-01

Flutes are elongated sediment ridges formed at the base of glaciers and ice sheets. In this paper, we show that flutes can be the product of a corkscrew-like spiral flow in basal ice that removes sediment from troughs between flutes and deposits it at their crests, as first suggested by Shaw and Freschauf. In order to generate the type of basal ice flow required for this mechanism, the viscous rheology of ice must allow for the generation of deviatoric normal stresses transverse to the main flow direction. This type of behavior, which is commonly observed in real nonlinearly viscous and viscoelastic fluids, can be described by a Reiner-Rivlin rheology. Here, we develop a mathematical model that describes the role of these transverse stresses in generating spiral flows in basal ice and investigate how these flows lead to the amplification of initially small basal topography and the eventual formation of assemblies of evenly spaced subglacial flutes.

Effects of heat and mass transfer on unsteady boundary layer flow of a chemical reacting Casson fluid

NASA Astrophysics Data System (ADS)

Khan, Kashif Ali; Butt, Asma Rashid; Raza, Nauman

2018-03-01

In this study, an endeavor is to observe the unsteady two-dimensional boundary layer flow with heat and mass transfer behavior of Casson fluid past a stretching sheet in presence of wall mass transfer by ignoring the effects of viscous dissipation. Chemical reaction of linear order is also invoked here. Similarity transformation have been applied to reduce the governing equations of momentum, energy and mass into non-linear ordinary differential equations; then Homotopy analysis method (HAM) is applied to solve these equations. Numerical work is done carefully with a well-known software MATHEMATICA for the examination of non-dimensional velocity, temperature, and concentration profiles, and then results are presented graphically. The skin friction (viscous drag), local Nusselt number (rate of heat transfer) and Sherwood number (rate of mass transfer) are discussed and presented in tabular form for several factors which are monitoring the flow model.

Sensitivity of grounding line dynamics to viscoelastic deformation of the solid Earth: Inferences from a fully coupled ice sheet - solid Earth model

NASA Astrophysics Data System (ADS)

Konrad, H.; Sasgen, I.; Thoma, M.; Klemann, V.; Grosfeld, K.; Martinec, Z.

2013-12-01

The interactions of ice sheets with the sea level and the solid Earth are important factors for the stability of the ice shelves and the tributary inland ice (e.g. Thomas and Bentley, 1978; Gomez et al, 2012). First, changes in ice extent and ice thickness induce viscoelastic deformation of the Earth surface and Earth's gravity field. In turn, global and local changes in sea level and bathymetry affect the grounding line and, subsequently, alter the ice dynamic behaviour. Here, we investigate these feedbacks for a synthetic ice sheet configuration as well as for the Antarctic ice sheet using a three-dimensional thermomechanical ice sheet and shelf model, coupled to a viscoelastic solid-Earth and gravitationally self-consistent sea-level model. The respective ice sheet undergoes a forcing from rising sea level, warming ocean, and/or changing surface mass balance. The coupling is realized by exchanging ice thickness, Earth surface deformation, and sea level periodically. We apply several sets of viscoelastic Earth parameters to our coupled model, e.g. simulating a low-viscous upper mantle present at the Antarctic Peninsula (Ivins et al., 2011). Special focus of our study lies on the evolution of Earth surface deformation and local sea level changes, as well as on the accompanying grounding line evolution. N. Gomez, D. Pollard, J. X. Mitrovica, P. Huybers, and P. U. Clark 2012. Evolution of a coupled marine ice sheet-sea level model, J. Geophys. Res., 117, F01013, doi:10.1029/2011JF002128. E. R. Ivins, M. M. Watkins, D.-N. Yuan, R. Dietrich, G. Casassa, and A. Rülke 2011. On-land ice loss and glacial isostatic adjustment at the Drake Passage: 2003-2009, J. Geophys. Res. 116, B02403, doi: 10.1029/2010JB007607 R. H. Thomas and C. R. Bentley 1978. A model for Holocene retreat of the West Antarctic Ice Sheet, Quaternary Research, 10 (2), pages 150-170, doi: 10.1016/0033-5894(78)90098-4.

Method for shaping sheet thermoplastic and the like

NASA Technical Reports Server (NTRS)

Akilian, Mireille K. (Inventor); Schattenburg, Mark L. (Inventor)

2011-01-01

Processes and apparati for shaping sheet glass or thermoplastic materials use force from a layer of a flowing fluid, such as air, between the sheet and a mandrel at close to the softening temperature of the thermoplastic. The shape is preserved by cooling. The shape of the air bearing mandrel and the pressure distribution of the fluid contribute to the final shape. A process can be conducted on one or two surfaces such that the force from the air layer is on one or two surfaces of the sheet. The gap size between the sheet and mandrel determines the pressure profile in the gap, which also determines the final sheet shape. In general, smaller gaps lead to larger viscous forces. The pressure profile depends on the shape of the mandrel, the size of the fluid gap and the sheet and the fluid supply pressure.

An analysis method for multi-component airfoils in separated flow

NASA Technical Reports Server (NTRS)

Rao, B. M.; Duorak, F. A.; Maskew, B.

1980-01-01

The multi-component airfoil program (Langley-MCARF) for attached flow is modified to accept the free vortex sheet separation-flow model program (Analytical Methods, Inc.-CLMAX). The viscous effects are incorporated into the calculation by representing the boundary layer displacement thickness with an appropriate source distribution. The separation flow model incorporated into MCARF was applied to single component airfoils. Calculated pressure distributions for angles of attack up to the stall are in close agreement with experimental measurements. Even at higher angles of attack beyond the stall, correct trends of separation, decrease in lift coefficients, and increase in pitching moment coefficients are predicted.

Near wakes of advanced turbopropellers

NASA Technical Reports Server (NTRS)

Hanson, D. B.; Patrick, W. P.

1989-01-01

The flow in the wake of a model single rotation Prop-Fan rotor operating in a wind tunnel was traversed with a hot-wire anemometer system designed to determine the 3 periodic velocity components. Special data acquisition and data reduction methods were required to deal with the high data frequency, narrow wakes, and large fluctuating air angles in the tip vortex region. The model tip helical Mach number was 1.17, simulating the cruise condition. Although the flow field is complex, flow features such as viscous velocity defects, vortex sheets, tip vortices, and propagating acoustic pulses are clearly identified with the aid of a simple analytical wake theory.

Intrusion of Magmatic Bodies Into the Continental Crust: 3-D Numerical Models

NASA Astrophysics Data System (ADS)

Gorczyk, Weronika; Vogt, Katharina

2018-03-01

Magma intrusion is a major material transfer process in the Earth's continental crust. Yet the mechanical behavior of the intruding magma and its host are a matter of debate. In this study we present a series of numerical thermomechanical simulations on magma emplacement in 3-D. Our results demonstrate the response of the continental crust to magma intrusion. We observe change in intrusion geometries between dikes, cone sheets, sills, plutons, ponds, funnels, finger-shaped and stock-like intrusions, and injection time. The rheology and temperature of the host are the main controlling factors in the transition between these different modes of intrusion. Viscous deformation in the warm and deep crust favors host rock displacement and plutons at the crust-mantle boundary forming deep-seated plutons or magma ponds in the lower to middle crust. Brittle deformation in the cool and shallow crust induces cone-shaped fractures in the host rock and enables emplacement of finger- or stock-like intrusions at shallow or intermediate depth. Here the passage of magmatic and hydrothermal fluids from the intrusion through the fracture pattern may result in the formation of ore deposits. A combination of viscous and brittle deformation forms funnel-shaped intrusions in the middle crust. Intrusion of low-density magma may more over result in T-shaped intrusions in cross section with magma sheets at the surface.

Modeling ocean wave propagation under sea ice covers

NASA Astrophysics Data System (ADS)

Zhao, Xin; Shen, Hayley H.; Cheng, Sukun

2015-02-01

Operational ocean wave models need to work globally, yet current ocean wave models can only treat ice-covered regions crudely. The purpose of this paper is to provide a brief overview of ice effects on wave propagation and different research methodology used in studying these effects. Based on its proximity to land or sea, sea ice can be classified as: landfast ice zone, shear zone, and the marginal ice zone. All ice covers attenuate wave energy. Only long swells can penetrate deep into an ice cover. Being closest to open water, wave propagation in the marginal ice zone is the most complex to model. The physical appearance of sea ice in the marginal ice zone varies. Grease ice, pancake ice, brash ice, floe aggregates, and continuous ice sheet may be found in this zone at different times and locations. These types of ice are formed under different thermal-mechanical forcing. There are three classic models that describe wave propagation through an idealized ice cover: mass loading, thin elastic plate, and viscous layer models. From physical arguments we may conjecture that mass loading model is suitable for disjoint aggregates of ice floes much smaller than the wavelength, thin elastic plate model is suitable for a continuous ice sheet, and the viscous layer model is suitable for grease ice. For different sea ice types we may need different wave ice interaction models. A recently proposed viscoelastic model is able to synthesize all three classic models into one. Under suitable limiting conditions it converges to the three previous models. The complete theoretical framework for evaluating wave propagation through various ice covers need to be implemented in the operational ocean wave models. In this review, we introduce the sea ice types, previous wave ice interaction models, wave attenuation mechanisms, the methods to calculate wave reflection and transmission between different ice covers, and the effect of ice floe breaking on shaping the sea ice morphology. Laboratory experiments, field measurements and numerical simulations supporting the fundamental research in wave-ice interaction models are discussed. We conclude with some outlook of future research needs in this field.

Observation of the magnetospheric ``sash'' and its implications relative to solar-wind/magnetospheric coupling: A multisatellite event analysis

NASA Astrophysics Data System (ADS)

Maynard, N. C.; Savin, S.; Erickson, G. M.; Kawano, H.; Němeček, Z.; Peterson, W. K.; Šafránoková, J.; Sandahl, I.; Scudder, J. D.; Siscoe, G. L.; Sonnerup, B. U. Ö.; Weimer, D. R.; White, W. W.; Wilson, G. R.

2001-04-01

Using a unique data set from the Wind, Polar, Interball 1, Magion 4, and Defense Meteorological Satellite Program (DMSP) F11 satellites, comparisons with the Integrated Space Weather Model (ISM) have provided validation of the global structure predicted by the ISM model, which in turn has allowed us to use the model to interpret the data to further understand boundary layers and magnetospheric processes. The comparisons have shown that the magnetospheric ``sash'' [White et al., 1998], a region of low magnetic field discovered by the MHD modeling which extends along the high-latitude flank of the magnetopause, is related to the turbulent boundary layer on the high-latitude magnetopause, recently mapped by Interball 1. The sash in the data and in the model has rotational discontinuity properties, expected for a reconnection site. At some point near or behind the terminator, the sash becomes a site for reconnection of open field lines, which were previously opened by merging on the dayside. This indicates that significant reconnection in the magnetotail occurs on the flanks. Polar mapped to the high-density extension of the sash into the tilted plasma sheet. The source of the magnetosheath plasma observed by Polar on closed field lines behind the terminator was plasma entry through the low field connection of the sash to the central plasma sheet. The Polar magnetic field line footprints in each hemisphere are moving in different directions. Above and below the tilted plasma sheet the flows in the model are consistent with the corresponding flows in the ionosphere. The turbulence in the plasma sheet allows the convection patterns from each hemisphere to adjust. The boundary layer in the equatorial plane on the flank for this interplanetary magnetic field BY condition, which is below the tilted central plasma sheet, is several RE thick and is on tailward flowing open field lines. This thick boundary layer shields the magnetopause from viscous forces and must be driven by magnetic tension. Above the plasma sheet the boundary layer is dominated by the sash, and the model indicates that the open region inside the sash is considerably thinner.

Study of velocity and temperature distributions in boundary layer flow of fourth grade fluid over an exponential stretching sheet

NASA Astrophysics Data System (ADS)

Khan, Najeeb Alam; Saeed, Umair Bin; Sultan, Faqiha; Ullah, Saif; Rehman, Abdul

2018-02-01

This study deals with the investigation of boundary layer flow of a fourth grade fluid and heat transfer over an exponential stretching sheet. For analyzing two heating processes, namely, (i) prescribed surface temperature (PST), and (ii) prescribed heat flux (PHF), the temperature distribution in a fluid has been considered. The suitable transformations associated with the velocity components and temperature, have been employed for reducing the nonlinear model equation to a system of ordinary differential equations. The flow and temperature fields are revealed by solving these reduced nonlinear equations through an effective analytical method. The important findings in this analysis are to observe the effects of viscoelastic, cross-viscous, third grade fluid, and fourth grade fluid parameters on the constructed analytical expression for velocity profile. Likewise, the heat transfer properties are studied for Prandtl and Eckert numbers.

The Effects of Thermal Radiation on an Unsteady MHD Axisymmetric Stagnation-Point Flow over a Shrinking Sheet in Presence of Temperature Dependent Thermal Conductivity with Navier Slip

PubMed Central

Mondal, Sabyasachi; Haroun, Nageeb A. H.; Sibanda, Precious

2015-01-01

In this paper, the magnetohydrodynamic (MHD) axisymmetric stagnation-point flow of an unsteady and electrically conducting incompressible viscous fluid in with temperature dependent thermal conductivity, thermal radiation and Navier slip is investigated. The flow is due to a shrinking surface that is shrunk axisymmetrically in its own plane with a linear velocity. The magnetic field is imposed normally to the sheet. The model equations that describe this fluid flow are solved by using the spectral relaxation method. Here, heat transfer processes are discussed for two different types of wall heating; (a) a prescribed surface temperature and (b) a prescribed surface heat flux. We discuss and evaluate how the various parameters affect the fluid flow, heat transfer and the temperature field with the aid of different graphical presentations and tabulated results. PMID:26414006

High-performance graphene-based supercapacitors made by a scalable blade-coating approach

NASA Astrophysics Data System (ADS)

Wang, Bin; Liu, Jinzhang; Mirri, Francesca; Pasquali, Matteo; Motta, Nunzio; Holmes, John W.

2016-04-01

Graphene oxide (GO) sheets can form liquid crystals (LCs) in their aqueous dispersions that are more viscous with a stronger LC feature. In this work we combine the viscous LC-GO solution with the blade-coating technique to make GO films, for constructing graphene-based supercapacitors in a scalable way. Reduced GO (rGO) films are prepared by wet chemical methods, using either hydrazine (HZ) or hydroiodic acid (HI). Solid-state supercapacitors with rGO films as electrodes and highly conductive carbon nanotube films as current collectors are fabricated and the capacitive properties of different rGO films are compared. It is found that the HZ-rGO film is superior to the HI-rGO film in achieving high capacitance, owing to the 3D structure of graphene sheets in the electrode. Compared to gelled electrolyte, the use of liquid electrolyte (H2SO4) can further increase the capacitance to 265 F per gram (corresponding to 52 mF per cm2) of the HZ-rGO film.

The electromagnetic field for an open magnetosphere

NASA Technical Reports Server (NTRS)

Heikkila, W. J.

1984-01-01

The boundary-layer-dominated models of the earth EM field developed by Heikkila (1975, 1978, 1982, and 1983) and Heikkila et al. (1979) to account for deficiencies in the electric-field descriptions of quasi-steady-state magnetic-field-reconnection models (such as that of Cowley, 1980) are characterized, reviewing the arguments and indicating the most important implications. The mechanisms of boundary-layer formation and field direction reversal are explained and illustrated with diagrams, and it is inferred that boundary-layer phenomena rather than magnetic reconnection may be the cause of large-scale magnetospheric circulation, convection, plasma-sheet formation and sunward convection, and auroras, the boundary layer acting basically as a viscous process mediating solar-wind/magnetosphere interactions.

The multi-millennial Antarctic commitment to future sea-level rise

NASA Astrophysics Data System (ADS)

Golledge, N. R.; Kowalewski, D. E.; Naish, T. R.; Levy, R. H.; Fogwill, C. J.; Gasson, E. G. W.

2015-10-01

Atmospheric warming is projected to increase global mean surface temperatures by 0.3 to 4.8 degrees Celsius above pre-industrial values by the end of this century. If anthropogenic emissions continue unchecked, the warming increase may reach 8-10 degrees Celsius by 2300 (ref. 2). The contribution that large ice sheets will make to sea-level rise under such warming scenarios is difficult to quantify because the equilibrium-response timescale of ice sheets is longer than those of the atmosphere or ocean. Here we use a coupled ice-sheet/ice-shelf model to show that if atmospheric warming exceeds 1.5 to 2 degrees Celsius above present, collapse of the major Antarctic ice shelves triggers a centennial- to millennial-scale response of the Antarctic ice sheet in which enhanced viscous flow produces a long-term commitment (an unstoppable contribution) to sea-level rise. Our simulations represent the response of the present-day Antarctic ice-sheet system to the oceanic and climatic changes of four representative concentration pathways (RCPs) from the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. We find that substantial Antarctic ice loss can be prevented only by limiting greenhouse gas emissions to RCP 2.6 levels. Higher-emissions scenarios lead to ice loss from Antarctic that will raise sea level by 0.6-3 metres by the year 2300. Our results imply that greenhouse gas emissions in the next few decades will strongly influence the long-term contribution of the Antarctic ice sheet to global sea level.

The multi-millennial Antarctic commitment to future sea-level rise.

PubMed

Golledge, N R; Kowalewski, D E; Naish, T R; Levy, R H; Fogwill, C J; Gasson, E G W

2015-10-15

Atmospheric warming is projected to increase global mean surface temperatures by 0.3 to 4.8 degrees Celsius above pre-industrial values by the end of this century. If anthropogenic emissions continue unchecked, the warming increase may reach 8-10 degrees Celsius by 2300 (ref. 2). The contribution that large ice sheets will make to sea-level rise under such warming scenarios is difficult to quantify because the equilibrium-response timescale of ice sheets is longer than those of the atmosphere or ocean. Here we use a coupled ice-sheet/ice-shelf model to show that if atmospheric warming exceeds 1.5 to 2 degrees Celsius above present, collapse of the major Antarctic ice shelves triggers a centennial- to millennial-scale response of the Antarctic ice sheet in which enhanced viscous flow produces a long-term commitment (an unstoppable contribution) to sea-level rise. Our simulations represent the response of the present-day Antarctic ice-sheet system to the oceanic and climatic changes of four representative concentration pathways (RCPs) from the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. We find that substantial Antarctic ice loss can be prevented only by limiting greenhouse gas emissions to RCP 2.6 levels. Higher-emissions scenarios lead to ice loss from Antarctic that will raise sea level by 0.6-3 metres by the year 2300. Our results imply that greenhouse gas emissions in the next few decades will strongly influence the long-term contribution of the Antarctic ice sheet to global sea level.

Magma emplacement in 3D

NASA Astrophysics Data System (ADS)

Gorczyk, W.; Vogt, K.

2017-12-01

Magma intrusion is a major material transfer process in Earth's continental crust. Yet, the mechanical behavior of the intruding magma and its host are a matter of debate. In this study, we present a series of numerical thermo-mechanical experiments on mafic magma emplacement in 3D.In our model, we place the magmatic source region (40 km diameter) at the base of the mantle lithosphere and connect it to the crust by a 3 km wide channel, which may have evolved at early stages of magmatism during rapid ascent of hot magmatic fluids/melts. Our results demonstrate continental crustal response due to magma intrusion. We observe change in intrusion geometries between dikes, cone-sheets, sills, plutons, ponds, funnels, finger-shaped and stock-like intrusions as well as injection time. The rheology and temperature of the host-rock are the main controlling factors in the transition between these different modes of intrusion. Viscous deformation in the warm and deep crust favours host rock displacement and magma pools along the crust-mantle boundary forming deep-seated plutons or magma ponds in the lower to middle-crust. Brittle deformation in the cool and shallow crust induces cone-shaped fractures in the host rock and enables emplacement of finger- or stock-like intrusions at shallow or intermediate depth. A combination of viscous and brittle deformation forms funnel-shaped intrusions in the middle-crust. Low-density source magma results in T-shaped intrusions in cross-section with magma sheets at the surface.

Igneous intrusion models for floor fracturing in lunar craters

NASA Technical Reports Server (NTRS)

Wichman, R. W.; Schultz, P. H.

1991-01-01

Lunar floor-fractured craters are primarily located near the maria and frequently contain ponded mare units and dark mantling deposits. Fracturing is confined to the crater interior, often producing a moat-like feature near the floor edge, and crater depth is commonly reduced by uplift of the crater floor. Although viscous relaxation of crater topography can produce such uplift, the close association of modification with surface volcanism supports a model linking floor fracture to crater-centered igneous intrusions. The consequences of two intrusion models for the lunar interior are quantitatively explored. The first model is based on terrestrial laccoliths and describes a shallow intrusion beneath the crater. The second model is based on cone sheet complexes where surface deformation results from a deeper magma chamber. Both models, their fit to observed crater modifications and possible implications for local volcanism are described.

Instability of the Antarctic Ross Sea Embayment as climate warms

NASA Astrophysics Data System (ADS)

Hughes, Terence; Zhao, Zihong; Hintz, Raymond; Fastook, James

2017-06-01

Collapse of the Antarctic Ice Sheet since the Last Glacial Maximum 18,000 years ago is most pronounced in the Ross Sea Embayment, which is partly ice-free during Antarctic summers, thereby breaching the O-ring of ice shelves and sea ice surrounding Antarctica that stabilizes the ice sheet. The O-ring may have vanished during Early Holocene (5000 to 3000 B.C.), Roman (1 to 400 A.D.), and Medieval (900 to 1300 A.D.) warm periods and reappeared during the Little Ice Age (1300 to 1900 A.D.). We postulate further collapse in the embayment during the post-1900 warming may be forestalled because East Antarctic outlet glaciers "nail" the Ross Ice Shelf to the Transantarctic Mountains so it can resist the push from West Antarctic ice streams. Our hypothesis is examined for Byrd Glacier and a static ice shelf using three modeling experiments having plastic, viscous, and viscoplastic solutions as more data and improved modeling became available. Observed crevasse patterns were not reproduced. A new research study is needed to model a dynamic Ross Ice Shelf with all its feeder ice streams, outlet glaciers, and ice calving dynamics in three dimensions over time to fully test our hypothesis. The required model must allow accelerated calving if further warming melts sea ice and discerps the ice shelf. Calving must then successively pull the outlet glacier "nails" so collapse of the marine West Antarctic Ice Sheet proceeds to completion.

Fundamental Processes of Atomization in Fluid-Fluid Flows

NASA Technical Reports Server (NTRS)

McCready, M. J.; Chang, H.-C.; Leighton, D. T.

2001-01-01

This report outlines the major results of the grant "Fundamental Processes of Atomization in Fluid-Fluid Flows." These include: 1) the demonstration that atomization in liquid/liquid shear flow is driven by a viscous shear instability that triggers the formation of a long thin sheet; 2) discovery of a new mode of interfacial instability for oscillatory two-layer systems whereby a mode that originates within the less viscous liquid phase causes interfacial deformation as the oscillation proceeds; 3) the demonstration that rivulet formation from gravity front occurs because the local front shape specified by gravity and surface tension changes from a nose to a wedge geometry, thus triggering a large increase in viscous resistance; and 4) extension of the studies on nonlinear wave evolution on falling films and in stratified flow, particularly the evolution towards large-amplitude solitary waves that tend to generate drops.

Spray formation processes of impinging jet injectors

NASA Technical Reports Server (NTRS)

Anderson, W. E.; Ryan, H. M.; Pal, S.; Santoro, R. J.

1993-01-01

A study examining impinging liquid jets has been underway to determine physical mechanisms responsible for combustion instabilities in liquid bi-propellant rocket engines. Primary atomization has been identified as an important process. Measurements of atomization length, wave structure, and drop size and velocity distribution were made under various ambient conditions. Test parameters included geometric effects and flow effects. It was observed that pre-impingement jet conditions, specifically whether they were laminar or turbulent, had the major effect on primary atomization. Comparison of the measurements with results from a two dimensional linear aerodynamic stability model of a thinning, viscous sheet were made. Measured turbulent impinging jet characteristics were contrary to model predictions; the structure of waves generated near the point of jet impingement were dependent primarily on jet diameter and independent of jet velocity. It has been postulated that these impact waves are related to pressure and momentum fluctuations near the impingement region and control the eventual disintegration of the liquid sheet into ligaments. Examination of the temporal characteristics of primary atomization (ligament shedding frequency) strongly suggests that the periodic nature of primary atomization is a key process in combustion instability.

Arctic landfast sea ice

NASA Astrophysics Data System (ADS)

Konig, Christof S.

Landfast ice is sea ice which forms and remains fixed along a coast, where it is attached either to the shore, or held between shoals or grounded icebergs. Landfast ice fundamentally modifies the momentum exchange between atmosphere and ocean, as compared to pack ice. It thus affects the heat and freshwater exchange between air and ocean and impacts on the location of ocean upwelling and downwelling zones. Further, the landfast ice edge is essential for numerous Arctic mammals and Inupiat who depend on them for their subsistence. The current generation of sea ice models is not capable of reproducing certain aspects of landfast ice formation, maintenance, and disintegration even when the spatial resolution would be sufficient to resolve such features. In my work I develop a new ice model that permits the existence of landfast sea ice even in the presence of offshore winds, as is observed in mature. Based on viscous-plastic as well as elastic-viscous-plastic ice dynamics I add tensile strength to the ice rheology and re-derive the equations as well as numerical methods to solve them. Through numerical experiments on simplified domains, the effects of those changes are demonstrated. It is found that the modifications enable landfast ice modeling, as desired. The elastic-viscous-plastic rheology leads to initial velocity fluctuations within the landfast ice that weaken the ice sheet and break it up much faster than theoretically predicted. Solving the viscous-plastic rheology using an implicit numerical method avoids those waves and comes much closer to theoretical predictions. Improvements in landfast ice modeling can only verified in comparison to observed data. I have extracted landfast sea ice data of several decades from several sources to create a landfast sea ice climatology that can be used for that purpose. Statistical analysis of the data shows several factors that significantly influence landfast ice distribution: distance from the coastline, ocean depth, as well as the strength of offshore winds during nine out of the twelve months each year. Additionally, I identify regions where landfast ice appearance has been increasing or decreasing over the observed time span.

The multi-millennial Antarctic commitment to future sea-level rise

NASA Astrophysics Data System (ADS)

Golledge, Nicholas R.; Kowalewski, Douglas E.; Naish, Timothy R.; Levy, Richard H.; Fogwill, Christopher J.; Gasson, Edward G. W.

2016-04-01

Atmospheric warming is projected to increase global mean surface temperatures by 0.3 to 4.8 degrees Celsius above present values by the end of this century (Collins et al., 2013). If anthropogenic emissions continue unchecked, the warming increase may reach 8-10 degrees Celsius by 2300 (Rogelj et al., 2012). The contribution that large ice sheets will make to sea-level rise under such warming scenarios is difficult to quantify because the equilibrium-response timescale of ice sheets is longer than those of the atmosphere or ocean. Here we use a coupled ice-sheet/ice-shelf model to show that if atmospheric warming exceeds 1.5 to 2 degrees Celsius above present, collapse of the major Antarctic ice shelves triggers a centennial- to millennial-scale response of the Antarctic ice sheet in which enhanced viscous flow produces a long-term commitment (an unstoppable contribution) to sea-level rise. Our simulations represent the response of the present-day Antarctic ice-sheet system to the oceanic and climatic changes of four representative concentration pathways (RCPs) from the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Collins et al., 2013). We find that substantial Antarctic ice loss can be prevented only by limiting greenhouse gas emissions to RCP 2.6 levels. Higher-emissions scenarios lead to ice loss from Antarctic that will raise sea level by 0.6-3 metres by the year 2300. Our results imply that greenhouse gas emissions in the next few decades will strongly influence the long-term contribution of the Antarctic ice sheet to global sea level.

Emergent cell and tissue dynamics from subcellular modeling of active biomechanical processes

NASA Astrophysics Data System (ADS)

Sandersius, S. A.; Weijer, C. J.; Newman, T. J.

2011-08-01

Cells and the tissues they form are not passive material bodies. Cells change their behavior in response to external biochemical and biomechanical cues. Behavioral changes, such as morphological deformation, proliferation and migration, are striking in many multicellular processes such as morphogenesis, wound healing and cancer progression. Cell-based modeling of these phenomena requires algorithms that can capture active cell behavior and their emergent tissue-level phenotypes. In this paper, we report on extensions of the subcellular element model to model active biomechanical subcellular processes. These processes lead to emergent cell and tissue level phenotypes at larger scales, including (i) adaptive shape deformations in cells responding to slow stretching, (ii) viscous flow of embryonic tissues, and (iii) streaming patterns of chemotactic cells in epithelial-like sheets. In each case, we connect our simulation results to recent experiments.

Nanoconfined ionic liquids: Disentangling electrostatic and viscous forces

NASA Astrophysics Data System (ADS)

Lhermerout, Romain; Perkin, Susan

2018-01-01

Recent reports of surface forces across nanoconfined ionic liquids have revealed the existence of an anomalously long-ranged interaction apparently of electrostatic origin. Ionic liquids are viscous, and therefore it is important to inspect rigorously whether the observed repulsive forces are indeed equilibrium forces or, rather, arise from the viscous force during drainage of the fluid between two confining surfaces. In this paper we present our direct measurements of surface forces between mica sheets approaching in the ionic liquid [C2C1Im ] [NTf2] , exploring three orders of magnitude in approach velocity. Trajectories are systematically fitted by solving the equation of motion, allowing us to disentangle the viscous and equilibrium contributions. First, we find that the drainage obeys classical hydrodynamics with a negative slip boundary condition in the range of the structural force, implying that a nanometer -thick portion of the liquid in the vicinity of the solid surface is composed of ordered molecules that do not contribute to the flow. Second, we show that a long-range static force must indeed be invoked, in addition to the viscous force, in order to describe the data quantitatively. This equilibrium interaction decays exponentially and with decay length in agreement with the screening length reported for the same system in previous studies. In those studies the decay was simply checked to be independent of velocity and measured at a low approach rate, rather than explicitly taking account of viscous effects: we explain why this gives indistinguishable outcomes for the screening length by noting that the viscous force is linear to very good approximation over a wide range of distances.

The fluid dynamics of the chocolate fountain

NASA Astrophysics Data System (ADS)

Townsend, Adam K.; Wilson, Helen J.

2016-01-01

We consider the fluid dynamics of the chocolate fountain. Molten chocolate is a mildly shear-thinning non-Newtonian fluid. Dividing the flow into three main domains—the pumped flow up the centre, the film flow over each dome, and the freely falling curtain flow between the domes—we generate a wide-ranging study of Newtonian and non-Newtonian fluid mechanics. The central pumped flow is a benchmark to elucidate the effects of shear-thinning. The dome flow can be modelled as a thin-film flow with the leading-order effects being a simple balance of gravity and viscosity. Finally, the curtain flow is analytically intractable but is related to the existing theory of water bells (both inviscid and viscous). In pipe flow, Newtonian fluids exhibit a parabolic velocity profile; shear-thinning makes the profile more blunted. In thin-film flow over the dome, gravitational and viscous effects balance and the dome shape is not important beyond the local slope. We find that the chocolate thins and slows down as it travels down the dome. Finally, in the curtain flow, we predict the shape of the falling sheet for an inviscid fluid, and compare this with the literature to predict the shape for a viscous fluid, having shown that viscous forces are too great to ignore. We also find that the primary effect driving the shape of the curtain (which falls inwards towards the axis of the fountain) is surface tension. We find that the three domains provide excellent introductions to non-Newtonian mechanics, the important mathematical technique of scaling, and how to manipulate existing data to make our own predictions. We also find that the topic generates interest among the public in our engagement work.

Postglacial Rebound From Space Geodesy

NASA Astrophysics Data System (ADS)

Argus, D. F.; Peltier, W. R.

2005-12-01

To study the viscous response of the earth to the unloading of the late Pleistocene ice sheets and, to a lesser extent, the elastic response of the earth to current changes in ice sheet mass, we integrate geodetic observations from VLBI over 24 years, from SLR over 23 years, from DORIS over 12 years, and from GPS over 11 years. The excellent geodetic velocity solutions upon which this study are based are from Chopo Ma and Dan MacMillan (Goddard Space Flight Center), Michael Heflin (Jet Propulsion Laboratory), John Ries and Richard Eanes (Center for Space Research, University of Texas at Austin), and Pascal Willis (Institut Geogpraphique National and Jet Propulsion Laboratory). The rates of uplift and subsidence we determine, which in places differ significantly from published studies, are constraining postglacial rebound models like that of Peltier [1994], Peltier [1996], Milne [2001], and Peltier [2004]. We find the following: Yellowknife is rising at 5.7 ±1.8 mm/yr (95% confidence limits), showing [Peltier 2002, 2004], with complimentary ground observations of gravity [Lambert et al. 2001], that the western part of the Laurentide ice sheet was thicker during the Last Glacial Maximum than previously [Peltier 1994, 1996] believed. Onsala (Sweden) is rising at 2.4 ±1.3 mm/yr and Algonquin Park (Ontario) is rising at 2.3 ±1.4 mm/yr, constraining the positions of the margins of the Fennoscandian and Laurentide ice sheets during the Last Glacial Maximum. The eastern United States is falling at ~1 mm/yr, suggesting that the area around the ancient Laurentide ice sheet is subsiding more slowly than predicted by the model of Peltier [2004]. Western and central Europe are falling at ~0.5 mm/yr, suggesting that the area around the ancient Fennoscandian ice sheet is hardly subsiding at all, consistent with the model of Peltier [2004]. Kellyville (Greenland) is falling insignificantly at 1.1 ±4.3 mm/yr, not requiring current loading of the ice sheet [Wahr et al. 2002, Tarasov and Peltier 2002]. Ny Alesund (Spitsbergen) is rising at 6.3 ±1.7 mm/yr, in elastic response to current unloading of a glacier [Hagedoorn and Wolf 2003]. Sites along the margins of the ancient ice sheets are moving horizontally away from the former ice centers at ~1 mm/yr, more slowly than predicted by the model of Peltier [2004]. Sites far from the ancient ice sheets are moving horizontally hardly at all, at insignificant speeds again slower than predicted by the model of Peltier [2004]. Much of the Antarctic, Eurasian, and North American plates are deforming horizontally hardly at all, allowing the angular velocities of the plates to be estimated meaningfully. Satisfying the geodetic observations is going to require revising the ice sheet thickness and the viscosity of the upper mantle (or the elastic lithosphere thickness) in the model of Peltier [2004].

A comprehensive study of the electrically conducting water based CuO and Al2O3 nanoparticles over coupled nanofluid-sheet interface

NASA Astrophysics Data System (ADS)

Ahmad, R.

2016-02-01

Many studies on nanofluid flow over a permeable/impermeable sheet prescribe the kinematics of the sheet and disregard the sheet’s mechanics. However, the current study is one of the infrequent contributions that anticipate the mechanics of both the electrically conducting nanofluid (a homogeneous mixture of nanoparticles and base fluid) and the sheet. Two types of nanoparticles, alumina and copper, with water as a base fluid over the sheet are considered. With the help of the similarity transformations, the corresponding partial differential equations for the coupled nanofluid-sheet interface are transformed into a system of ordinary differential equations. The simulations are done by using the experimentally verified results from the previous studies for viscosity and thermal conductivity. Self-similar solutions are attained by considering both analytical and numerical techniques. Dual skin friction coefficients are attained with different copper and alumina nanoparticles over both the stretching and viscous sheets. The influence of the Eckert number, magnetic and mass suction/blowing parameters on the dimensionless velocity, temperature, skin friction and heat transfer rates over the nanofluid-sheet interface are presented graphically as well as numerically. The obtained results are of potential benefit for studying nanofluid flow over various soft surfaces such as synthetic plastics, soft silicone sheet and soft synthetic rubber sheet. These surfaces are easily deformed by thermal fluctuations.

α-Helix Unwinding as Force Buffer in Spectrins.

PubMed

Takahashi, Hirohide; Rico, Felix; Chipot, Christophe; Scheuring, Simon

2018-03-27

Spectrins are cytoskeletal proteins located at the inner face of the plasma membrane, making connections between membrane anchors and the actin cortex, and between actin filaments. Spectrins share a common structure forming a bundle of 3 α-helices and play a major role during cell deformation. Here, we used high-speed force spectroscopy and steered molecular dynamics simulations to understand the mechanical stability of spectrin, revealing a molecular force buffering function. We find that spectrin acts as a soft spring at short extensions (70-100 Å). Under continuous external stretching, its α-helices unwind, leading to a viscous mechanical response over larger extensions (100-300 Å), represented by a constant-force plateau in force/extension curves. This viscous force buffering emerges from a quasi-equilibrium competition between disruption and re-formation of α-helical hydrogen bonds. Our results suggest that, in contrast to β-sheet proteins, which unfold in a catastrophic event, α-helical spectrins dominantly unwind, providing a viscous force buffer over extensions about 5 times their folded length.

Conditions for a steady ice sheet ice shelf junction

NASA Astrophysics Data System (ADS)

Nowicki, S. M. J.; Wingham, D. J.

2008-01-01

This paper investigates the conditions under which a marine ice sheet may adopt a steady profile. The ice is treated as a linear viscous fluid caused to flow from a rigid base to and over water, treated as a denser but inviscid fluid. The solutions in the region around the point of flotation, or 'transition' zone, are calculated numerically. In-flow and out-flow conditions appropriate to ice sheet and ice shelf flow are applied at the ends of the transition zone and the rigid base is specified; the flow and steady free surfaces are determined as part of the solutions. The basal stress upstream, and the basal deflection downstream, of the flotation point are examined to determine which of these steady solutions satisfy 'contact' conditions that would prevent (i) the steady downstream basal deflection contacting the downstream base, and (ii) the upstream ice commencing to float in the event it was melted at the base. In the case that the upstream bed is allowed to slide, we find only one mass flux that satisfies the contact conditions. When no sliding is allowed at the bed, however, we find a range of mass fluxes satisfy the contact conditions. The effect of 'backpressure' on the solutions is investigated, and is found to have no affect on the qualitative behaviour of the junctions. To the extent that the numerical, linearly viscous treatment may be applied to the case of ice flowing out over the ocean, we conclude that when sliding is present, Weertman's 'instability' hypothesis holds.

Toothpaste lava from the Barren Island volcano (Andaman Sea)

NASA Astrophysics Data System (ADS)

Sheth, Hetu C.; Ray, Jyotiranjan S.; Kumar, Alok; Bhutani, Rajneesh; Awasthi, Neeraj

2011-04-01

Toothpaste lava is a basaltic lava flow type transitional between pahoehoe and aa and has been described from Paricutin, Kilauea and Etna volcanoes. Here we describe a spectacular example of toothpaste lava, forming part of a recent (possibly 1994-95) aa flow on the active volcano of Barren Island (Andaman Sea). This flow of subalkalic basalt shows abundant squeeze-ups of viscous toothpasate lava near its entry into the sea. The squeeze-ups are sheets and slabs, up to several meters across and tens of centimeters thick, extruded from boccas. They are often prominently curved, have striated upper surfaces with close-spaced, en echelon linear ridges and grooves, broad wave-like undulations perpendicular to the striations, and sometimes, clefts. Textural, geochemical, and Sr-Nd isotopic data on the squeeze-ups and the exposed aa flow core indicate very crystal-rich, viscous, and isotopically very homogeneous lava. We envisage that a greatly reduced speed of this viscous flow at the coastline, possibly aided by a shallowing of the basal slope, led to lateral spreading of the flow, which caused tension in its upper parts. This, with continued (albeit dwindling) lava supply at the back, led to widespread tearing of the flow surface and extrusion of the squeeze-ups. The larger slabs, while extruding in a plastic condition, curved under their own weight, whereas their surfaces experienced brittle deformation, forming the en echelon grooves. The extruded, detached, and rotated sheets and slabs were carried forward for some distance atop the very slowly advancing aa core, before the flow solidified.

Entropy generation in a second grade magnetohydrodynamic nanofluid flow over a convectively heated stretching sheet with nonlinear thermal radiation and viscous dissipation

NASA Astrophysics Data System (ADS)

Sithole, Hloniphile; Mondal, Hiranmoy; Sibanda, Precious

2018-06-01

This study addresses entropy generation in magnetohydrodynamic flow of a second grade nanofluid over a convectively heated stretching sheet with nonlinear thermal radiation and viscous dissipation. The second grade fluid is assumed to be electrically conducting and is permeated by an applied non-uniform magnetic field. We further consider the impact on the fluid properties and the Nusselt number of homogeneous-heterogeneous reactions and a convective boundary condition. The mathematical equations are solved using the spectral local linearization method. Computations for skin-friction coefficient and local Nusselt number are carried out and displayed in a table. It is observed that the effects of the thermophoresis parameter is to increase the temperature distributions throughout the boundary layer. The entropy generation is enhanced by larger magnetic parameters and increasing Reynolds number. The aim of this manuscript is to pay more attention of entropy generation analysis with heat and fluid flow on second grade nanofluids to improve the system performance. Also the fluid velocity and temperature in the boundary layer region rise significantly for increasing the values of the second grade nanofluid parameter.

Geometry of thin liquid sheet flows

NASA Technical Reports Server (NTRS)

Chubb, Donald L.; Calfo, Frederick D.; Mcconley, Marc W.; Mcmaster, Matthew S.; Afjeh, Abdollah A.

1994-01-01

Incompresible, thin sheet flows have been of research interest for many years. Those studies were mainly concerned with the stability of the flow in a surrounding gas. Squire was the first to carry out a linear, invicid stability analysis of sheet flow in air and compare the results with experiment. Dombrowski and Fraser did an experimental study of the disintegration of sheet flows using several viscous liquids. They also detected the formulation of holes in their sheet flows. Hagerty and Shea carried out an inviscid stability analysis and calculated growth rates with experimental values. They compared their calculated growth rates with experimental values. Taylor studied extensively the stability of thin liquid sheets both theoretically and experimentally. He showed that thin sheets in a vacuum are stable. Brown experimentally investigated thin liquid sheet flows as a method of application of thin films. Clark and Dumbrowski carried out second-order stability analysis for invicid sheet flows. Lin introduced viscosity into the linear stability analysis of thin sheet flows in a vacuum. Mansour and Chigier conducted an experimental study of the breakup of a sheet flow surrounded by high-speed air. Lin et al. did a linear stability analysis that included viscosity and a surrounding gas. Rangel and Sirignano carried out both a linear and nonlinear invisid stability analysis that applies for any density ratio between the sheet liquid and the surrounding gas. Now there is renewed interest in sheet flows because of their possible application as low mass radiating surfaces. The objective of this study is to investigate the fluid dynamics of sheet flows that are of interest for a space radiator system. Analytical expressions that govern the sheet geometry are compared with experimental results. Since a space radiator will operate in a vacuum, the analysis does not include any drag force on the sheet flow.

Magnetohydrodynamic dissipative flow across the slendering stretching sheet with temperature dependent variable viscosity

NASA Astrophysics Data System (ADS)

Jayachandra Babu, M.; Sandeep, N.; Ali, M. E.; Nuhait, Abdullah O.

The boundary layer flow across a slendering stretching sheet has gotten awesome consideration due to its inexhaustible pragmatic applications in nuclear reactor technology, acoustical components, chemical and manufacturing procedures, for example, polymer extrusion, and machine design. By keeping this in view, we analyzed the two-dimensional MHD flow across a slendering stretching sheet within the sight of variable viscosity and viscous dissipation. The sheet is thought to be convectively warmed. Convective boundary conditions through heat and mass are employed. Similarity transformations used to change over the administering nonlinear partial differential equations as a group of nonlinear ordinary differential equations. Runge-Kutta based shooting technique is utilized to solve the converted equations. Numerical estimations of the physical parameters involved in the problem are calculated for the friction factor, local Nusselt and Sherwood numbers. Viscosity variation parameter and chemical reaction parameter shows the opposite impact to each other on the concentration profile. Heat and mass transfer Biot numbers are helpful to enhance the temperature and concentration respectively.

A k-Omega Turbulence Model for Quasi-Three-Dimensional Turbomachinery Flows

NASA Technical Reports Server (NTRS)

Chima, Rodrick V.

1995-01-01

A two-equation k-omega turbulence model has been developed and applied to a quasi-three-dimensional viscous analysis code for blade-to-blade flows in turbomachinery. the code includes the effects of rotation, radius change, and variable stream sheet thickness. The flow equations are given and the explicit runge-Kutta solution scheme is described. the k-omega model equations are also given and the upwind implicit approximate-factorization solution scheme is described. Three cases were calculated: transitional flow over a flat plate, a transonic compressor rotor, and transonic turbine vane with heat transfer. Results were compared to theory, experimental data, and to results using the Baldwin-Lomax turbulence model. The two models compared reasonably well with the data and surprisingly well with each other. Although the k-omega model behaves well numerically and simulates effects of transition, freestream turbulence, and wall roughness, it was not decisively better than the Baldwin-Lomax model for the cases considered here.

Peeling without precursors

NASA Astrophysics Data System (ADS)

Lister, John; Skinner, Dominic; Large, Tim

2017-11-01

The peeling by fluid injection of an elastic sheet away from a substrate is often regularised by invoking a thin prewetting film or a low-viscosity phase in the tip. Here we analyse fluid-driven peeling without such precursors, and consider an elastic sheet either bonded to, or simply laid on, an elastic substrate. To resolve the `elastic contact-line problem' that arises from viscous flow and beam theory, we determine the near-tip behaviour from lubrication theory coupled to the full equations of elasticity and fracture. The result is a law for the tip propagation speed in terms of the remote loading and the toughness of the sheet-substrate bonding (which might be zero). There are distinct modes of failure, according to whether there is slip ahead of the fluid front. The propagation-speed law gives rise to new similarity solutions for the spread of a fluid-filled blister in different regimes.

The stress field below the NE German Basin: effects induced by the Alpine collision

NASA Astrophysics Data System (ADS)

Marotta, A. M.; Bayer, U.; Scheck, M.; Thybo, H.

2001-02-01

We use a thin-sheet approach for a viscous lithosphere to investigate the effects induced by the Alpine collision on the vertical deformation and regional stress in northern Europe, focusing on the NE German Basin. New seismic studies indicate a flexural-type deep crustal structure under the basin, which may be induced by compressive forces transmitted from the south and related to Alpine tectonics. Finite element techniques are used to solve the vertical deformation and stress field for a viscous European lithosphere with horizontal rheological heterogeneities. Our results support the idea that a relatively strong lithosphere below the northern margin of the German Basin at the transition into the Baltic Shield may explain the characteristic regional stress field, especially the fan-like pattern that is observed within the region.

Impact of Beads and Drops on a Repellent Solid Surface: A Unified Description

NASA Astrophysics Data System (ADS)

Arora, S.; Fromental, J.-M.; Mora, S.; Phou, Ty; Ramos, L.; Ligoure, C.

2018-04-01

We investigate freely expanding sheets formed by ultrasoft gel beads, and liquid and viscoelastic drops, produced by the impact of the bead or drop on a silicon wafer covered with a thin layer of liquid nitrogen that suppresses viscous dissipation thanks to an inverse Leidenfrost effect. Our experiments show a unified behavior for the impact dynamics that holds for solids, liquids, and viscoelastic fluids and that we rationalize by properly taking into account elastocapillary effects. In this framework, the classical impact dynamics of solids and liquids, as far as viscous dissipation is negligible, appears as the asymptotic limits of a universal theoretical description. A novel material-dependent characteristic velocity that includes both capillary and bulk elasticity emerges from this unified description of the physics of impact.

An Exploration of Heating Mechanisms in a Supra-arcade Plasma Sheet Formed after a Coronal Mass Ejection

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

Reeves, Katharine K.; Freed, Michael S.; McKenzie, David E.

We perform a detailed analysis of the thermal structure of the region above the post-eruption arcade for a flare that occurred on 2011 October 22. During this event, a sheet of hot plasma is visible above the flare loops in the 131 Å bandpass of the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory . Supra-arcade downflows (SADs) are observed traveling sunward through the post-eruption plasma sheet. We calculate differential emission measures using the AIA data and derive an emission measure weighted average temperature in the supra-arcade region. In areas where many SADs occur, the temperature of the supra-arcademore » plasma tends to increase, while in areas where no SADs are observed, the temperature tends to decrease. We calculate the plane-of-sky velocities in the supra-arcade plasma and use them to determine the potential heating due to adiabatic compression and viscous heating. Of the 13 SADs studied, 10 have noticeable signatures in both the adiabatic and the viscous terms. The adiabatic heating due to compression of plasma in front of the SADs is on the order of 0.1–0.2 MK/s, which is similar in magnitude to the estimated conductive cooling rate. This result supports the notion that SADs contribute locally to the heating of plasma in the supra-arcade region. We also find that in the region without SADs, the plasma cools at a rate that is slower than the estimated conductive cooling, indicating that additional heating mechanisms may act globally to keep the plasma temperature high.« less

Large-scale deformation related to the collision of the Aleutian Arc with Kamchatka

USGS Publications Warehouse

Gesit, Eric L.; Scholl, David W.

1994-01-01

The far western Aleutian Island Arc is actively colliding with Kamchatka. Westward motion of the Aleutian Arc is brought about by the tangential relative motion of the Pacific plate transferred to major, right-lateral shear zones north and south of the arc. Early geologic mapping of Cape Kamchatka (a promontory of Kamchatka along strike with the Aleutian Arc) revealed many similarities to the geology of the Aleutian Islands. Later studies support the notion that Cape Kamchatka is the farthest west Aleutian “island” and that it has been accreted to Kamchatka by the process of arc-continent collision. Deformation associated with the collision onshore Kamchatka includes gravimetrically determined crustal thickening and formation of a narrow thrust belt of intensely deformed rocks directly west of Cape Kamchatka. The trend of the thrust faults is concave toward the collision zone, indicating a radial distribution of maximum horizontal compressive stress. Offshore, major crustal faults trend either oblique to the Kamchatka margin or parallel to major Aleutian shear zones. These offshore faults are complex, accommodating both strike-slip and thrust displacements as documented by focal mechanisms and seismic reflection data. Earthquake activity is much higher in the offshore region within a zone bounded to the north by the northernmost Aleutian shear zone and to the west by an apparent aseismic front. Analysis of focal mechanisms in the region indicate that the present-day arc-continent “contact zone” is located directly east of Cape Kamchatka. In modeling the dynamics of the collision zone using thin viscous sheet theory, the rheological parameters are only partially constrained to values of n (the effective power law exponent) ≥ 3 and Ar(the Argand number) ≤ 30. These values are consistent with a forearc thermal profile of Kamchatka, previously determined from heat flow modeling. The thin viscous sheet modeling also indicates that onshore thrust faulting is a consequence, not only of compressive stresses resulting from the west directed collision, but also of sediment-induced coupling of the subducting Pacific plate.

Three-dimensional morphogenesis of MDCK cells induced by cellular contractile forces on a viscous substrate

PubMed Central

Imai, Misako; Furusawa, Kazuya; Mizutani, Takeomi; Kawabata, Kazushige; Haga, Hisashi

2015-01-01

Substrate physical properties are essential for many physiological events such as embryonic development and 3D tissue formation. Physical properties of the extracellular matrix such as viscoelasticity and geometrical constraints are understood as factors that affect cell behaviour. In this study, we focused on the relationship between epithelial cell 3D morphogenesis and the substrate viscosity. We observed that Madin-Darby Canine Kidney (MDCK) cells formed 3D structures on a viscous substrate (Matrigel). The structures appear as a tulip hat. We then changed the substrate viscosity by genipin (GP) treatment. GP is a cross-linker of amino groups. Cells cultured on GP-treated-matrigel changed their 3D morphology in a substrate viscosity-dependent manner. Furthermore, to elucidate the spatial distribution of the cellular contractile force, localization of mono-phosphorylated and di-phosphorylated myosin regulatory light chain (P-MRLCs) was visualized by immunofluorescence. P-MRLCs localized along the periphery of epithelial sheets. Treatment with Y-27632, a Rho-kinase inhibitor, blocked the P-MRLCs localization at the edge of epithelial sheets and halted 3D morphogenesis. Our results indicate that the substrate viscosity, the substrate deformation, and the cellular contractile forces induced by P-MRLCs play crucial roles in 3D morphogenesis. PMID:26374384

Three-dimensional morphogenesis of MDCK cells induced by cellular contractile forces on a viscous substrate.

PubMed

Imai, Misako; Furusawa, Kazuya; Mizutani, Takeomi; Kawabata, Kazushige; Haga, Hisashi

2015-09-16

Substrate physical properties are essential for many physiological events such as embryonic development and 3D tissue formation. Physical properties of the extracellular matrix such as viscoelasticity and geometrical constraints are understood as factors that affect cell behaviour. In this study, we focused on the relationship between epithelial cell 3D morphogenesis and the substrate viscosity. We observed that Madin-Darby Canine Kidney (MDCK) cells formed 3D structures on a viscous substrate (Matrigel). The structures appear as a tulip hat. We then changed the substrate viscosity by genipin (GP) treatment. GP is a cross-linker of amino groups. Cells cultured on GP-treated-matrigel changed their 3D morphology in a substrate viscosity-dependent manner. Furthermore, to elucidate the spatial distribution of the cellular contractile force, localization of mono-phosphorylated and di-phosphorylated myosin regulatory light chain (P-MRLCs) was visualized by immunofluorescence. P-MRLCs localized along the periphery of epithelial sheets. Treatment with Y-27632, a Rho-kinase inhibitor, blocked the P-MRLCs localization at the edge of epithelial sheets and halted 3D morphogenesis. Our results indicate that the substrate viscosity, the substrate deformation, and the cellular contractile forces induced by P-MRLCs play crucial roles in 3D morphogenesis.

One-dimensional reduction of viscous jets. II. Applications

NASA Astrophysics Data System (ADS)

Pitrou, Cyril

2018-04-01

In a companion paper [Phys. Rev. E 97, 043115 (2018), 10.1103/PhysRevE.97.043115], a formalism allowing to describe viscous fibers as one-dimensional objects was developed. We apply it to the special case of a viscous fluid torus. This allows to highlight the differences with the basic viscous string model and with its viscous rod model extension. In particular, an elliptic deformation of the torus section appears because of surface tension effects, and this cannot be described by viscous string nor viscous rod models. Furthermore, we study the Rayleigh-Plateau instability for periodic deformations around the perfect torus, and we show that the instability is not sufficient to lead to the torus breakup in several droplets before it collapses to a single spherical drop. Conversely, a rotating torus is dynamically attracted toward a stationary solution, around which the instability can develop freely and split the torus in multiple droplets.

Sources of Local Time Asymmetries in Magnetodiscs

NASA Astrophysics Data System (ADS)

Arridge, C. S.; Kane, M.; Sergis, N.; Khurana, K. K.; Jackman, C. M.

2015-04-01

The rapidly rotating magnetospheres at Jupiter and Saturn contain a near-equatorial thin current sheet over most local times known as the magnetodisc, resembling a wrapped-up magnetotail. The Pioneer, Voyager, Ulysses, Galileo, Cassini and New Horizons spacecraft at Jupiter and Saturn have provided extensive datasets from which to observationally identify local time asymmetries in these magnetodiscs. Imaging in the infrared and ultraviolet from ground- and space-based instruments have also revealed the presence of local time asymmetries in the aurora which therefore must map to local time asymmetries in the magnetosphere. Asymmetries are found in (i) the configuration of the magnetic field and magnetospheric currents, where a thicker disc is found in the noon and dusk sectors; (ii) plasma flows where the plasma flow has local time-dependent radial components; (iii) a thicker plasma sheet in the dusk sector. Many of these features are also reproduced in global MHD simulations. Several models have been developed to interpret these various observations and typically fall into two groups: ones which invoke coupling with the solar wind (via reconnection or viscous processes) and ones which invoke internal rotational processes operating inside an asymmetrical external boundary. In this paper we review these observational in situ findings, review the models which seek to explain them, and highlight open questions and directions for future work.

Viscous cosmology in new holographic dark energy model and the cosmic acceleration

NASA Astrophysics Data System (ADS)

Singh, C. P.; Srivastava, Milan

2018-03-01

In this work, we study a flat Friedmann-Robertson-Walker universe filled with dark matter and viscous new holographic dark energy. We present four possible solutions of the model depending on the choice of the viscous term. We obtain the evolution of the cosmological quantities such as scale factor, deceleration parameter and transition redshift to observe the effect of viscosity in the evolution. We also emphasis upon the two independent geometrical diagnostics for our model, namely the statefinder and the Om diagnostics. In the first case we study new holographic dark energy model without viscous and obtain power-law expansion of the universe which gives constant deceleration parameter and statefinder parameters. In the limit of the parameter, the model approaches to Λ CDM model. In new holographic dark energy model with viscous, the bulk viscous coefficient is assumed as ζ =ζ 0+ζ 1H, where ζ 0 and ζ 1 are constants, and H is the Hubble parameter. In this model, we obtain all possible solutions with viscous term and analyze the expansion history of the universe. We draw the evolution graphs of the scale factor and deceleration parameter. It is observed that the universe transits from deceleration to acceleration for small values of ζ in late time. However, it accelerates very fast from the beginning for large values of ζ . By illustrating the evolutionary trajectories in r-s and r-q planes, we find that our model behaves as an quintessence like for small values of viscous coefficient and a Chaplygin gas like for large values of bulk viscous coefficient at early stage. However, model has close resemblance to that of the Λ CDM cosmology in late time. The Om has positive and negative curvatures for phantom and quintessence models, respectively depending on ζ . Our study shows that the bulk viscosity plays very important role in the expansion history of the universe.

A mechanism for crustal recycling on Venus

NASA Technical Reports Server (NTRS)

Lenardic, A.; Kaula, W. M.; Bindschadler, D. L.

1993-01-01

Entrainment of lower crust by convective mantle downflows is proposed as a crustal recycling mechanism on Venus. The mechanism is characterized by thin sheets of crust being pulled into the mantle by viscous flow stresses. Finite element models of crust/mantle interaction are used to explore tectonic conditions under which crustal entrainment may occur. The recycling scenarios suggested by the numerical models are analogous to previously studied problems for which analytic and experimental relationships assessing entrainment rates have been derived. We use these relationships to estimate crustal recycling rates on Venus. Estimated rates are largely determined by (1) strain rate at the crust/mantle interface (higher strain rate leads to greater entrainment); and (2) effective viscosity of the lower crust (viscosity closer to that of mantle lithosphere leads to greater entrainment). Reasonable geologic strain rates and available crustal flow laws suggest entrainment can recycle approximately equal 1 cu km of crust per year under favorable conditions.

An experimental study of non-isothermal miscible displacements in a Hele-Shaw cell

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

Nagatsu, Yuichiro; Fujita, Norihito; Kato, Yoshihito

Non-isothermal miscible displacements in a radial Hele-Shaw cell were experimentally investigated using a scheme in which room temperature liquids of relatively high viscosity were displaced by high-temperature (80 C), less-viscous liquids. Fundamental characteristics have been presented regarding how the effect of a non-isothermal field on miscible displacement patterns varies in terms of factors such as the viscosity ratio of the more- and less-viscous liquids at 20 C, M{sub 20}, the rate of an increase in the pattern's area, R, and the gap width of the cell, b. The concept of area density was used to quantitatively evaluate the effect ofmore » the non-isothermal fields on the patterns. We have found that the effect of the non-isothermal field on the patterns does not monotonically vary with M{sub 20} and b. In contrast, it increases with R in the present experimental condition. The experimental results can be explained by introducing an assumption in which heat is transferred mainly to the plates of the cell, in other words, the temperature of the more-viscous liquid remains constant, whereas that of the less-viscous liquid spatiotemporally decreases and the viscosity of it increases along with the temperature decrease. Visualization of non-isothermal field in the cell has been done by means of a thermo sheet and the results support the assumption mentioned above. (author)« less

The evolution of impact basins - Viscous relaxation of topographic relief. [for lunar surface modeling

NASA Technical Reports Server (NTRS)

Solomon, S. C.; Comer, R. P.; Head, J. W.

1982-01-01

A topographic profile of the young large lunar basin, Orientale, is presented in order to examine the effects of viscous relaxation on basin topography. Analytical models for viscous flow are considered, showing a wavelength-dependence of time constants for viscous decay on the decrease in viscosity with depth and on the extent of the isostatic compensation of the initial topography. Lunar rheological models which are developed include a half-space model for uniform Newtonian viscosity, density, and gravitational acceleration, a layer over inviscid half space model with material inviscid over geological time scales, and a layer with isostatic compensation where a uniformly viscous layer overlies an inviscid half space of higher density. Greater roughness is concluded, and has been observed, on the moon's dark side due to continued lower temperatures since the time of heavy bombardment.

New thermodynamics of entropy generation minimization with nonlinear thermal radiation and nanomaterials

NASA Astrophysics Data System (ADS)

Hayat, T.; Khan, M. Ijaz; Qayyum, Sumaira; Alsaedi, A.; Khan, M. Imran

2018-03-01

This research addressed entropy generation for MHD stagnation point flow of viscous nanofluid over a stretching surface. Characteristics of heat transport are analyzed through nonlinear radiation and heat generation/absorption. Nanoliquid features for Brownian moment and thermophoresis have been considered. Fluid in the presence of constant applied inclined magnetic field is considered. Flow problem is mathematically modeled and governing expressions are changed into nonlinear ordinary ones by utilizing appropriate transformations. The effects of pertinent variables on velocity, nanoparticle concentration and temperature are discussed graphically. Furthermore Brownian motion and thermophoresis effects on entropy generation and Bejan number have been examined. Total entropy generation is inspected through various flow variables. Consideration is mainly given to the convergence process. Velocity, temperature and mass gradients at the surface of sheet are calculated numerically.

On ice rifts and the stability of non-Newtonian extensional flows on a sphere

NASA Astrophysics Data System (ADS)

Sayag, Roiy

2017-11-01

Rifts that form at the fronts of floating ice shelves that spread into the ocean can trigger major calving events in the ice. The deformation of ice can be modeled as a thin viscous film driven by buoyancy. The front of such a viscous film that propagates over a flat surface with no-slip basal conditions is known to have stable axisymmetric solutions. In contrast, when the fluid propagates under free-slip conditions at the substrate, the front can become unstable to small perturbations if the fluid is sufficiently strain-rate softening. Consequently, the front will develop tongues with a characteristic wavelength that coarsens over time, a pattern that is reminiscent of ice rifts. Here we investigate the stability of a spherical sheet of power-law fluids under free-slip basal conditions. The fluid is discharged at constant flux and axisymmetrically with respect to the pole, and propagates towards the equator. The propagating front in such a situation may become unstable due to its failure to sustain large extensional forces, resulting in the formation of rifts. This study has implications to understanding the cause of patterns that are observed on shells of floating ice in a range of planetary objects, and whether open rifts that sustain life were feasible in snowball earth. Israel Science Foundation 1368/16.

Rotorcraft acoustic radiation prediction based on a refined blade-vortex interaction model

NASA Astrophysics Data System (ADS)

Rule, John Allen

1997-08-01

The analysis of rotorcraft aerodynamics and acoustics is a challenging problem, primarily due to the fact that a rotorcraft continually flies through its own wake. The generation mechanism for a rotorcraft wake, which is dominated by strong, concentrated blade-tip trailing vortices, is similar to that in fixed wing aerodynamics. However, following blades encounter shed vortices from previous blades before they are swept downstream, resulting in sharp, impulsive loading on the blades. The blade/wake encounter, known as Blade-Vortex Interaction, or BVI, is responsible for a significant amount of vibratory loading and the characteristic rotorcraft acoustic signature in certain flight regimes. The present work addressed three different aspects of this interaction at a fundamental level. First, an analytical model for the prediction of trailing vortex structure is discussed. The model as presented is the culmination of a lengthy research effort to isolate the key physical mechanisms which govern vortex sheet rollup. Based on the Betz model, properties of the flow such as mass flux, axial momentum flux, and axial flux of angular momentum are conserved on either a differential or integral basis during the rollup process. The formation of a viscous central core was facilitated by the assumption of a turbulent mixing process with final vortex velocity profiles chosen to be consistent with a rotational flow mixing model and experimental observation. A general derivation of the method is outlined, followed by a comparison of model predictions with experimental vortex measurements, and finally a viscous blade drag model to account for additional effects of aerodynamic drag on vortex structure. The second phase of this program involved the development of a new formulation of lifting surface theory with the ultimate goal of an accurate, reduced order hybrid analytical/numerical model for fast rotorcraft load calculations. Currently, accurate rotorcraft airload analyses are limited by the massive computational power required to capture the small time scale events associated with BVI. This problem has two primary facets: accurate knowledge of the wake geometry, and accurate resolution of the impulsive loading imposed by a tip vortex on a blade. The present work addressed the second facet, providing a mathematical framework for solving the impulsive loading problem analytically, then asymptotically matching this solution to a low-resolution numerical calculation. A method was developed which uses continuous sheets of integrated boundary elements to model the lifting surface and wake. Special elements were developed to capture local behavior in high-gradient regions of the flow, thereby reducing the burden placed on the surrounding numerical method. Unsteady calculations for several classical cases were made in both frequency and time domain to demonstrate the performance of the method. Finally, a new unsteady, compressible boundary element method was applied to the problem of BVI acoustic radiation prediction. This numerical method, combined with the viscous core trailing vortex model, was used to duplicate the geometry and flight configuration of a detailed experimental BVI study carried out at NASA Ames Research Center. Blade surface pressure and near- and far-field acoustic radiation calculations were made. All calculations were shown to compare favorably with experimentally measured values. The linear boundary element method with non-linear corrections proved sufficient over most of the rotor azimuth, and particular in the region of the blade vortex interaction, suggesting that full non-linear CFD schemes are not necessary for rotorcraft noise prediction.

Self-propulsion of a planar electric or magnetic microbot immersed in a polar viscous fluid

NASA Astrophysics Data System (ADS)

Felderhof, B. U.

2011-05-01

A planar sheet immersed in an electrically polar liquid like water can propel itself by means of a plane wave charge density propagating in the sheet. The corresponding running electric wave polarizes the fluid and causes an electrical torque density to act on the fluid. The sheet is convected by the fluid motion resulting from the conversion of rotational particle motion, generated by the torque density, into translational fluid motion by the mechanism of friction and spin diffusion. Similarly, a planar sheet immersed in a magnetic ferrofluid can propel itself by means of a plane wave current density in the sheet and the torque density acting on the fluid corresponding to the running wave magnetic field and magnetization. The effect is studied on the basis of the micropolar fluid equations of motion and Maxwell’s equations of electrostatics or magnetostatics, respectively. An analytic expression is derived for the velocity of the sheet by perturbation theory to second order in powers of the amplitude of the driving charge or current density. Under the assumption that the equilibrium magnetic equation of state may be used in linearized form and that higher harmonics than the first may be neglected, a set of self-consistent integral equations is derived which can be solved numerically by iteration. In typical situations the second-order perturbation theory turns out to be quite accurate.

ISSM-SESAW v1.0: mesh-based computation of gravitationally consistent sea-level and geodetic signatures caused by cryosphere and climate driven mass change

NASA Astrophysics Data System (ADS)

Adhikari, Surendra; Ivins, Erik R.; Larour, Eric

2016-03-01

A classical Green's function approach for computing gravitationally consistent sea-level variations associated with mass redistribution on the earth's surface employed in contemporary sea-level models naturally suits the spectral methods for numerical evaluation. The capability of these methods to resolve high wave number features such as small glaciers is limited by the need for large numbers of pixels and high-degree (associated Legendre) series truncation. Incorporating a spectral model into (components of) earth system models that generally operate on a mesh system also requires repetitive forward and inverse transforms. In order to overcome these limitations, we present a method that functions efficiently on an unstructured mesh, thus capturing the physics operating at kilometer scale yet capable of simulating geophysical observables that are inherently of global scale with minimal computational cost. The goal of the current version of this model is to provide high-resolution solid-earth, gravitational, sea-level and rotational responses for earth system models operating in the domain of the earth's outer fluid envelope on timescales less than about 1 century when viscous effects can largely be ignored over most of the globe. The model has numerous important geophysical applications. For example, we compute time-varying computations of global geodetic and sea-level signatures associated with recent ice-sheet changes that are derived from space gravimetry observations. We also demonstrate the capability of our model to simultaneously resolve kilometer-scale sources of the earth's time-varying surface mass transport, derived from high-resolution modeling of polar ice sheets, and predict the corresponding local and global geodetic signatures.

Mesh-free distributed point source method for modeling viscous fluid motion between disks vibrating at ultrasonic frequency.

PubMed

Wada, Yuji; Kundu, Tribikram; Nakamura, Kentaro

2014-08-01

The distributed point source method (DPSM) is extended to model wave propagation in viscous fluids. Appropriate estimation on attenuation and boundary layer formation due to fluid viscosity is necessary for the ultrasonic devices used for acoustic streaming or ultrasonic levitation. The equations for DPSM modeling in viscous fluids are derived in this paper by decomposing the linearized viscous fluid equations into two components-dilatational and rotational components. By considering complex P- and S-wave numbers, the acoustic fields in viscous fluids can be calculated following similar calculation steps that are used for wave propagation modeling in solids. From the calculations reported the precision of DPSM is found comparable to that of the finite element method (FEM) for a fundamental ultrasonic field problem. The particle velocity parallel to the two bounding surfaces of the viscous fluid layer between two rigid plates (one in motion and one stationary) is calculated. The finite element results agree well with the DPSM results that were generated faster than the transient FEM results.

Continuation of tailored composite structures of ordered staple thermoplastic material

NASA Technical Reports Server (NTRS)

Santare, Michael H.; Pipes, R. Byron

1992-01-01

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Kaluza-Klein Bulk Viscous Fluid Cosmological Models and the Validity of the Second Law of Thermodynamics in f(R, T) Gravity

NASA Astrophysics Data System (ADS)

Samanta, Gauranga Charan; Myrzakulov, Ratbay; Shah, Parth

2017-04-01

The authors considered the bulk viscous fluid in f(R, T) gravity within the framework of Kaluza-Klein space time. The bulk viscous coefficient (ξ) expressed as ξ = {ξ_0} + {ξ_1}{{\\dot a} \\over a} + {ξ_2}{{\\ddot a} \\over {\\dot a}}, where ξ0, ξ1, and ξ2 are positive constants. We take p=(γ-1)ρ, where 0≤γ≤2 as an equation of state for perfect fluid. The exact solutions to the corresponding field equations are given by assuming a particular model of the form of f(R, T)=R+2f(T), where f(T)=λT, λ is constant. We studied the cosmological model in two stages, in first stage: we studied the model with no viscosity, and in second stage: we studied the model involve with viscosity. The cosmological model involve with viscosity is studied by five possible scenarios for bulk viscous fluid coefficient (ξ). The total bulk viscous coefficient seems to be negative, when the bulk viscous coefficient is proportional to {ξ _2}{{\\ddot a} \\over {\\dot a}}, hence, the second law of thermodynamics is not valid; however, it is valid with the generalised second law of thermodynamics. The total bulk viscous coefficient seems to be positive, when the bulk viscous coefficient is proportional to ξ = {ξ _1}{{\\dot a} \\over a} + {ξ _2}{{\\ddot a} \\over {\\dot a}} and ξ = {ξ _0} + {ξ _1}{{\\dot a} \\over a} + {ξ _2}{{\\ddot a} \\over {\\dot a}}, so the second law of thermodynamics and the generalised second law of thermodynamics is satisfied throughout the evolution. We calculate statefinder parameters of the model and observed that it is different from the ∧CDM model. Finally, some physical and geometrical properties of the models are discussed.

An analysis of the viscous flow through a compact radial turbine by the average passage approach

NASA Technical Reports Server (NTRS)

Heidmann, James D.; Beach, Timothy A.

1990-01-01

A steady, three-dimensional viscous average passage computer code is used to analyze the flow through a compact radial turbine rotor. The code models the flow as spatially periodic from blade passage to blade passage. Results from the code using varying computational models are compared with each other and with experimental data. These results include blade surface velocities and pressures, exit vorticity and entropy contour plots, shroud pressures, and spanwise exit total temperature, total pressure, and swirl distributions. The three computational models used are inviscid, viscous with no blade clearance, and viscous with blade clearance. It is found that modeling viscous effects improves correlation with experimental data, while modeling hub and tip clearances further improves some comparisons. Experimental results such as a local maximum of exit swirl, reduced exit total pressures at the walls, and exit total temperature magnitudes are explained by interpretation of the flow physics and computed secondary flows. Trends in the computed blade loading diagrams are similarly explained.

Diverse Geological Applications For Basil: A 2d Finite-deformation Computational Algorithm

NASA Astrophysics Data System (ADS)

Houseman, Gregory A.; Barr, Terence D.; Evans, Lynn

Geological processes are often characterised by large finite-deformation continuum strains, on the order of 100% or greater. Microstructural processes cause deformation that may be represented by a viscous constitutive mechanism, with viscosity that may depend on temperature, pressure, or strain-rate. We have developed an effective com- putational algorithm for the evaluation of 2D deformation fields produced by Newto- nian or non-Newtonian viscous flow. With the implementation of this algorithm as a computer program, Basil, we have applied it to a range of diverse applications in Earth Sciences. Viscous flow fields in 2D may be defined for the thin-sheet case or, using a velocity-pressure formulation, for the plane-strain case. Flow fields are represented using 2D triangular elements with quadratic interpolation for velocity components and linear for pressure. The main matrix equation is solved by an efficient and compact conjugate gradient algorithm with iteration for non-Newtonian viscosity. Regular grids may be used, or grids based on a random distribution of points. Definition of the prob- lem requires that velocities, tractions, or some combination of the two, are specified on all external boundary nodes. Compliant boundaries may also be defined, based on the idea that traction is opposed to and proportional to boundary displacement rate. In- ternal boundary segments, allowing fault-like displacements within a viscous medium have also been developed, and we find that the computed displacement field around the fault tip is accurately represented for Newtonian and non-Newtonian viscosities, in spite of the stress singularity at the fault tip. Basil has been applied by us and colleagues to problems that include: thin sheet calculations of continental collision, Rayleigh-Taylor instability of the continental mantle lithosphere, deformation fields around fault terminations at the outcrop scale, stress and deformation fields in and around porphyroblasts, and deformation of the subducted oceanic slab. Application of Basil to a diverse range of topics is facilitated by the use of command syntax input files that allow most aspects of the calculation to be controlled easily, together with a post-processing package, Sybil, for display and interpretation of the results. Sybil uses a menu-driven graphical interface to access a powerful combination of commands, to- gether with log files that allow repetitive tasks to be more automated

Symposium of Naval Hydrodynamics (14th) held at Ann Arbor, Michigan on August 23-27, 1982,

DTIC Science & Technology

1982-01-01

Chahine -Viscous Effects on the Stability of Cavitating Line Vortices -. 195 Jaakko V. Pylkknen Nuclei and Cavitation 215 Jean -Pierre Le G9fu and Yves...the sectional area of the sheet cavity at • ,this position. . .~ . . * , ILv ’. - ,’ 4 4- ,*-. . ... . 4.4 K"% Nuclei and Cavitation Jean -Pierre Le Goff...the experiments, with analysing .- the results and with running computer programs. Thanks are also due to U H Pinto who developed a substantial part of

Quasi-Static Viscoelastic Finite Element Model of an Aircraft Tire

NASA Technical Reports Server (NTRS)

Johnson, Arthur R.; Tanner, John A.; Mason, Angela J.

1999-01-01

An elastic large displacement thick-shell mixed finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at the element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses is included in the mixed variational functional. The nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the nonlinear elastic equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to computationally simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Deterministic Stress Modeling of Hot Gas Segregation in a Turbine

NASA Technical Reports Server (NTRS)

Busby, Judy; Sondak, Doug; Staubach, Brent; Davis, Roger

1998-01-01

Simulation of unsteady viscous turbomachinery flowfields is presently impractical as a design tool due to the long run times required. Designers rely predominantly on steady-state simulations, but these simulations do not account for some of the important unsteady flow physics. Unsteady flow effects can be modeled as source terms in the steady flow equations. These source terms, referred to as Lumped Deterministic Stresses (LDS), can be used to drive steady flow solution procedures to reproduce the time-average of an unsteady flow solution. The goal of this work is to investigate the feasibility of using inviscid lumped deterministic stresses to model unsteady combustion hot streak migration effects on the turbine blade tip and outer air seal heat loads using a steady computational approach. The LDS model is obtained from an unsteady inviscid calculation. The LDS model is then used with a steady viscous computation to simulate the time-averaged viscous solution. Both two-dimensional and three-dimensional applications are examined. The inviscid LDS model produces good results for the two-dimensional case and requires less than 10% of the CPU time of the unsteady viscous run. For the three-dimensional case, the LDS model does a good job of reproducing the time-averaged viscous temperature migration and separation as well as heat load on the outer air seal at a CPU cost that is 25% of that of an unsteady viscous computation.

Role of Viscous Dissipative Processes on the Wetting of Textured Surfaces

PubMed Central

Grewal, H. S.; Nam Kim, Hong; Cho, Il-Joo; Yoon, Eui-Sung

2015-01-01

We investigate the role of viscous forces on the wetting of hydrophobic, semi-hydrophobic, and hydrophilic textured surfaces as second-order effects. We show that during the initial contact, the transition from inertia- to viscous-dominant regime occurs regardless of their surface topography and chemistry. Furthermore, we demonstrate the effect of viscosity on the apparent contact angle under quasi-static conditions by modulating the ratio of a water/glycerol mixture and show the effect of viscosity, especially on the semi-hydrophobic and hydrophobic textured substrates. The reason why the viscous force does not affect the apparent contact angle of the hydrophilic surface is explained based on the relationship between the disjoining pressure and surface chemistry. We further propose a wetting model that can predict the apparent contact angle of a liquid drop on a textured substrate by incorporating a viscous force component in the force balance equation. This model can predict apparent contact angles on semi-hydrophobic and hydrophobic textured surfaces exhibiting Wenzel state more accurately than the Wenzel model, indicating the importance of viscous forces in determining the apparent contact angle. The modified model can be applied for estimating the wetting properties of arbitrary engineered surfaces. PMID:26390958

A discrete geometric approach for simulating the dynamics of thin viscous threads

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

Audoly, B., E-mail: audoly@lmm.jussieu.fr; Clauvelin, N.; Brun, P.-T.

We present a numerical model for the dynamics of thin viscous threads based on a discrete, Lagrangian formulation of the smooth equations. The model makes use of a condensed set of coordinates, called the centerline/spin representation: the kinematic constraints linking the centerline's tangent to the orientation of the material frame is used to eliminate two out of three degrees of freedom associated with rotations. Based on a description of twist inspired from discrete differential geometry and from variational principles, we build a full-fledged discrete viscous thread model, which includes in particular a discrete representation of the internal viscous stress. Consistencymore » of the discrete model with the classical, smooth equations for thin threads is established formally. Our numerical method is validated against reference solutions for steady coiling. The method makes it possible to simulate the unsteady behavior of thin viscous threads in a robust and efficient way, including the combined effects of inertia, stretching, bending, twisting, large rotations and surface tension.« less

On the membrane approximation in isothermal film casting

NASA Astrophysics Data System (ADS)

Hagen, Thomas

2014-08-01

In this work, a one-dimensional model for isothermal film casting is studied. Film casting is an important engineering process to manufacture thin films and sheets from a highly viscous polymer melt. The model equations account for variations in film width and film thickness, and arise from thinness and kinematic assumptions for the free liquid film. The first aspect of our study is a rigorous discussion of the existence and uniqueness of stationary solutions. This objective is approached via the argument principle, exploiting the homotopy invariance of a family of analytic functions. As our second objective, we analyze the linearization of the governing equations about stationary solutions. It is shown that solutions for the associated boundary-initial value problem are given by a strongly continuous semigroup of bounded linear operators. To reach this result, we cast the relevant Cauchy problem in a more accessible form. These transformed equations allow us insight into the regularity of the semigroup, thus yielding the validity of the spectral mapping theorem for the semigroup and the spectrally determined growth property.

The effect of convective boundary condition on MHD mixed convection boundary layer flow over an exponentially stretching vertical sheet

NASA Astrophysics Data System (ADS)

Isa, Siti Suzilliana Putri Mohamed; Arifin, Norihan Md.; Nazar, Roslinda; Bachok, Norfifah; Ali, Fadzilah Md

2017-12-01

A theoretical study that describes the magnetohydrodynamic mixed convection boundary layer flow with heat transfer over an exponentially stretching sheet with an exponential temperature distribution has been presented herein. This study is conducted in the presence of convective heat exchange at the surface and its surroundings. The system is controlled by viscous dissipation and internal heat generation effects. The governing nonlinear partial differential equations are converted into ordinary differential equations by a similarity transformation. The converted equations are then solved numerically using the shooting method. The results related to skin friction coefficient, local Nusselt number, velocity and temperature profiles are presented for several sets of values of the parameters. The effects of the governing parameters on the features of the flow and heat transfer are examined in detail in this study.

Current Sheets in Pulsar Magnetospheres and Winds: Particle Acceleration and Pulsed Gamma Ray Emission

NASA Astrophysics Data System (ADS)

Arons, Jonathan

The research proposed addresses understanding of the origin of non-thermal energy in the Universe, a subject beginning with the discovery of Cosmic Rays and continues, including the study of relativistic compact objects - neutron stars and black holes. Observed Rotation Powered Pulsars (RPPs) have rotational energy loss implying they have TeraGauss magnetic fields and electric potentials as large as 40 PetaVolts. The rotational energy lost is reprocessed into particles which manifest themselves in high energy gamma ray photon emission (GeV to TeV). Observations of pulsars from the FERMI Gamma Ray Observatory, launched into orbit in 2008, have revealed 130 of these stars (and still counting), thus demonstrating the presence of efficient cosmic accelerators within the strongly magnetized regions surrounding the rotating neutron stars. Understanding the physics of these and other Cosmic Accelerators is a major goal of astrophysical research. A new model for particle acceleration in the current sheets separating the closed and open field line regions of pulsars' magnetospheres, and separating regions of opposite magnetization in the relativistic winds emerging from those magnetopsheres, will be developed. The currents established in recent global models of the magnetosphere will be used as input to a magnetic field aligned acceleration model that takes account of the current carrying particles' inertia, generalizing models of the terrestrial aurora to the relativistic regime. The results will be applied to the spectacular new results from the FERMI gamma ray observatory on gamma ray pulsars, to probe the physics of the generation of the relativistic wind that carries rotational energy away from the compact stars, illuminating the whole problem of how compact objects can energize their surroundings. The work to be performed if this proposal is funded involves extending and developing concepts from plasma physics on dissipation of magnetic energy in thin sheets of electric current that separate regions of differing magnetization into the domain of highly relativistic magnetic fields - those with energy density large compared to the rest mass energy of the charged particles - the plasma - caught in that field. The investigators will create theoretical and computational models of the magnetic dissipation - a form of viscous flow in the thin sheets of electric current that form in the magnetized regions around the rotating stars - using Particle in-Cell plasma simulations. These simulations use a large computer to solve the equations of motion of many charged particles - millions to billions in the research that will be pursued - to unravel the dissipation of those fields and the acceleration of beams of particles in the thin sheets. The results will be incorporated into macroscopic MHD models of the magnetic structures around the stars which determine the location and strength of the current sheets, so as to model and analyze the pulsed gamma ray emission seen from hundreds of Rotation Powered Pulsars. The computational models will be assisted by ``pencil and paper'' theoretical modeling designed to motivate and interpret the computer simulations, and connect them to the observations.

Investigating Different Patterns of Slab Deformation in the Lower Mantle

NASA Astrophysics Data System (ADS)

Zhang, J.; McNamara, A. K.

2017-12-01

The geometry of slabs within the upper mantle have been relatively well-imaged by tomography and regional seismic studies; however, the style of slab deformation in the lower mantle remains poorly understood. Although tomography models reveal that the lower mantle beneath paleo-subduction regions are faster-than-average, the resolution is not high enough to resolve how slabs are actually deforming there. Slabs have long been hypothesized as viscous, tabular sheets that subduct at the surface, descend through the mantle, and impinge on the core-mantle boundary (CMB). Geodynamical studies have shown a wide range of possible deformational behaviors, ranging from stiff, buckling slabs to more-ductile masses of accumulating slab material undergoing pure shear. Of particular interest is how rheology and 3D spherical geometry control the shape and deformational style of slabs as they descend deeper into the mantle. We performed high resolution 3D spherical calculations to explore slab deformation in deep mantle as a function of slab strength. In our model, kinematic velocity boundary conditions are imposed on the surface to simulate a moving plate which guides the formation of a subducting slab. In addition, a viscosity jump at the transition zone is applied. We find that although a slab subducts as a large tabular sheet from the surface, it doesn't always maintain such geometry. Instead, it typically breaks apart into a few smaller and narrower sheets which can even turn into cylindrical-shaped downwelling after subducting into deep mantle. Since seismic anisotropy is hypothesized to originate from crystal preferred orientation (CPO) in a slab when it impinges on the CMB and is predicted with significant help of time-dependent deformation information from the geodynamic models, our findings on lower mantle slab deformation patterns may enhance the understanding towards the cause of characteristic patterns of predicted seismic anisotropy.

A depth-of-field limited particle image velocimetry technique applied to oscillatory boundary layer flow over a porous bed

NASA Astrophysics Data System (ADS)

Lara, J. L.; Cowen, E. A.; Sou, I. M.

2002-06-01

Boundary layer flows are ubiquitous in the environment, but their study is often complicated by their thinness, geometric irregularity and boundary porosity. In this paper, we present an approach to making laboratory-based particle image velocimetry (PIV) measurements in these complex flow environments. Clear polycarbonate spheres were used to model a porous and rough bed. The strong curvature of the spheres results in a diffuse volume illuminated region instead of the more traditional finite and thin light sheet illuminated region, resulting in the imaging of both in-focus and significantly out-of-focus particles. Results of a traditional cross-correlation-based PIV-type analysis of these images demonstrate that the mean and turbulent features of an oscillatory boundary layer driven by a free-surface wave over an irregular-shaped porous bed can be robustly measured. Measurements of the mean flow, turbulent intensities, viscous and turbulent stresses are presented and discussed. Velocity spectra have been calculated showing an inertial subrange confirming that the PIV analysis is sufficiently robust to extract turbulence. The presented technique is particularly well suited for the study of highly dynamic free-surface flows that prevent the delivery of the light sheet from above the bed, such as swash flows.

The role of membrane-like stresses in determining the stability and sensitivity of the Antarctic ice sheets: back pressure and grounding line motion.

PubMed

Hindmarsh, Richard C A

2006-07-15

Membrane stresses act along thin bodies which are relatively well lubricated on both surfaces. They operate in ice sheets because the bottom is either sliding, or is much less viscous than the top owing to stress and heat softening of the basal ice. Ice streams flow over very well lubricated beds, and are restrained at their sides. The ideal of the perfectly slippery bed is considered in this paper, and the propagation of mechanical effects along an ice stream considered by applying spatially varying horizontal body forces. Propagation distances depend sensitively on the rheological index, and can be very large for ice-type rheologies.A new analytical solution for ice-shelf profiles and grounded tractionless stream profiles is presented, which show blow up of the profile in a finite distance upstream at locations where the flux is non-zero. This is a feature of an earlier analytical solution for a floating shelf.The length scale of decay of membrane stresses from the grounding line is investigated through scale analysis. In ice sheets, such effects decay over distances of several tens of kilometres, creating a vertical boundary layer between sheet flow and shelf flow, where membrane stresses adjust. Bounded, physically reasonable steady surface profiles only exist conditionally in this boundary layer. Where bounded steady profiles exist, adjacent profile equilibria for the whole ice sheet corresponding to different grounded areas occur (neutral equilibrium). If no solution in the boundary layer can exist, the ice-sheet profile must change.The conditions for existence can be written in terms of whether the basal rate factor (sliding or internal deformation) is too large to permit a steady solution. The critical value depends extremely sensitively on ice velocity and the back stress applied at the grounding line. High ice velocity and high stress both favour the existence of solutions and stability. Changes in these parameters can cause the steady solution existence criterion to be traversed, and the ice-sheet dynamics to change.A finite difference model which represents both neutral equilibrium and the dynamical transition is presented, and preliminary investigations into its numerical sensitivity are carried out. Evidence for the existence of a long wavelength instability is presented through the solution of a numerical eigenproblem, which will hamper predictability.

Swimming in an anisotropic fluid: How speed depends on alignment angle

NASA Astrophysics Data System (ADS)

Shi, Juan; Powers, Thomas R.

2017-12-01

Orientational order in a fluid affects the swimming behavior of flagellated microorganisms. For example, bacteria tend to swim along the director in lyotropic nematic liquid crystals. To better understand how anisotropy affects propulsion, we study the problem of a sheet supporting small-amplitude traveling waves, also known as the Taylor swimmer, in a nematic liquid crystal. For the case of weak anchoring of the nematic director at the swimmer surface and in the limit of a minimally anisotropic model, we calculate the swimming speed as a function of the angle between the swimmer and the nematic director. The effect of the anisotropy can be to increase or decrease the swimming speed, depending on the angle of alignment. We also show that elastic torque dominates the viscous torque for small-amplitude waves and that the torque tends to align the swimmer along the local director.

Theory of buttressed marine ice sheet dynamics and its application to the assessment of tipping-point conditions

NASA Astrophysics Data System (ADS)

Pegler, S.

2017-12-01

Understanding the fate of the West Antarctic Ice Sheet is constrained by difficulties of resolving the buttressing effect of ice shelves and its dynamic response to grounding-line migration. This effect may be responsible for protecting a large majority of outlet glaciers in Antarctica against surging into the ocean. I present a theoretical methodology for assessing the positions and stability of grounding lines that incorporates closed-form, dynamic descriptions of ice-shelf buttressing and extensional (stretching) viscous stresses. The method is applied to assess the conditions for grounding-line tipping points. Such points are shown to produce abrupt `cliff-edge' transitions to runaway retreat, representing the so-called marine ice sheet instability. Depending on the bed profile, melt and calving rates, a tipping point can either lie very near to a local maximum in the bed topography or potentially far upstream of it, along a reverse bed. The model predictions for both wide and narrow embayments are validated by numerical simulations and laboratory experiments. A case study of Pine Island Glacier indicates the possibility for long-term stabilisation, with the analytical method affording an extensive exploration of scenarios. The theory also elucidates a mode of grounding-line migration controlled entirely by the determinants of the ice-shelf buttressing force, with a loss of sensitivity to basal conditions, contrasting with the conclusion from one-dimensional theory that the ice shelf is irrelevant. The results provide an interpretive framework for understanding grounding-line dynamics, its coupling with ice-shelf dynamics, an efficient exploration of parameter variation, and a complement to large-scale simulation.

Magma Fragmentation

NASA Astrophysics Data System (ADS)

Gonnermann, Helge M.

2015-05-01

Magma fragmentation is the breakup of a continuous volume of molten rock into discrete pieces, called pyroclasts. Because magma contains bubbles of compressible magmatic volatiles, decompression of low-viscosity magma leads to rapid expansion. The magma is torn into fragments, as it is stretched into hydrodynamically unstable sheets and filaments. If the magma is highly viscous, resistance to bubble growth will instead lead to excess gas pressure and the magma will deform viscoelastically by fracturing like a glassy solid, resulting in the formation of a violently expanding gas-pyroclast mixture. In either case, fragmentation represents the conversion of potential energy into the surface energy of the newly created fragments and the kinetic energy of the expanding gas-pyroclast mixture. If magma comes into contact with external water, the conversion of thermal energy will vaporize water and quench magma at the melt-water interface, thus creating dynamic stresses that cause fragmentation and the release of kinetic energy. Lastly, shear deformation of highly viscous magma may cause brittle fractures and release seismic energy.

From viscous to elastic sheets: Dynamics of smectic freely floating films

NASA Astrophysics Data System (ADS)

Harth, Kirsten; May, Kathrin; Trittel, Torsten; Stannarius, Ralf

2015-03-01

Oscillations and rupture of bubbles, composed of an inner fluid separated from an outer fluid by a membrane, represent an old but still immensely active field of research. Membrane properties except surface tension are often neglected for simple fluid films (e.g. soap bubbles), whereas they govern the dynamics in systems with more complex membranes (e.g. vesicles). Due to their layered phase structure, smectic liquid crystals can form stable, uniform and easy-to handle fluid films of immense aspect ratios. Recently, freely floating bubbles detached from a support were prepared. We analyze the relaxation from strongly non-spherical shapes and the rupture dynamics of such bubbles using high-speed video recordings. Peculiar dynamics intermediate between those of simple viscous fluid films and an elastic response emerge: Oscillations, slowed relaxation and even the formation of wrinkles and extrusions. We characterize these phenomena and propose explanations. We acknowledge funding by the German Aerospace Center DLR within Project OASIS-CO and German Science Foundation Project STA 425-28.

In-vivo viscous properties of the heel pad by stress-relaxation experiment based on a spherical indentation.

PubMed

Suzuki, Ryo; Ito, Kohta; Lee, Taeyong; Ogihara, Naomichi

2017-12-01

Identifying the viscous properties of the plantar soft tissue is crucial not only for understanding the dynamic interaction of the foot with the ground during locomotion, but also for development of improved footwear products and therapeutic footwear interventions. In the present study, the viscous and hyperelastic material properties of the plantar soft tissue were experimentally identified using a spherical indentation test and an analytical contact model of the spherical indentation test. Force-relaxation curves of the heel pads were obtained from the indentation experiment. The curves were fit to the contact model incorporating a five-element Maxwell model to identify the viscous material parameters. The finite element method with the experimentally identified viscoelastic parameters could successfully reproduce the measured force-relaxation curves, indicating the material parameters were correctly estimated using the proposed method. Although there are some methodological limitations, the proposed framework to identify the viscous material properties may facilitate the development of subject-specific finite element modeling of the foot and other biological materials. Copyright © 2017 IPEM. Published by Elsevier Ltd. All rights reserved.

Dynamic analysis of submerged microscale plates: the effects of acoustic radiation and viscous dissipation

PubMed Central

Ma, Xianghong

2016-01-01

The aim of this paper is to study the dynamic characteristics of micromechanical rectangular plates used as sensing elements in a viscous compressible fluid. A novel modelling procedure for the plate–fluid interaction problem is developed on the basis of linearized Navier–Stokes equations and no-slip conditions. Analytical expression for the fluid-loading impedance is obtained using a double Fourier transform approach. This modelling work provides us an analytical means to study the effects of inertial loading, acoustic radiation and viscous dissipation of the fluid acting on the vibration of microplates. The numerical simulation is conducted on microplates with different boundary conditions and fluids with different viscosities. The simulation results reveal that the acoustic radiation dominates the damping mechanism of the submerged microplates. It is also proved that microplates offer better sensitivities (Q-factors) than the conventional beam type microcantilevers being mass sensing platforms in a viscous fluid environment. The frequency response features of microplates under highly viscous fluid loading are studied using the present model. The dynamics of the microplates with all edges clamped are less influenced by the highly viscous dissipation of the fluid than the microplates with other types of boundary conditions. PMID:27118914

Dynamic analysis of submerged microscale plates: the effects of acoustic radiation and viscous dissipation.

PubMed

Wu, Zhangming; Ma, Xianghong

2016-03-01

The aim of this paper is to study the dynamic characteristics of micromechanical rectangular plates used as sensing elements in a viscous compressible fluid. A novel modelling procedure for the plate-fluid interaction problem is developed on the basis of linearized Navier-Stokes equations and no-slip conditions. Analytical expression for the fluid-loading impedance is obtained using a double Fourier transform approach. This modelling work provides us an analytical means to study the effects of inertial loading, acoustic radiation and viscous dissipation of the fluid acting on the vibration of microplates. The numerical simulation is conducted on microplates with different boundary conditions and fluids with different viscosities. The simulation results reveal that the acoustic radiation dominates the damping mechanism of the submerged microplates. It is also proved that microplates offer better sensitivities (Q-factors) than the conventional beam type microcantilevers being mass sensing platforms in a viscous fluid environment. The frequency response features of microplates under highly viscous fluid loading are studied using the present model. The dynamics of the microplates with all edges clamped are less influenced by the highly viscous dissipation of the fluid than the microplates with other types of boundary conditions.

Linear instabilities of a planar liquid sheet in a static electric field for intermediate relaxation and convection of surface charges

NASA Astrophysics Data System (ADS)

Yoshinaga, Takao

2018-04-01

Linear temporal instabilities of a two-dimensional planar liquid sheet in a static electric field are investigated when the relaxation and convection of surface electric charges are considered. Both viscous sheet liquid and inviscid surrounding liquid are placed between two parallel sheath walls, on which an external electric field is imposed. In particular, effects of the electric Peclet number {Pe} (charge relaxation time/convection time) and the electric Euler number Λ (electric pressure/liquid inertial) on the instabilities are emphasized for the symmetric and antisymmetric deformations of the sheet. It is found that the unstable mode is composed of the aerodynamic and electric modes, which are merged with each other for the symmetric deformation and separated for the antisymmetric deformation. For the symmetric deformation, the combined mode is more destabilized with the decrease of {Pe} and the increase of Λ. On the other hand, for the antisymmetric deformation, the electric mode is more destabilized and the aerodynamic mode is left unchanged with the decrease of {Pe}, while the electric mode is more destabilized but the aerodynamic mode is more stabilized with the increase of Λ. It is also found for both symmetric and antisymmetric deformations that the instabilities are most suppressed when {σ }R≃ 1/{ε }P ({σ }R: conductivity ratio of the surrounding to the sheet liquid, {ε }P: permittivity ratio of the sheet to the surrounding liquid), whose trend of the instabilities is more enhanced with the decrease of {Pe} except for vanishingly small {Pe}.

Computational studies of an impulsively started viscous flow

NASA Technical Reports Server (NTRS)

Davis, Sanford S.

1988-01-01

Progress in validating incompressible Navier-Stokes codes is described using a predictor/corrector scheme. The flow field under study is the impulsive start of a circular cylinder and the unsteady evolution of the separation bubble. In the current code, a uniform asymptotic expansion is used as an initial condition in order to correctly capture the initial growth of the vortex sheet. Volocity fields at selected instants of time are decomposed into vectorial representations of Navier-Stokes equations which are then used to analyze dominant contributions in the boundary-layer region.

Impact of an interbedded viscous décollement on the structural and kinematic coupling in fold-and-thrust belts: Insights from analogue modeling

NASA Astrophysics Data System (ADS)

Borderie, Sandra; Graveleau, Fabien; Witt, César; Vendeville, Bruno C.

2018-01-01

Fold-and-thrust belts (FTBs) can be segmented both across and along strike because of various factors including tectonic and stratigraphic inheritance. In this study, we investigated along/across-strike structural interactions in a FTB propagating toward a foreland which displays contrasted lithological sequences. A set of analogue models was performed in a compressional box where a single viscous level of varying width was interbedded within a frictional series. The tectonic interaction between the viscous and the frictional provinces was tested both along and across strike. Results indicate that a frictional province influences the along-strike tectonic evolution of an adjacent viscous province. This influence decreases when the width of the viscous province increases. The frictional provinces control the taper, structural style, obliquity of the structures' trend and kinematics of the shallow deformation front of the viscous province. Results evidence how far a frictional province can impact the deformation of an adjacent viscous province. For frictional-viscous wedges, it appears that the critical taper theory, which is generally applied in 2-D, should be likely considered in terms of 3-D. Moreover, the kinematics of the deep deformation front shows mutual influences between the adjacent viscous and frictional provinces. Experimental results are compared to natural examples in the Kuqa Basin (Southern Tian Shan, China) and the Salt Range (Pakistan), and give an insight to a better understanding of the dynamics of fold-and-thrust belts bearing a viscous décollement, such as salt.

Flow-induced flutter in a wall-bounded elastic sheet

NASA Astrophysics Data System (ADS)

Weidman, M. S.; Argentina, M.; Hosoi, A. E.; Mahadevan, L.

2004-11-01

Inspired by voice production in natural and artificial systems, we consider the flow between a long but finite flexible elastic sheet and a rigid wall close to it. We derive evolution equations for the coupled dynamics of the fluid and solid in two limits corresponding to the viscously dominated and inertially dominated regimes of the flow. In both situations, the inertia of the solid remains important. We show that a long wavelength instability via a 1:1 resonance mechanism arises in both situations when the flow rate is increased beyond a critical threshold. We also compare the results of our analytical, numerical and scaling calculations with those of simple experiments. Finally we comment on the rich nonlinear dynamics of these systems which suggest that at least some aspects of voice and song production may be more a manifestation of physics rather than neurophysiology.

Global existence of the three-dimensional viscous quantum magnetohydrodynamic model

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

Yang, Jianwei, E-mail: yangjianwei@ncwu.edu.cn; Ju, Qiangchang, E-mail: qiangchang-ju@yahoo.com

2014-08-15

The global-in-time existence of weak solutions to the viscous quantum Magnetohydrodynamic equations in a three-dimensional torus with large data is proved. The global existence of weak solutions to the viscous quantum Magnetohydrodynamic equations is shown by using the Faedo-Galerkin method and weak compactness techniques.

Derivation of revised formulae for eddy viscous forces used in the ocean general circulation model

NASA Technical Reports Server (NTRS)

Chou, Ru Ling

1988-01-01

Presented is a re-derivation of the eddy viscous dissipation tensor commonly used in present oceanographic general circulation models. When isotropy is imposed, the currently-used form of the tensor fails to return to the laplacian operator. In this paper, the source of this error is identified in a consistent derivation of the tensor in both rectangular and earth spherical coordinates, and the correct form of the eddy viscous tensor is presented.

Viscoelastic responses of a hard transition zone - Effects on postglacial uplifts and rotational signatures

NASA Technical Reports Server (NTRS)

Spada, Giorgio; Sabadini, Roberto; Yuen, David A.

1991-01-01

A five-layer viscoelastic spherical model is used to calculate the transient displacements of postglacial rebound, the induced polar motions, and the temporal variations of the geopotential up to degree 8 of the zonal coefficients. Two models - one with two viscoelastic layers separated at 670 km, and the other with three layers in which a hard garnet layer lies between the upper and lower mantle - are compared. Forward modeling shows that it may be possible to discern the presence of a hard garnet layer with a viscosity of at least ten times greater than the upper mantle, on the basis of uplift data near the center of the former Laurentide ice-sheet and from polar wander and j2 data. Temporal variations of higher gravity harmonics, such as j6 and j8, can potentially place even tighter constraints on the rheological properties of the hard transition zone. A lower mantle viscosity between 2 and 4 x 10 to the 22nd Pa is generally preferred in models with a garnet layer which may be as large as 50 times more viscous than the upper mantle.

Magnetic and velocity fields in a dynamo operating at extremely small Ekman and magnetic Prandtl numbers

NASA Astrophysics Data System (ADS)

Šimkanin, Ján; Kyselica, Juraj

2017-12-01

Numerical simulations of the geodynamo are becoming more realistic because of advances in computer technology. Here, the geodynamo model is investigated numerically at the extremely low Ekman and magnetic Prandtl numbers using the PARODY dynamo code. These parameters are more realistic than those used in previous numerical studies of the geodynamo. Our model is based on the Boussinesq approximation and the temperature gradient between upper and lower boundaries is a source of convection. This study attempts to answer the question how realistic the geodynamo models are. Numerical results show that our dynamo belongs to the strong-field dynamos. The generated magnetic field is dipolar and large-scale while convection is small-scale and sheet-like flows (plumes) are preferred to a columnar convection. Scales of magnetic and velocity fields are separated, which enables hydromagnetic dynamos to maintain the magnetic field at the low magnetic Prandtl numbers. The inner core rotation rate is lower than that in previous geodynamo models. On the other hand, dimensional magnitudes of velocity and magnetic fields and those of the magnetic and viscous dissipation are larger than those expected in the Earth's core due to our parameter range chosen.

Experiment study of mud to the moving process influent about viscous debris flow along slope

NASA Astrophysics Data System (ADS)

Jun, JiXian; Ying, Liang; Li, Pan Hua; Qiang, OuGuo

2018-01-01

Mud is the main component of viscous debris flow. The physical model experiments of viscous debris flow were carried out through the mixing mud with different density and fixed components of coarse particles. The width, longitudinal movement distance and motion velocity were recorded by video cameras during experiment. Through viscous debris flow physical model experiments, the influence of mud to transverse width, longitudinal movement distance and motion velocity was discussed. The physical model experiment results show that the motion forms change from inviscid particle flow to viscous debris flow and to the whole mass sliding with the increase of mud density; the width and the length along the slope decrease with mud density increasing; the movement process has classified phenomena about viscous debris flow composed by different mud densities: the velocity increases rapidly with time and the change gradient is steady when the density of mud is lower than 1.413g/cm3; the movement process can be divided into two stages when the density of mud is higher than 1.413g/cm3: the movement velocity is lower and the gradient change is small in the initial stage; but in the second stage, the movement velocity increases quickly, and the gradient is higher than the first stage, and with steady value.

A new Eulerian model for viscous and heat conducting compressible flows

NASA Astrophysics Data System (ADS)

Svärd, Magnus

2018-09-01

In this article, a suite of physically inconsistent properties of the Navier-Stokes equations, associated with the lack of mass diffusion and the definition of velocity, is presented. We show that these inconsistencies are consequences of the Lagrangian derivation that models viscous stresses rather than diffusion. A new model for compressible and diffusive (viscous and heat conducting) flows of an ideal gas, is derived in a purely Eulerian framework. We propose that these equations supersede the Navier-Stokes equations. A few numerical experiments demonstrate some differences and similarities between the new system and the Navier-Stokes equations.

Zonal flow evolution and overstability in accretion discs

NASA Astrophysics Data System (ADS)

Vanon, R.; Ogilvie, G. I.

2017-04-01

This work presents a linear analytical calculation on the stability and evolution of a compressible, viscous self-gravitating (SG) Keplerian disc with both horizontal thermal diffusion and a constant cooling time-scale when an axisymmetric structure is present and freely evolving. The calculation makes use of the shearing sheet model and is carried out for a range of cooling times. Although the solutions to the inviscid problem with no cooling or diffusion are well known, it is non-trivial to predict the effect caused by the introduction of cooling and of small diffusivities; this work focuses on perturbations of intermediate wavelengths, therefore representing an extension to the classical stability analysis on thermal and viscous instabilities. For density wave modes, the analysis can be simplified by means of a regular perturbation analysis; considering both shear and thermal diffusivities, the system is found to be overstable for intermediate and long wavelengths for values of the Toomre parameter Q ≲ 2; a non-SG instability is also detected for wavelengths ≳18H, where H is the disc scale-height, as long as γ ≲ 1.305. The regular perturbation analysis does not, however, hold for the entropy and potential vorticity slow modes as their ideal growth rates are degenerate. To understand their evolution, equations for the axisymmetric structure's amplitudes in these two quantities are analytically derived and their instability regions obtained. The instability appears boosted by increasing the value of the adiabatic index and of the Prandtl number, while it is quenched by efficient cooling.

On spontaneous formation of current sheets: Untwisted magnetic fields

NASA Astrophysics Data System (ADS)

Bhattacharyya, R.; Low, B. C.; Smolarkiewicz, P. K.

2010-11-01

This is a study of the spontaneous formation of electric current sheets in an incompressible viscous fluid with perfect electrical conductivity, governed by the magnetohydrodynamic Navier-Stokes equations. Numerical solutions to two initial value problems are presented for a three-dimensional, periodic, untwisted magnetic field evolving, with no change in magnetic topology under the frozen-in condition and at characteristic fluid Reynolds numbers of the order of 500, from a nonequilibrium initial state with the fluid at rest. The evolution converts magnetic free energy into kinetic energy to be all dissipated away by viscosity so that the field settles into a minimum-energy, static equilibrium. The solutions demonstrate that, as a consequence of the frozen-in condition, current sheets must form during the evolution despite the geometric simplicity of the prescribed initial fields. In addition to the current sheets associated with magnetic neutral points and field reversal layers, other sheets not associated with such magnetic features are also in evidence. These current sheets form on magnetic flux surfaces. This property is used to achieve a high degree of the frozen-in condition in the simulations, by describing the magnetic field entirely in terms of the advection of its flux surfaces and integrating the resulting governing equations with a customized version of a general-purpose high-resolution (viz., nonoscillatory) hydrodynamical simulation code EULAG [J. M. Prusa et al., Comput. Fluids 37, 1193 (2008)]. Incompressibility imposes the additional global constraint that the flux surfaces must evolve with no change in the spatial volumes they enclose. In this approach, current sheet formation is demonstrated graphically by the progressive pressing together of suitably selected flux surfaces until their separation has diminished below the minimal resolved distance on a fixed grid. The frozen-in condition then fails in the simulation as the field reconnects through an effecting numerical resistivity. The principal results are related to the Parker theory of current-sheet formation and dissipation in the solar corona.

A Viscoelastic Hybrid Shell Finite Element

NASA Technical Reports Server (NTRS)

Johnson, Arthur

1999-01-01

An elastic large displacement thick-shell hybrid finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at he element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses in included in the mixed variational functional. Nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to numerically simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

From viscous to elastic sheets: Dynamics of smectic freely floating films

NASA Astrophysics Data System (ADS)

Stannarius, Ralf; Harth, Kirsten; May, Kathrin; Trittel, Torsten

The dynamics of droplets and bubbles, particularly on microscopic scales, are of considerable importance in biological, environmental, and technical contexts. Soap bubbles, vesicles and components of biological cells are well known examples where the dynamic behavior is significantly influenced by the properties of thin membranes enclosed by fluids. Two-dimensional membrane motions couple to 3D shape transformations. Smectic liquid crystal mesogens form phases with internal molecular layer order. Free-standing films are easily prepared from this class of materials. They represent simple model systems for membrane dynamics and pattern formation in a quasi two-dimensional fluid. These films are usually spanned over a frame, and they can be inflated to bubbles on a support. Recently, closed microscopic shells of liquid-crystalline materials suspended in an outer fluid without contact to a solid support have been introduced and studied. With a special technique, we prepare millimetre to centimetre sized smectic bubbles in air (similar to soap bubbles). Their distinct feature is the fact that any change of surface area is coupled to a restructuring of the layers in the membrane. High-speed cameras are used to observe the shape transformations of freely floating bubbles from a distorted initial shape to a sphere. Bursting dynamics are recorded and compared to models. Most strikingly, an unpreceded cross-over from inviscid to viscous and elastic behaviour with increasing thickness of the membrane is found: Whereas thin bubbles behave almost like inviscid fluids, the relaxation dynamics slows down considerably for larger film thicknesses. Surface wrinkling and formation of extrusions are observed. We will present a characterization and an expalantion for the above phenomena.

Entropy Analysis in Mixed Convection MHD flow of Nanofluid over a Non-linear Stretching Sheet

NASA Astrophysics Data System (ADS)

Matin, Meisam Habibi; Nobari, Mohammad Reza Heirani; Jahangiri, Pouyan

This article deals with a numerical study of entropy analysis in mixed convection MHD flow of nanofluid over a non-linear stretching sheet taking into account the effects of viscous dissipation and variable magnetic field. The nanofluid is made of such nano particles as SiO2 with pure water as a base fluid. To analyze the problem, at first the boundary layer equations are transformed into non-linear ordinary equations using a similarity transformation. The resultant equations are then solved numerically using the Keller-Box scheme based on the implicit finite-difference method. The effects of different non-dimensional governing parameters such as magnetic parameter, nanoparticles volume fraction, Nusselt, Richardson, Eckert, Hartman, Brinkman, Reynolds and entropy generation numbers are investigated in details. The results indicate that increasing the nano particles to the base fluids causes the reduction in shear forces and a decrease in stretching sheet heat transfer coefficient. Also, decreasing the magnetic parameter and increasing the Eckert number result in improves heat transfer rate. Furthermore, the surface acts as a strong source of irreversibility due to the higher entropy generation number near the surface.

Universal Rim Thickness in Unsteady Sheet Fragmentation.

PubMed

Wang, Y; Dandekar, R; Bustos, N; Poulain, S; Bourouiba, L

2018-05-18

Unsteady fragmentation of a fluid bulk into droplets is important for epidemiology as it governs the transport of pathogens from sneezes and coughs, or from contaminated crops in agriculture. It is also ubiquitous in industrial processes such as paint, coating, and combustion. Unsteady fragmentation is distinct from steady fragmentation on which most theoretical efforts have been focused thus far. We address this gap by studying a canonical unsteady fragmentation process: the breakup from a drop impact on a finite surface where the drop fluid is transferred to a free expanding sheet of time-varying properties and bounded by a rim of time-varying thickness. The continuous rim destabilization selects the final spray droplets, yet this process remains poorly understood. We combine theory with advanced image analysis to study the unsteady rim destabilization. We show that, at all times, the rim thickness is governed by a local instantaneous Bond number equal to unity, defined with the instantaneous, local, unsteady rim acceleration. This criterion is found to be robust and universal for a family of unsteady inviscid fluid sheet fragmentation phenomena, from impacts of drops on various surface geometries to impacts on films. We discuss under which viscous and viscoelastic conditions the criterion continues to govern the unsteady rim thickness.

Finite Element Solution of Unsteady Mixed Convection Flow of Micropolar Fluid over a Porous Shrinking Sheet

PubMed Central

Gupta, Diksha; Singh, Bani

2014-01-01

The objective of this investigation is to analyze the effect of unsteadiness on the mixed convection boundary layer flow of micropolar fluid over a permeable shrinking sheet in the presence of viscous dissipation. At the sheet a variable distribution of suction is assumed. The unsteadiness in the flow and temperature fields is caused by the time dependence of the shrinking velocity and surface temperature. With the aid of similarity transformations, the governing partial differential equations are transformed into a set of nonlinear ordinary differential equations, which are solved numerically, using variational finite element method. The influence of important physical parameters, namely, suction parameter, unsteadiness parameter, buoyancy parameter and Eckert number on the velocity, microrotation, and temperature functions is investigated and analyzed with the help of their graphical representations. Additionally skin friction and the rate of heat transfer have also been computed. Under special conditions, an exact solution for the flow velocity is compared with the numerical results obtained by finite element method. An excellent agreement is observed for the two sets of solutions. Furthermore, to verify the convergence of numerical results, calculations are conducted with increasing number of elements. PMID:24672310

Universal Rim Thickness in Unsteady Sheet Fragmentation

NASA Astrophysics Data System (ADS)

Wang, Y.; Dandekar, R.; Bustos, N.; Poulain, S.; Bourouiba, L.

2018-05-01

Unsteady fragmentation of a fluid bulk into droplets is important for epidemiology as it governs the transport of pathogens from sneezes and coughs, or from contaminated crops in agriculture. It is also ubiquitous in industrial processes such as paint, coating, and combustion. Unsteady fragmentation is distinct from steady fragmentation on which most theoretical efforts have been focused thus far. We address this gap by studying a canonical unsteady fragmentation process: the breakup from a drop impact on a finite surface where the drop fluid is transferred to a free expanding sheet of time-varying properties and bounded by a rim of time-varying thickness. The continuous rim destabilization selects the final spray droplets, yet this process remains poorly understood. We combine theory with advanced image analysis to study the unsteady rim destabilization. We show that, at all times, the rim thickness is governed by a local instantaneous Bond number equal to unity, defined with the instantaneous, local, unsteady rim acceleration. This criterion is found to be robust and universal for a family of unsteady inviscid fluid sheet fragmentation phenomena, from impacts of drops on various surface geometries to impacts on films. We discuss under which viscous and viscoelastic conditions the criterion continues to govern the unsteady rim thickness.

Calculation of unsteady transonic flows with mild separation by viscous-inviscid interaction

NASA Technical Reports Server (NTRS)

Howlett, James T.

1992-01-01

This paper presents a method for calculating viscous effects in two- and three-dimensional unsteady transonic flow fields. An integral boundary-layer method for turbulent viscous flow is coupled with the transonic small-disturbance potential equation in a quasi-steady manner. The viscous effects are modeled with Green's lag-entrainment equations for attached flow and an inverse boundary-layer method for flows that involve mild separation. The boundary-layer method is used stripwise to approximate three-dimensional effects. Applications are given for two-dimensional airfoils, aileron buzz, and a wing planform. Comparisons with inviscid calculations, other viscous calculation methods, and experimental data are presented. The results demonstrate that the present technique can economically and accurately calculate unsteady transonic flow fields that have viscous-inviscid interactions with mild flow separation.

Advanced Space Propulsion System Flowfield Modeling

NASA Technical Reports Server (NTRS)

Smith, Sheldon

1998-01-01

Solar thermal upper stage propulsion systems currently under development utilize small low chamber pressure/high area ratio nozzles. Consequently, the resulting flow in the nozzle is highly viscous, with the boundary layer flow comprising a significant fraction of the total nozzle flow area. Conventional uncoupled flow methods which treat the nozzle boundary layer and inviscid flowfield separately by combining the two calculations via the influence of the boundary layer displacement thickness on the inviscid flowfield are not accurate enough to adequately treat highly viscous nozzles. Navier Stokes models such as VNAP2 can treat these flowfields but cannot perform a vacuum plume expansion for applications where the exhaust plume produces induced environments on adjacent structures. This study is built upon recently developed artificial intelligence methods and user interface methodologies to couple the VNAP2 model for treating viscous nozzle flowfields with a vacuum plume flowfield model (RAMP2) that is currently a part of the Plume Environment Prediction (PEP) Model. This study integrated the VNAP2 code into the PEP model to produce an accurate, practical and user friendly tool for calculating highly viscous nozzle and exhaust plume flowfields.

Learning from jellyfish: Fluid transport in muscular pumps at intermediate Reynolds numbers

NASA Astrophysics Data System (ADS)

Nawroth, Janna; Dabiri, John

2010-11-01

Biologically inspired hydrodynamic propulsion and maneuvering strategies promise the advancement of medical implants and minimally invasive clinical tools. We have chosen juvenile jellyfish as a model system for investigating fluid dynamics and morphological properties underlying fluid transport by a muscular pump at intermediate Reynolds numbers. Recently we have described how natural variations in viscous forces are balanced by changes in jellyfish body shape (phenotypic plasticity), to the effect of facilitating efficient body-fluid interaction. Complementing these studies in our live model organisms, we are also engaged in engineering an artificial jellyfish, that is, a jellyfish-inspired construct of a flexible plastic sheet actuated by a monolayer of rat cardiomyocytes. The main challenges here are (1) to derive a body shape and deformation suitable for effective fluid transport under physiological conditions, (2) to understand the mechanical properties of the muscular film and derive a design capable of the desired deformation, (3) to master the proper alignment and timely contraction of the muscle component needed to achieve the desired deformation, and (4) to evaluate the performance of the design.

Theoretical modeling of electron mobility in superfluid {sup 4}He

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

Aitken, Frédéric; Bonifaci, Nelly; Haeften, Klaus von

The Orsay-Trento bosonic density functional theory model is extended to include dissipation due to the viscous response of superfluid {sup 4}He present at finite temperatures. The viscous functional is derived from the Navier-Stokes equation by using the Madelung transformation and includes the contribution of interfacial viscous response present at the gas-liquid boundaries. This contribution was obtained by calibrating the model against the experimentally determined electron mobilities from 1.2 K to 2.1 K along the saturated vapor pressure line, where the viscous response is dominated by thermal rotons. The temperature dependence of ion mobility was calculated for several different solvation cavitymore » sizes and the data are rationalized in the context of roton scattering and Stokes limited mobility models. Results are compared to the experimentally observed “exotic ion” data, which provides estimates for the corresponding bubble sizes in the liquid. Possible sources of such ions are briefly discussed.« less

Modeling of brittle-viscous flow using discrete particles

NASA Astrophysics Data System (ADS)

Thordén Haug, Øystein; Barabasch, Jessica; Virgo, Simon; Souche, Alban; Galland, Olivier; Mair, Karen; Abe, Steffen; Urai, Janos L.

2017-04-01

Many geological processes involve both viscous flow and brittle fractures, e.g. boudinage, folding and magmatic intrusions. Numerical modeling of such viscous-brittle materials poses challenges: one has to account for the discrete fracturing, the continuous viscous flow, the coupling between them, and potential pressure dependence of the flow. The Discrete Element Method (DEM) is a numerical technique, widely used for studying fracture of geomaterials. However, the implementation of viscous fluid flow in discrete element models is not trivial. In this study, we model quasi-viscous fluid flow behavior using Esys-Particle software (Abe et al., 2004). We build on the methodology of Abe and Urai (2012) where a combination of elastic repulsion and dashpot interactions between the discrete particles is implemented. Several benchmarks are presented to illustrate the material properties. Here, we present extensive, systematic material tests to characterize the rheology of quasi-viscous DEM particle packing. We present two tests: a simple shear test and a channel flow test, both in 2D and 3D. In the simple shear tests, simulations were performed in a box, where the upper wall is moved with a constant velocity in the x-direction, causing shear deformation of the particle assemblage. Here, the boundary conditions are periodic on the sides, with constant forces on the upper and lower walls. In the channel flow tests, a piston pushes a sample through a channel by Poisseuille flow. For both setups, we present the resulting stress-strain relationships over a range of material parameters, confining stress and strain rate. Results show power-law dependence between stress and strain rate, with a non-linear dependence on confining force. The material is strain softening under some conditions (which). Additionally, volumetric strain can be dilatant or compactant, depending on porosity, confining pressure and strain rate. Constitutive relations are implemented in a way that limits the range of viscosities. For identical pressure and strain rate, an order of magnitude range in viscosity can be investigated. The extensive material testing indicates that DEM particles interacting by a combination of elastic repulsion and dashpots can be used to model viscous flows. This allows us to exploit the fracturing capabilities of the discrete element methods and study systems that involve both viscous flow and brittle fracturing. However, the small viscosity range achievable using this approach does constraint the applicability for systems where larger viscosity ranges are required, such as folding of viscous layers of contrasting viscosities. References: Abe, S., Place, D., & Mora, P. (2004). A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics. PAGEOPH, 161(11-12), 2265-2277. http://doi.org/10.1007/s00024-004-2562-x Abe, S., and J. L. Urai (2012), Discrete element modeling of boudinage: Insights on rock rheology, matrix flow, and evolution of geometry, JGR., 117, B01407, doi:10.1029/2011JB00855

Locomotion in a liquid crystal near a wall

NASA Astrophysics Data System (ADS)

Powers, Thomas; Krieger, Madison; Spagnolie, Saverio

2015-11-01

Recent observations of bacteria swimming in nematic liquid crystal solution motivate the theoretical study of how swimming speed depends on liquid crystal properties. We consider the Taylor sheet near a wall, in which propulsion is achieved by the propagation of traveling waves along the length of the swimmer. Using the lubrication approximation, we determine how swimming speed depends on the Ericksen number, which is the ratio of elastic to viscous stresses. We also study the effect of anchoring strength, at the surface of the swimmer and the surface of the wall. Supported by NSF-CBET 1437195.

Can Satellite Geodesy Disentangle Holocene Rebound and Present-Day Glacier Balance Signatures?

NASA Technical Reports Server (NTRS)

Irvins, E.; James, T.; Yoder, C.

1995-01-01

The secular drift of the precession of the ascending node of the LAGOES -1 satellite is apparently linked to the Earth s paleoclimate through the slow viscous response of the mantle to ice sheet/ocean mass transfer during the last great continental deglaciation . The secular node acceleration is particularly sensitive to the longest wavelengths of the paleo -surface loading that have been memorized by the mantle glacio -isostatic flow. Tide gauge records for the last 130 years show a post-glacial rebound-corrected sea-level rise of 2.4 n 0.9 mm yr-1.

Three-dimensional real-time imaging of bi-phasic flow through porous media

NASA Astrophysics Data System (ADS)

Sharma, Prerna; Aswathi, P.; Sane, Anit; Ghosh, Shankar; Bhattacharya, S.

2011-11-01

We present a scanning laser-sheet video imaging technique to image bi-phasic flow in three-dimensional porous media in real time with pore-scale spatial resolution, i.e., 35 μm and 500 μm for directions parallel and perpendicular to the flow, respectively. The technique is illustrated for the case of viscous fingering. Using suitable image processing protocols, both the morphology and the movement of the two-fluid interface, were quantitatively estimated. Furthermore, a macroscopic parameter such as the displacement efficiency obtained from a microscopic (pore-scale) analysis demonstrates the versatility and usefulness of the method.

Application of a Third Order Upwind Scheme to Viscous Flow over Clean and Iced Wings

NASA Technical Reports Server (NTRS)

Bangalore, A.; Phaengsook, N.; Sankar, L. N.

1994-01-01

A 3-D compressible Navier-Stokes solver has been developed and applied to 3-D viscous flow over clean and iced wings. This method uses a third order accurate finite volume scheme with flux difference splitting to model the inviscid fluxes, and second order accurate symmetric differences to model the viscous terms. The effects of turbulence are modeled using a Kappa-epsilon model. In the vicinity of the sold walls the kappa and epsilon values are modeled using Gorski's algebraic model. Sampling results are presented for surface pressure distributions, for untapered swept clean and iced wings made of NACA 0012 airfoil sections. The leading edge of these sections is modified using a simulated ice shape. Comparisons with experimental data are given.

The controlling effect of viscous dissipation on magma flow in silicic conduits

USGS Publications Warehouse

Mastin, L.G.

2005-01-01

Nearly all volcanic conduit models assume that flow is Newtonian and isothermal. Such models predict that, during high-flux silicic eruptions, gradients in pressure with depth increase upward as magma accelerates and becomes more viscous, leading to extremely low pressure and fragmentation at a depth of kilometers below the surface. In this paper I show that shear heating, also known as viscous dissipation, dramatically reduces the pressure gradient required for flow and concentrates shear in narrow zones along the conduit margin. The reduction in friction may eliminate the zone of low pressure predicted by isothermal models and move the fragmentation level up to the surface.

Lifecycle of miscible viscous fingering: onset to shutdown

NASA Astrophysics Data System (ADS)

Nijjer, Japinder S.; Hewitt, Duncan R.; Neufeld, Jerome A.

2017-11-01

When a viscous fluid is injected into a porous medium or Hele-Shaw cell that is initially saturated with a more viscous fluid, the flow can be unstable to viscous fingering. We investigate the long-time dynamics of miscible viscous fingering in a homogeneous, planar, two-dimensional porous medium using high-resolution numerical simulations. At late times, we identify a new flow regime which consists of a pair of counter-propagating fingers that diffuse and slow, leaving a linearly well-mixed interior. We derive an analytic solution for this regime, and show that, in contrast to previous suggestions, the flow always evolves to this regime irrespective of the viscosity ratio and Peclet number. As a consequence, we find the instability can only ever generate a finite amount of advective mixing. We also describe the full life-cycle of miscible viscous fingering, which can be partitioned into three regimes: an early-time linearly unstable regime, an intermediate-time non-linear regime, and a late-time exchange-flow regime. We identify, using linear stability theory, a critical Peclet number below which the flow is always stable, and derive a model for the evolution of the transversely averaged concentration in the intermediate-time regime, which extends previous empirical models.

Application of viscous and Iwan modal damping models to experimental measurements from bolted structures

DOE PAGES

Deaner, Brandon J.; Allen, Matthew S.; Starr, Michael James; ...

2015-01-20

Measurements are presented from a two-beam structure with several bolted interfaces in order to characterize the nonlinear damping introduced by the joints. The measurements (all at force levels below macroslip) reveal that each underlying mode of the structure is well approximated by a single degree-of-freedom (SDOF) system with a nonlinear mechanical joint. At low enough force levels, the measurements show dissipation that scales as the second power of the applied force, agreeing with theory for a linear viscously damped system. This is attributed to linear viscous behavior of the material and/or damping provided by the support structure. At larger forcemore » levels, the damping is observed to behave nonlinearly, suggesting that damping from the mechanical joints is dominant. A model is presented that captures these effects, consisting of a spring and viscous damping element in parallel with a four-parameter Iwan model. As a result, the parameters of this model are identified for each mode of the structure and comparisons suggest that the model captures the stiffness and damping accurately over a range of forcing levels.« less

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

Zhang, Na; Zhang, Peng; Kang, Wei

Multiscale simulations of fluids such as blood represent a major computational challenge of coupling the disparate spatiotemporal scales between molecular and macroscopic transport phenomena characterizing such complex fluids. In this paper, a coarse-grained (CG) particle model is developed for simulating blood flow by modifying the Morse potential, traditionally used in Molecular Dynamics for modeling vibrating structures. The modified Morse potential is parameterized with effective mass scales for reproducing blood viscous flow properties, including density, pressure, viscosity, compressibility and characteristic flow dynamics of human blood plasma fluid. The parameterization follows a standard inverse-problem approach in which the optimal micro parameters aremore » systematically searched, by gradually decoupling loosely correlated parameter spaces, to match the macro physical quantities of viscous blood flow. The predictions of this particle based multiscale model compare favorably to classic viscous flow solutions such as Counter-Poiseuille and Couette flows. It demonstrates that such coarse grained particle model can be applied to replicate the dynamics of viscous blood flow, with the advantage of bridging the gap between macroscopic flow scales and the cellular scales characterizing blood flow that continuum based models fail to handle adequately.« less

Viscoacoustic model for near-field ultrasonic levitation.

PubMed

Melikhov, Ivan; Chivilikhin, Sergey; Amosov, Alexey; Jeanson, Romain

2016-11-01

Ultrasonic near-field levitation allows for contactless support and transportation of an object over vibrating surface. We developed an accurate model predicting pressure distribution in the gap between the surface and levitating object. The formulation covers a wide range of the air flow regimes: from viscous squeezed flow dominating in small gap to acoustic wave propagation in larger gap. The paper explains derivation of the governing equations from the basic fluid dynamics. The nonreflective boundary conditions were developed to properly define air flow at the outlet. Comparing to direct computational fluid dynamics modeling our approach allows achieving good accuracy while keeping the computation cost low. Using the model we studied the levitation force as a function of gap distance. It was shown that there are three distinguished flow regimes: purely viscous, viscoacoustic, and acoustic. The regimes are defined by the balance of viscous and inertial forces. In the viscous regime the pressure in the gap is close to uniform while in the intermediate viscoacoustic and the acoustic regimes the pressure profile is wavy. The model was validated by a dedicated levitation experiment and compared to similar published results.

Viscoacoustic model for near-field ultrasonic levitation

NASA Astrophysics Data System (ADS)

Melikhov, Ivan; Chivilikhin, Sergey; Amosov, Alexey; Jeanson, Romain

2016-11-01

Ultrasonic near-field levitation allows for contactless support and transportation of an object over vibrating surface. We developed an accurate model predicting pressure distribution in the gap between the surface and levitating object. The formulation covers a wide range of the air flow regimes: from viscous squeezed flow dominating in small gap to acoustic wave propagation in larger gap. The paper explains derivation of the governing equations from the basic fluid dynamics. The nonreflective boundary conditions were developed to properly define air flow at the outlet. Comparing to direct computational fluid dynamics modeling our approach allows achieving good accuracy while keeping the computation cost low. Using the model we studied the levitation force as a function of gap distance. It was shown that there are three distinguished flow regimes: purely viscous, viscoacoustic, and acoustic. The regimes are defined by the balance of viscous and inertial forces. In the viscous regime the pressure in the gap is close to uniform while in the intermediate viscoacoustic and the acoustic regimes the pressure profile is wavy. The model was validated by a dedicated levitation experiment and compared to similar published results.

Clast-fabric development in a shearing granular material: Implications for subglacial till and fault gouge

USGS Publications Warehouse

Hooyer, T.S.; Iverson, N.R.

2000-01-01

Elongate clasts in subglacial till and in fault gouge align during shearing, but the relation between clast-fabric strength and cumulative shear strain for such materials is effectively unknown. This relation was explored in experiments with a large ring-shear device in which a till and a viscous putty that contained isolated clasts were sheared to high strains. As expected, rotation of clasts in the putty is closely approximated by the theory of G.B. Jeffrey, who derived the orbits of rigid ellipsoids in a slowly shearing fluid. Clast rotation in the till, however, is strikingly different. Rather than orbiting through the shear plane as predicted by Jeffery, most clasts rotate into the shear plane and remain there, resulting in strong fabrics regardless of the aspect ratios and initial orientations of clasts. This divergent behavior is likely due to slip of the till matrix along the surfaces of clasts, which is a natural expectation in a granular material but violates the no-slip condition of Jeffery's model. These results do not support the widespread belief that subglacial till deformation results in weak clast fabrics. Thus, many tills with weak fabrics thought to have been sheared subglacially to high strains, like many basal tills of the Laurentide Ice Sheet, may have been sheared only slightly with little effect on either ice-sheet dynamics or sediment transport. In addition, these results indicate that in simple shear the rotation of clasts in till and in fault gouge is best analyzed with the model of A. March, who treated inclusions as passive markers.

Investigation into Behavior of a Steam-Water Mixture Flow Through Holes in a Submerged Perforated Sheet at High Void Fractions

NASA Astrophysics Data System (ADS)

Melikhov, V. I.; Melikhov, O. I.; Nerovnov, A. A.; Nikonov, S. M.

2018-01-01

Processing of experimental data on the pressure difference across a submerged perforated sheet (SPS) revealed that, at sufficiently high void fractions under SPS, the pressure difference across it became less than the pressure difference for the pure steam stream with the same flowrate. To find the cause of this, the effect of a liquid film, which can be formed on the SPS upstream surface as a result of water droplets' impact and can smooth over sharp edges of holes in SDS, was examined. This can decrease the pressure drop across the sharp edges of holes. This assumption was checked through numerical solution to several model problems in the axisymmetric formulation for a steam flow in a round pipe with an orifice. The flow of steam and water was modeled using the viscous incompressible liquid approximation, while the turbulence was described by the k-ɛ model. The evolution of the interfacial area was modeled using the VOF model. The following model problems of steam flow through an orifice were studied: a single-phase flow, a flow through the orifice with a liquid film on its upstream surface, a flow through a chamfered hole, and a flow through the orifice with a liquid film on its upstream surface without liquid supply to the film. The predictions demonstrate that even the approximate account of the liquid film effect on the steam flow yields a considerable decrease in the pressure drop across the hole (from 8 to 24%) due to smoothing its sharp outlet edges over. This makes it possible to make a conclusion that the cause of a decrease in the pressure drop across SPS observed in the experiments at high void fractions is the formation of a liquid film, which smooths the sharp edges of the hole.

New holographic dark energy model with constant bulk viscosity in modified f(R,T) gravity theory

NASA Astrophysics Data System (ADS)

Srivastava, Milan; Singh, C. P.

2018-06-01

The aim of this paper is to study new holographic dark energy (HDE) model in modified f(R,T) gravity theory within the framework of a flat Friedmann-Robertson-Walker model with bulk viscous matter content. It is thought that the negative pressure caused by the bulk viscosity can play the role of dark energy component, and drive the accelerating expansion of the universe. This is the motive of this paper to observe such phenomena with bulk viscosity. In the specific model f(R,T)=R+λ T, where R is the Ricci scalar, T the trace of the energy-momentum tensor and λ is a constant, we find the solution for non-viscous and viscous new HDE models. We analyze new HDE model with constant bulk viscosity, ζ =ζ 0= const. to explain the present accelerated expansion of the universe. We classify all possible scenarios (deceleration, acceleration and their transition) with possible positive and negative ranges of λ over the constraint on ζ 0 to analyze the evolution of the universe. We obtain the solutions of scale factor and deceleration parameter, and discuss the evolution of the universe. We observe the future finite-time singularities of type I and III at a finite time under certain constraints on λ . We also investigate the statefinder and Om diagnostics of the viscous new HDE model to discriminate with other existing dark energy models. In late time the viscous new HDE model approaches to Λ CDM model. We also discuss the thermodynamics and entropy of the model and find that it satisfies the second law of thermodynamics.

Crustal fingering: solidification on a viscously unstable interface

NASA Astrophysics Data System (ADS)

Fu, Xiaojing; Jimenez-Martinez, Joaquin; Cueto-Felgueroso, Luis; Porter, Mark; Juanes, Ruben

2017-11-01

Motivated by the formation of gas hydrates in seafloor sediments, here we study the volumetric expansion of a less viscous gas pocket into a more viscous liquid when the gas-liquid interfaces readily solidify due to hydrate formation. We first present a high-pressure microfluidic experiment to study the depressurization-controlled expansion of a Xenon gas pocket in a water-filled Hele-Shaw cell. The evolution of the pocket is controlled by three processes: (1) volumetric expansion of the gas; (2) rupturing of existing hydrate films on the gas-liquid interface; and (3) formation of new hydrate films. These result in gas fingering leading to a complex labyrinth pattern. To reproduce these observations, we propose a phase-field model that describes the formation of hydrate shell on viscously unstable interfaces. We design the free energy of the three-phase system to rigorously account for interfacial effects, gas compressibility and phase transitions. We model the hydrate shell as a highly viscous fluid with shear-thinning rheology to reproduce shell-rupturing behavior. We present high-resolution numerical simulations of the model, which illustrate the emergence of complex crustal fingering patterns as a result of gas expansion dynamics modulated by hydrate growth at the interface.

Spiral pattern in a radial displacement in a Hele-Shaw cell

NASA Astrophysics Data System (ADS)

Ban, Mitsumasa; Nagatsu, Yuichiro; Hayashi, Atsushi; Kato, Yoshihiro; Tada, Yutaka

2008-11-01

When a reactive and miscible less-viscous liquid displaces a more-viscous liquid in a Hele-Shaw cell, reactive miscible viscous fingering takes place. We have experimentally shown that the pattern created by the displacement of a more-viscous fluid by a less-viscous one in a radial Hele-Shaw cell develops not radially but spirally when a more-viscous sodium polyacrylate solution is displaced by a less-viscous trivalent iron ion (Fe^3+) solution with a sufficiently high concentration of Fe^3+. Another experiment in order to investigate the mechanism of spiral pattern formation revealed that an instantaneous chemical reaction takes place between the two fluids and at high Fe^3+ concentrations it produces a film of the gel at the contact plane. The gel is formed by three-dimensional network structures between the polyacrylate solution and the trivalent iron ion (Fe^3+) solution. We have proposed a physical model that the gel's film is responsible for the form of the spiral pattern.

Implications of Grain Size Evolution for the Effective Stress Exponent in Ice

NASA Astrophysics Data System (ADS)

Behn, M. D.; Goldsby, D. L.; Hirth, G.

2016-12-01

Viscous flow in ice has typically been described by the Glen law—a non-Newtonian, power-law relationship between stress and strain-rate with a stress exponent n 3. The Glen law is attributed to grain-size-insensitive dislocation creep; however, laboratory and field studies demonstrate that deformation in ice is strongly dependent on grain size. This has led to the hypothesis that at sufficiently low stresses, ice flow is controlled by grain boundary sliding [1], which explicitly incorporates the grain-size dependence of ice rheology. Yet, neither dislocation creep (n 4), nor grain boundary sliding (n 1.8), have stress exponents that match the value of n 3 for the Glen law. Thus, although the Glen law provides an approximate description of ice flow in glaciers and ice sheets, its functional form cannot be explained by a single deformation mechanism. Here we seek to understand the origin of the n 3 dependence of the Glen law through a new model for grain-size evolution in ice. In our model, grain size evolves in response to the balance between dynamic recrystallization and grain growth. To simulate these processes we adapt the "wattmeter" [2], originally developed within the solid-Earth community to quantify grain size in crustal and mantle rocks. The wattmeter posits that grain size is controlled by a balance between the mechanical work required for grain growth and dynamic grain size reduction. The evolution of grain size in turn controls the relative contributions of dislocation creep and grain boundary sliding, and thus the effective stress exponent for ice flow. Using this approach, we first benchmark our grain size evolution model on experimental data and then calculate grain size in two end-member scenarios: (1) as a function of depth within an ice-sheet, and (2) across an ice-stream margin. We show that the calculated grain sizes match ice core observations for the interior of ice sheets. Furthermore, owing to the influence of grain size on strain rate, the variation in grain size with deformation conditions results in an effective stress exponent intermediate between grain boundary sliding and dislocation creep. [1] Goldsby & Kohlstedt, JGR, 2001; [2] Austin & Evans, Geology, 1997

Including fluid shear viscosity in a structural acoustic finite element model using a scalar fluid representation

PubMed Central

Cheng, Lei; Li, Yizeng; Grosh, Karl

2013-01-01

An approximate boundary condition is developed in this paper to model fluid shear viscosity at boundaries of coupled fluid-structure system. The effect of shear viscosity is approximated by a correction term to the inviscid boundary condition, written in terms of second order in-plane derivatives of pressure. Both thin and thick viscous boundary layer approximations are formulated; the latter subsumes the former. These approximations are used to develop a variational formation, upon which a viscous finite element method (FEM) model is based, requiring only minor modifications to the boundary integral contributions of an existing inviscid FEM model. Since this FEM formulation has only one degree of freedom for pressure, it holds a great computational advantage over the conventional viscous FEM formulation which requires discretization of the full set of linearized Navier-Stokes equations. The results from thick viscous boundary layer approximation are found to be in good agreement with the prediction from a Navier-Stokes model. When applicable, thin viscous boundary layer approximation also gives accurate results with computational simplicity compared to the thick boundary layer formulation. Direct comparison of simulation results using the boundary layer approximations and a full, linearized Navier-Stokes model are made and used to evaluate the accuracy of the approximate technique. Guidelines are given for the parameter ranges over which the accurate application of the thick and thin boundary approximations can be used for a fluid-structure interaction problem. PMID:23729844

Including fluid shear viscosity in a structural acoustic finite element model using a scalar fluid representation.

PubMed

Cheng, Lei; Li, Yizeng; Grosh, Karl

2013-08-15

An approximate boundary condition is developed in this paper to model fluid shear viscosity at boundaries of coupled fluid-structure system. The effect of shear viscosity is approximated by a correction term to the inviscid boundary condition, written in terms of second order in-plane derivatives of pressure. Both thin and thick viscous boundary layer approximations are formulated; the latter subsumes the former. These approximations are used to develop a variational formation, upon which a viscous finite element method (FEM) model is based, requiring only minor modifications to the boundary integral contributions of an existing inviscid FEM model. Since this FEM formulation has only one degree of freedom for pressure, it holds a great computational advantage over the conventional viscous FEM formulation which requires discretization of the full set of linearized Navier-Stokes equations. The results from thick viscous boundary layer approximation are found to be in good agreement with the prediction from a Navier-Stokes model. When applicable, thin viscous boundary layer approximation also gives accurate results with computational simplicity compared to the thick boundary layer formulation. Direct comparison of simulation results using the boundary layer approximations and a full, linearized Navier-Stokes model are made and used to evaluate the accuracy of the approximate technique. Guidelines are given for the parameter ranges over which the accurate application of the thick and thin boundary approximations can be used for a fluid-structure interaction problem.

PDF modeling of near-wall turbulent flows

NASA Astrophysics Data System (ADS)

Dreeben, Thomas David

1997-06-01

Pdf methods are extended to include modeling of wall- bounded turbulent flows. For flows in which resolution of the viscous sublayer is desired, a Pdf near-wall model is developed in which the Generalized Langevin model is combined with an exact model for viscous transport. Durbin's method of elliptic relaxation is used to incorporate the wall effects into the governing equations without the use of wall functions or damping functions. Close to the wall, the Generalized Langevin model provides an analogy to the effect of the fluctuating continuity equation. This enables accurate modeling of the near-wall turbulent statistics. Demonstrated accuracy for fully-developed channel flow is achieved with a Pdf/Monte Carlo simulation, and with its related Reynolds-stress closure. For flows in which the details of the viscous sublayer are not important, a Pdf wall- function method is developed with the Simplified Langevin model.

A Hermite-based lattice Boltzmann model with artificial viscosity for compressible viscous flows

NASA Astrophysics Data System (ADS)

Qiu, Ruofan; Chen, Rongqian; Zhu, Chenxiang; You, Yancheng

2018-05-01

A lattice Boltzmann model on Hermite basis for compressible viscous flows is presented in this paper. The model is developed in the framework of double-distribution-function approach, which has adjustable specific-heat ratio and Prandtl number. It contains a density distribution function for the flow field and a total energy distribution function for the temperature field. The equilibrium distribution function is determined by Hermite expansion, and the D3Q27 and D3Q39 three-dimensional (3D) discrete velocity models are used, in which the discrete velocity model can be replaced easily. Moreover, an artificial viscosity is introduced to enhance the model for capturing shock waves. The model is tested through several cases of compressible flows, including 3D supersonic viscous flows with boundary layer. The effect of artificial viscosity is estimated. Besides, D3Q27 and D3Q39 models are further compared in the present platform.

Collisional damping of the geodesic acoustic mode with toroidal rotation. I. Viscous damping

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

Gong, Xueyu; Xie, Baoyi; Chen, You

2016-03-15

With the dispersion relation derived for the geodesic acoustic mode in toroidally rotating tokamak plasmas using the fluid model, the effect of the toroidal rotation on the collisional viscous damping of the geodesic acoustic mode is investigated. It is found that the collisional viscous damping of the geodesic acoustic mode has weak increase with respect to the toroidal Mach number.

Three mechanisms model of shale gas in real state transport through a single nanopore

NASA Astrophysics Data System (ADS)

Li, Dongdong; Zhang, Yanyu; Sun, Xiaofei; Li, Peng; Zhao, Fengkai

2018-02-01

At present, the apparent permeability models of shale gas consider only the viscous flow and Knudsen diffusion of free gas, but do not take into account the influence of surface diffusion. Moreover, it is assumed that shale gas is in ideal state. In this paper, shale gas is assumed in real state, a new apparent permeability model for shale gas transport through a single nanopore is developed that captures many important migration mechanisms, such as viscous flow and Knudsen diffusion of free gas, surface diffusion of adsorbed gas. According to experimental data, the accuracy of apparent permeability model was verified. What’s more, the effects of pressure and pore radius on apparent permeability, and the effects on the permeability fraction of viscous flow, Knudsen diffusion and surface diffusion were analysed, separately. Finally, the results indicate that the error of the developed model in this paper was 3.02%, which is less than the existing models. Pressure and pore radius seriously affect the apparent permeability of shale gas. When the pore radius is small or pressure is low, the surface diffusion cannot be ignored. When the pressure and the pore radius is big, the viscous flow occupies the main position.

Inviscid Wall-Modeled Large Eddy Simulations for Improved Efficiency

NASA Astrophysics Data System (ADS)

Aikens, Kurt; Craft, Kyle; Redman, Andrew

2015-11-01

The accuracy of an inviscid flow assumption for wall-modeled large eddy simulations (LES) is examined because of its ability to reduce simulation costs. This assumption is not generally applicable for wall-bounded flows due to the high velocity gradients found near walls. In wall-modeled LES, however, neither the viscous near-wall region or the viscous length scales in the outer flow are resolved. Therefore, the viscous terms in the Navier-Stokes equations have little impact on the resolved flowfield. Zero pressure gradient flat plate boundary layer results are presented for both viscous and inviscid simulations using a wall model developed previously. The results are very similar and compare favorably to those from another wall model methodology and experimental data. Furthermore, the inviscid assumption reduces simulation costs by about 25% and 39% for supersonic and subsonic flows, respectively. Future research directions are discussed as are preliminary efforts to extend the wall model to include the effects of unresolved wall roughness. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1053575. Computational resources on TACC Stampede were provided under XSEDE allocation ENG150001.

Parameterizing the Morse Potential for Coarse-Grained Modeling of Blood Plasma

PubMed Central

Zhang, Na; Zhang, Peng; Kang, Wei; Bluestein, Danny; Deng, Yuefan

2014-01-01

Multiscale simulations of fluids such as blood represent a major computational challenge of coupling the disparate spatiotemporal scales between molecular and macroscopic transport phenomena characterizing such complex fluids. In this paper, a coarse-grained (CG) particle model is developed for simulating blood flow by modifying the Morse potential, traditionally used in Molecular Dynamics for modeling vibrating structures. The modified Morse potential is parameterized with effective mass scales for reproducing blood viscous flow properties, including density, pressure, viscosity, compressibility and characteristic flow dynamics of human blood plasma fluid. The parameterization follows a standard inverse-problem approach in which the optimal micro parameters are systematically searched, by gradually decoupling loosely correlated parameter spaces, to match the macro physical quantities of viscous blood flow. The predictions of this particle based multiscale model compare favorably to classic viscous flow solutions such as Counter-Poiseuille and Couette flows. It demonstrates that such coarse grained particle model can be applied to replicate the dynamics of viscous blood flow, with the advantage of bridging the gap between macroscopic flow scales and the cellular scales characterizing blood flow that continuum based models fail to handle adequately. PMID:24910470

Contributions of Kinetic Energy and Viscous Dissipation to Airway Resistance in Pulmonary Inspiratory and Expiratory Airflows in Successive Symmetric Airway Models With Various Bifurcation Angles.

PubMed

Choi, Sanghun; Choi, Jiwoong; Lin, Ching-Long

2018-01-01

The aim of this study was to investigate and quantify contributions of kinetic energy and viscous dissipation to airway resistance during inspiration and expiration at various flow rates in airway models of different bifurcation angles. We employed symmetric airway models up to the 20th generation with the following five different bifurcation angles at a tracheal flow rate of 20 L/min: 15 deg, 25 deg, 35 deg, 45 deg, and 55 deg. Thus, a total of ten computational fluid dynamics (CFD) simulations for both inspiration and expiration were conducted. Furthermore, we performed additional four simulations with tracheal flow rate values of 10 and 40 L/min for a bifurcation angle of 35 deg to study the effect of flow rate on inspiration and expiration. Using an energy balance equation, we quantified contributions of the pressure drop associated with kinetic energy and viscous dissipation. Kinetic energy was found to be a key variable that explained the differences in airway resistance on inspiration and expiration. The total pressure drop and airway resistance were larger during expiration than inspiration, whereas wall shear stress and viscous dissipation were larger during inspiration than expiration. The dimensional analysis demonstrated that the coefficients of kinetic energy and viscous dissipation were strongly correlated with generation number. In addition, the viscous dissipation coefficient was significantly correlated with bifurcation angle and tracheal flow rate. We performed multiple linear regressions to determine the coefficients of kinetic energy and viscous dissipation, which could be utilized to better estimate the pressure drop in broader ranges of successive bifurcation structures.

Gravitational waves from remnant massive neutron stars of binary neutron star merger: Viscous hydrodynamics effects

NASA Astrophysics Data System (ADS)

Shibata, Masaru; Kiuchi, Kenta

2017-06-01

Employing a simplified version of the Israel-Stewart formalism of general-relativistic shear-viscous hydrodynamics, we explore the evolution of a remnant massive neutron star of binary neutron star merger and pay special attention to the resulting gravitational waveforms. We find that for the plausible values of the so-called viscous alpha parameter of the order 10-2 the degree of the differential rotation in the remnant massive neutron star is significantly reduced in the viscous time scale, ≲5 ms . Associated with this, the degree of nonaxisymmetric deformation is also reduced quickly, and as a consequence, the amplitude of quasiperiodic gravitational waves emitted also decays in the viscous time scale. Our results indicate that for modeling the evolution of the merger remnants of binary neutron stars we would have to take into account magnetohydrodynamics effects, which in nature could provide the viscous effects.

Bianchi Type VI1 Viscous Fluid Cosmological Model in Wesson´s Theory of Gravitation

NASA Astrophysics Data System (ADS)

Khadekar, G. S.; Avachar, G. R.

2007-03-01

Field equations of a scale invariant theory of gravitation proposed by Wesson [1, 2] are obtained in the presence of viscous fluid with the aid of Bianchi type VIh space-time with the time dependent gauge function (Dirac gauge). It is found that Bianchi type VIh (h = 1) space-time with viscous fluid is feasible in this theory, whereas Bianchi type VIh (h = -1, 0) space-times are not feasible in this theory, even in the presence of viscosity. For the feasible case, by assuming a relation connecting viscosity and metric coefficient, we have obtained a nonsingular-radiating model. We have discussed some physical and kinematical properties of the models.

Ion and electron heating characteristics of magnetic reconnection in tokamak plasma merging experiments

NASA Astrophysics Data System (ADS)

Ono, Y.; Tanabe, H.; Yamada, T.; Inomoto, M.; T, Ii; Inoue, S.; Gi, K.; Watanabe, T.; Gryaznevich, M.; Scannell, R.; Michael, C.; Cheng, C. Z.

2012-12-01

Recently, the TS-3 and TS-4 tokamak merging experiments revealed significant plasma heating during magnetic reconnection. A key question is how and where ions and electrons are heated during magnetic reconnection. Two-dimensional measurements of ion and electron temperatures and plasma flow made clear that electrons are heated inside the current sheet mainly by the Ohmic heating and ions are heated in the downstream areas mainly by the reconnection outflows. The outflow kinetic energy is thermalized by the fast shock formation and viscous damping. The magnetic reconnection converts the reconnecting magnetic field energy mostly to the ion thermal energy in the outflow region whose size is much larger than the current sheet size for electron heating. The ion heating energy is proportional to the square of the reconnection magnetic field component B_p^2 . This scaling of reconnection heating indicates the significant ion heating effect of magnetic reconnection, which leads to a new high-field reconnection heating experiment for fusion plasmas.

Rippling Instability of a Collapsing Bubble

NASA Technical Reports Server (NTRS)

daSilveira, Rava; Chaieb, Sahraoui; Mahadevan, L.

1999-01-01

The rippling instability of a liquid sheet was first observed by Debregeas, de Gennes, an Brochard-Wyart [Science 279, 1704 (1998)] on a hemispherical bubble resting on a free surface. Unlike a soap bubble, it collapses under its own weight while bursting, and folds into a wavy structure which breaks the original axisymmetry. In fact, this effect occurs for both purely elastic and purely viscous (liquid) sheets, and an analogy can be made between the two mechanisms. We present a theory for the onset of the instability in both cases, in which the growth of the corrugation out of an inextensible initial condition is governed by the competition between gravitational and bending (shearing) forces. The instability occurs for a range of densities, stiffnesses (viscosities), and sizes, a result which arises less from dynamics than from geometry, suggesting a wide validity. We further obtain a quantitative expression for the number of ripples. Finally, we present the results of experiments, which are consistent with our predictions.

The role of zonal flows in disc gravito-turbulence

NASA Astrophysics Data System (ADS)

Vanon, R.

2018-07-01

The work presented here focuses on the role of zonal flows in the self-sustenance of gravito-turbulence in accretion discs. The numerical analysis is conducted using a bespoke pseudo-spectral code in fully compressible, non-linear conditions. The disc in question, which is modelled using the shearing sheet approximation, is assumed to be self-gravitating, viscous, and thermally diffusive; a constant cooling time-scale is also considered. Zonal flows are found to emerge at the onset of gravito-turbulence and they remain closely linked to the turbulent state. A cycle of zonal flow formation and destruction is established, mediated by a slow mode instability (which allows zonal flows to grow) and a non-axisymmetric instability (which disrupts the zonal flow), which is found to repeat numerous times. It is in fact the disruptive action of the non-axisymmetric instability to form new leading and trailing shearing waves, allowing energy to be extracted from the background flow and ensuring the self-sustenance of the gravito-turbulent regime.

The role of zonal flows in disc gravito-turbulence

NASA Astrophysics Data System (ADS)

Vanon, R.

2018-04-01

The work presented here focuses on the role of zonal flows in the self-sustenance of gravito-turbulence in accretion discs. The numerical analysis is conducted using a bespoke pseudo-spectral code in fully compressible, non-linear conditions. The disc in question, which is modelled using the shearing sheet approximation, is assumed to be self-gravitating, viscous, and thermally diffusive; a constant cooling timescale is also considered. Zonal flows are found to emerge at the onset of gravito-turbulence and they remain closely linked to the turbulent state. A cycle of zonal flow formation and destruction is established, mediated by a slow mode instability (which allows zonal flows to grow) and a non-axisymmetric instability (which disrupts the zonal flow), which is found to repeat numerous times. It is in fact the disruptive action of the non-axisymmetric instability to form new leading and trailing shearing waves, allowing energy to be extracted from the background flow and ensuring the self-sustenance of the gravito-turbulent regime.

Formation of lobate debris aprons on Mars: Assessment of regional ice sheet collapse and debris-cover armoring

NASA Astrophysics Data System (ADS)

Fastook, James L.; Head, James W.; Marchant, David R.

2014-01-01

Lobate debris aprons (LDA) are lobate-shaped aprons surrounding scarps and isolated massifs that are concentrated in the vicinity of the northern Dichotomy Boundary on Mars. LDAs have been interpreted as (1) ice-cemented talus aprons undergoing viscous flow, (2) local debris-covered alpine-like glaciers, or (3) remnants of the collapse of a regional retreating ice sheet. We investigate the plausibility that LDAs are remnants of a more extensive regional ice sheet by modeling this process. We find that as a regional ice sheet collapses, the surface drops below cliff and massif bedrock margins, exposing bedrock and regolith, and initiating debris deposition on the surface of a cold-based glacier. Reduced sublimation due to debris-cover armoring of the proto-LDA surface produces a surface slope and consequent ice flow that carries the armoring debris away from the rock outcrops. As collapse and ice retreat continue the debris train eventually reaches the substrate surface at the front of the glacier, leaving the entire LDA armored by debris cover. Using a simplified ice flow model we are able to characterize the temperature and sublimation rate that would be necessary to produce LDAs with a wide range of specified lateral extents and thicknesses. We then apply this method to a database of documented LDA parameters (height, lateral extent) from the Dichotomy Boundary region, and assess the implications for predicted climate conditions during their formation and the range of formation times implied by the model. We find that for the population examined here, typical temperatures are in the range of -85 to -40 °C and typical sublimation rates lie in the range of 6-14 mm/a. Lobate debris apron formation times (from the point of bedrock exposure to complete debris cover) cluster near 400-500 ka. These results show that LDA length and thickness characteristics are consistent with climate conditions and a formation scenario typical of the collapse of a regional retreating ice sheet and exposure of bedrock cliffs. This scenario helps resolve many of the unusual characteristics of lobate debris aprons (LDA) and lineated valley fill (LVF). For example, the distribution of LVF is very consistent with extensive flow of glacial ice from plateau icefields, and the acquisition of a debris cover in the waning stages of retreat of the regional cover as the bedrock scarps are exposed. The typical concentric development of LDA around massifs is much more consistent with ice sheet retreat than insolation-related local accumulation and flow. We thus conclude that the retreating ice-sheet model is robust and should be investigated and tested in more detail. In addition, these results clearly show that the lobate debris aprons in the vicinity of the Dichotomy Boundary could not have attained temperatures near or above the ice melting point and retained their current shape, a finding that supports subzero temperatures for the last several hundred million years, the age of the LDA surfaces. A further implication is that the LDA ice has been preserved for at least several hundred million years, and could potentially contain the record of the climate of Mars, preserved since that time below a sublimation lag deposit.

One-dimensional reduction of viscous jets. I. Theory

NASA Astrophysics Data System (ADS)

Pitrou, Cyril

2018-04-01

We build a general formalism to describe thin viscous jets as one-dimensional objects with an internal structure. We present in full generality the steps needed to describe the viscous jets around their central line, and we argue that the Taylor expansion of all fields around that line is conveniently expressed in terms of symmetric trace-free tensors living in the two dimensions of the fiber sections. We recover the standard results of axisymmetric jets and we report the first and second corrections to the lowest order description, also allowing for a rotational component around the axis of symmetry. When applied to generally curved fibers, the lowest order description corresponds to a viscous string model whose sections are circular. However, when including the first corrections, we find that curved jets generically develop elliptic sections. Several subtle effects imply that the first corrections cannot be described by a rod model since it amounts to selectively discard some corrections. However, in a fast rotating frame, we find that the dominant effects induced by inertial and Coriolis forces should be correctly described by rod models. For completeness, we also recover the constitutive relations for forces and torques in rod models and exhibit a missing term in the lowest order expression of viscous torque. Given that our method is based on tensors, the complexity of all computations has been beaten down by using an appropriate tensor algebra package such as xAct, allowing us to obtain a one-dimensional description of curved viscous jets with all the first order corrections consistently included. Finally, we find a description for straight fibers with elliptic sections as a special case of these results, and recover that ellipticity is dynamically damped by surface tension. An application to toroidal viscous fibers is presented in the companion paper [Pitrou, Phys. Rev. E 97, 043116 (2018), 10.1103/PhysRevE.97.043116].

Implicit marching solution of compressible viscous subsonic flow in planar and axisymmetric ducts. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Towne, C. E.; Hoffman, J. D.

1982-01-01

A new streamwise marching procedure was developed and coded for compressible viscous subsonic flow in planar or axisymmetric ducts with or without centerbodies. The continuity, streamwise momentum, cross-flow momentum, and energy equations are written in generalized orthogonal curvilinear coordinates. To allow the use of a marching procedure, second derivatives in the streamwise momentum equation are written as the sum of a known two dimensional imposed pressure field and an unknown one dimensional viscous correction. For turbulent flow, the Reynolds stress and heat flux terms are modeled using two-layer eddy viscosity turbulence models.

Topography associated with crustal flow in continental collisions, with application to Tibet

NASA Astrophysics Data System (ADS)

Bendick, R.; McKenzie, D.; Etienne, J.

2008-10-01

Collision between an undeformable indenter and a viscous region generates isostatically compensated topography by solid-state flow. We model this process numerically, using a finite element scheme. The slope, amplitude and symmetry of the topographic signal depend on the indenter size and the Argand number of the viscous region, a dimensionless ratio of gravitational body forces to viscous forces. When applied to convergent continental settings, these scaling rules provide estimates of the position of an indenter at depth and the mechanical properties of the viscous region, especially effective viscosity. In Tibet, forward modelling suggests that some elevated, low relief topography within the northern plateau may be attributed to lower crustal flow, stimulated by a crustal indenter, possibly Indian lithosphere. The best-fit model constrains the northernmost limit of this indenter to 33.7°N and the maximum effective viscosity of Eurasian middle and lower crust to 1 × 1020 +/- 0.3 × 1020 Pa s.

Wettability and Flow Rate Impacts on Immiscible Displacement: A Theoretical Model

NASA Astrophysics Data System (ADS)

Hu, Ran; Wan, Jiamin; Yang, Zhibing; Chen, Yi-Feng; Tokunaga, Tetsu

2018-04-01

When a more viscous fluid displaces a less viscous one in porous media, viscous pressure drop stabilizes the displacement front against capillary pressure fluctuation. For this favorable viscous ratio conditions, previous studies focused on the front instability under slow flow conditions but did not address competing effects of wettability and flow rate. Here we study how this competition controls displacement patterns. We propose a theoretical model that describes the crossover from fingering to stable flow as a function of invading fluid contact angle θ and capillary number Ca. The phase diagram predicted by the model shows that decreasing θ stabilizes the displacement for θ≥45° and the critical contact angle θc increases with Ca. The boundary between corner flow and cooperative filling for θ < 45° is also described. This work extends the classic phase diagram and has potential applications in predicting CO2 capillary trapping and manipulating wettability to enhance gas/oil displacement efficiency.

Exponential Stellar Disks in Low Surface Brightness Galaxies: A Critical Test of Viscous Evolution

NASA Astrophysics Data System (ADS)

Bell, Eric F.

2002-12-01

Viscous redistribution of mass in Milky Way-type galactic disks is an appealing way of generating an exponential stellar profile over many scale lengths, almost independent of initial conditions, requiring only that the viscous timescale and star formation timescale are approximately equal. However, galaxies with solid-body rotation curves cannot undergo viscous evolution. Low surface brightness (LSB) galaxies have exponential surface brightness profiles, yet have slowly rising, nearly solid-body rotation curves. Because of this, viscous evolution may be inefficient in LSB galaxies: the exponential profiles, instead, would give important insight into initial conditions for galaxy disk formation. Using star formation laws from the literature and tuning the efficiency of viscous processes to reproduce an exponential stellar profile in Milky Way-type galaxies, I test the role of viscous evolution in LSB galaxies. Under the conservative and not unreasonable condition that LSB galaxies are gravitationally unstable for at least a part of their lives, I find that it is impossible to rule out a significant role for viscous evolution. This type of model still offers an attractive way of producing exponential disks, even in LSB galaxies with slowly rising rotation curves.

A Weakly Nonlinear Model for the Damping of Resonantly Forced Density Waves in Dense Planetary Rings

NASA Astrophysics Data System (ADS)

Lehmann, Marius; Schmidt, Jürgen; Salo, Heikki

2016-10-01

In this paper, we address the stability of resonantly forced density waves in dense planetary rings. Goldreich & Tremaine have already argued that density waves might be unstable, depending on the relationship between the ring’s viscosity and the surface mass density. In the recent paper Schmidt et al., we have pointed out that when—within a fluid description of the ring dynamics—the criterion for viscous overstability is satisfied, forced spiral density waves become unstable as well. In this case, linear theory fails to describe the damping, but nonlinearity of the underlying equations guarantees a finite amplitude and eventually a damping of the wave. We apply the multiple scale formalism to derive a weakly nonlinear damping relation from a hydrodynamical model. This relation describes the resonant excitation and nonlinear viscous damping of spiral density waves in a vertically integrated fluid disk with density dependent transport coefficients. The model consistently predicts density waves to be (linearly) unstable in a ring region where the conditions for viscous overstability are met. Sufficiently far away from the Lindblad resonance, the surface mass density perturbation is predicted to saturate to a constant value due to nonlinear viscous damping. The wave’s damping lengths of the model depend on certain input parameters, such as the distance to the threshold for viscous overstability in parameter space and the ground state surface mass density.

Modeling and Simulation of Viscous Electro-Active Polymers

PubMed Central

Vogel, Franziska; Göktepe, Serdar; Steinmann, Paul; Kuhl, Ellen

2014-01-01

Electro-active materials are capable of undergoing large deformation when stimulated by an electric field. They can be divided into electronic and ionic electro-active polymers (EAPs) depending on their actuation mechanism based on their composition. We consider electronic EAPs, for which attractive Coulomb forces or local re-orientation of polar groups cause a bulk deformation. Many of these materials exhibit pronounced visco-elastic behavior. Here we show the development and implementation of a constitutive model, which captures the influence of the electric field on the visco-elastic response within a geometrically non-linear finite element framework. The electric field affects not only the equilibrium part of the strain energy function, but also the viscous part. To adopt the familiar additive split of the strain from the small strain setting, we formulate the governing equations in the logarithmic strain space and additively decompose the logarithmic strain into elastic and viscous parts. We show that the incorporation of the electric field in the viscous response significantly alters the relaxation and hysteresis behavior of the model. Our parametric study demonstrates that the model is sensitive to the choice of the electro-viscous coupling parameters. We simulate several actuator structures to illustrate the performance of the method in typical relaxation and creep scenarios. Our model could serve as a design tool for micro-electro-mechanical systems, microfluidic devices, and stimuli-responsive gels such as artificial skin, tactile displays, or artificial muscle. PMID:25267881

Computational Modeling of Shape Memory Polymer Origami that Responds to Light

NASA Astrophysics Data System (ADS)

Mailen, Russell William

Shape memory polymers (SMPs) transform in response to external stimuli, such as infrared (IR) light. Although SMPs have many applications, this investigation focuses on their use as actuators in self-folding origami structures. Ink patterned on the surface of the SMP sheet absorbs thermal energy from the IR light, which produces localized heating. The material shrinks wherever the activation temperature is exceeded and can produce out-of-plane deformation. The time and temperature dependent response of these SMPs provides unique opportunities for developing complex three-dimensional (3D) structures from initially flat sheets through self-folding origami, but the application of this technique requires predicting accurately the final folded or deformed shape. Furthermore, current computational approaches for SMPs do not fully couple the thermo-mechanical response of the material. Hence, a proposed nonlinear, 3D, thermo-viscoelastic finite element framework was formulated to predict deformed shapes for different self-folding systems and compared to experimental results for self-folding origami structures. A detailed understanding of the shape memory response and the effect of controllable design parameters, such as the ink pattern, pre-strain conditions, and applied thermal and mechanical fields, allows for a predictive understanding and design of functional, 3D structures. The proposed modeling framework was used to obtain a fundamental understanding of the thermo-mechanical behavior of SMPs and the impact of the material behavior on hinged self-folding. These predictions indicated how the thermal and mechanical conditions during pre-strain significantly affect the shrinking and folding response of the SMP. Additionally, the externally applied thermal loads significantly influenced the folding rate and maximum bending angle. The computational framework was also adapted to understand the effects of fully coupling the thermal and mechanical response of the material. This updated framework accounted for external heat sources, such as ambient temperature and incident surface heat flux, as well as internal temperature changes due to conduction and viscous heat generation. Viscous heating during the pre-strain sequence affected the residual stresses after cooling due to accelerated viscoelastic relaxation. This resulted in a delayed shrinking and folding response. Other factors that affected the folding response include sheet thickness, hinge width, degree of pre-strain, and hinge temperature. The predicted results indicated that the maximum bending angle can be increased for a folded structure by increasing the hinge width, degree of pre-strain, and hinge surface temperature. Folding time can be reduced by decreasing the sheet thickness, increasing the hinge width, and increasing the hinge temperature. The coupled thermo-mechanical approach was also extended to investigate both curved and folded structures by varying the ink pattern and the substrate geometry. With this approach, two continuous curvature mechanisms were obtained. One was an indirect curvature mechanism which resulted from internal stresses that evolved from the shrinking of activated regions of the material relative to unactivated regions. The second was a direct curvature mechanism that resulted from ink distributed in gradients across the surface of the material. Furthermore, the effects of hinge orientation, proximity of multiple hinges, sheet aspect ratio, and axisymmetric ink patterns were characterized for other shapes, such as rectangles and discs. The findings of this investigation clearly indicate that this validated computational approach can be used to predict and understand the myriad mechanisms of self-folding origami structures. By varying the location of ink on the polymer surface and making changes to the substrate geometry, complex 3D structures can be obtained. The developed thermo-mechanical framework can be used to design optimized origami structures for biomedical devices, space telescopes, and functional, engineered origami devices.

LRS Bianchi Type-I Bulk Viscous Cosmological Models in f( R, T) Gravity

NASA Astrophysics Data System (ADS)

Sahoo, P.; Reddy, R.

2018-03-01

We have studied the locally rotationally symmetric (LRS) Bianchi type-I cosmological model in f ( R, T) gravity (R is the Ricci scalar and T is the trace of the stress energy tensor) with bulk viscous fluid as matter content. The model is constructed for the linear form f ( R, T) = R + 2 f ( T). The exact solution of the field equations is obtained by using a time varying deceleration parameter q for a suitable choice of the function f ( T). In this case, the bulk viscous pressure \\overline{p} is found to be negative and the energy density ρ is found to be positive. The obtained model is anisotropic, accelerating, and compatible with the results of astronomical observations. Also, some important features of physical parameters of this model have been discussed.

Vibrational modes of hydraulic fractures: Inference of fracture geometry from resonant frequencies and attenuation

NASA Astrophysics Data System (ADS)

Lipovsky, Bradley P.; Dunham, Eric M.

2015-02-01

Oscillatory seismic signals arising from resonant vibrations of hydraulic fractures are observed in many geologic systems, including volcanoes, glaciers and ice sheets, and hydrocarbon and geothermal reservoirs. To better quantify the physical dimensions of fluid-filled cracks and properties of the fluids within them, we study wave motion along a thin hydraulic fracture waveguide. We present a linearized analysis, valid at wavelengths greater than the fracture aperture, that accounts for quasi-static elastic deformation of the fracture walls, as well as fluid viscosity, inertia, and compressibility. In the long-wavelength limit, anomalously dispersed guided waves known as crack or Krauklis waves propagate with restoring force from fracture wall elasticity. At shorter wavelengths, the waves become sound waves within the fluid channel. Wave attenuation in our model is due to fluid viscosity, rather than seismic radiation from crack tips or fracture wall roughness. We characterize viscous damping at both low frequencies, where the flow is always fully developed, and at high frequencies, where the flow has a nearly constant velocity profile away from viscous boundary layers near the fracture walls. Most observable seismic signals from resonating fractures likely arise in the boundary layer crack wave limit, where fluid-solid coupling is pronounced and attenuation is minimal. We present a method to estimate the aperture and length of a resonating hydraulic fracture using both the seismically observed quality factor and characteristic frequency. Finally, we develop scaling relations between seismic moment and characteristic frequency that might be useful when interpreting the statistics of hydraulic fracture events.

Shock-induced bubble jetting into a viscous fluid with application to tissue injury in shock-wave lithotripsy.

PubMed

Freund, J B; Shukla, R K; Evan, A P

2009-11-01

Shock waves in liquids are known to cause spherical gas bubbles to rapidly collapse and form strong re-entrant jets in the direction of the propagating shock. The interaction of these jets with an adjacent viscous liquid is investigated using finite-volume simulation methods. This configuration serves as a model for tissue injury during shock-wave lithotripsy, a medical procedure to remove kidney stones. In this case, the viscous fluid provides a crude model for the tissue. It is found that for viscosities comparable to what might be expected in tissue, the jet that forms upon collapse of a small bubble fails to penetrate deeply into the viscous fluid "tissue." A simple model reproduces the penetration distance versus viscosity observed in the simulations and leads to a phenomenological model for the spreading of injury with multiple shocks. For a reasonable selection of a single efficiency parameter, this model is able to reproduce in vivo observations of an apparent 1000-shock threshold before wide-spread tissue injury occurs in targeted kidneys and the approximate extent of this injury after a typical clinical dose of 2000 shock waves.

Shock-induced bubble jetting into a viscous fluid with application to tissue injury in shock-wave lithotripsy

PubMed Central

Freund, J. B.; Shukla, R. K.; Evan, A. P.

2009-01-01

Shock waves in liquids are known to cause spherical gas bubbles to rapidly collapse and form strong re-entrant jets in the direction of the propagating shock. The interaction of these jets with an adjacent viscous liquid is investigated using finite-volume simulation methods. This configuration serves as a model for tissue injury during shock-wave lithotripsy, a medical procedure to remove kidney stones. In this case, the viscous fluid provides a crude model for the tissue. It is found that for viscosities comparable to what might be expected in tissue, the jet that forms upon collapse of a small bubble fails to penetrate deeply into the viscous fluid “tissue.” A simple model reproduces the penetration distance versus viscosity observed in the simulations and leads to a phenomenological model for the spreading of injury with multiple shocks. For a reasonable selection of a single efficiency parameter, this model is able to reproduce in vivo observations of an apparent 1000-shock threshold before wide-spread tissue injury occurs in targeted kidneys and the approximate extent of this injury after a typical clinical dose of 2000 shock waves. PMID:19894850

Investigating the viscous interaction and its role in generating the ionospheric potential during the Whole Heliosphere Interval

NASA Astrophysics Data System (ADS)

Bruntz, R.; Lopez, R. E.; Bhattarai, S. K.; Pham, K. H.; Deng, Y.; Huang, Y.; Wiltberger, M.; Lyon, J. G.

2012-07-01

The Whole Heliosphere Interval (WHI), comprising March 20-April 16, 2008 (DOY 80-107), is a single Carrington Rotation (2068) designated for intense study through observations and simulations. We used solar wind data from the WHI to run the Coupled Magnetosphere-Ionosphere-Thermosphere (CMIT) and stand-alone Lyon-Fedder-Mobarry (LFM) models. The LFM model was also run with the WHI solar wind plasma parameters but with zero interplanetary magnetic field (IMF). With no IMF, we expect that the cross-polar cap potential (CPCP) is due entirely to the viscous interaction. Comparing the LFM runs with and without the IMF, we found that during strong driving with southward IMF Bz, the viscous potential could be a significant fraction of the total CPCP. During times of northward IMF Bz, the CPCP was generally lower than the CPCP value from the IMF=0 run. LFM tends to produce high polar cap potentials, but by using the Bruntz et al. (2012) viscous potential formula (ΦV=μn0.439V1.33, where μ=0.00431) and the IMF=0 LFM run, we calculated a scaling factor γ=1.54, which can be used to scale the LFM CPCP during the WHI down to realistic values. The Newell et al. (2008) viscous merging term can similarly be used to predict the viscous potential using the formula: ΦV=νn1/2V2, where the value ν=6.39×10-5 was also found using the zero IMF run. Both formulas were found to perform better when V (solar wind)=Vx, rather than Vtotal, yielding similar, accurate predictions of the LFM viscous potential, with R2>0.91 for both formulas. The γ factor was also used to scale down the LFM CPCP from the full solar wind run, with most of the resultant values matching the CPCP from the Weimer05 model well, even though γ was derived independent of the Weimer05 model or the full LFM data. We interpret this to be an indication that the conductivity model in LFM is producing values that are too low, thus elevating the CPCP values.

Lubrication pressure and fractional viscous damping effects on the spring-block model of earthquakes

NASA Astrophysics Data System (ADS)

Tanekou, G. B.; Fogang, C. F.; Kengne, R.; Pelap, F. B.

2018-04-01

We examine the dynamical behaviours of the "single mass-spring" model for earthquakes considering lubrication pressure effects on pre-existing faults and viscous fractional damping. The lubrication pressure supports a part of the load, thereby reducing the normal stress and the associated friction across the gap. During the co-seismic phase, all of the strain accumulated during the inter-seismic duration does not recover; a fraction of this strain remains as a result of viscous relaxation. Viscous damping friction makes it possible to study rocks at depth possessing visco-elastic behaviours. At increasing depths, rock deformation gradually transitions from brittle to ductile. The fractional derivative is based on the properties of rocks, including information about previous deformation events ( i.e., the so-called memory effect). Increasing the fractional derivative can extend or delay the transition from stick-slip oscillation to a stable equilibrium state and even suppress it. For the single block model, the interactions of the introduced lubrication pressure and viscous damping are found to give rise to oscillation death, which corresponds to aseismic fault behaviour. Our result shows that the earthquake occurrence increases with increases in both the damping coefficient and the lubrication pressure. We have also revealed that the accumulation of large stresses can be controlled via artificial lubrication.

Drinking with a hairy tongue: viscous entrainment by dipping hairy surfaces

NASA Astrophysics Data System (ADS)

Nasto, Alice; Brun, Pierre-Thomas; Alvarado, José; Bush, John; Hosoi, Anette

2016-11-01

Nectar-drinking bats have tongues covered with hair-like papillae, enhancing their ability to take up viscous nectar by dipping. Using a combination of model experiments and theory reminiscent of Landau-Levich-Derjaguin dip coating, we rationalize this mechanism of viscous entrainment in a hairy texture. For the model experiments, hairy surfaces are fabricated using laser cut molds and casting samples with PDMS elastomer. Modeling the liquid trapped within the texture using a Darcy-Brinkman like approach, we derive the drainage flow solution. The amount of fluid that is entrained is dependent on the viscosity of the fluid, the density of the hairs, and the dipping speed. We find that there is an optimal hair density to maximize fluid uptake.

Hydrodynamic characteristics of viscous fluid flow in screw channels formed by two ribs

NASA Astrophysics Data System (ADS)

Kadyirov, A. I.; Abaydullin, B. R.; Vachagina, E. K.

2018-03-01

The mathematical model of laminar viscous flows in screw channels, formed by two ribs, is developed using the helical coordinates. The numerical results of the flow with helical symmetry are presented.

Modeling and Control of a Tethered Rotorcraft

DTIC Science & Technology

2010-07-30

viscous damper with damping coefficient Cv. Visco-elastic line force is written in terms of components Δx, Δy, and Δz, of the difference vector formed...tether drag coefficient CS = tether damping coefficient Cv = viscous damping coefficient d = diameter of the tether En = n x n identity matrix FA...matrix consisting of Iyy and Izz k = rotor head stiffness KLAT, KLON = steady state flapping gains Ks, Kv = static and viscous stiffness Lj

Kevlar reinforced neoprene composites

NASA Technical Reports Server (NTRS)

Penn, B. G.; Daniels, J. G.; White, W. T.; Thompson, L. M.; Clemons, L. M.

1985-01-01

Kevlar/neoprene composites were prepared by two techniques. One method involved the fabrication of a composite from a rubber prepreg prepared by coating Kevlar with viscous neoprene solution and then allowing the solvent to evaporate (solution impregnation technique). The second method involved heating a stack of Kevlar/neoprene sheets at a temperature sufficient to cause polymer flow (melt flow technique). There was no significant difference in the breaking strength and percent elongation for samples obtained by the two methods; however the shear strength obtained for samples fabricated by the solution impregnation technique (275 psi) was significantly higher than that found for the melt flow fabricated samples (110 psi).

Viscous entrainment on hairy surfaces

NASA Astrophysics Data System (ADS)

Nasto, Alice; Brun, P.-T.; Hosoi, A. E.

2018-02-01

Nectar-drinking bats and honeybees have tongues covered with hairlike structures, enhancing their ability to take up viscous nectar by dipping. Using a combination of model experiments and theory, we explore the physical mechanisms that govern viscous entrainment in a hairy texture. Hairy surfaces are fabricated using laser cut molds and casting samples with polydimethylsiloxane (PDMS) elastomer. We model the liquid trapped within the texture using a Darcy-Brinkmann-like approach and derive the drainage flow solution. The amount of fluid that is entrained is dependent on the viscosity of the fluid, the density of the hairs, and the withdrawal speed. Both experiments and theory reveal an optimal hair density to maximize fluid uptake.

A soft porous drop in linear flows

NASA Astrophysics Data System (ADS)

Young, Yuan-Nan; Miksis, Michael; Mori, Yoichiro; Shelley, Michael

2017-11-01

The cellular cytoplasm consists a viscous fluid filled with fibrous networks that also have their own dynamics. Such fluid-structure interactions have been modeled as a soft porous material immersed in a viscous fluid. In this talk we focus on the hydrodynamics of a viscous drop filled with soft porous material inside. Suspended in a Stokes flow, such a porous viscous drop is allowed to deform, both the drop interface and the porous structures inside. Special focus is on the deformation dynamics of both the porosity and the shape of the drop under simple flows such as a uniform streaming flow and linear flows. We examine the effects of flow boundary conditions at interface between the porous drop and the surrounding viscous fluid. We also examine the dynamics of a porous drop with active stress from the porous network.

Alveolar Thin Layer Flows and Surfactant Dynamics

NASA Astrophysics Data System (ADS)

Roumie, Ahmad; Jbaily, Abdulrahman; Szeri, Andrew J.

2017-11-01

Pulmonary surfactants play a vital role in everyday respiration. They regulate surface tension in the lungs by diffusing through the hypophase, a liquid layer that lines the interior surface of the alveoli, and adsorbing to the existing air-fluid interface. This decreases the equilibrium surface tension value by as much as a factor of 3, minimizing breathing effort and preventing lung collapse at the end of exhalation. Given that the hypophase thickness h lies within the range 0.1 μm < h <0.5 μm , and that the average alveolar radius R is 100 μm , for some purposes the hypophase may usefully be modeled as a fluid layer on a flat sheet representing the alveolar wall. Moreover, because of the large aspect ratio, the lubrication approximation can be applied. The aim of the present work is to study the interaction between the straining of the alveolar wall and the fluid flow in the hypophase. The analysis is governed by the relative magnitudes of the time scales of surfactant diffusion, adsorption, desorption, viscous dissipation and sheet straining. Cases of particular interest include non-uniform surfactant concentration at the interface, leading to Marangoni flows and a non-uniform hypophase thickness profile. The analytical formulation and numerical simulations are presented. This work is motivated by a need to understand alveolar deformation during breathing, and to do so in a way that derives from improved understanding of the fluid mechanics of the problem.

Generation of region 1 current by magnetospheric pressure gradients

NASA Technical Reports Server (NTRS)

Yang, Y. S.; Spiro, R. W.; Wolf, R. A.

1994-01-01

The Rice Convection Model (RCM) is used to illustrate theoretical possibilities for generating region 1 Birkeland currents by pressure gradients on closed field lines in the Earth's magnetosphere. Inertial effects and viscous forces are neglected. The RCM is applied to idealized cases, to emphasize the basic physical ideas rather than realistic representation of the actual magnetosphere. Ionospheric conductance is taken to be uniform, and the simplest possible representations of the magnetospheric plasma are used. Three basic cases are considered: (1) the case of pure northward Interplanetary Magnetic Field (IMF), with cusp merging assumed to create new closed field lines near the nose of the magnetosphere, following the suggestion by Song and Russell (1992); (2) the case where Dungey-type reconnection occurs at the nose, but magnetosheath plasma somehow enters closed field lines on the dawnside and duskside of the merging region, causing a pressure-driven low-latitude boundary layer; and (3) the case where Dungey-type reconnection occurs at the nose, but region 1 currents flow on sunward drifting plasma sheet field lines. In case 1, currents of region 1 sense are generated by pressure gradients, but those currents do not supply the power for ionospheric convection. Results for case 2 suggest that pressure gradients at the inner edge of the low-latitude boundary layer might generate a large fraction of the region 1 Birkeland currents that drive magnetospheric convection. Results for case 3 indicate that pressure gradients in the plasma sheet could provide part of the region 1 current.

Procedures for the computation of unsteady transonic flows including viscous effects

NASA Technical Reports Server (NTRS)

Rizzetta, D. P.

1982-01-01

Modifications of the code LTRAN2, developed by Ballhaus and Goorjian, which account for viscous effects in the computation of planar unsteady transonic flows are presented. Two models are considered and their theoretical development and numerical implementation is discussed. Computational examples employing both models are compared with inviscid solutions and with experimental data. Use of the modified code is described.

Protoplanetary disc `isochrones' and the evolution of discs in the M˙-Md plane

NASA Astrophysics Data System (ADS)

Lodato, Giuseppe; Scardoni, Chiara E.; Manara, Carlo F.; Testi, Leonardo

2017-12-01

In this paper, we compare simple viscous diffusion models for the disc evolution with the results of recent surveys of the properties of young protoplanetary discs. We introduce the useful concept of 'disc isochrones' in the accretion rate-disc mass plane and explore a set of Monte Carlo realization of disc initial conditions. We find that such simple viscous models can provide a remarkable agreement with the available data in the Lupus star forming region, with the key requirement that the average viscous evolutionary time-scale of the discs is comparable to the cluster age. Our models produce naturally a correlation between mass accretion rate and disc mass that is shallower than linear, contrary to previous results and in agreement with observations. We also predict that a linear correlation, with a tighter scatter, should be found for more evolved disc populations. Finally, we find that such viscous models can reproduce the observations in the Lupus region only in the assumption that the efficiency of angular momentum transport is a growing function of radius, thus putting interesting constraints on the nature of the microscopic processes that lead to disc accretion.

A model for cytoplasmic rheology consistent with magnetic twisting cytometry.

PubMed

Butler, J P; Kelly, S M

1998-01-01

Magnetic twisting cytometry is gaining wide applicability as a tool for the investigation of the rheological properties of cells and the mechanical properties of receptor-cytoskeletal interactions. Current technology involves the application and release of magnetically induced torques on small magnetic particles bound to or inside cells, with measurements of the resulting angular rotation of the particles. The properties of purely elastic or purely viscous materials can be determined by the angular strain and strain rate, respectively. However, the cytoskeleton and its linkage to cell surface receptors display elastic, viscous, and even plastic deformation, and the simultaneous characterization of these properties using only elastic or viscous models is internally inconsistent. Data interpretation is complicated by the fact that in current technology, the applied torques are not constant in time, but decrease as the particles rotate. This paper describes an internally consistent model consisting of a parallel viscoelastic element in series with a parallel viscoelastic element, and one approach to quantitative parameter evaluation. The unified model reproduces all essential features seen in data obtained from a wide variety of cell populations, and contains the pure elastic, viscoelastic, and viscous cases as subsets.

46 CFR 160.060-1 - Incorporation by reference.

Code of Federal Regulations, 2010 CFR

2010-10-01

...-1: Sheet 1—Cutting Pattern and General Arrangement, Model AY. Sheet 2—Cutting Pattern and General Arrangement, Model CYM. Sheet 3—Cutting Pattern and General Arrangement, Model CYS. Sheet 4—Insert Pattern, Model AY. Sheet 5—Insert Pattern, Model CYM. Sheet 6—Insert Pattern, Model CYS. (c) Copies on file...

Viscothermal Coupling Effects on Sound Attenuation in Concentrated Colloidal Dispersions.

NASA Astrophysics Data System (ADS)

Han, Wei

1995-11-01

This thesis describes a Unified Coupled Phase Continuum (UCPC) model to analyze sound propagation through aerosols, emulsions and suspensions in terms of frequency dependent attenuation coefficient and sound speed. Expressions for the viscous and thermal coupling coefficients explicitly account for the effects of particle size, shape factor, orientation as well as concentration and the sound frequency. The UCPC model also takes into account the intrinsic acoustic absorption within the fluid medium due to its viscosity and heat conductivity. The effective complex wave number as a function of frequency is derived. A frequency- and concentration-dependent complex Nusselt number for the interfacial thermal coupling coefficient is derived using an approximate similarity between the 'viscous skin drag' and 'heat conduction flux' associated with the discontinuous suspended phase, on the basis of a cell model. The theoretical predictions of attenuation spectra provide satisfactory agreement with reported experimental data on two concentrated suspensions (polystyrene latex and kaolin pigment), two concentrated emulsions (toluene -in-water, n-hexadecane-in-water), and two aerosols (oleic acid droplets-in-nitrogen, alumina-in-air), covering a wide range of relative magnitudes (from 10^ {-3} to 10^{3}) of thermal versus viscous contributions, for dispersed phase volume fractions as high as 50%. The relative differences between the additive result of separate viscous and thermal loss estimates and combined viscothermal absorption results are also presented. Effects of particle shape on viscous attenuation of sound in concentrated suspensions of non-spherical clay particles are studied. Attenuation spectra for 18 frequencies from 3 to 100 MHz are measured and analyzed for eleven kaolin clay slurries with solid concentrations ranging from 0.6% to 35% (w/w). A modified viscous drag coefficient that considers frequency, concentration, particle size, shape and orientation of spheroids, is developed and applied to estimate the viscous attenuation coefficients. With incorporation of particle size and shape distributions (PSSD), predictions agree quantitatively with observed attenuation coefficients. The effects of particle aspect ratio and orientation become more evident as particle concentrations and frequencies are increased. The UCPC model combined with the ultrasonic spectroscopy techniques can provide for theoretical and experimental frameworks in characterization of concentrated colloidal dispersions.

A critical assessment of viscous models of trench topography and corner flow

NASA Technical Reports Server (NTRS)

Zhang, J.; Hager, B. H.; Raefsky, A.

1984-01-01

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

An assessment of viscous effects in computational simulation of benign and burst vortex flows on generic fighter wind-tunnel models using TEAM code

NASA Technical Reports Server (NTRS)

Kinard, Tim A.; Harris, Brenda W.; Raj, Pradeep

1995-01-01

Vortex flows on a twin-tail and a single-tail modular transonic vortex interaction (MTVI) model, representative of a generic fighter configuration, are computationally simulated in this study using the Three-dimensional Euler/Navier-Stokes Aerodynamic Method (TEAM). The primary objective is to provide an assessment of viscous effects on benign (10 deg angle of attack) and burst (35 deg angle of attack) vortex flow solutions. This study was conducted in support of a NASA project aimed at assessing the viability of using Euler technology to predict aerodynamic characteristics of aircraft configurations at moderate-to-high angles of attack in a preliminary design environment. The TEAM code solves the Euler and Reynolds-average Navier-Stokes equations on patched multiblock structured grids. Its algorithm is based on a cell-centered finite-volume formulation with multistage time-stepping scheme. Viscous effects are assessed by comparing the computed inviscid and viscous solutions with each other and experimental data. Also, results of Euler solution sensitivity to grid density and numerical dissipation are presented for the twin-tail model. The results show that proper accounting of viscous effects is necessary for detailed design and optimization but Euler solutions can provide meaningful guidelines for preliminary design of flight vehicles which exhibit vortex flows in parts of their flight envelope.

Viscous anisotropy of textured olivine aggregates: 2. Micromechanical model

NASA Astrophysics Data System (ADS)

Hansen, Lars N.; Conrad, Clinton P.; Boneh, Yuval; Skemer, Philip; Warren, Jessica M.; Kohlstedt, David L.

2016-10-01

The significant viscous anisotropy that results from crystallographic alignment (texture) of olivine grains in deformed upper mantle rocks strongly influences a large variety of geodynamic processes. Our ability to explore the effects of anisotropic viscosity in simulations of these processes requires a mechanical model that can predict the magnitude of anisotropy and its evolution. Unfortunately, existing models of olivine textural evolution and viscous anisotropy are calibrated for relatively small deformations and simple strain paths, making them less general than desired for many large-scale geodynamic scenarios. Here we develop a new set of micromechanical models to describe the mechanical behavior and textural evolution of olivine through a large range of strains and complex strain histories. For the mechanical behavior, we explore two extreme scenarios, one in which each grain experiences the same stress tensor (Sachs model) and one in which each grain undergoes a strain rate as close as possible to the macroscopic strain rate (pseudo-Taylor model). For the textural evolution, we develop a new model in which the director method is used to control the rate of grain rotation and the available slip systems in olivine are used to control the axis of rotation. Only recently has enough laboratory data on the deformation of olivine become available to calibrate these models. We use these new data to conduct inversions for the best parameters to characterize both the mechanical and textural evolution models. These inversions demonstrate that the calibrated pseudo-Taylor model best reproduces the mechanical observations. Additionally, the pseudo-Taylor textural evolution model can reasonably reproduce the observed texture strength, shape, and orientation after large and complex deformations. A quantitative comparison between our calibrated models and previously published models reveals that our new models excel in predicting the magnitude of viscous anisotropy and the details of the textural evolution. In addition, we demonstrate that the mechanical and textural evolution models can be coupled and used to reproduce mechanical evolution during large-strain torsion tests. This set of models therefore provides a new geodynamic tool for incorporating viscous anisotropy into large-scale numerical simulations.

Numerical study of Free Convective Viscous Dissipative flow along Vertical Cone with Influence of Radiation using Network Simulation method

NASA Astrophysics Data System (ADS)

Kannan, R. M.; Pullepu, Bapuji; Immanuel, Y.

2018-04-01

A two dimensional mathematical model is formulated for the transient laminar free convective flow with heat transfer over an incompressible viscous fluid past a vertical cone with uniform surface heat flux with combined effects of viscous dissipation and radiation. The dimensionless boundary layer equations of the flow which are transient, coupled and nonlinear Partial differential equations are solved using the Network Simulation Method (NSM), a powerful numerical technique which demonstrates high efficiency and accuracy by employing the network simulator computer code Pspice. The velocity and temperature profiles have been investigated for various factors, namely viscous dissipation parameter ε, Prandtl number Pr and radiation Rd are analyzed graphically.

The thinning of viscous liquid threads.

NASA Astrophysics Data System (ADS)

Castrejon-Pita, J. Rafael; Castrejon-Pita, Alfonso A.; Hutchings, Ian M.

2012-11-01

The thinning neck of dripping droplets is studied experimentally for viscous Newtonian fluids. High speed imaging is used to measure the minimum neck diameter in terms of the time τ to breakup. Mixtures of water and glycerol with viscosities ranging from 20 to 363 mPa s are used to model the Newtonian behavior. The results show the transition from potential to inertial-viscous regimes occurs at the predicted values of ~Oh2. Before this transition the neck contraction rate follows the inviscid scaling law ~τ 2 / 3 . After the transition, the neck thinning tends towards the linear viscous scaling law ~ τ . Project supported by the EPSRC-UK (EP/G029458/1) and Cambridge-KACST.

Formation of magnetic discontinuities through viscous relaxation

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

Kumar, Sanjay; Bhattacharyya, R.; Smolarkiewicz, P. K.

2014-05-15

According to Parker's magnetostatic theorem, tangential discontinuities in magnetic field, or current sheets (CSs), are generally unavoidable in an equilibrium magnetofluid with infinite electrical conductivity and complex magnetic topology. These CSs are due to a failure of a magnetic field in achieving force-balance everywhere and preserving its topology while remaining in a spatially continuous state. A recent work [Kumar, Bhattacharyya, and Smolarkiewicz, Phys. Plasmas 20, 112903 (2013)] demonstrated this CS formation utilizing numerical simulations in terms of the vector magnetic field. The magnetohydrodynamic simulations presented here complement the above work by demonstrating CS formation by employing a novel approach ofmore » describing the magnetofluid evolution in terms of magnetic flux surfaces instead of the vector magnetic field. The magnetic flux surfaces being the possible sites on which CSs develop, this approach provides a direct visualization of the CS formation, helpful in understanding the governing dynamics. The simulations confirm development of tangential discontinuities through a favorable contortion of magnetic flux surfaces, as the magnetofluid undergoes a topology-preserving viscous relaxation from an initial non-equilibrium state with twisted magnetic field. A crucial finding of this work is in its demonstration of CS formation at spatial locations away from the magnetic nulls.« less

From viscous to elastic sheets: Dynamics of smectic bubbles

NASA Astrophysics Data System (ADS)

Harth, Kirsten; Trittel, Torsten; van der Meer, Devaraj; Stannarius, Ralf

2015-11-01

Oscillations and rupture of bubbles composed of an inner fluid separated from an outer fluid by a membrane, represent an old but still immensely active field of research. Membrane properties apart from surface tension are often neglected for fluids (e.g. soap bubbles), whereas they govern the dynamics in systems with a rigid membrane (e.g. vesicles). Due to their layered phase structure, smectic liquid crystals can form stable, uniform and easy-to-handle fluid films of immense aspect ratios. Only recently, freely floating bubbles detached from a support could be prepared. We analyze their relaxation from strongly non-spherical shapes and the rupture using high-speed video recordings. Peculiar dynamics intermediate between simple viscous fluid films and an elastic response are observed: Fast oscillations, slowed relaxation and even the reversible formation of wrinkles and extrusions. Bubble rupture deviates qualitatively from previously observed behavior of simple Newtonian and other complex fluids. It becomes retarded by at least two orders of magnitude compared to the predictions of Taylor and Culick. A transition between fluid-like and elastic behavior is seen with increasing thickness. We give experimental results, an intuitive explanation and a novel hydrodynamic description.

On a viscous critical-stress model of martensitic phase transitions

NASA Astrophysics Data System (ADS)

Weatherwax, John; Vaynblat, Dimitri; Bruno, Oscar; Rosales, Ruben

2007-09-01

The solid-to-solid phase transitions that result from shock loading of certain materials, such as the graphite-to-diamond transition and the α-ɛ transition in iron, have long been subjects of a substantial theoretical and experimental literature. Recently a model for such transitions was introduced which, based on a CS condition (CS) and without use of fitting parameters, accounts quantitatively for existing observations in a number of systems [Bruno and Vaynblat, Proc. R. Soc. London, Ser. A 457, 2871 (2001)]. While the results of the CS model match the main features of the available experimental data, disagreements in some details between the predictions of this model and experiment, attributable to an ideal character of the CS model, do exist. In this article we present a version of the CS model, the viscous CS model (vCS), as well as a numerical method for its solution. This model and the corresponding solver results in a much improved overall CS modeling capability. The innovations we introduce include: (1) Enhancement of the model by inclusion of viscous phase-transition effects; as well as a numerical solver that allows for a fully rigorous treatment of both, the (2) Rarefaction fans (which had previously been approximated by "rarefaction discontinuities"), and (3) viscous phase-transition effects, that are part of the vCS model. In particular we show that the vCS model accounts accurately for well known "gradual" rises in the α-ɛ transition which, in the original CS model, were somewhat crudely approximated as jump discontinuities.

Effects of viscous pressure on warm inflationary generalized cosmic Chaplygin gas model

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

Sharif, M.; Saleem, Rabia, E-mail: msharif.math@pu.edu.pk, E-mail: rabiasaleem1988@yahoo.com

This paper is devoted to study the effects of bulk viscous pressure on an inflationary generalized cosmic Chaplygin gas model using FRW background. The matter contents of the universe are assumed to be inflaton and imperfect fluid. We evaluate inflaton fields, potentials and entropy density for variable as well as constant dissipation and bulk viscous coefficients in weak as well as high dissipative regimes during intermediate era. In order to discuss inflationary perturbations, we evaluate entropy density, scalar (tensor) power spectra, their corresponding spectral indices, tensor-scalar ratio and running of spectral index in terms of inflaton which are constrained usingmore » recent Planck, WMAP7 and Bicep2 probes.« less

Three-dimensional viscous rotor flow calculations using a viscous-inviscid interaction approach

NASA Technical Reports Server (NTRS)

Chen, Ching S.; Bridgeman, John O.

1990-01-01

A three-dimensional viscous-inviscid interaction analysis was developed to predict the performance of rotors in hover and in forward flight at subsonic and transonic tip speeds. The analysis solves the full-potential and boundary-layer equations by finite-difference numerical procedures. Calculations were made for several different model rotor configurations. The results were compared with predictions from a two-dimensional integral method and with experimental data. The comparisons show good agreement between predictions and test data.

Extension of the dielectric breakdown model for simulation of viscous fingering at finite viscosity ratios.

PubMed

Doorwar, Shashvat; Mohanty, Kishore K

2014-07-01

Immiscible displacement of viscous oil by water in a petroleum reservoir is often hydrodynamically unstable. Due to similarities between the physics of dielectric breakdown and immiscible flow in porous media, we extend the existing dielectric breakdown model to simulate viscous fingering patterns for a wide range of viscosity ratios (μ(r)). At low values of power-law index η, the system behaves like a stable Eden growth model and as the value of η is increased to unity, diffusion limited aggregation-like fractals appear. This model is compared with our two-dimensional (2D) experiments to develop a correlation between the viscosity ratio and the power index, i.e., η = 10(-5)μ(r)(0.8775). The 2D and three-dimensional (3D) simulation data appear scalable. The fingering pattern in 3D simulations at finite viscosity ratios appear qualitatively similar to the few experimental results published in the literature.

Sensor Phenomenology and Feature Development for Improved Sonar-based Detection & Classification of Underwater UXO

DTIC Science & Technology

2009-12-01

the validity of approximating poroelastic media with acoustic or acoustic /elastic models , and to characterize how scattering physics will differ for...elastic buried object (yellow rectangle in the figure) in three types of environments: • (1) Model 1: acoustic layer on top of a poroelastic medium with a...porosity gradient and no viscous damping. • (2) Model 2: acoustic layer on top of a poroelastic medium with a porosity gradient and viscous damping

Microfluidic System Simulation Including the Electro-Viscous Effect

NASA Technical Reports Server (NTRS)

Rojas, Eileen; Chen, C. P.; Majumdar, Alok

2007-01-01

This paper describes a practical approach using a general purpose lumped-parameter computer program, GFSSP (Generalized Fluid System Simulation Program) for calculating flow distribution in a network of micro-channels including electro-viscous effects due to the existence of electrical double layer (EDL). In this study, an empirical formulation for calculating an effective viscosity of ionic solutions based on dimensional analysis is described to account for surface charge and bulk fluid conductivity, which give rise to electro-viscous effect in microfluidics network. Two dimensional slit micro flow data was used to determine the model coefficients. Geometry effect is then included through a Poiseuille number correlation in GFSSP. The bi-power model was used to calculate flow distribution of isotropically etched straight channel and T-junction microflows involving ionic solutions. Performance of the proposed model is assessed against experimental test data.

Constraining Lithospheric Stress on Mars From Mars Global Surveyor (MGS) Topography, Gravity, and Crustal Thickness

NASA Astrophysics Data System (ADS)

Dimitrova, L. L.; Holt, W. E.; Haines, A. J.

2004-12-01

The quantification of lithospheric dynamics on Mars is of fundamental importance to the understanding of Martian geologic history and surface morphology. The global stress field associated with gravitational potential energy differences (GPE) constitutes a significant fraction of the total stress field. We have obtained 0.25 × 0.25 degree data sets of MOLA topography and crustal thickness from Zuber et al 2000. We calculate the GPE associated with the topography by vertically integrating density to a given lithospheric depth using either Airy isostasy assumptions or the crustal thickness model of Zuber et al 2000. Using a finite-element thin sheet method we solve the full 3-D force-balance equations for stress magnitudes and orientations within the lithosphere associated with the horizontal gradients in GPE. We assume ρ crust = 2900~ kg m-3, ~ρ mantle = 3500~ kg m-3, ~g = 3.7~ ms-2 and various lithospheric thicknesses. We explore both viscous and elastic rheologies. It is interesting to note that our solutions depend on the rheology only through the ratio of shear and bulk moduli for the elastic case or the shear and bulk viscosities for the viscous case. Finally, we also calculate the expected style and orientation of the associated elastic strain. To first order, all stress field solutions are consistent with a tectonically inactive region relaxing due to excess or deficit of mass. Thus topographic highs e.g., Tharsis Mons, Olympus Mons, and Alba Patera, are in deviatoric extension, while topographic lows e.g., Valles Marineras and impact basins, are in deviatoric compression. At short wavelengths, however, features are regionally supported and differences in the stress field solutions occur between the Airy isostasy and the crustal thickness models. For example, several lowlands e.g., Isidis and Argyre Planitae, are in deviatoric extension in the crustal thickness model. Increasing the lithospheric depth increases the magnitudes of stresses but does not change the deformation style. Comparisons of viscous and elastic solutions with the same ratio of the shear and bulk viscosities for the viscous case and shear and bulk moduli for the elastic case show small variations in stress magnitudes and orientations among the solutions. We compare our results with inferred fault styles and orientations e.g., Mege 2001, Montessi and Zuber 2003. The GPE variations appear to match the radial graben features on Tharsis, as well as the wrinkle ridges at lower elevations. However, they do not predict the wrinkle ridges on the flanks of Tharsis and some grabens around Alba Patera extend into the lower topography regions beyond the area where the GPE associated stresses predict extension. Calculations of Banerdt and Glombek 2000 show the importance of flexural loading on the stress field. Their flexural loading solution provides compression at higher elevations on Tharsis that may explain some of the wrinkle ridges not predicted by GPE variations. Therefore, it appears that a combined solution of GPE inferred vertically integrated deviatoric stresses and the flexurally produced stresses will provide a better fit to the fault styles and orientations observed on Mars.

Continuum calculations of continental deformation in transcurrent environments

NASA Technical Reports Server (NTRS)

Sonder, L. J.; England, P. C.; Houseman, G. A.

1986-01-01

A thin viscous sheet approximation is used to investigate continental deformation near a strike-slip boundary. The vertically averaged velocity field is calculated for a medium characterized by a power law rheology with stress exponent n. Driving stresses include those applied along boundaries of the sheet and those arising from buoyancy forces related to lateral differences in crustal thickness. Exact and approximate analytic solutions for a region with a sinusoidal strike-slip boundary condition are compared with solutions for more geologically relevant boundary conditions obtained using a finite element technique. The across-strike length scale of the deformation is approximately 1/4pi x sq rt n times the dominant wavelength of the imposed strike-slip boundary condition for both the analytic and the numerical solutions; this result is consistent with length scales observed in continental regions of large-scale transcurrent faulting. An approximate, linear relationship between displacement and rotation is found that depends only on the deformation length scale and the rheology. Calculated displacements, finite rotations, and distribution of crustal thicknesses are consistent with those observed in the region of the Pacific-North America plate boundary in California.

Modeling of viscoelastic properties of nonpermeable porous rocks saturated with highly viscous fluid at seismic frequencies at the core scale

NASA Astrophysics Data System (ADS)

Wang, Zizhen; Schmitt, Douglas R.; Wang, Ruihe

2017-08-01

A core scale modeling method for viscoelastic properties of rocks saturated with viscous fluid at low frequencies is developed based on the stress-strain method. The elastic moduli dispersion of viscous fluid is described by the Maxwell's spring-dash pot model. Based on this modeling method, we numerically test the effects of frequency, fluid viscosity, porosity, pore size, and pore aspect ratio on the storage moduli and the stress-strain phase lag of saturated rocks. And we also compared the modeling results to the Hashin-Shtrikman bounds and the coherent potential approximation (CPA). The dynamic moduli calculated from the modeling are lower than the predictions of CPA, and both of these fall between the Hashin-Shtrikman bounds. The modeling results indicate that the frequency and the fluid viscosity have similar effects on the dynamic moduli dispersion of fully saturated rocks. We observed the Debye peak in the phase lag variation with the change of frequency and viscosity. The pore structure parameters, such as porosity, pore size, and aspect ratio affect the rock frame stiffness and result in different viscoelastic behaviors of the saturated rocks. The stress-strain phase lags are larger with smaller stiffness contrasts between the rock frame and the pore fluid. The viscoelastic properties of saturated rocks are more sensitive to aspect ratio compared to other pore structure parameters. The results suggest that significant seismic dispersion (at about 50-200 Hz) might be expected for both compressional and shear waves passing through rocks saturated with highly viscous fluids.