Search results for "Numerical"
showing 10 items of 2002 documents
Initial strain effects in multilayer composite laminates
2001
A boundary integral formulation for the analysis of stress fields induced in composite laminates by initial strains, such as may be due to temperature changes and moisture absorption is presented. The study is formulated on the basis of the theory of generalized orthotropic thermo-elasticity and the governing integral equations are directly deduced through the generalized reciprocity theorem. A suitable expression of the problem fundamental solutions is given for use in computations. The resulting linear system of algebraic equations is obtained by the boundary element method and stress interlaminar distributions in the boundary-layer are calculated by using a boundary only discretization. …
BIEM-based variational principles for elastoplasticity with unilateral contact boundary conditions
1998
The structural step problem for elastic-plastic internal-variable materials is addressed in the presence of frictionless unilateral contact conditions. Basing on the BIEM (boundary integral equation method) and making use of deformation-theory plasticity (through the backward-difference method of computational plasticity), two variational principles are shown to characterize the solution to the step problem: one is a stationarity principle having as unknowns all the problem variables, the other is a saddle-point principle having as unknowns the increments of the boundary tractions and displacements, along with the plastic strain increments in the domain. The discretization by boundary and i…
A Mesh-free Particle Method for Transient Full-wave Simulation
2007
A mesh-free particle method is presented for electromagnetic (EM) transient simulation. The basic idea is to obtain numerical solutions for the partial differential equations describing the EM problem in time domain, by using a set of particles, considered as spatial interpolation points of the field variables, arbitrarily placed in the problem domain and by avoiding the use of a regular mesh. Irregular problems geometry with diffused non-homogeneous media can be modeled only with an initial set of arbitrarily distributed particles. The time dependence is accounted for with an explicit finite difference scheme. Moreover the particle discretization can be improved during the process time ste…
A new solver for incompressible non-isothermal flows in natural and mixed convection over unstructured grids
2022
Abstract In the present paper we propose a new numerical methodology for the solution of 2D non-isothermal incompressible flows for natural and mixed convection in irregular geometries. The governing equations are the Incompressible Navier-Stokes Equations and the Energy Conservation Equation. Fluid velocity and temperature are coupled in the buoyancy term of the momentum equations according to the Oberbeck–Boussinesq approximation. The governing equations are discretized over unstructured triangular meshes satisfying the Delaunay property. Thanks to the Oberbeck–Boussinesq hypothesis, the flow and energy problems are solved in an uncoupled way, and two fractional time step procedures are s…
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…
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…
High-fidelity analysis of multilayered shells with cut-outs via the discontinuous Galerkin method
2021
Abstract A novel numerical method for the analysis of multilayered shells with cut-outs is presented. In the proposed approach, the shell geometry is represented via either analytical functions or NURBS parametrizations , while generally-shaped cut-outs are defined implicitly within the shell modelling domain via a level set function . The multilayered shell problem is addressed via the Equivalent-Single-Layer approach whereby high-order polynomial functions are employed to approximate the covariant components of the displacement field throughout the shell thickness. The shell governing equations are then derived from the Principle of Virtual Displacements of three-dimensional elasticity an…
Guaranteed and computable error bounds for approximations constructed by an iterative decoupling of the Biot problem
2021
The paper is concerned with guaranteed a posteriori error estimates for a class of evolutionary problems related to poroelastic media governed by the quasi-static linear Biot equations. The system is decoupled by employing the fixed-stress split scheme, which leads to an iteratively solved semi-discrete system. The error bounds are derived by combining a posteriori estimates for contractive mappings with functional type error control for elliptic partial differential equations. The estimates are applicable to any approximation in the admissible functional space and are independent of the discretization method. They are fully computable, do not contain mesh-dependent constants, and provide r…
ADI schemes for valuing European options under the Bates model
2018
Abstract This paper is concerned with the adaptation of alternating direction implicit (ADI) time discretization schemes for the numerical solution of partial integro-differential equations (PIDEs) with application to the Bates model in finance. Three different adaptations are formulated and their (von Neumann) stability is analyzed. Ample numerical experiments are provided for the Bates PIDE, illustrating the actual stability and convergence behaviour of the three adaptations.
Dynamic shakedown of structures with variable appended masses and subjected to repeated excitations
1996
Elastic shakedown for discrete, or finite-element discretized, structures subjected to combinations of static and time-variable loads is addressed in the hypothesis of elastic-perfectly plastic material behavior. The static load is conceived as the weight of an additional mass appended to the structure, whereas the time-variable load is conceived as an unknown sequence of excitations belonging to a specified domain, with intervals between subsequent excitations during which the structure is considered as being motionless. It is shown that, in the plane of the static and time-variable load parameters, the structure's dynamic shakedown domain is nonconvex and that its boundary curve generally…