A collocation method for solving singular integro-differential equations

This paper describes a collocation method for numerically solving Cauchy-type linear singular integro-differential equations. The numerical method is based on the transformation of the integro-differential equation into an integral equation, and then applying a collocation method to solve the latter. The collocation points are chosen as the Chebyshev nodes. Uniform convergence of the resulting method is then discussed. Numerical examples are presented and solved by the numerical techniques.

A numerical method for solving the dynamic three-dimensional Ericksen-Leslie equations for nematic liquid crystals subject to a strong magnetic field

A finite difference technique, based on a projection method, is developed for solving the dynamic three-dimensional Ericksen-Leslie equations for nematic liquid crystals subject to a strong magnetic field. The governing equations in this situation are derived using primitive variables and are solved using the ideas behind the GENSMAC methodology (Tome and McKee [32]; Tome et al. [34]). The resulting numerical technique is then validated by comparing the numerical solution against an analytic solution for steady three-dimensional flow between two-parallel plates subject to a strong magnetic field. The validated code is then employed to solve channel flow for which there is no analytic solution. (C) 2009 Elsevier B.V. All rights reserved.

An approach for solving the lot-sizing problem of a market-driven foundry

Foundries can be found all over Brazil and they are very important to its economy. In 2008, a mixed integer-programming model for small market-driven foundries was published, attempting to minimize delivery delays. We undertook a study of that model. Here, we present a new approach based on the decomposition of the problem into two sub-problems: production planning of alloys and production planning of items. Both sub-problems are solved using a Lagrangian heuristic based on transferences. An important aspect of the proposed heuristic is its ability to take into account a secondary practice objective solution: the furnace waste. Computational tests show that the approach proposed here is able to generate good quality solutions that outperform prior results. Journal of the Operational Research Society (2010) 61, 108-114. doi:10.1057/jors.2008.151

An implicit technique for solving 3D low Reynolds number moving free surface flows

This paper describes the development of an implicit finite difference method for solving transient three-dimensional incompressible free surface flows. To reduce the CPU time of explicit low-Reynolds number calculations, we have combined a projection method with an implicit technique for treating the pressure on the free surface. The projection method is employed to uncouple the velocity and the pressure fields, allowing each variable to be solved separately. We employ the normal stress condition on the free surface to derive an implicit technique for calculating the pressure at the free surface. Numerical results demonstrate that this modification is essential for the construction of methods that are more stable than those provided by discretizing the free surface explicitly. In addition, we show that the proposed method can be applied to viscoelastic fluids. Numerical results include the simulation of jet buckling and extrudate swell for Reynolds numbers in the range [0.01, 0.5]. (C) 2008 Elsevier Inc. All rights reserved.

A QCD sum rule calculation for the Y(4140) narrow structure

We use the QCD sum rules to evaluate the mass of a possible scalar mesonic state that couples to a molecular D(s)*(D) over bar (s)* current. We find a mass m(Ds)*(Ds)* = (4.14 +/- 0.09) GeV, which is in an excellent agreement with the recently observed Y(4140) charmonium state. We consider the contributions of condensates up to dimension-eight, we work at leading order in alpha(s) and we keep terms which are linear in the strange quark mass m(s). We also consider a molecular D*(D) over bar* current and we obtain m m(D)*(D)* = (4.13 +/- 0.10), around 200 MeV above the mass of the Y(3930) charmonium state. We conclude that it is possible to describe the Y(4140) structure as a D(s)*(D) over bar (s)* molecular state or even as a mixture of D(s)*(D) over bar (s)* and D*(D) over bar* molecular states. (C) 2009 Elsevier B.V. All rights reserved.

QCD sum rule calculation for the charmonium-like structures in the J/psi phi and J/psi omega invariant mass spectra

Using the QCD sum rules we test if the charmonium-like structure Y(4274), observed in the J/psi phi invariant mass spectrum, can be described with a D(s)(D) over bar (s0)(2317)+ h.c. molecular current with J(PC) = 0(-+). We consider the contributions of condensates up to dimension ten and we work at leading order in alpha(s). We keep terms which are linear in the strange quark mass m(s). The mass obtained for such state is mD(s)D(s0) = (4.78 +/- 0.54) GeV. We also consider a molecular 0(-+) D (D) over bar (0)(2400)+ h.c. current and we obtain m(DD0) = (4.55 +/- 0.49) GeV. Our study shows that the newly observed Y(4274) in the J/psi phi invariant mass spectrum can be, considering the uncertainties, described using a molecular charmonium current. (C) 2011 Elsevier B.V. All rights reserved.

The degree of polydispersity in micellar dispersions: A sum rule approach

Statistical properties of a two-dimensional ideal dispersion of polydisperse micelles are derived by analyzing the convergence properties of a sum rule set by mass conservation. Internal micellar degrees of freedom are accounted for by a microscopic model describing small displacements of the constituting amphiphiles with respect to their equilibrium positions. The transfer matrix (TM) method is employed to compute internal micelle partition function. We show that the conditions under which the sum rule is saturated by the largest eigenvalue of the TM determine the value of amphiphile concentration above which the dispersion becomes highly polydisperse and micelle sizes approach a Schultz distribution. (C) 2011 Elsevier B.V. All rights reserved.

A Bayesian Semiparametric Approach for Solving the Discrete Calibration Problem

In this article, we introduce a semi-parametric Bayesian approach based on Dirichlet process priors for the discrete calibration problem in binomial regression models. An interesting topic is the dosimetry problem related to the dose-response model. A hierarchical formulation is provided so that a Markov chain Monte Carlo approach is developed. The methodology is applied to simulated and real data.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.