LaCàN - Laboratori de Càlcul Numèric
http://hdl.handle.net/2117/2072
Sun, 29 May 2016 02:18:20 GMT2016-05-29T02:18:20ZLaCàN - Laboratori de Càlcul Numèrichttp://upcommons.upc.edu/bitstream/id/22442/logoLaCaN128.gif
http://hdl.handle.net/2117/2072
Spectral analysis of the Beznar dam accelerogram. Comparison with results in the northeastern of Spain
http://hdl.handle.net/2117/87083
Spectral analysis of the Beznar dam accelerogram. Comparison with results in the northeastern of Spain
Canas Torres, José Antonio; Barbat Barbat, Horia Alejandro; Pujades Beneit, Lluís; Egozcue Rubí, Juan José; Sarrate Ramos, Josep
Numerical analysis applied to the Beznar dam accelerograms (longitudinal and transversal components) shows that the Fourier amplitude spectrum is a good approximation to the pseudo-velocity spectra. Thís fact makes possible to use displacement ar velocíty seismograrns to generate acceleration of the ground in places
where acceleration data are not available.
Maximae accelerations determined near and at the Beznar dam are compared wíth the predicted acceleration values far the regían compressed by the Pyrenees Mountains, the lberic System and the Catalonia Coastal Mountains usíng a established theoretical formula
far this region. Theoretical and observad values are consistent among them.
Tue, 17 May 2016 07:55:41 GMThttp://hdl.handle.net/2117/870832016-05-17T07:55:41ZCanas Torres, José AntonioBarbat Barbat, Horia AlejandroPujades Beneit, LluísEgozcue Rubí, Juan JoséSarrate Ramos, JosepNumerical analysis applied to the Beznar dam accelerograms (longitudinal and transversal components) shows that the Fourier amplitude spectrum is a good approximation to the pseudo-velocity spectra. Thís fact makes possible to use displacement ar velocíty seismograrns to generate acceleration of the ground in places
where acceleration data are not available.
Maximae accelerations determined near and at the Beznar dam are compared wíth the predicted acceleration values far the regían compressed by the Pyrenees Mountains, the lberic System and the Catalonia Coastal Mountains usíng a established theoretical formula
far this region. Theoretical and observad values are consistent among them.Boundary element method based on preliminary discretization
http://hdl.handle.net/2117/87019
Boundary element method based on preliminary discretization
Poblet-Puig, Jordi; Valyaev, Valery; Shanin, Andrey
A new numerical method for solving wave diffraction problems is given. The method is based on the concept of boundary elements; i.e., the unknown values are the field values on the surface of the scatterer. An analog of a boundary element method rather than a numerical approximation of the initial (continuous) problem is constructed for an approximate statement of the problem on the discrete lattice. Although it reduces the accuracy of the method, it helps to simplify the implementation significantly since the Green functions of the problem are no longer singular. In order to ensure the solution to the diffraction problem is unique (i.e., to suppress fictitious resonances), a new method is constructed similarly to the CFIE approach developed for the classical boundary element method.
The final publication is available at Springer via http://dx.doi.org/10.1134/S2070048214020082
Thu, 12 May 2016 13:41:06 GMThttp://hdl.handle.net/2117/870192016-05-12T13:41:06ZPoblet-Puig, JordiValyaev, ValeryShanin, AndreyA new numerical method for solving wave diffraction problems is given. The method is based on the concept of boundary elements; i.e., the unknown values are the field values on the surface of the scatterer. An analog of a boundary element method rather than a numerical approximation of the initial (continuous) problem is constructed for an approximate statement of the problem on the discrete lattice. Although it reduces the accuracy of the method, it helps to simplify the implementation significantly since the Green functions of the problem are no longer singular. In order to ensure the solution to the diffraction problem is unique (i.e., to suppress fictitious resonances), a new method is constructed similarly to the CFIE approach developed for the classical boundary element method.An error estimator for separated representations of highly multidimensional models
http://hdl.handle.net/2117/87000
An error estimator for separated representations of highly multidimensional models
Asaad Ammaar, Mouafk; Chinesta, F.; Díez, F.; Huerta, Antonio
Thu, 12 May 2016 10:49:40 GMThttp://hdl.handle.net/2117/870002016-05-12T10:49:40ZAsaad Ammaar, MouafkChinesta, F.Díez, F.Huerta, AntonioNumerical modeling of undersea acoustics using a partition of unity method with plane waves enrichment
http://hdl.handle.net/2117/86965
Numerical modeling of undersea acoustics using a partition of unity method with plane waves enrichment
Hospital Bravo, Raúl; Sarrate Ramos, Josep; Díez, Pedro
A new 2D numerical model to predict the underwater acoustic propagation is obtained by exploring the potential of the Partition of Unity Method (PUM) enriched with plane waves. The aim of the work is to obtain sound pressure level distributions when multiple operational noise sources are present, in order to assess the acoustic impact over the marine fauna. The model takes advantage of the suitability of the PUM for solving the Helmholtz equation, especially for the practical case of large domains and medium frequencies. The seawater acoustic absorption and the acoustic reflectance of the sea surface and sea bottom are explicitly considered, and perfectly matched layers (PML) are placed at the lateral artificial boundaries to avoid spurious reflexions. The model includes semi-analytical integration rules which are adapted to highly oscillatory integrands with the aim of reducing the computational cost of the integration step. In addition, we develop a novel strategy to mitigate the ill-conditioning of the elemental and global system matrices. Specifically, we compute a low-rank approximation of the local space of solutions, which in turn reduces the number of degrees of freedom, the CPU time and the memory footprint. Numerical examples are presented to illustrate the capabilities of the model and to assess its accuracy.
The final publication is available at Springer via http://dx.doi.org/10.1007/s00466-015-1257-8
Wed, 11 May 2016 17:52:15 GMThttp://hdl.handle.net/2117/869652016-05-11T17:52:15ZHospital Bravo, RaúlSarrate Ramos, JosepDíez, PedroA new 2D numerical model to predict the underwater acoustic propagation is obtained by exploring the potential of the Partition of Unity Method (PUM) enriched with plane waves. The aim of the work is to obtain sound pressure level distributions when multiple operational noise sources are present, in order to assess the acoustic impact over the marine fauna. The model takes advantage of the suitability of the PUM for solving the Helmholtz equation, especially for the practical case of large domains and medium frequencies. The seawater acoustic absorption and the acoustic reflectance of the sea surface and sea bottom are explicitly considered, and perfectly matched layers (PML) are placed at the lateral artificial boundaries to avoid spurious reflexions. The model includes semi-analytical integration rules which are adapted to highly oscillatory integrands with the aim of reducing the computational cost of the integration step. In addition, we develop a novel strategy to mitigate the ill-conditioning of the elemental and global system matrices. Specifically, we compute a low-rank approximation of the local space of solutions, which in turn reduces the number of degrees of freedom, the CPU time and the memory footprint. Numerical examples are presented to illustrate the capabilities of the model and to assess its accuracy.Sliding contact conditions using the master-slave approach with application on geometrically non-linear beams
http://hdl.handle.net/2117/86751
Sliding contact conditions using the master-slave approach with application on geometrically non-linear beams
Muñoz Romero, José; Jelenic, Gordan
Frictionless sliding conditions between two bodies are usually defined using either the method of Lagrangian multipliers or by prescribing an artificial (penalty) stiffness which resists the penetration at the contact point. Both of these methods impose the condition that the contact force should be normal to the contact surface, with the Lagrangian multiplier or the penalty parameter serving as a measure of this force. In this work, an alternative approach is undertaken: the frictionless sliding condition is defined through a relationship between nodal parameters of the virtual displacements of a discretised principle of virtual work. This method, which does not involve additional force parameters or degrees of freedom, is known as the master–slave or the minimum-set method and is particularly convenient for displacement-based finite-element implementation. The method is analysed in detail in context of bilateral sliding constraints characteristic of prismatic and cylindrical joints in flexible beam assemblies undergoing large overall motion. Two numerical examples are presented and assessed against the results in the literature.
Mon, 09 May 2016 09:26:06 GMThttp://hdl.handle.net/2117/867512016-05-09T09:26:06ZMuñoz Romero, JoséJelenic, GordanFrictionless sliding conditions between two bodies are usually defined using either the method of Lagrangian multipliers or by prescribing an artificial (penalty) stiffness which resists the penetration at the contact point. Both of these methods impose the condition that the contact force should be normal to the contact surface, with the Lagrangian multiplier or the penalty parameter serving as a measure of this force. In this work, an alternative approach is undertaken: the frictionless sliding condition is defined through a relationship between nodal parameters of the virtual displacements of a discretised principle of virtual work. This method, which does not involve additional force parameters or degrees of freedom, is known as the master–slave or the minimum-set method and is particularly convenient for displacement-based finite-element implementation. The method is analysed in detail in context of bilateral sliding constraints characteristic of prismatic and cylindrical joints in flexible beam assemblies undergoing large overall motion. Two numerical examples are presented and assessed against the results in the literature.Fast r-adaptivity for multiple queries of heterogeneous stochastic material fields
http://hdl.handle.net/2117/86637
Fast r-adaptivity for multiple queries of heterogeneous stochastic material fields
Cottereau, Régis; Díez, Pedro
We present an r-adaptivity approach for boundary value problems with randomly fluctuating material parameters solved through the Monte Carlo or stochastic collocation methods. This approach tailors a specific mesh for each sample of the problem. It only requires the computation of the solution of a single deterministic problem with the same geometry and the average parameter, whose numerical cost becomes marginal for large number of samples. Starting from the mesh used to solve that deterministic problem, the nodes are moved depending on the particular sample of mechanical parameter field. The reduction in the error is small for each sample but sums up to reduce the overall bias on the statistics estimated through the Monte Carlo scheme. Several numerical examples in 2D are presented.
Thu, 05 May 2016 12:08:08 GMThttp://hdl.handle.net/2117/866372016-05-05T12:08:08ZCottereau, RégisDíez, PedroWe present an r-adaptivity approach for boundary value problems with randomly fluctuating material parameters solved through the Monte Carlo or stochastic collocation methods. This approach tailors a specific mesh for each sample of the problem. It only requires the computation of the solution of a single deterministic problem with the same geometry and the average parameter, whose numerical cost becomes marginal for large number of samples. Starting from the mesh used to solve that deterministic problem, the nodes are moved depending on the particular sample of mechanical parameter field. The reduction in the error is small for each sample but sums up to reduce the overall bias on the statistics estimated through the Monte Carlo scheme. Several numerical examples in 2D are presented.AAR-based decomposition method for lower-bound limit analysis
http://hdl.handle.net/2117/86620
AAR-based decomposition method for lower-bound limit analysis
Muñoz Romero, José; Rabiei, Syednima
Despite the recent progress in optimisation techniques, finite-element stability analysis of realistic three-dimensional problems is still hampered by the size of the resulting optimisation problem. Current solvers may take a prohibitive computational time, if they give a solution at all. The possible remedies to this are the design of adaptive de-remeshing techniques, decomposition of the system of equations or of the optimisation problem. This paper concentrates on the last approach, and presents an algorithm especially suited for limit analysis. Optimisation problems in limit analysis are in general convex but non-linear. This fact renders the design of decomposition techniques specially challenging. The efficiency of general approaches such as Benders or Dantzig–Wolfe is not always satisfactory, and strongly depends on the structure of the optimisation problem. This work presents a new method that is based on rewriting the feasibility region of the global optimisation problem as the intersection of two subsets. By resorting to the averaged alternating reflections (AAR) method in order to find the distance between the sets, the optimisation problem is successfully solved in a decomposed manner. Some representative examples illustrate the application of the method and its efficiency with respect to other well-known decomposition algorithms.
Thu, 05 May 2016 08:50:31 GMThttp://hdl.handle.net/2117/866202016-05-05T08:50:31ZMuñoz Romero, JoséRabiei, SyednimaDespite the recent progress in optimisation techniques, finite-element stability analysis of realistic three-dimensional problems is still hampered by the size of the resulting optimisation problem. Current solvers may take a prohibitive computational time, if they give a solution at all. The possible remedies to this are the design of adaptive de-remeshing techniques, decomposition of the system of equations or of the optimisation problem. This paper concentrates on the last approach, and presents an algorithm especially suited for limit analysis. Optimisation problems in limit analysis are in general convex but non-linear. This fact renders the design of decomposition techniques specially challenging. The efficiency of general approaches such as Benders or Dantzig–Wolfe is not always satisfactory, and strongly depends on the structure of the optimisation problem. This work presents a new method that is based on rewriting the feasibility region of the global optimisation problem as the intersection of two subsets. By resorting to the averaged alternating reflections (AAR) method in order to find the distance between the sets, the optimisation problem is successfully solved in a decomposed manner. Some representative examples illustrate the application of the method and its efficiency with respect to other well-known decomposition algorithms.Constructive and destructive interplay between piezoelectricity and flexoelectricity in flexural sensors and actuators
http://hdl.handle.net/2117/86474
Constructive and destructive interplay between piezoelectricity and flexoelectricity in flexural sensors and actuators
Abdollahi Hosnijeh, Amir; Arias Vicente, Irene
Flexoelectricity is an electromechanical effect coupling polarization to strain gradients. It fundamentally differs from piezoelectricity because of its size-dependence and symmetry. Flexoelectricity is generally perceived as a small effect noticeable only at the nanoscale. Since ferroelectric ceramics have a particularly high flexoelectric coefficient, however, it may play a significant role as piezoelectric transducers shrink to the submicrometer scale. We examine this issue with a continuum model self-consistently treating piezo- and flexoelectricity. We show that in piezoelectric device configurations that induce strain gradients and at small but technologically relevant scales, the electromechanical coupling may be dominated by flexoelectricity. More importantly, depending on the device design flexoelectricity may enhance or reduce the effective piezoelectric effect. Focusing on bimorph configurations, we show that configurations that are equivalent at large scales exhibit dramatically different behavior for thicknesses below 100¿nm for typical piezoelectric materials. Our results suggest flexoelectric-aware designs for small-scale piezoelectric bimorph transducers.
Mon, 02 May 2016 12:12:55 GMThttp://hdl.handle.net/2117/864742016-05-02T12:12:55ZAbdollahi Hosnijeh, AmirArias Vicente, IreneFlexoelectricity is an electromechanical effect coupling polarization to strain gradients. It fundamentally differs from piezoelectricity because of its size-dependence and symmetry. Flexoelectricity is generally perceived as a small effect noticeable only at the nanoscale. Since ferroelectric ceramics have a particularly high flexoelectric coefficient, however, it may play a significant role as piezoelectric transducers shrink to the submicrometer scale. We examine this issue with a continuum model self-consistently treating piezo- and flexoelectricity. We show that in piezoelectric device configurations that induce strain gradients and at small but technologically relevant scales, the electromechanical coupling may be dominated by flexoelectricity. More importantly, depending on the device design flexoelectricity may enhance or reduce the effective piezoelectric effect. Focusing on bimorph configurations, we show that configurations that are equivalent at large scales exhibit dramatically different behavior for thicknesses below 100¿nm for typical piezoelectric materials. Our results suggest flexoelectric-aware designs for small-scale piezoelectric bimorph transducers.Cell-centred model for the simulation of curved cellular monolayers
http://hdl.handle.net/2117/86189
Cell-centred model for the simulation of curved cellular monolayers
Mosafa, Payman; Asadipour, Nina; Millán, Raúl Daniel; Rodríguez Ferran, Antonio; Muñoz Romero, José
This paper presents a cell-centred model for the simulation of planar and curved multicellular soft tissues. We propose a computational model that includes stress relaxation due to cell reorganisation (intercellular connectivity changes) and cytoskeleton remodelling (intracellular changes). Cells are represented by their cell centres, and their mechanical interaction is modelled through active non-linear elastic laws with a dynamically changing resting length. Special attention is paid to the handling of connectivity changes between cells, and the relaxation that the tissues exhibit under these topological changes. Cell--cell connectivity is computed by resorting to a Delaunay triangulation, which is combined with a mapping technique in order to obtain triangulations on curved manifolds. Our numerical results show that even a linear elastic cell--cell interaction model may induce a global non-linear response due to the reorganisation of the cell connectivity. This plastic-like behaviour is combined with a non-linear rheological law where the resting length depends on the elastic strain, mimicking the global visco-elastic response of tissues. The model is applied to simulate the elongation of planar and curved monolayers.
Tue, 26 Apr 2016 10:28:30 GMThttp://hdl.handle.net/2117/861892016-04-26T10:28:30ZMosafa, PaymanAsadipour, NinaMillán, Raúl DanielRodríguez Ferran, AntonioMuñoz Romero, JoséThis paper presents a cell-centred model for the simulation of planar and curved multicellular soft tissues. We propose a computational model that includes stress relaxation due to cell reorganisation (intercellular connectivity changes) and cytoskeleton remodelling (intracellular changes). Cells are represented by their cell centres, and their mechanical interaction is modelled through active non-linear elastic laws with a dynamically changing resting length. Special attention is paid to the handling of connectivity changes between cells, and the relaxation that the tissues exhibit under these topological changes. Cell--cell connectivity is computed by resorting to a Delaunay triangulation, which is combined with a mapping technique in order to obtain triangulations on curved manifolds. Our numerical results show that even a linear elastic cell--cell interaction model may induce a global non-linear response due to the reorganisation of the cell connectivity. This plastic-like behaviour is combined with a non-linear rheological law where the resting length depends on the elastic strain, mimicking the global visco-elastic response of tissues. The model is applied to simulate the elongation of planar and curved monolayers.Delay assignment optimization strategies at pre-tactical and tactical levels
http://hdl.handle.net/2117/85616
Delay assignment optimization strategies at pre-tactical and tactical levels
Villardi de Montlaur, Adeline de; Delgado Muñoz, Luis
—This paper compares different optimization strategies
for the minimization of flight and passenger delays at
two levels: pre-tactical, with on-ground delay at origin, and
tactical, with airborne delay close to the destination airport.
The optimization model is based on the ground holding problem
and uses various cost functions. The scenario considered takes
place in a busy European airport and includes realistic values of
traffic. Uncertainty is introduced in the model for the passenger
allocation, minimum time required for turnaround and tactical
uncertainty. Performance of the various optimization processes is
presented and compared to ratio by schedule results.
Wed, 13 Apr 2016 13:05:20 GMThttp://hdl.handle.net/2117/856162016-04-13T13:05:20ZVillardi de Montlaur, Adeline deDelgado Muñoz, Luis—This paper compares different optimization strategies
for the minimization of flight and passenger delays at
two levels: pre-tactical, with on-ground delay at origin, and
tactical, with airborne delay close to the destination airport.
The optimization model is based on the ground holding problem
and uses various cost functions. The scenario considered takes
place in a busy European airport and includes realistic values of
traffic. Uncertainty is introduced in the model for the passenger
allocation, minimum time required for turnaround and tactical
uncertainty. Performance of the various optimization processes is
presented and compared to ratio by schedule results.