Time domain numerical modeling of wave propagation in 2D heterogeneous porous media

This paper deals with the numerical modeling of wave propagation in porous media described by Biot’s theory. The viscous efforts between the fluid and the elastic skeleton are assumed to be a linear function of the relative velocity, which is valid in the low-frequency range. The coexistence of propagating fast compressional wave and shear wave, and of a diffusive slow compressional wave, makes numerical modeling tricky. To avoid restrictions on the time step, the Biot’s system is splitted into two parts: the propagative part is discretized by a fourth-order ADER scheme, while the diffusive part is solved analytically. Near the material interfaces, a space–time mesh refinement is implemented to capture the small spatial scales related to the slow compressional wave. The jump conditions along the interfaces are discretized by an immersed interface method. Numerical experiments and comparisons with exact solutions confirm the accuracy of the numerical modeling. The efficiency of the approach is illustrated by simulations of multiple scattering.     

Numerical solution of an extended White–Metzner model for eccentric Taylor–Couette flow

In this study, we have developed a new numerical approach to solve differential-type viscoelastic fluid models for a commonly used benchmark problem, namely, the steady Taylor—Couette flow between eccentric cylinders. The proposed numerical approach is special in that the nonlinear system of discretized algebraic flow equations is solved iteratively using a Newton–Krylov method along with an inverse-based incomplete lower-upper preconditioner. The numerical approach has been validated by solving the benchmark problem for the upper-convected Maxwell model at a large Deborah number. Excellent agreement with the numerical data reported in the literature has been found. In addition, a parameter study was performed for an extended White–Metzner model. A large eccentricity ratio was chosen for the cylinder system in order to allow flow recirculation to occur. We detected several interesting phenomena caused by the large eccentricity ratio of the cylinder system and by the viscoelastic nature of the fluid. Encouraged by the results of this study, we intend to investigate other polymeric fluids having a more complex microstructure in an eccentric annular flow field.     

对CELS算法在多重网格下的数值试验

求解不可压流体耦合方程的直接解法（CELS算法）在许多计算问题中体现了其优越性．多重网格法在加速数值计算的收敛速度方面也体现出了其有效性．本文把CELS算法在多重网格下实施,使用原始变量法,对圆管突扩通道和外掠后台阶的层流流动进行了计算,并与Armaly等人的实验数据进行了比较．结果表明,这种结合是很有效的．     

采用非交错网格压力修正法计算不可压及可压亚、跨、超音速流动

采用非交错网格压力修正法计算不可压及可压亚、跨、超音速流动朱斌，蔡睿贤，蒋洪德（中国科学院工程热物理研究所北京１０００８０）关键词：压力修正法，亚、跨、超音速流动一、前言源于计算不可压流动的压力修正方法，由于具有计算时间短及不受马赫数限制的优点，近年...     

Modeling of a rotor speed transient response with radial rubbing

A rotor-stator model of a turbogenerator is introduced in order to investigate speed transients with rotor-to-stator rubbing caused by an accidental blade-off imbalance. In order to assess the angular deceleration of the rotor due to rubbing, the angular position of its cross-section is considered as an unknown of the problem. Displacement fields are discretized through a finite element formulation. The highly nonlinear equations due to contact conditions are solved through an explicit prediction-correction time-marching procedure combined with the Lagrange multiplier approach dealing with a node-to-line contact strategy. The developed numerical tool is suitable for analyzing rotor-stator interactions in turbomachines as the system passes through critical speeds during an accidental shutdown. The sensitivity of the system response to modeling, physical and numerical parameters is investigated. The results highlight the significant role of the friction coefficient together with the diaphragm modeling, from rigid to fully flexible, in the interaction phenomenon. Rigid models have the advantage of simplicity and provide reasonable estimations of the overall response of the turbine. A flexible model, however, may be more computationally intensive but is more appropriate in order to accurately capture quantities of interest such as shaft eccentricity and bearing loads.     

A composite grid solver for conjugate heat transfer in fluid–structure systems

We describe a numerical method for modeling temperature-dependent fluid flow coupled to heat transfer in solids. This approach to conjugate heat transfer can be used to compute transient and steady state solutions to a wide range of fluid–solid systems in complex two- and three-dimensional geometry. Fluids are modeled with the temperature-dependent incompressible Navier–Stokes equations using the Boussinesq approximation. Solids with heat transfer are modeled with the heat equation. Appropriate interface equations are applied to couple the solutions across different domains. The computational region is divided into a number of sub-domains corresponding to fluid domains and solid domains. There may be multiple fluid domains and multiple solid domains. Each fluid or solid sub-domain is discretized with an overlapping grid. The entire region is associated with a composite grid which is the union of the overlapping grids for the sub-domains. A different physics solver (fluid solver or solid solver) is associated with each sub-domain. A higher-level multi-domain solver manages the entire solution process.     

3D phase-field simulations of interfacial dynamics in Newtonian and viscoelastic fluids

This work presents a three-dimensional finite-element algorithm, based on the phase-field model, for computing interfacial flows of Newtonian and complex fluids. A 3D adaptive meshing scheme produces fine grid covering the interface and coarse mesh in the bulk. It is key to accurate resolution of the interface at manageable computational costs. The coupled Navier–Stokes and Cahn–Hilliard equations, plus the constitutive equation for non-Newtonian fluids, are solved using second-order implicit time stepping. Within each time step, Newton iteration is used to handle the nonlinearity, and the linear algebraic system is solved by preconditioned Krylov methods. The phase-field model, with a physically diffuse interface, affords the method several advantages in computing interfacial dynamics. One is the ease in simulating topological changes such as interfacial rupture and coalescence. Another is the capability of computing contact line motion without invoking ad hoc slip conditions. As validation of the 3D numerical scheme, we have computed drop deformation in an elongational flow, relaxation of a deformed drop to the spherical shape, and drop spreading on a partially wetting substrate. The results are compared with numerical and experimental results in the literature as well as our own axisymmetric computations where appropriate. Excellent agreement is achieved provided that the 3D interface is adequately resolved by using a sufficiently thin diffuse interface and refined grid. Since our model involves several coupled partial differential equations and we use a fully implicit scheme, the matrix inversion requires a large memory. This puts a limit on the scale of problems that can be simulated in 3D, especially for viscoelastic fluids.     

Metachronal propulsion of non-Newtonian viscoelastic mucus in an axisymmetric tube with ciliated walls

Cilia-induced flow of viscoelastic mucus through an idealized two-dimensional model of the human trachea is presented.The cilia motion is simulated by a metachronal wave pattern which enables the mobilization of highly viscous mucus even at nonzero Reynolds numbers.The viscoelastic mucus is analyzed with the upper convected Maxwell viscoelastic formulation which features a relaxation time and accurately captures normal stress generation in shear flows.The governing equations are transformed from fixed to wave(laboratory)frame with appropriate variables and resulting differential equations are perturbed about wave number.The trachea is treated as an axisymmetric ciliated tube.Radial and axial distributions in axial velocity are calculated via the regular perturbation method and pressure rise is computed with numerical integration using symbolic software MATHEMATICA^‘TM’.The influence of selected parameters which is cilia length,and Maxwell viscoelastic material parameter i.e.relaxation time for prescribed values of wave number are visualized graphically.Pressure rise is observed to increase considerably with elevation in both cilia length and relaxation time whereas the axial velocity is markedly decelerated.The simulations provide some insight into viscous-dominated cilia propulsion of rheological mucus and also serve as a benchmark for more advanced modeling.     

Numerical investigation of electro hydrodynamics (EHD) enhanced water evaporation using Large Eddy Simulation turbulent model

In this paper, the electric field effect on enhancement of the water evaporation rate in a channel is numerically investigated. The coupled equations of electrical field, flow field, temperature field and species concentration fields are discretized using Finite Volume Method (FVM) and the implicit/hybrid difference form. This equations are solved via SIMPLE algorithm and the Kaptsov hypothesis. The turbulent flow is modelled by Large Eddy Simulation (LES). The numerical results show that the water evaporation rate is increased with the presence of electric field, but the effect of electric field is diminished at high Reynolds number.     

AN EXACT METHOD FOR THE FREE VIBRATION ANALYSIS OF TIMOSHENKO-KELVIN BEAMS WITH OSCILLATORS

In this paper, a modern exact method is proposed for solving the problem of free vibrations of a Timoshenko-type viscoelastic beam with discrete rigid bodies, connected to the beam by means of viscoelastic constraints. The phenomenon of free vibrations of this discrete-continuous system is described by a set of three partial and two subsystem ordinary differential equations with generalized boundary conditions and initial conditions. Vector notation of the equations allows one to identify the self-adjoint linear operators of inertia, stiffness and damping. In this case, these operators are not homothetic hence a separation of variables in this set of equations is possible only in a complex Hilbert space. Such separation of variables leads to ordinary differential equations of motion with respect to time as well as to a set of three ordinary differential equations with respect to a spatial variable and two subsystem algebraical equations. Solution of the boundary-value problem was carried out in the classical way, but its results are complex conjugated. Using these results and the fundamental principle, describing the orthogonality property of complex eigenvectors, the problem of free vibrations of the system with arbitrary initial conditions has been finally solved exactly.     

Numerical simulation of Gaussian beam scattering by complex particles of arbitrary shape and structure

An efficient numerical method based on the surface integral equations is introduced to simulate the scattering of Gaussian beam by complex particles that consist of an arbitrarily shaped host particle and multiple internal inclusions of arbitrary shape. In particular, the incident focused Gaussian beam is described by the Davis fifth-order approximate expressions in combination with rotation defined by Euler angles. The established surface integral equations are discretized with the method of moments, where the unknown equivalent electric and magnetic currents induced on the surfaces of the host particle and the internal inclusions are expanded using the Rao–Wilton–Glisson (RWG) basis functions. The resultant matrix equations are solved by using the parallel conjugate gradient method. The proposed numerical method is validated and its capability illustrated in several characteristic examples.     

Two-dimensional oblique stagnation-point flow towards a stretching surface in a viscoelastic fluid

In this paper, the steady two-dimensional stagnation-point flow of a viscoelastic Walters’ B’ fluid over a stretching surface is examined. It is assumed that the fluid impinges on the wall obliquely. Using similarity variables, the governing partial differential equations are transformed into a set of two non-dimensional ordinary differential equations. These equations are then solved numerically using the shooting method with a finite-difference technique.     

Solution of Hartree-Fock equations in coordinate space for axially symmetric nuclei

As an alternative to the conventional deformed harmonic oscillator basis expansion, a method is developed in which the coupled integro-differential Hartree-Fock equations are solved directly in coordinate space for a simple effective interaction. The single particle wave functions are obtained on a finite mesh by minimizing a discretized energy functional and solving the resulting finite difference equations using the Lanczos algorithm. Expressions to correct the total energy to second order in the mesh spacing are derived, and the accuracy of the method is demonstrated by numerical comparison with spherical results. Applications and advantages of this new technique are briefly discussed.     

Nonlinear dynamics of higher-dimensional system for an axially accelerating viscoelastic beam with in-plane and out-of-plane vibrations

In this paper, the bifurcations and chaotic motions of higher-dimensional nonlinear systems are investigated for the nonplanar nonlinear vibrations of an axially accelerating moving viscoelastic beam. The Kelvin viscoelastic model is chosen to describe the viscoelastic property of the beam material. Firstly, the nonlinear governing equations of nonplanar motion for an axially accelerating moving viscoelastic beam are established by using the generalized Hamilton’s principle for the first time. Then, based on the Galerkin’s discretization, the governing equations of nonplanar motion are simplified to a six-degree-of-freedom nonlinear system and a three-degree-of-freedom nonlinear system with parametric excitation, respectively. At last, numerical simulations, including the Poincare map, phase portrait and Lyapunov exponents are used to analyze the complex nonlinear dynamic behaviors of the axially accelerating moving viscoelastic beam. The bifurcation diagrams for the in-plane and out-of-plane displacements via the mean axial velocity, the amplitude of velocity fluctuation and the frequency of velocity fluctuation are respectively presented when other parameters are fixed. The Lyapunov exponents are calculated to identify the existence of the chaotic motions. From the numerical results, it is indicated that the periodic, quasi-periodic and chaotic motions occur for the nonplanar nonlinear vibrations of the axially accelerating moving viscoelastic beam. Observing the in-plane nonlinear vibrations of the axially accelerating moving viscoelastic beam from the numerical results, it is found that the nonlinear responses of the six-degree-of-freedom nonlinear system are much different from that of the three-degree-of-freedom nonlinear system when all parameters are same.     

Comparison of the generalized Riemann solver and the gas-kinetic scheme for inviscid compressible flow simulations

The generalized Riemann problem (GRP) scheme for the Euler equations and gas-kinetic scheme (GKS) for the Boltzmann equation are two high resolution shock capturing schemes for fluid simulations. The difference is that one is based on the characteristics of the inviscid Euler equations and their wave interactions, and the other is based on the particle transport and collisions. The similarity between them is that both methods can use identical MUSCL-type initial reconstructions around a cell interface, and the spatial slopes on both sides of a cell interface involve in the gas evolution process and the construction of a time-dependent flux function. Although both methods have been applied successfully to the inviscid compressible flow computations, their performances have never been compared. Since both methods use the same initial reconstruction, any difference is solely coming from different underlying mechanism in their flux evaluation. Therefore, such a comparison is important to help us to understand the correspondence between physical modeling and numerical performances. Since GRP is so faithfully solving the inviscid Euler equations, the comparison can be also used to show the validity of solving the Euler equations itself. The numerical comparison shows that the GRP exhibits a slightly better computational efficiency, and has comparable accuracy with GKS for the Euler solutions in 1D case, but the GKS is more robust than GRP. For the 2D high Mach number flow simulations, the GKS is absent from the shock instability and converges to the steady state solutions faster than the GRP. The GRP has carbuncle phenomena, likes a cloud hanging over exact Riemann solvers. The GRP and GKS use different physical processes to describe the flow motion starting from a discontinuity. One is based on the assumption of equilibrium state with infinite number of particle collisions, and the other starts from the non-equilibrium free transport process to evolve into an equilibrium one through particle collisions. The different mechanism in the flux evaluation deviates their numerical performance. Through this study, we may conclude scientifically that it may NOT be valid to use the Euler equations as governing equations to construct numerical fluxes in a discretized space with limited cell resolution. To adapt the Navier–Stokes (NS) equations is NOT valid either because the NS equations describe the flow behavior on the hydrodynamic scale and have no any corresponding physics starting from a discontinuity. This fact alludes to the consistency of the Euler and Navier–Stokes equations with the continuum assumption and the necessity of a direct modeling of the physical process in the discretized space in the construction of numerical scheme when modeling very high Mach number flows. The development of numerical algorithm is similar to the modeling process in deriving the governing equations, but the control volume here cannot be shrunk to zero.     

Algebraic Splitting for Incompressible Navier–Stokes Equations

Fully discretized incompressible Navier–Stokes equations are solved by splitting the algebraic system with an approximate factorization. This splitting affects the temporal convergence order of velocity and pressure. The splitting error is proportional to the pressure variable, and a simple analysis shows that the original convergence order of the time-integration scheme can be retained by solving for incremental pressure. The combination of splitting and incremental pressure is shown to be equivalent to an error-correcting method using the full pressure. In numerical experiments employing a third-order time-integration scheme and various orders for the pressure increment, the splitting error is shown to control the convergence order, and the full order of the scheme is recaptured for both velocity and pressure. The difference between perturbing the momentum or the continuity equation is also explored.     

Four Spacetime Dimensional Simulation of Rheological Waves in Solids and the Merits of Thermodynamics

The recent results attained from a thermodynamically conceived numerical scheme applied on wave propagation in viscoelastic/rheological solids are generalized here, both in the sense that the scheme is extended to four spacetime dimensions and in the aspect of the virtues of a thermodynamical approach. Regarding the scheme, the arrangement of which quantity is represented where in discretized spacetime, including the question of appropriately realizing the boundary conditions, is nontrivial. In parallel, placing the problem in the thermodynamical framework proves to be beneficial in regards to monitoring and controlling numerical artefacts—instability, dissipation error, and dispersion error. This, in addition to the observed preciseness, speed, and resource-friendliness, makes the thermodynamically extended symplectic approach that is presented here advantageous above commercial finite element software solutions.     

A meshless method based on moving Kriging interpolation for a two-dimensional time-fractional diffusion equation

Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.     

Nonlinear mechanical response of supercooled melts under applied forces

We review recent progress on a microscopic theoretical approach to describe the nonlinear response of glass-forming colloidal dispersions under strong external forcing leading to homogeneous and inhomogeneous flow. Using mode-coupling theory (MCT), constitutive equations for the rheology of viscoelastic shear-thinning fluids are obtained. These are, in suitably simplified form, employed in continuum fluid dynamics, solved by a hybrid-Lattice Boltzmann (LB) algorithm that was developed to deal with long-lasting memory effects. The combined microscopic theoretical and mesoscopic numerical approach captures a number of phenomena far from equilibrium, including the yielding of metastable states, process-dependent mechanical properties, and inhomogeneous pressure-driven channel flow.     

Direct pore-level modeling of incompressible fluid flow in porous media

We present a dynamic particle-based model for direct pore-level modeling of incompressible viscous fluid flow in disordered porous media. The model is capable of simulating flow directly in three-dimensional high-resolution micro-CT images of rock samples. It is based on moving particle semi-implicit (MPS) method. We modify this technique in order to improve its stability for flow in porous media problems. Using the micro-CT image of a rock sample, the entire medium, i.e., solid and fluid, is discretized into particles. The incompressible Navier–Stokes equations are then solved for each particle using the MPS summations. The model handles highly irregular fluid–solid boundaries effectively. An algorithm to split and merge fluid particles is also introduced. To handle the computational load, we present a parallel version of the model that runs on distributed memory computer clusters. The accuracy of the model is validated against the analytical, numerical, and experimental data available in the literature. The validated model is then used to simulate both unsteady- and steady-state flow of an incompressible fluid directly in a representative elementary volume (REV) size micro-CT image of a naturally-occurring sandstone with 3.398 μm resolution. We analyze the quality and consistency of the predicted flow behavior and calculate absolute permeability using the steady-state flow rate.