Mathematical simulation of skidding

An approach to formulate a mathematical model of variable structure is proposed to describe the motion of a vehicle in different cases of loss of wheel contact with the road. This model is formed by a set of mathematical models of various dimensions and by some conditions of transition from one model to another. The effect of steering on the occurrence of a skid is studied.     

Numerical simulation of flow-induced bi-directional oscillations

Flow-induced vibration (FIV) by vortex shedding behind a submerged cylinder can lead to damage of nuclear components. With respect to such a serious scenario, various experiments and numerical simulations have been conducted to predict the vibration phenomena. Especially in simulation, the immersed finite element method (IFEM) is a promising approach to solve fluid-structure interaction problems because it needs less computational resources. In this paper, two-dimensional motions of cylinders are simulated by using IFEM to obtain their vibration characteristics. Three benchmark tests such as flow past a fixed circular cylinder, in-line oscillation of a circular cylinder and flow-induced vibration with uni-directional motion are performed to verify the proposed numerical method. Furthermore, bi-directional motions of two horizontally and vertically arranged cylinders as well as that of a single cylinder in fluid flow are analyzed, and then key findings are fully discussed.     

Coupled CFD-DEM simulation of hydrodynamic bridging at constrictions

This paper presents a coupled CFD-DEM approach to simulate the flow of particulate suspensions in the intermediate concentration regime where solid volume concentration is 1% < ϕ < 50%. In particular, hydrodynamic multi-particle bridging during flow through a single constriction in a rectangular channel is studied. It is shown that for neutrally buoyant, monodispersed particulate suspensions, the probability of jamming increases with the particle concentration. There also exists a critical particle concentration (ϕ*) for spontaneous bridging, which depends on the ratio of pore size to particle size, the flow velocity, the particle-fluid density contrast, and the flow geometry leading to the constriction. The ϕ* has a strong dependence on the outlet-to-particle relative size (R_o). For 1.5 ≤ R_o ≤ 2.5, a direct transition from a flowing state to a jammed state was observed. For R_o ≥ 3, the flowing state typically transitioned to a dense state characterized by the accumulation of particles near the constriction before jamming. Increasing the inlet-to-particle relative size (R_ip) lowers ϕ* by increasing the number of particles arriving at the constriction simultaneously. The effect of changing R_ip is more pronounced at high R_o when the probability of bridging is lower. A high fluid velocity increases particle interactions near the constriction and accelerates the onset of bridging. However, no distinct effect of velocity on ϕ* was observed in this study. A higher particle-to-fluid density ratio (ρ_p/ρ_f) increases the probability of bridging and leads to a lower ϕ* in a given constriction geometry. The effect saturates at higher ρ_p/ρ_fwhen gravitational forces completely dominate over viscous drag forces. ϕ* is also found to decrease with increasing angle of constriction convergence (θ) for θ < 30°, but increases beyond that at $\theta ={60}^{\circ }$.     

带自由表面水流与离散颗粒耦合模型研究

基于考察泥沙运动的细观行为特征,采用离散单元法(DEM)模拟泥沙颗粒运动,结合带自由表面的水动力学计算模型,建立了CFD-DEM耦合数值模型。计算程序开发基于Fortran语言来实现。耦合模型中实现了硬球模型和软球模型两种颗粒碰撞模型,应用范围较广。作为自由表面水流与泥沙颗粒流数值模型的初步研究,在模型建立的基础上,对模型做了基本的验证。分别通过单颗粒静水沉降和混合颗粒群分选两个计算工况,验证了模型的正确性及模拟精度。该耦合模型可进一步丰富带自由表面水流条件下泥沙运动的研究手段。     

Hardware-in-the-loop simulation of wave energy converters based on dielectric elastomer generators

Dielectric elastomer generators (DEGs) are soft electrostatic generators based on low-cost electroactive polymer materials. These devices have attracted the attention of the marine energy community as a promising solution to implement economically viable wave energy converters (WECs). This paper introduces a hardware-in-the-loop (HIL) simulation framework for a class of WECs that combines the concept of the oscillating water columns (OWCs) with the DEGs. The proposed HIL system replicates in a laboratory environment the realistic operating conditions of an OWC/DEG plant, while drastically reducing the experimental burden compared to wave tank or sea tests. The HIL simulator is driven by a closed-loop real-time hydrodynamic model that is based on a novel coupling criterion which allows rendering a realistic dynamic response for a diversity of scenarios, including large scale DEG plants, whose dimensions and topologies are largely different from those available in the HIL setup. A case study is also introduced, which simulates the application of DEGs on an OWC plant installed in a mild real sea laboratory test-site. Comparisons with available real sea-test data demonstrated the ability of the HIL setup to effectively replicate a realistic operating scenario. The insights gathered on the promising performance of the analysed OWC/DEG systems pave the way to pursue further sea trials in the future.

     

Mathematical simulation of biogenic processes in oil strata

An attempt is made to construct a mathematical model of the development of bacteria in an oil stratum and investigate on the basis of this model the main features of the process.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 84–91, October–December, 1981.     

Mathematical modeling and computer simulation of fire phenomena

An approach to the study of gas phase combustion and convection processes in fires using a combination of mathematical analysis and computer simulation is presented. It seeks to solve the governing equations directly (if approximately) by decomposing the fire into a large-scale convective and radiative transport problem coupled to a small-scale thermal-element model of combustion and radiative emission. The thermal-element model solves the combustion equations in a local Lagrangian coordinate system convected by the large-scale motion, which in turn is driven by the heat released by the combustion processes. The large-scale flow is studied using finite-difference techniques to solve large-eddy simulations of the Navier-Stokes equations. The basic theory behind the methodology is outlined and sample results of both large- and small-scale phenomena are presented. An analytical model of a large eddy is used to show how the simulation can be assembled to yield radiation feedback from a fire plume to a target surface.     

Mathematical study of the small oscillations of a pendulum containing a viscous liquid

The authors study the small oscillations of a pendulum containing a viscous liquid.Using the principle of virtual work, they determine the equations of motions and the variational formulation of the problem. Renouncing a little long method of the Russian scientists, they obtain directly from the variational equation, an operator equation in a Hilbert space.The corresponding bundle of operators belongs to a known class, as that the authors can study the spectrum of the problem and specify the set of eigenvalues according to the coefficient of viscosity.     

Dynamic simulation of the hydrodynamic interaction among immersed particles in Stokes flow

A numerical method for the dynamic simulation of the hydrodynamic interaction among particles in Stokes flow is developed. The method couples the quasi-static Stokes equations for the fluid with the equilibrium equations for the particles. The boundary element method is used to represent the velocity at a general field point in terms of surface velocities and stresses. However, neither the stresses nor the velocities are assumed to be known on the surface of the particles. Kinematic equations relating the linear and angular velocities at the centroids of the particles to the surface velocities are combined with the discretized boundary element equations and the equilibrium equations to generate a system of linear equations. The associated coefficient matrix is correspondent to the grand resistance matrix which relates the velocities of the particles to a given geometry.     

A Stokes model with cavitation for the numerical simulation of hydrodynamic lubrication

We present a cavitation model based on the Stokes equation and formulate adaptive finite element methods for its numerical solution. A posteriori error estimates and adaptive algorithms are derived, and numerical examples illustrating the theory are supplied, in particular with comparison to the simplified Reynolds model of lubrication. Copyright © 2010 John Wiley & Sons, Ltd.     

Mathematical simulation of avalanche-streamer transition in strong electric fields

Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 10–15, January–February, 1989.     

Numerical simulation of transonic limit cycle oscillations using high-order low-diffusion schemes

This paper simulates the NLR7301 airfoil limit cycle oscillation (LCO) caused by fluid–structure interaction (FSI) using Reynolds averaged Navier–Stokes equations (RANS) coupled with Spalart–Allmaras (S–A) one-equation turbulence model. A low diffusion E-CUSP (LDE) scheme with 5th order weighted essentially nonoscillatory scheme (WENO) is employed to calculate the inviscid fluxes. A fully conservative 4th order central differencing is used for the viscous terms. A fully coupled fluid–structural interaction model is employed. For the case computed in this paper, the predicted LCO frequency, amplitudes, averaged lift and moment, all agree excellently with the experiment performed by Schewe et al. The solutions appear to have bifurcation and are dependent on the initial fields or initial perturbation. The developed computational fluid dynamics (CFD)/computational structure dynamics (CSD) simulation is able to capture the LCO with very small amplitudes measured in the experiment. This is attributed to the high order low diffusion schemes, fully coupled FSI model, and the turbulence model used. This research appears to be the first time that a numerical simulation of LCO matches the experiment. The simulation confirms several observations of the experiment.     

Mathematical simulation of micellar-polymer displacement of oil from flooded strata

The micellar-polymer method of increasing the oil recovery from strata [1] is currently regarded as promising. The method consists of injecting into an oil stratum, which has previously undergone ordinary flooding, a relatively small amount, a slug, of micellar solution (5–10% of of the pore volume), which is propelled through the stratum by slugs of a highly viscous buffer fluid (aqueous solution of a polymer). In turn, the system of slugs is propelled from the injection points to the extraction wells by the water used for ordinary flooding. The displacement of the oil that remains after flooding in the stratum is achieved by a decrease in the coefficient of surface tension at the boundaries of the micellar solution with the oil and the water to the value 10^–2-10^–3 dyn/cm, which leads to a decrease in the amount of fixed oil and also to a control of the mobility of the fluids, which is achieved by varying the concentrations of the components of the micellar solution and the buffer fluid. The main components of micellar solutions are: a hydrocarbon fluid (oil or its fractions), water, surface-active substances. The relationships between the main components, and also the addition of salts and alcohol to the water component have a strong influence on the interaction between the solution and the stratal oil and water [2]. The micellar solution considered in the present paper dissolves oil but does not mix with water; the relationships between the components in it are characteristic of the solutions used to increase oil recovery from strata.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 84–93, November–December, 1982.     

Numerical simulation and hydrodynamic visualization of transient viscous flow around an oscillating aerofoil

Unsteady viscous flow around a large-amplitude and high-frequency oscillating aerofoil is examined in this paper by numerical simulation and experimental visualization. The numerical method is based on the combination of a fourth-order Hermitian finite difference scheme for the stream function equation and a classical second-order scheme to solve the vorticity transport equation. Experiments are carried out by a traditional visualization method using solid tracers suspended in water. The comparison between numerical and experimental results is found to be satisfactory. Time evolutions of the flow structure are presented for Reynolds numbers of 3 × 10³ and 10⁴. The influence of the amplitude and frequency of the oscillating motion on the dynamic stall is analysed.     

Predicting cortical oscillations with bidirectional LSTM network: a simulation study

It has been stated that up-down-state (UDS) cortical oscillation levels between excitatory and inhibitory neurons play a fundamental role in brain network construction. Predicting the time series behaviors of neurons in periodic and chaotic regimes can help in improving diseases, higher-order human activities, and memory consolidation. Predicting the time series is usually done by machine learning methods. In paper, the deep bidirectional long short-term memory (DBLSTM) network is employed to predict the time evolution of regular, large-scale UDS oscillations produced by a previously developed neocortical network model. In noisy time-series prediction tasks, we compared the DBLSTM performance with two other variants of deep LSTM networks: standard LSTM, LSTM projected, and gated recurrent unit (GRU) cells. We also applied the classic seasonal autoregressive integrated moving average (SARIMA) time-series prediction method as an additional baseline. The results are justified through qualitative resemblance between the bifurcation diagrams of the actual and predicted outputs and quantitative error analyses of the network performance. The results of extensive simulations showed that the DBLSTM network provides accurate short and long-term predictions in both periodic and chaotic behavioral regimes and offers robust solutions in the presence of the corruption process.

     

Mathematical simulation of unsteady flow past a wing with jets

The various approximate approaches to the investigation of the unsteady aerodynamic characteristics of an airfoil with jet flap [1–3] are applicable only for an airfoil, low jet intensity, and low oscillation frequencies. In the present paper, the method of discrete vortices [4] is generalized to the case of unsteady flow past a wing with jets and arbitrary shape in plan. The problem is solved in the linear formulation; the conditions used are standard: no flow through the wing and jet, finite velocities at the trailing edges where there is no jet, and also a dynamical condition on the jet. The wing and jet are assumed to be thin and the medium inviscid and incompressible.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 139–144, May–June, 1982.     

Mathematical simulation of mixed convection in a vertical cylindrical vessel

Mixed laminar convection in a vertical cylindrical vessel is investigated. A flow model is developed and the temperature distribution in the core is obtained for a given specific heat flux on the lateral surface. The results of the analytic simulation are compared with numerical solutions of the problem for a specific heat flux which is constant on the lateral wall or varies linearly with height. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 9–17, November–December, 1998.     

A numerical simulation of hydrodynamic forces of ground-effect problem using Lagrange's equation of motion

On the basis of the potential flow theory, Lagrange's equation of motion is used to study the unsteady ground-effect problem. The forces and moments acting on the moving body are solved in terms of the derivatives of added masses in which the generalized Taylor's formulae are applied. The singular integral equations used to solve the surface source intensities and their derivatives are regularized by the Gauss flux theorem and are therefore amenable to the direct use of the Gaussian quadrature formula. In illustration, the condition of a prolate spheroid moving in the fore-and-aft direction at constant speed past a flat ground with a protrusion is considered. The hydrodynamic forces and moments acting on the moving spheroid are investigated systematically by varying the size of the protrusion and the cruising height of the spheroid. © 1998 John Wiley & Sons, Ltd.     

Numerical simulation of plunging wave breaking by the weakly compressible smoothed particle hydrodynamic method

In this paper, a numerical meshless method called the weakly compressible smoothed particle hydrodynamic (WCSPH) method, which is a two-dimensionalmodel of a weakly compressible fluid, is applied to simulate the plunging wave breaking process. This model solves the viscous fluid equations to obtain the velocity and density fields and also solves the equation of state to obtain the pressure field. The WCSPH method is demonstrated to have a higher computational efficiency than the basic SPH model. To simulate the turbulent behavior of the fluid flow in the wave breaking procedure, a sub particle scale (SPS) model is used, which is obtained from the Large eddy simulation (LES) theory. To consider the accuracy of the standard WCSPH model (WCSPH model without considering the turbulent effect), a dam break test is performed, and model results are compared with available experimental data.