Viscosity-adjusted estimation of pressure head and pump flow with quasi-pulsatile modulation of rotary blood pump for a total artificial heart.

PubMed

Yurimoto, Terumi; Hara, Shintaro; Isoyama, Takashi; Saito, Itsuro; Ono, Toshiya; Abe, Yusuke

2016-09-01

Estimation of pressure and flow has been an important subject for developing implantable artificial hearts. To realize real-time viscosity-adjusted estimation of pressure head and pump flow for a total artificial heart, we propose the table estimation method with quasi-pulsatile modulation of rotary blood pump in which systolic high flow and diastolic low flow phased are generated. The table estimation method utilizes three kinds of tables: viscosity, pressure and flow tables. Viscosity is estimated from the characteristic that differential value in motor speed between systolic and diastolic phases varies depending on viscosity. Potential of this estimation method was investigated using mock circulation system. Glycerin solution diluted with salty water was used to adjust viscosity of fluid. In verification of this method using continuous flow data, fairly good estimation could be possible when differential pulse width modulation (PWM) value of the motor between systolic and diastolic phases was high. In estimation under quasi-pulsatile condition, inertia correction was provided and fairly good estimation was possible when the differential PWM value was high, which was not different from the verification results using continuous flow data. In the experiment of real-time estimation applying moving average method to the estimated viscosity, fair estimation could be possible when the differential PWM value was high, showing that real-time viscosity-adjusted estimation of pressure head and pump flow would be possible with this novel estimation method when the differential PWM value would be set high.

Discontinuous Galerkin finite element method for the nonlinear hyperbolic problems with entropy-based artificial viscosity stabilization

NASA Astrophysics Data System (ADS)

Zingan, Valentin Nikolaevich

This work develops a discontinuous Galerkin finite element discretization of non- linear hyperbolic conservation equations with efficient and robust high order stabilization built on an entropy-based artificial viscosity approximation. The solutions of equations are represented by elementwise polynomials of an arbitrary degree p > 0 which are continuous within each element but discontinuous on the boundaries. The discretization of equations in time is done by means of high order explicit Runge-Kutta methods identified with respective Butcher tableaux. To stabilize a numerical solution in the vicinity of shock waves and simultaneously preserve the smooth parts from smearing, we add some reasonable amount of artificial viscosity in accordance with the physical principle of entropy production in the interior of shock waves. The viscosity coefficient is proportional to the local size of the residual of an entropy equation and is bounded from above by the first-order artificial viscosity defined by a local wave speed. Since the residual of an entropy equation is supposed to be vanishingly small in smooth regions (of the order of the Local Truncation Error) and arbitrarily large in shocks, the entropy viscosity is almost zero everywhere except the shocks, where it reaches the first-order upper bound. One- and two-dimensional benchmark test cases are presented for nonlinear hyperbolic scalar conservation laws and the system of compressible Euler equations. These tests demonstrate the satisfactory stability properties of the method and optimal convergence rates as well. All numerical solutions to the test problems agree well with the reference solutions found in the literature. We conclude that the new method developed in the present work is a valuable alternative to currently existing techniques of viscous stabilization.

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.

An analysis of artificial viscosity effects on reacting flows using a spectral multi-domain technique

NASA Technical Reports Server (NTRS)

Macaraeg, M. G.; Streett, C. L.; Hussaini, M. Y.

1987-01-01

Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.

Measurement of the dynamic viscosity of hybrid engine oil -Cuo-MWCNT nanofluid, development of a practical viscosity correlation and utilizing the artificial neural network

NASA Astrophysics Data System (ADS)

Aghaei, Alireza; Khorasanizadeh, Hossein; Sheikhzadeh, Ghanbar Ali

2018-01-01

The main objectives of this study have been measurement of the dynamic viscosity of CuO-MWCNTs/SAE 5w-50 hybrid nanofluid, utilization of artificial neural networks (ANN) and development of a new viscosity model. The new nanofluid has been prepared by a two-stage procedure with volume fractions of 0.05, 0.1, 0.25, 0.5, 0.75 and 1%. Then, utilizing a Brookfield viscometer, its dynamic viscosity has been measured for temperatures of 5, 15, 25, 35, 45, 55 °C. The experimental results demonstrate that the viscosity increases by increasing the nanoparticles volume fraction and decreases by increasing temperature. Based on the experimental data the maximum and minimum nanofluid viscosity enhancements, when the volume fraction increases from 0.05 to 1, are 35.52% and 12.92% for constant temperatures of 55 and 15 °C, respectively. The higher viscosity of oil engine in higher temperatures is an advantage, thus this result is important. The measured nanofluid viscosity magnitudes in various shear rates show that this hybrid nanofluid is Newtonian. An ANN model has been employed to predict the viscosity of the CuO-MWCNTs/SAE 5w-50 hybrid nanofluid and the results showed that the ANN can estimate the viscosity efficiently and accurately. Eventually, for viscosity estimation a new temperature and volume fraction based third-degree polynomial empirical model has been developed. The comparison shows that this model is in good agreement with the experimental data.

Convergence of the Critical Cooling Rate for Protoplanetary Disk Fragmentation Achieved: The Key Role of Numerical Dissipation of Angular Momentum

NASA Astrophysics Data System (ADS)

Deng, Hongping; Mayer, Lucio; Meru, Farzana

2017-09-01

We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics (SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code. Our aim is to understand the cause of the nonconvergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementations and compare them with MFM, which does not employ any artificial viscosity. With MFM we demonstrate convergence of the critical cooling timescale for fragmentation at {β }{crit}≈ 3. Nonconvergence persists in SPH codes. We show how the nonconvergence problem is caused by artificial fragmentation triggered by excessive dissipation of angular momentum in domains with large velocity derivatives. With increased resolution, such domains become more prominent. Vorticity lags behind density, due to numerical viscous dissipation in these regions, promoting collapse with longer cooling times. Such effect is shown to be dominant over the competing tendency of artificial viscosity to diminish with increasing resolution. When the initial conditions are first relaxed for several orbits, the flow is more regular, with lower shear and vorticity in nonaxisymmetric regions, aiding convergence. Yet MFM is the only method that converges exactly. Our findings are of general interest, as numerical dissipation via artificial viscosity or advection errors can also occur in grid-based codes. Indeed, for the FARGO code values of {β }{crit} significantly higher than our converged estimate have been reported in the literature. Finally, we discuss implications for giant planet formation via disk instability.

Artificial fluid properties for large-eddy simulation of compressible turbulent mixing

NASA Astrophysics Data System (ADS)

Cook, Andrew W.

2007-05-01

An alternative methodology is described for large-eddy simulation (LES) of flows involving shocks, turbulence, and mixing. In lieu of filtering the governing equations, it is postulated that the large-scale behavior of a LES fluid, i.e., a fluid with artificial properties, will be similar to that of a real fluid, provided the artificial properties obey certain constraints. The artificial properties consist of modifications to the shear viscosity, bulk viscosity, thermal conductivity, and species diffusivity of a fluid. The modified transport coefficients are designed to damp out high wavenumber modes, close to the resolution limit, without corrupting lower modes. Requisite behavior of the artificial properties is discussed and results are shown for a variety of test problems, each designed to exercise different aspects of the models. When combined with a tenth-order compact scheme, the overall method exhibits excellent resolution characteristics for turbulent mixing, while capturing shocks and material interfaces in a crisp fashion.

Modified Method of Adaptive Artificial Viscosity for Solution of Gas Dynamics Problems on Parallel Computer Systems

NASA Astrophysics Data System (ADS)

Popov, Igor; Sukov, Sergey

2018-02-01

A modification of the adaptive artificial viscosity (AAV) method is considered. This modification is based on one stage time approximation and is adopted to calculation of gasdynamics problems on unstructured grids with an arbitrary type of grid elements. The proposed numerical method has simplified logic, better performance and parallel efficiency compared to the implementation of the original AAV method. Computer experiments evidence the robustness and convergence of the method to difference solution.

High-order centered difference methods with sharp shock resolution

NASA Technical Reports Server (NTRS)

Gustafsson, Bertil; Olsson, Pelle

1994-01-01

In this paper we consider high-order centered finite difference approximations of hyperbolic conservation laws. We propose different ways of adding artificial viscosity to obtain sharp shock resolution. For the Riemann problem we give simple explicit formulas for obtaining stationary one and two-point shocks. This can be done for any order of accuracy. It is shown that the addition of artificial viscosity is equivalent to ensuring the Lax k-shock condition. We also show numerical experiments that verify the theoretical results.

How Artificial Should the Treatment of a Plasma's Viscosity Be?

NASA Astrophysics Data System (ADS)

Whitney, K. G.; Velikovich, A. L.; Thornhill, J. W.; Davis, J.

1999-11-01

Electron viscosity dominates over ion viscosity and is important in describing the generation of shock fronts in highly ionizable plasmas. The sizes of shock front jumps in electron and ion temperature are determined from the magnitudes of the heat flow vector and pressure tensor, which, in turn, acquire non-negligible nonlinear contributions from the temperature and density gradients when these gradients are large. Thus, a consistent treatment of steep gradient formation in plasmas must come from investigations that include the effects of these nonlinear contributions to heat and momentum transport. Coefficients for each of five nonlinear contributions to the pressure tensor for an (r,z) Z-pinch geometry are presented and discussed in this talk. Hydrodynamic code calculations generally are not designed to provide a testbed for directly evaluating the kinetic energy dissipation that occurs at shock fronts; therefore, the strength of these nonlinear pressure tensor terms will be estimated by post-processing a Z-pinch hydrodynamics calculation and a steady-state planar shock wave calculation.

Grid Effect on Spherical Shallow Water Jets Using Continuous and Discontinuous Galerkin Methods

DTIC Science & Technology

2013-01-01

The high-order Legendre -Gauss-Lobatto (LGL) points are added to the linear grid by projecting the linear elements onto the auxiliary gnomonic space...mapping, the triangles are subdivided into smaller ones by a Lagrange polynomial of order nI . The number of quadrilateral elements and grid points of...of the acceleration of gravity and the vertical height of the fluid), ν∇2 is the artificial viscosity term of viscous coefficient ν = 1× 105 m2 s−1

GASOLINE: Smoothed Particle Hydrodynamics (SPH) code

NASA Astrophysics Data System (ADS)

N-Body Shop

2017-10-01

Gasoline solves the equations of gravity and hydrodynamics in astrophysical problems, including simulations of planets, stars, and galaxies. It uses an SPH method that features correct mixing behavior in multiphase fluids and minimal artificial viscosity. This method is identical to the SPH method used in the ChaNGa code (ascl:1105.005), allowing users to extend results to problems requiring >100,000 cores. Gasoline uses a fast, memory-efficient O(N log N) KD-Tree to solve Poisson's Equation for gravity and avoids artificial viscosity in non-shocking compressive flows.

Swallowing Mechanics Associated with Artificial Airways, Bolus Properties, and Penetration-Aspiration Status in Trauma Patients

ERIC Educational Resources Information Center

Dietsch, Angela M.; Rowley, Christopher B.; Solomon, Nancy Pearl; Pearson, William G., Jr.

2017-01-01

Purpose: Artificial airway procedures such as intubation and Sare common in the treatment of traumatic injuries, and bolus modifications may be implemented to help manage swallowing disorders. This study assessed artificial airway status, bolus properties (volume and viscosity), and the occurrence of laryngeal penetration and/or aspiration in…

Vanishing Viscosity Approach to the Compressible Euler Equations for Transonic Nozzle and Spherically Symmetric Flows

NASA Astrophysics Data System (ADS)

Chen, Gui-Qiang G.; Schrecker, Matthew R. I.

2018-04-01

We are concerned with globally defined entropy solutions to the Euler equations for compressible fluid flows in transonic nozzles with general cross-sectional areas. Such nozzles include the de Laval nozzles and other more general nozzles whose cross-sectional area functions are allowed at the nozzle ends to be either zero (closed ends) or infinity (unbounded ends). To achieve this, in this paper, we develop a vanishing viscosity method to construct globally defined approximate solutions and then establish essential uniform estimates in weighted L p norms for the whole range of physical adiabatic exponents γ\\in (1, ∞) , so that the viscosity approximate solutions satisfy the general L p compensated compactness framework. The viscosity method is designed to incorporate artificial viscosity terms with the natural Dirichlet boundary conditions to ensure the uniform estimates. Then such estimates lead to both the convergence of the approximate solutions and the existence theory of globally defined finite-energy entropy solutions to the Euler equations for transonic flows that may have different end-states in the class of nozzles with general cross-sectional areas for all γ\\in (1, ∞) . The approach and techniques developed here apply to other problems with similar difficulties. In particular, we successfully apply them to construct globally defined spherically symmetric entropy solutions to the Euler equations for all γ\\in (1, ∞).

Experimental investigation, model development and sensitivity analysis of rheological behavior of ZnO/10W40 nano-lubricants for automotive applications

NASA Astrophysics Data System (ADS)

Hemmat Esfe, Mohammad; Saedodin, Seyfolah; Rejvani, Mousa; Shahram, Jalal

2017-06-01

In the present study, rheological behavior of ZnO/10W40 nano-lubricant is investigated by an experimental approach. Firstly, ZnO nanoparticles of 10-30 nm were dispersed in 10W40 engine oil with solid volume fractions of 0.25-2%, then the viscosity of the composed nano-lubricant was measured in temperature ranges of 5-55 °C and in various shear rates. From analyzing the results, it was revealed that both of the base oil and nano-lubricants are non-Newtonian fluids which exhibit shear thinning behavior. Sensitivity of viscosity to the solid volume fraction enhancement was calculated by a new correlation which was proposed in terms of solid volume fraction and temperature. In order to attain an accurate model by which experimental data are predicted, an artificial neural network (ANN) with a hidden layer and 5 neurons was designed. This model was considerably accurate in predicting experimental data of dynamic viscosity as R-squared and average absolute relative deviation (AARD %) were respectively 0.9999 and 0.0502.

Artificial viscosity to cure the carbuncle phenomenon: The three-dimensional case

NASA Astrophysics Data System (ADS)

Rodionov, Alexander V.

2018-05-01

The carbuncle phenomenon (also known as the shock instability) has remained a serious computational challenge since it was first noticed and described [1,2]. In [3] the author presented a summary on this subject and proposed a new technique for curing the problem. Its idea is to introduce some dissipation in the form of right-hand sides of the Navier-Stokes equations into the basic method of solving Euler equations; in so doing, the molecular viscosity coefficient is replaced by the artificial viscosity coefficient. The new cure for the carbuncle flaw was tested and tuned for the case of using first-order schemes in two-dimensional simulations. Its efficiency was demonstrated on several well-known test problems. In this paper we extend the technique of [3] to the case of three-dimensional simulations.

The Overshoot Phenomenon in Geodynamics Codes

NASA Astrophysics Data System (ADS)

Kommu, R. K.; Heien, E. M.; Kellogg, L. H.; Bangerth, W.; Heister, T.; Studley, E. H.

2013-12-01

The overshoot phenomenon is a common occurrence in numerical software when a continuous function on a finite dimensional discretized space is used to approximate a discontinuous jump, in temperature and material concentration, for example. The resulting solution overshoots, and undershoots, the discontinuous jump. Numerical simulations play an extremely important role in mantle convection research. This is both due to the strong temperature and stress dependence of viscosity and also due to the inaccessibility of deep earth. Under these circumstances, it is essential that mantle convection simulations be extremely accurate and reliable. CitcomS and ASPECT are two finite element based mantle convection simulations developed and maintained by the Computational Infrastructure for Geodynamics. CitcomS is a finite element based mantle convection code that is designed to run on multiple high-performance computing platforms. ASPECT, an adaptive mesh refinement (AMR) code built on the Deal.II library, is also a finite element based mantle convection code that scales well on various HPC platforms. CitcomS and ASPECT both exhibit the overshoot phenomenon. One attempt at controlling the overshoot uses the Entropy Viscosity method, which introduces an artificial diffusion term in the energy equation of mantle convection. This artificial diffusion term is small where the temperature field is smooth. We present results from CitcomS and ASPECT that quantify the effect of the Entropy Viscosity method in reducing the overshoot phenomenon. In the discontinuous Galerkin (DG) finite element method, the test functions used in the method are continuous within each element but are discontinuous across inter-element boundaries. The solution space in the DG method is discontinuous. FEniCS is a collection of free software tools that automate the solution of differential equations using finite element methods. In this work we also present results from a finite element mantle convection simulation implemented in FEniCS that investigates the effect of using DG elements in reducing the overshoot problem.

Spreading properties of cosmetic emollients: Use of synthetic skin surface to elucidate structural effect.

PubMed

Douguet, Marine; Picard, Céline; Savary, Géraldine; Merlaud, Fabien; Loubat-Bouleuc, Nathalie; Grisel, Michel

2017-06-01

The study focuses on the impact of structural and physicochemical properties of emollients on their spreadability. Fifty-three emollients, among which esters, silicones, vegetable and mineral oils, have been characterized. Their viscosity, surface tension, density and spreadability have been measured. Vitro-skin ® , an artificial skin substitute, was used as an artificial porous substrate to measure spreadability. Two different methods have been selected to characterize spreadability, namely contact angle and spreading value. Dynamic contact angle measurements showed that emollient spreadability is first governed by spontaneous spreading and that, in a second phase, absorption and migration into the porous substrate becomes the driver of the extension of the spreading area. Statistical analysis of physicochemical and spreading value data revealed that viscosity has a major impact on the spreading behavior of emollients whatever their chemical type. A special emphasis was placed on the ester family in which chemical diversity is very wide. The results highlighted a difference between "high viscosity esters" for which viscosity is the main factor impacting spreadability and "low viscosity esters" for which structural variations (mono/diester, saturated/unsaturated chain, linear/branched chain) have to be considered in addition to viscosity. Linear regressions were used to express spreading value as a function of viscosity for each of the four emollient families tested (esters, silicones, vegetable and mineral oils). These regressions allowed the development of reliable predictive models as a powerful tool for formulators to forecast spreadability of emollients. Copyright © 2017 Elsevier B.V. All rights reserved.

Computation of turbulent pipe and duct flow using third order upwind scheme

NASA Technical Reports Server (NTRS)

Kawamura, T.

1986-01-01

The fully developed turbulence in a circular pipe and in a square duct is simulated directly without using turbulence models in the Navier-Stokes equations. The utilized method employs a third-order upwind scheme for the approximation to the nonlinear term and the second-order Adams-Bashforth method for the time derivative in the Navier-Stokes equation. The computational results appear to capture the large-scale turbulent structures at least qualitatively. The significance of the artificial viscosity inherent in the present scheme is discussed.

Simulation of separated flow past a bluff body using Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Ghia, K. N.; Ghia, U.; Osswald, G. A.; Liu, C. A.

1987-01-01

Two-dimensional flow past a bluff body is presently simulated on the basis of an analysis that employs the incompressible, unsteady Navier-Stokes equations in terms of vorticity and stream function. The fully implicit, time-marching, alternating-direction, implicit-block Gaussian elimination used is a direct method with second-order spatial accuracy; this allows it to avoid the introduction of any artificial viscosity. Attention is given to the simulation of flow past a circular cylinder with and without symmetry, requiring the use of either the half or the full cylinder, respectively.

Experiment and Artificial Neural Network Prediction of Thermal Conductivity and Viscosity for Alumina-Water Nanofluids

PubMed Central

Zhao, Ningbo; Li, Zhiming

2017-01-01

To effectively predict the thermal conductivity and viscosity of alumina (Al2O3)-water nanofluids, an artificial neural network (ANN) approach was investigated in the present study. Firstly, using a two-step method, four Al2O3-water nanofluids were prepared respectively by dispersing different volume fractions (1.31%, 2.72%, 4.25%, and 5.92%) of nanoparticles with the average diameter of 30 nm. On this basis, the thermal conductivity and viscosity of the above nanofluids were analyzed experimentally under various temperatures ranging from 296 to 313 K. Then a radial basis function (RBF) neural network was constructed to predict the thermal conductivity and viscosity of Al2O3-water nanofluids as a function of nanoparticle volume fraction and temperature. The experimental results showed that both nanoparticle volume fraction and temperature could enhance the thermal conductivity of Al2O3-water nanofluids. However, the viscosity only depended strongly on Al2O3 nanoparticle volume fraction and was increased slightly by changing temperature. In addition, the comparative analysis revealed that the RBF neural network had an excellent ability to predict the thermal conductivity and viscosity of Al2O3-water nanofluids with the mean absolute percent errors of 0.5177% and 0.5618%, respectively. This demonstrated that the ANN provided an effective way to predict the thermophysical properties of nanofluids with limited experimental data. PMID:28772913

Experiment and Artificial Neural Network Prediction of Thermal Conductivity and Viscosity for Alumina-Water Nanofluids.

PubMed

Zhao, Ningbo; Li, Zhiming

2017-05-19

To effectively predict the thermal conductivity and viscosity of alumina (Al₂O₃)-water nanofluids, an artificial neural network (ANN) approach was investigated in the present study. Firstly, using a two-step method, four Al₂O₃-water nanofluids were prepared respectively by dispersing different volume fractions (1.31%, 2.72%, 4.25%, and 5.92%) of nanoparticles with the average diameter of 30 nm. On this basis, the thermal conductivity and viscosity of the above nanofluids were analyzed experimentally under various temperatures ranging from 296 to 313 K. Then a radial basis function (RBF) neural network was constructed to predict the thermal conductivity and viscosity of Al₂O₃-water nanofluids as a function of nanoparticle volume fraction and temperature. The experimental results showed that both nanoparticle volume fraction and temperature could enhance the thermal conductivity of Al₂O₃-water nanofluids. However, the viscosity only depended strongly on Al₂O₃ nanoparticle volume fraction and was increased slightly by changing temperature. In addition, the comparative analysis revealed that the RBF neural network had an excellent ability to predict the thermal conductivity and viscosity of Al₂O₃-water nanofluids with the mean absolute percent errors of 0.5177% and 0.5618%, respectively. This demonstrated that the ANN provided an effective way to predict the thermophysical properties of nanofluids with limited experimental data.

Sucralose Destabilization of Protein Structure.

PubMed

Chen, Lee; Shukla, Nimesh; Cho, Inha; Cohn, Erin; Taylor, Erika A; Othon, Christina M

2015-04-16

Sucralose is a commonly employed artificial sweetener that behaves very differently than its natural disaccharide counterpart, sucrose, in terms of its interaction with biomolecules. The presence of sucralose in solution is found to destabilize the native structure of two model protein systems: the globular protein bovine serum albumin and an enzyme staphylococcal nuclease. The melting temperature of these proteins decreases as a linear function of sucralose concentration. We correlate this destabilization to the increased polarity of the molecule. The strongly polar nature is manifested as a large dielectric friction exerted on the excited-state rotational diffusion of tryptophan using time-resolved fluorescence anisotropy. Tryptophan exhibits rotational diffusion proportional to the measured bulk viscosity for sucrose solutions over a wide range of concentrations, consistent with a Stokes-Einstein model. For sucralose solutions, however, the diffusion is dependent on the concentration, strongly diverging from the viscosity predictions, and results in heterogeneous rotational diffusion.

Influence of numerical dissipation in computing supersonic vortex-dominated flows

NASA Technical Reports Server (NTRS)

Kandil, O. A.; Chuang, A.

1986-01-01

Steady supersonic vortex-dominated flows are solved using the unsteady Euler equations for conical and three-dimensional flows around sharp- and round-edged delta wings. The computational method is a finite-volume scheme which uses a four-stage Runge-Kutta time stepping with explicit second- and fourth-order dissipation terms. The grid is generated by a modified Joukowski transformation. The steady flow solution is obtained through time-stepping with initial conditions corresponding to the freestream conditions, and the bow shock is captured as a part of the solution. The scheme is applied to flat-plate and elliptic-section wings with a leading edge sweep of 70 deg at an angle of attack of 10 deg and a freestream Mach number of 2.0. Three grid sizes of 29 x 39, 65 x 65 and 100 x 100 have been used. The results for sharp-edged wings show that they are consistent with all grid sizes and variation of the artificial viscosity coefficients. The results for round-edged wings show that separated and attached flow solutions can be obtained by varying the artificial viscosity coefficients. They also show that the solutions are independent of the way time stepping is done. Local time-stepping and global minimum time-steeping produce same solutions.

A comparison of performance of several artificial intelligence methods for predicting the dynamic viscosity of TiO2/SAE 50 nano-lubricant

NASA Astrophysics Data System (ADS)

Hemmat Esfe, Mohammad; Tatar, Afshin; Ahangar, Mohammad Reza Hassani; Rostamian, Hossein

2018-02-01

Since the conventional thermal fluids such as water, oil, and ethylene glycol have poor thermal properties, the tiny solid particles are added to these fluids to increase their heat transfer improvement. As viscosity determines the rheological behavior of a fluid, studying the parameters affecting the viscosity is crucial. Since the experimental measurement of viscosity is expensive and time consuming, predicting this parameter is the apt method. In this work, three artificial intelligence methods containing Genetic Algorithm-Radial Basis Function Neural Networks (GA-RBF), Least Square Support Vector Machine (LS-SVM) and Gene Expression Programming (GEP) were applied to predict the viscosity of TiO2/SAE 50 nano-lubricant with Non-Newtonian power-law behavior using experimental data. The correlation factor (R2), Average Absolute Relative Deviation (AARD), Root Mean Square Error (RMSE), and Margin of Deviation were employed to investigate the accuracy of the proposed models. RMSE values of 0.58, 1.28, and 6.59 and R2 values of 0.99998, 0.99991, and 0.99777 reveal the accuracy of the proposed models for respective GA-RBF, CSA-LSSVM, and GEP methods. Among the developed models, the GA-RBF shows the best accuracy.

Plasma viscosity with mass transport in spherical inertial confinement fusion implosion simulations

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

Vold, E. L.; Molvig, K.; Joglekar, A. S.

2015-11-15

The effects of viscosity and small-scale atomic-level mixing on plasmas in inertial confinement fusion (ICF) currently represent challenges in ICF research. Many current ICF hydrodynamic codes ignore the effects of viscosity though recent research indicates viscosity and mixing by classical transport processes may have a substantial impact on implosion dynamics. We have implemented a Lagrangian hydrodynamic code in one-dimensional spherical geometry with plasma viscosity and mass transport and including a three temperature model for ions, electrons, and radiation treated in a gray radiation diffusion approximation. The code is used to study ICF implosion differences with and without plasma viscosity andmore » to determine the impacts of viscosity on temperature histories and neutron yield. It was found that plasma viscosity has substantial impacts on ICF shock dynamics characterized by shock burn timing, maximum burn temperatures, convergence ratio, and time history of neutron production rates. Plasma viscosity reduces the need for artificial viscosity to maintain numerical stability in the Lagrangian formulation and also modifies the flux-limiting needed for electron thermal conduction.« less

Numerical dissipation vs. subgrid-scale modelling for large eddy simulation

NASA Astrophysics Data System (ADS)

Dairay, Thibault; Lamballais, Eric; Laizet, Sylvain; Vassilicos, John Christos

2017-05-01

This study presents an alternative way to perform large eddy simulation based on a targeted numerical dissipation introduced by the discretization of the viscous term. It is shown that this regularisation technique is equivalent to the use of spectral vanishing viscosity. The flexibility of the method ensures high-order accuracy while controlling the level and spectral features of this purely numerical viscosity. A Pao-like spectral closure based on physical arguments is used to scale this numerical viscosity a priori. It is shown that this way of approaching large eddy simulation is more efficient and accurate than the use of the very popular Smagorinsky model in standard as well as in dynamic version. The main strength of being able to correctly calibrate numerical dissipation is the possibility to regularise the solution at the mesh scale. Thanks to this property, it is shown that the solution can be seen as numerically converged. Conversely, the two versions of the Smagorinsky model are found unable to ensure regularisation while showing a strong sensitivity to numerical errors. The originality of the present approach is that it can be viewed as implicit large eddy simulation, in the sense that the numerical error is the source of artificial dissipation, but also as explicit subgrid-scale modelling, because of the equivalence with spectral viscosity prescribed on a physical basis.

A computational fluid dynamics simulation of the hypersonic flight of the Pegasus(TM) vehicle using an artificial viscosity model and a nonlinear filtering method. M.S. Thesis

NASA Technical Reports Server (NTRS)

Mendoza, John Cadiz

1995-01-01

The computational fluid dynamics code, PARC3D, is tested to see if its use of non-physical artificial dissipation affects the accuracy of its results. This is accomplished by simulating a shock-laminar boundary layer interaction and several hypersonic flight conditions of the Pegasus(TM) launch vehicle using full artificial dissipation, low artificial dissipation, and the Engquist filter. Before the filter is applied to the PARC3D code, it is validated in one-dimensional and two-dimensional form in a MacCormack scheme against the Riemann and convergent duct problem. For this explicit scheme, the filter shows great improvements in accuracy and computational time as opposed to the nonfiltered solutions. However, for the implicit PARC3D code it is found that the best estimate of the Pegasus experimental heat fluxes and surface pressures is the simulation utilizing low artificial dissipation and no filter. The filter does improve accuracy over the artificially dissipative case but at a computational expense greater than that achieved by the low artificial dissipation case which has no computational time penalty and shows better results. For the shock-boundary layer simulation, the filter does well in terms of accuracy for a strong impingement shock but not as well for weaker shock strengths. Furthermore, for the latter problem the filter reduces the required computational time to convergence by 18.7 percent.

Evaluation of sensor arrays for engine oils using artificial oil alteration

NASA Astrophysics Data System (ADS)

Sen, Sedat; Schneidhofer, Christoph; Dörr, Nicole; Vellekoop, Michael J.

2011-06-01

With respect to varying operation conditions, only sensors directly installed in the engine can detect the current oil condition hence enabling to get the right time for the oil change. Usually, only one parameter is not sufficient to obtain reliable information about the current oil condition. For this reason, appropriate sensor principles were evaluated for the design of sensor arrays for the measurement of critical lubricant parameters. In this contribution, we report on the development of a sensor array for engine oils using laboratory analyses of used engine oils for the correlation with sensor signals. The sensor array comprises the measurement of conductivity, permittivity, viscosity and temperature as well as oil corrosiveness as a consequence of acidification of the lubricant. As a key method, rapid evaluation of the sensors was done by short term simulation of entire oil change intervals based on artificial oil alteration. Thereby, the compatibility of the sensor array to the lubricant and the oil deterioration during the artificial alteration process was observed by the sensors and confirmed by additional laboratory analyses of oil samples take.

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 3: Programmer's reference

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This User's Guide describes the program's features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Biofluid lubrication for artificial joints

NASA Astrophysics Data System (ADS)

Pendleton, Alice Mae

This research investigated biofluid lubrication related to artificial joints using tribological and rheological approaches. Biofluids studied here represent two categories of fluids, base fluids and nanostructured biofluids. Base fluids were studied through comparison of synthetic fluids (simulated body fluid and hyaluronic acid) as well as natural biofluids (from dogs, horses, and humans) in terms of viscosity and fluid shear stress. The nano-structured biofluids were formed using molecules having well-defined shapes. Understanding nano-structured biofluids leads to new ways of design and synthesis of biofluids that are beneficial for artificial joint performance. Experimental approaches were utilized in the present research. This includes basic analysis of biofluids' property, such as viscosity, fluid shear stress, and shear rate using rheological experiments. Tribological investigation and surface characterization were conducted in order to understand effects of molecular and nanostructures on fluid lubrication. Workpiece surface structure and wear mechanisms were investigated using a scanning electron microscope and a transmission electron microscope. The surface topography was examined using a profilometer. The results demonstrated that with the adding of solid additives, such as crown ether or fullerene acted as rough as the other solids in the 3-body wear systems. In addition, the fullerene supplied low friction and low wear, which designates the lubrication purpose of this particular particle system. This dissertation is constructed of six chapters. The first chapter is an introduction to body fluids, as mentioned earlier. After Chapter II, it examines the motivation and approach of the present research, Chapter III discusses the experimental approaches, including materials, experimental setup, and conditions. In Chapter IV, lubrication properties of various fluids are discussed. The tribological properties and performance nanostructured biofluids are discussed in Chapter V, followed by summary and conclusions in Chapter VI.

Viscous regularization of the full set of nonequilibrium-diffusion Grey Radiation-Hydrodynamic equations

DOE PAGES

Delchini, Marc O.; Ragusa, Jean C.; Ferguson, Jim

2017-02-17

A viscous regularization technique, based on the local entropy residual, was proposed by Delchini et al. (2015) to stabilize the nonequilibrium-diffusion Grey Radiation-Hydrodynamic equations using an artificial viscosity technique. This viscous regularization is modulated by the local entropy production and is consistent with the entropy minimum principle. However, Delchini et al. (2015) only based their work on the hyperbolic parts of the Grey Radiation-Hydrodynamic equations and thus omitted the relaxation and diffusion terms present in the material energy and radiation energy equations. Here in this paper, we extend the theoretical grounds for the method and derive an entropy minimum principlemore » for the full set of nonequilibrium-diffusion Grey Radiation-Hydrodynamic equations. This further strengthens the applicability of the entropy viscosity method as a stabilization technique for radiation-hydrodynamic shock simulations. Radiative shock calculations using constant and temperature-dependent opacities are compared against semi-analytical reference solutions, and we present a procedure to perform spatial convergence studies of such simulations.« less

Strawberry-flavored yogurts and whey beverages: What is the sensory profile of the ideal product?

PubMed

Janiaski, D R; Pimentel, T C; Cruz, A G; Prudencio, S H

2016-07-01

This study aimed to evaluate the sensory profile and Brazilian consumers' liking of strawberry-flavored yogurts and whey beverages (fermented or nonfermented) with different fat contents that were sweetened with sugar or nonsugar sweeteners. We also determined the influence of sensory attributes on consumer preferences and the profile of the ideal product. Nonfermented whey beverages (NFWB) and "light" yogurt were less liked. The NFWB were less acidic, less viscous, and with lower smoothness of mouthcoating, sweeter and with a more intense artificial strawberry aroma (ASA) than the fermented products. Low-fat yogurts were more liked, more viscous, and had higher smoothness of mouthcoating than nonfat yogurts. Fermented-whey beverages were as liked as yogurts. Viscosity and smoothness of mouthcoating positively influenced consumer liking. The ideal product had higher levels of brightness, artificial strawberry taste, artificial strawberry aroma, and sweet taste; intermediate smoothness of mouthcoating, color, and viscosity; and low particles, acid taste, and aroma. Copyright © 2016 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.

Lectures series in computational fluid dynamics

NASA Technical Reports Server (NTRS)

Thompson, Kevin W.

1987-01-01

The lecture notes cover the basic principles of computational fluid dynamics (CFD). They are oriented more toward practical applications than theory, and are intended to serve as a unified source for basic material in the CFD field as well as an introduction to more specialized topics in artificial viscosity and boundary conditions. Each chapter in the test is associated with a videotaped lecture. The basic properties of conservation laws, wave equations, and shock waves are described. The duality of the conservation law and wave representations is investigated, and shock waves are examined in some detail. Finite difference techniques are introduced for the solution of wave equations and conservation laws. Stability analysis for finite difference approximations are presented. A consistent description of artificial viscosity methods are provided. Finally, the problem of nonreflecting boundary conditions are treated.

Upwind methods for the Baer–Nunziato equations and higher-order reconstruction using artificial viscosity

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

Fraysse, F., E-mail: francois.fraysse@rs2n.eu; E. T. S. de Ingeniería Aeronáutica y del Espacio, Universidad Politécnica de Madrid, Madrid; Redondo, C.

This article is devoted to the numerical discretisation of the hyperbolic two-phase flow model of Baer and Nunziato. A special attention is paid on the discretisation of intercell flux functions in the framework of Finite Volume and Discontinuous Galerkin approaches, where care has to be taken to efficiently approximate the non-conservative products inherent to the model equations. Various upwind approximate Riemann solvers have been tested on a bench of discontinuous test cases. New discretisation schemes are proposed in a Discontinuous Galerkin framework following the criterion of Abgrall and the path-conservative formalism. A stabilisation technique based on artificial viscosity is appliedmore » to the high-order Discontinuous Galerkin method and compared against classical TVD-MUSCL Finite Volume flux reconstruction.« less

The effects of rotational flow, viscosity, thickness, and shape on transonic flutter dip phenomena

NASA Technical Reports Server (NTRS)

Reddy, T. S. R.; Srivastava, Rakesh; Kaza, Krishna Rao V.

1988-01-01

The transonic flutter dip phenomena on thin airfoils, which are employed for propfan blades, is investigated using an integrated Euler/Navier-Stokes code and a two degrees of freedom typical section structural model. As a part of the code validation, the flutter characteristics of the NACA 64A010 airfoil are also investigated. In addition, the effects of artificial dissipation models, rotational flow, initial conditions, mean angle of attack, viscosity, airfoil thickness and shape on flutter are investigated. The results obtained with a Euler code for the NACA 64A010 airfoil are in reasonable agreement with published results obtained by using transonic small disturbance and Euler codes. The two artificial dissipation models, one based on the local pressure gradient scaled by a common factor and the other based on the local pressure gradient scaled by a spectral radius, predicted the same flutter speeds except in the recovery region for the case studied. The effects of rotational flow, initial conditions, mean angle of attack, and viscosity for the Reynold's number studied seem to be negligible or small on the minima of the flutter dip.

Shock Capturing with PDE-Based Artificial Viscosity for an Adaptive, Higher-Order Discontinuous Galerkin Finite Element Method

DTIC Science & Technology

2008-06-01

Geometry Interpolation The function space , VpH , consists of discontinuous, piecewise-polynomials. This work used a polynomial basis for VpH such...between a piecewise-constant and smooth variation of viscosity in both a one- dimensional and multi- dimensional setting. Before continuing with the ...inviscid, transonic flow past a NACA 0012 at zero angle of attack and freestream Mach number of M∞ = 0.95. The

Heat Generation in Axial and Centrifugal Flow Left Ventricular Assist Devices.

PubMed

Yost, Gardner; Joseph, Christine Rachel; Royston, Thomas; Tatooles, Antone; Bhat, Geetha

Despite increasing use of left ventricular assist devices (LVADs) as a surgical treatment for advanced heart failure in an era of improved outcomes with LVAD support, the mechanical interactions between these pumps and the cardiovascular system are not completely understood. We utilized an in vitro mock circulatory loop to analyze the heat production incurred by operation of an axial flow and centrifugal flow LVAD. A HeartMate II and a HeartWare HVAD were connected to an abbreviated flow loop and were implanted in a viscoelastic gel. Temperature was measured at the surface of each LVAD. Device speed and fluid viscosity were altered and, in the HeartMate II, as artificial thrombi were attached to the inflow stator, impeller, and outflow stator. The surface temperatures of both LVADs increased in all trials and reached a plateau within 80 minutes of flow initiation. Rate of heat generation and maximum system temperature were greater when speed was increased, when viscosity was increased, and when artificial thrombi were attached to the HeartMate II impeller. Normal operation of these two widely utilized LVADs results in appreciable heat generation in vitro. Increased pump loading resulted in more rapid heat generation, which was particularly severe when a large thrombus was attached to the impeller of the HeartMate II. While heat accumulation in vivo is likely minimized by greater dissipation in the blood and soft tissues, focal temperature gains with the pump housing of these two devices during long-term operation may have negative hematological consequences.

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 2: User's guide

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the User's Guide, and describes the program's features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Experimental study of the effects of surface mucus viscosity on the glottic cycle.

PubMed

Ayache, Stéphane; Ouaknine, Maurice; Dejonkere, Philippe; Prindere, Pierre; Giovanni, Antoine

2004-03-01

Numerous clinical findings indicate that viscosity of laryngeal mucosa is a crucial factor in glottal perfomance. Experience using experimental test benches has shown the importance of humidifying air stream used to induce vibration in excised larynges. Nevertheless, there is a lack of knowledge particularly regarding the physicochemical properties of laryngeal mucus. The purpose of this study was to research vocal fold vibration in excised larynges using artificial mucus of precisely known viscosity. Eight freshly harvested porcine larynges were examined. Parameters measured were Fo and vocal fold contact time. Measurements were performed under three conditions: basal (no fluid application on vocal cord surface), after application of a fluid of 60cP viscosity (Visc60), and after application of a fluid of 100cP viscosity (Visc100). Electroglottographic measurements were performed at two different times for each condition: 1 s after airflow onset (T1) and 6 seconds after airflow onset (T2). Statistical analysis consisted of comparing data obtained under each condition at T1 and T2. The results showed a significant decrease in Fo after application of Visc60 and Visc100 fluids and a decrease in Fo at T2. Closure time was significantly higher under Visc60 conditions and under Visc100 conditions than under basal conditions. Application of artificial mucus to the mucosa of the vocal folds lowered vibratory frequency and prolonged the contact phase. Our interpretation of this data is that the presence of mucus on the surface of the vocal folds generated superficial tension and caused adhesion, which is a source of nonlinearity in vocal vibration.

Improved High Resolution Models of Subduction Dynamics: Use of transversely isotropic viscosity with a free-surface

NASA Astrophysics Data System (ADS)

Liu, X.; Gurnis, M.; Stadler, G.; Rudi, J.; Ratnaswamy, V.; Ghattas, O.

2017-12-01

Dynamic topography, or uncompensated topography, is controlled by internal dynamics, and provide constraints on the buoyancy structure and rheological parameters in the mantle. Compared with other surface manifestations such as the geoid, dynamic topography is very sensitive to shallower and more regional mantle structure. For example, the significant dynamic topography above the subduction zone potentially provides a rich mine for inferring the rheological and mechanical properties such as plate coupling, flow, and lateral viscosity variations, all critical in plate tectonics. However, employing subduction zone topography in the inversion study requires that we have a better understanding of the topography from forward models, especially the influence of the viscosity formulation, numerical resolution, and other factors. One common approach to formulating a fault between the subducted slab and the overriding plates in viscous flow models assumes a thin weak zone. However, due to the large lateral variation in viscosity, topography from free-slip numerical models typically has artificially large magnitude as well as high-frequency undulations over subduction zone, which adds to the difficulty in making comparisons between model results and observations. In this study, we formulate a weak zone with the transversely isotropic viscosity (TI) where the tangential viscosity is much smaller than the viscosity in the normal direction. Similar with isotropic weak zone models, TI models effectively decouple subducted slabs from the overriding plates. However, we find that the topography in TI models is largely reduced compared with that in weak zone models assuming an isotropic viscosity. Moreover, the artificial `tooth paste' squeezing effect observed in isotropic weak zone models vanishes in TI models, although the difference becomes less significant when the dip angle is small. We also implement a free-surface condition in our numerical models, which has a smoothing effect on the topography. With the improved model configuration, we can use the adjoint inversion method in a high-resolution model and employ topography in addition to other observables such as the plate motion to infer critical mechanical and rheological parameters in the subduction zone.

Multigrid calculation of three-dimensional viscous cascade flows

NASA Technical Reports Server (NTRS)

Arnone, A.; Liou, M.-S.; Povinelli, L. A.

1991-01-01

A 3-D code for viscous cascade flow prediction was developed. The space discretization uses a cell-centered scheme with eigenvalue scaling to weigh the artificial dissipation terms. Computational efficiency of a four stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full multigrid method. The Baldwin-Lomax eddy viscosity model is used for turbulence closure. A zonal, nonperiodic grid is used to minimize mesh distortion in and downstream of the throat region. Applications are presented for an annular vane with and without end wall contouring, and for a large scale linear cascade. The calculation is validated by comparing with experiments and by studying grid dependency.

Multigrid calculation of three-dimensional viscous cascade flows

NASA Technical Reports Server (NTRS)

Arnone, A.; Liou, M.-S.; Povinelli, L. A.

1991-01-01

A three-dimensional code for viscous cascade flow prediction has been developed. The space discretization uses a cell-centered scheme with eigenvalue scaling to weigh the artificial dissipation terms. Computational efficiency of a four-stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full-multigrid method. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. A zonal, nonperiodic grid is used to minimize mesh distortion in and downstream of the throat region. Applications are presented for an annular vane with and without end wall contouring, and for a large-scale linear cascade. The calculation is validated by comparing with experiments and by studying grid dependency.

Plasma viscosity with mass transport in spherical inertial confinement fusion implosion simulations

DOE PAGES

Vold, Erik Lehman; Joglekar, Archis S.; Ortega, Mario I.; ...

2015-11-20

The effects of viscosity and small-scale atomic-level mixing on plasmas in inertial confinement fusion(ICF) currently represent challenges in ICF research. Many current ICF hydrodynamic codes ignore the effects of viscosity though recent research indicates viscosity and mixing by classical transport processes may have a substantial impact on implosion dynamics. In this paper, we have implemented a Lagrangian hydrodynamic code in one-dimensional spherical geometry with plasmaviscosity and mass transport and including a three temperature model for ions, electrons, and radiation treated in a gray radiation diffusion approximation. The code is used to study ICF implosion differences with and without plasmaviscosity andmore » to determine the impacts of viscosity on temperature histories and neutron yield. It was found that plasmaviscosity has substantial impacts on ICF shock dynamics characterized by shock burn timing, maximum burn temperatures, convergence ratio, and time history of neutron production rates. Finally, plasmaviscosity reduces the need for artificial viscosity to maintain numerical stability in the Lagrangian formulation and also modifies the flux-limiting needed for electron thermal conduction.« less

Smoothed Particle Hydrodynamics Simulations of Ultrarelativistic Shocks with Artificial Viscosity

NASA Astrophysics Data System (ADS)

Siegler, S.; Riffert, H.

2000-03-01

We present a fully Lagrangian conservation form of the general relativistic hydrodynamic equations for perfect fluids with artificial viscosity in a given arbitrary background spacetime. This conservation formulation is achieved by choosing suitable Lagrangian time evolution variables, from which the generic fluid variables of rest-mass density, 3-velocity, and thermodynamic pressure have to be determined. We present the corresponding equations for an ideal gas and show the existence and uniqueness of the solution. On the basis of the Lagrangian formulation we have developed a three-dimensional general relativistic smoothed particle hydrodynamics (SPH) code using the standard SPH formalism as known from nonrelativistic fluid dynamics. One-dimensional simulations of a shock tube and a wall shock are presented together with a two-dimensional test calculation of an inclined shock tube. With our method we can model ultrarelativistic fluid flows including shocks with Lorentz factors of even 1000.

A constitutive relation for the viscous flow of an oriented fiber assembly

NASA Technical Reports Server (NTRS)

Pipes, R. B.; Hearle, J. W. S.; Beaussart, A. J.; Sastry, A. M.; Okine, R. K.

1991-01-01

A constitutive relation for an equivalent, homogeneous fluid is developed for the anisotropic viscous flow of an oriented assembly of discontinuous fibers suspended in a viscous fluid. The anisotropic viscous compliance matrix can be expressed in terms of three constants by assuming the equivalent fluid to be incompressible and the microstructure to have axial symmetry (transversely isotropic). By means of a micromechanics analysis, the three terms of the constitutive relation are expressed in terms of the viscosity of the matrix fluid, the fiber aspect ratio, and the fiber volume fraction. A comparison of the viscosity terms reveals that the elongational viscosity in the fiber direction varies as the square of the fiber aspect ratio and a complex function of the fiber volume fraction. Furthermore, the ratio of the axial elongational viscosity to the transverse elongational viscosity and both axial and transverse shear viscosities was shown to be 10 exp 4 - 10 exp 6 for fiber aspect ratio of 100-1000, except at extreme values of the fiber volume fraction.

Artificial neural network based particle size prediction of polymeric nanoparticles.

PubMed

Youshia, John; Ali, Mohamed Ehab; Lamprecht, Alf

2017-10-01

Particle size of nanoparticles and the respective polydispersity are key factors influencing their biopharmaceutical behavior in a large variety of therapeutic applications. Predicting these attributes would skip many preliminary studies usually required to optimize formulations. The aim was to build a mathematical model capable of predicting the particle size of polymeric nanoparticles produced by a pharmaceutical polymer of choice. Polymer properties controlling the particle size were identified as molecular weight, hydrophobicity and surface activity, and were quantified by measuring polymer viscosity, contact angle and interfacial tension, respectively. A model was built using artificial neural network including these properties as input with particle size and polydispersity index as output. The established model successfully predicted particle size of nanoparticles covering a range of 70-400nm prepared from other polymers. The percentage bias for particle prediction was 2%, 4% and 6%, for the training, validation and testing data, respectively. Polymer surface activity was found to have the highest impact on the particle size followed by viscosity and finally hydrophobicity. Results of this study successfully highlighted polymer properties affecting particle size and confirmed the usefulness of artificial neural networks in predicting the particle size and polydispersity of polymeric nanoparticles. Copyright © 2017 Elsevier B.V. All rights reserved.

Dynamics and structures of transitional viscoelastic turbulence in channel flow

NASA Astrophysics Data System (ADS)

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

2017-11-01

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

Numerical Viscosity and the Survival of Gas Giant Protoplanets in Disk Simulations

NASA Astrophysics Data System (ADS)

Pickett, Megan K.; Durisen, Richard H.

2007-01-01

We present three-dimensional hydrodynamic simulations of a gravitationally unstable protoplanetary disk model under the condition of local isothermality. Ordinarily, local isothermality precludes the need for an artificial viscosity (AV) scheme to mediate shocks. Without AV, the disk evolves violently, shredding into dense (although short-lived) clumps. When we introduce our AV treatment in the momentum equation, but without heating due to irreversible compression, our grid-based simulations begin to resemble smoothed particle hydrodynamics (SPH) calculations, where clumps are more likely to survive many orbits. In fact, the standard SPH viscosity appears comparable in strength to the AV that leads to clump longevity in our code. This sensitivity to one numerical parameter suggests extreme caution in interpreting simulations by any code in which long-lived gaseous protoplanetary bodies appear.

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Benson, Thomas J.; Suresh, Ambady

1990-01-01

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Discontinuous Galerkin (DG) Method for solving time dependent convection-diffusion type temperature equation : Demonstration and Comparison with Other Methods in the Mantle Convection Code ASPECT

NASA Astrophysics Data System (ADS)

He, Y.; Puckett, E. G.; Billen, M. I.; Kellogg, L. H.

2016-12-01

For a convection-dominated system, like convection in the Earth's mantle, accurate modeling of the temperature field in terms of the interaction between convective and diffusive processes is one of the most common numerical challenges. In the geodynamics community using Finite Element Method (FEM) with artificial entropy viscosity is a popular approach to resolve this difficulty, but introduce numerical diffusion. The extra artificial viscosity added into the temperature system will not only oversmooth the temperature field where the convective process dominates, but also change the physical properties by increasing the local material conductivity, which will eventually change the local conservation of energy. Accurate modeling of temperature is especially important in the mantle, where material properties are strongly dependent on temperature. In subduction zones, for example, the rheology of the cold sinking slab depends nonlinearly on the temperature, and physical processes such as slab detachment, rollback, and melting all are sensitively dependent on temperature and rheology. Therefore methods that overly smooth the temperature may inaccurately represent the physical processes governing subduction, lithospheric instabilities, plume generation and other aspects of mantle convection. Here we present a method for modeling the temperature field in mantle dynamics simulations using a new solver implemented in the ASPECT software. The new solver for the temperature equation uses a Discontinuous Galerkin (DG) approach, which combines features of both finite element and finite volume methods, and is particularly suitable for problems satisfying the conservation law, and the solution has a large variation locally. Furthermore, we have applied a post-processing technique to insure that the solution satisfies a local discrete maximum principle in order to eliminate the overshoots and undershoots in the temperature locally. To demonstrate the capabilities of this new method we present benchmark results (e.g., falling sphere), and a simple subduction models with kinematic surface boundary condition. To evaluate the trade-offs in computational speed and solution accuracy we present results for the same benchmarks using the Finite Element entropy viscosity method available in ASPECT.

Stability analysis of Eulerian-Lagrangian methods for the one-dimensional shallow-water equations

USGS Publications Warehouse

Casulli, V.; Cheng, R.T.

1990-01-01

In this paper stability and error analyses are discussed for some finite difference methods when applied to the one-dimensional shallow-water equations. Two finite difference formulations, which are based on a combined Eulerian-Lagrangian approach, are discussed. In the first part of this paper the results of numerical analyses for an explicit Eulerian-Lagrangian method (ELM) have shown that the method is unconditionally stable. This method, which is a generalized fixed grid method of characteristics, covers the Courant-Isaacson-Rees method as a special case. Some artificial viscosity is introduced by this scheme. However, because the method is unconditionally stable, the artificial viscosity can be brought under control either by reducing the spatial increment or by increasing the size of time step. The second part of the paper discusses a class of semi-implicit finite difference methods for the one-dimensional shallow-water equations. This method, when the Eulerian-Lagrangian approach is used for the convective terms, is also unconditionally stable and highly accurate for small space increments or large time steps. The semi-implicit methods seem to be more computationally efficient than the explicit ELM; at each time step a single tridiagonal system of linear equations is solved. The combined explicit and implicit ELM is best used in formulating a solution strategy for solving a network of interconnected channels. The explicit ELM is used at channel junctions for each time step. The semi-implicit method is then applied to the interior points in each channel segment. Following this solution strategy, the channel network problem can be reduced to a set of independent one-dimensional open-channel flow problems. Numerical results support properties given by the stability and error analyses. ?? 1990.

Finite-Difference Lattice Boltzmann Scheme for High-Speed Compressible Flow: Two-Dimensional Case

NASA Astrophysics Data System (ADS)

Gan, Yan-Biao; Xu, Ai-Guo; Zhang, Guang-Cai; Zhang, Ping; Zhang, Lei; Li, Ying-Jun

2008-07-01

Lattice Boltzmann (LB) modeling of high-speed compressible flows has long been attempted by various authors. One common weakness of most of previous models is the instability problem when the Mach number of the flow is large. In this paper we present a finite-difference LB model, which works for flows with flexible ratios of specific heats and a wide range of Mach number, from 0 to 30 or higher. Besides the discrete-velocity-model by Watari [Physica A 382 (2007) 502], a modified Lax Wendroff finite difference scheme and an artificial viscosity are introduced. The combination of the finite-difference scheme and the adding of artificial viscosity must find a balance of numerical stability versus accuracy. The proposed model is validated by recovering results of some well-known benchmark tests: shock tubes and shock reflections. The new model may be used to track shock waves and/or to study the non-equilibrium procedure in the transition between the regular and Mach reflections of shock waves, etc.

Lowered pH Alters Decay but Not Speed of Tectorial Membrane Waves

NASA Astrophysics Data System (ADS)

Farrahi, Shirin; Ghaffari, Roozbeh; Freeman, Dennis M.

2011-11-01

Tectorial membrane (TM) traveling waves and mechanical shear impedances were measured in artificial endolymph baths at neutral and acidic pHs. Lowering pH from 7 to 4 significantly decreases the spatial extent of TM waves but has a relatively minor effect on wave speed. At pH 4, the imaginary component of TM shear impedance, which relates to the shear modulus, drops significantly; whereas, the real component, which relates to viscosity, is reduced less. These results suggest that shear modulus, and not viscosity, controls the extent of TM waves at lower pH.

Effect of viscosity on tear drainage and ocular residence time.

PubMed

Zhu, Heng; Chauhan, Anuj

2008-08-01

An increase in residence time of dry eye medications including artificial tears will likely enhance therapeutic benefits. The drainage rates and the residence time of eye drops depend on the viscosity of the instilled fluids. However, a quantitative understanding of the dependence of drainage rates and the residence time on viscosity is lacking. The current study aims to develop a mathematical model for the drainage of Newtonian fluids and also for power-law non-Newtonian fluids of different viscosities. This study is an extension of our previous study on the mathematical model of tear drainage. The tear drainage model is modified to describe the drainage of Newtonian fluids with viscosities higher than the tear viscosity and power-law non-Newtonian fluids with rheological parameters obtained from fitting experimental data in literature. The drainage rate through canaliculi was derived from the modified drainage model and was incorporated into a tear mass balance to calculate the transients of total solute quantity in ocular fluids and the bioavailability of instilled drugs. For Newtonian fluids, increasing the viscosity does not affect the drainage rate unless the viscosity exceeds a critical value of about 4.4 cp. The viscosity has a maximum impact on drainage rate around a value of about 100 cp. The trends are similar for shear thinning power law fluids. The transients of total solute quantity, and the residence time agrees at least qualitatively with experimental studies. A mathematical model has been developed for the drainage of Newtonian fluids and power-law fluids through canaliculi. The model can quantitatively explain different experimental observations on the effect of viscosity on the residence of instilled fluids on the ocular surface. The current study is helpful for understanding the mechanism of fluid drainage from the ocular surface and for improving the design of dry eye treatments.

Implicit Space-Time Conservation Element and Solution Element Schemes

NASA Technical Reports Server (NTRS)

Chang, Sin-Chung; Himansu, Ananda; Wang, Xiao-Yen

1999-01-01

Artificial numerical dissipation is in important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate implicit numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The new schemes presented are two highly accurate implicit solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to convection-dominated equations with very small viscosity. The stability and consistency of the schemes are analysed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme.

An artificial nonlinear diffusivity method for supersonic reacting flows with shocks

NASA Astrophysics Data System (ADS)

Fiorina, B.; Lele, S. K.

2007-03-01

A computational approach for modeling interactions between shocks waves, contact discontinuities and reactions zones with a high-order compact scheme is investigated. To prevent the formation of spurious oscillations around shocks, artificial nonlinear viscosity [A.W. Cook, W.H. Cabot, A high-wavenumber viscosity for high resolution numerical method, J. Comput. Phys. 195 (2004) 594-601] based on high-order derivative of the strain rate tensor is used. To capture temperature and species discontinuities a nonlinear diffusivity based on the entropy gradient is added. It is shown that the damping of 'wiggles' is controlled by the model constants and is largely independent of the mesh size and the shock strength. The same holds for the numerical shock thickness and allows a determination of the L2 error. In the shock tube problem, with fluids of different initial entropy separated by the diaphragm, an artificial diffusivity is required to accurately capture the contact surface. Finally, the method is applied to a shock wave propagating into a medium with non-uniform density/entropy and to a CJ detonation wave. Multi-dimensional formulation of the model is presented and is illustrated by a 2D oblique wave reflection from an inviscid wall, by a 2D supersonic blunt body flow and by a Mach reflection problem.

Applications of Taylor-Galerkin finite element method to compressible internal flow problems

NASA Technical Reports Server (NTRS)

Sohn, Jeong L.; Kim, Yongmo; Chung, T. J.

1989-01-01

A two-step Taylor-Galerkin finite element method with Lapidus' artificial viscosity scheme is applied to several test cases for internal compressible inviscid flow problems. Investigations for the effect of supersonic/subsonic inlet and outlet boundary conditions on computational results are particularly emphasized.

Localized Artificial Viscosity Stabilization of Discontinuous Galerkin Methods for Nonhydrostatic Mesoscale Atmospheric Modeling

DTIC Science & Technology

2014-01-01

with a Riemann flux !"#! (!! ,!!!,), where !!! denotes the solution outside the current element !. Various (approximate) Riemann solvers ...can be used to calculate the Riemann flux, and the Rusanov Riemann solver is adopted in this paper. Then Eq. (7) can be rewritten as !

Relaxation and approximate factorization methods for the unsteady full potential equation

NASA Technical Reports Server (NTRS)

Shankar, V.; Ide, H.; Gorski, J.

1984-01-01

The unsteady form of the full potential equation is solved in conservation form, using implicit methods based on approximate factorization and relaxation schemes. A local time linearization for density is introduced to enable solution to the equation in terms of phi, the velocity potential. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity, to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi obtained from requirements of density continuity. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. Results are presented for flows over airfoils, cylinders, and spheres. Comparisons are made with available Euler and full potential results.

Possible benefits of catheters with lateral holes in coronary thrombus aspiration: a computational study for different clot viscosities and vacuum pressures.

PubMed

Soleimani, Sajjad; Dubini, Gabriele; Pennati, Giancarlo

2014-10-01

According to a number of clinical studies, coronary aspiration catheters are useful tools to remove a thrombus (blood clot) blocking a coronary artery. However, these thrombectomy devices may fail to remove the blood clot entirely. Few studies have been devoted to a systematic analysis of factors affecting clot aspiration. The geometric characteristics of the aspiration catheter, the physical properties of the thrombus, and the applied vacuum pressure are crucial parameters. In this study, the aspiration of a blood clot blocking a coronary bifurcation is computationally simulated. The clot is modeled as a highly viscous fluid, and a two-phase (blood and clot) problem is solved. The effects of geometric variations in the tip of the coronary catheter, including lateral hole size and location, are investigated considering different aspiration pressures and clot viscosities. A Bird-Carreau model is adopted for blood viscosity, while a power law model is used to describe the clot rheology. Computational results for blood clot aspiration show that the presence of holes in the lateral part of the tip of the catheter can be beneficial depending on clot viscosity, hole features, and applied aspiration pressure. In general, the holes are beneficial when the clot viscosity is low, while aspiration catheters without any extra lateral holes exhibit better performance for higher clot viscosity. However, when higher aspiration pressures are applied, the catheters tend to behave relatively similarly in removing clots with various viscosities, reducing the role of the clot viscosity. Copyright © 2014 International Center for Artificial Organs and Transplantation and Wiley Periodicals, Inc.

StarSmasher: Smoothed Particle Hydrodynamics code for smashing stars and planets

NASA Astrophysics Data System (ADS)

Gaburov, Evghenii; Lombardi, James C., Jr.; Portegies Zwart, Simon; Rasio, F. A.

2018-05-01

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian particle method that approximates a continuous fluid as discrete nodes, each carrying various parameters such as mass, position, velocity, pressure, and temperature. In an SPH simulation the resolution scales with the particle density; StarSmasher is able to handle both equal-mass and equal number-density particle models. StarSmasher solves for hydro forces by calculating the pressure for each particle as a function of the particle's properties - density, internal energy, and internal properties (e.g. temperature and mean molecular weight). The code implements variational equations of motion and libraries to calculate the gravitational forces between particles using direct summation on NVIDIA graphics cards. Using a direct summation instead of a tree-based algorithm for gravity increases the accuracy of the gravity calculations at the cost of speed. The code uses a cubic spline for the smoothing kernel and an artificial viscosity prescription coupled with a Balsara Switch to prevent unphysical interparticle penetration. The code also implements an artificial relaxation force to the equations of motion to add a drag term to the calculated accelerations during relaxation integrations. Initially called StarCrash, StarSmasher was developed originally by Rasio.

Study on Solution Properties of Binary Mixtures of Some Industrially Important Solvents with Cyclohexylamine and Cyclohexanone at 298.15 K

NASA Astrophysics Data System (ADS)

Roy, Mahendra Nath; Das, Rajesh Kumar; Chanda, Riju

2010-03-01

Densities and viscosities were measured for the binary mixtures of cyclohexylamine and cyclohexanone with butyl acetate, butanone, butylamine, tert-butylamine, and 2-butoxyethanol at 298.15 K over the entire composition range. From density data, the values of the excess molar volume ( V E) have been calculated. The experimental viscosity data were correlated by means of the equation of Grunberg-Nissan. The density and viscosity data have been analyzed in terms of some semiempirical viscosity models. The results are discussed in terms of molecular interactions and structural effects. The excess molar volume is found to be either negative or positive depending on the molecular interactions and the nature of the liquid mixtures and is discussed in terms of molecular interactions and structural changes.

Determining a membrane's shear modulus, independent of its area-dilatation modulus, via capsule flow in a converging micro-capillary.

PubMed

Dimitrakopoulos, P; Kuriakose, S

2015-04-14

Determination of the elastic properties of the membrane of artificial capsules is essential for the better design of the various devices that are utilized in their engineering and biomedical applications. However this task is complicated owing to the combined effects of the shear and area-dilatation moduli on the capsule deformation. Based on computational investigation, we propose a new methodology to determine a membrane's shear modulus, independent of its area-dilatation modulus, by flowing strain-hardening capsules in a converging micro-capillary of comparable size under Stokes flow conditions, and comparing the experimental measurements of the capsule elongation overshooting with computational data. The capsule prestress, if any, can also be determined with the same methodology. The elongation overshooting is practically independent of the viscosity ratio for low and moderate viscosity ratios, and thus a wide range of capsule fluids can be employed. Our proposed experimental device can be readily produced via glass fabrication while owing to the continuous flow in the micro-capillary, the characterization of a large number of artificial capsules is possible.

[The management of Dysphagia. Part 1: diagnostics].

PubMed

Schönweiler, R

2004-12-01

Patients with chronic dysphagia are often in need of artificial nutrition; though being well balanced in terms of energy and vitamins, patients are at a high risk for the loss of resistance and body weight. Dysphagia also causes a severe drawback of the overall quality of life. This paper gives an overview of the present management of dysphagia from the point of view of otolaryngologists, head-neck-surgeons, phoniatricians, and medical speech-language-voice-pathologists. The physiology of swallowing and typical symptoms of dysphagia are first explained. Then the current most important diagnostic procedures as orofacial and laryngeal function analysis, video-endoscopy, and quantitative assessments, are discussed (part 1). This also includes considerations on bolus viscosity variation, postures, swallowing maneuvers, and sensory enhancement procedures, while actual options like botulinum toxin, passy-muir speaking valve, electromyographic biofeedback, and electrostimulation are also mentioned (part 2).

Experimental study on the pressure and pulse wave propagation in viscoelastic vessel tubes-effects of liquid viscosity and tube stiffness.

PubMed

Ikenaga, Yuki; Nishi, Shohei; Komagata, Yuka; Saito, Masashi; Lagrée, Pierre-Yves; Asada, Takaaki; Matsukawa, Mami

2013-11-01

A pulse wave is the displacement wave which arises because of ejection of blood from the heart and reflection at vascular bed and distal point. The investigation of pressure waves leads to understanding the propagation characteristics of a pulse wave. To investigate the pulse wave behavior, an experimental study was performed using an artificial polymer tube and viscous liquid. A polyurethane tube and glycerin solution were used to simulate a blood vessel and blood, respectively. In the case of the 40 wt% glycerin solution, which corresponds to the viscosity of ordinary blood, the attenuation coefficient of a pressure wave in the tube decreased from 4.3 to 1.6 dB/m because of the tube stiffness (Young's modulus: 60 to 200 kPa). When the viscosity of liquid increased from approximately 4 to 10 mPa·s (the range of human blood viscosity) in the stiff tube, the attenuation coefficient of the pressure wave changed from 1.6 to 3.2 dB/m. The hardening of the blood vessel caused by aging and the increase of blood viscosity caused by illness possibly have opposite effects on the intravascular pressure wave. The effect of the viscosity of a liquid on the amplitude of a pressure wave was then considered using a phantom simulating human blood vessels. As a result, in the typical range of blood viscosity, the amplitude ratio of the waves obtained by the experiments with water and glycerin solution became 1:0.83. In comparison with clinical data, this value is much smaller than that seen from blood vessel hardening. Thus, it can be concluded that the blood viscosity seldom affects the attenuation of a pulse wave.

A comparative study of high-viscosity cement percutaneous vertebroplasty vs. low-viscosity cement percutaneous kyphoplasty for treatment of osteoporotic vertebral compression fractures.

PubMed

Sun, Kai; Liu, Yang; Peng, Hao; Tan, Jun-Feng; Zhang, Mi; Zheng, Xian-Nian; Chen, Fang-Zhou; Li, Ming-Hui

2016-06-01

The clinical effects of two different methods-high-viscosity cement percutaneous vertebroplasty (PVP) and low-viscosity cement percutaneous kyphoplasty (PKP) in the treatment of osteoporotic vertebral compression fractures (OVCFs) were investigated. From June 2010 to August 2013, 98 cases of OVCFs were included in our study. Forty-six patients underwent high-viscosity PVP and 52 patients underwent low-viscosity PKP. The occurrence of cement leakage was observed. Pain relief and functional activity were evaluated using the Visual Analog Scale (VAS) and Oswestry Disability Index (ODI), respectively. Restoration of the vertebral body height and angle of kyphosis were assessed by comparing preoperative and postoperative measurements of the anterior heights, middle heights and the kyphotic angle of the fractured vertebra. Nine out of the 54 vertebra bodies and 11 out of the 60 vertebra bodies were observed to have cement leakage in the high-viscosity PVP and low-viscosity PKP groups, respectively. The rate of cement leakage, correction of anterior vertebral height and kyphotic angles showed no significant differences between the two groups (P>0.05). Low-viscosity PKP had significant advantage in terms of the restoration of middle vertebral height as compared with the high-viscosity PVP (P<0.05). Both groups showed significant improvements in pain relief and functional capacity status after surgery (P<0.05). It was concluded that high-viscosity PVP and low-viscosity PKP have similar clinical effects in terms of the rate of cement leakage, restoration of the anterior vertebral body height, changes of kyphotic angles, functional activity, and pain relief. Low-viscosity PKP is better than high-viscosity PVP in restoring the height of the middle vertebra.

A new flux-limited, two-dimensional, nonsymmetric tensor shock viscosity for DYNA2D: Progress report

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

Benson, D.J.

1988-08-15

A new artificial viscosity that is being developed for DYNA2D is described in this progress report. It incorporates the one- dimensional, flux-limited viscosity developed by Randy Christensen into a two-dimensional setting. In addition, it tries to overcome some of the difficulties associated with traditional bulk viscosity formulations by using an nonsymmetric viscosity tensor. A different nonsymmetric approach (and a symmetrized version of it) is being used in SHALE by Len Margolin with great success. In a later section of this paper, the results of DYNA2D and SHALE for the ''Saltzman problem'' are compared. The general organization of this report ismore » as follows: First it presents the one-dimensional, flux-limited viscosity developed by Christensen. An nonsymmetric tensor viscosity is then developed and the addition of flux-limiting to it is discussed. The issue of boundary conditions turns out to be crucial, and some issues associated with them are unresolved. The current boundary conditions and the issues associated with them are discussed. Example calculations are also given. The current work is summarized and future work is outlined in the first chapter. This brief progress report is the basis of a future journal paper. The work presented in here is unfinished, and some of the results and algorithms will undoubtedly be changed before the final paper. No attempt has been made to survey the results of others in this report, but a limited survey will be incorporated in the journal paper. 7 refs., 39 figs.« less

Time evolution of the eddy viscosity in two-dimensional navier-stokes flow

PubMed

Chaves; Gama

2000-02-01

The time evolution of the eddy viscosity associated with an unforced two-dimensional incompressible Navier-Stokes flow is analyzed by direct numerical simulation. The initial condition is such that the eddy viscosity is isotropic and negative. It is shown by concrete examples that the Navier-Stokes dynamics stabilizes negative eddy viscosity effects. In other words, this dynamics moves monotonically the initial negative eddy viscosity to positive values before relaxation due to viscous term occurs.

Nonlinear second order evolution inclusions with noncoercive viscosity term

NASA Astrophysics Data System (ADS)

Papageorgiou, Nikolaos S.; Rădulescu, Vicenţiu D.; Repovš, Dušan D.

2018-04-01

In this paper we deal with a second order nonlinear evolution inclusion, with a nonmonotone, noncoercive viscosity term. Using a parabolic regularization (approximation) of the problem and a priori bounds that permit passing to the limit, we prove that the problem has a solution.

Renormalizing a viscous fluid model for large scale structure formation

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

Führer, Florian; Rigopoulos, Gerasimos, E-mail: fuhrer@thphys.uni-heidelberg.de, E-mail: gerasimos.rigopoulos@ncl.ac.uk

2016-02-01

Using the Stochastic Adhesion Model (SAM) as a simple toy model for cosmic structure formation, we study renormalization and the removal of the cutoff dependence from loop integrals in perturbative calculations. SAM shares the same symmetry with the full system of continuity+Euler equations and includes a viscosity term and a stochastic noise term, similar to the effective theories recently put forward to model CDM clustering. We show in this context that if the viscosity and noise terms are treated as perturbative corrections to the standard eulerian perturbation theory, they are necessarily non-local in time. To ensure Galilean Invariance higher ordermore » vertices related to the viscosity and the noise must then be added and we explicitly show at one-loop that these terms act as counter terms for vertex diagrams. The Ward Identities ensure that the non-local-in-time theory can be renormalized consistently. Another possibility is to include the viscosity in the linear propagator, resulting in exponential damping at high wavenumber. The resulting local-in-time theory is then renormalizable to one loop, requiring less free parameters for its renormalization.« less

Color change of composite resins subjected to accelerated artificial aging.

PubMed

Tornavoi, Denise Cremonezzi; Agnelli, José Augusto Marcondes; Panzeri, Heitor; Dos Reis, Andréa Cândido

2013-01-01

The aim of this study was to evaluate the influence of accelerated artificial aging (AAA) on the color change of composite resins used in dentistry. Three composite resins were evaluated: Two microhybrids and one hybrid of higher viscosity, with different amounts and sizes of filler particles, shades C2 and B2. A total of 54 specimens were obtained (18 for each composite resin), made of a Teflon matrix (15 mm in diameter and 2 mm in height). The color measurements were obtained with a Spectrophotometer, (PCB 6807 BYK Gardner) before and after AAA. Data were submitted to the Kolmogorov-Smirnov test (α >0.05), ANOVA and Tukey test (α <0.05). After statistical analysis, the color difference among composite resins with the same shades was analyzed. All composite resins showed unacceptable color changes after AAA (ΔE > 3). Considering the variable ∆E, it was observed that the color tone C2 was already statistically different for the microhybrid composite resin prior to AAA (P < 0.05) and in shade B2 for hybrid of higher viscosity and microhybrid with barium glass fluoride aluminum and silica dioxide (P < 0.01). After this process, a statistically significant difference was observed only for shade B2 between microhybrid composite resins (P < 0.01) and for hybrid of higher viscosity and microhybrid with barium glass fluoride aluminum and silica dioxide (P < 0.05). Regarding the color difference within a same composite resin group, before aging the composite resin hybrid of higher viscosity B2 showed the highest color variation rate and microhybrid with zirconium/silica C2 showed the lowest. All composite resins presented unacceptable color changes after 382 h of aging and different composite resins with same hue, presented different colors before being subjected to the aging process (B2 and C2) and after (B2). It was also observed color difference within a group of the same composite resin and same hue.

Design of supercritical swept wings

NASA Technical Reports Server (NTRS)

Garabedian, P.; Mcfadden, G.

1982-01-01

Computational fluid dynamics are used to discuss problems inherent to transonic three-dimensional flow past supercritical swept wings. The formulation for a boundary value problem for the flow past the wing is provided, including consideration of weak shock waves and the use of parabolic coordinates. A swept wing code is developed which requires a mesh of 152 x 10 x 12 points and 200 time cycles. A formula for wave drag is calculated, based on the idea that the conservation form of the momentum equation becomes an entropy inequality measuring the drag, expressible in terms of a small-disturbance equation for a potential function in two dimensions. The entropy inequality has been incorporated in a two-dimensional code for the analysis of transonic flow over airfoils. A method of artificial viscosity is explored for optimum pressure distributions with design, and involves a free boundary problem considering speed over only a portion of the wing.

Reference Correlation for the Viscosity of Carbon Dioxide

NASA Astrophysics Data System (ADS)

Laesecke, Arno; Muzny, Chris D.

2017-03-01

A comprehensive database of experimental and computed data for the viscosity of carbon dioxide (CO2) was compiled and a new reference correlation was developed. Literature results based on an ab initio potential energy surface were the foundation of the correlation of the viscosity in the limit of zero density in the temperature range from 100 to 2000 K. Guided symbolic regression was employed to obtain a new functional form that extrapolates correctly to 0 and to 10 000 K. Coordinated measurements at low density made it possible to implement the temperature dependence of the Rainwater-Friend theory in the linear-in-density viscosity term. The residual viscosity could be formulated with a scaling term ργ/T, the significance of which was confirmed by symbolic regression. The final viscosity correlation covers temperatures from 100 to 2000 K for gaseous CO2 and from 220 to 700 K with pressures along the melting line up to 8000 MPa for compressed and supercritical liquid states. The data representation is more accurate than with the previous correlations, and the covered pressure and temperature range is significantly extended. The critical enhancement of the viscosity of CO2 is included in the new correlation.

L1-Based Approximations of PDEs and Applications

DTIC Science & Technology

2012-09-05

the analysis of the Navier-Stokes equations. The early versions of artificial vis- cosities being overly dissipative, the interest for these technique ...Guermond, and B. Popov. Stability analysis of explicit en- tropy viscosity methods for non-linear scalar conservation equations. Math. Comp., 2012... methods for solv- ing mathematical models of nonlinear phenomena such as nonlinear conservation laws, surface/image/data reconstruction problems

[The influence of sodium hyaluronate solution and artificial tears on higher-order aberrations].

PubMed

Yamashita, Tsutomu; Ochi, Shintarou; Inoue, Yasushi; Miki, Atsushi; Kiryu, Junichi; Tabuchi, Akio

2013-12-01

To investigate the influence of sodium hyaluronate solution (HA) and artificial tears (AT) on higher-order aberrations (HOAs). Twenty four eyes of 24 normal subjects and 11 eyes of 11 dry eye patients were examined. Cornea and ocular wavefront aberrations (total, spherical-like and coma-like) were measured with a Hartmann-Shack wavefront aberrometer before and after 0.1% or 0.3% HA, AT. The consecutively obtained data of the cornea and ocular HOAs were analyzed in the central 4-mm diameter for coma-like, spherical-like and total HOAs. Average HOAs, as well as fluctuation index (FI) and stability index (SI) of the HOAs over time were compared between the two groups. In normal subjects, the AVE of all aberration parameters and FI showed an increase depending on viscosity of the HA (p < 0.001). After AT and 0.1% HA treatment the cornea aberration of the dry eye patients changed from a sawtooth pattern to a stable pattern. Cornea HOAs decreased, and the optical characteristics showed improvement after AT and 0.1% HA in the dry eye patients. HOAs increased depending on the viscosity of the HA, and optical stability worsened.

Quantum hydrodynamics: capturing a reactive scattering resonance.

PubMed

Derrickson, Sean W; Bittner, Eric R; Kendrick, Brian K

2005-08-01

The hydrodynamic equations of motion associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a meshless method based upon a moving least-squares approach. An arbitrary Lagrangian-Eulerian frame of reference and a regridding algorithm which adds and deletes computational points are used to maintain a uniform and nearly constant interparticle spacing. The methodology also uses averaged fields to maintain unitary time evolution. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. A new and more robust artificial viscosity algorithm is presented which gives accurate scattering results and is capable of capturing quantum resonances. The methodology is applied to a one-dimensional model chemical reaction that is known to exhibit a quantum resonance. The correlation function approach is used to compute the reactive scattering matrix, reaction probability, and time delay as a function of energy. Excellent agreement is obtained between the scattering results based upon the quantum hydrodynamic approach and those based upon standard quantum mechanics. This is the first clear demonstration of the ability of moving grid approaches to accurately and robustly reproduce resonance structures in a scattering system.

Preparation of ion-activated in situ gel systems of scopolamine hydrobromide and evaluation of its antimotion sickness efficacy.

PubMed

Cao, Shi-lei; Zhang, Qi-zhi; Jiang, Xin-guo

2007-04-01

To develop a novel, in situ gel system for nasal delivery of scopolamine hydrobromide (SCOP) and study its efficacy on motion sickness. SCOP in situ gels at 0.2%, 0.5%, and 1.0% gellan gum concentration (w/v) were prepared, respectively, and characterized in terms of viscosity, in vitro release, and nasal ciliotoxicity. Single photon emission computing tomography technique was used to evaluate the nasal residence time of gel containing (99m)Tc tracer. The antimotion sickness efficacy produced by the in situ gel formulation was investigated in rats and compared with those achieved after subcutaneous and oral administration. The viscosity of the gellan gum formulations either in solution or in gel increased with increasing concentrations of gellan gum. Its release in vitro was moderate in artificial nasal fluid. The micrographic results showed that in situ gels were safe, without nasal ciliotoxicity. In comparison with phosphate buffer saline, a prolonged radioactivity of (99m)Tc in the rabbit nasal cavity was observed after administration of the gellan gum formulation. Intranasal SCOP in situ gel at a dose of 100 microg/kg decreased symptoms of motion sickness significantly in comparison with subcutaneous and oral administration (P<0.01). SCOP nasal in situ gel is a safe and promising therapeutic alternative to existing medications for motion sickness.

Evaluation of Future Fuels in a High Pressure Common Rail System. Part 3. John Deere 4.5L Powertech Plus

DTIC Science & Technology

2013-01-01

within the fuel injectors regardless of viscosity or lubricity levels. 15. SUBJECT TERMS Fuel Lubricity, Viscosity, HPCR, Synthetic Jet Fuel, Denso... injectors regardless of viscosity or lubricity levels. UNCLASSIFIED v UNCLASSIFIED FOREWORD/ACKNOWLEDGMENTS The U.S. Army TARDEC...40 3.10 UPPER INJECTOR CONNECTING PIN

Reference Correlation for the Viscosity of Carbon Dioxide1

PubMed Central

Laesecke, Arno; Muzny, Chris D.

2017-01-01

A comprehensive database of experimental and computed data for the viscosity of carbon dioxide (CO2) was compiled and a new reference correlation was developed. Literature results based on an ab initio potential energy surface were the foundation of the correlation of the viscosity in the limit of zero density in the temperature range from 100 K to 2000 K. Guided symbolic regression was employed to obtain a new functional form that extrapolates correctly to T → 0 K and to 10 000 K. Coordinated measurements at low density made it possible to implement the temperature dependence of the Rainwater-Friend theory in the linear-in-density viscosity term. The residual viscosity could be formulated with a scaling term ργ/T the significance of which was confirmed by symbolic regression. The final viscosity correlation covers temperatures from 100 K to 2000 K for gaseous CO2, and from 220 K to 700 K with pressures along the melting line up to 8000 MPa for compressed and supercritical liquid states. The data representation is more accurate than with the previous correlations, and the covered pressure and temperature range is significantly extended. The critical enhancement of the viscosity of CO2 is included in the new correlation. PMID:28736460

Entropy-Based Approach To Nonlinear Stability

NASA Technical Reports Server (NTRS)

Merriam, Marshal L.

1991-01-01

NASA technical memorandum suggests schemes for numerical solution of differential equations of flow made more accurate and robust by invoking second law of thermodynamics. Proposes instead of using artificial viscosity to suppress such unphysical solutions as spurious numerical oscillations and nonlinear instabilities, one should formulate equations so that rate of production of entropy within each cell of computational grid be nonnegative, as required by second law.

Dynamical influences on the moment of inertia tensor from lateral viscosity variations inferred from seismic tomographic models

NASA Technical Reports Server (NTRS)

Zhang, Shuxia; Yuen, David A.

1994-01-01

We have investigated the influences of lateral variations of viscosity on the moment of inertia tensor from viscous flows due to the density anomalies in the mantle inferred from seismic tomographic models. The scaling relations between the density and the seismic anomalies is taken as either a constant or a function increasing with depth in accord with the recent high-pressure experimental studies. The viscosity is taken as an exponential function of the 3D density anomaly. In models with an isoviscous background, the effects on the perturbed moment of inertia tensor from the lateral viscosity variations are smaller than those due to variations in the radial viscosity profiles. In mantle models with a background viscosity increasing with depth, the influences of the lateral viscosity variations are significant. The most striking feature in the latter case is that the two off-diagonal elements delta I(sub xz) and delta I(sub yz) in the inertia tensor exhibit greatest sensitivity to lateral variations of the viscosity. While the other elements of the inertia change by only about a few tens of percent in the range of lateral viscosity contrast considered (less than 300), delta I(sub xz) and delta I(sub yz) can vary up to 40 times even with a change in sign, depending on the radial viscosity stratification and the location of the strongest lateral variations. The increase in the velocity-density scaling relation with depth can reduce the influences of the lateral viscosity variations, but it does not change the overall sensitive nature of delta I(sub xz) and delta I(sub yz). This study demonstrates clearly that the lateral viscosity variations, especially in the upper mantle, must be considered in the determination of long-term polar wander, since the variations in the delta I(sub xz) and delta I(sub yz) terms are directly responsible for exciting rotational movements.

The viscous lee wave problem and its implications for ocean modelling

NASA Astrophysics Data System (ADS)

Shakespeare, Callum J.; Hogg, Andrew McC.

2017-05-01

Ocean circulation models employ 'turbulent' viscosity and diffusivity to represent unresolved sub-gridscale processes such as breaking internal waves. Computational power has now advanced sufficiently to permit regional ocean circulation models to be run at sufficiently high (100 m-1 km) horizontal resolution to resolve a significant part of the internal wave spectrum. Here we develop theory for boundary generated internal waves in such models, and in particular, where the waves dissipate their energy. We focus specifically on the steady lee wave problem where stationary waves are generated by a large-scale flow acting across ocean bottom topography. We generalise the energy flux expressions of [Bell, T., 1975. Topographically generated internal waves in the open ocean. J. Geophys. Res. 80, 320-327] to include the effect of arbitrary viscosity and diffusivity. Applying these results for realistic parameter choices we show that in the present generation of models with O(1) m2s-1 horizontal viscosity/diffusivity boundary-generated waves will inevitably dissipate the majority of their energy within a few hundred metres of the boundary. This dissipation is a direct consequence of the artificially high viscosity/diffusivity, which is not always physically justified in numerical models. Hence, caution is necessary in comparing model results to ocean observations. Our theory further predicts that O(10-2) m2s-1 horizontal and O(10-4) m2s-1 vertical viscosity/diffusivity is required to achieve a qualitatively inviscid representation of internal wave dynamics in ocean models.

A new traffic model with a lane-changing viscosity term

NASA Astrophysics Data System (ADS)

Ko, Hung-Tang; Liu, Xiao-He; Guo, Ming-Min; Wu, Zheng

2015-09-01

In this paper, a new continuum traffic flow model is proposed, with a lane-changing source term in the continuity equation and a lane-changing viscosity term in the acceleration equation. Based on previous literature, the source term addresses the impact of speed difference and density difference between adjacent lanes, which provides better precision for free lane-changing simulation; the viscosity term turns lane-changing behavior to a “force” that may influence speed distribution. Using a flux-splitting scheme for the model discretization, two cases are investigated numerically. The case under a homogeneous initial condition shows that the numerical results by our model agree well with the analytical ones; the case with a small initial disturbance shows that our model can simulate the evolution of perturbation, including propagation, dissipation, cluster effect and stop-and-go phenomenon. Project supported by the National Natural Science Foundation of China (Grant Nos. 11002035 and 11372147) and Hui-Chun Chin and Tsung-Dao Lee Chinese Undergraduate Research Endowment (Grant No. CURE 14024).

High speed inviscid compressible flow by the finite element method

NASA Technical Reports Server (NTRS)

Zienkiewicz, O. C.; Loehner, R.; Morgan, K.

1984-01-01

The finite element method and an explicit time stepping algorithm which is based on Taylor-Galerkin schemes with an appropriate artificial viscosity is combined with an automatic mesh refinement process which is designed to produce accurate steady state solutions to problems of inviscid compressible flow in two dimensions. The results of two test problems are included which demonstrate the excellent performance characteristics of the proposed procedures.

Structural properties of aqueous metoprolol succinate solutions. Density, viscosity, and refractive index at 311 K

NASA Astrophysics Data System (ADS)

Deosarkar, S. D.; Kalyankar, T. M.

2013-06-01

Density, viscosity and refractive index of aqueous solutions of metoprolol succinate of different concentrations (0.005-0.05 mol dm-3) were measured at 38°C. Apparent molar volume of resultant solutions were calculated and fitted to the Masson's equation and apparent molar volume at infinite dilution was determined graphically. Viscosity data of solutions has been fitted to the Jone-Dole equation and viscosity A- and B-coefficients were determined graphically. Physicochemical data obtained were discussed in terms of molecular interactions.

Damping of drop oscillations by surfactants and surface viscosity

NASA Technical Reports Server (NTRS)

Rush, Brian M.; Nadim, Ali

1999-01-01

An energy equation is derived for the general case of a viscous drop suspended in a viscous medium with surfactants contaminating the interface. It contains terms that clearly identify dissipation contributions from the viscous effects in the bulk fluids, surface shear and dilatational viscosity effects at the interface, and surfactant transport. An efficient boundary integral method is developed which incorporates the effects of a constant surface dilatational viscosity in simulations of an oscillating two-dimensional inviscid drop. Surface dilatational viscosity is shown to have a significant damping effect on the otherwise undamped inviscid oscillations.

Inflationary universe in terms of a van der Waals viscous fluid

NASA Astrophysics Data System (ADS)

Brevik, I.; Elizalde, E.; Odintsov, S. D.; Timoshkin, A. V.

The inflationary expansion of our early-time universe is considered in terms of the van der Waals equation, as equation of state for the cosmic fluid, where a bulk viscosity contribution is assumed to be present. The corresponding gravitational equations for the energy density in a homogeneous and isotropic Friedmann-Lemaître-Robertson-Walker universe are solved, and an analytic expression for the scale factor is obtained. Attention is paid, specifically, to the role of the viscosity term in the accelerated expansion; the values of the slow-roll parameters, the spectral index, and the tensor-to-scalar ratio for the van der Waals model are calculated and compared with the most recent astronomical data from the Planck satellite. By imposing reasonable restrictions on the parameters of the van der Waals equation, in the presence of viscosity, it is shown to be possible for this model to comply quite precisely with the observational data. One can therefore conclude that the inclusion of viscosity in the theory of the inflationary epoch may definitely improve the cosmological models.

Surface functionalisation with viscosity-sensitive BODIPY molecular rotor

NASA Astrophysics Data System (ADS)

Vyšniauskas, Aurimas; Lopez-Duarte, Ismael; Thompson, Alex J.; Bull, James A.; Kuimova, Marina K.

2018-07-01

Surface functionalisation with viscosity sensitive dyes termed ‘molecular rotors’ can potentially open up new opportunities in sensing, for example for non-invasive biological viscosity imaging, in studying the effect of shear stress on lipid membranes and in cells, and in imaging contacts between surfaces upon applied pressure. We have functionalised microscope slides with BODIPY-based molecular rotor capable of viscosity sensing via its fluorescence lifetime. We have optimised functionalisation conditions and prepared the slides with the BODIPY rotor attached directly to the surface of glass slides and through polymer linkers of 5 kDa and 40 kDa in mass. The slides were characterised for their sensitivity to viscosity, and used to measure viscosity of supported lipid bilayers during photooxidation, and of giant unilamellar vesicles lying on the surface of the slide. We conclude that our functionalised slides show promise for a variety of viscosity sensing applications.

Time accurate application of the MacCormack 2-4 scheme on massively parallel computers

NASA Technical Reports Server (NTRS)

Hudson, Dale A.; Long, Lyle N.

1995-01-01

Many recent computational efforts in turbulence and acoustics research have used higher order numerical algorithms. One popular method has been the explicit MacCormack 2-4 scheme. The MacCormack 2-4 scheme is second order accurate in time and fourth order accurate in space, and is stable for CFL's below 2/3. Current research has shown that the method can give accurate results but does exhibit significant Gibbs phenomena at sharp discontinuities. The impact of adding Jameson type second, third, and fourth order artificial viscosity was examined here. Category 2 problems, the nonlinear traveling wave and the Riemann problem, were computed using a CFL number of 0.25. This research has found that dispersion errors can be significantly reduced or nearly eliminated by using a combination of second and third order terms in the damping. Use of second and fourth order terms reduced the magnitude of dispersion errors but not as effectively as the second and third order combination. The program was coded using Thinking Machine's CM Fortran, a variant of Fortran 90/High Performance Fortran, and was executed on a 2K CM-200. Simple extrapolation boundary conditions were used for both problems.

Simulation of 2-dimensional viscous flow through cascades using a semi-elliptic analysis and hybrid C-H grids

NASA Technical Reports Server (NTRS)

Ramamurti, R.; Ghia, U.; Ghia, K. N.

1988-01-01

A semi-elliptic formulation, termed the interacting parabolized Navier-Stokes (IPNS) formulation, is developed for the analysis of a class of subsonic viscous flows for which streamwise diffusion is neglible but which are significantly influenced by upstream interactions. The IPNS equations are obtained from the Navier-Stokes equations by dropping the streamwise viscous-diffusion terms but retaining upstream influence via the streamwise pressure-gradient. A two-step alternating-direction-explicit numerical scheme is developed to solve these equations. The quasi-linearization and discretization of the equations are carefully examined so that no artificial viscosity is added externally to the scheme. Also, solutions to compressible as well as nearly compressible flows are obtained without any modification either in the analysis or in the solution process. The procedure is applied to constricted channels and cascade passages formed by airfoils of various shapes. These geometries are represented using numerically generated curilinear boundary-oriented coordinates forming an H-grid. A hybrid C-H grid, more appropriate for cascade of airfoils with rounded leading edges, was also developed. Satisfactory results are obtained for flows through cascades of Joukowski airfoils.

Transport Coefficients in weakly compressible turbulence

NASA Technical Reports Server (NTRS)

Rubinstein, Robert; Erlebacher, Gordon

1996-01-01

A theory of transport coefficients in weakly compressible turbulence is derived by applying Yoshizawa's two-scale direct interaction approximation to the compressible equations of motion linearized about a state of incompressible turbulence. The result is a generalization of the eddy viscosity representation of incompressible turbulence. In addition to the usual incompressible eddy viscosity, the calculation generates eddy diffusivities for entropy and pressure, and an effective bulk viscosity acting on the mean flow. The compressible fluctuations also generate an effective turbulent mean pressure and corrections to the speed of sound. Finally, a prediction unique to Yoshizawa's two-scale approximation is that terms containing gradients of incompressible turbulence quantities also appear in the mean flow equations. The form these terms take is described.

Time-asymptotic solutions of the Navier-Stokes equation for free shear flows using an alternating-direction implicit method

NASA Technical Reports Server (NTRS)

Rudy, D. H.; Morris, D. J.

1976-01-01

An uncoupled time asymptotic alternating direction implicit method for solving the Navier-Stokes equations was tested on two laminar parallel mixing flows. A constant total temperature was assumed in order to eliminate the need to solve the full energy equation; consequently, static temperature was evaluated by using algebraic relationship. For the mixing of two supersonic streams at a Reynolds number of 1,000, convergent solutions were obtained for a time step 5 times the maximum allowable size for an explicit method. The solution diverged for a time step 10 times the explicit limit. Improved convergence was obtained when upwind differencing was used for convective terms. Larger time steps were not possible with either upwind differencing or the diagonally dominant scheme. Artificial viscosity was added to the continuity equation in order to eliminate divergence for the mixing of a subsonic stream with a supersonic stream at a Reynolds number of 1,000.

Viscous compressible flow direct and inverse computation and illustrations

NASA Technical Reports Server (NTRS)

Yang, T. T.; Ntone, F.

1986-01-01

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

An artificial viscosity method for the design of supercritical airfoils

NASA Technical Reports Server (NTRS)

Mcfadden, G. B.

1979-01-01

A numerical technique is presented for the design of two-dimensional supercritical wing sections with low wave drag. The method is a design mode of the analysis code H which gives excellent agreement with experimental results and is widely used in the aircraft industry. Topics covered include the partial differential equations of transonic flow, the computational procedure and results; the design procedure; a convergence theorem; and description of the code.

Partial molar volumes and viscosities of aqueous hippuric acid solutions containing LiCl and MnCl2 · 4H2O at 303.15 K

NASA Astrophysics Data System (ADS)

Deosarkar, S. D.; Tawde, P. D.; Zinjade, A. B.; Shaikh, A. I.

2015-09-01

Density (ρ) and viscosity (η) of aqueous hippuric acid (HA) solutions containing LiCl and MnCl2 · 4H2O have been studied at 303.15 K in order to understand volumetric and viscometric behavior of these systems. Apparent molar volume (φv) of salts were calculated from density data and fitted to Massons relation and partial molar volumes (φ{v/0}) at infinite dilution were determined. Relative viscosity data has been used to determine viscosity A and B coefficients using Jones-Dole relation. Partial molar volume and viscosity coefficients have been discussed in terms of ion-solvent interactions and overall structural fittings in solution.

Swallowing Mechanics Associated With Artificial Airways, Bolus Properties, and Penetration-Aspiration Status in Trauma Patients.

PubMed

Dietsch, Angela M; Rowley, Christopher B; Solomon, Nancy Pearl; Pearson, William G

2017-09-18

Artificial airway procedures such as intubation and tracheotomy are common in the treatment of traumatic injuries, and bolus modifications may be implemented to help manage swallowing disorders. This study assessed artificial airway status, bolus properties (volume and viscosity), and the occurrence of laryngeal penetration and/or aspiration in relation to mechanical features of swallowing. Coordinates of anatomical landmarks were extracted at minimum and maximum hyolaryngeal excursion from 228 videofluoroscopic swallowing studies representing 69 traumatically injured U.S. military service members with dysphagia. Morphometric canonical variate and regression analyses examined associations between swallowing mechanics and bolus properties based on artificial airway and penetration-aspiration status. Significant differences in swallowing mechanics were detected between extubated versus tracheotomized (D = 1.32, p < .0001), extubated versus decannulated (D = 1.74, p < .0001), and decannulated versus tracheotomized (D = 1.24, p < .0001) groups per post hoc discriminant function analysis. Tracheotomy-in-situ and decannulated subgroups exhibited increased head/neck extension and posterior relocation of the larynx. Swallowing mechanics associated with (a) penetration-aspiration status and (b) bolus properties were moderately related for extubated and decannulated subgroups, but not the tracheotomized subgroup, per morphometric regression analysis. Specific differences in swallowing mechanics associated with artificial airway status and certain bolus properties may guide therapeutic intervention in trauma-based dysphagia.

Variational approach to the volume viscosity of fluids

NASA Astrophysics Data System (ADS)

Zuckerwar, Allan J.; Ash, Robert L.

2006-04-01

The variational principle of Hamilton is applied to develop an analytical formulation to describe the volume viscosity in fluids. The procedure described here differs from those used in the past in that a dissipative process is represented by the chemical affinity and progress variable (sometimes called "order parameter") of a reacting species. These state variables appear in the variational integral in two places: first, in the expression for the internal energy, and second, in a subsidiary condition accounting for the conservation of the reacting species. As a result of the variational procedure, two dissipative terms appear in the Navier-Stokes equation. The first is the traditional volume viscosity term, proportional to the dilatational component of velocity; the second term is proportional to the material time derivative of the pressure gradient. Values of the respective volume viscosity coefficients are determined by applying the resulting volume-viscous Navier-Stokes equation to the case of acoustical propagation and then comparing expressions for the dispersion and absorption of sound. The formulation includes the special case of equilibration of the translational degrees of freedom. As examples, values are tabulated for dry and humid air, argon, and sea water.

Smoothed particle hydrodynamics simulations of black hole accretion: a step to model black hole feedback in galaxies

NASA Astrophysics Data System (ADS)

Barai, Paramita; Proga, Daniel; Nagamine, Kentaro

2011-11-01

We test how accurately the smoothed particle hydrodynamics (SPH) numerical technique can follow spherically symmetric Bondi accretion. Using the 3D SPH code GADGET-3, we perform simulations of gas accretion on to a central supermassive black hole of mass 108 M⊙ within the radial range of 0.1-200 pc. We carry out simulations without and with radiative heating by a central X-ray corona and radiative cooling. For an adiabatic case, the radial profiles of hydrodynamical properties match the Bondi solution, except near the inner and outer radius of the computational domain. The deviation from the Bondi solution close to the inner radius is caused by the combination of numerical resolution, artificial viscosity and our inner boundary condition. Near the outer radius (≤200 pc), we observe either an outflow or development of a non-spherical inflow unless the outer boundary conditions are very stringently implemented. Despite these issues related to the boundary conditions, we find that adiabatic Bondi accretion can be reproduced for durations of a few dynamical times at the Bondi radius, and for longer times if the outer radius is increased. In particular, the mass inflow rate at the inner boundary, which we measure, is within 3-4 per cent of the Bondi accretion rate. With radiative heating and cooling included, the spherically accreting gas takes a longer time to reach a steady state than the adiabatic Bondi accretion runs, and in some cases does not reach a steady state even within several hundred dynamical times. We find that artificial viscosity causes excessive heating near the inner radius, making the thermal properties of the gas inconsistent with a physical solution. This overheating occurs typically only in the supersonic part of the flow, so that it does not affect the mass accretion rate. We see that increasing the X-ray luminosity produces a lower central mass inflow rate, implying that feedback due to radiative heating is operational in our simulations. With a sufficiently high X-ray luminosity, the inflowing gas is radiatively heated up, and an outflow develops. We conclude that the SPH simulations can capture the gas dynamics needed to study radiative feedback, provided artificial viscosity alters only highly supersonic part of the inflow.

Note: A portable automatic capillary viscometer for transparent and opaque liquids

NASA Astrophysics Data System (ADS)

Soltani Ghalehjooghi, A.; Minaei, S.; Gholipour Zanjani, N.; Beheshti, B.

2017-07-01

A portable automatic capillary viscometer, equipped with an AVR microcontroller, was designed and developed. The viscometer was calibrated with Certified Reference Material (CRM) s200 and utilized for measurement of kinematic viscosity. A quadratic equation was developed for calibration of the instrument at various temperatures. Also, a model was developed for viscosity determination in terms of the viscometer dimensions. Development of the portable viscometer provides for on-site monitoring of engine oil viscosity.

Reversing the direction of galvanotaxis with controlled increases in boundary layer viscosity

NASA Astrophysics Data System (ADS)

Kobylkevich, Brian M.; Sarkar, Anyesha; Carlberg, Brady R.; Huang, Ling; Ranjit, Suman; Graham, David M.; Messerli, Mark A.

2018-05-01

Weak external electric fields (EFs) polarize cellular structure and direct most migrating cells (galvanotaxis) toward the cathode, making it a useful tool during tissue engineering and for healing epidermal wounds. However, the biophysical mechanisms for sensing weak EFs remain elusive. We have reinvestigated the mechanism of cathode-directed water flow (electro-osmosis) in the boundary layer of cells, by reducing it with neutral, viscous polymers. We report that increasing viscosity with low molecular weight polymers decreases cathodal migration and promotes anodal migration in a concentration dependent manner. In contrast, increased viscosity with high molecular weight polymers does not affect directionality. We explain the contradictory results in terms of porosity and hydraulic permeability between the polymers rather than in terms of bulk viscosity. These results provide the first evidence for controlled reversal of galvanotaxis using viscous agents and position the field closer to identifying the putative electric field receptor, a fundamental, outside-in signaling receptor that controls cellular polarity for different cell types.

Reversing the direction of galvanotaxis with controlled increases in boundary layer viscosity.

PubMed

Kobylkevich, Brian M; Sarkar, Anyesha; Carlberg, Brady R; Huang, Ling; Ranjit, Suman; Graham, David M; Messerli, Mark A

2018-03-09

Weak external electric fields (EFs) polarize cellular structure and direct most migrating cells (galvanotaxis) toward the cathode, making it a useful tool during tissue engineering and for healing epidermal wounds. However, the biophysical mechanisms for sensing weak EFs remain elusive. We have reinvestigated the mechanism of cathode-directed water flow (electro-osmosis) in the boundary layer of cells, by reducing it with neutral, viscous polymers. We report that increasing viscosity with low molecular weight polymers decreases cathodal migration and promotes anodal migration in a concentration dependent manner. In contrast, increased viscosity with high molecular weight polymers does not affect directionality. We explain the contradictory results in terms of porosity and hydraulic permeability between the polymers rather than in terms of bulk viscosity. These results provide the first evidence for controlled reversal of galvanotaxis using viscous agents and position the field closer to identifying the putative electric field receptor, a fundamental, outside-in signaling receptor that controls cellular polarity for different cell types.

Freeze-dried, mucoadhesive system for vaginal delivery of the HIV microbicide, dapivirine: optimisation by an artificial neural network.

PubMed

Woolfson, A David; Umrethia, Manish L; Kett, Victoria L; Malcolm, R Karl

2010-03-30

Dapivirine mucoadhesive gels and freeze-dried tablets were prepared using a 3x3x2 factorial design. An artificial neural network (ANN) with multi-layer perception was used to investigate the effect of hydroxypropyl-methylcellulose (HPMC): polyvinylpyrrolidone (PVP) ratio (X1), mucoadhesive concentration (X2) and delivery system (gel or freeze-dried mucoadhesive tablet, X3) on response variables; cumulative release of dapivirine at 24h (Q(24)), mucoadhesive force (F(max)) and zero-rate viscosity. Optimisation was performed by minimising the error between the experimental and predicted values of responses by ANN. The method was validated using check point analysis by preparing six formulations of gels and their corresponding freeze-dried tablets randomly selected from within the design space of contour plots. Experimental and predicted values of response variables were not significantly different (p>0.05, two-sided paired t-test). For gels, Q(24) values were higher than their corresponding freeze-dried tablets. F(max) values for freeze-dried tablets were significantly different (2-4 times greater, p>0.05, two-sided paired t-test) compared to equivalent gels. Freeze-dried tablets having lower values for X1 and higher values for X2 components offered the best compromise between effective dapivirine release, mucoadhesion and viscosity such that increased vaginal residence time was likely to be achieved. Copyright (c) 2009 Elsevier B.V. All rights reserved.

Towards a versatile technique for tracking nanoparticle-mucus interaction: a step on the road

NASA Astrophysics Data System (ADS)

Nafee, N.; Schneider, M.

2014-02-01

Respiratory mucus is one of the main barriers for nanoparticle-based pulmonary delivery systems. This holds true especially for lung diseases like cystic fibrosis, where a very tenacious thick mucus layer hinders particle diffusion to the lung epithelium or the target area. Typically, mean square displacement of particles is used for mobility evaluation. In contrast, our objective is to develop a feasible technique to track directed particle penetration as a prerequisite for efficient pulmonary nanotherapy. Therefore, particle diffusion in artificial mucus was monitored based on confocal laser scanning microscopy (CLSM) and particle-mucus interaction was observed. As pharmaceutical relevant and benign materials, solid lipid nanoparticles (SLNs) were prepared by hot-melt emulsification using glyceryl behenate and different stabilizing agents such as poloxamer-407, tween-80, and polyvinyl alcohol (PVA). The diffusion of labeled SLNs in stained artificial sputum representing CF-patient sputum was verified by 3D time laps imaging. Thus, the effect of coating, particle size and mucus viscosity on nanoparticle diffusion was studied. Using image analysis software "Image J", the total fluorescent signal after 30 min in case of poloxamer-coated SLNs was 5 and 100 folds higher than tween- and PVA-coated SLNs, respectively. Nevertheless, increasing mucus viscosity reduced the diffusion of tweencoated SLNs by a factor of 10. Studying particle-mucus interaction by CLSM can be considered a promising and versatile technique.

Silt and gas accumulation beneath an artificial recharge spreading basin, Southwestern Utah, U.S.A.

USGS Publications Warehouse

Heilweil, V.M.; Solomon, D.K.; Ortiz, G.

2009-01-01

Sand Hollow Reservoir in southwestern Utah, USA, is operated for both surface-water storage and artificial recharge to the underlying Navajo Sandstone. The total volume of estimated artificial recharge between 2002 and 2007 is 85 million cubic meters (69,000 acre-feet). Since 2002, artificial recharge rates have generally been declining and are inversely correlated with the increasing surface area of the reservoir. Permeability testing of core samples retrieved from beneath the reservoir indicates that this decline may not be due to silt accumulation. Artificial recharge rates also show much seasonal variability. Calculations of apparent intrinsic permeability show that these variations can only partly be explained by variation in water viscosity associated with seasonal changes in water temperature. Sporadic seasonal trends in recharge rates and intrinsic permeability during 2002-2004 could be associated with the large fluctuations in reservoir elevation and wetted area. From 2005 through 2007, the reservoir was mostly full and there has been a more consistent seasonal pattern of minimum recharge rates during the summer and maximum rates during the autumn. Total dissolved-gas pressure measurements indicate the presence of biogenic gas bubbles in the shallow sediments beneath the shallower parts of Sand Hollow Reservoir when the water is warmer. Permeability reduction associated with this gas clogging may contribute to the decrease in artificial recharge rates during the spring and summer, with a subsequently increasing recharge rates in the autumn associated with a decline in volume of gas bubbles. Other possible causes for seasonal variation in artificial recharge rates require further investigation.

Third prize: the impact of fluid environment manipulation on shockwave lithotripsy artificial calculi fragmentation rates.

PubMed

Méndez-Probst, Carlos E; Fernadez, Alfonso; Erdeljan, Petar; Vanjecek, Maaike; Cadieux, Peter A; Razvi, Hassan

2011-03-01

Studies have suggested that shockwave lithotripsy (SWL) stone fragmentation rates can be affected by characteristics of the fluid media surrounding the stone, although evidence to implicate the impact of urine specific gravity (SG) is limited and inconclusive. Our aim is to further explore the impact fluid media and SGs have on stone fragmentation using a variable focus lithotripter. Artificial stones were presoaked for 24 hours in urine and then shocked in various fluid media including artificial urine (SG 1.010 control, 1.020, and 1.07), human pooled urine (HPU), degassed HPU, Pentastarch, 100% and 30% contrast, degassed 30% contrast, 100% ethanol, deionized water (dH(2)O), degassed dH(2)O, 5% glucose, Ringer lactate, 0.9% saline, glycerol, whole blood, and lubricating gel. After soaking, SWL using the Modulith SLX-F2 electromagnetic lithotripter was performed. Fragments were dried and sieved using a 4-mm diameter opening grid. Fragments >4 mm were weighed and fragmentation coefficients (FCs) calculated (pre-SWL weight - post-SWL weight)/(pre-SWL weight) × 100. Fifteen stones were shocked for each fluid group. Fluid type, viscosity, and degassing all significantly impacted stone fragmentation. While the solutions' SG, per se, did not appear to affect stone fragmentation, the use of degassed 30% contrast significantly improved stone destruction over the SG 1.010 artificial urine control (95.3% vs 71.4, P < 0.01). Furthermore, degassing improved comminution rates by increasing the number of completely fragmented stones (FC = 100%). Using degassed 30% contrast, 12/15 stones were completely fragmented, compared with only 2/15 in the control group (P = 0.007). Among the whole blood, glycerol, and lubricating gel groups, only 1/15, 0/15, and 1/15 stones reached 100% FC respectively in the narrow focus, possibly because of the detrimental impact of increased viscosity. Different fluid media can significantly affect FC in vitro. Among the various fluids tested, degassed 30% contrast significantly increased the FC and total number of completely fragmented stones.

Scheduling the blended solution as industrial CO2 absorber in separation process by back-propagation artificial neural networks.

PubMed

Abdollahi, Yadollah; Sairi, Nor Asrina; Said, Suhana Binti Mohd; Abouzari-lotf, Ebrahim; Zakaria, Azmi; Sabri, Mohd Faizul Bin Mohd; Islam, Aminul; Alias, Yatimah

2015-11-05

It is believe that 80% industrial of carbon dioxide can be controlled by separation and storage technologies which use the blended ionic liquids absorber. Among the blended absorbers, the mixture of water, N-methyldiethanolamine (MDEA) and guanidinium trifluoromethane sulfonate (gua) has presented the superior stripping qualities. However, the blended solution has illustrated high viscosity that affects the cost of separation process. In this work, the blended fabrication was scheduled with is the process arranging, controlling and optimizing. Therefore, the blend's components and operating temperature were modeled and optimized as input effective variables to minimize its viscosity as the final output by using back-propagation artificial neural network (ANN). The modeling was carried out by four mathematical algorithms with individual experimental design to obtain the optimum topology using root mean squared error (RMSE), R-squared (R(2)) and absolute average deviation (AAD). As a result, the final model (QP-4-8-1) with minimum RMSE and AAD as well as the highest R(2) was selected to navigate the fabrication of the blended solution. Therefore, the model was applied to obtain the optimum initial level of the input variables which were included temperature 303-323 K, x[gua], 0-0.033, x[MDAE], 0.3-0.4, and x[H2O], 0.7-1.0. Moreover, the model has obtained the relative importance ordered of the variables which included x[gua]>temperature>x[MDEA]>x[H2O]. Therefore, none of the variables was negligible in the fabrication. Furthermore, the model predicted the optimum points of the variables to minimize the viscosity which was validated by further experiments. The validated results confirmed the model schedulability. Accordingly, ANN succeeds to model the initial components of the blended solutions as absorber of CO2 capture in separation technologies that is able to industries scale up. Copyright © 2015 Elsevier B.V. All rights reserved.

Electroosmotic flow in capillary channels filled with nonconstant viscosity electrolytes: exact solution of the Navier-Stokes equation.

PubMed

Otevrel, Marek; Klepárník, Karel

2002-10-01

The partial differential equation describing unsteady velocity profile of electroosmotic flow (EOF) in a cylindrical capillary filled with a nonconstant viscosity electrolyte was derived. Analytical solution, based on the general Navier-Stokes equation, was found for constant viscosity electrolytes using the separation of variables (Fourier method). For the case of a nonconstant viscosity electrolyte, the steady-state velocity profile was calculated assuming that the viscosity decreases exponentially in the direction from the wall to the capillary center. Since the respective equations with nonconstant viscosity term are not solvable in general, the method of continuous binding conditions was used to solve this problem. In this method, an arbitrary viscosity profile can be modeled. The theoretical conclusions show that the relaxation times at which an EOF approaches the steady state are too short to have an impact on a separation process in any real systems. A viscous layer at the wall affects EOF significantly, if it is thicker than the Debye length of the electric double layer. The presented description of the EOF dynamics is applicable to any microfluidic systems.

Large Eddy Simulation ... Where Do We Stand? International Workshop Held in St. Petersburg Beach, Florida on 19-21 December 1990.

DTIC Science & Technology

1990-01-01

S. Orszag, Chairman 1. P. Moin Some Issues in Computation of Turbulent Flows. 2. M. Lesieur, P. Comte, X. Normand, 0. Metais and A. Silveira Spectral...Richtmeyer’s computational experience with one-dimensional shock waves (1950) indicated the value of a non-linear artificial viscosity. Charney and... computer architecture and the advantages of semi-Lagrangian advective schemes may lure large-scale atmospheric modelers back to finite-difference

Hemoglobin-vesicle, a cellular artificial oxygen carrier that fulfils the physiological roles of the red blood cell structure.

PubMed

Sakai, Hiromi; Sou, Keitaro; Horinouchi, Hirohisa; Kobayashi, Koichi; Tsuchida, Eishun

2010-01-01

Hb-vesicles (HbV) are artificial O(2) carriers encapsulating concentrated Hb solution (35 g/dL) with a phospholipid bilayer membrane (liposome). The concentration of the HbV suspension is extremely high ([Hb] = 10 g/dL) and it has an O(2) carrying capacity that is comparable to that of blood. HbV is much smaller than RBC (250 vs. 8000 nm), but it recreates the functions of RBCs; (i) the slower rate of O(2) unloading than Hb solution; (ii) colloid osmotic pressure is zero; (iii) the viscosity of a HbV suspension is adjustable to that of blood; (iv) HbV is finally captured by and degraded in RES; (v) co-encapsulation of an allosteric effector to regulate O(2) affinity; (vi) the lipid bilayer membrane prevents direct contact of Hb and vasculature; (vii) NO-binding is retarded to some extent by an intracellular diffusion barrier, and HbV does not induce vasoconstriction. (viii) Both RBC and HbV can be a carrier of not only O(2) but also exogenous CO. However, HbV has limitations such as a shorter functional half-life when compared with RBCs. On the other hand, the advantages of HbV are that it is pathogen-free and blood-type-antigen-free; moreover, it can withstand long-term storage of a few years, none of which can be achieved by the RBC transfusion systems.

Contribution to viscosity from the structural relaxation via the atomic scale Green-Kubo stress correlation function.

PubMed

Levashov, V A

2017-11-14

We studied the connection between the structural relaxation and viscosity for a binary model of repulsive particles in the supercooled liquid regime. The used approach is based on the decomposition of the macroscopic Green-Kubo stress correlation function into the correlation functions between the atomic level stresses. Previously we used the approach to study an iron-like single component system of particles. The role of vibrational motion has been addressed through the demonstration of the relationship between viscosity and the shear waves propagating over large distances. In our previous considerations, however, we did not discuss the role of the structural relaxation. Here we suggest that the contribution to viscosity from the structural relaxation can be taken into account through the consideration of the contribution from the atomic stress auto-correlation term only. This conclusion, however, does not mean that only the auto-correlation term represents the contribution to viscosity from the structural relaxation. Previously the role of the structural relaxation for viscosity has been addressed through the considerations of the transitions between inherent structures and within the mode-coupling theory by other authors. In the present work, we study the structural relaxation through the considerations of the parent liquid and the atomic level stress correlations in it. The comparison with the results obtained on the inherent structures also is made. Our current results suggest, as our previous observations, that in the supercooled liquid regime, the vibrational contribution to viscosity extends over the times that are much larger than the Einstein's vibrational period and much larger than the times that it takes for the shear waves to propagate over the model systems. Besides addressing the atomic level shear stress correlations, we also studied correlations between the atomic level pressure elements.

The effect of ilmenite viscosity on the dynamics and evolution of an overturned lunar cumulate mantle

NASA Astrophysics Data System (ADS)

Zhang, Nan; Dygert, Nick; Liang, Yan; Parmentier, E. M.

2017-07-01

Lunar cumulate mantle overturn and the subsequent upwelling of overturned mantle cumulates provide a potential framework for understanding the first-order thermochemical evolution of the Moon. Upwelling of ilmenite-bearing cumulates (IBCs) after the overturn has a dominant influence on the dynamics and long-term thermal evolution of the lunar mantle. An important parameter determining the stability and convective behavior of the IBC is its viscosity, which was recently constrained through rock deformation experiments. To examine the effect of IBC viscosity on the upwelling of overturned lunar cumulate mantle, here we conduct three-dimensional mantle convection models with an evolving core superposed by an IBC-rich layer, which resulted from mantle overturn after magma ocean solidification. Our modeling shows that a reduction of mantle viscosity by 1 order of magnitude, due to the presence of ilmenite, can dramatically change convective planform and long-term lunar mantle evolution. Our model results suggest a relatively stable partially molten IBC layer that has surrounded the lunar core to the present day.