Hybrid mesh generation with embedded surfaces using a multiple marching direction approach

This paper proposes a method for the creation of hybrid meshes with embedded surfaces for viscous flow simulations as an extension of the multiple marching direction approach (AIAA J. 2007; 45 (1):162–167). The multiple marching direction approach enables to place semi‐structured elements around singular points, where valid semi‐structured elements cannot be placed using conventional hybrid mesh generation methods. This feature is discussed first with a couple of examples. Elements sometimes need to be clustered inside a computational domain to obtain more accurate results. For example, solution features, such as shocks, vortex cores and wake regions, can be extracted during the process of adaptive mesh generation. These features can be represented as surface meshes embedded in a computational domain. Semi‐structured elements can be placed around the embedded surface meshes using the multiple marching direction approach with a pretreatment method. Tetrahedral elements can be placed easily instead. A couple of results are presented to demonstrate the capability of the mesh generation method. Copyright © 2008 John Wiley & Sons, Ltd.     

On a discretization process of fractional-order Logistic differential equation

In this work we introduce a discretization process to discretize fractional-order differential equations. First of all, we consider the fractional-order Logistic differential equation then, we consider the corresponding fractional-order Logistic differential equation with piecewise constant arguments and we apply the proposed discretization on it. The stability of the fixed points of the resultant dynamical system and the Lyapunov exponent are investigated. Finally, we study some dynamic behavior of the resultant systems such as bifurcation and chaos.     

Phase diagram of the 3D bimodal random-field Ising model

The one-parametric Wang-Landau (WL) method is implemented together with an extrapolation scheme to yield approximations of the two-dimensional (exchange-energy, field-energy) density of states (DOS) of the 3D bimodal random-field Ising model (RFIM). The present approach generalizes our earlier WL implementations, by handling the final stage of the WL process as an entropic sampling scheme, appropriate for the recording of the required two-parametric histograms. We test the accuracy of the proposed extrapolation scheme and then apply it to study the size-shift behavior of the phase diagram of the 3D bimodal RFIM. We present a finite-size converging approach and a well-behaved sequence of estimates for the critical disorder strength. Their asymptotic shift-behavior yields the critical disorder strength and the associated correlation length's exponent, in agreement with previous estimates from ground-state studies of the model.     

Infinitely Many Critical Points of Non-Differentiable Functions and Applications to a Neumann-Type Problem Involving the p-Laplacian

For a family of functionals in a Banach space, which are possibly non-smooth and depend also on a positive real parameter, the existence of a sequence of critical points (according to Motreanu and Panagiotopoulos (“Minimax Theorems and Qualitative Properties of the Solutions of Hemivariational Inequalities,” Nonconvex Optimization Applications, Vol. 29, Kluwer, Dordrecht, 1998, Chap. 3)) is established by mainly adapting a new technique due to Ricceri (2000, J. Comput. Appl. Math.113, 401-410). Two applications are then presented. Both of them treat the Neumann problem for an elliptic variational-hemivariational inequality with p-Laplacian.     

On dynamics of vertices of locally connected polynomial Julia sets

Let be a polynomial whose Julia set is locally connected. Then a non-preperiodic non-precritical vertex of must have the limit set which coincides with the limit set of an appropriately chosen recurrent critical point of . In particular, if all critical points of are non-recurrent then all vertices of are preperiodic or precritical.

     

Positivity preserving finite element approximation

We consider finite element operators defined on ``rough' functions in a bounded polyhedron in . Insisting on preserving positivity in the approximations, we discover an intriguing and basic difference between approximating functions which vanish on the boundary of and approximating general functions which do not. We give impossibility results for approximation of general functions to more than first order accuracy at extreme points of . We also give impossibility results about invariance of positive operators on finite element functions. This is in striking contrast to the well-studied case without positivity.

     

Hausdorff Dimension of Sets of Generic Points for Gibbs Measures

For a translation invariant Gibbs measure on the configuration space X of a lattice finite spin system, we consider the set X of generic points. Using a Breiman type convergence theorem on the set X of generic points of an arbitrary translation invariant probability measure on X, we evaluate the Hausdorff dimension of the set X with respect to any metric out of a wide class of scale metrics on X (including Billingsley metrics generated by Gibbs measures).     

On the Stability of the O(N)-Invariant and the Cubic-Invariant Three-Dimensional N-Component Renormalization-Group Fixed Points in the Hierarchical Approximation

We compute renormalization-group fixed points and their spectrum in an ultralocal approximation. We study a case of two competing nontrivial fixed points for a three-dimensional real N-component field: the O(N)-invariant fixed point vs. the cubic-invariant fixed point. We compute the critical value N _c of the cubic ⁴-perturbation at the O(N)-fixed point. The O(N)-fixed point is stable under a cubic ⁴-perturbation below N _c; above N _c it is unstable. The Critical value comes out as 2.219435<N _c<2.219436 in the ultralocal approximation. We also compute the critical value of N at the cubic invariant fixed point. Within the accuracy of our computations, the two values coincide.     

On a Primal-Dual Analytic Center Cutting Plane Method for Variational Inequalities

We present an algorithm for variational inequalities VI( , Y) that uses a primal-dual version of the Analytic Center Cutting Plane Method. The point-to-set mapping is assumed to be monotone, or pseudomonotone. Each computation of a new analytic center requires at most four Newton iterations, in theory, and in practice one or sometimes two. Linear equalities that may be included in the definition of the set Y are taken explicitly into account.We report numerical experiments on several well—known variational inequality problems as well as on one where the functional results from the solution of large subproblems. The method is robust and competitive with algorithms which use the same information as this one.     

Numerical algorithms for uniform Airy-type asymptotic expansions

Airy-type asymptotic representations of a class of special functions are considered from a numerical point of view. It is well known that the evaluation of the coefficients of the asymptotic series near the transition point is a difficult problem. We discuss two methods for computing the asymptotic series. One method is based on expanding the coefficients of the asymptotic series in Maclaurin series. In the second method we consider auxiliary functions that can be computed more efficiently than the coefficients in the first method, and we do not need the tabulation of many coefficients. The methods are quite general, but the paper concentrates on Bessel functions, in particular on the differential equation of the Bessel functions, which has a turning point character when order and argument of the Bessel functions are equal.