Tsunamis - a model of their generation and propagation

A general solution is presented for water waves generated by an arbitrary movement of the bed (in space and time) in a two-dimensional fluid domain with a uniform depth. The integral solution which is developed is based on a linearized approximation to the complete (nonlinear) set of governing equations. The general solution is evaluated for the specific case of a uniform upthrust or downthrow of a block section of the bed; two time-displacement histories of the bed movement are considered.

An integral solution (based on a linear theory) is also developed for a three-dimensional fluid domain of uniform depth for a class of bed movements which are axially symmetric. The integral solution is evaluated for the specific case of a block upthrust or downthrow of a section of the bed, circular in planform, with a time-displacement history identical to one of the motions used in the two-dimensional model.

Since the linear solutions are developed from a linearized approximation of the complete nonlinear description of wave behavior, the applicability of these solutions is investigated. Two types of non-linear effects are found which limit the applicability of the linear theory: (1) large nonlinear effects which occur in the region of generation during the bed movement, and (2) the gradual growth of nonlinear effects during wave propagation.

A model of wave behavior, which includes, in an approximate manner, both linear and nonlinear effects is presented for computing wave profiles after the linear theory has become invalid due to the growth of nonlinearities during wave propagation.

An experimental program has been conducted to confirm both the linear model for the two-dimensional fluid domain and the strategy suggested for determining wave profiles during propagation after the linear theory becomes invalid. The effect of a more general time-displacement history of the moving bed than those employed in the theoretical models is also investigated experimentally.

The linear theory is found to accurately approximate the wave behavior in the region of generation whenever the total displacement of the bed is much less than the water depth. Curves are developed and confirmed by the experiments which predict gross features of the lead wave propagating from the region of generation once the values of certain nondimensional parameters (which characterize the generation process) are known. For example, the maximum amplitude of the lead wave propagating from the region of generation has been found to never exceed approximately one-half of the total bed displacement. The gross features of the tsunami resulting from the Alaskan earthquake of 27 March 1964 can be estimated from the results of this study.

Seafloor-riser interaction model

Fatigue stresses associated with extreme storms, vessel movements, and vortex-induced vibrations are critical to the performance of steel catenary risers. The critical location for fatigue damage often occurs within the touchdown zone, where cyclic interaction of the riser with the seabed occurs. Developing a model for seabed stiffness requires characterization of a number of complex nonlinear processes including trench formation, nonlinear soil stiffness, soil suction, and breakaway of the riser from the seafloor. The analytical framework utilized in this research considers the riser-seafloor interaction problem in terms of a pipe resting on a bed of springs, the stiffness characteristics of which are described by nonlinear load-deflection (P-y) curves. The P-y model allows for first penetration and uplift, as well as repenetration and small range motions within the bounding loop defined by extreme loading. The backbone curve is constructed from knowledge of the soil strength, the rate of strength increase with depth, trench width, and two additional parameters, while three parameters are necessary for the cyclic response. © ASCE 2009.

Model categories for orthogonal calculus

We restate the notion of orthogonal calculus in terms of model categories. This provides a cleaner set of results and makes the role of O(n)-equivariance clearer. Thus we develop model structures for the category of n-polynomial and n-homogeneous functors, along with Quillen pairs relating them. We then classify n-homogeneous functors, via a zig-zag of Quillen equivalences, in terms of spectra with an O(n)-action. This improves upon the classification theorem of Weiss. As an application, we develop a variant of orthogonal calculus by replacing topological spaces with orthogonal spectra.

An heuristic approach for optimal location of gas supply units in transportation system

The best places to locate the Gas Supply Units (GSUs) on a natural gas systems and their optimal allocation to loads are the key factors to organize an efficient upstream gas infrastructure. The number of GSUs and their optimal location in a gas network is a decision problem that can be formulated as a linear programming problem. Our emphasis is on the formulation and use of a suitable location model, reflecting real-world operations and constraints of a natural gas system. This paper presents a heuristic model, based on lagrangean approach, developed for finding the optimal GSUs location on a natural gas network, minimizing expenses and maximizing throughput and security of supply.The location model is applied to the Iberian high pressure natural gas network, a system modelised with 65 demand nodes. These nodes are linked by physical and virtual pipelines – road trucks with gas in liquefied form. The location model result shows the best places to locate, with the optimal demand allocation and the most economical gas transport mode: by pipeline or by road truck.

Natural gas system operation: lagrangean optimization techniques

To comply with natural gas demand growth patterns and Europe´s import dependency, the gas industry needs to organize an efficient upstream infrastructure. The best location of Gas Supply Units – GSUs and the alternative transportation mode – by phisical or virtual pipelines, are the key of a successful industry. In this work we study the optimal location of GSUs, as well as determining the most efficient allocation from gas loads to sources, selecting the best transportation mode, observing specific technical restrictions and minimizing system total costs. For the location of GSUs on system we use the P-median problem, for assigning gas demands nodes to source facilities we use the classical transportation problem. The developed model is an optimisation-based approach, based on a Lagrangean heuristic, using Lagrangean relaxation for P-median problems – Simple Lagrangean Heuristic. The solution of this heuristic can be improved by adding a local search procedure - the Lagrangean Reallocation Heuristic. These two heuristics, Simple Lagrangean and Lagrangean Reallocation, were tested on a realistic network - the primary Iberian natural gas network, organized with 65 nodes, connected by physical and virtual pipelines. Computational results are presented for both approaches, showing the location gas sources and allocation loads arrangement, system total costs and gas transportation mode.

Model predictions of the occurrence of non-Maxwellian plasmas, and analysis of their effects on EISCAT data

The recent identification of non-thermal plasmas using EISCAT data has been made possible by their occurrence during large, short-lived flow bursts. For steady, yet rapid, ion convection the only available signature is the shape of the spectrum, which is unreliable because it is open to distortion by noise and sampling uncertainty and can be mimicked by other phenomena. Nevertheless, spectral shape does give an indication of the presence of non-thermal plasma, and the characteristic shape has been observed for long periods (of the order of an hour or more) in some experiments. To evaluate this type of event properly one needs to compare it to what would be expected theoretically. Predictions have been made using the coupled thermosphere-ionosphere model developed at University College London and the University of Sheffield to show where and when non-Maxwellian plasmas would be expected in the auroral zone. Geometrical and other factors then govern whether these are detectable by radar. The results are applicable to any incoherent scatter radar in this area, but the work presented here concentrates on predictions with regard to experiments on the EISCAT facility.

How well can we model stream phosphorus concentrations in agricultural catchments?

Mechanistic catchment-scale phosphorus models appear to perform poorly where diffuse sources dominate. We investigate the reasons for this for one model, INCA-P, testing model output against 18 months of daily data in a small Scottish catchment. We examine key model processes and provide recommendations for model improvement and simplification. Improvements to the particulate phosphorus simulation are especially needed. The model evaluation procedure is then generalised to provide a checklist for identifying why model performance may be poor or unreliable, incorporating calibration, data, structural and conceptual challenges. There needs to be greater recognition that current models struggle to produce positive Nash–Sutcliffe statistics in agricultural catchments when evaluated against daily data. Phosphorus modelling is difficult, but models are not as useless as this might suggest. We found a combination of correlation coefficients, bias, a comparison of distributions and a visual assessment of time series a better means of identifying realistic simulations.

Abordagens complementares para problemas de p-medianas

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A kinetic model of phosphorus metabolism in growing goats

The effect of increasing phosphorus (P) intake on P utilization was investigated in balance experiments using 12 Saanen goats, 4 to 5 mo of age and weighing 20 to 30 kg. The goats were given similar diets with various concentrations of P, and 32P was injected to trace the movement of P in the body. A P metabolism model with four pools was developed to compute P exchanges in the system. The results showed that P absorption, bone resorption, and excretion of urinary P and endogenous and fecal P all play a part in the homeostatic control of P. Endogenous fecal output was positively correlated to P intake (P < .01). Bone resorption of P was not influenced by intake of P, and P recycling from tissues to the blood pool was lesser for low P intake. Endogenous P loss occurred even in animals fed an inadequate P diet, resulting in a negative P balance. The extrapolated minimum endogenous loss in feces was .067 g of P/d. The minimum P intake for maintenance in Saanen goats was calculated to be .61 g of P/ d or .055 g of P/(kg.75·d) at 25 kg BW. Model outputs indicate greater P flow from the blood pool to the gut and vice versa as P intake increased. Intake of P did not significantly affect P flow from bone and soft tissue to blood. The kinetic model and regressions could be used to estimate P requirement and the fate of P in goats and could also be extrapolated to both sheep and cattle.

Einfluss der Modulation von P-Glykoprotein (P-gp) auf die Verfügbarkeit und Wirksamkeit der pharmakologischen Behandlung psychiatrischer Erkrankungen im Tiermodell

P-Glykoprotein (P-gp) ist ein ATP-verbrauchender Transporter, der in Organschranken exprimiert wird, um Fremdstoffe auszuschleusen, darunter auch Psychopharmaka. Im Rahmen dieser Arbeit wurde im Tiermodell der Maus untersucht, welche pharmakokinetischen und pharmakodynamischen Konsequenzen sich bei Verabreichung von Risperidon als P-gp Modellsubstrat ergeben, wenn die Expression von P-gp induziert wird. Als potenzielle Induktoren wurden Dexamethason, Rifampicin, Quercetin, 5-Pregnen-3ß-ol-20-on-16α-Carbonitril (PCN) und Acitretin geprüft. Es konnte gezeigt werden, dass alle Substanzen die Verteilung von Risperidon und seinem aktiven Metaboliten 9-Hydroxyrisperidon beeinflussten. Während sich für Quercetin und Acitretin leichte P-gp inhibitorische Eigenschaften ergaben, die an Hand von erhöhten Konzentrationen von Risperidon und 9-Hydroxyrisperidon gezeigt werden konnten, führten die bekannten P-gp Induktoren Rifampicin, Dexamethason und PCN zu verringerten Konzentrationen im Vergleich zur Kontrollgruppe. Durch Western Blot Untersuchungen wurde bestätigt, dass die Induktoren die P-gp Expression im Hirngewebe tendenziell steigerten. Dies sprach dafür, dass bei Verabreichung einer Komedikation, die P-gp induziert, mit einer veränderten Verteilung von P-gp Substraten zu rechnen ist. Darüber hinaus konnte nachgewiesen werden, dass durch eine Hemmung bzw. Induktion von P-gp nicht nur die Pharmakokinetik, sondern auch die Pharmakodynamik von Risperidon und 9-Hydroxyrisperidon verändert wird. Dies wurde durch verhaltenspharmakologische Untersuchungen gezeigt. Durch Risperidon induzierte motorische Effekte auf dem RotaRod waren nach Induktion von P-gp abgeschwächt. Dies zeigte sich auch für Haloperidol, welches kein Substrat ist. Da P-gp abhängige Effekte in diesem Fall keine bedeutende Rolle spielen, ist davon auszugehen, dass neben der Induktion von P-gp an der Blut-Hirn Schranke auch andere Mechanismen wie z.B. eine Induktion von Enzymen der CYP-Familie an den beobachteten Effekten beteiligt sind. Bei Untersuchungen von kognitiven Leistungen in der Barnes Maze konnte gezeigt werden, dass Haloperidol im Gegensatz zu Risperidon das Lernverhalten negativ beeinflussen kann. Eine P-gp Induktion schien jedoch keinen deutlichen Einfluss auf das Lernverhalten unter Antipsychotika-Gabe zu haben und sprach vielmehr für substanzabhängige Effekte der einzelnen Antipsychotika bzw. P-gp Modulatoren. Zusatzuntersuchungen zur Hirngängigkeit von Acitretin, einem synthetischen Retinoid, welches derzeit als potenzielles Antidementivum geprüft wird, konnten belegen, dass es die Blut-Hirn Schranke überwindet. Bereits 1h nach Injektion war Acitretin in hoher Konzentration im Gehirn nachweisbar. Durch die Analyse zur Verteilung von Acitretin in Hirngewebe und Serum von P-gp Wildtyp und P-gp doppel knockout Mäusen konnte belegt werden, dass Acitretin nicht P-gp abhängig transportiert wird. Die Daten insgesamt betrachtet, lassen den Schluss zu, dass durch Verabreichung von Medikamenten, die P-gp Modulatoren sind, bei Antipsychotika mit pharmakokinetischen Interaktionen zu rechnen ist, welche die Wirksamkeit der Medikamente einschränken können.

Métodos acotados de aplicación de reglas de evolución para la implementación de Sistemas P de Transición

La característica fundamental de la Computación Natural se basa en el empleo de conceptos, principios y mecanismos del funcionamiento de la Naturaleza. La Computación Natural -y dentro de ésta, la Computación de Membranas- surge como una posible alternativa a la computación clásica y como resultado de la búsqueda de nuevos modelos de computación que puedan superar las limitaciones presentes en los modelos convencionales. En concreto, la Computación de Membranas se originó como un intento de formular un nuevo modelo computacional inspirado en la estructura y el funcionamiento de las células biológicas: los sistemas basados en este modelo constan de una estructura de membranas que actúan a la vez como separadores y como canales de comunicación, y dentro de esa estructura se alojan multiconjuntos de objetos que evolucionan de acuerdo a unas determinadas reglas de evolución. Al conjunto de dispositivos contemplados por la Computación de Membranas se les denomina genéricamente como Sistemas P. Hasta el momento los Sistemas P sólo han sido estudiados a nivel teórico y no han sido plenamente implementados ni en medios electrónicos, ni en medios bioquímicos, sólo han sido simulados o parcialmente implementados. Por tanto, la implantación de estos sistemas es un reto de investigación abierto. Esta tesis aborda uno de los problemas que debe ser resuelto para conseguir la implantación de los Sistemas P sobre plataformas hardware. El problema concreto se centra en el modelo de los Sistemas P de Transición y surge de la necesidad de disponer de algoritmos de aplicación de reglas que, independientemente de la plataforma hardware sobre la que se implementen, cumplan los requisitos de ser no deterministas, masivamente paralelos y además su tiempo de ejecución esté estáticamente acotado. Como resultado se ha obtenido un conjunto de algoritmos (tanto para plataformas secuenciales, como para plataformas paralelas) que se adecúan a las diferentes configuraciones de los Sistemas P. ABSTRACT The main feature of Natural Computing is the use of concepts, principles and mechanisms inspired by Nature. Natural Computing and within it, Membrane Computing emerges as an potential alternative to conventional computing and as from the search for new models of computation that may overcome the existing limitations in conventional models. Specifically, Membrane Computing was created to formulate a new computational paradigm inspired by the structure and functioning of biological cells: it consists of a membrane structure, which acts as separators as well as communication channels, and within this structure are stored multisets of objects that evolve according to certain evolution rules. The set of computing devices addressed by Membrane Computing are generically known P systems. Up to now, no P systems have been fully implemented yet in electronic or biochemical means. They only have been studied in theory, simulated or partially implemented. Therefore, the implementation of these systems is an open research challenge. This thesis addresses one of the problems to be solved in order to deploy P systems on hardware platforms. This specific problem is focused on the Transition P System model and emerges from the need of providing application rules algorithms that independently on the hardware platform on which they are implemented, meets the requirements of being nondeterministic, massively parallel and runtime-bounded. As a result, this thesis has developed a set of algorithms for both platforms, sequential and parallel, adapted to all possible configurations of P systems.

Strong Consistency of the Conditional Least Squares Estimator for a Nonstationary Process. Example of the Garch Model

2000 Mathematics Subject Classification: Primary: 62M10, 62J02, 62F12, 62M05, 62P05, 62P10; secondary: 60G46, 60F15.

Model updating using genetic algorithms with sequential niche technique

The authors would like to express their gratitude to organizations and people that supported this research. Piotr Omenzetter’s work within the Lloyd’s Register Foundation Centre for Safety and Reliability Engineering at the University of Aberdeen is supported by Lloyd’s Register Foundation. The Foundation helps to protect life and property by supporting engineering-related education, public engagement and the application of research. Ben Ryder of Aurecon and Graeme Cummings of HEB Construction assisted in obtaining access to the bridge and information for modelling. Luke Williams and Graham Bougen, undergraduate research students, assisted with testing.

Model updating of a full-scale bridge structure using particle swarm optimization

Peer reviewed