Fast numerical methods for mixed, singular Helmholtz boundary value problems and Laplace eigenvalue problems - with applications to antenna design, sloshing, electromagnetic scattering and spectral geometry

This thesis presents a novel class of algorithms for the solution of scattering and eigenvalue problems on general two-dimensional domains under a variety of boundary conditions, including non-smooth domains and certain "Zaremba" boundary conditions - for which Dirichlet and Neumann conditions are specified on various portions of the domain boundary. The theoretical basis of the methods for the Zaremba problems on smooth domains concern detailed information, which is put forth for the first time in this thesis, about the singularity structure of solutions of the Laplace operator under boundary conditions of Zaremba type. The new methods, which are based on use of Green functions and integral equations, incorporate a number of algorithmic innovations, including a fast and robust eigenvalue-search algorithm, use of the Fourier Continuation method for regularization of all smooth-domain Zaremba singularities, and newly derived quadrature rules which give rise to high-order convergence even around singular points for the Zaremba problem. The resulting algorithms enjoy high-order convergence, and they can tackle a variety of elliptic problems under general boundary conditions, including, for example, eigenvalue problems, scattering problems, and, in particular, eigenfunction expansion for time-domain problems in non-separable physical domains with mixed boundary conditions.

Prediction of Multi-Target Networks of Neuroprotective Compounds with Entropy Indices and Synthesis, Assay, and Theoretical Study of New Asymmetric 1,2-Rasagiline Carbamates

In a multi-target complex network, the links (L-ij) represent the interactions between the drug (d(i)) and the target (t(j)), characterized by different experimental measures (K-i, K-m, IC50, etc.) obtained in pharmacological assays under diverse boundary conditions (c(j)). In this work, we handle Shannon entropy measures for developing a model encompassing a multi-target network of neuroprotective/neurotoxic compounds reported in the CHEMBL database. The model predicts correctly >8300 experimental outcomes with Accuracy, Specificity, and Sensitivity above 80%-90% on training and external validation series. Indeed, the model can calculate different outcomes for >30 experimental measures in >400 different experimental protocolsin relation with >150 molecular and cellular targets on 11 different organisms (including human). Hereafter, we reported by the first time the synthesis, characterization, and experimental assays of a new series of chiral 1,2-rasagiline carbamate derivatives not reported in previous works. The experimental tests included: (1) assay in absence of neurotoxic agents; (2) in the presence of glutamate; and (3) in the presence of H2O2. Lastly, we used the new Assessing Links with Moving Averages (ALMA)-entropy model to predict possible outcomes for the new compounds in a high number of pharmacological tests not carried out experimentally.

Feedback on the Rate and Depth of Chest Compressions during Cardiopulmonary Resuscitation Using Only Accelerometers

Background Quality of cardiopulmonary resuscitation (CPR) is key to increase survival from cardiac arrest. Providing chest compressions with adequate rate and depth is difficult even for well-trained rescuers. The use of real-time feedback devices is intended to contribute to enhance chest compression quality. These devices are typically based on the double integration of the acceleration to obtain the chest displacement during compressions. The integration process is inherently unstable and leads to important errors unless boundary conditions are applied for each compression cycle. Commercial solutions use additional reference signals to establish these conditions, requiring additional sensors. Our aim was to study the accuracy of three methods based solely on the acceleration signal to provide feedback on the compression rate and depth. Materials and Methods We simulated a CPR scenario with several volunteers grouped in couples providing chest compressions on a resuscitation manikin. Different target rates (80, 100, 120, and 140 compressions per minute) and a target depth of at least 50 mm were indicated. The manikin was equipped with a displacement sensor. The accelerometer was placed between the rescuer's hands and the manikin's chest. We designed three alternatives to direct integration based on different principles (linear filtering, analysis of velocity, and spectral analysis of acceleration). We evaluated their accuracy by comparing the estimated depth and rate with the values obtained from the reference displacement sensor. Results The median (IQR) percent error was 5.9% (2.8-10.3), 6.3% (2.9-11.3), and 2.5% (1.2-4.4) for depth and 1.7% (0.0-2.3), 0.0% (0.0-2.0), and 0.9% (0.4-1.6) for rate, respectively. Depth accuracy depended on the target rate (p < 0.001) and on the rescuer couple (p < 0.001) within each method. Conclusions Accurate feedback on chest compression depth and rate during CPR is possible using exclusively the chest acceleration signal. The algorithm based on spectral analysis showed the best performance. Despite these encouraging results, further research should be conducted to asses the performance of these algorithms with clinical data.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Subdivision-stabilised immersed b-spline finite elements for moving boundary flows

An immersed finite element method is presented to compute flows with complex moving boundaries on a fixed Cartesian grid. The viscous, incompressible fluid flow equations are discretized with b-spline basis functions. The two-scale relation for b-splines is used to implement an elegant and efficient technique to satisfy the LBB condition. On non-grid-aligned fluid domains and at moving boundaries, the boundary conditions are enforced with a consistent penalty method as originally proposed by Nitsche. In addition, a special extrapolation technique is employed to prevent the loss of numerical stability in presence of arbitrarily small cut-cells. The versatility and accuracy of the proposed approach is demonstrated by means of convergence studies and comparisons with previous experimental and computational investigations.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

Static analysis of Euler-Bernoulli beams with multiple unilateral cracks under combined axial and transverse loads

The paper deals with the static analysis of pre-damaged Euler-Bernoulli beams with any number of unilateral cracks and subjected to tensile or compression forces combined with arbitrary transverse loads. The mathematical representation of cracks with a bilateral behaviour (i.e. always open) via Dirac delta functions is extended by introducing a convenient switching variable, which allows each crack to be open or closed depending on the sign of the axial strain at the crack centre. The proposed model leads to analytical solutions, which depend on four integration constants (to be computed by enforcing the boundary conditions) along with the Boolean switching variables associated with the cracks (whose role is to turn on and off the additional flexibility due to the presence of the cracks). An efficient computational procedure is also presented and numerically validated. For this purpose, the proposed approach is applied to two pre-damaged beams, with different damage and loading conditions, and the results so obtained are compared against those given by a standard finite element code (in which the correct opening of the cracks is pre-assigned), always showing a perfect agreement. © 2013 Elsevier Ltd. All rights reserved.

潮波模型的优化开边界方法研究

数值模式是潮波研究的一种有利手段，但在研究中会面临各种具体问题，包括开边界条件的确定、底摩擦系数和耗散系数的选取等。数据同化是解决这些问题的一种途径，即利用有限数量的潮汐观测资料对潮波进行最优估计，其根本目的是迫使模型预报值逼近观测值，使模式不要偏离实际情况太远。本文采用了一种优化开边界方法，沿着数值模型的开边界优化潮汐水位信息，目的是设法使数值解在动力约束的意义下接近观测值，获得研究区域的潮汐结果。边界值由指定优化问题的解来定，以提高模拟区域的潮汐精度，最优问题的解是基于通过开边界的能量通量的变化，处理开边界处的观测值与计算值之差的最小化。这里提供了辐射型边界条件，由Reid 和Bodine（本文简称为RB）推导，我们将采用的优化后的RB方法（称为ORB）是优化开边界的特殊情况。 本文对理想矩形海域（ E- E， N- N， 分辨率 ）进行了潮波模拟，有东部开边界，模式采用ECOM3D模式。对数据结果的误差分析采用，振幅平均偏差，平均绝对偏差，平均相对误差和均方根偏差四个值来衡量模拟结果的好坏程度。 需要优化入开边界的解析潮汐值本文采用的解析解由方国洪《海湾的潮汐与潮流》（1966年）方法提供，为验证本文所做的解析解和方文的一致，本文做了其第一个例子的关键值a,b,z,结果与其结果吻合的相当好。但略有差别，分析的可能原因是两法在具体迭代方案和计算机保留小数上有区别造成微小误差。另外，我们取m=20，得到更精确的数值，我们发现对前十项的各项参数值，取m=10,m=20各项参数略有改进。当然我们可以获得m更大的各项参数值。 同时为了检验解析解的正确性讨论m和l变化对边界值的影响，结果指出，增大m，m=20时，u的模最大在本身u1或u2的模的6%；m=100时，u的模最大在本身u1或u2的模的4%；m再增大，m=1000时，u的模最大在本身u1或u2的模的4%，改变不大。当l<1时， =0处u的模最大为2。当l=1时， =0处u的模最大为0.1，当l>1时，l越大，u的模越小，当l=10时，u的模最大为0.001，可以认为为0。 为检验该优化方法的应用情况，我们对理想矩形区域进行模拟，首先将本文所采用的优化开边界方法应用于30m的情况，在开边界优化入开边界得出模式解，所得模拟结果与解析解吻合得相当好，该模式解和解析解在整个区域上，振幅平均绝对偏差为9.9cm,相位平均绝对偏差只有4.0 ，均方根偏差只有13.3cm，说明该优化方法在潮波模型中有效。 为验证该优化方法在各种条件下的模拟结果情况，在下面我们做了三类敏感性试验： 第一类试验：为证明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，我们来比较ORB条件与RB条件的优劣，我们模拟用了两个不同的摩擦系数，k分别为：0,0.00006。 结果显示，针对不同摩擦系数，显示在开边界上使用ORB条件的解比使用RB条件的解无论是振幅还是相位都有显著改善，两个试验均方根偏差优化程度分别为84.3%，83.7%。说明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，大大提高了模拟水平。上述的两个试验得出， k=0.00006优化结果比k=0的好。 第二类试验，使用ORB条件确定优化开边界情况下，在东西边界加入出入流的情况，流考虑线性和非线性情况，结果显示，加入流的情况，潮汐模拟的效果降低不少，流为1Sv的情况要比5Sv的情况均方根偏差相差20cm,而不加流的情况只有0.2cm。线性流和非线性流情况两者模式解相差不大，振幅，相位各项指数都相近, 说明流的线性与否对结果影响不大。 第三类试验，不仅在开边界使用ORB条件，在模式内部也使用ORB条件，比较了内部优化和不优化情况与解析解的偏差。结果显示，选用不同的k,振幅都能得到很好的模拟，而相位相对较差。另外，在内部优化的情况下，考虑不同的k的模式解， 我们选用了与解析解相近的6个模式解的k,结果显示，不同的k,振幅都能得到很好的模拟，而相位较差。 总之，在开边界使用ORB条件比使用RB条件好，振幅相位都有大幅度改进，在加入出入流情况下，流的大小对模拟结果有影响，但线形流和非线性流差别不大。内部优化的结果显示，模式采用不同的k都能很好模拟解析解的振幅。

Boundary Element Method (BEM) Analysis for Galvanic Corrosion of Hot Dip Galvanized Steel Immersed in Seawater

A numerical analysis of galvanic corrosion of hot-dip galvanized steel immersed in seawater was presented. The analysis was based on the boundary element methods (BEMs) coupled with Newton-Raphson iterative technique to treat the nonlinear boundary conditions, which were determined by the experimental polarization curves. Results showed that galvanic current density concentrates on the boundary of steel substrate and zinc coating, and the sacrificial protection of zinc coating to steel substrate results in overprotection of steel cathode. Not only oxygen reduction but also hydrogen reduction could occur as cathode reactions, which probably led up to the adsorption and absorption of hydrogen atoms. Flat galvanized steel tensile sample shows a brittle behavior similar to hydrogen embrittlement according to the SSRT (show strain rate test) in seawater.

孔隙介质中油气二次运移动力学机制研究

Study of dynamical mechanism of hydrocarbon secondary migration is the key research project of China Petroleum and Chemical Corporation in the ninth "Five-Year Plan", and this research is the hot and difficult spot and frontline in the domain of reservoir forming study in recent years. It is a common recognition that the dynamical mechanism of hydrocarbon's secondary migration is the important factor to control the reservoir type, distribution and oil/gas abundance. Therefore, to study this mechanism and establish the modes of hydrocarbon's migration and accumulation in different reservoirs under different conditions are of great theoretical meaningfulness and practical value on both developing the theory and method of hydrocarbon migration/accumulation dynamics in terrestrial rift-subsidence lacustrine basins and guiding the exploration and production. A laboratory for physical simulation of hydrocarbon's secondary migration/accumulation mechanism has been build up. 12 types of physical simulation tests to determine the volume of oil/gas migration and accumulation within these 3 series of plentiful hydrocarbon sources, different hydrocarbon abundance and pore level have been carried out under the guide of multidisciplinary theories, applying various methods and techniques, and 24 migration/accumulation modes have been established. The innovative results and recognition are as follows: 1, The oil/gas migration and accumulation modes for sandstones of moderate, fine grain size and silt in these six paleo depositional environments of shallow lake, fluvial, lacustrine, fluvial-deltaic, turbidite-delta, and salty-semi salty lake have been established. A new view has been put forward that the oil/gas volumetric increment during their migration and accumulation in different porous media of different rocks has similar features and evolution history. 2. During oil/gas migration and accumulation in different grain-sized porous media or different reservoirs, all the volumetric increment had experienced three period of rapid increasing, balanced and slower increasing and limited increasing. The dynamical process of oil/gas secondary migration and accumulation has been expounded. 3 The two new concepts of "source supply abundance" and "source supply intensity" have been proposed for the first time, and the physical simulation for hydrocarbon's migration, accumulation and forming a reservoir has been realized. 4, Source supply abundance is the important factor to control the accumulated volume of oil phase in the porous media. It is impossible to accumulate large amount of hydrocarbon volume in an open boundary system when the source supply abundance is low, i.e. impossible to form reservoirs of high productivity. 5 The above 12 types of physical simulation tests indicated that enough energy (pressure) of the oil sources is the decisive factor to ensure hydrocarbon's entering, flowing and accumulating through porous media, and both oil and gas phase will accumulate into the favorable places nearest to the oil sources. 6 The theory, method and related techniques for physical simulation of hydrocarbon's secondary migration/accumulation mechanism have been formed and applied to the E&P of Shengtuo rollover anticline and Niuzhuang turbidite lithological reservoirs. 7 This study developed the theory and method of hydrocarbon migration/accumulation dynamics in terrestrial rift-subsidence lacustrine basins, and the benefits and social effect are remarkable.

Boundary triplets and M-functions for non-selfadjoint operators, with applications to elliptic PDEs and block operator matrices

B.M. Brown, M. Marletta, S. Naboko, I. Wood: Boundary triplets and M-functions for non-selfadjoint operators, with applications to elliptic PDEs and block operator matrices, J. London Math. Soc., June 2008; 77: 700-718. The full text of this article will be made available in this repository in June 2009 Sponsorship: EPSRC,INTAS

A computational model of coupled heat and moisture transfer with phase change in granular sugar during varying environmental conditions

As part of a comprehensive effort to predict the development of caking in granular materials, a mathematical model is introduced to model simultaneous heat and moisture transfer with phase change in porous media when undergoing temperature oscillations/cycling. The resulting model partial differential equations were solved using finite-volume procedures in the context of the PHYSICA framework and then applied to the analysis of sugar in storage. The influence of temperature on absorption/desorption and diffusion coefficients is coupled into the transport equations. The temperature profile, the depth of penetration of the temperature oscillation into the bulk solid, and the solids moisture content distribution were first calculated, and these proved to be in good agreement with experimental data. Then, the influence of temperature oscillation on absolute humidity, moisture concentration, and moisture migration for different parameters and boundary conditions was examined. As expected, the results show that moisture near boundary regions responds faster than farther away from them with surface temperature changes. The moisture absorption and desorption in materials occurs mainly near boundary regions (where interactions with the environment are more pronounced). Small amounts of solids moisture content, driven by both temperature and vapour concentration gradients, migrate between boundary and center with oscillating temperature.

Airborne transport of Saharan dust to the Mediterranean and to the Atlantic

The Sahara desert is a significant source of particulate pollution not only to the Mediterranean region, but also to the Atlantic and beyond. In this paper, PM 10 exceedences recorded in the UK and the island of Crete are studied and their source investigated, using Lagrangian Particle Dispersion (LPD) methods. Forward and inverse simulations identify Saharan dust storms as the primary source of these episodes. The methodology used allows comparison between this primary source and other possible candidates, for example large forest fires or volcanic eruptions. Two LPD models are used in the simulations, namely the open source code FLEXPART and the proprietary code HYSPLIT. Driven by the same meteorological fields (the ECMWF MARS archive and the PSU/NCAR Mesoscale model, known as MM5) the codes produce similar, but not identical predictions. This inter-model comparison enables a critical assessment of the physical modelling assumptions employed in each code, plus the influence of boundary conditions and solution grid density. The outputs, in the form of particle concentrations evolving in time, are compared against satellite images and receptor data from multiple ground-based sites. Quantitative comparisons are good, especially in predicting the time of arrival of the dust plume in a particular location.

Correlated electron transport in molecular electronics

Theoretical and experimental values to date for the resistances of single molecules commonly disagree by orders of magnitude. By reformulating the transport problem using boundary conditions suitable for correlated many-electron systems, we approach electron transport across molecules from a new standpoint. Application of our correlated formalism to benzene-dithiol gives current-voltage characteristics close to experimental observations. The method can solve the open system quantum many-body problem accurately, treats spin exactly, and is valid beyond the linear response regime.

Dynamical simulation of inelastic quantum transport

A method for correlated quantum electron-ion dynamics is combined with a method for electronic open boundaries to simulate in real time the heating, and eventual equilibration at an elevated vibrational energy, of a quantum ion under current flow in an atomic wire, together with the response of the current to the ionic heating. The method can also be used to extract inelastic current voltage corrections under steady-state conditions. However, in its present form the open-boundary method contains an approximation that limits the resolution of current-voltage features. The results of the simulations are tested against analytical results from scattering theory. Directions for the improvement of the method are summarized at the end.