A comparison of methods for the coupled analysis of floating structures

The dynamic analysis of a deepwater floating platform and the associated mooring/riser system should ideally be fully coupled to ensure a reliable response prediction. It is generally held that a time domain analysis is the only means of capturing the various coupling and nonlinear effects accurately. However, in recent work it has been found that for an ultra-deepwater floating system (2000m water depth), the highly efficient frequency domain approach can provide highly accurate response predictions. One reason for this is the accuracy of the drag linearization procedure over both first and second order motions, another reason is the minimal geometric nonlinearity displayed by the mooring lines in deepwater. In this paper, the aim is to develop an efficient analysis method for intermediate water depths, where both mooring/vessel coupling and geometric nonlinearity are of importance. It is found that the standard frequency domain approach is not so accurate for this case and two alternative methods are investigated. In the first, an enhanced frequency domain approach is adopted, in which line nonlinearities are linearized in a systematic way. In the second, a hybrid approach is adopted in which the low frequency motion is solved in the time domain while the high frequency motion is solved in the frequency domain; the two analyses are coupled by the fact that (i) the low frequency motion affects the mooring line geometry for the high frequency motion, and (ii) the high frequency motion affects the drag forces which damp the low frequency motion. The accuracy and efficiency of each of the methods are systematically compared. Copyright © 2007 by ASME.

Fluid mechanics of pulse pressure-gain combustors

The fluid dynamic operation of a valveless pulse combustor has been studied experimentally and numerically. Through phase-locked chemiluminescence and pressure measurements it is shown that mechanical energy is created periodically in the flame surface, with an efficiency of 1.6%. This mechanical energy leaves the pulse combustor through unsteady jets at the aerovalve inlet and the tailpipe exit stations. Two thermodynamically distinct flows are identified: a flow that is transported from inlet to exit and participates in combustion along the way, and a flow that is ingested and then ejected from the combustor without undergoing combustion. It is the latter of these two flows which has the greatest quantity of net work done on it. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc.

Fluid Dynamics in Physics, Engineering and Environmental Applications

The book contains invited lectures and selected contributions presented at the Enzo Levi and XVII Annual Meeting of the Fluid Dynamic Division of the Mexican Physical Society in 2011.

Swimming dynamics of bidirectional artificial flagella

We study magnetic artificial flagella whose swimming speed and direction can be controlled using light and magnetic field as external triggers. The dependence of the swimming velocity on the system parameters (e.g., length, stiffness, fluid viscosity, and magnetic field) is explored using a computational framework in which the magnetostatic, fluid dynamic, and solid mechanics equations are solved simultaneously. A dimensionless analysis is carried out to obtain an optimal combination of system parameters for which the swimming velocity is maximal. The swimming direction reversal is addressed by incorporating photoresponsive materials, which in the photoactuated state can mimic natural mastigonemes. © 2013 American Physical Society.

Swimming dynamics of bidirectional artificial flagella.

We study magnetic artificial flagella whose swimming speed and direction can be controlled using light and magnetic field as external triggers. The dependence of the swimming velocity on the system parameters (e.g., length, stiffness, fluid viscosity, and magnetic field) is explored using a computational framework in which the magnetostatic, fluid dynamic, and solid mechanics equations are solved simultaneously. A dimensionless analysis is carried out to obtain an optimal combination of system parameters for which the swimming velocity is maximal. The swimming direction reversal is addressed by incorporating photoresponsive materials, which in the photoactuated state can mimic natural mastigonemes.

An important pitfall of pseudo-static finite element analysis

Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified methods of dynamic analysis and time domain analysis. The pseudo-static FE approach can accurately model the in situ, stresses prior to seismic loading (when it follows a static analysis simulating the construction sequence) is relatively simple and not as computationally expensive as the time domain approach. However this method should be used with caution as the results can be sensitive to the choice of the mesh dimensions. In this paper two simple examples of pseudo-static finite element analysis are examined parametrically, a homogeneous slope and a cantilever retaining wall, exploring the sensitivity of the pseudo-static analysis results on the adopted mesh size. The mesh dependence was found to be more pronounced for problems with high critical seismic coefficients values (e.g. gentle slopes or small walls), as in these cases a generalised layer failure mechanism is developed simultaneously with the slope or wall mechanism. In general the mesh width was found not to affect notably the predicted value of critical seismic coefficient but to have a major impact on the predicted movements. © 2012 Elsevier Ltd.

Analysis of current induced by long internal solitary waves in stratified ocean

An approximate theoretical expression for the current induced by long internal solitary waves is presented when the ocean is continuously or two-layer stratified. Particular attention is paid to characterizing velocity fields in terms of magnitude, flow components, and their temporal evolution/spatial distribution. For the two-layer case, the effects of the upper/lower layer depths and the relative layer density difference upon the induced current are further studied. The results show that the horizontal components are basically uniform in each layer with a shear at the interface. In contrast, the vertical counterparts vary monotonically in the direction of the water depth in each layer while they change sign across the interface or when the wave peak passes through. In addition, though the vertical components are generally one order of magnitude smaller than the horizontal ones, they can never be neglected in predicting the heave response of floating platforms in gravitationally neutral balance. Comparisons are made between the partial theoretical results and the observational field data. Future research directions regarding the internal wave induced flow field are also indicated.

Dynamic Measurement System for Dipoles of RCS for CS

A dynamic measurement system was developed by the Institute of Modern Physics (IMP) for the dipole prototype of Rapid Cycle Synchrotron (RCS) of China Spallation Neutron Source (CSNS). The repetition frequency of RCS is 25 Hz. The probe is a moving arc searching-coil, and the data acquisition system is based on the dynamic analysis modular of National Instrument. To get the error of high order harmonics of the field at basic frequency, the hardware integrator is replaced by a high speed ADC with software filter and integrator. A series of harmonic coefficients of field are used to express the varieties of dynamic fields in space and time simultaneously. The measurement system has been tested in Institute of High Energy Physics (IHEP), and the property of the dipole prototype of RCS has been measured. Some measurement results and the repeatability of system are illustrated in this paper.

Influence of the earth rotation on the interfacial wave solutions and wave-induced tangential stress in a two-layer fluid system

Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer fluid with a top free surface and a flat bottom. The solutions were deduced from the general form of linear fluid dynamic equations of two-layer fluid under the f-plane approximation, and wave-induced tangential stress were estimated based on the solutions obtained. As expected; the solutions derived from the present work include as special cases those obtained by Sun et al. (2004. Science in China, Set. D, 47(12): 1147-1154) for geostrophic small amplitude surface wave solutions and wave-induced tangential stress if tire density of the upper layer is much smaller than that of the lower layer. The results show that the interface and the surface will oscillate synchronously, and the influence of the earth's rotation both on the surface wave solutions and the interfacial wave solutions should be considered.

三维油藏数值模拟正反演算法研究

The dynamic prediction of complex reservoir development is one of the important research contents of dynamic analysis of oil and gas development. With the increase development of time, the permeabilities and porosities of reservoirs and the permeability of block reservoir at its boundaries are dynamically changing. How to track the dynamic change of permeability and porosity and make certain the permeability of block reservoir at its boundary is an important practical problem. To study developing dynamic prediction of complex reservoir, the key problem of research of dynamic prediction of complex reservoir development is realizing inversion of permeability and porosity. To realize the inversion, first of all, the fast forward and inverse method of 3-dimension reservoir simulation must be studied. Although the inversion has been widely applied to exploration and logging, it has not been applied to3-dimension reservoir simulation. Therefore, the study of fast forward and inverse method of 3-dimension reservoir simulation is a cutting-edge problem, takes on important realistic signification and application value. In this dissertation, 2-dimension and 3-dimension fluid equations in porous media are discretized by finite difference, obtaining finite difference equations to meet the inner boundary conditions by Peaceman's equations, giving successive over relaxation iteration of 3-dimension fluid equations in porous media and the dimensional analysis. Several equation-solving methods are compared in common use, analyzing its convergence and convergence rate. The alternating direction implicit procedure of 2-dimension has been turned into successive over relaxation iteration of alternating direction implicit procedure of 3-dimension fluid equations in porous media, which possesses the virtues of fast computing speed, needing small memory of computer, good adaptability for heterogeneous media and fast convergence rate. The geological model of channel-sandy reservoir has been generated with the help of stochastic simulation technique, whose cross sections of channel-sandy reservoir are parabolic shapes. This method makes the hard data commendably meet, very suit for geological modeling of containing complex boundary surface reservoir. To verify reliability of the method, theoretical solution and numerical solution are compared by simplifying model of 3-dimension fluid equations in porous media, whose results show that the only difference of the two pressure curves is that the numerical solution is lower than theoretical at the wellbore in the same space. It proves that using finite difference to solve fluid equations in porous media is reliable. As numerical examples of 3-dimension heterogeneous reservoir of the single-well and multi-well, the pressure distributions have been computed respectively, which show the pressure distributions there are clearly difference as difference of the permeabilities is greater than one order of magnitude, otherwise there are no clearly difference. As application, the pressure distribution of the channel-sandy reservoir have been computed, which indicates that the space distribution of pressure strongly relies on the direction of permeability, and is sensitive for space distributions of permeability. In this dissertation, the Peaceman's equations have been modified into solving vertical well problem and horizontal well problem simultaneously. In porous media, a 3D layer reservoir in which contain vertical wells and horizontal wells has been calculated with iteration. For channel-sandy reservoir in which there are also vertical wells and horizontal wells, a 3D transient heterogeneous fluid equation has been discretized. As an example, the space distribution of pressure has been calculated with iteration. The results of examples are accord with the fact, which shows the modification of Peaceman's equation is correct. The problem has been solved in the space where there are vertical and horizontal wells. In the dissertation, the nonuniform grid permeability integration equation upscaling method, the nonuniform grid 2D flow rate upscaling method and the nonuniform grid 3D flow rate upscaling method have been studied respectively. In those methods, they enhance computing speed greatly, but the computing speed of 3D flow rate upscaling method is faster than that of 2D flow rate upscaling method, and the precision of 3D flow rate upscaling method is better than that of 2D flow rate upscaling method. The results also show that the solutions of upscaling method are very approximating to that of fine grid blocks. In this paper, 4 methods of fast adaptive nonuniform grid upscaling method of 3D fluid equations in porous media have been put forward, and applied to calculate 3D heterogeneous reservoir and channel-sandy reservoir, whose computing results show that the solutions of nonuniform adaptive upscaling method of 3D heterogeneous fluid equations in porous media are very approximating to that of fine grid blocks in the regions the permeability or porosity being abnormity and very approximating to that of coarsen grid blocks in the other region, however, the computing speed of adaptive upscaling method is 100 times faster than that of fine grid block method. The formula of sensitivity coefficients are derived from initial boundary value problems of fluid equations in porous media by Green's reciprocity principle. The sensitivity coefficients of wellbore pressure to permeability parameters are given by Peaceman's equation and calculated by means of numerical calculation method of 3D transient anisotropic fluid equation in porous media and verified by direct method. The computing results are in excellent agreement with those obtained by the direct method, which shows feasibility of the method. In the dissertation, the calculating examples are also given for 3D reservoir, channel-sandy reservoir and 3D multi-well reservoir, whose numerical results indicate: around the well hole, the value of the sensitivity coefficients of permeability is very large, the value of the sensitivity coefficients of porosity is very large too, but the sensitivity coefficients of porosity is much less than the sensitivity coefficients of permeability, so that the effect of the sensitivity coefficients of permeability for inversion of reservoir parameters is much greater than that of the sensitivity coefficients of porosity. Because computing the sensitivity coefficients needs to call twice the program of reservoir simulation in one iteration, realizing inversion of reservoir parameters must be sustained by the fast forward method. Using the sensitivity coefficients of permeability and porosity, conditioned on observed valley erosion thickness in wells (hard data), the inversion of the permeabilities and porosities in the homogeneous reservoir, homogeneous reservoir only along the certain direction and block reservoir are implemented by Gauss-Newton method or conjugate gradient method respectively. The results of our examples are very approximating to the real data of permeability and porosity, but the convergence rate of conjugate gradient method is much faster than that of Gauss-Newton method.

Nonlinear MHD stability of aluminium reduction cells

An industrial electrolysis cell used to produce primary aluminium is sensitive to waves at the interface of liquid aluminium and electrolyte. The interface waves are similar to stratified sea layers [1], but the penetrating electric current and the associated magnetic field are intricately involved in the oscillation process, and the observed wave frequencies are shifted from the purely hydrodynamic ones [2]. The interface stability problem is of great practical importance because the electrolytic aluminium production is a major electrical energy consumer, and it is related to environmental pollution rate. The stability analysis was started in [3] and a short summary of the main developments is given in [2]. Important aspects of the multiple mode interaction have been introduced in [4], and a widely used linear friction law first applied in [5]. In [6] a systematic perturbation expansion is developed for the fluid dynamics and electric current problems permitting reduction of the three-dimensional problem to a two dimensional one. The procedure is more generally known as “shallow water approximation” which can be extended for the case of weakly non-linear and dispersive waves. The Boussinesq formulation permits to generalise the problem for non-unidirectionally propagating waves accounting for side walls and for a two fluid layer interface [1]. Attempts to extend the electrolytic cell wave modelling to the weakly nonlinear case have started in [7] where the basic equations are derived, including the nonlinearity and linear dispersion terms. An alternative approach for the nonlinear numerical simulation for an electrolysis cell wave evolution is attempted in [8 and references there], yet, omitting the dispersion terms and without a proper account for the dissipation, the model can predict unstable waves growth only. The present paper contains a generalisation of the previous non linear wave equations [7] by accounting for the turbulent horizontal circulation flows in the two fluid layers. The inclusion of the turbulence model is essential in order to explain the small amplitude self-sustained oscillations of the liquid metal surface observed in real cells, known as “MHD noise”. The fluid dynamic model is coupled to the extended electromagnetic simulation including not only the fluid layers, but the whole bus bar circuit and the ferromagnetic effects [9].

Comparison of Experimental PIV Data and CFD Simulations for Flow in a Diesel Particulate Filter Inlet Diffuser

Flow maldistribution of the exhaust gas entering a Diesel Particulate Filter (DPF) can cause uneven soot distribution during loading and excessive temperature gradients during the regeneration phase. Minimising the magnitude of this maldistribution is therefore an important consideration in the design of the inlet pipe and diffuser, particularly in situations where packaging constraints dictate bends in the inlet pipe close to the filter, or a sharp diffuser angle. This paper describes the use of Particle Image Velocimetry (PIV) to validate a Computational Fluid Dynamic (CFD) model of the flow within the inlet diffuser of a DPF so that CFD can be used with confidence as a tool to minimise this flow maldistribution. PIV is used to study the flow of gas into a DPF over a range of steady state flow conditions. The distribution of flow approaching the front face of the substrate was of particular interest to this study. Optically clear diffusing cones were designed and placed between pipe and substrate to allow PIV analysis to take place. Stereoscopic PIV was used to eliminate any error produced by the optical aberrations caused by looking through the curved wall of the inlet cone. In parallel to the experiments, numerical analysis was carried out using a CFD program with an incorporated DPF model. Boundary conditions for the CFD simulations were taken from the experimental data, allowing an experimental validation of the numerical results. The CFD model incorporated a DPF model, the cement layers seen in segmented filters and the intumescent matting that is commonly used to pack the filter into a metal casing. The mesh contained approximately 580,000 cells and used the realizable ?-e turbulence model. The CFD simulation predicted both pressure drop across the DPF and the velocity field within the cone and at the DPF face with reasonable accuracy, providing confidence in the use the CFD in future work to design new, more efficient cones.

Finite element modeling of debonding failures in FRP-strengthened RC beams: A dynamic approach

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton-Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results.

Finite Element Analysis of the Enhanced Strength of Laterally Restrained RC Slabs

This paper presents the numerical simulation of the ultimate behaviour of 85 one-way and two-way spanning laterally restrained concrete slabs of variable thickness, span, reinforcement ratio, strength and boundary conditions reported in literature by different authors. The developed numerical model was described and all the assumptions were illustrated. ABAQUS, a Finite Element Analysis suite of software, was employed. Non-linear implicit static general analysis method offered by ABAQUS was used. Other analysis methods were also discussed in general in terms of application such as Explicit Dynamic Analysis and Riks method. The aim is to demonstrate the ability and efficacy of FEA to simulate the ultimate load behaviour of slabs considering different material properties and boundary conditions. The authors intended to present a numerical model that provides consistent predictions of the ultimate behaviour of laterally restrained slabs that could be used as an alternative for expensive real life testing as well as for the design and assessment of new and existing structures respectively. The enhanced strength of laterally-restrained slabs compared with conventional design methods predictions is believed to be due to compressive membrane action (CMA). CMA is an inherent phenomenon of laterally restrained concrete beams/slabs. The numerical predictions obtained from the developed model were in good correlation with the experimental results and with those obtained from the CMA method developed at the Queen’s University Belfast, UK.

Elastodynamic Modeling and Analysis for an Exechon Parallel Kinematic Machine

As a newly invented parallel kinematic machine (PKM), Exechon has attracted intensive attention from both academic and industrial fields due to its conceptual high performance. Nevertheless, the dynamic behaviors of Exechon PKM have not been thoroughly investigated because of its structural and kinematic complexities. To identify the dynamic characteristics of Exechon PKM, an elastodynamic model is proposed with the substructure synthesis technique in this paper. The Exechon PKM is divided into a moving platform subsystem, a fixed base subsystem and three limb subsystems according to its structural features. Differential equations of motion for the limb subsystem are derived through finite element (FE) formulations by modeling the complex limb structure as a spatial beam with corresponding geometric cross sections. Meanwhile, revolute, universal, and spherical joints are simplified into virtual lumped springs associated with equivalent stiffnesses and mass at their geometric centers. Differential equations of motion for the moving platform are derived with Newton's second law after treating the platform as a rigid body due to its comparatively high rigidity. After introducing the deformation compatibility conditions between the platform and the limbs, governing differential equations of motion for Exechon PKM are derived. The solution to characteristic equations leads to natural frequencies and corresponding modal shapes of the PKM at any typical configuration. In order to predict the dynamic behaviors in a quick manner, an algorithm is proposed to numerically compute the distributions of natural frequencies throughout the workspace. Simulation results reveal that the lower natural frequencies are strongly position-dependent and distributed axial-symmetrically due to the structure symmetry of the limbs. At the last stage, a parametric analysis is carried out to identify the effects of structural, dimensional, and stiffness parameters on the system's dynamic characteristics with the purpose of providing useful information for optimal design and performance improvement of the Exechon PKM. The elastodynamic modeling methodology and dynamic analysis procedure can be well extended to other overconstrained PKMs with minor modifications.