2D numerical simulation of submerged hydraulic jumps

Hydraulic jumps are usually used to dissipate energy in hydraulic engineering. In this paper, the turbulent submerged hydraulic jumps are simulated by solving the unsteady Reynolds averaged Navier–Stokes equations along with the continuity equation and the standard k–? equations for turbulence modeling. The Lagrangian moving grid method is employed for the simulation of the free surface. In the developed model, kinematic free-surface boundary condition is solved simultaneously with the momentum and continuity equations, so that the water elevation can be obtained along with velocity and pressure fields as part of the solution. Computational results are presented for Froude numbers ranging from 3.2 to 8.2 and submergence factors ranging from 0.24 to 0.85. Comparisons with experimental measurements show that numerical model can simulate the velocity field, variation of free surface, maximum velocity, Reynolds shear and normal stresses at various stations with reasonable accuracy.     

3D numerical simulation on fluid-structure interaction of structure subjected to underwater explosion with cavitation

In the underwater-shock environment,cavitation occurs near the structural surface.The dynamic response of fluid-structure interactions is influenced seriously by the cavitation effects.It is also the d...     

Coupling strategies for the numerical simulation of blood flow in deformable arteries by 3D and 1D models

The fluid structure interaction mechanism in vascular dynamics can be described by either 3D or 1D models, depending on the level of detail of the flow and pressure patterns needed for analysis. A successful strategy that has been proposed in the past years is the so-called geometrical multiscale approach, which consists of coupling both 3D and 1D models so as to use the former only in those regions where details of the fluid flow are needed and describe the remaining part of the vascular network by the simplified 1D model.In this paper we review recently proposed strategies to couple the 3D and 1D models, and within the 3D model, to couple the fluid and structure sub-problems. The 3D/1D coupling strategy relies on the imposition of the continuity of flow rate and total normal stress at the interface. On the other hand, the fluid–structure coupling strategy employs Robin transmission conditions. We present some numerical results and show the effectiveness of the new approaches.     

Non-Boussinesq turbulent buoyant jet of a low-density gas leaks into high-density ambient

In this article, we study the problem of low-density gas jet injected into high-density ambient numerically which is important in applications such as fuel injection and leaks. It is assumed that the local rate of entrainment is consisted of two components; one is the component of entrainment due to jet momentum while the other is the component of entrainment due to buoyancy. The integral models of the mass, momentum and concentration fluxes are obtained and transformed to a set of ordinary differential equations using some similarity transformations. The resulting system is solved to determine the centerline quantities which are used to get the mean axial velocity, mean concentration and mean density of the jet. Therefore, the centerline and mean quantities are used together with the governing equation to determine some important turbulent quantities such as, cross-stream velocity, Reynolds stress, velocity-concentration correlation, turbulent eddy viscosity and turbulent eddy diffusivity. Throughout this paper the developed model is verified by comparing the present results with experimental results and jet/plume theory from the literature.     

3D mathematical model and numerical simulation for laying marine cable along prescribed trajectory on seabed

In this paper, we establish a mathematical model, that is, an initial-boundary value problem of partial differential equations (PDE), to describe the three dimensional (3D) motion of a marine cable being laid onto the seabed of varying depth. Based on this PDE model, a numerical method for simulating dynamically the moving process of marine cable is developed, which can be used to determine the cable laying operations in practice such that the marine cable can be laid on the prescribed position on seabed.     

Regularized shallow water equations and an efficient method for numerical simulation of shallow water flows

Regularized shallow water equations are derived as based on a regularization of the Navier-Stokes equations in the form of quasi-gasdynamic and quasi-hydrodynamic equations. Efficient finite-difference algorithms based on the regularized shallow water equations are proposed for the numerical simulation of shallow water flows. The capabilities of the model are examined by computing a test Riemann problem, the flow over an obstacle, and asymmetric dam break.     

Regularized shallow water equations for numerical simulation of flows with a moving shoreline

A numerical algorithm for simulating free-surface flows based on regularized shallow water equations is adapted to flows involving moving dry-bed areas. Well-balanced versions of the algorithm are constructed. Test computations of flows with dry-bed areas in the cases of water runup onto a plane beach and a constant-slope beach are presented. An example of tsunami simulation is given.     

ALE-based 3D FE simulation of electromagnetic forming

Electromagnetic forming is a contact-free high-speed forming process. The deformation of the work piece is driven by the Lorentz force which results from the interaction of a pulsed magnetic field with eddy currents induced in the work piece by the field itself. The purpose of this work is to present a fully-coupled three-dimensional simulation of this process. For the mechanical structure, a thermoelastic, viscoplastic, electromagnetic material model is relevant, which is incorporated in a large-deformation dynamic formulation. The electromagnetic fields are governed by Maxwell's equations under quasistatic conditions. To consider their reduced regularity at material interfaces Nédélec elements are applied. Coupling takes the form of the Lorentz force, the electromotive intensity and the current geometry of the work piece. A staggered solution scheme based on a Lagrangian mesh for the work piece and an ALE formulation for the electromagnetic field is employed. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Coupled numerical simulation and sensitivity assessment for quality modelling for water distribution systems

Sensitivity analysis is an important tool which can be used to investigate the stability of a process perturbed by parameter changes and uncertainty impacts. In this work the unsteady sensitivity equations for complex looped pipe networks are solved. Special attention is focused on the coupled version of these equations, with the direct problem. For this purpose a splitting method using a Total Variation Diminishing (TVD) scheme with very good quality of stability is set up and validated on a benchmark pipe network.     

A numerical study of singular stress field of 3D cracks

We investigate the stress distribution and the variation of the mode I stress intensity factor along a straight three-dimensional (3D) crack by the finite element method. The results are checked against plane strain theory near the mid-crack and against the 3D theory of Zhu at the free surface. Although Zhu's formulation is not perfect and has some typographical errors. The surface stress distribution of his results are in line with the present study by the finite element method. The stress intensity factors at the free surface are found to be much lower than that at the mid-crack.     

Temporally regularized direct numerical simulation

Experience with fluid-flow simulation suggests that, in some instances, under-resolved direct numerical simulation (DNS), without a residual-stress model per se but with artificial damping of small scales to account for energy lost in the cascade from resolved to unresolved scales, may be as reliable as simulations based on more complex models of turbulence. One efficient and versatile manner to selectively damp under-resolved spatial scales is by a relaxation regularization, e.g. Stolz and Adams [S. Stolz, N.A. Adams, An approximate deconvolution procedure for large eddy simulation, Phys. Fluids II (1999) 1699-1701]. We consider the analogous approach based on time scales, time filtering and damping of under-resolved temporal features. The paper explores theoretical and practical aspects of temporally damped fluid-flow simulations. We prove existence of solutions to the resulting continuum model. We also establish the effect of the damping of under-resolved temporal features as the energy balance and dissipation and prove that the time fluctuations → 0 in a precise sense. The method is then demonstrated to obtain both steady-state and time-dependent coarse-grid solutions of the Navier-Stokes equations.     

On numerical simulation of tsunami run-up

The improved adaptive grid method is presented for the numerical modelling of propagation of long surface waves in large water areas and their interaction with coasts. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The skim of balance theory of 3D garment simulation

The balance theory of 3D garment simulation tries to integrally study the complex system of human body-fabric-environment. The theory includes the geometrical balance, dynamic balance, pressure balance, and heat and moisture transfer balance. In this paper, the ICAD-balance equations are established based on the framework of balance theory and the heat and moisture transfer model of garment. A finite volume method is employed to solve the equations. Experiment results show that ICAD-balance equations can excellently simulate real situation. Therefore, our balance theory provides a good tool for customized clothing design, electronic mirror for clothing-test-platform, and comfortable clothing evaluation research.     

High-resolution numerical simulation of 2D nonlinear wave structures in electromagnetic fluids with absorbing boundary conditions

Here we show how the full set of governing equations for the dynamics of charged-particle fluids in an electromagnetic field may be solved numerically in order to model nonlinear wave structures propagating in two dimensions. We employ a source-term adaptation and two-fluid extension of the second-order high-resolution central scheme of Balbas et al. (2004) [1]. The model employed is a 2D extension of that used by Baboolal and Bharuthram (2007) [5] in studies of 1D shocks and solitons in a two-fluid plasma under 3D electromagnetic fields. Further, we outline the use of free-flow boundary conditions to obtain stable wave structures over sufficiently long modelling times. As illustrative results, we examine the formation and evolution of shock-like and soliton structures of the magnetosonic mode.     

电离层层析成像及数值模拟

利用正则化反演法,对电离层层析成像技术进行了数值模拟.模拟结果表明,正则化反演法能够准确反演出电离层赤道异常的电子密度分布.对正则化反演误差的原因进行了分析,定义了代表重建结果精度的参数并研究其随信号数量的变化,结果显示造成反演误差的主要原因是信息量不足,在提供充足的数据条件下,正则化反演法能够得到很高精度的重建结果.数值模拟实验进一步证明了电离层层析技术在反演电离层电子密度分布中有很大的应用价值.     

稠油开采数值模拟研究

针对稠油油藏使用水平井并且注蒸汽的开采方式,建立了相应的数学模型,模型考虑了油、气、水三个相态及所有组分,和由于稠油流体特性而引起的启动压力,对数学模型使用九点差分进行离散,用全隐式方法将数值模型线性化,给出了模型等效化的处理方法,将计算结果与模拟区块的实际生产数据做动态历史拟合,并提出了开发调整方案,取得了良好效果.该模型与其它方法做了比较,模型所得结果与实际产量的绝对值偏差为1.9％,是所有方法中最小的.模型与求解方法适用于稠油油藏注蒸汽水平并开采数值模拟.