Search results for "Numerical Analysis"
showing 10 items of 883 documents
Generalized wave propagation problems and discrete exterior calculus
2018
We introduce a general class of second-order boundary value problems unifying application areas such as acoustics, electromagnetism, elastodynamics, quantum mechanics, and so on, into a single framework. This also enables us to solve wave propagation problems very efficiently with a single software system. The solution method precisely follows the conservation laws in finite-dimensional systems, whereas the constitutive relations are imposed approximately. We employ discrete exterior calculus for the spatial discretization, use natural crystal structures for three-dimensional meshing, and derive a “discrete Hodge” adapted to harmonic wave. The numerical experiments indicate that the cumulat…
Residual a posteriori error estimation for frictional contact with Nitsche method
2023
We consider frictional contact problems in small strain elasticity discretized with finite elements and Nitsche method. Both bilateral and unilateral contact problems are taken into account, as well as both Tresca and Coulomb models for the friction. We derive residual a posteriori error estimates for each friction model, following [Chouly et al, IMA J. Numer. Anal. (38) 2018, pp. 921-954]. For the incomplete variant of Nitsche, we prove an upper bound for the dual norm of the residual, for Tresca and Coulomb friction, without any extra regularity and without a saturation assumption. Numerical experiments allow to assess the accuracy of the estimates and their interest for adaptive meshing …
Rocking behaviour of multi-block columns subjected to pulse-typeground motion accelerations
2016
Ancient columns, made with a variety of materials such as marble, granite, stone or masonry are an important part of theEuropean cultural heritage. In particular columns of ancient temples in Greece and Sicily which support only the architrave arecharacterized by small axial load values. This feature together with the slenderness typical of these structural members clearlyhighlights as the evaluation of the rocking behaviour is a key aspect of their safety assessment and maintenance. It has to be notedthat the rocking response of rectangular cross-sectional columns modelled as monolithic rigid elements, has been widely investigatedsince the first theoretical study carried out by Housner (19…
Improving the Downwind Sail Design Process by Means of a Novel FSI Approach
2021
The process of designing a sail can be a challenging task because of the difficulties in predicting the real aerodynamic performance. This is especially true in the case of downwind sails, where the evaluation of the real shapes and aerodynamic forces can be very complex because of turbulent and detached flows and the high-deformable behavior of structures. Of course, numerical methods are very useful and reliable tools to investigate sail performances, and their use, also as a result of the exponential growth of computational resources at a very low cost, is spreading more and more, even in not highly competitive fields. This paper presents a new methodology to support sail designers in ev…
Numerical methods for a nonlinear impact model: A comparative study with closed-form corrections
2011
A physically based impact model-already known and exploited in the field of sound synthesis-is studied using both analytical tools and numerical simulations. It is shown that the Hamiltonian of a physical system composed of a mass impacting on a wall can be expressed analytically as a function of the mass velocity during contact. Moreover, an efficient and accurate approximation for the mass outbound velocity is presented, which allows to estimate the Hamiltonian at the end of the contact. Analytical results are then compared to numerical simulations obtained by discretizing the system with several numerical methods. It is shown that, for some regions of the parameter space, the trajectorie…
Reverse Catmull-Clark Subdivision
2006
Reverse subdivision consists in constructing a coarse mesh of a model from a finer mesh of this same model. In this paper, we give formulas for reverse Catmull-Clark subdivision. These formulas allow the constructing of a coarse mesh for almost all meshes. The condition for being able to apply these formulas is that the mesh to be reversed must be generated by the subdivision of a coarse mesh. Except for this condition, the mesh can be arbitrary. Vertices can be regular or extraordinary and the mesh itself can be arbitrary (triangular, quadrilateral…).
Numerical Recovery of Source Singularities via the Radiative Transfer Equation with Partial Data
2013
The inverse source problem for the radiative transfer equation is considered, with partial data. Here we demonstrate numerical computation of the normal operator $X_{V}^{*}X_{V}$ where $X_{V}$ is the partial data solution operator to the radiative transfer equation. The numerical scheme is based in part on a forward solver designed by F. Monard and G. Bal. We will see that one can detect quite well the visible singularities of an internal optical source $f$ for generic anisotropic $k$ and $\sigma$, with or without noise added to the accessible data $X_{V}f$. In particular, we use a truncated Neumann series to estimate $X_{V}$ and $X_{V}^{*}$, which provides a good approximation of $X_{V}^{*…
A posteriori error estimates for time-dependent reaction-diffusion problems based on the Payne-Weinberger inequality
2015
We consider evolutionary reaction-diffusion problem with mixed Dirichlet--Robin boundary conditions. For this class of problems, we derive two-sided estimates of the distance between any function in the admissible energy space and exact solution of the problem. The estimates (majorants and minorants) are explicitly computable and do not contain unknown functions or constants. Moreover, it is proved that the estimates are equivalent to the energy norm of the deviation from the exact solution.
IMEX schemes for pricing options under jump–diffusion models
2014
We propose families of IMEX time discretization schemes for the partial integro-differential equation derived for the pricing of options under a jump-diffusion process. The schemes include the families of IMEX-midpoint, IMEX-CNAB and IMEX-BDF2 schemes. Each family is defined by a convex combination parameter [email protected]?[0,1], which divides the zeroth-order term due to the jumps between the implicit and explicit parts in the time discretization. These IMEX schemes lead to tridiagonal systems, which can be solved extremely efficiently. The schemes are studied through Fourier stability analysis and numerical experiments. It is found that, under suitable assumptions and time step restric…
LOCAL CONTROL OF SOUND IN STOCHASTIC DOMAINS BASED ON FINITE ELEMENT MODELS
2011
A numerical method for optimizing the local control of sound in a stochastic domain is developed. A three-dimensional enclosed acoustic space, for example, a cabin with acoustic actuators in given locations is modeled using the finite element method in the frequency domain. The optimal local noise control signals minimizing the least square of the pressure field in the silent region are given by the solution of a quadratic optimization problem. The developed method computes a robust local noise control in the presence of randomly varying parameters such as variations in the acoustic space. Numerical examples consider the noise experienced by a vehicle driver with a varying posture. In a mod…