(Table 1) Stable isotope record of Middle Cretaceous benthic foraminifera

Preservation of planktic foraminiferal calcite has received widespread attention in recent years, but the taphonomy of benthic foraminiferal calcite and its influence on the deep-sea palaeotemperature record have gone comparatively unreported. Numerical modeling indicates that the carbonate recrystallization histories of deep-sea sections are dominated by events in their early burial history, meaning that the degree of exchange between sediments and pore fluids during the early postburial phase holds the key to determining the palaeotemperature significance of diagenetic alteration of benthic foraminifera. Postburial sedimentation rate and lithology are likely to be important determinants of the paleoceanographic significance of this sediment-pore fluid interaction. Here we report an investigation of the impact of extreme change in sedimentation rate (a prolonged and widespread Upper Cretaceous hiatus in the North Atlantic Ocean) on the preservation and d18O of benthic foraminifera of Middle Cretaceous age (nannofossil zone NC10, uppermost Albian/lowermost Cenomanian, ~99 Ma ago) from multiple drill sites. At sites where this hiatus immediately overlies NC10, benthic foraminifera appear to display at least moderate preservation of the whole test. However, on closer inspection, these tests are shown to be extremely poorly preserved internally and yield d18O values substantially higher than those from contemporaneous better preserved benthic foraminifera at sites without an immediately overlying hiatus. These high d18O values are interpreted to indicate alteration close to the seafloor in cooler waters during the Late Cretaceous hiatus. Intersite differences in lithology modulate the diagenetic impact of this extreme change in sedimentation rate. Our results highlight the importance of thorough examination of benthic foraminiferal wall structures and lend support to the view that sedimentation rate and lithology are key factors controlling the paleoceanographic significance of diagenetic alteration of biogenic carbonates.

U, FE, AND LI ISOTOPE RATIOS IN THE MINERALIZED AND BARREN AREAS OF THE KOMBOLGIE BASIN

The Paleo- to Meso-Proterozoic Jabiluka unconformity related uranium mine is located within the Alligator River Uranium Field, found in the Northern Territories, Australia. The uranium ore is hosted in the late middle Paleoproterozoic Cahill Formation, which is unconformably overlain by a group of unmetamorphosed conglomerates known as the Kombolgie subgroup. The Kombolgie subgroup provided the source for oxidized basinal brines, carrying U as the mobile form U(VI), which interacted with reducing lithologies in the Cahill formation, thus reducing U(VI) to the solid U(IV), and leading to the precipitation of uraninite (UO2). In order to characterize fluid interaction with the ore body and compare that to areas without mineralization, several isotopic tracers were studied on a series of clay samples from drill core at Jabiluka as well as in barren areas throughout the ARUF. Among the potential tracers, three were selected: U (redox sensitive and recent fluid mobilization), Fe (redox sensitive), and Li (fractionated by hydrothermal fluids and adsorption reactions). δ238U values were found to be closely linked to the mineralogy, with samples with higher K/Al ratios (indicating high illite and low chlorite concentrations) having higher δ238U values. This demonstrates that 235U preferentially absorbs onto the surface of chlorite during hydrothermal circulation. In addition, δ234U values lie far from secular equilibrium (δ234U of 30‰), indicating there was addition or removal of 234U from the surface of the samples from recent (<2.5Ma) interactions of mobile fluids. δ57Fe values were found to be related to lithology and spatially to known uranium deposits. Decreasing δ57Fe values were found with increasing depth to the unconformity in a drill hole directly above the ore zone, but not in drill holes in the barren area. Similarly to δ238U, δ7Li is found to correlate with mineralogy, with higher δ7Li values associated with samples with more chlorite. In addition, higher δ7Li values are found at greater depth throughout the basin, indicating that the direction of the mineralizing fluid circulation was upwards from the Cahill formation to the Kombolgie subgroup.

Fluid-structure interaction analysis of full scale vehicle-barrier impact using coupled SPH-FEA

Portable water-filled barriers (PWFB) are roadside structures used to separate moving traffic from work-zones. Numerical PWFB modelling is preferred in the design stages prior to actual testing. This paper aims to study the fluid-structure interaction of PWFB under vehicular impact using several methods. The strategy to treat water as non-structural mass was proposed and the errors were investigated. It was found that water can be treated with the FEA-NSM model for velocities higher than 80kmh-1. However, full SPH/FEA model is still the best treatment for water and necessary for lower impact velocities. The findings in this paper can be used as guidelines for modelling and designing PWFB.

Performance of buried tunnels subjected to surface blast incorporating fluid structure interaction

This paper uses finite element techniques to investigate the performance of buried tunnels subjected to surface blasts incorporating fully coupled Fluid Structure Interaction and appropriate material models which simulate strain rate effects. Modelling techniques are first validated against existing experimental results and then used to treat the blast induced shock wave propagation and tunnel response in dry and saturated sands. Results show that the tunnel buried in saturated sand responds earlier than that in dry sand. Tunnel deformations decrease with distance from explosive in both sands, as expected. In the vicinity of the explosive, the tunnel buried in saturated sand suffered permanent deformation in both axial and circumferential directions, whereas the tunnel buried in dry sand recovered from most of the axial deformation. Overall, response of the tunnel in saturated sand is more severe for a given blast event and shows the detrimental effect of pore water on the blast response of buried tunnels. The validated modelling techniques developed in this paper can be used to investigate the blast response of tunnels buried in dry and saturated sands.

Fluid-structure interaction analysis by coupled FE-SPH model based on a novel searching algorithm

Fluid–Structure Interaction (FSI) problem is significant in science and engineering, which leads to challenges for computational mechanics. The coupled model of Finite Element and Smoothed Particle Hydrodynamics (FE-SPH) is a robust technique for simulation of FSI problems. However, two important steps of neighbor searching and contact searching in the coupled FE-SPH model are extremely time-consuming. Point-In-Box (PIB) searching algorithm has been developed by Swegle to improve the efficiency of searching. However, it has a shortcoming that efficiency of searching can be significantly affected by the distribution of points (nodes in FEM and particles in SPH). In this paper, in order to improve the efficiency of searching, a novel Striped-PIB (S-PIB) searching algorithm is proposed to overcome the shortcoming of PIB algorithm that caused by points distribution, and the two time-consuming steps of neighbor searching and contact searching are integrated into one searching step. The accuracy and efficiency of the newly developed searching algorithm is studied on by efficiency test and FSI problems. It has been found that the newly developed model can significantly improve the computational efficiency and it is believed to be a powerful tool for the FSI analysis.

Stress analysis of carotid atheroma in a transient ischaemic attack patient using the MRI-based fluid-structure interaction method

Rupture of atherosclerotic plaque is a major cause of mortality. Plaque stress analysis, based on patient-specific multisequence in vivo MRI, can provide critical information for the understanding of plaque rupture and could eventually lead to plaque rupture prediction. However, the direct link between stress and plaque rupture is not fully understood. In the present study, the plaque from a patient who recently experienced a transient ischaemic attack (TIA) was studied using a fluid-structure interaction method to quantify stress distribution in the plaque region based on in vivo MR images. The results showed that wall shear stress is generally low in the artery with a slight increase at the plaque throat owing to minor luminal narrowing. The oscillatory shear index is much higher in the proximal part of the plaque. Both local wall stress concentrations and the relative stress variation distribution during a cardiac cycle indicate that the actual plaque rupture site is collocated with the highest rupture risk region in the studied patient.

Fluid-structure interaction study on human arterial plaque with patient specific geometry and boundary conditions

Atherosclerotic plaque rupture has been extensively considered as the leading cause of death in western countries. It is believed that high stresses within plaque can be an important factor on triggering the rupture of the plaque. Stress analysis in the coronary and carotid arteries with plaque have been developed by many researchers from 2D to 3-D models, from structure analysis only to the Fluid-Structure Interaction (FSI) models[1].

A monolithic strategy for fluid-structure interaction problems

In this work, we present a new monolithic strategy for solving fluid-structure interaction problems involving incompressible fluids, within the context of the finite element method. This strategy, similar to the continuum dynamics, conserves certain properties, and thus provides a rational basis for the design of the time-stepping strategy; detailed proofs of the conservation of these properties are provided. The proposed algorithm works with displacement and velocity variables for the structure and fluid, respectively, and introduces no new variables to enforce velocity or traction continuity. Any existing structural dynamics algorithm can be used without change in the proposed method. Use of the exact tangent stiffness matrix ensures that the algorithm converges quadratically within each time step. An analytical solution is presented for one of the benchmark problems used in the literature, namely, the piston problem. A number of benchmark problems including problems involving free surfaces such as sloshing and the breaking dam problem are used to demonstrate the good performance of the proposed method. Copyright (C) 2010 John Wiley & Sons, Ltd.

Exciting forces due to interaction of water waves with a submerged sphere in an ice-covered two-layer fluid of finite depth

Scattering of water waves by a sphere in a two-layer fluid, where the upper layer has an ice-cover modelled as an elastic plate of very small thickness, while the lower one has a rigid horizontal bottom surface, is investigated within the framework of linearized water wave theory. The effects of surface tension at the surface of separation is neglected. There exist two modes of time-harmonic waves - the one with lower wave number propagating along the ice-cover and the one with higher wave number along the interface. Method of multipole expansions is used to find the particular solution for the problem of wave scattering by a submerged sphere placed in either of the layers. The exciting forces for vertical and horizontal directions are derived and plotted against different values of the wave number for different submersion depths of the sphere and flexural rigidity of the ice-cover. When the flexural rigidity and the density of the ice-cover are taken to be zero, the numerical results for the exciting forces for the problem with free surface are recovered as particular cases. (C) 2011 Elsevier Ltd. All rights reserved.

Fluid-Solid Interaction in Particle-Laden Flows

Liu Qingquan, Singh V.P

Combination Of Cfd And Csd Packages For Fluid-Structure Interaction

In this article the UDF script file in the Fluent software was rewritten as the "connecting file" for the Fluent and the ANSYS/ABAQUS in order that the joined file can be used to do aero-elastic computations. In this way the fluid field is computed by solving the Navier-Stokes equations and the structure movement is integrated by the dynamics directly. An analysis of the computed results shows that this coupled method designed for simulating aero-elastic systems is workable and can be used for the other fluid-structure interaction problems.

Dynamic interaction of an elastic container with fluid

The vibration analysis of an elastic container with partially filled fluid was investigated in this paper. The container is made of a thin cylinder and two circular plates at the ends. The axis of the cylinder is in the horizontal direction. It is difficult to solve this problem because the complex system is not axially symmetric. The equations of motion for this system were derived. An incompressible and ideal fluid model is used in the present work. Solutions of the equations were obtained by the generalized variational method. The solution was expressed in a series of normalized generalized Fourier's functions. This series converged rapidly, and so its approximate solution was obtained with high precision. The agreement of the calculated values with the experimental result is good. It should be mentioned that with our method, the computer time is less than that with the finite-element method.

The interaction of two pairs of solitary waves in a two-fluid system

This paper deals with the interaction of solitary waves in a two-fluid system which consistsof two superimposed incompressible inviscid fluids with a free surface and a horizontal rigidbottom. Under the assumption of shallow water wave, we first derive the basic equationssuitable for the model considered, a generalized form of the Boussinesq equations, then usingthe PLK method and the reductive perturbation method, obtain the second-order approximatesolution for the head-on collision between two pairs of interface and surface solitary waves,and give their maximum amplitudes during the collision and the nonuniform phase shiftsafter the collision which lead to the distortion of the wave profiles.

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment