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.

Parameters identification by a single point free response measurement

When a thin rectangular plate is restrained on the two long edges and free on the remaining edges, the equivalent stiffness of the restraining joints can be identified by the order of the natural frequencies obtained using the free response of the plate at a single location. This work presents a method to identify the equivalent stiffness of the restraining joints, being represented as simply supporting the plate but elastically restraining it in rotation. An integral transform is used to map the autospectrum of the free response from the frequency domain to the stiffness domain in order to identify the equivalent torsional stiffness of the restrained edges of the plate and also the order of natural frequencies. The kernel of the integral transform is built interpolating data from a finite element model of the plate. The method introduced in this paper can also be applied to plates or shells with different shapes and boundary conditions. © 2011 Elsevier Ltd. All rights reserved.

Vibrations of coupled cylindrical-conical shells

The free vibrational characteristics of coupled conical-cylindrical shells is presented. The equations of motion for the cylindrical shell are solved using a wave approach while the equations of motion for the conical shells are solved using a power series solution. The use of both Donnell-Mushtari and Flügge equations of motion are investigated and their limitations are discussed. Results are presented in terms of natural frequencies for different boundary conditions and the purely torsional mode solution is described. The results from the analytical model presented are compared with those obtained from a finite element model solved with Nastran and other data available in literature.

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.

Guidance of collective cell migration by substrate geometry.

Collective behavior refers to the emergence of complex migration patterns over scales larger than those of the individual elements constituting a system. It plays a pivotal role in biological systems in regulating various processes such as gastrulation, morphogenesis and tissue organization. Here, by combining experimental approaches and numerical modeling, we explore the role of cell density ('crowding'), strength of intercellular adhesion ('cohesion') and boundary conditions imposed by extracellular matrix (ECM) proteins ('constraints') in regulating the emergence of collective behavior within epithelial cell sheets. Our results show that the geometrical confinement of cells into well-defined circles induces a persistent, coordinated and synchronized rotation of cells that depends on cell density. The speed of such rotating large-scale movements slows down as the density increases. Furthermore, such collective rotation behavior depends on the size of the micropatterned circles: we observe a rotating motion of the overall cell population in the same direction for sizes of up to 200 μm. The rotating cells move as a solid body, with a uniform angular velocity. Interestingly, this upper limit leads to length scales that are similar to the natural correlation length observed for unconfined epithelial cell sheets. This behavior is strongly altered in cells that present a downregulation of adherens junctions and in cancerous cell types. We anticipate that our system provides a simple and easy approach to investigate collective cell behavior in a well-controlled and systematic manner.

Bayesian parameter identification for vibro-acoustic models with nonparametric uncertainty

Vibration and acoustic analysis at higher frequencies faces two challenges: computing the response without using an excessive number of degrees of freedom, and quantifying its uncertainty due to small spatial variations in geometry, material properties and boundary conditions. Efficient models make use of the observation that when the response of a decoupled vibro-acoustic subsystem is sufficiently sensitive to uncertainty in such spatial variations, the local statistics of its natural frequencies and mode shapes saturate to universal probability distributions. This holds irrespective of the causes that underly these spatial variations and thus leads to a nonparametric description of uncertainty. This work deals with the identification of uncertain parameters in such models by using experimental data. One of the difficulties is that both experimental errors and modeling errors, due to the nonparametric uncertainty that is inherent to the model type, are present. This is tackled by employing a Bayesian inference strategy. The prior probability distribution of the uncertain parameters is constructed using the maximum entropy principle. The likelihood function that is subsequently computed takes the experimental information, the experimental errors and the modeling errors into account. The posterior probability distribution, which is computed with the Markov Chain Monte Carlo method, provides a full uncertainty quantification of the identified parameters, and indicates how well their uncertainty is reduced, with respect to the prior information, by the experimental data. © 2013 Taylor & Francis Group, London.

One-dimensional quantum waveguide theory of Rashba electrons

The ballistic spin transport in one-dimensional waveguides with the Rashba effect is studied. Due to the Rashba effect, there are two electron states with different wave vectors for the same energy. The wave functions of two Rashba electron states are derived, and it is found that their phase depend on the direction of the circuit and the spin directions of two states are perpendicular to the circuit, with the +pi/2 and -pi/2 angles, respectively. The boundary conditions of the wave functions and their derivatives at the intersection of circuits are given, which can be used to investigate the waveguide transport properties of Rashba spin electron in circuits of any shape and structure. The eigenstates of the closed circular and square loops are studied by using the transfer matrix method. The transfer matrix M(E) of a circular arc is obtained by dividing the circular arc into N segments and multiplying the transfer matrix of each straight segment. The energies of eigenstates in the closed loop are obtained by solving the equation det[M(E)-I]=0. For the circular ring, the eigenenergies obtained with this method are in agreement with those obtained by solving the Schrodinger equation. For the square loop, the analytic formula of the eigenenergies is obtained first The transport properties of the AB ring and AB square loop and double square loop are studied using the boundary conditions and the transfer matrix method In the case of no magnetic field, the zero points of the reflection coefficients are just the energies of eigenstates in closed loops. In the case of magnetic field, the transmission and reflection coefficients all oscillate with the magnetic field; the oscillating period is Phi(m)=hc/e, independent of the shape of the loop, and Phi(m) is the magnetic flux through the loop. For the double loop the oscillating period is Phi(m)=hc/2e, in agreement with the experimental result. At last, we compared our method with Koga's experiment. (C) 2009 American Institute of Physics. [doi: 10.1063/1.3253752]

Micromixing enhanced by pulsating flows

We consider the micromixing enhancement by pulsating flows. Dimensionless governing equations and boundary conditions were developed for T-type micromixers with two inlet pulsating flows. The problem involves a set of parameters. Three key dimensionless parameters are identified: the Reynolds number, the Strouhal number, and the disturbance amplitude. Suitable Strouhal number or disturbance amplitude causes symmetrical meniscus-shape mixing interfaces, separating the whole mixing channel into a set of segments. Thus uniform exit species concentration can be reached. Too large or too small Strouhal number or disturbance amplitude yields the meniscus-shape mixing interfaces deviating from the centerline of the mixing channel, deteriorating the mixing performance. The optimized disturbance amplitude is increased with increases in Strouhal numbers. Low Reynolds number needs larger disturbance amplitude.

Ground state energy of He isoelectronic sequence treated variationally via Hylleraas-like wavefunction

In this study, the energy for the ground state of helium and a few helium-like ions (Z=1-6) is computed variationally by using a Hylleraas-like wavefunction. A four-parameters wavefunction, satisfying boundary conditions for coalescence points, is combined with a Hylleraas-like basis set which explicitly incorporates r12 interelectronic distance. The main contribution of this work is the introduction of modified correlation terms leading to the definition of integral transforms which provide the calculation of expectation value of energy to be done analytically over single-particle coordinates instead of Hylleraas coordinates.

黄土区考虑滞后作用的坡地水量转化模型

建立了考虑滞后作用的坡地水量转化数学模型，并对模型不同情况下的定解条件进行了分析。利用陕西长武王东沟小流域坡地天然降雨－入渗－产流的实测资料对模型进行了系统检验，结果表明该模型具有较高精度，考虑滞后作用时的模拟值更加接近实测值。在降雨入渗过程中，滞后作用对土壤含水率、入渗速率、径流过程、产流量都有相当大的影响。

A theoretical framework of tracer methods for marine sediment dynamics

A new theoretical framework of tracer methods is proposed in the present contribution, on the basis of mass conservation. This model is applicable for both artificial and natural tracers. It can be used to calculate the spatial distribution patterns of sediment transport rate, thus providing independent information and verification for the results derived from empirical formulae. For the procedures of the calculation, first, the tracer concentration and topographic maps of two times are obtained. Then, the spatial and temporal changes in the concentration and seabed elevation are calculated, and the boundary conditions required are determined by field observations (such as flow and bedform migration measurements). Finally, based upon eqs. (1) and (13), the transport rate is calculated and expressed as a function of the position over the study area. Further, appropriate modifications to the model may allow the tracer to have different densities and grain size distributions from the bulk sediment.

Observations on the nocturnal activity and feeding behavior of Anguilla japonica glass eels under laboratory conditions

Glass eels of the temperate anguillid species, Anguilla japonica, clearly showed a nocturnal activity rhythm under laboratory conditions. Light-dark cycle was a determinant factor affecting their photonegative behavior, nocturnal locomotor activity, and feeding behavior. Under natural light conditions, glass eels remained in shelters with little daytime feeding, but came out to forage during darkness. They moved and foraged actively in the following dark, and then their activity gradually declined possibly because of food satiation. They finally buried in the sand or stayed in tubes immediately after the lights came on. Under constant light, glass eels often came out of the shelters to forage in the lights but spent little time moving outside the shelters (e.g. swimming or crawling on the sand). Glass eels took shelter to avoid light and preferred tubes to sand for shelter possibly because tubes were much easier for them to take refuge in than sand. Feeding and locomotor activities of the glass eels were nocturnal and well synchronized. They appeared to depend on olfaction rather than vision to detect and capture prey in darkness. Feeding was the driving force for glass eels to come out of sand under constant light. However, in the dark, some glass eels swam or crept actively on sand even when they were fully fed. The lunar cycles of activity rhythms of glass eels that have been observed in some estuarine areas were not detected under these laboratory conditions.

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.

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

Fixational Instability and Natural Image Statistics: Implications for Early Visual Representations

Under natural viewing conditions small movements of the eye, head, and body prevent the maintenance of a steady direction of gaze. It is known that stimuli tend to fade when they a restabilized on the retina for several seconds. However; it is unclear whether the physiological motion of the retinal image serves a visual purpose during the brief periods of natural visual fixation. This study examines the impact of fixational instability on the statistics of the visua1 input to the retina and on the structure of neural activity in the early visual system. We show that fixational instability introduces a component in the retinal input signals that in the presence of natural images, lacks spatial correlations. This component strongly influences neural activity in a model of the LGN. It decorrelates cell responses even if the contrast sensitivity functions of simulated cells arc not perfectly tuned to counterbalance the power-law spectrum of natural images. A decorrelation of neural activity at the early stages of the visual system has been proposed to be beneficial for discarding statistical redundancies in the input signals. The results of this study suggest that fixational instability might contribute to establishing efficient representations of natural stimuli.