Search results for "Numerical"

A degenerating convection-diffusion system modelling froth flotation with drainage

2022

Abstract Froth flotation is a common unit operation used in mineral processing. It serves to separate valuable mineral particles from worthless gangue particles in finely ground ores. The valuable mineral particles are hydrophobic and attach to bubbles of air injected into the pulp. This creates bubble-particle aggregates that rise to the top of the flotation column where they accumulate to a froth or foam layer that is removed through a launder for further processing. At the same time, the hydrophilic gangue particles settle and are removed continuously. The drainage of liquid due to capillarity is essential for the formation of a stable froth layer. This effect is included into a previous…

A generalized Newton iteration for computing the solution of the inverse Henderson problem

2020

We develop a generalized Newton scheme IHNC for the construction of effective pair potentials for systems of interacting point-like particles.The construction is made in such a way that the distribution of the particles matches a given radial distribution function. The IHNC iteration uses the hypernetted-chain integral equation for an approximate evaluation of the inverse of the Jacobian of the forward operator. In contrast to the full Newton method realized in the Inverse Monte Carlo (IMC) scheme, the IHNC algorithm requires only a single molecular dynamics computation of the radial distribution function per iteration step, and no further expensive cross-correlations. Numerical experiments…

Approximation and quasicontinuity of Besov and Triebel–Lizorkin functions

2016

We show that, for $0<s<1$, $0<p<\infty$, $0<q<\infty$, Haj\l asz-Besov and Haj\l asz-Triebel-Lizorkin functions can be approximated in the norm by discrete median convolutions. This allows us to show that, for these functions, the limit of medians, \[ \lim_{r\to 0}m_u^\gamma(B(x,r))=u^*(x), \] exists quasieverywhere and defines a quasicontinuous representative of $u$. The above limit exists quasieverywhere also for Haj\l asz functions $u\in M^{s,p}$, $0<s\le 1$, $0<p<\infty$, but approximation of $u$ in $M^{s,p}$ by discrete (median) convolutions is not in general possible.

On a global superconvergence of the gradient of linear triangular elements

1987

Abstract We study a simple superconvergent scheme which recovers the gradient when solving a second-order elliptic problem in the plane by the usual linear elements. The recovered gradient globally approximates the true gradient even by one order of accuracy higher in the L 2 -norm than the piecewise constant gradient of the Ritz—Galerkin solution. A superconvergent approximation to the boundary flux is presented as well.

A strain-difference-based nonlocal elasticity model

2004

Abstract A two-component local/nonlocal constitutive model for (macroscopically) inhomogeneous linear elastic materials (but constant internal length) is proposed, in which the stress is the sum of the local stress and a nonlocal-type stress expressed in terms of the strain difference field, hence identically vanishing in the case of uniform strain. Attention is focused upon the particular case of piecewise homogeneous material. The proposed model is thermodynamically consistent with a suitable free energy potential. It constitutes an improved form of the Vermeer and Brinkgreve [A new effective nonlocal strain measure for softening plasticity. In: Chambon, R., Desrues, J., Vardulakis, I. (E…

Energy dissipative characteristic schemes for the diffusive Oldroyd-B viscoelastic fluid

2015

Domain decomposition in the symmetric boundary element analysis

2002

Recent developments in the symmetric boundary element method (SBEM) have shown a clear superiority of this formulation over the collocation method. Its competitiveness has been tested in comparison to the finite element method (FEM) and is manifested in several engineering problems in which internal boundaries are present, i.e. those in which the body shows a jump in the physical characteristics of the material and in which an appropriate study of the response must be used. When we work in the ambit of the SBE formulation, the body is subdivided into macroelements characterized by some relations which link the interface boundary unknowns to the external actions. These relations, valid for e…

Elastic plastic analysis iterative solution

1998

The step-by-step analysis of finite element elastic plastic structures subjected to an assigned (quasi-static) loading history, is considered; it identifies with the well-known sequence of linear complementarity problems. An iterative technique devoted to solve the relevant linear complementarity problem is presented. It is based on the recursive solution of a suitable linear complementarity problem, deduced from the relevant one and easier than it. The procedure convergency is proved. Some noticing particular cases are examined. The physical meaning of the procedure is shown to be a plastic relaxation. The suitable numerical ranges for some check parameter values, to be utilized in the app…

Multidomain boundary integral formulation for piezoelectric materials fracture mechanics

2001

Abstract A boundary element method and its numerical implementation for the analysis of piezoelectric materials are presented with the aim to exploit their features in linear electroelastic fracture mechanics. The problem is formulated employing generalized displacements, that is displacements and electric potential, and generalized tractions, that is tractions and electric displacement. The generalized displacements boundary integral equation is obtained by using the closed form of the piezoelasticity fundamental solutions. These are derived through a displacement based modified Lekhnitskii’s functions approach. The multidomain boundary element technique is implemented to achieve the numer…

Numerical Algorithms Based on Characteristic Domain Decomposition for Obstacle Problems

1997

A new numerical solution algorithm for obstacle problems is proposed, where the characteristic domain decomposition into active and inactive subdomains separated by the free boundary is approximated by a Schwarz method. Such an approach gives an opportunity to apply fast linear system solvers to genuinely non-linear obstacle problems. Other solution algorithms, like projected relaxation methods and active set strategies, are compared to the new solution algorithm. Numerical experiments related to the elastoplastic torsion problem are included showing the efficiency of the new approach.