Application of the log-conformation tensor to three-dimensional time-dependent free surface flows

The numerical simulation of flows of highly elastic fluids has been the subject of intense research over the past decades with important industrial applications. Therefore, many efforts have been made to improve the convergence capabilities of the numerical methods employed to simulate viscoelastic fluid flows. An important contribution for the solution of the High-Weissenberg Number Problem has been presented by Fattal and Kupferman [J. Non-Newton. Fluid. Mech. 123 (2004) 281-285] who developed the matrix-logarithm of the conformation tensor technique, henceforth called log-conformation tensor. Its advantage is a better approximation of the large growth of the stress tensor that occur in some regions of the flow and it is doubly beneficial in that it ensures physically correct stress fields, allowing converged computations at high Weissenberg number flows. In this work we investigate the application of the log-conformation tensor to three-dimensional unsteady free surface flows. The log-conformation tensor formulation was applied to solve the Upper-Convected Maxwell (UCM) constitutive equation while the momentum equation was solved using a finite difference Marker-and-Cell type method. The resulting developed code is validated by comparing the log-conformation results with the analytic solution for fully developed pipe flows. To illustrate the stability of the log-conformation tensor approach in solving three-dimensional free surface flows, results from the simulation of the extrudate swell and jet buckling phenomena of UCM fluids at high Weissenberg numbers are presented. (C) 2012 Elsevier B.V. All rights reserved.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids. (C) 2010 Elsevier B.V. All rights reserved.

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.

Water regulation of actinomycin-D binding to DNA: the interplay among drug affinity, DNA long-range conformation, and hydration

Actiaomycin-D (actD) binds to natural DNA at two different classes of binding sites, weak and strong. The affinity for these sites is highly dependent on DNA se(sequence and solution conditions, and the interaction appears to be purely entropic driven Although the entropic character of this reaction has been attributed to the release of water molecules upon drug to DNA complex formation, the mechanism by which hydration regulates actD binding and discrimination between different classes of binding sites on natural DNA is still unknown. In this work, we investigate the role of hydration on this reaction using the osmotic stress method. We skew that the decrease of solution water activity, due to the addition of sucrose, glycerol ethylene glycol, and betaine, favors drug binding to the strong binding sites on DNA by increasing both the apparent binding affinity Delta G, and the number of DNA base pairs apparently occupied by the bound drug n(bp/actD). These binding parameters vary linearly with the logarithm of the molar fraction of water in solution log(X-w), which indicates the contribution of water binding to the energetic of the reaction. It is demonstrated that the hydration change measured upon binding increases proportionally to the apparent size of the binding site n(bp/uctD). This indicates that n(bp/actD) measured from the Scatchard plod is a measure of the size of the DNA molecule changing conformation due to ligand binding. We also find that the contribution of DNA deformation, gauged by n(bp/act) to the total free energy of binding Delta G, is given by Delta G = Delta G(local) + n(bp/actD) x delta G(DNA), where Delta G(local), = -8020 +/- 51 cal/mol of actD bound and delta G(DNa) = -24.1 +/- 1.7cal/mol of base pair at 25 degrees C. We interpret Delta G(local), as the energetic contribution due to the direct interactions of actD with the actual tetranucleotide binding site, and it n(bp/actB) X delta G(DNA) as that due to change inconformation, induced by binding, of it n(bp/actD) DNA base pairs flanking the local site. This interpretation is supported by the agreement found between the value of delta G(DNA) and the torsional free energy change measured independently. We conclude suggesting an allosteric model for ligand binding to DNA, such that the increase in binding affinity is achieved by increasing the relaxation of the unfavorable free energy of binding storage at the local site through a larger number of DNA base pairs. The new aspect on this model is that the size of the complex is not fixed but determined by solutions conditions, such as water activity, which modulate the energetic barrier to change helix conformation. These results may suggest that long-range allosteric transitions of duplex DNA are involved in the inhibition of RNA synthesis by actD, and more generally, in the regulation of transcription. (C) 2000 John Wiley & Sons, Inc.

Effects Of High Pressure Homogenization On The Activity, Stability, Kinetics And Three-dimensional Conformation Of A Glucose Oxidase Produced By Aspergillus Niger.

High pressure homogenization (HPH) is a non-thermal method, which has been employed to change the activity and stability of biotechnologically relevant enzymes. This work investigated how HPH affects the structural and functional characteristics of a glucose oxidase (GO) from Aspergillus niger. The enzyme was homogenized at 75 and 150 MPa and the effects were evaluated with respect to the enzyme activity, stability, kinetic parameters and molecular structure. The enzyme showed a pH-dependent response to the HPH treatment, with reduction or maintenance of activity at pH 4.5-6.0 and a remarkable activity increase (30-300%) at pH 6.5 in all tested temperatures (15, 50 and 75°C). The enzyme thermal tolerance was reduced due to HPH treatment and the storage for 24 h at high temperatures (50 and 75°C) also caused a reduction of activity. Interestingly, at lower temperatures (15°C) the activity levels were slightly higher than that observed for native enzyme or at least maintained. These effects of HPH treatment on function and stability of GO were further investigated by spectroscopic methods. Both fluorescence and circular dichroism revealed conformational changes in the molecular structure of the enzyme that might be associated with the distinct functional and stability behavior of GO.

Aspects of a noncommutative scalar/tensor duality

We study the noncommutative massless Kalb-Ramond gauge field coupled to a dynamical U(1) gauge field in the adjoint representation together with a compensating vector field. We derive the Seiberg-Witten map and obtain the corresponding mapped action to first order in theta. The (emergent) gravity structure found in other situations is not present here. The off-shell dual scalar theory is derived and it does not coincide with the Seiberg-Witten mapped scalar theory. Dispersion relations are also discussed. The p-form generalization of the Seiberg-Witten map to order theta is also derived.

One-loop energy-momentum tensor in QED with electriclike background

We have obtained nonperturbative one-loop expressions for the mean-energy-momentum tensor and current density of Dirac's field on a constant electriclike back-round. One of the goals of this calculation is to give a consistent description of backreaction in such a theory. Two cases of initial states are considered: the vacuum state and the thermal equilibrium state. First, we perform calculations for the vacuum initial state. In the obtained expressions, we separate the contributions due to particle creation and vacuum polarization. The latter contribution,, are related to the Heisenberg-Euler Lagrangian. Then, we Study the case of the thermal initial state. Here, we separate the contributions due to particle creation, vacuum polarization, and the contributions due to the work of the external field on the particles at the initial state. All these contributions are studied in detail, in different regimes of weak and strong fields and low and high temperatures. The obtained results allow us to establish restrictions on the electric field and its duration under which QED with a strong constant electric field is consistent. Under such restrictions, one can neglect the backreaction of particles created by the electric field. Some of the obtained results generalize the calculations of Heisenberg-Euler for energy density to the case of arbitrary strong electric fields.

INFERENCE FOR EIGENVALUES AND EIGENVECTORS OF GAUSSIAN SYMMETRIC MATRICES

This article presents maximum likelihood estimators (MLEs) and log-likelihood ratio (LLR) tests for the eigenvalues and eigenvectors of Gaussian random symmetric matrices of arbitrary dimension, where the observations are independent repeated samples from one or two populations. These inference problems are relevant in the analysis of diffusion tensor imaging data and polarized cosmic background radiation data, where the observations are, respectively, 3 x 3 and 2 x 2 symmetric positive definite matrices. The parameter sets involved in the inference problems for eigenvalues and eigenvectors are subsets of Euclidean space that are either affine subspaces, embedded submanifolds that are invariant under orthogonal transformations or polyhedral convex cones. We show that for a class of sets that includes the ones considered in this paper, the MLEs of the mean parameter do not depend on the covariance parameters if and only if the covariance structure is orthogonally invariant. Closed-form expressions for the MLEs and the associated LLRs are derived for this covariance structure.

Residuals for log-Burr XII regression models in survival analysis

In this paper, we compare three residuals to assess departures from the error assumptions as well as to detect outlying observations in log-Burr XII regression models with censored observations. These residuals can also be used for the log-logistic regression model, which is a special case of the log-Burr XII regression model. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and the empirical distribution of each residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended to the modified martingale-type residual in log-Burr XII regression models with censored data.

A log-extended Weibull regression model

A bathtub-shaped failure rate function is very useful in survival analysis and reliability studies. The well-known lifetime distributions do not have this property. For the first time, we propose a location-scale regression model based on the logarithm of an extended Weibull distribution which has the ability to deal with bathtub-shaped failure rate functions. We use the method of maximum likelihood to estimate the model parameters and some inferential procedures are presented. We reanalyze a real data set under the new model and the log-modified Weibull regression model. We perform a model check based on martingale-type residuals and generated envelopes and the statistics AIC and BIC to select appropriate models. (C) 2009 Elsevier B.V. All rights reserved.

The generalized log-gamma mixture model with covariates: local influence and residual analysis

In a sample of censored survival times, the presence of an immune proportion of individuals who are not subject to death, failure or relapse, may be indicated by a relatively high number of individuals with large censored survival times. In this paper the generalized log-gamma model is modified for the possibility that long-term survivors may be present in the data. The model attempts to separately estimate the effects of covariates on the surviving fraction, that is, the proportion of the population for which the event never occurs. The logistic function is used for the regression model of the surviving fraction. Inference for the model parameters is considered via maximum likelihood. Some influence methods, such as the local influence and total local influence of an individual are derived, analyzed and discussed. Finally, a data set from the medical area is analyzed under the log-gamma generalized mixture model. A residual analysis is performed in order to select an appropriate model.

Log-modified Weibull regression models with censored data: Sensitivity and residual analysis

This paper proposes a regression model considering the modified Weibull distribution. This distribution can be used to model bathtub-shaped failure rate functions. Assuming censored data, we consider maximum likelihood and Jackknife estimators for the parameters of the model. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and we also present some ways to perform global influence. Besides, for different parameter settings, sample sizes and censoring percentages, various simulations are performed and the empirical distribution of the modified deviance residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for a martingale-type residual in log-modified Weibull regression models with censored data. Finally, we analyze a real data set under log-modified Weibull regression models. A diagnostic analysis and a model checking based on the modified deviance residual are performed to select appropriate models. (c) 2008 Elsevier B.V. All rights reserved.

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

Sudden onset of log-periodicity and superdiffusion in non-Markovian random walks with amnestically induced persistence: exact results

Random walks can undergo transitions from normal diffusion to anomalous diffusion as some relevant parameter varies, for instance the L,vy index in L,vy flights. Here we derive the Fokker-Planck equation for a two-parameter family of non-Markovian random walks with amnestically induced persistence. We investigate two distinct transitions: one order parameter quantifies log-periodicity and discrete scale invariance in the first moment of the propagator, whereas the second order parameter, known as the Hurst exponent, describes the growth of the second moment. We report numerical and analytical results for six critical exponents, which together completely characterize the properties of the transitions. We find that the critical exponents related to the diffusion-superdiffusion transition are identical in the positive feedback and negative feedback branches of the critical line, even though the former leads to classical superdiffusion whereas the latter gives rise to log-periodic superdiffusion.

An evaluation of the time dependence of the anisotropy of the water diffusion tensor in acute human ischemia

We have performed MRI examinations to determine the water diffusion tensor in the brain of six patients who were admitted to the hospital within 12 h after the onset of cerebral ischemic symptoms. The examinations have been carried out immediately after admission, and thereafter at varying intervals up to 90 days post admission. Maps of the trace of the diffusion tensor, the fractional anisotropy and the lattice index, as well as maps of cerebral blood perfusion parameters, were generated to quantitatively assess the character of the water diffusion tensor in the infarcted area. In patients with significant perfusion deficits and substantial lesion volume changes, four of six cases, our measurements show a monotonic and significant decrease in the diffusion anisotropy within the ischemic lesion as a function of time. We propose that retrospective analysis of this quantity, in combination with brain tissue segmentation and cerebral perfusion maps, may be used in future studies to assess the severity of the ischemic event. (C) 1999 Elsevier Science Inc.