Search results for "Discretization"
showing 10 items of 237 documents
Anisotropic potential of velocity fields in real fluids: Application to the MAST solution of shallow water equations
2013
In the present paper it is first shown that, due to their structure, the general governing equations of uncompressible real fluids can be regarded as an "anisotropic" potential flow problem and closed streamlines cannot occur at any time. For a discretized velocity field, a fast iterative procedure is proposed to order the computational elements at the beginning of each time level, allowing a sequential solution element by element of the advection problem. Some closed circuits could appear due to the discretization error and the elements involved in these circuits could not be ordered. We prove in the paper that the total flux of these not ordered elements goes to zero by refining the compu…
A new discretization for the polarizable continuum model within the domain decomposition paradigm
2016
International audience; We present a new algorithm to solve the polarizable continuum model equation in a framework compatible with the strategy previously developed by us for the conductor-like screening model based on Schwarz’s domain decomposition method (ddCOSMO). The new discretization is systematically improvable and is fully consistent with ddCOSMO so that it reproduces ddCOSMO results for large dielectric constants.
On the equivalence between the Scheduled Relaxation Jacobi method and Richardson's non-stationary method
2017
The Scheduled Relaxation Jacobi (SRJ) method is an extension of the classical Jacobi iterative method to solve linear systems of equations ($Au=b$) associated with elliptic problems. It inherits its robustness and accelerates its convergence rate computing a set of $P$ relaxation factors that result from a minimization problem. In a typical SRJ scheme, the former set of factors is employed in cycles of $M$ consecutive iterations until a prescribed tolerance is reached. We present the analytic form for the optimal set of relaxation factors for the case in which all of them are different, and find that the resulting algorithm is equivalent to a non-stationary generalized Richardson's method. …
A numerical method to calculate the muon relaxation function in the presence of diffusion
2014
We present an accurate and efficient method to calculate the effect of random fluctuations of the local field at the muon, for instance in the case muon diffusion, within the framework of the strong collision approximation. The method is based on a reformulation of the Markovian process over a discretized time base, leading to a summation equation for the muon polarization function which is solved by discrete Fourier transform. The latter is formally analogous, though not identical, to the integral equation of the original continuous-time model, solved by Laplace transform. With real-case parameter values, the solution of the discrete-time strong collision model is found to approximate the …
The multiple slope discontinuity beam element for nonlinear analysis of RC framed structures
2018
The seismic nonlinear response of reinforced concrete structures permits to identify critical zones of an existing structure and to better plan its rehabilitation process. It is obtained by performing finite element analysis using numerical models classifiable into two categories: lumped plasticity models and distributed plasticity models. The present work is devoted to the implementation, in a finite element environment, of an elastoplastic Euler–Bernoulli beam element showing possible slope discontinuities at any position along the beam span, in the framework of a modified lumped plasticity. The differential equation of an Euler–Bernoulli beam element under static loads in presence of mul…
Fast Direct Solver for a Time-harmonic Electromagnetic Problem with an Application
2003
A fast direct solution of a periodic problem derived from the time-harmonic Maxwell’s equations is considered. The problem is discretized by low order hexahedral finite elements proposed by Nedelec. The solver is based on the application of FFT, and it has the computational cost O(N log N). An application to scattering of an electromagnetic wave by a periodic structure is presented.
Harmonic solution of semiconductor transport equations for microwave and millimetre-wave device modelling
2004
The transport equations for charges in a semiconductor have been solved for a periodic voltage excitation by means of a harmonic approach, for modelling of microwave and millimetre-wave active devices. The solution is based on the expansion of the unknown physical quantities in Fourier series in the time domain, and on the discretisation in the space domain. A Waveform-Balance technique in the time domain is used to solve the resulting non-linear equations system. In this way the time step is determined only by Nyquist's sampling requirements at the operating frequency, irrespective of the relaxation times of the semiconductor. This approach allows for a longer time step, and therefore a sh…
Electrodynamics in complex systems
1996
This paper discusses recent theoretical efforts to develop a general and flexible method for the calculation of the field distributions around and inside complex optical systems involving both dielectric and metallic materials. Starting from the usual light-matter coupling Hamiltonian, we derive a self-consistent equation for the optical field in arbitrary optical systems composed of N different subdomains. We show that an appropriate solving procedure based on the real-space discretization of each subdomain raises the present approach to the rank of an accurate predictive numerical scheme. In order to illustrate its applicability, we use this formalism to address challenging problems relat…
Efficient Iterative Solution of Time-harmonic Scattering by Objects in Layered Fluid
2008
We consider the computation of time-harmonic acoustic scattering by sound-soft or elastic objects in layered media. An example of such problem is the scattering by a mine buried in sediment. The computational domain can be tens or hundreds of meters long while the target requires modeling of details smaller than one centimeter. A discretized problem can have several billion degrees of freedom.
A Consistent Boundary/Interior Element Method for Evolutive Elastic Plastic Structural Analysis
1993
A symmetric/sign-definite formulation of the BEM to address the evolutive elastic plastic analysis of structures is presented. A wide class of material models with internal variables and thermodynamic potential is considered. Different energy methods—namely the boundary min-max principle, the Helmholtz free energy and the maximum intrinsic dissipation theorem—axe employed in order to provide the discretization operations by boundary elements and cell elements with inherent variational consistency. The resulting space-discretized equations can be solved by a step-by-step procedure and a predictor/corrector iteration scheme, with corrections operated locally cell-by-cell, just as with the FEM…