Fluctuating hydrodynamic methods for fluid-structure interactions in confined channel geometries

We develop computational teractions subject to thermal fluctuations geometry. The methods take into account methods for the study of fluid-structure in- when confined within channels with slit-like the hydrodynamic coupling and diffusivity of microstructures when influenced by their proximity to no-slip walls. We develop stochas- tic numerical methods subject to no-slip boundary conditions using a staggered finite volume discretization. We introduce techniques for discretizing stochastic systems in a manner that ensures results consistent with statistical mechanics. We show how an exact fluctuation-dissipation condition can be used for this purpose to discretize the stochastic driving fields and combined with an exact projection method to enforce incompressibil- ity. We demonstrate our computational methods by investigating how the proximity of ellipsoidal colloids to the channel wall affects their active hydrodynamic responses and passive diffusivity. We also study for a large number of interacting particles collective drift-diffusion dynamics and associated correlation h/actions. We expect the introduced stochastic computational methods to be broadly applicable to applications in which con- finement effects play an important role in the dynamics of microstructures subject to hydrodynamic coupling and thermal fluctuations.

Investigation on charge parameters of underwater contact explosion based on axisymmetric SPH method

This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.

Numerical Simulation of Water Mitigation Effects on Shock Wave with SPH Method

The water mitigation effect on the propagation of shock wave was investigated numerically. The traditional smoothed particle hydrodynamics (SPH) method was modified based on Riemann solution. The comparison of numerical results with the analytical solution indicated that the modified SPH method has more advantages than the traditional SPH method. Using the modified SPH algorithm, a series of one-dimensional planar wave propagation problems were investigated, focusing on the influence of the air-gap between the high-pressure air and water and the thickness of water. The numerical results showed that water mitigation effect is significant. Up to 60% shock wave pressure reduction could be achieved with the existence of water, and the shape of shock wave was also changed greatly. It is seemly that the small air-gap between the high-pressure air and water has more influence on water mitigation effect.

Modeling of fluid-induced vibrations and identification of hydrodynamic forces on flow control valves

Dynamics and vibration of control valves under flow-induced vibration are analyzed. Hydrodynamic load characteristics and structural response under flow-induced vibration are mainly influenced by inertia, damping, elastic, geometric characteristics and hydraulic parameters. The purpose of this work is to investigate the dynamic behavior of control valves in the response to self-excited fluid flow. An analytical and numerical method is developed to simulate the dynamic and vibrational behavior of sliding dam valves, in response to flow excitation. In order to demonstrate the effectiveness of proposed model, the simulation results are validated with experimental ones. Finally, to achieve the optimal valve geometry, numerical results for various shapes of valves are compared. Rounded valve with the least amount of flow turbulence obtains lower fluctuations and vibration amplitude compared with the flat and steep valves. Simulation results demonstrate that with the optimal design requirements of valves, vibration amplitude can be reduced by an average to 30%.

Numerical analysis of submarine landslides using a smoothed particle hydrodynamics depth integral model

Submarine landslides can cause severe damage to marine engineering structures. Their sliding velocity and runout distance are two major parameters for quantifying and analyzing the risk of submarine landslides.Currently, commercial calculation programs such as BING have limitations in simulating underwater soil movements. All of these processes can be consistently simulated through a smoothed particle hydrodynamics（SPH） depth integrated model. The basis of the model is a control equation that was developed to take into account the effects of soil consolidation and erosion. In this work, the frictional rheological mode has been used to perform a simulation study of submarine landslides. Time-history curves of the sliding body＇s velocity, height,and length under various conditions of water depth, slope gradient, contact friction coefficient, and erosion rate are compared; the maximum sliding distance and velocity are calculated; and patterns of variation are discussed.The findings of this study can provide a reference for disaster warnings and pipeline route selection.

Prediction for asphalt pavement water film thickness based on artificial neural network

In order to study the variation o f the asphalt pavement water film thickness influenced by multi-factors,anew method for predicting water film thickness was developed by the combination o f the artificial neural network(ANN)a d two-dimensional shallow water equations based on hydrodynamic theory.Multi-factors included the rainfall intensity,pavement width,cross slope,longitudinal slope a d pavement roughness coefficient.The two-dimensional hydrodynamic method was validated by a natural rainfall event.Based on the design scheme o f Shen-Sha expressway engineering project,the limited training data obtained by the two-dimensional hydrodynamic simulation model was used to predict water film thickness.Furthermore,the distribution of the water film thickness influenced by multi-factors on the pavement was analyzed.The accuracy o f the ANN model was verified by the18sets o f data with a precision o f0.991.The simulation results indicate that the water film thickness increases from the median strip to the edge o f the pavement.The water film thickness variation is obviously influenced by rainfall intensity.Under the condition that the pavement width is20m and t e rainfall intensity is3m m/h,t e water film thickness is below10mm in the fast lane and20mm in t e lateral lane.Athough there is fluctuation due to the amount oftraining data,compared with the calculation on the basis o f the existing criterion and theory,t e ANN model exhibits a better performance for depicting the macroscopic distribution of the asphalt pavement water film.

Numerical modeling of sediment transport based on unsteady and steady flows by incompressible smoothed particle hydrodynamics method

The purpose of the present paper is to introduce a simple two-part multi-phase model for the sediment transport problems based on the incompressible smoothed particle hydrodynamics（ISPH） method. The proposed model simulates the movement of sediment particles in two parts. The sediment particles are classified into three categories, including the motionless particles, moving particles behave like a rigid body, and moving particles with a pseudo fluid behavior. The criterion for the classification of sediment particles is the Bingham rheological model. Verification of the present model is performed by simulation of the dam break waves on movable beds with different conditions and the bed scouring under steady flow condition. Comparison of the present model results, the experimental data and available numerical results show that it has good ability to simulate flow pattern and sediment transport.

Dependence of the length of solar filament threads on the magnetic configuration

High-resolution Ha observations indicate that filaments consist of an as- sembly of thin threads. In quiescent filaments, the threads are generally short, whereas in active region filaments, the threads are generally long. In order to explain these observational features, we performed one-dimensional radiative hydrodynamic sim- ulations of filament formation along a dipped magnetic flux tube in the framework of the chromospheric evaporation-coronal condensation model. The geometry of a dipped magnetic flux tube is characterized by three parameters, i.e., the depth （D）, the half-width （w） and the altitude （h） of the magnetic dip. A survey of the parame- ters in numerical simulations shows that when allowing the filament thread to grow in 5 days, the maximum length （Lth） of the filament thread increases linearly with w, and decreases linearly with D and h. The dependence is fitted into a linear function Lth = 0.84w --0.88D - 2.78h ＋ 17.31 （Mm）. Such a relation can qualitatively explain why quiescent filaments have shorter threads and active region filaments have longer threads.

An explanation for the peculiar periphery of supernova remnant G309.2–0.6

Supernova remnant（SNR） G309.2–0.6 has a peculiar radio morphology with two bright ears to the southwest and northeast, although the main shell outside the ears is roughly circular. Based on an earlier proposal that the supernova ejecta has a jet component with extra energy, the dynamical evolution of the remnant is solved using 3 D hydrodynamical（HD） simulation to investigate the formation of the periphery of the remnant. Assuming the ejecta with a kinetic energy of 10^51 erg and a mass of 3 M⊙evolved in a uniform ambient medium for a time of-4000 yr and the jet component has cylindrical symmetry with a half open angle of 10°, the result indicates that the energy contained in the jet is about10%–15% of the kinetic energy of the entire ejecta to reproduce the detected profile. This study supports that the remnant originated from a jet-driven core-collapse supernova.

HYDRODYNAMIC ANALYSIS AND DESIGN METHOD OF OWC WAVE ENERGY CONVERTORS

In this paper,using potential flow theory and assuming that the pressure in the air cabin is proportional to the vertical velocity of the water column,we establish a hydrodynamic model for OWC wave energy convertors,resulting in a collection of illustrative plates,from which the behaviour of an OWC and the relation between its parameters are discussed. Optimization theory is used to design an optimum convertor working in irregular waves. It is found that the numerical results fit well the experimental data.

Hydrodynamics characterization of a choanoid fluidized bed bioreactor used in the bioartificial liver system： Fully resolved simulation with a fctitious domain method

Choanoid fluidized bed bioreactors （CFBBs） are newly developed core devices used in bioartificial liver- support systems to detoxify blood plasma of patients with microencapsulated liver cells. Direct numerical simulations （DNS） with a direct-forcing/fictitious domain （DF/FD） method were conducted to study the hydrodynamic performance of a CFBB. The effects of particle-fluid density ratio, particle number, and fil- ter screens preventing particles flowing out of the reactor were investigated. Depending on density ratio, two flow patterns are evident： the circulation mode in which the suspension rises along one sidewall and descends along the other sidewall, and the non-circulation mode in which the whole suspension roughly flows upward. The circulation mode takes place under non-neutral-buoyancy where the particle sedimentation dominates, whereas the non-circulation mode occurs under pure or near-neutral buoy- ancy with particle-fluid density ratios of unity or near unity. With particle-fluid density ratio of 1.01, the bioartificial liver reactor performs optimally as the significant particle accumulation existing in the non-circulation mode and the large shear forces on particles in the circulation mode are avoided. At higher particle volume fractions, more particles accumulate at the filter screens and a secondary counter circulation to the primary flow is observed at the top of the bed. Modelled as porous media, the filter screens play a negative role on particle fluidization velocities; without screens, particles are fluidized faster because of the higher fluid velocities in the jet center region. This work extends the DF/FD-based DNS to a fluidized bed and accounts for effects from inclined side walls and porous media, providing some hydrodynamics insight that is important for CFBB design and operation optimization.

Efficient turbulent compressible convection in the deep stellar atmosphere

We report on an application of gas-kinetic BGK scheme to the computation of turbulent compressible convection in the stellar interior. After incorporating the Sub-grid Scale （SGS） turbulence model into the BGK scheme, we tested the effects of numerical parameters on the quantitative relationships among the thermodynamic variables, their fluctuations and correlations in a very deep, initially gravity-stratified stellar atmosphere. Comparison indicates that the thermal properties and dynamic properties are dominated by different aspects of numerical models separately. An adjustable Deardorff constant in the SGS model cu, = 0.25 and an amplitude of artificial viscosity in the gas-kinetic BGK scheme C2 = 0 are appropriate for the current study. We also calculated the densityweighted auto- and cross-correlation functions in Xiong＇s turbulent stellar convection theory based on which the gradient type of models of the non-local transport and the anisotropy of the turbulence were preliminarily studied. No universal relations or constant parameters were found for these models.

Numerical Simulation of Interaction of a Current with a Circular Cylinder near a Rigid Bed

The numerical modeling of a 2D flow around a horizontal cylinder near a rigid bed with the gap ratios G/D = 0, 0.2, 0.3 at Reynolds numbers 840, 4500, 9000, and 9500 is investigated by using weakly compressible smoothed particle hydrodynamics. The velocity field and the separation angles from the present simulations are compared with those obtained from the experimental measurements and are in a good agreement. The results show that the maximum value of shear stress on the bed increases as the cylinder closes the bed and suddenly decreases when the cylinder contact the wall.

A One-Dimensional Relativistic Shock Model for the Light Curve of Gamma-ray Bursts

We investigate the forming of gamma-ray burst pulses with a simple onedimensional relativistic shock model. The mechanism is that a ＂central engine＂ drives forward the nearby plasma inside the fireball to generate a series of pressure waves. We give a relativistic geometric recurrence formula that connects the time when the pressure waves are produced and the time when the corresponding shocks occurred. This relation enables us to relate the pulse magnitude with the observation time. Our analysis shows that the evolution of the pressure waves leads to a fast rise and an exponential decay pulses. In determining the width of the pulses, the acceleration time is more important than that of the deceleration.

Numerical modeling of dam-break flood through intricate city layouts including underground spaces using GPU-based SPH method

This paper applies the meshfree Smoothed Particle Hydrodynamics （SPH） method with Graphical Processing Unit （GPU） parallel computing technique to investigate the highly complex 3-D dam-break flow in urban areas including underground spaces. Taking the advantage of GPUs parallel computing techniques, simulations involving more than 107 particles can be achieved. We use a virtual geometric plane boundary to handle the outermost solid wall in order to save considerable video card memory for the GPU computing. To evaluate the accuracy of the new GPU-based SPH model, qualitative and quantitative comparison to a real flooding experiment is performed and the results of a numerical model based on Shallow Water Equations （SWEs） is given with good accuracy. With the new GPU-based SPH model, the effects of the building layouts and underground spaces on the propagation of dambreak flood through an intricate city layout are examined.

TRACKING METHODS FOR FREE SURFACE AND SIMULATION OF A LIQUID DROPLET IMPACTING ON A SOLID SURFACE BASED ON SPH

With some popular tracking methods for free surface, simulations of several typical examples are carried out under various flow field conditions. The results show that the Smoothed Particle Hydrodynamics （SPH） method is very suitable in simulating the flow problems with a free surface. A viscous liquid droplet with an initial velocity impacting on a solid surface is simulated based on the SPH method, and the surface tension is considered by searching the free surface particles, the initial impact effect is considered by using the artificial viscosity method, boundary virtual particles and image virtual particles are introduced to deal with the boundary problem, and the boundary defect can be identified quite well. The comparisons of simulated results and experimental photographs show that the SPH method can not only exactly simulate the spreading process and the rebound process of a liquid droplet impacting on a solid surface but also accurately track the free surface particles, simulate the free-surface flow and determine the shape of the free surface due to its particle nature.

FLUID-STRUCTURE INTERACTION OF HYDRODYNAMIC DAMPER DURING THE RUSH INTO THE WATER CHANNEL

The hydrodynamic damper is a device to decrease the motion of armament carrier by use of the water resistance. When hydrodynamic damper rushes into the water channel with high velocity, it is a complicated flow phenomenon with fluid-structure interaction, free surface and moving interface. Numerical simulation using the Smoothed Particle Hydrodynamics （SPH） method coupled with the Finite Element （FE） method was successfully conducted to predict the dynamic characteristics of hydrodynamic damper. The water resistance, the pressure in the interface and the stress of structure were investigated, and the relationship among the peak of water resistance, initial velocity and actual draught was also discussed. The empirical formula was put forward to predict the water resistance. And it is found that the resistance coefficient is commonly in the range of 0.3 ≤ C ≤ 0.5, when the initial velocity is larger than 50 m/s. It can be seen that the SPH method coupled with the FE method has many obvious advantages over other numerical methods for this complicated flow problem with fluid-structure interaction.

Numerical simulation of submarine landslide tsunamis using particle based methods

This paper presents the simulation of tsunamis due to rigid and deformable landslides with consideration of submerged conditions by using particle methods. The smoothed particle hydrodynamics（SPH）, as a particle based method, is for solving problems of fast moving boundaries in the field of continuum mechanics. Other particle based methods, like the discrete element method（DEM）, are suitable for modeling the displacement and the collision related to the rigid landslides. In the present work, we use the SPH and the DEM to simulate tsunamis generated by rigid and deformable landslides with consideration of submerged conditions. The viscous free-surface flows are solved by a weakly compressible SPH and the displacement and the rotation of the rigid body slides are calculated using a multi-sphere DEM allowing for modeling solids of arbitrarily complex shapes. The fluid-solid interactions are simulated by coupling the SPH and the DEM. A rheology model combining the Papanastasiou and the Herschel-Bulkley models is applied to represent the viscoplastic behavior of the non-Newtonian flow in the submarine deformable landslide cases. Submarine landslide tsunamis due to rigid and deformable landslides are both simulated as typical landslide cases in this investigation. Our simulated results and the previous experimental results in the literatures are in good agreement, which shows that the proposed particle based methods are capable of modeling the submarine landslide tsunamis.

Simplified hydrodynamic models for the analysis of marine propellers in a wake-field

This paper presents a comparison among different hydrodynamic models for the analysis of the unsteady loads delivered by a marine propeller working in an axial, non-uniform inflow. Specifically, for a propeller subjected to a wake-field dominated by local high-frequency changes in space, the unsteady hydroloads predicted by the Nakatake formulation are compared with those given by the Theodorsen and Sears theories, respectively. Drawbacks and potentialities of these approaches are highlighted to assess a computationally efficient hydrodynamic solver for the analysis of operating conditions where propeller blades are significantly perturbed by a multi-harmonic onset-flow. Guidelines coming from this investigation may drive the choice of a fast and reliable unsteady propeller modeling that represents a good trade-off between accuracy of simulation and cost of computation within implementation in Computational Fluid Dynamics （CFD） solvers. The hydrodynamic formulations herein proposed are validated through numerical comparisons with the （accurate but computationally expensive） propeller loads predicted by a fully 3-D panel-method Boundary Element Method （BEM） solver, suited for the analysis of propellers operating in a complex hydrodynamic environment.

Sensitivity Analysis of the Hydrodynamic Coefficients in 4 Degrees of Freedom Ship Manoeuvring Mathematical Model

The S-type test is simulated based on a ship manoeuvring mathematical model of 4 degrees of freedom(4-DOF);simultaneously,sensitivity analysis of the hydrodynamic coefficients in the mathematical model is implemented by using an indirect method.The mathematical model is simplified by omitting the coefficients of smaller sensitivity according to the results of sensitivity analysis.The 10°/10° zigzag test and 35° turning circle manoeuvre are simulated with the original and the simplified mathematical models.The comparison of the simulation results shows the effectiveness of the sensitivity analysis and the validity of the simplified model.