Search results for "Discretization"
showing 10 items of 237 documents
Controllability method for acoustic scattering with spectral elements
2007
We formulate the Helmholtz equation as an exact controllability problem for the time-dependent wave equation. The problem is then discretized in time domain with central finite difference scheme and in space domain with spectral elements. This approach leads to high accuracy in spatial discretization. Moreover, the spectral element method results in diagonal mass matrices, which makes the time integration of the wave equation highly efficient. After discretization, the exact controllability problem is reformulated as a least-squares problem, which is solved by the conjugate gradient method. We illustrate the method with some numerical experiments, which demonstrate the significant improveme…
On an iterative method for a class of integral equations of the first kind
1987
In this paper, we investigate an iterative method which has been proposed [1] for the numerical solution of a special class of integral equations of the first kind, where one of the essential assumptions is the positivity of the kernel and the given right-hand side. Integral equations of this special type occur in experimental physics, astronomy, medical tomography and other fields where density functions cannot be measured directly, but are related to observable functions via integral equations. In order to take into account the non-negativity of density functions, the proposed iterative scheme was defined in such a way that only non-negative solutions can be approximated. The first part o…
The symmetric boundary element method for unilateral contact problems
2008
Abstract On the basis of the boundary integral equation method, in its symmetric formulation, the frictionless unilateral contact between two elastic bodies has been studied. A boundary discretization by boundary elements leads to an algebraic formulation in the form of a linear complementarity problem. In this paper the process of contact or detachment is obtained through a step by step analysis by using generalized (weighted) quantities as the check elements: the detachment or the contact phenomenon may happen when the weighted traction or the weighted displacement is greater than the weighted cohesion or weighted minimum reference gap, respectively. The applications are performed by usin…
The Line Element-less Method Analysis of orthotropic beam for the De Saint Venant torsion problem
2010
Abstract This paper deals with the extension of a novel numerical technique, labelled line element-less method (LEM), in order to provide approximate solutions of the De Saint Venant torsion problem for orthotropic beams having simply and multiply connected cross-section. A suitable transformation of coordinates allows to take full advantage of the theory of analytic complex functions as in the isotropic case. A complex potential function analytic in all the transformed domain whose real and imaginary parts are related to the shear stress components and to the orthotropic ratio is introduced and expanded in the double-ended Laurent series involving harmonic polynomials. An element-free weak…
A Boundary/Interior Element Discretization Method for the Analysis of Two- and Three-Dimensional Elastic-Plastic Structures
1992
A coupled boundary/interior element method is presented for the analysis of elastic-plastic structures with material models endowed of dual internal variables. The domain field modelling is limited to the only plastic strains and strain-like internal variables, represented by their node values at a set of strain points in each interior element. The formulation, based on a Galerkin-type approach, is variationally consistent and leads to a fully symmetric-definite equation system. The backward difference method is adopted for the step-by-step integration procedure, and each step is addressed by an iterative predictor/corrector solution scheme. The analysis method is expected to be most approp…
Approximation of Elliptic Hemivariational Inequalities
1999
From the previous chapter we know that there exist many important problems in mechanics in which constitutive laws are expressed by means of nonmonotone, possibly multivalued relations (nonmonotone multivalued stress-strain or reaction-displacement relations,e.g). The resulting mathematical model leads to an inclusion type problem involving multivalued nonmonotone mappings or to a substationary type problem for a nonsmooth, nonconvex superpotential expressed in terms of calculus of variation. It is the aim of this chapter to give a detailed study of a discretization of such a type of problems including the convergence analysis. Here we follow closely Miettinen and Haslinger, 1995, Miettinen…
High Order Compact Finite Difference Schemes for A Nonlinear Black-Scholes Equation
2001
A nonlinear Black-Scholes equation which models transaction costs arising in the hedging of portfolios is discretized semi-implicitly using high order compact finite difference schemes. A new compact scheme, generalizing the compact schemes of Rigal [29], is derived and proved to be unconditionally stable and non-oscillatory. The numerical results are compared to standard finite difference schemes. It turns out that the compact schemes have very satisfying stability and non-oscillatory properties and are generally more efficient than the considered classical schemes.
Method of Lines and Finite Difference Schemes with Exact Spectrum for Solving Some Linear Problems of Mathematical Physics
2013
In this paper linear initial-boundary-value problems of mathematical physics with different type boundary conditions BCs and periodic boundary conditions PBCs are studied. The finite difference scheme FDS and the finite difference scheme with exact spectrum FDSES are used for the space discretization. The solution in the time is obtained analytically and numerically, using the method of lines and continuous and discrete Fourier methods.
Perimeter symmetrization of some dynamic and stationary equations involving the Monge-Ampère operator
2017
We apply the perimeter symmetrization to a two-dimensional pseudo-parabolic dynamic problem associated to the Monge-Ampere operator as well as to the second order elliptic problem which arises after an implicit time discretization of the dynamical equation. Curiously, the dynamical problem corresponds to a third order operator but becomes a singular second order parabolic equation (involving the 3-Laplacian operator) in the class of radially symmetric convex functions. Using symmetrization techniques some quantitative comparison estimates and several qualitative properties of solutions are given.
Line element-less method (LEM) for beam torsion solution (truly no-mesh method)
2008
In this paper a new numerical method for finding approximate solutions of the torsion problem is proposed. The method takes full advantage of the theory of analytic complex function. A new potential function directly in terms of shear stresses is proposed and expanded in the double-ended Laurent series involving harmonic polynomials. A novel element-free weak form procedure, labelled Line Element-Less Method (LEM), has been developed imposing that the square of the net flux across the border is minimum with respect to coefficients expansion. Numerical implementation of the LEM results in systems of linear algebraic equations involving symmetric and positive-definite matrices without resorti…