Improving the accuracy of central difference schemes

NASA Technical Reports Server (NTRS)

Turkel, Eli

1988-01-01

General difference approximations to the fluid dynamic equations require an artificial viscosity in order to converge to a steady state. This artificial viscosity serves two purposes. One is to suppress high frequency noise which is not damped by the central differences. The second purpose is to introduce an entropy-like condition so that shocks can be captured. These viscosities need a coefficient to measure the amount of viscosity to be added. In the standard scheme, a scalar coefficient is used based on the spectral radius of the Jacobian of the convective flux. However, this can add too much viscosity to the slower waves. Hence, it is suggested that a matrix viscosity be used. This gives an appropriate viscosity for each wave component. With this matrix valued coefficient, the central difference scheme becomes closer to upwind biased methods.

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.

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 1: Analysis description

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 Analysis Description, and presents the equations and solution procedure. 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 are described in detail.

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

NASA Technical Reports Server (NTRS)

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

1993-01-01

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal 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 Analysis Description, and presents the equations and solution procedure. It 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.

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.

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

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.

The reality of artificial viscosity

DOE PAGES

Margolin, L. G.

2018-02-24

Artificial viscosity is used in the computer simulation of high Reynolds number flows and is one of the oldest numerical artifices. In this work, I will describe the origin and the interpretation of artificial viscosity as a physical phenomenon. The basis of this interpretation is the finite scale theory, which describes the evolution of integral averages of the fluid solution over finite (length) scales. I will outline the derivation of finite scale Navier–Stokes equations and highlight the particular properties of the equations that depend on the finite scales. Those properties include enslavement, inviscid dissipation, and a law concerning the partitionmore » of total flux of conserved quantities into advective and diffusive components.« less

The reality of artificial viscosity

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

Margolin, L. G.

Artificial viscosity is used in the computer simulation of high Reynolds number flows and is one of the oldest numerical artifices. In this work, I will describe the origin and the interpretation of artificial viscosity as a physical phenomenon. The basis of this interpretation is the finite scale theory, which describes the evolution of integral averages of the fluid solution over finite (length) scales. I will outline the derivation of finite scale Navier–Stokes equations and highlight the particular properties of the equations that depend on the finite scales. Those properties include enslavement, inviscid dissipation, and a law concerning the partitionmore » of total flux of conserved quantities into advective and diffusive components.« less

Accuracy versus convergence rates for a three dimensional multistage Euler code

NASA Technical Reports Server (NTRS)

Turkel, Eli

1988-01-01

Using a central difference scheme, it is necessary to add an artificial viscosity in order to reach a steady state. This viscosity usually consists of a linear fourth difference to eliminate odd-even oscillations and a nonlinear second difference to suppress oscillations in the neighborhood of steep gradients. There are free constants in these differences. As one increases the artificial viscosity, the high modes are dissipated more and the scheme converges more rapidly. However, this higher level of viscosity smooths the shocks and eliminates other features of the flow. Thus, there is a conflict between the requirements of accuracy and efficiency. Examples are presented for a variety of three-dimensional inviscid solutions over isolated wings.

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.

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.

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 !

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.

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.

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.

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

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.

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.

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.

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.

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.

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…

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.

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

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.

Discontinuous Galerkin method with Gaussian artificial viscosity on graphical processing units for nonlinear acoustics

NASA Astrophysics Data System (ADS)

Tripathi, Bharat B.; Marchiano, Régis; Baskar, Sambandam; Coulouvrat, François

2015-10-01

Propagation of acoustical shock waves in complex geometry is a topic of interest in the field of nonlinear acoustics. For instance, simulation of Buzz Saw Noice requires the treatment of shock waves generated by the turbofan through the engines of aeroplanes with complex geometries and wall liners. Nevertheless, from a numerical point of view it remains a challenge. The two main hurdles are to take into account the complex geometry of the domain and to deal with the spurious oscillations (Gibbs phenomenon) near the discontinuities. In this work, first we derive the conservative hyperbolic system of nonlinear acoustics (up to quadratic nonlinear terms) using the fundamental equations of fluid dynamics. Then, we propose to adapt the classical nodal discontinuous Galerkin method to develop a high fidelity solver for nonlinear acoustics. The discontinuous Galerkin method is a hybrid of finite element and finite volume method and is very versatile to handle complex geometry. In order to obtain better performance, the method is parallelized on Graphical Processing Units. Like other numerical methods, discontinuous Galerkin method suffers with the problem of Gibbs phenomenon near the shock, which is a numerical artifact. Among the various ways to manage these spurious oscillations, we choose the method of parabolic regularization. Although, the introduction of artificial viscosity into the system is a popular way of managing shocks, we propose a new approach of introducing smooth artificial viscosity locally in each element, wherever needed. Firstly, a shock sensor using the linear coefficients of the spectral solution is used to locate the position of the discontinuities. Then, a viscosity coefficient depending on the shock sensor is introduced into the hyperbolic system of equations, only in the elements near the shock. The viscosity is applied as a two-dimensional Gaussian patch with its shape parameters depending on the element dimensions, referred here as Element Centered Smooth Artificial Viscosity. Using this numerical solver, various numerical experiments are presented for one and two-dimensional test cases in homogeneous and quiescent medium. This work is funded by CEFIPRA (Indo-French Centre for the Promotion of Advance Research) and partially aided by EGIDE (Campus France).

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.

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.

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.

Entropy-based artificial viscosity stabilization for non-equilibrium Grey Radiation-Hydrodynamics

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

Delchini, Marc O., E-mail: delchinm@email.tamu.edu; Ragusa, Jean C., E-mail: jean.ragusa@tamu.edu; Morel, Jim, E-mail: jim.morel@tamu.edu

2015-09-01

The entropy viscosity method is extended to the non-equilibrium Grey Radiation-Hydrodynamic equations. The method employs a viscous regularization to stabilize the numerical solution. The artificial viscosity coefficient is modulated by the entropy production and peaks at shock locations. The added dissipative terms are consistent with the entropy minimum principle. A new functional form of the entropy residual, suitable for the Radiation-Hydrodynamic equations, is derived. We demonstrate that the viscous regularization preserves the equilibrium diffusion limit. The equations are discretized with a standard Continuous Galerkin Finite Element Method and a fully implicit temporal integrator within the MOOSE multiphysics framework. The methodmore » of manufactured solutions is employed to demonstrate second-order accuracy in both the equilibrium diffusion and streaming limits. Several typical 1-D radiation-hydrodynamic test cases with shocks (from Mach 1.05 to Mach 50) are presented to establish the ability of the technique to capture and resolve shocks.« less

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.

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.

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

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.

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.

A Dynamic Eddy Viscosity Model for the Shallow Water Equations Solved by Spectral Element and Discontinuous Galerkin Methods

NASA Astrophysics Data System (ADS)

Marras, Simone; Suckale, Jenny; Giraldo, Francis X.; Constantinescu, Emil

2016-04-01

We present the solution of the viscous shallow water equations where viscosity is built as a residual-based subgrid scale model originally designed for large eddy simulation of compressible [1] and stratified flows [2]. The necessity of viscosity for a shallow water model not only finds motivation from mathematical analysis [3], but is supported by physical reasoning as can be seen by an analysis of the energetics of the solution. We simulated the flow of an idealized wave as it hits a set of obstacles. The kinetic energy spectrum of this flow shows that, although the inviscid Galerkin solutions -by spectral elements and discontinuous Galerkin [4]- preserve numerical stability in spite of the spurious oscillations in the proximity of the wave fronts, the slope of the energy cascade deviates from the theoretically expected values. We show that only a sufficiently small amount of dynamically adaptive viscosity removes the unwanted high-frequency modes while preserving the overall sharpness of the solution. In addition, it yields a physically plausible energy decay. This work is motivated by a larger interest in the application of a shallow water model to the solution of tsunami triggered coastal flows. In particular, coastal flows in regions around the world where coastal parks made of mitigation hills of different sizes and configurations are considered as a means to deviate the power of the incoming wave. References [1] M. Nazarov and J. Hoffman (2013) "Residual-based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods" Int. J. Numer. Methods Fluids, 71:339-357 [2] S. Marras, M. Nazarov, F. X. Giraldo (2015) "Stabilized high-order Galerkin methods based on a parameter-free dynamic SGS model for LES" J. Comput. Phys. 301:77-101 [3] J. F. Gerbeau and B. Perthame (2001) "Derivation of the viscous Saint-Venant system for laminar shallow water; numerical validation" Discrete Contin. Dyn. Syst. Ser. B, 1:89?102 [4] F. X. Giraldo and M. Restelli (2010) "High-order semi-implicit time-integrators for a triangular discontinuous Galerkin oceanic shallow water model. Int. J. Numer. Methods Fluids, 63:1077-1102

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

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.

A computational modeling approach for the characterization of mechanical properties of 3D alginate tissue scaffolds.

PubMed

Nair, K; Yan, K C; Sun, W

2008-01-01

Scaffold guided tissue engineering is an innovative approach wherein cells are seeded onto biocompatible and biodegradable materials to form 3-dimensional (3D) constructs that, when implanted in the body facilitate the regeneration of tissue. Tissue scaffolds act as artificial extracellular matrix providing the environment conducive for tissue growth. Characterization of scaffold properties is necessary to understand better the underlying processes involved in controlling cell behavior and formation of functional tissue. We report a computational modeling approach to characterize mechanical properties of 3D gellike biomaterial, specifically, 3D alginate scaffold encapsulated with cells. Alginate inherent nonlinearity and variations arising from minute changes in its concentration and viscosity make experimental evaluation of its mechanical properties a challenging and time consuming task. We developed an in silico model to determine the stress-strain relationship of alginate based scaffolds from experimental data. In particular, we compared the Ogden hyperelastic model to other hyperelastic material models and determined that this model was the most suitable to characterize the nonlinear behavior of alginate. We further propose a mathematical model that represents the alginate material constants in Ogden model as a function of concentrations and viscosity. This study demonstrates the model capability to predict mechanical properties of 3D alginate scaffolds.

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.

Development of Artificial Neural Network Model for Diesel Fuel Properties Prediction using Vibrational Spectroscopy.

PubMed

Bolanča, Tomislav; Marinović, Slavica; Ukić, Sime; Jukić, Ante; Rukavina, Vinko

2012-06-01

This paper describes development of artificial neural network models which can be used to correlate and predict diesel fuel properties from several FTIR-ATR absorbances and Raman intensities as input variables. Multilayer feed forward and radial basis function neural networks have been used to rapid and simultaneous prediction of cetane number, cetane index, density, viscosity, distillation temperatures at 10% (T10), 50% (T50) and 90% (T90) recovery, contents of total aromatics and polycyclic aromatic hydrocarbons of commercial diesel fuels. In this study two-phase training procedures for multilayer feed forward networks were applied. While first phase training algorithm was constantly the back propagation one, two second phase training algorithms were varied and compared, namely: conjugate gradient and quasi Newton. In case of radial basis function network, radial layer was trained using K-means radial assignment algorithm and three different radial spread algorithms: explicit, isotropic and K-nearest neighbour. The number of hidden layer neurons and experimental data points used for the training set have been optimized for both neural networks in order to insure good predictive ability by reducing unnecessary experimental work. This work shows that developed artificial neural network models can determine main properties of diesel fuels simultaneously based on a single and fast IR or Raman measurement.

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.

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.

Effect of artificially lengthened vocal tract on vocal fold oscillation's fundamental frequency.

PubMed

Hanamitsu, Masakazu; Kataoka, Hideyuki

2004-06-01

The fundamental frequency of vocal fold oscillation (F(0)) is controlled by laryngeal mechanics and aerodynamic properties. F(0) change per unit change of transglottal pressure (dF/dP) using a shutter valve has been studied and found to have nonlinear, V-shaped relationship with F(0). On the other hand, the vocal tract is also known to affect vocal fold oscillation. This study examined the effect of artificially lengthened vocal tract length on dF/dP. dF/dP was measured in six men using two mouthpieces of different lengths. The dF/dP graph for the longer vocal tract was shifted leftward relative to the shorter one. Using the one-mass model, the nadir of the "V" on the dF/dP graph was strongly influenced by the resonance around the first formant frequency. However, a more precise model is needed to account for the effects of viscosity and turbulence.

Study of heat transfer on physiological driven movement with CNT nanofluids and variable viscosity.

PubMed

Akbar, Noreen Sher; Kazmi, Naeem; Tripathi, Dharmendra; Mir, Nazir Ahmed

2016-11-01

With ongoing interest in CNT nanofluids and materials in biotechnology, energy and environment, microelectronics, composite materials etc., the current investigation is carried out to analyze the effects of variable viscosity and thermal conductivity of CNT nanofluids flow driven by cilia induced movement through a circular cylindrical tube. Metachronal wave is generated by the beating of cilia and mathematically modeled as elliptical wave propagation by Blake (1971). The problem is formulated in the form of nonlinear partial differential equations, which are simplified by using the dimensional analysis to avoid the complicacy of dimensional homogeneity. Lubrication theory is employed to linearize the governing equations and it is also physically appropriate for cilia movement. Analytical solutions for velocity, temperature and pressure gradient and stream function are obtained. The analytical results are numerically simulated by using the Mathematica Software and plotted the graphs for velocity profile, temperature profile, pressure gradient and stream lines for better discussion and visualization. This model is applicable in physiological transport phenomena to explore the nanotechnology in engineering the artificial cilia and ciliated tube/pipe. Copyright © 2016 Elsevier Ireland Ltd. 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.

Experimental investigation of the hydraulic and heat-transfer properties of artificially fractured granite.

PubMed

Luo, Jin; Zhu, Yongqiang; Guo, Qinghai; Tan, Long; Zhuang, Yaqin; Liu, Mingliang; Zhang, Canhai; Xiang, Wei; Rohn, Joachim

2017-01-05

In this paper, the hydraulic and heat-transfer properties of two sets of artificially fractured granite samples are investigated. First, the morphological information is determined using 3D modelling technology. The area ratio is used to describe the roughness of the fracture surface. Second, the hydraulic properties of fractured granite are tested by exposing samples to different confining pressures and temperatures. The results show that the hydraulic properties of the fractures are affected mainly by the area ratio, with a larger area ratio producing a larger fracture aperture and higher hydraulic conductivity. Both the hydraulic apertureand the hydraulic conductivity decrease with an increase in the confining pressure. Furthermore, the fracture aperture decreases with increasing rock temperature, but the hydraulic conductivity increases owing to a reduction of the viscosity of the fluid flowing through. Finally, the heat-transfer efficiency of the samples under coupled hydro-thermal-mechanical conditions is analysed and discussed.

Experimental investigation of the hydraulic and heat-transfer properties of artificially fractured granite

PubMed Central

Luo, Jin; Zhu, Yongqiang; Guo, Qinghai; Tan, Long; Zhuang, Yaqin; Liu, Mingliang; Zhang, Canhai; Xiang, Wei; Rohn, Joachim

2017-01-01

In this paper, the hydraulic and heat-transfer properties of two sets of artificially fractured granite samples are investigated. First, the morphological information is determined using 3D modelling technology. The area ratio is used to describe the roughness of the fracture surface. Second, the hydraulic properties of fractured granite are tested by exposing samples to different confining pressures and temperatures. The results show that the hydraulic properties of the fractures are affected mainly by the area ratio, with a larger area ratio producing a larger fracture aperture and higher hydraulic conductivity. Both the hydraulic apertureand the hydraulic conductivity decrease with an increase in the confining pressure. Furthermore, the fracture aperture decreases with increasing rock temperature, but the hydraulic conductivity increases owing to a reduction of the viscosity of the fluid flowing through. Finally, the heat-transfer efficiency of the samples under coupled hydro-thermal-mechanical conditions is analysed and discussed. PMID:28054594

Incompressible viscous flow computations for the pump components and the artificial heart

NASA Technical Reports Server (NTRS)

Kiris, Cetin

1992-01-01

A finite-difference, three-dimensional incompressible Navier-Stokes formulation to calculate the flow through turbopump components is utilized. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Both steady and unsteady flow calculations can be performed using the current algorithm. In this work, the equations are solved in steadily rotating reference frames by using the steady-state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using an algebraic mixing-length turbulence model. Numerical results are compared with experimental measurements and a good agreement is found between the two. Included in the appendix is a paper on incompressible viscous flow through artificial heart devices with moving boundaries. Time-accurate calculations, such as impeller and diffusor interaction, will be reported in future work.

Extreme learning machine for reduced order modeling of turbulent geophysical flows.

PubMed

San, Omer; Maulik, Romit

2018-04-01

We investigate the application of artificial neural networks to stabilize proper orthogonal decomposition-based reduced order models for quasistationary geophysical turbulent flows. An extreme learning machine concept is introduced for computing an eddy-viscosity closure dynamically to incorporate the effects of the truncated modes. We consider a four-gyre wind-driven ocean circulation problem as our prototype setting to assess the performance of the proposed data-driven approach. Our framework provides a significant reduction in computational time and effectively retains the dynamics of the full-order model during the forward simulation period beyond the training data set. Furthermore, we show that the method is robust for larger choices of time steps and can be used as an efficient and reliable tool for long time integration of general circulation models.

Extreme learning machine for reduced order modeling of turbulent geophysical flows

NASA Astrophysics Data System (ADS)

San, Omer; Maulik, Romit

2018-04-01

We investigate the application of artificial neural networks to stabilize proper orthogonal decomposition-based reduced order models for quasistationary geophysical turbulent flows. An extreme learning machine concept is introduced for computing an eddy-viscosity closure dynamically to incorporate the effects of the truncated modes. We consider a four-gyre wind-driven ocean circulation problem as our prototype setting to assess the performance of the proposed data-driven approach. Our framework provides a significant reduction in computational time and effectively retains the dynamics of the full-order model during the forward simulation period beyond the training data set. Furthermore, we show that the method is robust for larger choices of time steps and can be used as an efficient and reliable tool for long time integration of general circulation models.

Assessment of Beer Quality Based on a Robotic Pourer, Computer Vision, and Machine Learning Algorithms Using Commercial Beers.

PubMed

Gonzalez Viejo, Claudia; Fuentes, Sigfredo; Torrico, Damir D; Howell, Kate; Dunshea, Frank R

2018-05-01

Sensory attributes of beer are directly linked to perceived foam-related parameters and beer color. The aim of this study was to develop an objective predictive model using machine learning modeling to assess the intensity levels of sensory descriptors in beer using the physical measurements of color and foam-related parameters. A robotic pourer (RoboBEER), was used to obtain 15 color and foam-related parameters from 22 different commercial beer samples. A sensory session using quantitative descriptive analysis (QDA ® ) with trained panelists was conducted to assess the intensity of 10 beer descriptors. Results showed that the principal component analysis explained 64% of data variability with correlations found between foam-related descriptors from sensory and RoboBEER such as the positive and significant correlation between carbon dioxide and carbonation mouthfeel (R = 0.62), correlation of viscosity to sensory, and maximum volume of foam and total lifetime of foam (R = 0.75, R = 0.77, respectively). Using the RoboBEER parameters as inputs, an artificial neural network (ANN) regression model showed high correlation (R = 0.91) to predict the intensity levels of 10 related sensory descriptors such as yeast, grains and hops aromas, hops flavor, bitter, sour and sweet tastes, viscosity, carbonation, and astringency. This paper is a novel approach for food science using machine modeling techniques that could contribute significantly to rapid screenings of food and brewage products for the food industry and the implementation of Artificial Intelligence (AI). The use of RoboBEER to assess beer quality showed to be a reliable, objective, accurate, and less time-consuming method to predict sensory descriptors compared to trained sensory panels. Hence, this method could be useful as a rapid screening procedure to evaluate beer quality at the end of the production line for industry applications. © 2018 Institute of Food Technologists®.

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.

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.

The Effects of Thermal Energetics on Three-dimensional Hydrodynamic Instabilities in Massive Protostellar Disks. II. High-Resolution and Adiabatic Evolutions

NASA Astrophysics Data System (ADS)

Pickett, Brian K.; Cassen, Patrick; Durisen, Richard H.; Link, Robert

2000-02-01

In this paper, the effects of thermal energetics on the evolution of gravitationally unstable protostellar disks are investigated by means of three-dimensional hydrodynamic calculations. The initial states for the simulations correspond to stars with equilibrium, self-gravitating disks that are formed early in the collapse of a uniformly rotating, singular isothermal sphere. In a previous paper (Pickett et al.), it was shown that the nonlinear development of locally isentropic disturbances can be radically different than that of locally isothermal disturbances, even though growth in the linear regime may be similar. When multiple low-order modes grew rapidly in the star and inner disk region and saturated at moderate nonlinear levels in the isentropic evolution, the same modes in the isothermal evolution led to shredding of the disk into dense arclets and ejection of material. In this paper, we (1) examine the fate of the shredded disk with calculations at higher spatial resolution than the previous simulations had and (2) follow the evolution of the same initial state using an internal energy equation rather than the assumption of locally isentropic or locally isothermal conditions. Despite the complex structure of the nonlinear features that developed in the violently unstable isothermal disk referred to above, our previous calculation produced no gravitationally independent, long-lived stellar or planetary companions. The higher resolution calculations presented here confirm this result. When the disk of this model is cooled further, prompting even more violent instabilities, the end result is qualitatively the same--a shredded disk. At least for the disks studied here, it is difficult to produce condensations of material that do not shear away into fragmented spirals. It is argued that the ultimate fate of such fragments depends on how readily local internal energy is lost. On the other hand, if a dynamically unstable disk is to survive for very long times without shredding, then some mechanism must mitigate and control any violent phenomena that do occur. The prior simulations demonstrated a marked difference in final outcome, depending upon the efficiency of disk cooling under two different, idealized thermal conditions. We have here incorporated an internal energy equation that allows for arbitrary heating and cooling. Simulations are presented for adiabatic models with and without artificial viscosity. The artificial viscosity accounts for dissipation and heating due to shocks in the code physics. The expected nonaxisymmetric instabilities occur and grow as before in these energy equation evolutions. When artificial viscosity is not present, the model protostar displays behavior between the locally isentropic and locally isothermal cases of the last paper; a strong two-armed spiral grows to nonlinear amplitudes and saturates at a level higher than in the locally isentropic case. Since the amplitude of the spiral disturbance is large, it is expected that continued transport of material and angular momentum will occur well after the end of the calculation at nearly four outer rotation periods. The spiral is not strong enough, however, to disrupt the disk as in the locally isothermal case. When artificial viscosity is present, the same disturbances reach moderate nonlinear amplitude, then heat the gas, which in turn greatly reduces their strength and effects on the disk. Additional heating in the low-density regions of the disk also leads to a gentle flow of material vertically off the computational grid. The energy equation and high-resolution isothermal calculations are used to discuss the importance and relevance of the different thermal regimes so far examined, with particular attention to applications to star and planet formation.

Spermatozoa scattering by a microchannel feature: an elastohydrodynamic model

PubMed Central

Montenegro-Johnson, T. D.; Gadêlha, H.; Smith, D. J.

2015-01-01

Sperm traverse their microenvironment through viscous fluid by propagating flagellar waves; the waveform emerges as a consequence of elastic structure, internal active moments and low Reynolds number fluid dynamics. Engineered microchannels have recently been proposed as a method of sorting and manipulating motile cells; the interaction of cells with these artificial environments therefore warrants investigation. A numerical method is presented for large-amplitude elastohydrodynamic interaction of active swimmers with domain features. This method is employed to examine hydrodynamic scattering by a model microchannel backstep feature. Scattering is shown to depend on backstep height and the relative strength of viscous and elastic forces in the flagellum. In a ‘high viscosity’ parameter regime corresponding to human sperm in cervical mucus analogue, this hydrodynamic contribution to scattering is comparable in magnitude to recent data on contact effects, being of the order of 5°–10°. Scattering can be positive or negative depending on the relative strength of viscous and elastic effects, emphasizing the importance of viscosity on the interaction of sperm with their microenvironment. The modulation of scattering angle by viscosity is associated with variations in flagellar asymmetry induced by the elastohydrodynamic interaction with the boundary feature. PMID:26064617

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.

The pros and cons of code validation

NASA Technical Reports Server (NTRS)

Bobbitt, Percy J.

1988-01-01

Computational and wind tunnel error sources are examined and quantified using specific calculations of experimental data, and a substantial comparison of theoretical and experimental results, or a code validation, is discussed. Wind tunnel error sources considered include wall interference, sting effects, Reynolds number effects, flow quality and transition, and instrumentation such as strain gage balances, electronically scanned pressure systems, hot film gages, hot wire anemometers, and laser velocimeters. Computational error sources include math model equation sets, the solution algorithm, artificial viscosity/dissipation, boundary conditions, the uniqueness of solutions, grid resolution, turbulence modeling, and Reynolds number effects. It is concluded that, although improvements in theory are being made more quickly than in experiments, wind tunnel research has the advantage of the more realistic transition process of a right turbulence model in a free-transition test.

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.

An NPARC Turbulence Module with Wall Functions

NASA Technical Reports Server (NTRS)

Zhu, J.; Shih, T.-H.

1997-01-01

The turbulence module recently developed for the NPARC code has been extended to include wall functions. The Van Driest transformation is used so that the wall functions can be applied to both incompressible and compressible flows. The module is equipped with three two-equation K-epsilon turbulence models: Chien, Shih-Lumley and CMOTR models. Details of the wall functions as well as their numerical implementation are reported. It is shown that the inappropriate artificial viscosity in the near-wall region has a big influence on the solution of the wall function approach. A simple way to eliminate this influence is proposed, which gives satisfactory results during the code validation. The module can be easily linked to the NPARC code for practical applications.

Conservative, special-relativistic smoothed particle hydrodynamics

NASA Astrophysics Data System (ADS)

Rosswog, Stephan

2010-11-01

We present and test a new, special-relativistic formulation of smoothed particle hydrodynamics (SPH). Our approach benefits from several improvements with respect to earlier relativistic SPH formulations. It is self-consistently derived from the Lagrangian of an ideal fluid and accounts for the terms that stem from non-constant smoothing lengths, usually called “grad-h terms”. In our approach, we evolve the canonical momentum and the canonical energy per baryon and thus circumvent some of the problems that have plagued earlier formulations of relativistic SPH. We further use a much improved artificial viscosity prescription which uses the extreme local eigenvalues of the Euler equations and triggers selectively on (a) shocks and (b) velocity noise. The shock trigger accurately monitors the relative density slope and uses it to fine-tune the amount of artificial viscosity that is applied. This procedure substantially sharpens shock fronts while still avoiding post-shock noise. If not triggered, the viscosity parameter of each particle decays to zero. None of these viscosity triggers is specific to special relativity, both could also be applied in Newtonian SPH.The performance of the new scheme is explored in a large variety of benchmark tests where it delivers excellent results. Generally, the grad-h terms deliver minor, though worthwhile, improvements. As expected for a Lagrangian method, it performs close to perfect in supersonic advection tests, but also in strong relativistic shocks, usually considered a particular challenge for SPH, the method yields convincing results. For example, due to its perfect conservation properties, it is able to handle Lorentz factors as large as γ = 50,000 in the so-called wall shock test. Moreover, we find convincing results in a rarely shown, but challenging test that involves so-called relativistic simple waves and also in multi-dimensional shock tube tests.

The upwind control volume scheme for unstructured triangular grids

NASA Technical Reports Server (NTRS)

Giles, Michael; Anderson, W. Kyle; Roberts, Thomas W.

1989-01-01

A new algorithm for the numerical solution of the Euler equations is presented. This algorithm is particularly suited to the use of unstructured triangular meshes, allowing geometric flexibility. Solutions are second-order accurate in the steady state. Implementation of the algorithm requires minimal grid connectivity information, resulting in modest storage requirements, and should enhance the implementation of the scheme on massively parallel computers. A novel form of upwind differencing is developed, and is shown to yield sharp resolution of shocks. Two new artificial viscosity models are introduced that enhance the performance of the new scheme. Numerical results for transonic airfoil flows are presented, which demonstrate the performance of the algorithm.

Effect of Qingnao tablet on blood viscosity of rat model of blood stasis induced by epinephrine

NASA Astrophysics Data System (ADS)

Xie, Guoqi; Hao, Shaojun; Ma, Zhenzhen; Liu, Xiaobin; Li, Jun; Li, Wenjun; Zhang, Zhengchen

2018-04-01

To establish a rat model of blood stasis with adrenaline (Adr) subcutaneous injection and ice bath stimulation. The effects of different doses on the blood viscosity of blood stasis model rats were observed. The rats were randomly divided into 6 groups: blank control group (no model), model group, positive control group, high, middle and low dose group. The whole blood viscosity and plasma viscosity were detected by blood viscosity instrument. Compared with the blank group, model group, high shear, low shear whole blood viscosity and plasma viscosity were significantly increased, TT PT significantly shortened, APTT was significantly prolonged, FIB increased significantly, indicating that the model was successful. Compared with the model group, can significantly reduce the Naoluotong group (cut, low cut). Qingnaopian high dose group (low cut), middle dose group (cut, low shear blood viscosity) (P<0.01), Can significantly reduce Naoluotong qingnaopian group, high dose group (P<0.01), plasma viscosity decreased qingnaopian plasma viscosity in low dose group (P<0.05). Conclusion: qingnaopian could improve the blood rheology of blood stasis mice abnormal index, decrease the blood viscosity, blood stasis has certain hemostatic effect.

CAS2D: FORTRAN program for nonrotating blade-to-blade, steady, potential transonic cascade flows

NASA Technical Reports Server (NTRS)

Dulikravich, D. S.

1980-01-01

An exact, full-potential-equation (FPE) model for the steady, irrotational, homentropic and homoenergetic flow of a compressible, homocompositional, inviscid fluid through two dimensional planar cascades of airfoils was derived, together with its appropriate boundary conditions. A computer program, CAS2D, was developed that numerically solves an artificially time-dependent form of the actual FPE. The governing equation was discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field was discretized by providing a boundary-fitted, nonuniform computational mesh. The mesh was generated by using a sequence of conforming mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the FPE was solved iteratively by using successive line overrelaxation. The possible isentropic shocks were correctly captured by adding explicitly an artificial viscosity in a conservative form. In addition, a three-level consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two dimensional cascade flows.

An Arrhenius-type viscosity function to model sintering using the Skorohod Olevsky viscous sintering model within a finite element code.

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

Ewsuk, Kevin Gregory; Arguello, Jose Guadalupe, Jr.; Reiterer, Markus W.

2006-02-01

The ease and ability to predict sintering shrinkage and densification with the Skorohod-Olevsky viscous sintering (SOVS) model within a finite-element (FE) code have been improved with the use of an Arrhenius-type viscosity function. The need for a better viscosity function was identified by evaluating SOVS model predictions made using a previously published polynomial viscosity function. Predictions made using the original, polynomial viscosity function do not accurately reflect experimentally observed sintering behavior. To more easily and better predict sintering behavior using FE simulations, a thermally activated viscosity function based on creep theory was used with the SOVS model. In comparison withmore » the polynomial viscosity function, SOVS model predictions made using the Arrhenius-type viscosity function are more representative of experimentally observed viscosity and sintering behavior. Additionally, the effects of changes in heating rate on densification can easily be predicted with the Arrhenius-type viscosity function. Another attribute of the Arrhenius-type viscosity function is that it provides the potential to link different sintering models. For example, the apparent activation energy, Q, for densification used in the construction of the master sintering curve for a low-temperature cofire ceramic dielectric has been used as the apparent activation energy for material flow in the Arrhenius-type viscosity function to predict heating rate-dependent sintering behavior using the SOVS model.« less

GIA models with composite rheology and 3D viscosity: effect on GRACE mass balance in Antarctica

NASA Astrophysics Data System (ADS)

van der Wal, Wouter; Whitehouse, Pippa; Schrama, Ernst

2014-05-01

Most Glacial Isostatic Adjustment (GIA) models that have been used to correct GRACE data for the influence of GIA assume a radial stratification of viscosity in the Earth's mantle (1D viscosity). Seismic data in Antarctica indicate that there are large viscosity variations in the horizontal direction (3D viscosity). The purpose of this research is to determine the effect of 3D viscosity on GIA model output, and hence mass balance estimates in Antarctica. We use a GIA model with 3D viscosity and composite rheology in combination with ice loading histories ICE-5G and W12a. From comparisons with uplift and sea-level data in Fennoscandia and North America three preferred viscosity models are selected. For two of the 3D viscosity models the maximum gravity rate due to ICE-5G forcing is located over the Ronne-Filchner ice shelf. This is in contrast with the results obtained using a 1D model, in which the maximum gravity rate due to ICE-5G forcing is always located over the Ross ice shelf. This demonstrates that not all 3D viscosity models can be approximated with a 1D viscosity model. Using CSR release 5 GRACE data from February 2003 to June 2013 mass balance estimates for the three preferred viscosity models are -131 to -171 Gt/year for the ICE-5G model, and -48 to -57 Gt/year for the W12a model. The range due to Earth model uncertainty is larger than the error bar for GRACE (10 Gt/year), but smaller than the range resulting from the difference in ice loading histories.

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.

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.

Computational Fluid Dynamics of Whole-Body Aircraft

NASA Astrophysics Data System (ADS)

Agarwal, Ramesh

1999-01-01

The current state of the art in computational aerodynamics for whole-body aircraft flowfield simulations is described. Recent advances in geometry modeling, surface and volume grid generation, and flow simulation algorithms have led to accurate flowfield predictions for increasingly complex and realistic configurations. As a result, computational aerodynamics has emerged as a crucial enabling technology for the design and development of flight vehicles. Examples illustrating the current capability for the prediction of transport and fighter aircraft flowfields are presented. Unfortunately, accurate modeling of turbulence remains a major difficulty in the analysis of viscosity-dominated flows. In the future, inverse design methods, multidisciplinary design optimization methods, artificial intelligence technology, and massively parallel computer technology will be incorporated into computational aerodynamics, opening up greater opportunities for improved product design at substantially reduced costs.

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.

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.

Estimation of the Viscosities of Liquid Sn-Based Binary Lead-Free Solder Alloys

NASA Astrophysics Data System (ADS)

Wu, Min; Li, Jinquan

2018-01-01

The viscosity of a binary Sn-based lead-free solder alloy was calculated by combining the predicted model with the Miedema model. The viscosity factor was proposed and the relationship between the viscosity and surface tension was analyzed as well. The investigation result shows that the viscosity of Sn-based lead-free solders predicted from the predicted model shows excellent agreement with the reported values. The viscosity factor is determined by three physical parameters: atomic volume, electronic density, and electro-negativity. In addition, the apparent correlation between the surface tension and viscosity of the binary Sn-based Pb-free solder was obtained based on the predicted model.

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.

An entropy correction method for unsteady full potential flows with strong shocks

NASA Technical Reports Server (NTRS)

Whitlow, W., Jr.; Hafez, M. M.; Osher, S. J.

1986-01-01

An entropy correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to account for entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Comparisons between the present method and solutions of the Euler equations and between the present method and experimental data are presented. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation.

Incompressible viscous flow computations for the pump components and the artificial heart

NASA Technical Reports Server (NTRS)

Kiris, Cetin

1992-01-01

A finite difference, three dimensional incompressible Navier-Stokes formulation to calculate the flow through turbopump components is utilized. The solution method is based on the pseudo compressibility approach and uses an implicit upwind differencing scheme together with the Gauss-Seidel line relaxation method. Both steady and unsteady flow calculations can be performed using the current algorithm. Here, equations are solved in steadily rotating reference frames by using the steady state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using an algebraic mixing-length turbulence model. Numerical results are compared with experimental measurements and a good agreement is found between the two.

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.

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

Inference of viscosity jump at 670 km depth and lower mantle viscosity structure from GIA observations

NASA Astrophysics Data System (ADS)

Nakada, Masao; Okuno, Jun'ichi; Irie, Yoshiya

2018-03-01

A viscosity model with an exponential profile described by temperature (T) and pressure (P) distributions and constant activation energy (E_{{{um}}}^{{*}} for the upper mantle and E_{{{lm}}}^* for the lower mantle) and volume (V_{{{um}}}^{{*}} and V_{{{lm}}}^*) is employed in inferring the viscosity structure of the Earth's mantle from observations of glacial isostatic adjustment (GIA). We first construct standard viscosity models with an average upper-mantle viscosity ({\\bar{η }_{{{um}}}}) of 2 × 1020 Pa s, a typical value for the oceanic upper-mantle viscosity, satisfying the observationally derived three GIA-related observables, GIA-induced rate of change of the degree-two zonal harmonic of the geopotential, {\\dot{J}_2}, and differential relative sea level (RSL) changes for the Last Glacial Maximum sea levels at Barbados and Bonaparte Gulf in Australia and for RSL changes at 6 kyr BP for Karumba and Halifax Bay in Australia. Standard viscosity models inferred from three GIA-related observables are characterized by a viscosity of ˜1023 Pa s in the deep mantle for an assumed viscosity at 670 km depth, ηlm(670), of (1 - 50) × 1021 Pa s. Postglacial RSL changes at Southport, Bermuda and Everglades in the intermediate region of the North American ice sheet, largely dependent on its gross melting history, have a crucial potential for inference of a viscosity jump at 670 km depth. The analyses of these RSL changes based on the viscosity models with {\\bar{η }_{{{um}}}} ≥ 2 × 1020 Pa s and lower-mantle viscosity structures for the standard models yield permissible {\\bar{η }_{{{um}}}} and ηlm (670) values, although there is a trade-off between the viscosity and ice history models. Our preferred {\\bar{η }_{{{um}}}} and ηlm (670) values are ˜(7 - 9) × 1020 and ˜1022 Pa s, respectively, and the {\\bar{η }_{{{um}}}} is higher than that for the typical value of oceanic upper mantle, which may reflect a moderate laterally heterogeneous upper-mantle viscosity. The mantle viscosity structure adopted in this study depends on temperature distribution and activation energy and volume, and it is difficult to discuss the impact of each quantity on the inferred lower-mantle viscosity model. We conclude that models of smooth depth variation in the lower-mantle viscosity following η ( z ) ∝ {{ exp}}[ {( {E_{{{lm}}}^* + P( z )V_{{{lm}}}^*} )/{{R}}T( z )} ] with constant E_{{{lm}}}^* and V_{{{lm}}}^* are consistent with the GIA observations.

Inference of mantle viscosity for depth resolutions of GIA observations

NASA Astrophysics Data System (ADS)

Nakada, Masao; Okuno, Jun'ichi

2016-11-01

Inference of the mantle viscosity from observations for glacial isostatic adjustment (GIA) process has usually been conducted through the analyses based on the simple three-layer viscosity model characterized by lithospheric thickness, upper- and lower-mantle viscosities. Here, we examine the viscosity structures for the simple three-layer viscosity model and also for the two-layer lower-mantle viscosity model defined by viscosities of η670,D (670-D km depth) and ηD,2891 (D-2891 km depth) with D-values of 1191, 1691 and 2191 km. The upper-mantle rheological parameters for the two-layer lower-mantle viscosity model are the same as those for the simple three-layer one. For the simple three-layer viscosity model, rate of change of degree-two zonal harmonics of geopotential due to GIA process (GIA-induced J˙2) of -(6.0-6.5) × 10-11 yr-1 provides two permissible viscosity solutions for the lower mantle, (7-20) × 1021 and (5-9) × 1022 Pa s, and the analyses with observational constraints of the J˙2 and Last Glacial Maximum (LGM) sea levels at Barbados and Bonaparte Gulf indicate (5-9) × 1022 Pa s for the lower mantle. However, the analyses for the J˙2 based on the two-layer lower-mantle viscosity model only require a viscosity layer higher than (5-10) × 1021 Pa s for a depth above the core-mantle boundary (CMB), in which the value of (5-10) × 1021 Pa s corresponds to the solution of (7-20) × 1021 Pa s for the simple three-layer one. Moreover, the analyses with the J˙2 and LGM sea level constraints for the two-layer lower-mantle viscosity model indicate two viscosity solutions: η670,1191 > 3 × 1021 and η1191,2891 ˜ (5-10) × 1022 Pa s, and η670,1691 > 1022 and η1691,2891 ˜ (5-10) × 1022 Pa s. The inferred upper-mantle viscosity for such solutions is (1-4) × 1020 Pa s similar to the estimate for the simple three-layer viscosity model. That is, these analyses require a high viscosity layer of (5-10) × 1022 Pa s at least in the deep mantle, and suggest that the GIA-based lower-mantle viscosity structure should be treated carefully in discussing the mantle dynamics related to the viscosity jump at ˜670 km depth. We also preliminarily put additional constraints on these viscosity solutions by examining typical relative sea level (RSL) changes used to infer the lower-mantle viscosity. The viscosity solution inferred from the far-field RSL changes in the Australian region is consistent with those for the J˙2 and LGM sea levels, and the analyses for RSL changes at Southport and Bermuda in the intermediate region for the North American ice sheets suggest the solution of η670,D > 1022, ηD,2891 ˜ (5-10) × 1022 Pa s (D = 1191 or 1691 km) and upper-mantle viscosity higher than 6 × 1020 Pa s.

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.

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.

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, ∞).

Numerical Simulations of Hypersonic Boundary Layer Transition

NASA Astrophysics Data System (ADS)

Bartkowicz, Matthew David

Numerical schemes for supersonic flows tend to use large amounts of artificial viscosity for stability. This tends to damp out the small scale structures in the flow. Recently some low-dissipation methods have been proposed which selectively eliminate the artificial viscosity in regions which do not require it. This work builds upon the low-dissipation method of Subbareddy and Candler which uses the flux vector splitting method of Steger and Warming but identifies the dissipation portion to eliminate it. Computing accurate fluxes typically relies on large grid stencils or coupled linear systems that become computationally expensive to solve. Unstructured grids allow for CFD solutions to be obtained on complex geometries, unfortunately, it then becomes difficult to create a large stencil or the coupled linear system. Accurate solutions require grids that quickly become too large to be feasible. In this thesis a method is proposed to obtain more accurate solutions using relatively local data, making it suitable for unstructured grids composed of hexahedral elements. Fluxes are reconstructed using local gradients to extend the range of data used. The method is then validated on several test problems. Simulations of boundary layer transition are then performed. An elliptic cone at Mach 8 is simulated based on an experiment at the Princeton Gasdynamics Laboratory. A simulated acoustic noise boundary condition is imposed to model the noisy conditions of the wind tunnel and the transitioning boundary layer observed. A computation of an isolated roughness element is done based on an experiment in Purdue's Mach 6 quiet wind tunnel. The mechanism for transition is identified as an instability in the upstream separation region and a comparison is made to experimental data. In the CFD a fully turbulent boundary layer is observed downstream.

Modeling the viscosity of polydisperse suspensions: Improvements in prediction of limiting behavior

NASA Astrophysics Data System (ADS)

Mwasame, Paul M.; Wagner, Norman J.; Beris, Antony N.

2016-06-01

The present study develops a fully consistent extension of the approach pioneered by Farris ["Prediction of the viscosity of multimodal suspensions from unimodal viscosity data," Trans. Soc. Rheol. 12, 281-301 (1968)] to describe the viscosity of polydisperse suspensions significantly improving upon our previous model [P. M. Mwasame, N. J. Wagner, and A. N. Beris, "Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions," J. Rheol. 60, 225-240 (2016)]. The new model captures the Farris limit of large size differences between consecutive particle size classes in a suspension. Moreover, the new model includes a further generalization that enables its application to real, complex suspensions that deviate from ideal non-colloidal suspension behavior. The capability of the new model to predict the viscosity of complex suspensions is illustrated by comparison against experimental data.

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

The viscosity of magmatic silicate liquids: A model for calculation

NASA Technical Reports Server (NTRS)

Bottinga, Y.; Weill, D. F.

1971-01-01

A simple model has been designed to allow reasonably accurate calculations of viscosity as a function of temperature and composition. The problem of predicting viscosities of anhydrous silicate liquids has been investigated since such viscosity numbers are applicable to many extrusive melts and to nearly dry magmatic liquids in general. The fluidizing action of water dissolved in silicate melts is well recognized and it is now possible to predict the effect of water content on viscosity in a semiquantitative way. Water was not incorporated directly into the model. Viscosities of anhydrous compositions were calculated, and, where necessary, the effect of added water and estimated. The model can be easily modified to incorporate the effect of water whenever sufficient additional data are accumulated.

Fractional calculus model of articular cartilage based on experimental stress-relaxation

NASA Astrophysics Data System (ADS)

Smyth, P. A.; Green, I.

2015-05-01

Articular cartilage is a unique substance that protects joints from damage and wear. Many decades of research have led to detailed biphasic and triphasic models for the intricate structure and behavior of cartilage. However, the models contain many assumptions on boundary conditions, permeability, viscosity, model size, loading, etc., that complicate the description of cartilage. For impact studies or biomimetic applications, cartilage can be studied phenomenologically to reduce modeling complexity. This work reports experimental results on the stress-relaxation of equine articular cartilage in unconfined loading. The response is described by a fractional calculus viscoelastic model, which gives storage and loss moduli as functions of frequency, rendering multiple advantages: (1) the fractional calculus model is robust, meaning that fewer constants are needed to accurately capture a wide spectrum of viscoelastic behavior compared to other viscoelastic models (e.g., Prony series), (2) in the special case where the fractional derivative is 1/2, it is shown that there is a straightforward time-domain representation, (3) the eigenvalue problem is simplified in subsequent dynamic studies, and (4) cartilage stress-relaxation can be described with as few as three constants, giving an advantage for large-scale dynamic studies that account for joint motion or impact. Moreover, the resulting storage and loss moduli can quantify healthy, damaged, or cultured cartilage, as well as artificial joints. The proposed characterization is suited for high-level analysis of multiphase materials, where the separate contribution of each phase is not desired. Potential uses of this analysis include biomimetic dampers and bearings, or artificial joints where the effective stiffness and damping are fundamental parameters.

Predicting the viscosity of solids using steady-state creep behavior of the fibrous composites semi-theoretically

NASA Astrophysics Data System (ADS)

Monfared, Vahid

A semi-analytical formulation is presented for obtaining the viscosity of solids (such as metals) using the steady state creep model of the short-fiber composites. For achieving this aim, fluid mechanics theory is used for determining the viscosity. Sometimes, obtaining the viscosity is experimentally difficult and intricate. So, the present model may be beneficial to obtain the viscosity of metals.

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

Matyas, Josef; Cooley, Scott K.; Sundaram, S. K.

Slags of low viscosity readily penetrate the refractory lining in slagging gasifiers, causing rapid and severe corrosion called spalling. In addition, a low-viscosity slag that flows down the gasifier wall forms a relatively thin layer of slag on the refractory surface, allowing the corrosive gases in the gasifier to participate in the chemical reactions between the refractory and the slag. In contrast, a slag viscosity of <25 Pa•s at 1400°C is necessary to minimize the possibility of plugging the slag tap. There is a need to predict and optimize slag viscosity so slagging gasifiers can operate continuously at temperatures rangingmore » from 1300 to 1650°C. The approach adopted in this work was to statistically design and prepare simulated slags, measure the viscosity as a function of temperature, and develop a model to predict slag viscosity based on slag composition and temperature. Statistical design software was used to select compositions from a candidate set of all possible vertices that will optimally represent the composition space for 10 main components. A total of 21 slag compositions were generated, including 5 actual coal slag compositions. The Arrhenius equation was applied to measured viscosity versus temperature data of tested slags, and the Arrhenius coefficients (A and B in ln(vis) = A + B/T) were expressed as linear functions of the slag composition. The viscosity model was validated using 1) data splitting approach, and 2) viscosity/temperature data of selected slag compositions from the literature that were formulated and melted at Pacific Northwest National Laboratory. The capability of the model to predict the viscosity of coal slags was compared with the model developed by Browning et al. because this model can predict the viscosity of slags from coal ash better than the most commonly used empirical models found in the literature.« less

Nonequilibrium viscosity of glass

NASA Astrophysics Data System (ADS)

Mauro, John C.; Allan, Douglas C.; Potuzak, Marcel

2009-09-01

Since glass is a nonequilibrium material, its properties depend on both composition and thermal history. While most prior studies have focused on equilibrium liquid viscosity, an accurate description of nonequilibrium viscosity is essential for understanding the low temperature dynamics of glass. Departure from equilibrium occurs as a glass-forming system is cooled through the glass transition range. The glass transition involves a continuous breakdown of ergodicity as the system gradually becomes trapped in a subset of the available configurational phase space. At very low temperatures a glass is perfectly nonergodic (or “isostructural”), and the viscosity is described well by an Arrhenius form. However, the behavior of viscosity during the glass transition range itself is not yet understood. In this paper, we address the problem of glass viscosity using the enthalpy landscape model of Mauro and Loucks [Phys. Rev. B 76, 174202 (2007)] for selenium, an elemental glass former. To study a wide range of thermal histories, we compute nonequilibrium viscosity with cooling rates from 10-12 to 1012K/s . Based on these detailed landscape calculations, we propose a simplified phenomenological model capturing the essential physics of glass viscosity. The phenomenological model incorporates an ergodicity parameter that accounts for the continuous breakdown of ergodicity at the glass transition. We show a direct relationship between the nonequilibrium viscosity parameters and the fragility of the supercooled liquid. The nonequilibrium viscosity model is validated against experimental measurements of Corning EAGLE XG™ glass. The measurements are performed using a specially designed beam-bending apparatus capable of accurate nonequilibrium viscosity measurements up to 1016Pas . Using a common set of parameters, the phenomenological model provides an accurate description of EAGLE XG™ viscosity over the full range of measured temperatures and fictive temperatures.

Constraints on mantle viscosity from convection models with plate motion history

NASA Astrophysics Data System (ADS)

Mao, W.; Zhong, S.

2017-12-01

The Earth's long-wavelength geoid and dynamic topography are mainly controlled by the mantle buoyancy and viscosity structure. Previous dynamical models for the geoid provide constraints on the 1-D mantle viscosity, using mantle buoyancy derived from seismic topography models. However, it is a challenge in these studies on how to convert seismic velocity to density anomalies and mantle buoyancy. Furthermore, these studies provide constraints only on relative viscosity variations but not on absolute magnitude of viscosity. In this study, we formulate time-dependent 3-D spherical mantle convection models with imposed plate motion history and seek constraints on mantle viscosity structure for both its radial relative variations and its absolute magnitude (i.e., Rayleigh number), using the geoid from the convection models. We found that the geoid at intermediate wavelengths of degrees 4-9 is mainly controlled by the subducted slabs in the upper mantle and the upper part of lower mantle that result from subduction from the last 50 Myr or the Cenozoic. To fit the degrees 4-9 geoid, we need viscosity contrast β defined as the ratio of the lower mantle viscosity and the asthenospheric viscosity to be larger than 2000 and Ra to be 1e8 (defined by the Earth's radius). The best fit model leads to 57% variance reduction and 76% correlation between the model and the observations. However, the long-wavelength geoid at degrees 2-3 is controlled by the lower mantle structure which requires much longer time scale to develop, as seen from our modeling. The preferred viscosity structure and Rayleigh number as constrained by the Cenozoic plate motion and the degrees 4-9 geoid no longer provide adequate fit to the geoid in models with the plate motion history for the last 450 Myr. The degrees 4-9 geoid amplitude is smaller for the models with longer plate motion history and a smaller Ra is required to fit the observation. In order to satisfy the relative amplitude between degrees 2-3 and degrees 4-9 geoid, either a gradually increase of viscosity in the upper part of lower mantle or larger thermal expansivity in the lower mantle is needed. We also consider thermo-chemical models to examine the effects of the African and Pacific thermochemical piles (i.e., LLSVPSs) on the geoid and the inferred mantle viscosity and Ra.

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.

Nonlinear Convective Models of RR Lyrae Stars

NASA Astrophysics Data System (ADS)

Feuchtinger, M.; Dorfi, E. A.

The nonlinear behavior of RR Lyrae pulsations is investigated using a state-of-the-art numerical technique solving the full time-dependent system of radiation hydrodynamics. Grey radiative transfer is included by a variable Eddington-factor method and we use the time-dependent turbulent convection model according to Kuhfuss (1986, A&A 160, 116) in the version of Wuchterl (1995, Comp. Phys. Comm. 89, 19). OPAL opacities extended by the Alexander molecule opacities at temperatures below 6000 K and an equation of state according to Wuchterl (1990, A&A 238, 83) close the system. The resulting nonlinear system is discretized on an adaptive mesh developed by Dorfi & Drury (1987, J. Comp. Phys. 69, 175), which is important to provide the necessary spatial resolution in critical regions like ionization zones and shock waves. Additionally, we employ a second order advection scheme, a time centered temporal discretizaton and an artificial tensor viscosity in order to treat discontinuities. We compute fundamental as well first overtone models of RR Lyrae stars for a grid of stellar parameters both with and without convective energy transport in order to give a detailed picture of the pulsation-convection interaction. In order to investigate the influence of the different features of the convection model calculations with and without overshooting, turbulent pressure and turbulent viscosity are performed and compared with each other. A standard Fourier decomposition is used to confront the resulting light and radial velocity variations with recent observations and we show that the well known RR Lyrae phase discrepancy problem (Simon 1985, ApJ 299, 723) can be resolved with these stellar pulsation computations.

Variability of hemodynamic parameters using the common viscosity assumption in a computational fluid dynamics analysis of intracranial aneurysms.

PubMed

Suzuki, Takashi; Takao, Hiroyuki; Suzuki, Takamasa; Suzuki, Tomoaki; Masuda, Shunsuke; Dahmani, Chihebeddine; Watanabe, Mitsuyoshi; Mamori, Hiroya; Ishibashi, Toshihiro; Yamamoto, Hideki; Yamamoto, Makoto; Murayama, Yuichi

2017-01-01

In most simulations of intracranial aneurysm hemodynamics, blood is assumed to be a Newtonian fluid. However, it is a non-Newtonian fluid, and its viscosity profile differs among individuals. Therefore, the common viscosity assumption may not be valid for all patients. This study aims to test the suitability of the common viscosity assumption. Blood viscosity datasets were obtained from two healthy volunteers. Three simulations were performed for three different-sized aneurysms, two using measured value-based non-Newtonian models and one using a Newtonian model. The parameters proposed to predict an aneurysmal rupture obtained using the non-Newtonian models were compared with those obtained using the Newtonian model. The largest difference (25%) in the normalized wall shear stress (NWSS) was observed in the smallest aneurysm. Comparing the difference ratio to the NWSS with the Newtonian model between the two Non-Newtonian models, the difference of the ratio was 17.3%. Irrespective of the aneurysmal size, computational fluid dynamics simulations with either the common Newtonian or non-Newtonian viscosity assumption could lead to values different from those of the patient-specific viscosity model for hemodynamic parameters such as NWSS.

Fluid dynamic modeling of nano-thermite reactions

NASA Astrophysics Data System (ADS)

Martirosyan, Karen S.; Zyskin, Maxim; Jenkins, Charles M.; Yuki Horie, Yasuyuki

2014-03-01

This paper presents a direct numerical method based on gas dynamic equations to predict pressure evolution during the discharge of nanoenergetic materials. The direct numerical method provides for modeling reflections of the shock waves from the reactor walls that generates pressure-time fluctuations. The results of gas pressure prediction are consistent with the experimental evidence and estimates based on the self-similar solution. Artificial viscosity provides sufficient smoothing of shock wave discontinuity for the numerical procedure. The direct numerical method is more computationally demanding and flexible than self-similar solution, in particular it allows study of a shock wave in its early stage of reaction and allows the investigation of "slower" reactions, which may produce weaker shock waves. Moreover, numerical results indicate that peak pressure is not very sensitive to initial density and reaction time, providing that all the material reacts well before the shock wave arrives at the end of the reactor.

Fluid dynamic modeling of nano-thermite reactions

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

Martirosyan, Karen S., E-mail: karen.martirosyan@utb.edu; Zyskin, Maxim; Jenkins, Charles M.

2014-03-14

This paper presents a direct numerical method based on gas dynamic equations to predict pressure evolution during the discharge of nanoenergetic materials. The direct numerical method provides for modeling reflections of the shock waves from the reactor walls that generates pressure-time fluctuations. The results of gas pressure prediction are consistent with the experimental evidence and estimates based on the self-similar solution. Artificial viscosity provides sufficient smoothing of shock wave discontinuity for the numerical procedure. The direct numerical method is more computationally demanding and flexible than self-similar solution, in particular it allows study of a shock wave in its early stagemore » of reaction and allows the investigation of “slower” reactions, which may produce weaker shock waves. Moreover, numerical results indicate that peak pressure is not very sensitive to initial density and reaction time, providing that all the material reacts well before the shock wave arrives at the end of the reactor.« less

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

USGS Publications Warehouse

Pollitz, Fred F.; Thatcher, Wayne R.

2010-01-01

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

The Thermal Regulation of Gravitational Instabilities in Protoplanetary Disks

NASA Astrophysics Data System (ADS)

Pickett, Brian K.; Mejía, Annie C.; Durisen, Richard H.; Cassen, Patrick M.; Berry, Donald K.; Link, Robert P.

2003-06-01

We present a series of high-resolution, three-dimensional hydrodynamics simulations of a gravitationally unstable solar nebula model. The influences of both azimuthal grid resolution and the treatment of thermal processes on the origin and evolution of gravitational instabilities are investigated. In the first set of simulations, we vary the azimuthal resolution for a locally isothermal simulation, doubling and quadrupling the resolution used in a previous study; the largest number of grid points is (256,256,64) in cylindrical coordinates (r,ϕ,z). At this resolution, the disk breaks apart into a dozen short-lived condensations. Although our previous calculations underresolved the number and growth rate of clumps in the disk, the overall qualitative, but fundamental, conclusion remains: fragmentation under the locally isothermal condition in numerical simulations does not in itself lead to the survival of clumps to become gaseous giant protoplanets. Since local isothermality represents an extreme assumption about thermal processes in the disk, we also present several extended simulations in which heating from an artificial viscosity scheme and cooling from a simple volumetric cooling function are applied to two different models of the solar nebula. The models are differentiated primarily by disk temperature: a high-Q model generated directly by our self-consistent field equilibrium code and a low-Q model generated by cooling the high-Q model in a two-dimensional version of our hydrodynamics code. Here, ``high-Q'' and ``low-Q'' refer to the minimum values of the Toomre stability parameter Q in each disk, Qmin=1.8 and 0.9, respectively. Previous simulations, by ourselves as well as others, have focused on initial states that are already gravitationally unstable, i.e., models similar to the low-Q model. This paper presents for the first time the numerical evolution of an essentially stable initial equilibrium state (the high-Q model) to a severely unstable one by cooling. The additional heating and cooling are applied to each model over the outer half of the disk or the entire disk. The models are subject to the rapid growth of a four-armed spiral instability; the subsequent evolution of the models depends on the thermal behavior of the disk. The cooling function tends to overwhelm the heating included in our artificial viscosity prescription, and as a result the spiral structure strengthens. The spiral disturbances transport mass at prodigious rates during the early nonlinear stages of development and significantly alter the disk's vertical surface. Although dense condensations of material can appear, their character depends on the extent of the volumetric cooling in the disk. In the simulation of the high-Q model with heating and cooling applied throughout the disk, thin, dense rings form at radii ranging from 1 to 3 AU and steadily increase in mass; later companion formation may occur in these rings as cooling drives them toward instability. When heating and cooling are applied only over the outer radial half of the disk, however, a succession of single condensations appears near 5 AU. Each clump has roughly the mass of Saturn, and some survive a complete orbit. Since the clumps form near the artificial boundary in the treatment of the disk gas physics, the production of a clump in this case is a numerical artifact. Nevertheless, radially abrupt transitions in disk gas characteristics, for example, in opacity, might mimic the artificial boundary effects in our simulations and favor the production of stable companions in actual protostellar and protoplanetary disks. The ultimate survival of condensations as eventual stellar or substellar companions to the central star is still largely an open question.

Emergent universe model with dissipative effects

NASA Astrophysics Data System (ADS)

Debnath, P. S.; Paul, B. C.

2017-12-01

Emergent universe model is presented in general theory of relativity with isotropic fluid in addition to viscosity. We obtain cosmological solutions that permit emergent universe scenario in the presence of bulk viscosity that are described by either Eckart theory or Truncated Israel Stewart (TIS) theory. The stability of the solutions are also studied. In this case, the emergent universe (EU) model is analyzed with observational data. In the presence of viscosity, one obtains emergent universe scenario, which however is not permitted in the absence of viscosity. The EU model is compatible with cosmological observations.

Diffusion coefficient and shear viscosity of rigid water models.

PubMed

Tazi, Sami; Boţan, Alexandru; Salanne, Mathieu; Marry, Virginie; Turq, Pierre; Rotenberg, Benjamin

2012-07-18

We report the diffusion coefficient and viscosity of popular rigid water models: two non-polarizable ones (SPC/E with three sites, and TIP4P/2005 with four sites) and a polarizable one (Dang-Chang, four sites). We exploit the dependence of the diffusion coefficient on the system size (Yeh and Hummer 2004 J. Phys. Chem. B 108 15873) to obtain the size-independent value. This also provides an estimate of the viscosity of all water models, which we compare to the Green-Kubo result. In all cases, a good agreement is found. The TIP4P/2005 model is in better agreement with the experimental data for both diffusion and viscosity. The SPC/E and Dang-Chang models overestimate the diffusion coefficient and underestimate the viscosity.

Searching for 3D Viscosity Models of Glacial Isostatic Adjustment in Support of the Global ICE-6G_C Ice History Model

NASA Astrophysics Data System (ADS)

LI, T., II; Wu, P.; Steffen, H.; Wang, H.

2017-12-01

The global ice history model ICE-6G_C was constructed based on the laterally homogeneous earth model VM5a. The combined model of glacial isostatic adjustment (GIA) called ICE-6G_C (VM5a) fits global observations of GIA simultaneously well. However, seismic and geological observations clearly show that the Earth's mantle is laterally heterogeneous. Our aim therefore is to search for the best laterally heterogeneous viscosity models with ICE-6G_C ice history that is able to fit the global relative sea-level (RSL) data, the peak uplift rates (from GNSS) and peak g-dot rates (from the GRACE satellite mission) in Laurentia and Fennoscandia simultaneously. The Coupled Laplace-Finite Element Method is used to compute gravitationally self-consistent sea levels with time dependent coastlines and rotational feedback in addition to changes in deformation, gravity and the state of stress. As a start, the VM5a Earth model is used as the radial background viscosity structure but other radial background viscosity models will also be investigated. Lateral mantle viscosity structure is obtained by the superposition of the radial background viscosity and the lateral viscosity perturbations logarithmically. The latter is inferred from a seismic tomography model using a scaling relationship that takes into account the effects of anharmonicity, anelasticity and non-thermal effects. We will show that several laterally heterogeneous mantle viscosity models can fit the global sea level, GPS and GRACE data better than laterally homogeneous models, provided that the scaling relationship for mantle heterogeneity under northern Europe is allowed to be different from that under Laurentia. In addition, the effects of laterally heterogeneous lithosphere, as inferred from seismic tomography, and the lateral changes in sub-lithospheric properties will also be presented.

Temperature-viscosity models reassessed.

PubMed

Peleg, Micha

2017-05-04

The temperature effect on viscosity of liquid and semi-liquid foods has been traditionally described by the Arrhenius equation, a few other mathematical models, and more recently by the WLF and VTF (or VFT) equations. The essence of the Arrhenius equation is that the viscosity is proportional to the absolute temperature's reciprocal and governed by a single parameter, namely, the energy of activation. However, if the absolute temperature in K in the Arrhenius equation is replaced by T + b where both T and the adjustable b are in °C, the result is a two-parameter model, which has superior fit to experimental viscosity-temperature data. This modified version of the Arrhenius equation is also mathematically equal to the WLF and VTF equations, which are known to be equal to each other. Thus, despite their dissimilar appearances all three equations are essentially the same model, and when used to fit experimental temperature-viscosity data render exactly the same very high regression coefficient. It is shown that three new hybrid two-parameter mathematical models, whose formulation bears little resemblance to any of the conventional models, can also have excellent fit with r 2 ∼ 1. This is demonstrated by comparing the various models' regression coefficients to published viscosity-temperature relationships of 40% sucrose solution, soybean oil, and 70°Bx pear juice concentrate at different temperature ranges. Also compared are reconstructed temperature-viscosity curves using parameters calculated directly from 2 or 3 data points and fitted curves obtained by nonlinear regression using a larger number of experimental viscosity measurements.

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.

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.

Transonic cascade flow calculations using non-periodic C-type grids

NASA Technical Reports Server (NTRS)

Arnone, Andrea; Liou, Meng-Sing; Povinelli, Louis A.

1991-01-01

A new kind of C-type grid is proposed for turbomachinery flow calculations. This grid is nonperiodic on the wake and results in minimum skewness for cascades with high turning and large camber. Euler and Reynolds averaged Navier-Stokes equations are discretized on this type of grid using a finite volume approach. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. Jameson's explicit Runge-Kutta scheme is adopted for the integration in time, and computational efficiency is achieved through accelerating strategies such as multigriding and residual smoothing. A detailed numerical study was performed for a turbine rotor and for a vane. A grid dependence analysis is presented and the effect of artificial dissipation is also investigated. Comparison of calculations with experiments clearly demonstrates the advantage of the proposed grid.

Structure–property reduced order model for viscosity prediction in single-component CO 2 -binding organic liquids

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

Cantu, David C.; Malhotra, Deepika; Koech, Phillip K.

2016-01-01

CO2 capture from power generation with aqueous solvents remains energy intensive due to the high water content of the current technology, or the high viscosity of non-aqueous alternatives. Quantitative reduced models, connecting molecular structure to bulk properties, are key for developing structure-property relationships that enable molecular design. In this work, we describe such a model that quantitatively predicts viscosities of CO2 binding organic liquids (CO2BOLs) based solely on molecular structure and the amount of bound CO2. The functional form of the model correlates the viscosity with the CO2 loading and an electrostatic term describing the charge distribution between the CO2-bearingmore » functional group and the proton-receiving amine. Molecular simulations identify the proton shuttle between these groups within the same molecule to be the critical indicator of low viscosity. The model, developed to allow for quick screening of solvent libraries, paves the way towards the rational design of low viscosity non-aqueous solvent systems for post-combustion CO2 capture. Following these theoretical recommendations, synthetic efforts of promising candidates and viscosity measurement provide experimental validation and verification.« less

Live Cell Imaging of Viscosity in 3D Tumour Cell Models.

PubMed

Shirmanova, Marina V; Shimolina, Lubov' E; Lukina, Maria M; Zagaynova, Elena V; Kuimova, Marina K

2017-01-01

Abnormal levels of viscosity in tissues and cells are known to be associated with disease and malfunction. While methods to measure bulk macroscopic viscosity of bio-tissues are well developed, imaging viscosity at the microscopic scale remains a challenge, especially in vivo. Molecular rotors are small synthetic viscosity-sensitive fluorophores in which fluorescence parameters are strongly correlated to the microviscosity of their immediate environment. Hence, molecular rotors represent a promising instrument for mapping of viscosity in living cells and tissues at the microscopic level. Quantitative measurements of viscosity can be achieved by recording time-resolved fluorescence decays of molecular rotor using fluorescence lifetime imaging microscopy (FLIM), which is also suitable for dynamic viscosity mapping, both in cellulo and in vivo. Among tools of experimental oncology, 3D tumour cultures, or spheroids, are considered a more adequate in vitro model compared to a cellular monolayer, and represent a less labour-intensive and more unified approach compared to animal tumour models. This chapter describes a methodology for microviscosity imaging in tumour spheroids using BODIPY-based molecular rotors and two photon-excited FLIM.

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.

Effect of interactions between multiple interfaces on the rheological characteristics of double emulsions

NASA Astrophysics Data System (ADS)

Choi, Se Bin; Park, Jae Yong; Moon, Ji Young; Lee, Joon Sang

2018-06-01

In this study, we analyzed the rheological characteristics of double emulsions by using a three-dimensional lattice Boltzmann model. Numerical simulations indicate that interactions between multiple interfaces play a vital role in determining the shear stress on interfaces and affect deformations, which influence the relative viscosity of double emulsions. The large shear stress induced by droplets in contact increases the relative viscosity for high volume fractions. The double emulsions also show shear-thinning behavior, which corresponds with the Carreau model. The interfacial interference between the core and the deforming shell cause the relative viscosity to increase with increasing core-droplet radius. Finally, we investigated the dependence of the double-emulsion viscosity on the core-droplet viscosity. At high shear rates, the relative viscosity increases with increasing core-droplet viscosity. However, the trend is opposite at low shear rates, which results from the high inward flow (Marangoni flow) at low core-droplet viscosity.

Evaluating the Sensitivity of Glacial Isostatic Adjustment to a Hydrous Melt at 410 km Depth

NASA Astrophysics Data System (ADS)

Hill, A. M.; Milne, G. A.; Ranalli, G.

2017-12-01

We present a sensitivity analysis aimed at testing whether observables related to GIA can support or refute the existence of a low viscosity partial melt layer located above the mantle transition zone, as required by the so-called "Transition Zone Water Filter" model (Bercovici and Karato 2003). In total, 400 model runs were performed sampling a range of melt layer thicknesses (1, 10 & 20 km) and viscosities (1015 - 1019 Pas) as well as plausible viscosity values in the upper and lower mantle. Comparing model output of postglacial decay times and j2, 18 of the considered viscosity models were found to be compatible with all of the observational constraints. Amongst these, only three `background' upper and lower mantle viscosities are permitted regardless of the properties of the melt layer: an upper mantle value of 3×1020 Pas and lower mantle values of 1022, 3×1022 and 5×1022 Pas. Concerning the properties of the melt layer itself, a thin (1 km) layer may have any of the investigated viscosities (1015 to 1019 Pas). For thicker melt layers, the viscosity must be ≥1018 Pas (20 km) or ≥1017 Pas (10 km). Our results indicate clear parameter trade-offs between the properties of the melt layer and the background viscosity structure. Given that the observations permit several values of lower mantle viscosity, we conclude that tightening constraints on this parameter would be valuable for future investigation of the type presented here. Furthermore, while decay times from both locations considered in this investigation (Ångerman River, Sweden; Richmond Gulf, Canada) offer meaningful constraints on viscosity structure, the value for Richmond Gulf is significantly more uncertain and so increasing its precision would likely result in improved viscosity constraints.

Subduction zone evolution and low viscosity wedges and channels

NASA Astrophysics Data System (ADS)

Manea, Vlad; Gurnis, Michael

2007-12-01

Dehydration of subducting lithosphere likely transports fluid into the mantle wedge where the viscosity is decreased. Such a decrease in viscosity could form a low viscosity wedge (LVW) or a low viscosity channel (LVC) on top of the subducting slab. Using numerical models, we investigate the influence of low viscosity wedges and channels on subduction zone structure. Slab dip changes substantially with the viscosity reduction within the LVWs and LVCs. For models with or without trench rollback, overthickening of slabs is greatly reduced by LVWs or LVCs. Two divergent evolutionary pathways have been found depending on the maximum depth extent of the LVW and wedge viscosity. Assuming a viscosity contrast of 0.1 with background asthenosphere, models with a LVW that extends down to 400 km depth show a steeply dipping slab, while models with an LVW that extends to much shallower depth, such as 200 km, can produce slabs that are flat lying beneath the overriding plate. There is a narrow range of mantle viscosities that produces flat slabs (5 to10 × 10 19 Pa s) and the slab flattening process is enhanced by trench rollback. Slab can be decoupled from the overriding plate with a LVC if the thickness is at least a few 10 s of km, the viscosity reduction is at least a factor of two and the depth extent of the LVC is several hundred km. These models have important implications for the geochemical and spatial evolution of volcanic arcs and the state of stress within the overriding plate. The models explain the poor correlation between traditional geodynamic controls, subducting plate age and convergence rates, on slab dip. We predict that when volcanic arcs change their distance from the trench, they could be preceded by changes in arc chemistry. We predict that there could be a larger volatile input into the wedge when arcs migrate toward the trench and visa-versa. The transition of a subduction zone into the flat-lying regime could be preceded by changes in the volatile budget such that the dehydration front moves to shallower depths. Our flat-slab models shed some light on puzzling flat subduction systems, like in Central Mexico, where there is no deformation within the overriding plate above the flat segment. The lack of in-plane compression in Central Mexico suggests the presence of a low viscosity shear zone above the flat slab.

A Study of Oil Viscosity Mental Model

NASA Astrophysics Data System (ADS)

Albaiti; Liliasari; Sumarna, Omay; Abdulkadir Martoprawiro, Muhamad

2017-02-01

There is no study regarding on how to learn viscosity of the liquid (e.g. oil) by interconnecting macroscopic, sub-microscopic and symbolic levels. Therefore, the purpose of this research was to study the mental model of the oil viscosity. Intermolecular attractive force of oil constituent on the sub-microscopic level is depicted in the form of mental models. In this research, the viscosity data for some types of oil was measured by using Hoppler method. Viscosity of mineral oil SAE 20W-50, mineral oil SAE 15W-40 and synthetic oil SAE 10W-40 were 1.75, 1.31, and 1.03 Pa s, and the densities of these oils were 908.64, 885.04, and 877.02 kg/m3, respectively. The results showed that the greater density of the mineral oil that is assumed to be composed of linear chains of hydrocarbons, the longer the chain of hydrocarbon linear. Consequently, there are stronger the London force and greater the oil viscosity. The density and viscosity of synthetic oil are lower than that of both mineral oils. Synthetic oil structurally forms polymers with large branching. This structure affects a lower synthetic oil viscosity. This study contributes to construct a mental model of pre-service chemistry teachers.

Measuring viscosity with a resonant magnetic perturbation in the MST RFP

NASA Astrophysics Data System (ADS)

Fridström, Richard; Munaretto, Stefano; Frassinetti, Lorenzo; Chapman, Brett; Brunsell, Per; Sarff, John; MST Team

2016-10-01

Application of an m = 1 resonant magnetic perturbation (RMP) causes braking and locking of naturally rotating m = 1 tearing modes (TMs) in the MST RFP. The experimental TM dynamics are replicated by a theoretical model including the interaction between the RMP and multiple TMs [Fridström PoP 23, 062504 (2016)]. The viscosity is the only free parameter in the model, and it is chosen such that model TM velocity evolution matches that of the experiment. The model does not depend on the means by which the natural rotation is generated. The chosen value of the viscosity, about 40 m2/s, is consistent with separate measurements in MST using a biased probe to temporarily spin up the plasma. This viscosity is about 100 times larger than the classical prediction, likely due to magnetic stochasticity in the core of these plasmas. Viscosity is a key parameter in visco-resistive MHD codes like NIMROD. The validation of these codes requires measurement of the viscosity over a broad parameter range, which will now be possible with the RMP technique that, unlike the biased probe, is not limited to low-energy-density plasmas. Estimation with the RMP technique of the viscosity in several MST discharges suggests that the viscosity decreases as the electron beta increases. Work supported by USDOE.

A method for matching the refractive index and kinematic viscosity of a blood analog for flow visualization in hydraulic cardiovascular models.

PubMed

Nguyen, T T; Biadillah, Y; Mongrain, R; Brunette, J; Tardif, J C; Bertrand, O F

2004-08-01

In this work, we propose a simple method to simultaneously match the refractive index and kinematic viscosity of a circulating blood analog in hydraulic models for optical flow measurement techniques (PIV, PMFV, LDA, and LIF). The method is based on the determination of the volumetric proportions and temperature at which two transparent miscible liquids should be mixed to reproduce the targeted fluid characteristics. The temperature dependence models are a linear relation for the refractive index and an Arrhenius relation for the dynamic viscosity of each liquid. Then the dynamic viscosity of the mixture is represented with a Grunberg-Nissan model of type 1. Experimental tests for acrylic and blood viscosity were found to be in very good agreement with the targeted values (measured refractive index of 1.486 and kinematic viscosity of 3.454 milli-m2/s with targeted values of 1.47 and 3.300 milli-m2/s).

Anomalous intrinsic viscosity of octadecylamine-functionalised carbon nanotubes in suspension.

PubMed

Donovan, K J; Scott, K

2013-06-28

Single walled carbon nanotubes, SWCNTs, are used as a model cylinder of nanoscopic dimensions for testing rheological theories of how addition of cylindrical particles affects the viscosity of a suspension of such particles. Using the rate of growth of the accompanying induced linear dichroism following application of an applied electric field, the dynamics of carbon nanotube alignment is studied in suspensions of octadecylamine functionalised single walled carbon nanotubes, ODA-SWCNTs, in 1,2 dichloroethane. From such measurements the viscosity of the suspension is measured as the concentration of the suspension is varied. While working within the dilute limit the viscosity is found to increase linearly with concentration and the intrinsic viscosity of the suspension is found to be 8000. This anomalously high intrinsic viscosity is compared with the predictions of various models for a rigid cylinder and found to be incompatible with any of the current models. Some suggestions are made as to the way this ODA-SWCNT result may be eventually accommodated within other models.

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.

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.

Swimming efficiency in a shear-thinning fluid

NASA Astrophysics Data System (ADS)

Nganguia, Herve; Pietrzyk, Kyle; Pak, On Shun

2017-12-01

Micro-organisms expend energy moving through complex media. While propulsion speed is an important property of locomotion, efficiency is another factor that may determine the swimming gait adopted by a micro-organism in order to locomote in an energetically favorable manner. The efficiency of swimming in a Newtonian fluid is well characterized for different biological and artificial swimmers. However, these swimmers often encounter biological fluids displaying shear-thinning viscosities. Little is known about how this nonlinear rheology influences the efficiency of locomotion. Does the shear-thinning rheology render swimming more efficient or less? How does the swimming efficiency depend on the propulsion mechanism of a swimmer and rheological properties of the surrounding shear-thinning fluid? In this work, we address these fundamental questions on the efficiency of locomotion in a shear-thinning fluid by considering the squirmer model as a general locomotion model to represent different types of swimmers. Our analysis reveals how the choice of surface velocity distribution on a squirmer may reduce or enhance the swimming efficiency. We determine optimal shear rates at which the swimming efficiency can be substantially enhanced compared with the Newtonian case. The nontrivial variations of swimming efficiency prompt questions on how micro-organisms may tune their swimming gaits to exploit the shear-thinning rheology. The findings also provide insights into how artificial swimmers should be designed to move through complex media efficiently.

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.

Tidal dissipation in a viscoelastic planet

NASA Technical Reports Server (NTRS)

Ross, M.; Schubert, G.

1986-01-01

Tidal dissipation is examined using Maxwell standard liner solid (SLS), and Kelvin-Voigt models, and viscosity parameters are derived from the models that yield the amount of dissipation previously calculated for a moon model with QW = 100 in a hypothetical orbit closer to the earth. The relevance of these models is then assessed for simulating planetary tidal responses. Viscosities of 10 exp 14 and 10 ex 18 Pa s for the Kelvin-Voigt and Maxwell rheologies, respectively, are needed to match the dissipation rate calculated using the Q approach with a quality factor = 100. The SLS model requires a short time viscosity of 3 x 10 exp 17 Pa s to match the Q = 100 dissipation rate independent of the model's relaxation strength. Since Q = 100 is considered a representative value for the interiors of terrestrial planets, it is proposed that derived viscosities should characterize planetary materials. However, it is shown that neither the Kelvin-Voigt nor the SLS models simulate the behavior of real planetary materials on long time scales. The Maxwell model, by contrast, behaves realistically on both long and short time scales. The inferred Maxwell viscosity, corresponding to the time scale of days, is several times smaller than the longer time scale (greater than or equal to 10 exp 14 years) viscosity of the earth's mantle.

An exploration of viscosity models in the realm of kinetic theory of liquids originated fluids

NASA Astrophysics Data System (ADS)

Hussain, Azad; Ghafoor, Saadia; Malik, M. Y.; Jamal, Sarmad

The preeminent perspective of this article is to study flow of an Eyring Powell fluid model past a penetrable plate. To find the effects of variable viscosity on fluid model, continuity, momentum and energy equations are elaborated. Here, viscosity is taken as function of temperature. To understand the phenomenon, Reynold and Vogel models of variable viscosity are incorporated. The highly non-linear partial differential equations are transfigured into ordinary differential equations with the help of suitable similarity transformations. The numerical solution of the problem is presented. Graphs are plotted to visualize the behavior of pertinent parameters on the velocity and temperature profiles.

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.

Development of Viscosity Model for Petroleum Industry Applications

NASA Astrophysics Data System (ADS)

Motahhari, Hamed reza

Heavy oil and bitumen are challenging to produce and process due to their very high viscosity, but their viscosity can be reduced either by heating or dilution with a solvent. Given the key role of viscosity, an accurate viscosity model suitable for use with reservoir and process simulators is essential. While there are several viscosity models for natural gases and conventional oils, a compositional model applicable to heavy petroleum and diluents is lacking. The objective of this thesis is to develop a general compositional viscosity model that is applicable to natural gas mixtures, conventional crudes oils, heavy petroleum fluids, and their mixtures with solvents and other crudes. The recently developed Expanded Fluid (EF) viscosity correlation was selected as a suitable compositional viscosity model for petroleum applications. The correlation relates the viscosity of the fluid to its density over a broad range of pressures and temperatures. The other inputs are pressure and the dilute gas viscosity. Each fluid is characterized for the correlation by a set of fluid-specific parameters which are tuned to fit data. First, the applicability of the EF correlation was extended to asymmetric mixtures and liquid mixtures containing dissolved gas components. A new set of mass-fraction based mixing rules was developed to calculate the fluid-specific parameters for mixtures. The EF correlation with the new set of mixing rules predicted the viscosity of over 100 mixtures of hydrocarbon compounds and carbon dioxide with overall average absolute relative deviations (AARD) of less than 10% either with measured densities or densities estimated by Advanced Peng-Robinson equation of state (APR EoS). To improve the viscosity predictions with APR EoS-estimated densities, general correlations were developed for non-zero viscosity binary interaction parameters. The EF correlation was extended to non-hydrocarbon compounds typically encountered in natural gas industry. It was demonstrated that the framework of the correlation is valid for these compounds, except for compounds with strong hydrogen bonding such as water. A temperature dependency was introduced into the correlation for strongly hydrogen bonding compounds. The EF correlation fit the viscosity data of pure non-hydrocarbon compounds with AARDs below 6% and predicted the viscosity of sour and sweet natural gases and aqueous solutions of organic alcohols with overall AARDs less than 9%. An internally consistent estimation method was also developed to calculate the fluid-specific parameters for hydrocarbons when no experimental viscosity data are available. The method correlates the fluid-specific parameters to the molecular weight and specific gravity. The method was evaluated against viscosity data of over 250 pure hydrocarbon compounds and petroleum distillations cuts. The EF correlation predictions were found to be within the same order of magnitude of the measurements with an overall AARD of 31%. A methodology was then proposed to apply the EF viscosity correlation to crude oils characterized as mixtures of the defined components and pseudo-components. The above estimation methods are used to calculate the fluid-specific parameters for pseudo-components. Guidelines are provided for tuning of the correlation to available viscosity data, calculating the dilute gas viscosities, and improving the densities calculated with the Peng-Robinson EoS. The viscosities of over 10 dead and live crude oils and bitumen were predicted within a factor of 3 of the measured values using the measured density of the oils as the input. It was shown that single parameter tuning of the model improved the viscosity prediction to within 30% of the measured values. Finally, the performance of the EF correlation was evaluated for diluted heavy oils and bitumens. The required density and viscosity data were collected for over 20 diluted dead and live bitumen mixtures using an in-house capillary viscometer also equipped with an in-line density-meter at temperatures and pressures up to 175 °C and 10 MPa. The predictions of the correlation were found within the same order of magnitude of the measured values with overall AARDs less than 20%. It was shown that the predictions of the correlation with generalized non-zero interaction parameters for the solvent-oil pairs were improved to overall AARDs less than 10%.

Viscosity models for pure hydrocarbons at extreme conditions: A review and comparative study

DOE PAGES

Baled, Hseen O.; Gamwo, Isaac K.; Enick, Robert M.; ...

2018-01-12

Here, viscosity is a critical fundamental property required in many applications in the chemical and oil industries. In this review the performance of seven select viscosity models, representative of various predictive and correlative approaches, is discussed and evaluated by comparison to experimental data of 52 pure hydrocarbons including straight-chain alkanes, branched alkanes, cycloalkanes, and aromatics. This analysis considers viscosity data to extremely high-temperature, high-pressure conditions up to 573 K and 300 MPa. Unsatisfactory results are found, particularly at high pressures, with the Chung-Ajlan-Lee-Starling, Pedersen-Fredenslund, and Lohrenz-Bray-Clark models commonly used for oil reservoir simulation. If sufficient experimental viscosity data are readilymore » available to determine model-specific parameters, the free volume theory and the expanded fluid theory models provide generally comparable results that are superior to those obtained with the friction theory, particularly at pressures higher than 100 MPa. Otherwise, the entropy scaling method by Lötgering-Lin and Gross is recommended as the best predictive model.« less

Viscosity models for pure hydrocarbons at extreme conditions: A review and comparative study

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

Baled, Hseen O.; Gamwo, Isaac K.; Enick, Robert M.

Here, viscosity is a critical fundamental property required in many applications in the chemical and oil industries. In this review the performance of seven select viscosity models, representative of various predictive and correlative approaches, is discussed and evaluated by comparison to experimental data of 52 pure hydrocarbons including straight-chain alkanes, branched alkanes, cycloalkanes, and aromatics. This analysis considers viscosity data to extremely high-temperature, high-pressure conditions up to 573 K and 300 MPa. Unsatisfactory results are found, particularly at high pressures, with the Chung-Ajlan-Lee-Starling, Pedersen-Fredenslund, and Lohrenz-Bray-Clark models commonly used for oil reservoir simulation. If sufficient experimental viscosity data are readilymore » available to determine model-specific parameters, the free volume theory and the expanded fluid theory models provide generally comparable results that are superior to those obtained with the friction theory, particularly at pressures higher than 100 MPa. Otherwise, the entropy scaling method by Lötgering-Lin and Gross is recommended as the best predictive model.« less

Effects of fluid viscosity on a moving sonoluminescing bubble.

PubMed

Sadighi-Bonabi, Rasoul; Mirheydari, Mona; Rezaee, Nastaran; Ebrahimi, Homa

2011-08-01

Based on the quasi-adiabatic model, the parameters of the bubble interior for a moving single bubble sonoluminescence in water, adiponitrile, and N-methylformamide are calculated for various fluid viscosities. By using a complete form of the hydrodynamic force, the bubble trajectory is calculated for a moving single bubble sonoluminescence (m-SBSL). It is found that as the fluid viscosity increases, the unique circular path changes to an ellipsoidal and then linear form and along this incrementally increase of viscosity the light intensity increases. By using the Bremsstrahlung model to describe the bubble radiation, gradual increase of the viscosity results in brighter emissions. It is found that in fluids with higher viscosity the light intensity decreases as time passes.

Rheology and tribology of lubricants with polymeric viscosity modifiers

NASA Astrophysics Data System (ADS)

Babak, LotfizadehDehkordi

Elastohydrodynamic lubrication (EHL) theory has been used to model the lubrication state of antifriction machine elements, where initial viscosity and pressure viscosity coefficients are essential parameters in film thickness modeling. Since the pressures of lubricants in the contact zone can be very high, it is important to know the rheological properties of lubricants in these pressure and temperature regimes. The characteristics of viscosity behavior as a function of pressure are also essential for a universal definition of the pressure viscosity coefficient in order to estimate film thickness in an EHL regime. In this study, viscosities and pressure-viscosity coefficients of ten commercial engine and gear oils and seventeen laboratory-produced oil/polymer viscosity modifiers (VM) additives are measured up to 1.3 GPa at 40, 75 and 100 °C. For the first time, a sharp increase in the viscosity and piezoviscous factor is observed in both mineral-based and synthetic-based oils with different VMs. Analysis of the experimental results indicates that sharp increase in viscosity observed in these experiments are believed to arise from physical changes in the VMs, that is liquid-solid phase transition. Evidence is offered that polymer properties such as molecular weight, concentration and structure influence the onset of the phase transitions. A modified Yasutomi model, which normally describes the pressure dependence of the viscosity of lubricants very well, fails to predict the viscosity of the specimens above the onset of sharp increase in viscosity. A design of experiment (DOE) analysis using Design-Expert software indicates that pressure and temperature are the most critical parameters in the viscosity variation. Tribological tests demonstrate that wear in the contact, zone occurs at temperatures and stresses that coincides with the VM phase transitions in both commercial and laboratory synthesized oil/VMs. Tribological results also indicate that the onset of the sharp increase in viscosity can have significant and unanticipated consequences on the elastohydrodynamic contact and can adversely affect EHL theory. The onset of the steep rise in viscosity may also affect the torque and power losses in a mechanical system. Hence, this previously unknown behavior of the lubricant with VMs should be seriously considered in the application of lubricant in mechanical system.

New Numerical Approaches To thermal Convection In A Compositionally Stratified Fluid

NASA Astrophysics Data System (ADS)

Puckett, E. G.; Turcotte, D. L.; Kellogg, L. H.; Lokavarapu, H. V.; He, Y.; Robey, J.

2016-12-01

Seismic imaging of the mantle has revealed large and small scale heterogeneities in the lower mantle; specifically structures known as large low shear velocity provinces (LLSVP) below Africa and the South Pacific. Most interpretations propose that the heterogeneities are compositional in nature, differing from the overlying mantle, an interpretation that would be consistent with chemical geodynamic models. The LLSVP's are thought to be very old, meaning they have persisted thoughout much of Earth's history. Numerical modeling of persistent compositional interfaces present challenges to even state-of-the-art numerical methodology. It is extremely difficult to maintain sharp composition boundaries which migrate and distort with time dependent fingering without compositional diffusion and / or artificial diffusion. The compositional boundary must persist indefinitely. In this work we present computations of an initial compositionally stratified fluid that is subject to a thermal gradient ΔT = T1 - T0 across the height D of a rectangular domain over a range of buoyancy numbers B and Rayleigh numbers Ra. In these computations we compare three numerical approaches to modeling the movement of two distinct, thermally driven, compositional fields; namely, a high-order Finte Element Method (FEM) that employs artifical viscosity to preserve the maximum and minimum values of the compositional field, a Discontinous Galerkin (DG) method with a Bound Preserving (BP) limiter, and a Volume-of-Fluid (VOF) interface tracking algorithm. Our computations demonstrate that the FEM approach has far too much numerical diffusion to yield meaningful results, the DGBP method yields much better resuts but with small amounts of each compositional field being (numerically) entrained within the other compositional field, while the VOF method maintains a sharp interface between the two compositions throughout the computation. In the figure we show a comparison of between the three methods for a computation made with B = 1.111 and Ra = 10,000 after the flow has reached 'steady state'. (R) the images computed with the standard FEM method (with artifical viscosity), (C) the images computed with the DGBP method (with no artifical viscosity or diffusion due to discretization errors) and (L) the images computed with the VOF algorithm.

Laparoscopy and tribology: the effect of laparoscopic gas on peritoneal fluid.

PubMed

Ott, D E

2001-02-01

To assess the changes in viscosity of peritoneal fluid during laparoscopic exposure to CO2 insufflation. Analysis and mathematic modeling of peritoneal fluid viscosity in vivo and in vitro as a result of exposure to unconditioned CO2 (Canadian Task Force classification II-2). Medical school university research laboratory and hospital. Peritoneal fluid from 45 women. Peritoneal fluid was obtained at laparoscopy before insufflation and tested for viscosity after exposure to currently used raw dry unconditioned CO2. Peritoneal fluid viscosity was tested by viscometric methods and mathematic modeling. Initial viscosity of peritoneal fluid before gas exposure was 1.425 centipoise (cP). Viscosity measurements were obtained at 20-second intervals for gas flows of 1 and 3 L/minute. Increases in viscosity occur rapidly, and by 200 seconds it was 59 cP and 98 cP for 1 and 3 L flow rates, respectively. Very dry CO2 for laparoscopy causes peritoneal fluid viscosity to increase dramatically. (J Am Assoc Gynecol Laparosc 8(1):117-123, 2001)

Sweetness and other sensory properties of model fruit drinks: Does viscosity have an impact?

PubMed

Brandenstein, Cai V S; Busch-Stockfisch, Mechthild; Fischer, Markus

2015-03-15

The impact of thickening agents and viscosity levels on sensory perception was studied in model fruit drinks. Four formulations were prepared that varied in the sweetener blend (erythritol, maltitol and/or steviol glycosides). Locust bean gum and its blends with either xanthan or carrageenan were used to adjust viscosity levels (20, 40, and 70 mPa s). The ranges of viscosity and sweetness level were selected to represent a typical concentration range in commercially available beverages. An increase in viscosity resulted in significant increases in pulpiness, sliminess and perceived viscosity (P-values ≤ 0.001), which were not dependent on sweeteners or hydrocolloid type. Taste perception remained largely unchanged irrespective of the hydrocolloid used. The impact of viscosity on sweetness and taste perception was much smaller in the concentrations used than has been generally reported. The effect of the type of hydrocolloid on the perception of taste attributes was greater than that of viscosity. © 2014 Society of Chemical Industry.

Secular rotational motions and the mechanical structure of a dynamical viscoelastic earth

NASA Technical Reports Server (NTRS)

Yuen, D. A.; Sabadini, R.

1984-01-01

A survey is presented of analytical methods for computing the linear responses of the rotational axis of a layered viscoelastic earth to surface loading. Theoretical research in this area is first summarized, and the differences between the mechanical boundary conditions to be applied at the interface separating the upper and lower mantles for an adiabatically and chemically stratified mantle are discussed. Some examples of polar wander and secular variation of the spin rate from glacial excitation are presented for various types of chemical and viscosity stratifications. The effects of an artificial density jump at the base of the lithosphere in models are examined, and certain issues concerning the fluid tidal Love number for different types of density stratification are addressed. The meaning of effective plate thickness over geological time scales for rotational dynamics is discussed.

Compositional dependence of lower crustal viscosity

NASA Astrophysics Data System (ADS)

Shinevar, William J.; Behn, Mark D.; Hirth, Greg

2015-10-01

We calculate the viscosity structure of the lower continental crust as a function of its bulk composition using multiphase mixing theory. We use the Gibbs free-energy minimization routine Perple_X to calculate mineral assemblages for different crustal compositions under pressure and temperature conditions appropriate for the lower continental crust. The effective aggregate viscosities are then calculated using a rheologic mixing model and flow laws for the major crust-forming minerals. We investigate the viscosity of two lower crustal compositions: (i) basaltic (53 wt % SiO2) and (ii) andesitic (64 wt % SiO2). The andesitic model predicts aggregate viscosities similar to feldspar and approximately 1 order of magnitude greater than that of wet quartz. The viscosity range calculated for the andesitic crustal composition (particularly when hydrous phases are stable) is most similar to independent estimates of lower crust viscosity in actively deforming regions based on postglacial isostatic rebound, postseismic relaxation, and paleolake shoreline deflection.

Fluid friction and wall viscosity of the 1D blood flow model.

PubMed

Wang, Xiao-Fei; Nishi, Shohei; Matsukawa, Mami; Ghigo, Arthur; Lagrée, Pierre-Yves; Fullana, Jose-Maria

2016-02-29

We study the behavior of the pulse waves of water into a flexible tube for application to blood flow simulations. In pulse waves both fluid friction and wall viscosity are damping factors, and difficult to evaluate separately. In this paper, the coefficients of fluid friction and wall viscosity are estimated by fitting a nonlinear 1D flow model to experimental data. In the experimental setup, a distensible tube is connected to a piston pump at one end and closed at another end. The pressure and wall displacements are measured simultaneously. A good agreement between model predictions and experiments was achieved. For amplitude decrease, the effect of wall viscosity on the pulse wave has been shown as important as that of fluid viscosity. Copyright © 2016 Elsevier Ltd. All rights reserved.

Hybrid modeling of microbial exopolysaccharide (EPS) production: The case of Enterobacter A47.

PubMed

Marques, Rodolfo; von Stosch, Moritz; Portela, Rui M C; Torres, Cristiana A V; Antunes, Sílvia; Freitas, Filomena; Reis, Maria A M; Oliveira, Rui

2017-03-20

Enterobacter A47 is a bacterium that produces high amounts of a fucose-rich exopolysaccharide (EPS) from glycerol residue of the biodiesel industry. The fed-batch process is characterized by complex non-linear dynamics with highly viscous pseudo-plastic rheology due to the accumulation of EPS in the culture medium. In this paper, we study hybrid modeling as a methodology to increase the predictive power of models for EPS production optimization. We compare six hybrid structures that explore different levels of knowledge-based and machine-learning model components. Knowledge-based components consist of macroscopic material balances, Monod type kinetics, cardinal temperature and pH (CTP) dependency and power-law viscosity models. Unknown dependencies are set to be identified by a feedforward artificial neural network (ANN). A semiparametric identification schema is applied resorting to a data set of 13 independent fed-batch experiments. A parsimonious hybrid model was identified that describes the dynamics of the 13 experiments with the same parameterization. The final model is specific to Enterobacter A47 but can be easily extended to other microbial EPS processes. Copyright © 2017 Elsevier B.V. All rights reserved.

Comparison of observed rheological properties of hard wheat flour dough with predictions of the Giesekus-Leonov, White-Metzner and Phan-Thien Tanner models

NASA Technical Reports Server (NTRS)

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

1999-01-01

The measured rheological behavior of hard wheat flour dough was predicted using three nonlinear differential viscoelastic models. The Phan-Thien Tanner model gave good zero shear viscosity prediction, but overpredicted the shear viscosity at higher shear rates and the transient and extensional properties. The Giesekus-Leonov model gave similar predictions to the Phan-Thien Tanner model, but the extensional viscosity prediction showed extension thickening. Using high values of the mobility factor, extension thinning behavior was observed but the predictions were not satisfactory. The White-Metzner model gave good predictions of the steady shear viscosity and the first normal stress coefficient but it was unable to predict the uniaxial extensional viscosity as it exhibited asymptotic behavior in the tested extensional rates. It also predicted the transient shear properties with moderate accuracy in the transient phase, but very well at higher times, compared to the Phan-Thien Tanner model and the Giesekus-Leonov model. None of the models predicted all observed data consistently well. Overall the White-Metzner model appeared to make the best predictions of all the observed data.

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.

Basics of Physical Modeling in Coastal and Hydraulic Engineering

DTIC Science & Technology

2013-09-01

gravity (Fg), viscosity (Fv), surface tension (Fs), and elasticity (Fe) must have the same ratios. This requirement arises from Newton’s Second Law which...they are relatively small. Viscosity can be neglected in most free-surface models if the model is not too (a) (b) ERDC/CHL CHETN-XIII-3 September... viscosity is to ensure that the model flow is in the turbulent range, which occurs for Re above approximately 104. The Reynolds number (Re) is defined

Capillary waves with surface viscosity

NASA Astrophysics Data System (ADS)

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

2017-11-01

Experiments over the last 50 years have suggested a correlation between the surface (shear) viscosity and the stability of a foam or emulsion. With recent techniques allowing more accurate measurements of the elusive surface viscosity, we examine this link theoretically using small-amplitude capillary waves in the presence of the Marangoni effect and surface viscosity modelled via the Boussinesq-Scriven model. The surface viscosity effect is found to contribute a damping effect on the amplitude of the capillary wave with subtle differences to the effect of the convective-diffusive Marangoni transport. The general wave dispersion is augmented to take into account the Marangoni and surface viscosity effects, and a first-order correction to the critical damping wavelength is derived. The authors acknowledge the financial support of the Shell University Technology Centre for fuels and lubricants.

Mantle viscosity structure constrained by joint inversions of seismic velocities and density

NASA Astrophysics Data System (ADS)

Rudolph, M. L.; Moulik, P.; Lekic, V.

2017-12-01

The viscosity structure of Earth's deep mantle affects the thermal evolution of Earth, the ascent of mantle upwellings, sinking of subducted oceanic lithosphere, and the mixing of compositional heterogeneities in the mantle. Modeling the long-wavelength dynamic geoid allows us to constrain the radial viscosity profile of the mantle. Typically, in inversions for the mantle viscosity structure, wavespeed variations are mapped into density variations using a constant- or depth-dependent scaling factor. Here, we use a newly developed joint model of anisotropic Vs, Vp, density and transition zone topographies to generate a suite of solutions for the mantle viscosity structure directly from the seismologically constrained density structure. The density structure used to drive our forward models includes contributions from both thermal and compositional variations, including important contributions from compositionally dense material in the Large Low Velocity Provinces at the base of the mantle. These compositional variations have been neglected in the forward models used in most previous inversions and have the potential to significantly affect large-scale flow and thus the inferred viscosity structure. We use a transdimensional, hierarchical, Bayesian approach to solve the inverse problem, and our solutions for viscosity structure include an increase in viscosity below the base of the transition zone, in the shallow lower mantle. Using geoid dynamic response functions and an analysis of the correlation between the observed geoid and mantle structure, we demonstrate the underlying reason for this inference. Finally, we present a new family of solutions in which the data uncertainty is accounted for using covariance matrices associated with the mantle structure models.

Stereoisomeric effects on dynamic viscosity versus pressure and temperature for the system cis- + trans-decalin

NASA Astrophysics Data System (ADS)

Miyake, Yasufumi; Boned, Christian; Baylaucq, Antoine; Bessières, David; Zéberg-Mikkelsen, Claus K.; Galliéro, Guillaume; Ushiki, Hideharu

2007-07-01

In order to study the influence of stereoisomeric effects on the dynamic viscosity, an extensive experimental study of the viscosity of the binary system composed of the two stereoisomeric molecular forms of decalin - cis and trans - has been carried out for five different mixtures at three temperatures (303.15, 323.15 and 343.15) K and six isobars up to 100 MPa with a falling-body viscometer (a total of 90 points). The experimental relative uncertainty is estimated to be 2%. The variations of dynamic viscosity versus composition are discussed with respect to their behavior due to stereoisomerism. Four different models with a physical and theoretical background are studied in order to investigate how they take the stereoisomeric effect into account through their required model parameters. The evaluated models are based on the hard-sphere scheme, the concepts of the free-volume and the friction theory, and a model derived from molecular dynamics. Overall, a satisfactory representation of the viscosity of this binary system is found for the different models within the considered ( T, p) range taken into account their simplicity. All the models are able to distinguish between the two stereoisomeric decalin compounds. Further, based on the analysis of the model parameters performed on the pure compounds, it has been found that the use of simple mixing rules without introducing any binary interaction parameters are sufficient in order to predict the viscosity of cis + trans-decalin mixtures with the same accuracy in comparison with the experimental values as obtained for the pure compounds. In addition to these models, a semi-empirical self-referencing model and the simple mixing laws of Grunberg-Nissan and Katti-Chaudhri are also applied in the representation of the viscosity behavior of these systems.

Development of a Physiologically Based Computational Kidney Model to Describe the Renal Excretion of Hydrophilic Agents in Rats

PubMed Central

Niederalt, Christoph; Wendl, Thomas; Kuepfer, Lars; Claassen, Karina; Loosen, Roland; Willmann, Stefan; Lippert, Joerg; Schultze-Mosgau, Marcus; Winkler, Julia; Burghaus, Rolf; Bräutigam, Matthias; Pietsch, Hubertus; Lengsfeld, Philipp

2013-01-01

A physiologically based kidney model was developed to analyze the renal excretion and kidney exposure of hydrophilic agents, in particular contrast media, in rats. In order to study the influence of osmolality and viscosity changes, the model mechanistically represents urine concentration by water reabsorption in different segments of kidney tubules and viscosity dependent tubular fluid flow. The model was established using experimental data on the physiological steady state without administration of any contrast media or drugs. These data included the sodium and urea concentration gradient along the cortico-medullary axis, water reabsorption, urine flow, and sodium as well as urea urine concentrations for a normal hydration state. The model was evaluated by predicting the effects of mannitol and contrast media administration and comparing to experimental data on cortico-medullary concentration gradients, urine flow, urine viscosity, hydrostatic tubular pressures and single nephron glomerular filtration rate. Finally the model was used to analyze and compare typical examples of ionic and non-ionic monomeric as well as non-ionic dimeric contrast media with respect to their osmolality and viscosity. With the computational kidney model, urine flow depended mainly on osmolality, while osmolality and viscosity were important determinants for tubular hydrostatic pressure and kidney exposure. The low diuretic effect of dimeric contrast media in combination with their high intrinsic viscosity resulted in a high viscosity within the tubular fluid. In comparison to monomeric contrast media, this led to a higher increase in tubular pressure, to a reduction in glomerular filtration rate and tubular flow and to an increase in kidney exposure. The presented kidney model can be implemented into whole body physiologically based pharmacokinetic models and extended in order to simulate the renal excretion of lipophilic drugs which may also undergo active secretion and reabsorption. PMID:23355822

The Larger the Viscosity, the Higher the Bounce

NASA Astrophysics Data System (ADS)

Stern, Menachem; Klein Schaarsberg, Martin; Peters, Ivo; Dodge, Kevin; Zhang, Wendy; Jaeger, Heinrich

A low-viscosity liquid drop can bounce upon impact onto a solid. A high-viscosity drop typically just flattens, i.e., it splats. Surprisingly, our experiments with a droplet made of densely packed glass beads in silicone oil display the opposite behavior: the low-viscosity oil suspension drop splats. The high-viscosity oil suspension bounces. Increasing solvent viscosity increases the rebound energy. To gain insight into the underlying mechanism, we model the suspension as densely packed elastic spheres experiencing viscous lubrication drag between neighbors. The model reproduces the observed trends. Plots of elastic compression and drag experienced by the particles show that rebounds are made possible by (1) a fraction of the impact energy being stored during initial contact via elastic compression, (2) a rapid broadening of local lubrication drag interactions at the initial impact site into a spatially uniform upward force throughout the drop. Including finite wall drag due to the presence of ambient air into the numerical model diminishes and eventually cuts off the rebound.

A dynamic model of Venus's gravity field

NASA Technical Reports Server (NTRS)

Kiefer, W. S.; Richards, M. A.; Hager, B. H.; Bills, B. G.

1984-01-01

Unlike Earth, long wavelength gravity anomalies and topography correlate well on Venus. Venus's admittance curve from spherical harmonic degree 2 to 18 is inconsistent with either Airy or Pratt isostasy, but is consistent with dynamic support from mantle convection. A model using whole mantle flow and a high viscosity near surface layer overlying a constant viscosity mantle reproduces this admittance curve. On Earth, the effective viscosity deduced from geoid modeling increases by a factor of 300 from the asthenosphere to the lower mantle. These viscosity estimates may be biased by the neglect of lateral variations in mantle viscosity associated with hot plumes and cold subducted slabs. The different effective viscosity profiles for Earth and Venus may reflect their convective styles, with tectonism and mantle heat transport dominated by hot plumes on Venus and by subducted slabs on Earth. Convection at degree 2 appears much stronger on Earth than on Venus. A degree 2 convective structure may be unstable on Venus, but may have been stabilized on Earth by the insulating effects of the Pangean supercontinental assemblage.

A fast, time-accurate unsteady full potential scheme

NASA Technical Reports Server (NTRS)

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

1985-01-01

The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. 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 by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.

Numerical simulation of exploding pusher targets

NASA Astrophysics Data System (ADS)

Atzeni, S.; Rosenberg, M. J.; Gatu Johnson, M.; Petrasso, R. D.

2017-10-01

Exploding pusher targets, i.e. gas-filled large aspect-ratio glass or plastic shells, driven by a strong laser-generated shock, are widely used as pulsed sources of neutrons and fast charged particles. Recent experiments on exploding pushers provided evidence for the transition from a purely fluid behavior to a kinetic one. Indeed, fluid models largely overpredict yield and temperature as the Knudsen number Kn (ratio of ion mean-free path to compressed gas radius) is comparable or larger than one. At Kn = 0.3 - 1, fluid codes reasonably estimate integral quantities as yield and neutron-averaged temperatures, but do not reproduce burn radii, burn profiles and DD/DHe3 yield ratio. This motivated a detailed simulation study of intermediate-Kn exploding pushers. We will show how simulation results depend on models for laser-interaction, electron conductivity (flux-limited local vs nonlocal), viscosity (physical vs artificial), and ion mixing. Work partially supported by Sapienza Project C26A15YTMA, Sapienza 2016 (n. 257584), and Eurofusion Project AWP17-ENR-IFE-CEA-01.

Tethered fleximags: a physical model for ciliary propulsion.

NASA Astrophysics Data System (ADS)

Du Roure, Olivia; Babataheri, Avin; Jenffer, Patrice; Fermigier, Marc; Goubault, Cecile

2007-11-01

Fleximags are linear colloidal structures made of micron-sized superparamagnetic particles. Permanent links between colloids are established through molecules grafted on the particles. The elasticity of the linker bestows a flexibility to the filament. The fleximags have already been used to make one of the first artificial microswimmers (Dreyfus et al. Nature 2005) resembling a spermatozoon. They can also be anchored to a glass substrate isolated or as arrays. Those arrays build up experimental models of the array of cilia on paramecium for studying physical aspects of the propulsion. Here we'll show our first studies concerning anchored flexiamgs submitted to time-dependent field. The actuation is controlled by three electromagnets and allowing all types of 3-D movements: (a)symmetric beating in a plane, rotation... We first study one single anchored fleximag when the field is rotating on a cone. Only a part of the filament is moving reflecting the competition between magnetic interactions, elasticity and viscosity. The length of this mobile fraction decreases with frequency. We also study the induced flow by PIV.

Control Theory based Shape Design for the Incompressible Navier-Stokes Equations

NASA Astrophysics Data System (ADS)

Cowles, G.; Martinelli, L.

2003-12-01

A design method for shape optimization in incompressible turbulent viscous flow has been developed and validated for inverse design. The gradient information is determined using a control theory based algorithm. With such an approach, the cost of computing the gradient is negligible. An additional adjoint system must be solved which requires the cost of a single steady state flow solution. Thus, this method has an enormous advantage over traditional finite-difference based algorithms. The method of artificial compressibility is utilized to solve both the flow and adjoint systems. An algebraic turbulence model is used to compute the eddy viscosity. The method is validated using several inverse wing design test cases. In each case, the program must modify the shape of the initial wing such that its pressure distribution matches that of the target wing. Results are shown for the inversion of both finite thickness wings as well as zero thickness wings which can be considered a model of yacht sails.

Reconciling postseismic and interseismic surface deformation around strike-slip faults: Earthquake-cycle models with finite ruptures and viscous shear zones

NASA Astrophysics Data System (ADS)

Hearn, E. H.

2013-12-01

Geodetic surface velocity data show that after an energetic but brief phase of postseismic deformation, surface deformation around most major strike-slip faults tends to be localized and stationary, and can be modeled with a buried elastic dislocation creeping at or near the Holocene slip rate. Earthquake-cycle models incorporating an elastic layer over a Maxwell viscoelastic halfspace cannot explain this, even when the earliest postseismic deformation is ignored or modeled (e.g., as frictional afterslip). Models with heterogeneously distributed low-viscosity materials or power-law rheologies perform better, but to explain all phases of earthquake-cycle deformation, Burgers viscoelastic materials with extreme differences between their Maxwell and Kelvin element viscosities seem to be required. I present a suite of earthquake-cycle models to show that postseismic and interseismic deformation may be reconciled for a range of lithosphere architectures and rheologies if finite rupture length is taken into account. These models incorporate high-viscosity lithosphere optionally cut by a viscous shear zone, and a lower-viscosity mantle asthenosphere (all with a range of viscoelastic rheologies and parameters). Characteristic earthquakes with Mw = 7.0 - 7.9 are investigated, with interseismic intervals adjusted to maintain the same slip rate (10, 20 or 40 mm/yr). I find that a high-viscosity lower crust/uppermost mantle (or a high viscosity per unit width viscous shear zone at these depths) is required for localized and stationary interseismic deformation. For Mw = 7.9 characteristic earthquakes, the shear zone viscosity per unit width in the lower crust and uppermost mantle must exceed about 10^16 Pa s /m. For a layered viscoelastic model the lower crust and uppermost mantle effective viscosity must exceed about 10^20 Pa s. The range of admissible shear zone and lower lithosphere rheologies broadens considerably for faults producing more frequent but smaller characteristic earthquakes. Thus, minimum lithosphere or shear zone effective viscosities inferred from interseismic GPS data and infinite-fault earthquake-cycle models may be too high. The finite-fault models show that relaxation of viscoelastic material in the mid crust (most likely along a viscous shear zone) may be consistent with near- to intermediate-field postseismic deformation typical of recent Mw = 7.4 to 7.9 earthquakes. This deformation is compatible with more localized and time-invariant deformation during most of the interseismic interval if (1) shear zone viscosity per unit width increases with depth or (2) the shear zone material has a Burgers viscoelastic rheology.

ΛCDM model with dissipative nonextensive viscous dark matter

NASA Astrophysics Data System (ADS)

Gimenes, H. S.; Viswanathan, G. M.; Silva, R.

2018-03-01

Many models in cosmology typically assume the standard bulk viscosity. We study an alternative interpretation for the origin of the bulk viscosity. Using nonadditive statistics proposed by Tsallis, we propose a bulk viscosity component that can only exist by a nonextensive effect through the nonextensive/dissipative correspondence (NexDC). In this paper, we consider a ΛCDM model for a flat universe with a dissipative nonextensive viscous dark matter component, following the Eckart theory of bulk viscosity, without any perturbative approach. In order to analyze cosmological constraints, we use one of the most recent observations of Type Ia Supernova, baryon acoustic oscillations and cosmic microwave background data.

Towards adjoint-based inversion of time-dependent mantle convection with nonlinear viscosity

NASA Astrophysics Data System (ADS)

Li, Dunzhu; Gurnis, Michael; Stadler, Georg

2017-04-01

We develop and study an adjoint-based inversion method for the simultaneous recovery of initial temperature conditions and viscosity parameters in time-dependent mantle convection from the current mantle temperature and historic plate motion. Based on a realistic rheological model with temperature-dependent and strain-rate-dependent viscosity, we formulate the inversion as a PDE-constrained optimization problem. The objective functional includes the misfit of surface velocity (plate motion) history, the misfit of the current mantle temperature, and a regularization for the uncertain initial condition. The gradient of this functional with respect to the initial temperature and the uncertain viscosity parameters is computed by solving the adjoint of the mantle convection equations. This gradient is used in a pre-conditioned quasi-Newton minimization algorithm. We study the prospects and limitations of the inversion, as well as the computational performance of the method using two synthetic problems, a sinking cylinder and a realistic subduction model. The subduction model is characterized by the migration of a ridge toward a trench whereby both plate motions and subduction evolve. The results demonstrate: (1) for known viscosity parameters, the initial temperature can be well recovered, as in previous initial condition-only inversions where the effective viscosity was given; (2) for known initial temperature, viscosity parameters can be recovered accurately, despite the existence of trade-offs due to ill-conditioning; (3) for the joint inversion of initial condition and viscosity parameters, initial condition and effective viscosity can be reasonably recovered, but the high dimension of the parameter space and the resulting ill-posedness may limit recovery of viscosity parameters.

The Drainage of Thin, Vertical, Model Polyurethane Liquid Films

NASA Astrophysics Data System (ADS)

Snow, Steven; Pernisz, Udo; Braun, Richard; Naire, Shailesh

1999-11-01

We have successfully measured the drainage rate of thin, vertically-aligned, liquid films prepared from model polyurethane foam formulations. The pattern of interference fringes in these films was consistent with a wedge-shaped film profile. The time evolution of this wedge shape (the ``collapsing wedge") obeyed a power law relationship between fringe density s and time t of s = k t^m. Experimentally, m ranged from -0.47 to -0.92. The lower bound for m represented a case where the surface viscosity of the film was very high (a ``rigid" surface). Theoretical modeling of this case yielded m = -0.5, in excellent agreement with experiment. Instantaneous film drainage rate (dV/dt) could be extracted from the ``Collapsing Wedge" model. As expected, dV/dt scaled inversely with bulk viscosity. As surfactant concentration was varied at constant bulk viscosity, dV/dt passed through a maximum value, consistent with a model where the rigidity of the surface was a function of both the intensity of surface tension gradients and the surface viscosity of the film. The influence of surface viscosity on dV/dt was also modeled theoretically.

Role of rheology in reconstructing slab morphology in global mantle models

NASA Astrophysics Data System (ADS)

Bello, Léa; Coltice, Nicolas; Tackley, Paul; Müller, Dietmar

2015-04-01

Reconstructing the 3D structure of the Earth's mantle has been a challenge for geodynamicists for about 40 years. Although numerical models and computational capabilities have incredibly progressed, parameterizations used for modeling convection forced by plate motions are far from being Earth-like. Among the set of parameters, rheology is fundamental because it defines in a non-linear way the dynamics of slabs and plumes, and the organization of the lithosphere. Previous studies have employed diverse viscosity laws, most of them being temperature and depth dependent with relatively small viscosity contrasts. In this study, we evaluate the role of the temperature dependence of viscosity (variations up to 6 orders of magnitude) on reconstructing slab evolution in 3D spherical models of convection driven by plate history models. We also investigate the importance of pseudo-plasticity in such models. We show that strong temperature dependence of viscosity combined with pseudo-plasticity produce laterally and vertically continuous slabs, and flat subduction where trench retreat is fast (North, Central and South America). Moreover, pseudo-plasticity allows a consistent coupling between imposed plate motions and global convection, which is not possible with temperature-dependent viscosity only. However, even our most sophisticated model is not able to reproduce unambiguously stagnant slabs probably because of the simplicity of material properties we use here. The differences between models employing different viscosity laws are very large, larger than the differences between two models with the same rheology but using two different plate reconstructions or initial conditions.

A study on the dependence of nuclear viscosity on temperature

NASA Astrophysics Data System (ADS)

Vardaci, E.; Di Nitto, A.; Nadtochy, P. N.; La Rana, G.; Cinausero, M.; Prete, G.; Gelli, N.; Ashaduzzaman, M.; Davide, F.; Pulcini, A.; Quero, D.; Kozulin, E. M.; Knyazheva, G. N.; Itkis, I. M.

2018-05-01

Nuclear viscosity is an irreplaceable ingredient of nuclear fission collective dynamical models. It drives the exchange of energy between the collective variables and the thermal bath of single particle degrees of freedom. Its dependence on the shape and temperature is a matter of controversy. By using systems of intermediate fissility we have demonstrated in a recent study that the viscosity parameters is larger for compact shapes, and decreases for larger deformations of the fissioning system, at variance with the conclusions of the statistical model modified to include empirically viscosity and time scales. In this contribution we propose an experimental scenario to highlight the possible dependence of the viscosity from the temperature.

The relationship between plate velocity and trench viscosity in Newtonian and power-law subduction calculations

NASA Technical Reports Server (NTRS)

King, Scott D.; Hager, Bradford H.

1990-01-01

The relationship between oceanic trench viscosity and oceanic plate velocity is studied using a Newtonian rheology by varying the viscosity at the trench. The plate velocity is a function of the trench viscosity for fixed Rayleigh number and plate/slab viscosity. Slab velocities for non-Newtonian rheology calculations are significantly different from slab velocities from Newtonian rheology calculations at the same effective Rayleigh number. Both models give reasonable strain rates for the slab when compared with estimates of seismic strain rate. Non-Newtonian rheology eliminates the need for imposed weak zones and provides a self-consistent fluid dynamical mechanism for subduction in numerical convection models.

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.

Viscosity and Structure of CaO-SiO2-P2O5-FetO System with Varying P2O5 and FeO Content

NASA Astrophysics Data System (ADS)

Diao, Jiang; Gu, Pan; Liu, De-Man; Jiang, Lu; Wang, Cong; Xie, Bing

2017-10-01

A rotary viscosimeter and Raman spectrum were employed to measure the viscosity and structural information of the CaO-SiO2-P2O5-FetO system at 1673 K. The experimental data have been compared with the calculated results using different viscosity models. It shows that the National Physical Laboratory (NPL) and Pal models fit the CaO-SiO2-P2O5-FeOt system better. With the P2O5 content increasing from 5% to 14%, the viscosity increases from 0.12 Pa s to 0.27 Pa s. With the FeO content increasing from 30% to 40%, the viscosity decreases from 0.21 Pa s to 0.12 Pa s. Increasing FeO content makes the complicated molten melts become simple, and increasing P2O5 content will complicate the molten melts. The linear relation between viscosity and structure parameter Q(Si + P) was obtained by regression analysis. The calculated viscosity by using the optimized NPL and Pal model are almost identical with the fitted values.

Quantitative characterization of the viscosity of a microemulsion

NASA Technical Reports Server (NTRS)

Berg, Robert F.; Moldover, Michael R.; Huang, John S.

1987-01-01

The viscosity of the three-component microemulsion water/decane/AOT has been measured as a function of temperature and droplet volume fraction. At temperatures well below the phase-separation temperature the viscosity is described by treating the droplets as hard spheres suspended in decane. Upon approaching the two-phase region from low temperature, there is a large (as much as a factor of four) smooth increase of the viscosity which may be related to the percolation-like transition observed in the electrical conductivity. This increase in viscosity is not completely consistent with either a naive electroviscous model or a simple clustering model. The divergence of the viscosity near the critical point (39 C) is superimposed upon the smooth increase. The magnitude and temperature dependence of the critical divergence are similar to that seen near the critical points of binary liquid mixtures.

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.

Temperature and pressure dependences of kimberlite melts viscosity (experimental-theoretical study)

NASA Astrophysics Data System (ADS)

Persikov, Eduard; Bykhtiyarov, Pavel; Cokol, Alexsander

2016-04-01

Experimental data on temperature and pressure dependences of viscosity of model kimberlite melts (silicate 82 + carbonate 18, wt. %, 100NBO/T = 313) have been obtained for the first time at 100 MPa of CO2 pressure and at the lithostatic pressures up to 7.5 GPa in the temperature range 1350 oC - 1950 oC using radiation high gas pressure apparatus and press free split-sphere multi - anvil apparatus (BARS). Experimental data obtained on temperature and pressure dependences of viscosity of model kimberlite melts at moderate and high pressures is compared with predicted data on these dependences of viscosity of basaltic melts (100NBO/T = 58) in the same T, P - range. Dependences of the viscosity of model kimberlite and basaltic melts on temperature are consistent to the exponential Arrenian equation in the T, P - range of experimental study. The correct values of activation energies of viscous flow of kimberlite melts have been obtained for the first time. The activation energies of viscous flow of model kimberlite melts exponentially increase with increasing pressure and are equal: E = 130 ± 1.3 kJ/mole at moderate pressure (P = 100 MPa) and E = 160 ± 1.6 kJ/mole at high pressure (P = 5.5 GPa). It has been established too that the viscosity of model kimberlite melts exponentially increases on about half order of magnitude with increasing pressures from 100 MPa to 7.5 GPa at the isothermal condition (1800 oC). It has been established that viscosity of model kimberlite melts at the moderate pressure (100 MPa) is lover on about one order of magnitude to compare with the viscosity of basaltic melts, but at high pressure range (5.5 - 7.5 GPa), on the contrary, is higher on about half order of magnitude at the same values of the temperatures. Here we use both a new experimental data on viscosity of kimberlite melts and our structural chemical model for calculation and prediction the viscosity of magmatic melts [1] to determine the fundamental features of viscosity of kimberlite and basaltic magmas at the T, P - parameters of the Earth's crust and upper mantle. The Russian Foundation for Basic Research (project 15-05-01318) and the Russian Science Foundation (project 14-27-00054) are acknowledged for the financial support. [1] Persikov, E.S. & Bukhtiyarov, P.G. (2009) Russian Geology & Geophysics, 50, No 12, 1079-1090.

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.

Dependence of Perpendicular Viscosity on Magnetic Fluctuations in a Stochastic Topology

NASA Astrophysics Data System (ADS)

Fridström, R.; Chapman, B. E.; Almagri, A. F.; Frassinetti, L.; Brunsell, P. R.; Nishizawa, T.; Sarff, J. S.

2018-06-01

In a magnetically confined plasma with a stochastic magnetic field, the dependence of the perpendicular viscosity on the magnetic fluctuation amplitude is measured for the first time. With a controlled, ˜ tenfold variation in the fluctuation amplitude, the viscosity increases ˜100 -fold, exhibiting the same fluctuation-amplitude-squared dependence as the predicted rate of stochastic field line diffusion. The absolute value of the viscosity is well predicted by a model based on momentum transport in a stochastic field, the first in-depth test of this model.

Viscosity studies of water based magnetite nanofluids

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

Anu, K.; Hemalatha, J.

2016-05-23

Magnetite nanofluids of various concentrations have been synthesized through co-precipitation method. The structural and topographical studies made with the X-Ray Diffractometer and Atomic Force Microscope are presented in this paper. The density and viscosity studies for the ferrofluids of various concentrations have been made at room temperature. The experimental viscosities are compared with theoretical values obtained from Einstein, Batchelor and Wang models. An attempt to modify the Rosensweig model is made and the modified Rosensweig equation is reported. In addition, new empirical correlation is also proposed for predicting viscosity of ferrofluid at various concentrations.

Shear Viscosity Coefficient of 5d Liquid Transition Metals

NASA Astrophysics Data System (ADS)

Thakor, P. B.; Sonvane, Y. A.; Gajjar, P. N.; Jani, A. R.

2011-07-01

In the present paper we have calculated shear viscosity coefficient (η) of 5 d liquid transition metals. To calculate effective pair potential ν(r) and pair distribution function g(r) we have used our own newly constructed model potential and Percus- Yevick hard sphere (PYHS) structure factor S(q) respectively. We have also investigated the effect of different correction function like Hartree (H), Taylor (T) and Sarkar et al. (S) on shear viscosity coefficient (η). Our newly constructed model potential successfully explains the shear viscosity coefficient (η) of 5 d liquid transition metals.

Effect of Qing Nao tablet on blood stasis model of mice

NASA Astrophysics Data System (ADS)

Kong, Xuejun; Hao, Shaojun; Wang, Hongyu; Liu, Xiaobin; Xie, Guoqi; Li, Wenjun; Zhang, Zhengchen

2018-04-01

To investigate the effect of Qing Nao tablet on mouse model of blood stasis syndrome, 60 mice, male and female, were randomly divided into 6 groups, were fed with large, small doses of Qing Nao tablet suspension, Naoluotong saline suspension and the same volume (group 2, 0.1ml/10g), administer 1 times daily, orally for 15 days. Intragastric administration for first days, in addition to the 1 group saline group every day in the hind leg intramuscular saline, the other 5 groups each rat day hind leg muscle injection of dexamethasone 0.8mg/kg intramuscular injection every day, 1 times, 15 days. 1 hour continuous intramuscular injection and intramuscular drug perfusion on the sixteenth day after mice. The eyeball blood, heparin after whole blood viscosity test. Compared with the control group, model group, high and low shear viscosity were significantly increased (P<0.01), indicating that the model was successful. Compared with the model group, high dose group and Qing Nao tablet Naoluotong group can significantly reduce the viscosity at high shear and (P<0.01), middle dose Qing Nao tablet group can significantly reduce high shear and shear viscosity (P<0.05); large, middle dose Qing Nao tablet group can significantly reduce the low shear viscosity (P<0.05), Naoluotong group can significantly reduce the low shear viscosity (P<0.01); low dose Qing Nao tablet group were lower high cut, low shear viscosity and trend The potential (P>0.05). The Qing Nao tablet has a good effect on the model of blood stasis in mice.

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

Viscosity of dilute suspensions of rigid bead arrays at low shear: accounting for the variation in hydrodynamic stress over the bead surfaces.

PubMed

Allison, Stuart A; Pei, Hongxia

2009-06-11

In this work, we examine the viscosity of a dilute suspension of irregularly shaped particles at low shear. A particle is modeled as a rigid array of nonoverlapping beads of variable size and geometry. Starting from a boundary element formalism, approximate account is taken of the variation in hydrodynamic stress over the surface of the individual beads. For a touching dimer of two identical beads, the predicted viscosity is lower than the exact value by 5.2%. The methodology is then applied to several other model systems including tetramers of variable conformation and linear strings of touching beads. An analysis is also carried out of the viscosity and translational diffusion of several dilute amino acids and diglycine in water. It is concluded that continuum hydrodynamic modeling with stick boundary conditions is unable to account for the experimental viscosity and diffusion data simultaneously. A model intermediate between "stick" and "slip" could possibly reconcile theory and experiment.

Gas Accretion onto a Supermassive Black Hole: A Step to Model AGN Feedback

NASA Astrophysics Data System (ADS)

Nagamine, K.; Barai, P.; Proga, D.

2012-08-01

We study gas accretion onto a supermassive black hole (SMBH) using the 3D SPH code GADGET-3 on scales of 0.1-200 pc. First we test our code with the spherically symmetric, adiabatic Bondi accretion problem. We find that our simulation can reproduce the expected Bondi accretion flow very well for a limited amount of time until the effect of the outer boundary starts to be visible. We also find artificial heating of gas near the inner accretion boundary due to the artificial viscosity of SPH. Second, we implement radiative cooling and heating due to X-rays, and examine the impact of thermal feedback by the central X-ray source. The accretion flow roughly follows the Bondi solution for low central X-ray luminosities; however, the flow starts to exhibit non-spherical fragmentation due to the thermal instability for a certain range of central LX, and a strong overall outflow develops for greater LX. The cold gas develops filamentary structures that fall into the central SMBH, whereas the hot gas tries to escape through the channels in between the cold filaments. Such fragmentation of accreting gas can assist in the formation of clouds around AGN, induce star-formation, and contribute to the observed variability of narrow-line regions.

The nonlinear effects of the eddy viscosity and the bottom friction on the Lagrangian residual velocity in a narrow model bay

NASA Astrophysics Data System (ADS)

Deng, Fangjing; Jiang, Wensheng; Feng, Shizuo

2017-09-01

The nonlinear effects of the eddy viscosity and the bottom friction on the Lagrangian residual velocity (LRV) are studied numerically in a narrow model bay. Three groups of the experiments with different eddy viscosity and different forms of the bottom friction are designed and are carried out in the three kinds of the topography. When the eddy viscosity is obtained from a two-equation turbulence closure model, the pattern of the LRV is more complex than that of the time invariant eddy viscosity case and the intensity is from more than 1.3 times to one order smaller than that of the linear eddy viscosity condition. The LRV are also acquired when the eddy viscosity varies from the flood-averaged one to the ebb-averaged one. It is found that when the flood-averaged eddy viscosity is bigger than the ebb-averaged eddy viscosity (flood-dominated asymmetry), the direction of the breadth-averaged LRV and the 3D LRV is nearly opposite to that when the eddy viscosity asymmetry is reverse (ebb-dominated asymmetry). However, the intensity of the LRV for the ebb-dominated case decreases toward the flood-dominated case as the ratio of the maximum depth in the deep channel and the minimum depth on the shoal increases. The different forms of the bottom friction also play a role in the LRV. The structures of the 3D LRV and the depth-integrated LRV are simpler, and the intensity of the LRV is two times smaller when the linear bottom friction is used than those when the quadratic bottom friction is used.

A semi-analytic model of magnetized liner inertial fusion

DOE PAGES

McBride, Ryan D.; Slutz, Stephen A.

2015-05-21

Presented is a semi-analytic model of magnetized liner inertial fusion (MagLIF). This model accounts for several key aspects of MagLIF, including: (1) preheat of the fuel (optionally via laser absorption); (2) pulsed-power-driven liner implosion; (3) liner compressibility with an analytic equation of state, artificial viscosity, internal magnetic pressure, and ohmic heating; (4) adiabatic compression and heating of the fuel; (5) radiative losses and fuel opacity; (6) magnetic flux compression with Nernst thermoelectric losses; (7) magnetized electron and ion thermal conduction losses; (8) end losses; (9) enhanced losses due to prescribed dopant concentrations and contaminant mix; (10) deuterium-deuterium and deuterium-tritium primarymore » fusion reactions for arbitrary deuterium to tritium fuel ratios; and (11) magnetized α-particle fuel heating. We show that this simplified model, with its transparent and accessible physics, can be used to reproduce the general 1D behavior presented throughout the original MagLIF paper [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)]. We also discuss some important physics insights gained as a result of developing this model, such as the dependence of radiative loss rates on the radial fraction of the fuel that is preheated.« less

A semi-analytic model of magnetized liner inertial fusion

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

McBride, Ryan D.; Slutz, Stephen A.

Presented is a semi-analytic model of magnetized liner inertial fusion (MagLIF). This model accounts for several key aspects of MagLIF, including: (1) preheat of the fuel (optionally via laser absorption); (2) pulsed-power-driven liner implosion; (3) liner compressibility with an analytic equation of state, artificial viscosity, internal magnetic pressure, and ohmic heating; (4) adiabatic compression and heating of the fuel; (5) radiative losses and fuel opacity; (6) magnetic flux compression with Nernst thermoelectric losses; (7) magnetized electron and ion thermal conduction losses; (8) end losses; (9) enhanced losses due to prescribed dopant concentrations and contaminant mix; (10) deuterium-deuterium and deuterium-tritium primarymore » fusion reactions for arbitrary deuterium to tritium fuel ratios; and (11) magnetized α-particle fuel heating. We show that this simplified model, with its transparent and accessible physics, can be used to reproduce the general 1D behavior presented throughout the original MagLIF paper [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)]. We also discuss some important physics insights gained as a result of developing this model, such as the dependence of radiative loss rates on the radial fraction of the fuel that is preheated.« less

Weak ductile shear zone beneath the western North Anatolian Fault Zone: inferences from earthquake cycle model constrained by geodetic observations

NASA Astrophysics Data System (ADS)

Yamasaki, T.; Wright, T. J.; Houseman, G. A.

2013-12-01

After large earthquakes, rapid postseismic transient motions are commonly observed. Later in the loading cycle, strain is typically focused in narrow regions around the fault. In simple two-layer models of the loading cycle for strike-slip faults, rapid post-seismic transients require low viscosities beneath the elastic layer, but localized strain later in the cycle implies high viscosities in the crust. To explain this apparent paradox, complex transient rheologies have been invoked. Here we test an alternative hypothesis in which spatial variations in material properties of the crust can explain the geodetic observations. We use a 3D viscoelastic finite element code to examine two simple models of periodic fault slip: a stratified model in which crustal viscosity decreases exponentially with depth below an upper elastic layer, and a block model in which a low viscosity domain centered beneath the fault is embedded in a higher viscosity background representing normal crust. We test these models using GPS data acquired before and after the 1999 Izmit/Duzce earthquakes on the North Anatolian Fault Zone (Turkey). The model with depth-dependent viscosity can show both high postseismic velocities, and preseismic localization of the deformation, if the viscosity contrast from top to bottom of layer exceeds a factor of about 104. However, with no lateral variations in viscosity, this model cannot explain the proximity to the fault of maximum postseismic velocities. In contrast, the model which includes a localized weak zone beneath the faulted elastic lid can explain all the observations, if the weak zone extends down to mid-crustal levels and outward to 10 or 20 km from the fault. The non-dimensional ratio of relaxation time to earthquake repeat time, τ/Δt, is the critical parameter in controlling the observed deformation. In the weak-zone model, τ/Δt should be in the range 0.005 to 0.01 in the weak domain, and larger than ~ 1.0 elsewhere. This implies a viscosity in the weak zone of ~ 1018×0.3 Pa s, and larger than ~ 1020 Pa s outside this region. Models with sharp boundaries to the weak zone fit the data better than those with a smooth increase of viscosity away from the fault. Thus abrupt changes in material properties, such as those that might result from grain-size reduction, may be required in addition to any effect from shear heating. Unlike some previous models, we do not require non-linear stress-dependent viscosities. Our models imply that geodetic strain rates decay to a quasi-steady state within about 10% of the inter-earthquake period (years or decades) and that interseismic geodetic observations can therefore be used to infer the long-term geological slip rate, provided there has not been a recent earthquake. Rheologies inferred from postseismic studies alone likely reflect the rheology of the weak zone beneath the fault, and should not be used to infer the strength profile of normal lithosphere.

The micromechanics model analysis of the viscosity regulation of ultra-high strength concrete with low viscosity

NASA Astrophysics Data System (ADS)

Zhu, M.; Wang, F. G.; Wang, F. Z.; Liu, Y. P.

2017-02-01

The plastic viscosity of mortar and concrete with different binder content, sand ratio, water-binder ratio, microbead dosage and different class and dosage of fly ash were tested and calculated according tomicromechanics model proposed by A. Ghanbari and B.L. Karihaloo, The correlations between these parameters and fresh concrete workability were also investigated, which showed i. high consistence with the objective reality. When binder content, microbead dosage, fly ash dosage or the water-binder ratio was increased or sand ratio was reduced, the fresh concrete viscosity would decrease correspondingly. However their effects were not that same. The relationships between T50 a, V-funnel and inverted slump time with fresh concrete viscosity were established, respectively.

Effects of Artificial Viscosity on the Accuracy of High-reynolds-number Kappa-epsilon Turbulence Model

NASA Technical Reports Server (NTRS)

Chitsomboon, Tawit

1994-01-01

Wall functions, as used in the typical high Reynolds number k-epsilon turbulence model, can be implemented in various ways. A least disruptive method (to the flow solver) is to directly solve for the flow variables at the grid point next to the wall while prescribing the values of k and epsilon. For the centrally-differenced finite-difference scheme employing artificial viscocity (AV) as a stabilizing mechanism, this methodology proved to be totally useless. This is because the AV gives rise to a large error at the wall due to too steep a velocity gradient resulting from the use of a coarse grid as required by the wall function methodology. This error can be eliminated simply by extrapolating velocities at the wall, instead of using the physical values of the no-slip velocities (i.e. the zero value). The applicability of the technique used in this paper is demonstrated by solving a flow over a flat plate and comparing the results with those of experiments. It was also observed that AV gives rise to a velocity overshoot (about 1 percent) near the edge of the boundary layer. This small velocity error, however, can yield as much as 10 percent error in the momentum thickness. A method which integrates the boundary layer up to only the edge of the boundary (instead of infinity) was proposed and demonstrated to give better results than the standard method.

ARRHENIUS MODEL FOR HIGH-TEMPERATURE GLASS VISCOSITY WITH A CONSTANT PRE-EXPONENTIAL FACTOR

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

Hrma, Pavel R.

2008-04-15

A simplified form of the Arrhenius equation, ln η = A + B(x)/T, where η is the viscosity, T the temperature, x the composition vector, and A and B the Arrhenius coefficients, was fitted to glass-viscosity data for the processing temperature range (the range at which the viscosity is within 1 to 103 Pa.s) while setting A = constant and treating B(x) as a linear function of mass fractions of major components. Fitting the Arrhenius equation to over 550 viscosity data of commercial glasses and approximately 1000 viscosity data of glasses for nuclear-waste glasses resulted in the A values ofmore » -11.35 and -11.48, respectively. The R2 value ranged from 0.92 to 0.99 for commercial glasses and was 0.98 for waste glasses. The Arrhenius models estimate viscosities for melts of commercial glasses containing 42 to 84 mass% SiO2 within the temperature range of 1100 to 1550°C and viscosity range of 5 to 400 Pa.s and for waste glasses containing 32 to 60 mass% SiO2 within the temperature range of 850 to 1450°C and viscosity range of 0.4 to 250 Pa.s.« less

Computational fluid dynamics of airfoils and wings

NASA Technical Reports Server (NTRS)

Garabedian, P.; Mcfadden, G.

1982-01-01

It is pointed out that transonic flow is one of the fields where computational fluid dynamics turns out to be most effective. Codes for the design and analysis of supercritical airfoils and wings have become standard tools of the aircraft industry. The present investigation is concerned with mathematical models and theorems which account for some of the progress that has been made. The most successful aerodynamics codes are those for the analysis of flow at off-design conditions where weak shock waves appear. A major breakthrough was achieved by Murman and Cole (1971), who conceived of a retarded difference scheme which incorporates artificial viscosity to capture shocks in the supersonic zone. This concept has been used to develop codes for the analysis of transonic flow past a swept wing. Attention is given to the trailing edge and the boundary layer, entropy inequalities and wave drag, shockless airfoils, and the inverse swept wing code.

Heating, Cooling, and Gravitational Instabilities in Protostellar and Protoplanetary Disks

NASA Astrophysics Data System (ADS)

Pickett, B. K.; Mejia, A. C.; Durisen, R. H.

2001-12-01

We present three-dimensional hydrodynamic simulations of protostellar disk models, in order to explore how the interplay between heating and cooling regulates significant gravitational instabilities. Artificial viscosity is used to treat irreversible heating, such as would occur in shocks; volumetric cooling at several different rates is also applied throughout a broad radial region of the disk. We study the evolution of a disk that is already unstable (due to the low value of the Toomre Q parameter), and a marginally unstable disk that is cooled towards instability. The evolutions have implications for the transport of mass and angular momentum in protostellar disks, the effects of gravitational instabilities on the vertical structure of the disks, and the formation of stellar and substellar companions on dynamic time scales due to disk instabilties. This work is supported by grants from the NASA Planetary Geology and Geophysics and Origins of Solar Systems Programs.

A new method for solving the quantum hydrodynamic equations of motion: application to two-dimensional reactive scattering.

PubMed

Pauler, Denise K; Kendrick, Brian K

2004-01-08

The de Broglie-Bohm hydrodynamic equations of motion are solved using a meshless method based on a moving least squares approach and an arbitrary Lagrangian-Eulerian frame of reference. A regridding algorithm adds and deletes computational points as needed in order to maintain a uniform interparticle spacing, and unitary time evolution is obtained by propagating the wave packet using averaged fields. 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. The methodology is applied to a two-dimensional model collinear reaction with an activation barrier. Reaction probabilities are computed as a function of both time and energy, and are in excellent agreement with those based on the quantum trajectory method. (c) 2004 American Institute of Physics

New Numerical Approaches for Modeling Thermochemical Convection in a Compositionally Stratified Fluid

NASA Astrophysics Data System (ADS)

Puckett, E. G.; Turcotte, D. L.; He, Y.; Lokavarapu, H. V.; Robey, J.; Kellogg, L. H.

2017-12-01

Geochemical observations of mantle-derived rocks favor a nearly homogeneous upper mantle, the source of mid-ocean ridge basalts (MORB), and heterogeneous lower mantle regions.Plumes that generate ocean island basalts are thought to sample the lower mantle regions and exhibit more heterogeneity than MORB.These regions have been associated with lower mantle structures known as large low shear velocity provinces below Africa and the South Pacific.The isolation of these regions is attributed to compositional differences and density stratification that, consequently, have been the subject of computational and laboratory modeling designed to determine the parameter regime in which layering is stable and understanding how layering evolves.Mathematical models of persistent compositional interfaces in the Earth's mantle may be inherently unstable, at least in some regions of the parameter space relevant to the mantle.Computing approximations to solutions of such problems presents severe challenges, even to state-of-the-art numerical methods.Some numerical algorithms for modeling the interface between distinct compositions smear the interface at the boundary between compositions, such as methods that add numerical diffusion or `artificial viscosity' in order to stabilize the algorithm. We present two new algorithms for maintaining high-resolution and sharp computational boundaries in computations of these types of problems: a discontinuous Galerkin method with a bound preserving limiter and a Volume-of-Fluid interface tracking algorithm.We compare these new methods with two approaches widely used for modeling the advection of two distinct thermally driven compositional fields in mantle convection computations: a high-order accurate finite element advection algorithm with entropy viscosity and a particle method.We compare the performance of these four algorithms on three problems, including computing an approximation to the solution of an initially compositionally stratified fluid at Ra = 105 with buoyancy numbers {B} that vary from no stratification at B = 0 to stratified flow at large B.

A priori testing of subgrid-scale models for the velocity-pressure and vorticity-velocity formulations

NASA Technical Reports Server (NTRS)

Winckelmans, G. S.; Lund, T. S.; Carati, D.; Wray, A. A.

1996-01-01

Subgrid-scale models for Large Eddy Simulation (LES) in both the velocity-pressure and the vorticity-velocity formulations were evaluated and compared in a priori tests using spectral Direct Numerical Simulation (DNS) databases of isotropic turbulence: 128(exp 3) DNS of forced turbulence (Re(sub(lambda))=95.8) filtered, using the sharp cutoff filter, to both 32(exp 3) and 16(exp 3) synthetic LES fields; 512(exp 3) DNS of decaying turbulence (Re(sub(Lambda))=63.5) filtered to both 64(exp 3) and 32(exp 3) LES fields. Gaussian and top-hat filters were also used with the 128(exp 3) database. Different LES models were evaluated for each formulation: eddy-viscosity models, hyper eddy-viscosity models, mixed models, and scale-similarity models. Correlations between exact versus modeled subgrid-scale quantities were measured at three levels: tensor (traceless), vector (solenoidal 'force'), and scalar (dissipation) levels, and for both cases of uniform and variable coefficient(s). Different choices for the 1/T scaling appearing in the eddy-viscosity were also evaluated. It was found that the models for the vorticity-velocity formulation produce higher correlations with the filtered DNS data than their counterpart in the velocity-pressure formulation. It was also found that the hyper eddy-viscosity model performs better than the eddy viscosity model, in both formulations.

A Non-Arrhenian Viscosity Model for Natural Silicate Melts with Applications to Volcanology

NASA Astrophysics Data System (ADS)

Russell, J. K.; Giordano, D.; Dingwell, D. B.

2005-12-01

Silicate melt viscosity is the most important physical property in volcanic systems. It governs styles and rates of flow, velocity distributions in flowing magma, rates of vesiculation, and, ultimately, sets limits on coherent(vs. fragmented or disrupted) flow. The prediction of melt viscosity over the range of conditions found on terrestrial planets remains a challenge. However, the extraordinary increase in number and quality of published measurements of melt viscosity suggests the possibility of new models. Here we review the attributes of previous models for silicate melt viscosity and, then, present a new predictive model natural silicate melts. The importance of silicate melt viscosity was recognized early [1] and culminated in 2 models for predicting silicate melt viscosity [2,3]. These models used an Arrhenian T-dependence; they were limited by a limited experimental database dominated by high-T measurements. Subsequent models have aimed to: i) extend the compositional range of Arrhenian T-dependent models [4,5]; ii) to develop non-Arrhenian models for limited ranges of composition [6,7,8], iii) to develop new strategies for modelling the composition and T-dependence of viscosity [9,10,11], and, finally, to create chemical models for the non-Arrhenian T-dependence of natural melts [12]. We present a multicomponent model for the compositional and T dependence of silicate melt viscosity based on data spanning a wide range of anhydrous melt compositions. The experimental data include micropenetration and concentric cylinder viscometry measurements covering a viscosity range of 10-1 to 1012 Pa s and a T-range from 700 to 1650°C. These published data provide a high- quality database comprising ~ 800 experimental data on 44 well-characterized melt compositions. Our model uses the Adam-Gibbs equation to capture T-dependence: log η = A + B/[T · log (T/C)] where A, B, and C are adjustable parameters that vary for different melt compositions. We assume that all silicate melts converge to a common, but unknown, high-T limit (e.g., A) and that all compositional dependence is accommodated for by B and C. We adopt a linear compositional dependence for B and C: B = σi=1..n [xi βi] C = σi=1..n [xi γi] where xi's are the mole fractions of oxide components (n=8) and βi and γi are adjustable parameters. The model, therefore, comprises 2 · n+1 adjustable parameters which are optimized for against the experimental database including a common value of A and compositional coefficeints for B and C. The new model reproduces the original database to within experimental uncertainty and can predict the viscosity of silicate melts across the full range of conditions found in Nature. References Cited: [1] Friedman et al., 1963. J Geophys Res 68, 6523-6535. [2] Bottinga Y & Weill D 1972. Am J Sci 272, 438- 475. [3] Shaw HR 1972. Am J Sci 272, 438- 475. [4] Persikov ES 1991. Adv Phys Geochem 9, 1-40. [5] Prusevich AA 1988. Geol Geofiz 29, 67-69. [6] Baker DR 1996. Am Min 81, 126-134. [7] Hess KU & Dingwell DB 1996. Am Min 81, 1297- 1300. [8] Zhang, et al. 2003. Am min 88, 1741- 1752. [9] Russell et al. 2002. Eur J Min 14, 417-428. [10] Russell et al. 2003. Am Min 8, 1390- 1394. [11] Russell JK & Giordano D In Press. Geochim Cosmochim Acta. [12] Giordano D & Dingwell DB 2003. Earth Planet. Sci. Lett. 208, 337-349.

Reynolds Stress Distributions and the Measurement and Calculation of Eddy Viscosity in Gravity Currents

NASA Astrophysics Data System (ADS)

Kelly, R. W.; Chalk, C.; Dorrell, R. M.; Peakall, J.; Burns, A. D.; Keevil, G. M.; Thomas, R. E.; Williams, G.

2016-12-01

In the natural environment, gravity currents transport large volumes of sediment great distances and are often considered one of the most important mechanisms for sediment transport in ocean basins. Deposits from many individual submarine gravity currents, turbidites, ultimately form submarine fan systems. These are the largest sedimentary systems on the planet and contain valuable hydrocarbon reserves. Moreover, the impact of these currents on submarine technologies and seafloor infrastructure can be devastating and therefore they are of significant interest to a wide range of industries. Here we present experimental, numerical and theoretical models of time-averaged turbulent shear stresses, i.e. Reynolds stresses. Reynolds stresses can be conceptually parameterised by an eddy viscosity parameter that relates chaotic fluid motion to diffusive type processes. As such, it is a useful parameter for indicating the extent of internal mixing and is used extensively in both numerical and analytical modelling of both open-channel and gravity driven flows. However, a lack of knowledge of the turbulent structure of gravity currents limits many hydro- and morphodynamic models. High resolution 3-dimensional experimental velocity data, gathered using acoustic Doppler profiling velocimetry, enabled direct calculation of stresses and eddy viscosity. Comparison of experimental data to CFD and analytical models allowed the testing of eddy viscosity-based turbulent mixing models. The calculated eddy viscosity profile is parabolic in nature in both the upper and lower shear layers. However, an apparent breakdown in the Boussinesq hypothesis (used to calculate the eddy viscosity and upon which many numerical models are based) is observed in the region of the current around the velocity maximum. With the use of accompanying density data it is suggested that the effect of stratification on eddy viscosity is significant and alternative formulations may be required.

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.

Viscosity and transient electric birefringence study of clay colloidal aggregation.

PubMed

Bakk, Audun; Fossum, Jon O; da Silva, Geraldo J; Adland, Hans M; Mikkelsen, Arne; Elgsaeter, Arnljot

2002-02-01

We study a synthetic clay suspension of laponite at different particle and NaCl concentrations by measuring stationary shear viscosity and transient electrically induced birefringence (TEB). On one hand the viscosity data are consistent with the particles being spheres and the particles being associated with large amount bound water. On the other hand the viscosity data are also consistent with the particles being asymmetric, consistent with single laponite platelets associated with a very few monolayers of water. We analyze the TEB data by employing two different models of aggregate size (effective hydrodynamic radius) distribution: (1) bidisperse model and (2) log-normal distributed model. Both models fit, in the same manner, fairly well to the experimental TEB data and they indicate that the suspension consists of polydisperse particles. The models also appear to confirm that the aggregates increase in size vs increasing ionic strength. The smallest particles at low salt concentrations seem to be monomers and oligomers.

Internally heated mantle convection and the thermal and degassing history of the earth

NASA Technical Reports Server (NTRS)

Williams, David R.; Pan, Vivian

1992-01-01

An internally heated model of parameterized whole mantle convection with viscosity dependent on temperature and volatile content is examined. The model is run for 4l6 Gyr, and temperature, heat flow, degassing and regassing rates, stress, and viscosity are calculated. A nominal case is established which shows good agreement with accepted mantle values. The effects of changing various parameters are also tested. All cases show rapid cooling early in the planet's history and strong self-regulation of viscosity due to the temperature and volatile-content dependence. The effects of weakly stress-dependent viscosity are examined within the bounds of this model and are found to be small. Mantle water is typically outgassed rapidly to reach an equilibrium concentration on a time scale of less than 200 Myr for almost all models, the main exception being for models which start out with temperatures well below the melting temperature.

A singularity free approach to post glacial rebound calculations

NASA Technical Reports Server (NTRS)

Fang, Ming; Hager, Bradford H.

1994-01-01

Calculating the post glacial response of a viscoelastic Earth model using the exponential decay normal mode technique leads to intrinsic singularities if viscosity varies continuously as a function of radius. We develop a numerical implementation of the Complex Real Fourier transform (CRFT) method as an accurate and stable procedure to avoid these singularities. Using CRFT, we investigate the response of a set of Maxwell Earth models to surface loading. We find that the effect of expanding a layered viscosity structure into a continuously varying structure is to destroy the modes associated with the boundary between layers. Horizontal motion is more sensitive than vertical motion to the viscosity structure just below the lithosphere. Horizontal motion is less sensitive to the viscosity of the lower mantle than the vertical motion is. When the viscosity increases at 670 km depth by a factor of about 60, the response of the lower mantle is close to its elastic limit. Any further increase of the viscosity contrast at 670 km depth or further increase of viscosity as a continuous function of depth starting from 670 km depth is unlikely to be resolved.

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.

Gibbs Energy Additivity Approaches in Estimation of Dynamic Viscosities of n-Alkane-1-ol

NASA Astrophysics Data System (ADS)

Phankosol, S.; Krisnangkura, K.

2017-09-01

Alcohols are solvents for organic and inorganic substances. Dynamic viscosity of liquid is important transport properties. In this study models for estimating n-alkan-1-ol dynamic viscosities are correlated to the Martin’s rule of free energy additivity. Data available in literatures are used to validate and support the proposed equations. The dynamic viscosities of n-alkan-1-ol can be easily estimated from its carbon numbers (nc) and temperatures (T). The bias, average absolute deviation and coefficient of determination (R2) in estimating of n-alkan-1-ol are -0.17%, 1.73% and 0.999, respectively. The dynamic viscosities outside temperature between 288.15 and 363.15 K may be possibly estimated by this model but accuracy may be lower.

Two-phase model for prediction of cell-free layer width in blood flow

PubMed Central

Namgung, Bumseok; Ju, Meongkeun; Cabrales, Pedro; Kim, Sangho

2014-01-01

This study aimed to develop a numerical model capable of predicting changes in the cell-free layer (CFL) width in narrow tubes with consideration of red blood cell aggregation effects. The model development integrates to empirical relations for relative viscosity (ratio of apparent viscosity to medium viscosity) and core viscosity measured on independent blood samples to create a continuum model that includes these two regions. The constitutive relations were derived from in vitro experiments performed with three different glass-capillary tubes (inner diameter = 30, 50 and 100 μm) over a wide range of pseudoshear rates (5-300 s−1). The aggregation tendency of the blood samples was also varied by adding Dextran 500 kDa. Our model predicted that the CFL width was strongly modulated by the relative viscosity function. Aggregation increased the width of CFL, and this effect became more pronounced at low shear rates. The CFL widths predicted in the present study at high shear conditions were in agreement with those reported in previous studies. However, unlike previous multi-particle models, our model did not require a high computing cost, and it was capable of reproducing results for a thicker CFL width at low shear conditions, depending on aggregating tendency of the blood. PMID:23116701

Chemoviscosity modeling for thermosetting resins

NASA Technical Reports Server (NTRS)

Tiwari, S. N.; Hou, T. H.; Bai, J. M.

1985-01-01

A chemoviscosity model, which describes viscosity rise profiles accurately under various cure cycles, and correlates viscosity data to the changes of physical properties associated with structural transformations of the thermosetting resin system during cure, was established. Work completed on chemoviscosity modeling for thermosetting resins is reported.

Harnessing Electrostatic Forces to Grow Bio-inspired Hierarchical Vascular Networks

NASA Astrophysics Data System (ADS)

Behler, Kristopher; Melrose, Zachary; Schott, Andrew; Wetzel, Eric

2012-02-01

Vascular networks provide a system for fluid distribution. Artificial vascular materials with enhanced properties are currently being developed that could ultimately be integrated into systems reliant upon fluid transport while retaining their structural properties. An uninterrupted and controllable supply of liquid is optimal for many applications such as continual self-healing materials, in-situ delivery of index matched fluids, thermal management and drug delivery systems could benefit from a bio-inspired vascular approach that combines complex network geometries with minimal processing parameters. Two such approaches to induce vascular networks are electrohydrodynamic viscous fingering (EHVF) and electrical treeing (ET). EHVF is a phenomenon that occurs when a low viscosity liquid is forced through a high viscosity fluid or matrix, resulting in branches due to capillary and viscous forces in the high viscosity material. By applying voltages of 0 -- 60 kV, finger diameter is reduced. ET is the result of partial discharges in a dielectric material. In the vicinity of a small diameter electrode, the local electric field is greater than the global dielectric strength, causing a localized, step-wise, breakdown to occur forming a highly branched interconnected structure. ET is a viable method to produce networks on a smaller, micron, scale than the products of the EHVF method.

Seismic Constraints on the Mantle Viscosity Structure beneath Antarctica

NASA Astrophysics Data System (ADS)

Wiens, Douglas; Heeszel, David; Aster, Richard; Nyblade, Andrew; Wilson, Terry

2015-04-01

Lateral variations in upper mantle viscosity structure can have first order effects on glacial isostatic adjustment. These variations are expected to be particularly large for the Antarctic continent because of the stark geological contrast between ancient cratonic and recent tectonically active terrains in East and West Antarctica, respectively. A large misfit between observed and predicted GPS rates for West Antarctica probably results in part from the use of a laterally uniform viscosity structure. Although not linked by a simple relationship, mantle seismic velocities can provide important constraints on mantle viscosity structure, as they are both largely controlled by temperature and water content. Recent higher resolution seismic models for the Antarctic mantle, derived from data acquired by new seismic stations deployed in the AGAP/GAMSEIS and ANET/POLENET projects, offer the opportunity to use the seismic velocity structure to place new constraints on the viscosity of the Antarctic upper mantle. We use an Antarctic shear wave velocity model derived from array analysis of Rayleigh wave phase velocities [Heeszel et al, in prep] and examine a variety of methodologies for relating seismic, thermal and rheological parameters to compute a suite of viscosity models for the Antarctic mantle. A wide variety of viscosity structures can be derived using various assumptions, but they share several robust common elements. There is a viscosity contrast of at least two orders of magnitude between East and West Antarctica at depths of 80-250 km, reflecting the boundary between cold cratonic lithosphere in East Antarctica and warm upper mantle in West Antarctica. The region beneath the Ellsworth-Whitmore Mtns and extending to the Pensacola Mtns. shows intermediate viscosity between the extremes of East and West Antarctica. There are also significant variations between different parts of West Antarctica, with the lowest viscosity occurring beneath the Marie Byrd Land (MBL). The MBL Dome and adjacent coastal areas show extremely low viscosity (~1018Pa-s) for most parameterizations, suggesting that low mantle viscosity may produce a very rapid response to ice mass loss in this region.

VNAP2: A Computer Program for Computation of Two-dimensional, Time-dependent, Compressible, Turbulent Flow

NASA Technical Reports Server (NTRS)

Cline, M. C.

1981-01-01

A computer program, VNAP2, for calculating turbulent (as well as laminar and inviscid), steady, and unsteady flow is presented. It solves the two dimensional, time dependent, compressible Navier-Stokes equations. The turbulence is modeled with either an algebraic mixing length model, a one equation model, or the Jones-Launder two equation model. The geometry may be a single or a dual flowing stream. The interior grid points are computed using the unsplit MacCormack scheme. Two options to speed up the calculations for high Reynolds number flows are included. The boundary grid points are computed using a reference plane characteristic scheme with the viscous terms treated as source functions. An explicit artificial viscosity is included for shock computations. The fluid is assumed to be a perfect gas. The flow boundaries may be arbitrary curved solid walls, inflow/outflow boundaries, or free jet envelopes. Typical problems that can be solved concern nozzles, inlets, jet powered afterbodies, airfoils, and free jet expansions. The accuracy and efficiency of the program are shown by calculations of several inviscid and turbulent flows. The program and its use are described completely, and six sample cases and a code listing are included.

Advancing parabolic operators in thermodynamic MHD models: Explicit super time-stepping versus implicit schemes with Krylov solvers

NASA Astrophysics Data System (ADS)

Caplan, R. M.; Mikić, Z.; Linker, J. A.; Lionello, R.

2017-05-01

We explore the performance and advantages/disadvantages of using unconditionally stable explicit super time-stepping (STS) algorithms versus implicit schemes with Krylov solvers for integrating parabolic operators in thermodynamic MHD models of the solar corona. Specifically, we compare the second-order Runge-Kutta Legendre (RKL2) STS method with the implicit backward Euler scheme computed using the preconditioned conjugate gradient (PCG) solver with both a point-Jacobi and a non-overlapping domain decomposition ILU0 preconditioner. The algorithms are used to integrate anisotropic Spitzer thermal conduction and artificial kinematic viscosity at time-steps much larger than classic explicit stability criteria allow. A key component of the comparison is the use of an established MHD model (MAS) to compute a real-world simulation on a large HPC cluster. Special attention is placed on the parallel scaling of the algorithms. It is shown that, for a specific problem and model, the RKL2 method is comparable or surpasses the implicit method with PCG solvers in performance and scaling, but suffers from some accuracy limitations. These limitations, and the applicability of RKL methods are briefly discussed.

Eulerian adaptive finite-difference method for high-velocity impact and penetration problems

NASA Astrophysics Data System (ADS)

Barton, P. T.; Deiterding, R.; Meiron, D.; Pullin, D.

2013-05-01

Owing to the complex processes involved, faithful prediction of high-velocity impact events demands a simulation method delivering efficient calculations based on comprehensively formulated constitutive models. Such an approach is presented herein, employing a weighted essentially non-oscillatory (WENO) method within an adaptive mesh refinement (AMR) framework for the numerical solution of hyperbolic partial differential equations. Applied widely in computational fluid dynamics, these methods are well suited to the involved locally non-smooth finite deformations, circumventing any requirement for artificial viscosity functions for shock capturing. Application of the methods is facilitated through using a model of solid dynamics based upon hyper-elastic theory comprising kinematic evolution equations for the elastic distortion tensor. The model for finite inelastic deformations is phenomenologically equivalent to Maxwell's model of tangential stress relaxation. Closure relations tailored to the expected high-pressure states are proposed and calibrated for the materials of interest. Sharp interface resolution is achieved by employing level-set functions to track boundary motion, along with a ghost material method to capture the necessary internal boundary conditions for material interactions and stress-free surfaces. The approach is demonstrated for the simulation of high velocity impacts of steel projectiles on aluminium target plates in two and three dimensions.

The effect of hydroxyl functional groups and molar mass on the viscosity of non-crystalline organic and organic-water particles

NASA Astrophysics Data System (ADS)

Grayson, James W.; Evoy, Erin; Song, Mijung; Chu, Yangxi; Maclean, Adrian; Nguyen, Allena; Upshur, Mary Alice; Ebrahimi, Marzieh; Chan, Chak K.; Geiger, Franz M.; Thomson, Regan J.; Bertram, Allan K.

2017-07-01

The viscosities of three polyols and three saccharides, all in the non-crystalline state, have been studied. Two of the polyols (2-methyl-1,4-butanediol and 1,2,3-butanetriol) were studied under dry conditions, the third (1,2,3,4-butanetetrol) was studied as a function of relative humidity (RH), including under dry conditions, and the saccharides (glucose, raffinose, and maltohexaose) were studied as a function of RH. The mean viscosities of the polyols under dry conditions range from 1.5 × 10-1 to 3.7 × 101 Pa s, with the highest viscosity being that of the tetrol. Using a combination of data determined experimentally here and literature data for alkanes, alcohols, and polyols with a C3 to C6 carbon backbone, we show (1) there is a near-linear relationship between log10 (viscosity) and the number of hydroxyl groups in the molecule, (2) that on average the addition of one OH group increases the viscosity by a factor of approximately 22 to 45, (3) the sensitivity of viscosity to the addition of one OH group is not a strong function of the number of OH functional groups already present in the molecule up to three OH groups, and (4) higher sensitivities are observed when the molecule has more than three OH groups. Viscosities reported here for 1,2,3,4-butanetetrol particles are lower than previously reported measurements using aerosol optical tweezers, and additional studies are required to resolve these discrepancies. For saccharide particles at 30 % RH, viscosity increases by approximately 2-5 orders of magnitude as molar mass increases from 180 to 342 g mol-1, and at 80 % RH, viscosity increases by approximately 4-5 orders of magnitude as molar mass increases from 180 to 991 g mol-1. These results suggest oligomerization of highly oxidized compounds in atmospheric secondary organic aerosol (SOA) could lead to large increases in viscosity, and may be at least partially responsible for the high viscosities observed in some SOA. Finally, two quantitative structure-property relationship models (Sastri and Rao, 1992; Marrero-Morejón and Pardillo-Fontdevila, 2000) were used to predict the viscosity of alkanes, alcohols, and polyols with a C3-C6 carbon backbone. Both models show reasonably good agreement with measured viscosities for the alkanes, alcohols, and polyols studied here except for the case of a hexol, the viscosity of which is underpredicted by 1-3 orders of magnitude by each of the models.

High-Resolution Lithosphere Viscosity and Dynamics Revealed by Magnetotelluric Imaging

NASA Astrophysics Data System (ADS)

Liu, L.; Hasterok, D. P.

2016-12-01

An accurate viscosity structure is critical to truthfully modeling continental lithosphere dynamics, especially at spatial scales of <200 km where active tectonic deformation and volcanism occur. However, the effective viscosity structure of the lithosphere remains a key challenge in geodynamics due to the intimate involvement of viscosity with time and its dependence on many factors including strain rate, plastic failure, composition, and grain size. Current efforts on inferring the detailed lithosphere viscosity structure are sparse and large uncertainties and discrepancies still exist. Here we report an attempt to infer the effective lithospheric viscosity from a high-resolution magnetotelluric (MT) survey across the western United States. The high sensitivity of MT fields to the presence of electrically conductive fluids makes it a promising proxy for determining mechanical strength variations throughout the lithosphere. We demonstrate how a viscosity structure, approximated from electrical resistivity, results in a geodynamic model that successfully predicts short-wavelength surface topography, lithospheric deformation, and mantle upwelling beneath recent volcanism. The results indicate that lithosphere viscosity structure rather than the buoyancy structure is the dominant controlling factor for short-wavelength topography and intra-plate deformation in tectonically active regions. We further show that this viscosity is consistent with and more effective than that derived from laboratory-based rheology. We therefore propose that MT imaging provides a practical observational constraint for quantifying the dynamic evolution of the continental lithosphere.

The role of viscosity in TATB hot spot ignition

NASA Astrophysics Data System (ADS)

Fried, Laurence E.; Zepeda-Ruis, Luis; Howard, W. Michael; Najjar, Fady; Reaugh, John E.

2012-03-01

The role of dissipative effects, such as viscosity, in the ignition of high explosive pores is investigated using a coupled chemical, thermal, and hydrodynamic model. Chemical reactions are tracked with the Cheetah thermochemical code coupled to the ALE3D hydrodynamic code. We perform molecular dynamics simulations to determine the viscosity of liquid TATB. We also analyze shock wave experiments to obtain an estimate for the shock viscosity of TATB. Using the lower bound liquid-like viscosities, we find that the pore collapse is hydrodynamic in nature. Using the upper bound viscosity from shock wave experiments, we find that the pore collapse is closest to the viscous limit.

Viscosity as related to dietary fiber: a review.

PubMed

Dikeman, Cheryl L; Fahey, George C

2006-01-01

Viscosity is a physicochemical property associated with dietary fibers, particularly soluble dietary fibers. Viscous dietary fibers thicken when mixed with fluids and include polysaccharides such as gums, pectins, psyllium, and beta-glucans. Although insoluble fiber particles may affect viscosity measurement, viscosity is not an issue regards insoluble dietary fibers. Viscous fibers have been credited for beneficial physiological responses in human, animal, and animal-alternative in vitro models. The following article provides a review of viscosity as related to dietary fiber including definitions and instrumentation, factors affecting viscosity of solutions, and effects of viscous polysaccharides on glycemic response, blood lipid attenuation, intestinal enzymatic activity, digestibility, and laxation.

Comment on "An alternative theory to explain the effects of coalescing oil drops on mouthfeel" by B. Le Reverend and J. Engmann, Soft Matter, 2015, 11, 7077.

PubMed

Xia, Qiuyang

2016-03-28

In a recent paper by B. Le Reverend and J. Engmann, they used a model to explain the change in the perceived viscosity by phase separation. We improved this model by adding the drop in viscosity in the aqueous phase to it and we show how this will significantly change the conclusion in the original paper. The increase in viscosity due to phase separation is highly unlikely to happen because the drop in viscosity due to loss of oil is faster at a high oil concentration.

Viscosity and Solvation

ERIC Educational Resources Information Center

Robertson, C. T.

1973-01-01

Discusses theories underlying the phenomena of solution viscosities, involving the Jones and Dole equation, B-coefficient determination, and flickering cluster model. Indicates that viscosity measurements provide a basis for the study of the structural effects of ions in aqueous solutions and are applicable in teaching high school chemistry. (CC)

Effect of viscosity on droplet-droplet collisional interaction

NASA Astrophysics Data System (ADS)

Finotello, Giulia; Padding, Johan T.; Deen, Niels G.; Jongsma, Alfred; Innings, Fredrik; Kuipers, J. A. M.

2017-06-01

A complete knowledge of the effect of droplet viscosity on droplet-droplet collision outcomes is essential for industrial processes such as spray drying. When droplets with dispersed solids are dried, the apparent viscosity of the dispersed phase increases by many orders of magnitude, which drastically changes the outcome of a droplet-droplet collision. However, the effect of viscosity on the droplet collision regime boundaries demarcating coalescence and reflexive and stretching separation is still not entirely understood and a general model for collision outcome boundaries is not available. In this work, the effect of viscosity on the droplet-droplet collision outcome is studied using direct numerical simulations employing the volume of fluid method. The role of viscous energy dissipation is analysed in collisions of droplets with different sizes and different physical properties. From the simulations results, a general phenomenological model depending on the capillary number (Ca, accounting for viscosity), the impact parameter (B), the Weber number (We), and the size ratio (Δ) is proposed.

Understanding the importance of the temperature dependence of viscosity on the crystallization dynamics in the Ge2Sb2Te5 phase-change material

NASA Astrophysics Data System (ADS)

Aladool, A.; Aziz, M. M.; Wright, C. D.

2017-06-01

The crystallization dynamics in the phase-change material Ge2Sb2Te5 is modelled using the more detailed Master equation method over a wide range of heating rates commensurate with published ultrafast calorimetry experiments. Through the attachment and detachment of monomers, the Master rate equation naturally traces nucleation and growth of crystallites with temperature history to calculate the transient distribution of cluster sizes in the material. Both the attachment and detachment rates in this theory are strong functions of viscosity, and thus, the value of viscosity and its dependence on temperature significantly affect the crystallization process. In this paper, we use the physically realistic Mauro-Yue-Ellison-Gupta-Allan viscosity model in the Master equation approach to study the role of the viscosity model parameters on the crystallization dynamics in Ge2Sb2Te5 under ramped annealing conditions with heating rates up to 4 × 104 K/s. Furthermore, due to the relatively low computational cost of the Master equation method compared to atomistic level computations, an iterative numerical approach was developed to fit theoretical Kissinger plots simulated with the Master equation system to experimental Kissinger plots from ultrafast calorimetry measurements at increasing heating rates. This provided a more rigorous method (incorporating both nucleation and growth processes) to extract the viscosity model parameters from the analysis of experimental data. The simulations and analysis revealed the strong coupling between the glass transition temperature and fragility index in the viscosity and crystallization models and highlighted the role of the dependence of the glass transition temperature on the heating rate for the accurate estimation of the fragility index of phase-change materials from the analysis of experimental measurements.

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.