Efficient computation of the pressures developed during high-speed train passing events

The paper resumes the results obtained applying various implementations of the direct boundary element method (BEM) to the solution of the Laplace Equation governing the potential flow problem during everyday service manoeuvres of high-speed trains. In particular the results of train passing events at three different speed combinations are presented. Some recommendations are given in order to reduce calculation times which as is demonstrated can be cut down to not exceed reasonable limits even when using nowadays office PCs. Thus the method is shown to be a very valuable tool for the design engineer.

Elementos Finitos no Lineales en Mecánica de Sólidos y Estructuras. Introducción básica y formulación del Método de los Elementos Finitos : Curso

Con el método de elementos finitos se ha desarrollado un procedimiento para el tratamiento de ecuaciones diferenciales y su resolución numérica. En este capítulo se va a realizar la aplicación a problemas de potencial regidos por la ecuación de Laplace. Se trata de una situación relativamente sencilla pero que permite comentar la mayoría de las características del método de elementos finitos en desplazamientos y sirve como introducción a situaciones más complicadas. Además, esta ecuación rige problemas de interés como la filtración, la torsión de cilindros, la distribución de termperaturas en régimen estacionario, etc. En lo que sigue la exposición se limita a problemas planos. Tras un primer ejemplo representativo del tema se realiza la formulación débil del problema de potencial. Acto seguido se introduce el concepto de discretización con toda su parafernalia de elementos, nudos, variables nodales, funciones de interpolación etc. El capítulo finaliza con la mención de algunos temas más especializados que el lector estudioso deberá proseguir en otro tipo de textos.

Improvements for B.E.M. implementation in micros

This paper presents a computer program developed to run in a micro I.B.M.-P.C. wich incorporates some features in order to optimize the number of operations needed to compute the solution of plane potential problems governed by Laplace's equation by using the Boundary Integral Equation Method (B.I.E.M.). Also incorporated is a routine to plot isolines inside the domain under study.

Elementos de contorno p-adaptables en problemas tridimensionales

Se presentan en este artículo los resultados obtenidos con la aplicación de la filosofía de p-convergencia a problemas tridimensionales controlados por la ecuación de Laplace. Se pretende con ello sintetizar las ventajas inherentes a la discretización de contorno y a la jerarquización de las funciones de interpolación.

The B.I.E.M. applied to flow through porous media

This paper refers to the numerical solution of the classical Darcy's problem of plane fluid through isotropic media. Regarding the numerical procedure,the Laplace equation, is a classical one in mathematical physics and several procedures have been devised in order to solve it. So as to show the capability of the method, the paper presents some exemples.

Guía de introducción al Método de los Elementos de Contorno

Entre la impresionante floración de procedimientos de cálculo, provocada por la aplicación intensiva del ordenador, el llamado Método de los Elementos de Contorno (Boundary Element Method o Boundary Integral Equation Method) parece afianzarse como una alternativa útil al omnipresente Método de los Elementos Finitos que ya ha sido incorporado, como una herramienta de trabajo más, al cotidiano quehacer de la ingeniería. En España, tras unos intentos precursores que se señalan en el texto, la actividad más acusada en su desarrollo y mejora se ha centrado alrededor del Departamento que dirige uno de los autores. Después de la tesis doctoral de J. Domínguez en 1977 que introdujo en España la técnica del llamado "método directo", se han producido numerosas aportaciones en forma de artículos o tesis de investigación que han permitido alcanzar un nivel de conocimientos notable. En esta obrita se pretende transmitir parte de la experiencia adquirida, siquiera sea a nivel elemental y en un campo limitado de aplicación. La filosofía es semejante a la del pequeño libro de Hinton y Owen "A simple guide to finite elements" (Pineridge Press, 1980) que tanta aceptación ha tenido entre los principiantes. El libro se articula alrededor de un sólo tema, la solución del problema de Laplace, y se limitan los desarrollos matemáticos al mínimo imprescindible para el fácil seguimiento de áquel. Tras unos capítulos iniciales de motivación y centrado se desarrolla la técnica para problemas planos, tridimensionales y axisimétricos, limitando los razonamientos a los elementos más sencillos de variación constante o lineal. Finalmente, se incluye un capítulo descriptivo donde se avizoran temas que pueden provocar un futuro interés del estudioso. Para completar la información se ha añadido un apéndice en el que se recoge un pequeño programa para microordenador, con el objetivo de que se contemple la sencillez de programación para el caso plano. El programa es mejorable en muchos aspectos pero creemos que, con ello, mantiene un nivel de legibilidad adecuado para que el lector ensaye sobre él las modificaciones que se indican en los ejercicios al final del capítulo y justamente la provocación de ese aprendizaje es nuestro objetivo final.

Discussion of “Initial Pore Pressure from Vertical Surface Loads” by Jacobo Bielak (September, 1982)

The author presents a very interesting application of the ideas developed by Scott to determine the initial pore pressure in excess of the hydrostatic pore pressure in linear, elastic, homogeneous and isotropic soil-skeleton. Scott demonstrates that under vertical surface loads the problem is governed by Laplace's equation. Nevertheless the writers' think that it could be interesting to state clearly the conditions under which this analogy can be applied.

P-adaptive boundary elements for three-dimensional potential problems

This paper deals with the boundary element method (BEM) p-convergence approach applied to three-dimensional problems governed by Laplace's equation. The advantages derived from the boundary discretization and hierarchical interpolation functions are collated in order to minimize human effort in preparation of input data and improve numerical results.

One-dimensional linear analysis of the liquid injection or removal in a liquid bridge

The filling-withdrawal process of a long liquid bridge is analyzed using a one-dimensional linearized model for the dynamics of the liquid column. To carry out this study, a well-known standard operational method (Laplace transform) has been used, and time variation of both liquid velocity field and interface shape are obtained.

A time-domain finite element method for seakeeping and wave resistance problems

El objetivo de la tesis es la investigación de algoritmos numéricos para el desarrollo de herramientas numéricas para la simulación de problemas tanto de comportamiento en la mar como de resistencia al avance de buques y estructuras flotantes. La primera herramienta desarrollada resuelve el problema de difracción y radiación de olas. Se basan en el método de los elementos finitos (MEF) para la resolución de la ecuación de Laplace, así como en esquemas basados en MEF, integración a lo largo de líneas de corriente, y en diferencias finitas desarrollados para la condición de superficie libre. Se han desarrollado herramientas numéricas para la resolución de la dinámica de sólido rígido en sistemas multicuerpos con ligaduras. Estas herramientas han sido integradas junto con la herramienta de resolución de olas difractadas y radiadas para la resolución de problemas de interacción de cuerpos con olas. También se han diseñado algoritmos de acoplamientos con otras herramientas numéricas para la resolución de problemas multifísica. En particular, se han realizado acoplamientos con una herramienta numérica basada de cálculo de estructuras con MEF para problemas de interacción fluido-estructura, otra de cálculo de líneas de fondeo, y con una herramienta numérica de cálculo de flujos en tanques internos para problemas acoplados de comportamiento en la mar con “sloshing”. Se han realizado simulaciones numéricas para la validación y verificación de los algoritmos desarrollados, así como para el análisis de diferentes casos de estudio con aplicaciones diversas en los campos de la ingeniería naval, oceánica, y energías renovables marinas. ABSTRACT The objective of this thesis is the research on numerical algorithms to develop numerical tools to simulate seakeeping problems as well as wave resistance problems of ships and floating structures. The first tool developed is a wave diffraction-radiation solver. It is based on the finite element method (FEM) in order to solve the Laplace equation, as well as numerical schemes based on FEM, streamline integration, and finite difference method tailored for solving the free surface boundary condition. It has been developed numerical tools to solve solid body dynamics of multibody systems with body links across them. This tool has been integrated with the wave diffraction-radiation solver to solve wave-body interaction problems. Also it has been tailored coupling algorithms with other numerical tools in order to solve multi-physics problems. In particular, it has been performed coupling with a MEF structural solver to solve fluid-structure interaction problems, with a mooring solver, and with a solver capable of simulating internal flows in tanks to solve couple seakeeping-sloshing problems. Numerical simulations have been carried out to validate and verify the developed algorithms, as well as to analyze case studies in the areas of marine engineering, offshore engineering, and offshore renewable energy.

Fibroblast Growth Factor (FGF) Soluble Receptor 1 Acts as a Natural Inhibitor of FGF2 Neurotrophic Activity during Retinal Degeneration

Fibroblast growth factors (FGF) 1 and 2 and their tyrosine kinase receptor (FGFR) are present throughout the adult retina. FGFs are potential mitogens, but adult retinal cells are maintained in a nonproliferative state unless the retina is damaged. Our work aims to find a modulator of FGF signaling in normal and pathological retina. We identified and sequenced a truncated FGFR1 form from rat retina generated by the use of selective polyadenylation sites. This 70-kDa form of soluble extracellular FGFR1 (SR1) was distributed mainly localized in the inner nuclear layer of the retina, whereas the full-length FGFR1 form was detected in the retinal Muller glial cells. FGF2 and FGFR1 mRNA levels greatly increased in light-induced retinal degeneration. FGFR1 was detected in the radial fibers of activated retinal Muller glial cells. In contrast, SR1 mRNA synthesis followed a biphasic pattern of down- and up-regulation, and anti-SR1 staining was intense in retinal pigmented epithelial cells. The synthesis of SR1 and FGFR1 specifically and independently regulated in normal and degenerating retina suggests that changes in the proportion of various FGFR forms may control the bioavailability of FGFs and thus their potential as neurotrophic factors. This was demonstrated in vivo during retinal degeneration when recombinant SR1 inhibited the neurotrophic activity of exogenous FGF2 and increased damaging effects of light by inhibiting endogenous FGF. This study highlights the significance of the generation of SR1 in normal and pathological conditions.

New infinite families of exact sums of squares formulas, Jacobi elliptic functions, and Ramanujan’s tau function

In this paper, we give two infinite families of explicit exact formulas that generalize Jacobi’s (1829) 4 and 8 squares identities to 4n2 or 4n(n + 1) squares, respectively, without using cusp forms. Our 24 squares identity leads to a different formula for Ramanujan’s tau function τ(n), when n is odd. These results arise in the setting of Jacobi elliptic functions, Jacobi continued fractions, Hankel or Turánian determinants, Fourier series, Lambert series, inclusion/exclusion, Laplace expansion formula for determinants, and Schur functions. We have also obtained many additional infinite families of identities in this same setting that are analogous to the η-function identities in appendix I of Macdonald’s work [Macdonald, I. G. (1972) Invent. Math. 15, 91–143]. A special case of our methods yields a proof of the two conjectured [Kac, V. G. and Wakimoto, M. (1994) in Progress in Mathematics, eds. Brylinski, J.-L., Brylinski, R., Guillemin, V. & Kac, V. (Birkhäuser Boston, Boston, MA), Vol. 123, pp. 415–456] identities involving representing a positive integer by sums of 4n2 or 4n(n + 1) triangular numbers, respectively. Our 16 and 24 squares identities were originally obtained via multiple basic hypergeometric series, Gustafson’s Cℓ nonterminating 6φ5 summation theorem, and Andrews’ basic hypergeometric series proof of Jacobi’s 4 and 8 squares identities. We have (elsewhere) applied symmetry and Schur function techniques to this original approach to prove the existence of similar infinite families of sums of squares identities for n2 or n(n + 1) squares, respectively. Our sums of more than 8 squares identities are not the same as the formulas of Mathews (1895), Glaisher (1907), Ramanujan (1916), Mordell (1917, 1919), Hardy (1918, 1920), Kac and Wakimoto, and many others.

Mathematical analysis of activation thresholds in enzyme-catalyzed positive feedbacks: application to the feedbacks of blood coagulation.

A hierarchy of enzyme-catalyzed positive feedback loops is examined by mathematical and numerical analysis. Four systems are described, from the simplest, in which an enzyme catalyzes its own formation from an inactive precursor, to the most complex, in which two sequential feedback loops act in a cascade. In the latter we also examine the function of a long-range feedback, in which the final enzyme produced in the second loop activates the initial step in the first loop. When the enzymes generated are subject to inhibition or inactivation, all four systems exhibit threshold properties akin to excitable systems like neuron firing. For those that are amenable to mathematical analysis, expressions are derived that relate the excitation threshold to the kinetics of enzyme generation and inhibition and the initial conditions. For the most complex system, it was expedient to employ numerical simulation to demonstrate threshold behavior, and in this case long-range feedback was seen to have two distinct effects. At sufficiently high catalytic rates, this feedback is capable of exciting an otherwise subthreshold system. At lower catalytic rates, where the long-range feedback does not significantly affect the threshold, it nonetheless has a major effect in potentiating the response above the threshold. In particular, oscillatory behavior observed in simulations of sequential feedback loops is abolished when a long-range feedback is present.

Variability of the ocean heat content during the last millennium - an assessment with the ECHO-g Model

Studies addressing climate variability during the last millennium generally focus on variables with a direct influence on climate variability, like the fast thermal response to varying radiative forcing, or the large-scale changes in atmospheric dynamics (e. g. North Atlantic Oscillation). The ocean responds to these variations by slowly integrating in depth the upper heat flux changes, thus producing a delayed influence on ocean heat content (OHC) that can later impact low frequency SST (sea surface temperature) variability through reemergence processes. In this study, both the externally and internally driven variations of the OHC during the last millennium are investigated using a set of fully coupled simulations with the ECHO-G (coupled climate model ECHAMA4 and ocean model HOPE-G) atmosphere-ocean general circulation model (AOGCM). When compared to observations for the last 55 yr, the model tends to overestimate the global trends and underestimate the decadal OHC variability. Extending the analysis back to the last one thousand years, the main impact of the radiative forcing is an OHC increase at high latitudes, explained to some extent by a reduction in cloud cover and the subsequent increase of short-wave radiation at the surface. This OHC response is dominated by the effect of volcanism in the preindustrial era, and by the fast increase of GHGs during the last 150 yr. Likewise, salient impacts from internal climate variability are observed at regional scales. For instance, upper temperature in the equatorial Pacific is controlled by ENSO (El Nino Southern Oscillation) variability from interannual to multidecadal timescales. Also, both the Pacific Decadal Oscillation (PDO) and the Atlantic Multidecadal Oscillation (AMO) modulate intermittently the interdecadal OHC variability in the North Pacific and Mid Atlantic, respectively. The NAO, through its influence on North Atlantic surface heat fluxes and convection, also plays an important role on the OHC at multiple timescales, leading first to a cooling in the Labrador and Irminger seas, and later on to a North Atlantic warming, associated with a delayed impact on the AMO.

On the time to reach a critical number of infections in epidemic models with infective and susceptible immigrants.

In this paper we examine the time T to reach a critical number K0 of infections during an outbreak in an epidemic model with infective and susceptible immigrants. The underlying process X, which was first introduced by Ridler-Rowe (1967), is related to recurrent diseases and it appears to be analytically intractable. We present an approximating model inspired from the use of extreme values, and we derive formulae for the Laplace-Stieltjes transform of T and its moments, which are evaluated by using an iterative procedure. Numerical examples are presented to illustrate the effects of the contact and removal rates on the expected values of T and the threshold K0, when the initial time instant corresponds to an invasion time. We also study the exact reproduction number Rexact,0 and the population transmission number Rp, which are random versions of the basic reproduction number R0.