Dynamic fluid-structure interaction using finite volume unstructured mesh procedures

A three dimensional finite volume, unstructured mesh method for dynamic fluid-structure interation is described. The broad approach is conventional in that the fluid and structure are solved sequentially. The pressure and viscous stresses from the flow algorithm provide load conditions for the solid algorithm, whilst at the fluid structure interface the deformed structure provides boundary condition from the structure to the fluid. The structure algorithm also provides the necessary mesh adaptation for the flow field, the effect of which is accounted for in the flow algorithm. The procedures described in this work have several novel features, namely: * a single mesh covering the entire domain. * a Navier Stokes flow. * a single FV-UM discretisation approach for both the flow and solid mechanics procedures. * an implicit predictor-corrector version of the Newmark algorithm. * a single code embedding the whole strategy. The procedure is illustrated for a three dimensional loaded cantilever in fluid flow.

Two fluid approach to high speed impact interaction amongst solid structures

The issues surrounding collision of projectiles with structures has gained a high profile since the events of 11th September 2001. In such collision problems, the projectile penetrates the stucture so that tracking the interface between one material and another becomes very complex, especially if the projectile is essentially a vessel containing a fluid, e.g. fuel load. The subsequent combustion, heat transfer and melting and re-solidification process in the structure render this a very challenging computational modelling problem. The conventional approaches to the analysis of collision processes involves a Lagrangian-Lagrangian contact driven methodology. This approach suffers from a number of disadvantages in its implementation, most of which are associated with the challenges of the contact analysis component of the calculations. This paper describes a 'two fluid' approach to high speed impact between solid structures, where the objective is to overcome the problems of penetration and re-meshing. The work has been carried out using the finite volume, unstructured mesh multi-physics code PHYSICA+, where the three dimensional fluid flow, free surface, heat transfer, combustion, melting and re-solidification algorithms are approximated using cell-centred finite volume, unstructured mesh techniques on a collocated mesh. The basic procedure is illustrated for two cases of Newtonian and non-Newtonian flow to test various of its component capabilities in the analysis of problems of industrial interest.

A finite volume unstructured mesh approach to dynamic fluid–structure interaction: an assessment of the challenge of predicting the onset of flutter

Computational modelling of dynamic fluid–structure interaction (DFSI) is a considerable challenge. Our approach to this class of problems involves the use of a single software framework for all the phenomena involved, employing finite volume methods on unstructured meshes in three dimensions. This method enables time and space accurate calculations in a consistent manner. One key application of DFSI simulation is the analysis of the onset of flutter in aircraft wings, where the work of Yates et al. [Measured and Calculated Subsonic and Transonic Flutter Characteristics of a 45° degree Sweptback Wing Planform in Air and Freon-12 in the Langley Transonic Dynamic Tunnel. NASA Technical Note D-1616, 1963] on the AGARD 445.6 wing planform still provides the most comprehensive benchmark data available. This paper presents the results of a significant effort to model the onset of flutter for the AGARD 445.6 wing planform geometry. A series of key issues needs to be addressed for this computational approach. • The advantage of using a single mesh, in order to eliminate numerical problems when applying boundary conditions at the fluid-structure interface, is counteracted by the challenge of generating a suitably high quality mesh in both the fluid and structural domains. • The computational effort for this DFSI procedure, in terms of run time and memory requirements, is very significant. Practical simulations require even finer meshes and shorter time steps, requiring parallel implementation for operation on large, high performance parallel systems. • The consistency and completeness of the AGARD data in the public domain is inadequate for use in the validation of DFSI codes when predicting the onset of flutter.

A generic time and space accurate numerical approach to closely coupled fluid-structure interaction problems

Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, flow in elastic pipes and blood vessels and extrusion of metals through dies. However a comprehensive computational model of these multi-physics phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply even to the extent in metal forming, for example, that the deformation of the die is totally ignored. More recently, strategies for solving the full coupling between the fluid and soild mechanics behaviour have developed. Conventionally, the computational modelling of fluid structure interaction is problematical since computational fluid dynamics (CFD) is solved using finite volume (FV) methods and computational structural mechanics (CSM) is based entirely on finite element (FE) methods. In the past the concurrent, but rather disparate, development paths for the finite element and finite volume methods have resulted in numerical software tools for CFD and CSM that are different in almost every respect. Hence, progress is frustrated in modelling the emerging multi-physics problem of fluid structure interaction in a consistent manner. Unless the fluid-structure coupling is either one way, very weak or both, transferring and filtering data from one mesh and solution procedure to another may lead to significant problems in computational convergence. Using a novel three phase technique the full interaction between the fluid and the dynamic structural response are represented. The procedure is demonstrated on some challenging applications in complex three dimensional geometries involving aircraft flutter, metal forming and blood flow in arteries.

Seasonal plankton dynamics along a cross-shelf gradient off northern Spain

The development of population models able to reproduce the dynamics of zooplankton is a central issue when trying to understand how a changing environment would affect zooplankton in the future. Using 10 years of monthly data on phytoplankton and zooplankton abundance in the Bay of Biscay from the IEO's RADIALES time-series programme, we built non-parametric Generalized Additive Models (GAMs) able to reproduce the dynamics of plankton on the basis of environmental factors (nutrients, temperature, upwelling and photoperiod). We found that the interaction between these two plankton components is approximately linear, whereas the effects of environmental factors are non-linear. With the inclusion of the environmental variability, the main seasonal and inter-annual dynamic patterns observed within the studied plankton assemblage indicate the prevalence of bottom-up regulatory control. The statistically deduced models were used to simulate the dynamics of the phytoplankton and zooplankton. A good agreement between observations and simulations was obtained, especially for zooplankton. We are presently developing spatio-temporal GAM models for the North Sea based on the Continuous Plankton Recorder database.

The eigenphase formulation of fine-structure transition cross sections

The eigenphase formulation of Blatt and Biedenharn is applied to fine-structure transitions in *P atoms colliding with ‘S perturbers. Consideration is given to the limit of weak spin-orbit interaction. If the eigenphases are equal to the phaseshifts for elastic scattering by the molecular potentials then the expression for the total cross section reduces to the expression derived in the elastic approximation. However, a numerical comparison for the Li(2p ’P) + He(’S) system shows that the elastic molecular phaseshifts are not good approximations to the eigenphases. Hence the elastic approximation cannot be reliable.

Cross-Kerr-based information transfer processes

The realization of nonclassical states is an important task for many applications of quantum information processing. Usually, properly tailored interactions, different from goal to goal, are considered in order to accomplish specific tasks within the general framework of quantum state engineering. In this paper, we remark on the flexibility of a cross-Kerr nonlinear coupling in hybrid systems as an important ingredient in the engineering of nonclassical states. The general scenario we consider is the implementation of high cross-Kerr nonlinearity in cavity-quantum electrodynamics. In this context, we discuss the possibility of performing entanglement transfer and swapping between matter qubits and light fields initially prepared in separable coherent states. The recently introduced concept of entanglement reciprocation is also considered and shown to be possible with our scheme. We reinterpret some of our results in terms of applications of a generalized Ising interaction to systems of different nature.

Measurement of charge exchange and X-ray emission cross sections for solar wind-comet interactions

X-ray emission from a comet was observed for the first time in 1996. One of the mechanisms believed to be contributing to this surprisingly strong emission is the interaction of highly charged solar wind ions with cometary gases. Reported herein are total absolute charge-exchange and normalized line-emission (X-ray) cross sections for collisions of high-charge state (+3 to +10) C, N, O, and Ne ions with the cometary species H2O and CO2. It is found that in several cases the double charge-exchange cross sections can be large, and in the case of C3+ they are equal to those for single charge exchange. Present results are compared to cross section values used in recent comet models. The importance of applying accurate cross sections, including double charge exchange, to obtain absolute line-emission intensities is emphasized.

Liquid monomer-powder particle interaction in acrylic bone cement

It is known that the method used to mix the liquid monomer and powder of PMMA bone cement influences the quality of the cement that is used in total joint replacements. Mixing theory indicates that the interaction between the liquid monomer and the powder is affected by a number of parameters, such as cement viscosity and degree of agitation, with this knowledge utilized in the design of cement mixing devices. Therefore, the objectives of this study were to: (i) obtain information on the interaction of the liquid monomer and the powder in the case of an PMMA bone cement, (ii) show how this knowledge can be applied to the design of an automated cement mixing device, and (iii) compare the porosity, bending modulus, and bending strength of one commercially-available cement prepared using the automated mixer and prepared using a conventional mixer that is in current clinical use. Experimental data indicated that increasing the velocity and decreasing the viscosity of the systems produced cement that improved mechanical properties, which may contribute to better mechanical integrity and, hence, reduced tendency for aseptic loosening, of cemented hip implants.

Measurements of absolute, single charge-exchange cross sections of H+, He+ and He2+ with H2O and CO2

Absolute measurements have been made of single-electron charge-exchange cross sections of H+, He+, and He2+ in H2O and CO2 in the energy range 0.3-7.5 keV amu(-1). Collisions of this type occur in the interaction of solar wind ions with cometary gases and have been observed by the Giotto spacecraft using the Ion Mass Spectrometer/High Energy Range Spectrometer (IMS/HERS) during a close encounter with comet Halley in 1986. Increases in the He+ ion density, and in the He2+ to H+ density ratio were reported by Shelley et al, and Fuselier et al. and were explained by charge exchange. However, the lack of reliable cross sections for this process made interpretation of the data difficult. New cross sections are presented and discussed in relation to the Giotto observations.

Generation of entangled coherent states via cross-phase-modulation in a double electromagnetically induced transparency regime

The generation of an entangled coherent state is one of the most important ingredients of quantum information processing using coherent states. Recently, numerous schemes to achieve this task have been proposed. In order to generate travelling-wave entangled coherent states, cross-phase-modulation, optimized by optical Kerr effect enhancement in a dense medium in an electromagnetically induced transparency (EIT) regime, seems to be very promising. In this scenario, we propose a fully quantized model of a double-EIT scheme recently proposed [D. Petrosyan and G. Kurizki, Phys. Rev. A 65, 33 833 (2002)]: the quantization step is performed adopting a fully Hamiltonian approach. This allows us to write effective equations of motion for two interacting quantum fields of light that show how the dynamics of one field depends on the photon-number operator of the other. The preparation of a Schrodinger cat state, which is a superposition of two distinct coherent states, is briefly exposed. This is based on nonlinear interaction via double EIT of two light fields (initially prepared in coherent states) and on a detection step performed using a 50:50 beam splitter and two photodetectors. In order to show the entanglement of an entangled coherent state, we suggest to measure the joint quadrature variance of the field. We show that the entangled coherent states satisfy the sufficient condition for entanglement based on quadrature variance measurement. We also show how robust our scheme is against a low detection efficiency of homodyne detectors.

Interaction between interleukin 3 and dystrobrevin-binding protein 1 in schizophrenia

Schizophrenia is a common psychotic mental disorder that is believed to result from the effects of multiple genetic and environmental factors. In this study, we explored gene-gene interactions and main effects in both case-control (657 cases and 411 controls) and family-based (273 families, 1350 subjects) datasets of English or Irish ancestry. Fifty three markers in 8 genes were genotyped in the family sample and 44 markers in 7 genes were genotyped in the case-control sample. The Multifactor Dimensionality Reduction Pedigree Disequilibrium Test (MDR-PDT) was used to examine epistasis in the family dataset and a 3-locus model was identified (permuted p=0.003). The 3-locus model involved the IL3 (rs2069803), RGS4 (rs2661319), and DTNBP1 (rs21319539) genes. We used MDR to analyze the case-control dataset containing the same markers typed in the RGS4, IL3 and DTNBP1 genes and found evidence of a joint effect between IL3 (rs31400) and DTNBP1 (rs760761) (cross-validation consistency 4/5, balanced prediction accuracy=56.84%, p=0.019). While this is not a direct replication, the results obtained from both the family and case-control samples collectively suggest that IL3 and DTNBP1 are likely to interact and jointly contribute to increase risk for schizophrenia. We also observed a significant main effect in DTNBP1, which survived correction for multiple comparisons, and numerous nominally significant effects in several genes. (C) 2008 Elsevier B.V. All rights reserved.

Structural basis for the platelet-collagen interaction - The smallest motif within collagen that recognizes and activates platelet glycoprotein VI contains two glycine-proline-hydroxyproline triplets

Collagen-related peptide is a selective agonist for the platelet collagen receptor Glycoprotein VI. The triple helical peptide contains ten GPO triplets/strand (single letter amino acid nomenclature, where O is hydroxyproline) and so over-represents GPO compared with native collagen sequence. To investigate the ability of Glycoprotein VI to recognize GPO triplets in a setting more representative of the collagens, we synthesized a set of triple helical peptides containing fewer GPO triplets, varying their number and spacing within an inert (GPP)(n) backbone. The adhesion of recombinant human Glycoprotein VI ectodomain, like that of human platelets, to these peptides increased with their GPO content, and platelet adhesion was abolished by the specific anti-Glycoprotein VI-blocking antibody, 10B12. Platelet aggregation and protein tyrosine phosphorylation were induced only by cross-linked peptides and only those that contained two or more GPO triplets. Such peptides were less potent than cross-linked collagen-related peptide. Our data suggest that both the sequences GPOGPO and GPO center dot center dot center dot center dot center dot center dot center dot center dot center dot GPO represent functional Glycoprotein VI recognition motifs within collagen. Furthermore, we propose that the (GPO)(4) motif can support simultaneous binding of two glycoprotein VI molecules, in either a parallel or anti-parallel stacking arrangement, which could play an important role in activation of signaling.

Parenting stress and children with cerebral palsy: A European cross-sectional survey

Aim
The aim of this study was to describe stress in the parents of children with cerebral palsy and investigate associations with very high stress.

Method
A cross-sectional survey was conducted of parents of 818 children aged 8 to 12 years from nine regions in Europe. Families were eligible to participate if they were living in one of the specified geographic areas. Parental stress was captured using the Parenting Stress Index Short Form, which has 36 items and takes 10 minutes to complete. Parents rate items on a 5-point Likert scale, with higher scores indicating higher stress. The Short Form yields scores on three subscales and a Total Stress score. A trained research associate administered the questionnaire in the child’s home and visits lasted 90 to 120 minutes. All data collected were reported by parents unless otherwise stated.

Results
The Total Stress score on the Parenting Stress Index was dichotomized into scores of less than 99 or 99 or more, the latter indicating ‘very high’ stress. Most respondents were mothers (94%), and 26% reported very high stress levels. The parents of children with communication impairment had higher odds for very high stress (odds ratio [OR] 1.9; 95% confidence interval [CI] 1.2–3.0) than those whose child had no such impairment; the parents of children with moderate or severe pain had higher odds for very high stress (OR 1.7 [95% CI 1.1–2.4] and 2.5 [95% CI 1.5–4.3] respectively) than those whose child had no pain; and the parents of children with an intellectual impairment had higher odds for very high stress (OR 1.8; 95% CI 1.2–2.9) than those whose child had none. There was no association between very high stress and motor impairment. The subscales ‘parent–child dysfunctional interaction’ and ‘difficult child’ contributed most to the Total Stress score.

Interpretation
Parents of children with communication difficulties, intellectual impairment, or pain are at very high risk of stress. The final model explained 12% of the observed variation in very high stress.