Simulation results and applications of an advection bounded scheme to practical flows

This paper reports experiments on the use of a recently introduced advection bounded upwinding scheme, namely TOPUS (Computers & Fluids 57 (2012) 208-224), for flows of practical interest. The numerical results are compared against analytical, numerical and experimental data and show good agreement with them. It is concluded that the TOPUS scheme is a competent, powerful and generic scheme for complex flow phenomena.

Lê–Greuel type formula for the Euler obstruction and applications

The Euler obstruction of a function f can be viewed as a generalization of the Milnor number for functions defined on singular spaces. In this work, using the Euler obstruction of a function, we establish several Lê–Greuel type formulas for germs f:(X,0)→(C,0) and g:(X,0)→(C,0). We give applications when g is a generic linear form and when f and g have isolated singularities.

Dirac fermions in strong electric field and quantum transport in graphene

Our previous results on the nonperturbative calculations of the mean current and of the energy-momentum tensor in QED with the T-constant electric field are generalized to arbitrary dimensions. The renormalized mean values are found, and the vacuum polarization contributions and particle creation contributions to these mean values are isolated in the large T limit; we also relate the vacuum polarization contributions to the one-loop effective Euler-Heisenberg Lagrangian. Peculiarities in odd dimensions are considered in detail. We adapt general results obtained in 2 + 1 dimensions to the conditions which are realized in the Dirac model for graphene. We study the quantum electronic and energy transport in the graphene at low carrier density and low temperatures when quantum interference effects are important. Our description of the quantum transport in the graphene is based on the so-called generalized Furry picture in QED where the strong external field is taken into account nonperturbatively; this approach is not restricted to a semiclassical approximation for carriers and does not use any statistical assumptions inherent in the Boltzmann transport theory. In addition, we consider the evolution of the mean electromagnetic field in the graphene, taking into account the backreaction of the matter field to the applied external field. We find solutions of the corresponding Dirac-Maxwell set of equations and with their help we calculate the effective mean electromagnetic field and effective mean values of the current and the energy-momentum tensor. The nonlinear and linear I-V characteristics experimentally observed in both low-and high-mobility graphene samples are quite well explained in the framework of the proposed approach, their peculiarities being essentially due to the carrier creation from the vacuum by the applied electric field. DOI: 10.1103/PhysRevD.86.125022

On the Milnor fibre and L numbers of semi-weighted homogeneous arrangements

Inthispaperwestudygermsofpolynomialsformedbytheproductofsemi-weighted homogeneous polynomials of the same type, which we call semi-weighted homogeneous arrangements. It is shown how the L numbers of such polynomials are computed using only their weights and degree of homogeneity. A key point of the main theorem is to find the number called polar ratio of this polynomial class. An important consequence is the description of the Euler characteristic of the Milnor fibre of such arrangements only depending on their weights and degree of homogeneity. The constancy of the L numbers in families formed by such arrangements is shown, with the deformed terms having weighted degree greater than the weighted degree of the initial germ. Moreover, using the results of Massey applied to families of function germs, we obtain the constancy of the homology of the Milnor fibre in this family of semi-weighted homogeneous arrangements.

A novel self consistent calculation approach for the capacitance-voltage characteristics of semiconductor quantum wire transistors based on a split-gate configuration

Purpose - The purpose of this paper is to develop an efficient numerical algorithm for the self-consistent solution of Schrodinger and Poisson equations in one-dimensional systems. The goal is to compute the charge-control and capacitance-voltage characteristics of quantum wire transistors. Design/methodology/approach - The paper presents a numerical formulation employing a non-uniform finite difference discretization scheme, in which the wavefunctions and electronic energy levels are obtained by solving the Schrodinger equation through the split-operator method while a relaxation method in the FTCS scheme ("Forward Time Centered Space") is used to solve the two-dimensional Poisson equation. Findings - The numerical model is validated by taking previously published results as a benchmark and then applying them to yield the charge-control characteristics and the capacitance-voltage relationship for a split-gate quantum wire device. Originality/value - The paper helps to fulfill the need for C-V models of quantum wire device. To do so, the authors implemented a straightforward calculation method for the two-dimensional electronic carrier density n(x,y). The formulation reduces the computational procedure to a much simpler problem, similar to the one-dimensional quantization case, significantly diminishing running time.

Automatic design of decision-tree induction algorithms tailored to flexible-receptor docking data

Background: This paper addresses the prediction of the free energy of binding of a drug candidate with enzyme InhA associated with Mycobacterium tuberculosis. This problem is found within rational drug design, where interactions between drug candidates and target proteins are verified through molecular docking simulations. In this application, it is important not only to correctly predict the free energy of binding, but also to provide a comprehensible model that could be validated by a domain specialist. Decision-tree induction algorithms have been successfully used in drug-design related applications, specially considering that decision trees are simple to understand, interpret, and validate. There are several decision-tree induction algorithms available for general-use, but each one has a bias that makes it more suitable for a particular data distribution. In this article, we propose and investigate the automatic design of decision-tree induction algorithms tailored to particular drug-enzyme binding data sets. We investigate the performance of our new method for evaluating binding conformations of different drug candidates to InhA, and we analyze our findings with respect to decision tree accuracy, comprehensibility, and biological relevance. Results: The empirical analysis indicates that our method is capable of automatically generating decision-tree induction algorithms that significantly outperform the traditional C4.5 algorithm with respect to both accuracy and comprehensibility. In addition, we provide the biological interpretation of the rules generated by our approach, reinforcing the importance of comprehensible predictive models in this particular bioinformatics application. Conclusions: We conclude that automatically designing a decision-tree algorithm tailored to molecular docking data is a promising alternative for the prediction of the free energy from the binding of a drug candidate with a flexible-receptor.

Insights into the morphology of the Serra da Cangalha impact structure from geophysical modeling

Forward modeling is commonly applied to gravity field data of impact structures to determine the main gravity anomaly sources. In this context, we have developed 2.5-D gravity models of the Serra da Cangalha impact structure for the purpose of investigating geological bodies/structures underneath the crater. Interpretation of the models was supported by ground magnetic data acquired along profiles, as well as by high resolution aeromagnetic data. Ground magnetic data reveal the presence of short-wavelength anomalies probably related to shallow magnetic sources that could have been emplaced during the cratering process. Aeromagnetic data show that the basement underneath the crater occurs at an average depth of about 1.9 km, whereas in the region beneath the central uplift it is raised to 0.51 km below the current surface. These depths are also supported by 2.5-D gravity models showing a gentle relief for the basement beneath the central uplift area. Geophysical data were used to provide further constraints for numeral modeling of crater formation that provided important information on the structural modification that affected the rocks underneath the crater, as well as on shock-induced modifications of target rocks. The results showed that the morphology is consistent with the current observations of the crater and that Serra da Cangalha was formed by a meteorite of approximately 1.4 km diameter striking at 12 km s-1.

URANS calculations for smooth circular cylinder flow in a wide range of Reynolds Numbers: solution verification and validation

The flow around circular smooth fixed cylinder in a large range of Reynolds numbers is considered in this paper. In order to investigate this canonical case, we perform CFD calculations and apply verification & validation (V&V) procedures to draw conclusions regarding numerical error and, afterwards, assess the modeling errors and capabilities of this (U)RANS method to solve the problem. Eight Reynolds numbers between Re = 10 and Re 5 x 10(5) will be presented with, at least, four geometrically similar grids and five discretization in time for each case (when unsteady), together with strict control of iterative and round-off errors, allowing a consistent verification analysis with uncertainty estimation. Two-dimensional RANS, steady or unsteady, laminar or turbulent calculations are performed. The original 1994 k - omega SST turbulence model by Menter is used to model turbulence. The validation procedure is performed by comparing the numerical results with an extensive set of experimental results compiled from the literature. [DOI: 10.1115/1.4007571]

Haemodynamic effects of aliskiren in decompensated severe heart failure

Aim: The renin-angiotensin-aldosterone system (RAAS) has dual pathways to angiotensin II production; therefore, multiple blockages may be useful in heart failure. In this study, we evaluated the short-term haemodynamic effects of aliskiren, a direct renin inhibitor, in patients with decompensated severe heart failure who were also taking angiotensin-converting enzyme ( ACE) inhibitors. Materials and methods: A total of 16 patients (14 men, two women, mean age: 60.3 years) were enrolled in the study. The inclusion criteria included hospitalisation due to decompensated heart failure, ACE inhibitor use, and an ejection fraction < 40% (mean: 21.9 +/- 6.7%). The exclusion criteria were: creatinine > 2.0 mg/dl, cardiac pacemaker, serum K+ > 5.5 mEq/l, and systolic blood pressure < 70 mmHg. Patients either received 150 mg/d aliskiren for 7 days (aliskiren group, n = 10) or did not receive aliskiren (control group, n = 6). Primary end points were systemic vascular resistance and cardiac index values. Repeated-measures analysis of variance (ANOVA) was used to assess variables before and after intervention. A two-sided p-value < 0.05 was considered statistically significant. Results: Compared to pre-intervention levels, systemic vascular resistance was reduced by 20.4% in aliskiren patients, but it increased by 2.9% in control patients (p = 0.038). The cardiac index was not significantly increased by 19.0% in aliskiren patients, but decreased by 8.4% in control patients (p = 0.127). No differences in the pulmonary capillary or systolic blood pressure values were observed between the groups. Conclusion: Aliskiren use reduced systemic vascular resistance in patients with decompensated heart failure taking ACE inhibitors.

Transiting exoplanets from the CoRoT space mission XXI. CoRoT-19b: a low density planet orbiting an old inactive F9V-star

Context. Observations of transiting extrasolar planets are of key importance to our understanding of planets because their mass, radius, and mass density can be determined. These measurements indicate that planets of similar mass can have very different radii. For low-density planets, it is generally assumed that they are inflated owing to their proximity to the host-star. To determine the causes of this inflation, it is necessary to obtain a statistically significant sample of planets with precisely measured masses and radii. Aims. The CoRoT space mission allows us to achieve a very high photometric accuracy. By combining CoRoT data with high-precision radial velocity measurements, we derive precise planetary radii and masses. We report the discovery of CoRoT-19b, a gas-giant planet transiting an old, inactive F9V-type star with a period of four days. Methods. After excluding alternative physical configurations mimicking a planetary transit signal, we determine the radius and mass of the planet by combining CoRoT photometry with high-resolution spectroscopy obtained with the echelle spectrographs SOPHIE, HARPS, FIES, and SANDIFORD. To improve the precision of its ephemeris and the epoch, we observed additional transits with the TRAPPIST and Euler telescopes. Using HARPS spectra obtained during the transit, we then determine the projected angle between the spin of the star and the orbit of the planet. Results. We find that the host star of CoRoT-19b is an inactive F9V-type star close to the end of its main-sequence life. The host star has a mass M-* = 1.21 +/- 0.05 M-circle dot and radius R-* = 1.65 +/- 0.04 R-circle dot. The planet has a mass of M-P = 1.11 +/- 0.06 M-Jup and radius of R-P = 1.29 +/- 0.03 R-Jup. The resulting bulk density is only rho = 0.71 +/- 0.06 g cm (3), which is much lower than that for Jupiter. Conclusions. The exoplanet CoRoT-19b is an example of a giant planet of almost the same mass as Jupiter but a approximate to 30% larger radius.

A marker-and-cell approach to free surface 2-D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front-tracking method. The velocity field is computed using a finite-difference discretization of a modification of the NavierStokes equations. These equations together with the continuity equation are solved for the two-dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright (c) 2012 John Wiley & Sons, Ltd.

The analytic torsion of a disc

In this article, we study the Reidemeister torsion and the analytic torsion of the m dimensional disc, with the Ray and Singer homology basis (Adv Math 7:145-210, 1971). We prove that the Reidemeister torsion coincides with a power of the volume of the disc. We study the additional terms arising in the analytic torsion due to the boundary, using generalizations of the Cheeger-Muller theorem. We use a formula proved by Bruning and Ma (GAFA 16:767-873, 2006) that predicts a new anomaly boundary term beside the known term proportional to the Euler characteristic of the boundary (Luck, J Diff Geom 37:263-322, 1993). Some of our results extend to the case of the cone over a sphere, in particular we evaluate directly the analytic torsion for a cone over the circle and over the two sphere. We compare the results obtained in the low dimensional cases. We also consider a different formula for the boundary term given by Dai and Fang (Asian J Math 4:695-714, 2000), and we compare the results. The results of these work were announced in the study of Hartmann et al. (BUMI 2:529-533, 2009).

EQUIVALENCE OF REAL MILNOR FIBRATIONS FOR QUASI-HOMOGENEOUS SINGULARITIES

We show that for real quasi-homogeneous singularities f : (R-m, 0) -> (R-2, 0) with isolated singular point at the origin, the projection map of the Milnor fibration S-epsilon(m-1) \ K-epsilon -> S-1 is given by f/parallel to f parallel to. Moreover, for these singularities the two versions of the Milnor fibration, on the sphere and on a Milnor tube, are equivalent. In order to prove this, we show that the flow of the Euler vector field plays and important role. In addition, we present, in an easy way, a characterization of the critical points of the projection (f/parallel to f parallel to) : S-epsilon(m-1) \ K-epsilon -> S-1.

A continuously differentiable upwinding scheme for the simulation of fluid flow problems

This paper deals with the numerical solution of complex fluid dynamics problems using a new bounded high resolution upwind scheme (called SDPUS-C1 henceforth), for convection term discretization. The scheme is based on TVD and CBC stability criteria and is implemented in the context of the finite volume/difference methodologies, either into the CLAWPACK software package for compressible flows or in the Freeflow simulation system for incompressible viscous flows. The performance of the proposed upwind non-oscillatory scheme is demonstrated by solving two-dimensional compressible flow problems, such as shock wave propagation and two-dimensional/axisymmetric incompressible moving free surface flows. The numerical results demonstrate that this new cell-interface reconstruction technique works very well in several practical applications. (C) 2012 Elsevier Inc. All rights reserved.

Numerical assessment of stability of interface discontinuous finite element pressure spaces

The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.