Analytical and numerical solutions of the problem of gravity-induced mixing of a light layer using the <Emphasis Type="Italic">k</Emphasis>–<Emphasis Type="Italic">ε</Emphasis> turbulence model

An analytical solution for the self-similar stage in the problem of gravity-induced turbulent mixing in a light (heavy) layer is obtained on the basis of the k–ε model equations. The solution obtained is compared with the results of a numerical investigation of the problem using both three-dimensional direct numerical simulation and the k–ε model. The calculations were performed using the two- and three-dimensional versions of the EGAK method. The results of all the calculations and the available experimental data are in reasonable agreement.     

Numerical simulations of shock wave propagation in condensed multiphase materials

We propose to find out numerical solutions of a travelling shock wave in condensed mixtures by using a direct numerical simulation. Condensed multiphase materials under shock wave conditions are mechanically characterized by a unique pressure and a unique velocity. In this study, the mixture is considered as a collection of grains separated by interface between each material: this problem of interfaces is solved by a diffuse interface method. The results will be compared to existing one-dimensional numerical models, analytical solutions and also experimental data. The volume fraction (or the phase temperature) is not measured in experiments and it is then important to verify the behaviour of a phase quantity through various methods. A non-monotonous evolution of the volume fraction is obtained with analytical solution as well as numerical simulation.      

Construction of an approximate solution of the problem of a cavitational flow around a disk

An approach to analytical analysis of a steady axisymmetric cavitation flow of a potential ideal incompressible fluid is developed. An approximate analytical solution of Ryabushinskii’s problem for a disk-shaped cavitator is constructed. All integral and local characteristics of the flow calculated by the proposed approximation formulas in a wide range of cavitation numbers are demonstrated to differ insignificantly from the numerical solution of the problem in the exact nonlinear formulation.     

Identification of complex flows in Taylor–Couette counter-rotating cavities

The transition in confined rotating flows is a topical problem with many industrial and fundamental applications. The purpose of this study is to investigate the Taylor–Couette flow in a finite-length cavity with counter-rotating walls, for two aspect ratios L=5 or L=6. Two complex regimes of wavy vortex and spirals are emphasized for the first time via direct numerical simulation, by using a three-dimensional spectral method. The spatio-temporal behavior of the solutions is analyzed and compared to the few data actually available.     

Analytical and numerical study of double diffusive convection in a vertical enclosure

This paper presents an analytical and numerical study of natural convection of a double-diffusive fluid contained in a rectangular slot subject to uniform heat and mass fluxes along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number, Ra _T ; buoyancy ratio, N; Lewis number, Le; Prandtl number, Pr and aspect ratio of the cavity, A. In the first part of the analytical study a scale analysis is applied to the two extreme cases of heat-transfer and mass-transfer-driven flows. In the second part, an analytical solution, based on the parallel flow approximation, is reported for tall enclosures (A?1). Solutions for the flow fields, temperature and concentration distributions and Nusselt and Sherwood numbers are obtained in terms of the governing parameters of the problem. In the limits of heat-driven and solute-driven flows a good agreement is obtained between the prediction of the scale analysis and those of the analytical solution. The numerical solutions are based on the complete governing equations for two-dimensional flows, and cover the range 1≤Ra _T ≤10⁷, 0≤N≤10⁵, 10^-3≤Le≤10³, 1≤A≤20 and Pr=7. A good agreement is found between the analytical predictions and the numerical simulation.     

Enhanced equilibrium distribution functions for simulating immiscible multiphase flows with variable density ratios in a class of lattice Boltzmann models

This research examines the behavior of a class of lattice Boltzmann (LB) models designed to simulate immiscible multiphase flows. There is some debate in the scientific literature as to whether or not the “color gradient” models, also known as the Rothman–Keller (RK) models, are able to simulate flow with density variation. In this paper, we show that it is possible, by modifying the original equilibrium distribution functions, to capture the discontinuity present in the analytical momentum profile of the two-layered Couette flow with variable density ratios. Investigations carried out earlier were not able to simulate such a flow correctly. Now, with the proposed approach, the new scheme is compatible with the analytical solution, and it is also possible to simulate the two-layered Couette flow with simultaneous density ratios of O(1000) and viscosity ratios of O(100). To test the model in a more complex flow situation, i.e. with non-zero surface tension and a curved interface, an unsteady simulation of an oscillating bubble with variable density ratio is undertaken. The numerical frequency of the bubble is compared to that of the analytical frequency. It is demonstrated that the proposed modification greatly increases the accuracy of the model compared to the original model, i.e. the error can be up to one order of magnitude lower with the proposed enhanced equilibrium distribution functions. The authors believe that this improvement can be made to other RK models as well, which will allow the range of validity of these models to be extended. This is, in fact, what the authors found for the method proposed in this article.     

Extension of a combined analytical/numerical initial value problem solver for unsteady periodic flow

Here we describe analytical and numerical modifications that extend the Differential Reduced Ejector/ mixer Analysis (DREA), a combined analytical/numerical, multiple species ejector/mixing code developed for preliminary design applications, to apply to periodic unsteady flow. An unsteady periodic flow modelling capability opens a range of pertinent simulation problems including pulse detonation engines (PDE), internal combustion engine ICE applications, mixing enhancement and more fundamental fluid dynamic unsteadiness, e.g. fan instability/vortex shedding problems. Although mapping between steady and periodic forms for a scalar equation is a classical problem in applied mathematics, we will show that extension to systems of equations and, moreover, problems with complex initial conditions are more challenging. Additionally, the inherent large gradient initial condition singularities that are characteristic of mixing flows and that have greatly influenced the DREA code formulation, place considerable limitations on the use of numerical solution methods. Fortunately, using the combined analytical–numerical form of the DREA formulation, a successful formulation is developed and described. Comparison of this method with experimental measurements for jet flows with excitation shows reasonable agreement with the simulation. Other flow fields are presented to demonstrate the capabilities of the model. As such, we demonstrate that unsteady periodic effects can be included within the simple, efficient, coarse grid DREA implementation that has been the original intent of the DREA development effort, namely, to provide a viable tool where more complex and expensive models are inappropriate. Copyright © 2002 John Wiley & Sons, Ltd.     

Axisymmetric membrane in adhesive contact with rigid substrates: Analytical solutions under large deformation

The large deformation of an elastic axisymmetric membrane in adhesive contact with a rigid flat punch is studied. Detachment of membrane is analyzed using a critical energy release rate criterion. Two types of incompressible hyperelastic material models are considered: neo-Hookean and a class of materials whose elastic energy density functions are independent of the trace of the Cauchy–Green tensor (I₂-based material). We also include pre-stretch in our formulation and study the stability of detachment process. Closed form analytical solutions for the membrane stresses, deformed profiles and energy release rate are obtained in the regime of large longitudinal stretch. For the I₂-based material, we discover an interesting “pinching” instability where the contact angle suddenly increases in a displacement controlled test. The region of validity of our analytical solutions is determined by comparing them with numerical solutions of the governing equations. We found that the accuracy of our solution improves with pre-stretch; for pre-stretch ratios greater than 1.3, our analytical solution also works well in the small deformation regime.     

The theoretical study of the efficiency of diffusion deposition of nanoaerosols in the extended range of the Peclet numbers

The efficiencies of the diffusion deposition of nanoaerosols for a single fiber for the models of aerosol filter and wire mesh screen are studied numerically in the extended range of the Peclet number Pe. The rectangular periodic cell model for fluid flow and convective-diffusive transport of small aerosol particles is used. Most of the previous theoretical and experimental studies of single fiber diffusion deposition efficiency were for the case of Pe > 1. The array with uniform square or chess grid of fibers and of a row of circular cylindrical fibers are considered as the filter and wire mesh screen models. The flow and particles transport equations are solved numerically using the Boundary Element Method.The obtained numerical data are used to derive the approximate formulas for the deposition efficiency in the entire range of the Peclet number for the various porosities of the filter medium or distances between fibers in a wire mesh screen. The derived dependencies take into account nonlinearity of the deposition efficiency at the low Peclet numbers. The obtained analytical dependencies compare well with the numerical and experimental data.     

Large-Eddy Simulation of the Flow Over a Circular Cylinder at Reynolds Number 3900 Using the OpenFOAM Toolbox

The flow over a circular cylinder at Reynolds number 3900 and Mach number 0.2 was predicted numerically using the technique of large-eddy simulation. The computations were carried out with an O-type curvilinear grid of size of 300 × 300 × 64. The numerical simulations were performed using a second-order finite-volume method with central-difference schemes for the approximation of convective terms. A conventional Smagorinsky and a dynamic k-equation eddy viscosity sub-grid scale models were applied. The integration time interval for data sampling was extended up to 150 vortex shedding periods for the purpose of obtaining a fully converged mean flow field. The present numerical results were found to be in good agreement with existing experimental data and previously obtained large-eddy simulation results. This gives an indication on the adequacy and accuracy of the selected large-eddy simulation technique implemented in the OpenFOAM toolbox.     

Modal Analysis of Monolithic and Jointed Type Cantilever Beams with Non-Uniform Section

Modal analysis of non-uniform bolted structures are of significance in modeling many complex mechanical structures. There are vast literatures available related with the analytical as well as numerical modeling of bolted joint. However, most of the analytical model discuss about the modeling of first mode of uniform structures with single bolted joint. In this paper, we present the modeling of single as well as bolted non-uniform beams using approximate mode shapes. To develop the model, we first carry out experiments to measures the modal frequencies and shapes of the test structures. Subsequently, we also do numerical modeling of non-uniform beams in ANSYS to verify the validity of the Euler-Bernoulli beam theory in developing the analytical models. Finally, using the Euler-Bernoulli beam theory, we obtain the analytical values of frequencies using the approximate the mode shapes. The analytical results are found to be closer to the experimental results with a maximum percentage error of about 15 %. The model presented in the paper can be extended to the mechanical structures with many non-uniform sections with or without bolted joints.     

Flow past a suddenly cooled horizontal plate

The problem of unsteady laminar, incompressible free convection above a horizontal semi-infinite flat plate is studied theoretically. It is assumed that for timet<0 the plate is hotter than its surroundings and at timet=0 the plate is suddenly cooled to the same temperature of its surroundings. Three solutions of the momentum and energy equations are obtained, namely

1. an analytical solution which is valid for small time,
2. an asymptotic analytical solution which is valid for large time, and
3. a numerical solution which matches these two limiting analytical solutions.

It is found that the numerical solution matches the small and large time solutions accurately. Finally, the variation of the velocity, temperature, skin friction and heat transfer on the plate with time are discussed.     

The L-pseudo-solution using stochastic algorithm of Landweber

In this paper, we consider a linear equation Ax=u. A is an operator with an unbounded inverse in a Hilbert space. The right side u does not belong to the range of A. Obviously, a solution in classical sense does not exist and A ^?1 u does not have a sense. To solve this problem arising from many experimental fields of science, where the second member u stems from measurements, we propose a recurrent procedure which converges almost completely and in quadratic mean to L-pseudo-solution and for which we build up a confidence interval. To check the validity of our results, a numerical example which is standard in rheology is proposed.     

Modeling particle sedimentation in viscous fluids with a coupled immersed boundary method and discrete element method

Numerical techniques have increasingly been used to model fluid–particle two-phase flows. Coupling the immersed boundary method (IBM) and discrete element method (DEM) is one promising approach for modeling particulate flows. In this study, IBM was coupled with DEM to improve the reliability and accuracy of IBM for determining the positions of particles during the sedimentation process within viscous fluids. The required ratio of the particle diameter to the grid size (D/dx) was determined by comparing the simulation results with the analytical solution and experimental data. A dynamic mesh refinement model was utilised in the IBM model to refine the computational fluid dynamics grid near the particles. In addition, an optimum coupling interval between the IBM and DEM models was determined based on the experimental results of a single particle sedimentation within silicon oil at a Reynolds number of 1.5. The experimental results and the analytical solution were then utilised to validate the IBM–DEM model at Reynolds numbers of 4.1, 11.6, and 31.9. Finally, the validated model was utilised to investigate the sedimentation process for more than one particle by modeling the drafting-kissing-tumbling process and the Boycott phenomenon. Benchmark tests showed that the IBM–DEM technique preserves the advantages of DEM for tracking a group of particles, while the IBM provides a reliable and accurate approach for modeling the particle–fluid interaction.     

Series solution of the temperature distribution in the Falkner–Skan wedge flow by the homotopy analysis method

A new analytical method, namely the homotopy analysis method (HAM), has been applied to investigate the temperature field associated with the Falkner–Skan boundary-layer problem, and a series solution is provided in this paper. The results of the present work show agreement with those of numerical solutions in a large range of Prandtl numbers (0 < Pr ≤ 100), which demonstrates the validity of the present analysis.     

Dislocation shielding of a cohesive crack

Dislocation interaction with a cohesive crack is of increasing importance to computational modelling of crack nucleation/growth and related toughening mechanisms in confined structures and under cyclic fatigue conditions. Here, dislocation shielding of a Dugdale cohesive crack described by a rectangular traction-separation law is studied. The shielding is completely characterized by three non-dimensional parameters representing the effective fracture toughness, the cohesive strength, and the distance between the dislocations and the crack tip. A closed form analytical solution shows that, while the classical singular crack model predicts that a dislocation can shield or anti-shield a crack depending on the sign of its Burgers vector, at low cohesive strengths a dislocation always shields the cohesive crack irrespective of the Burgers vector. A numerical study shows the transition in shielding from the classical solution of Lin and Thomson (1986) in the high strength limit to the solution in the low strength limit. An asymptotic analysis yields an approximate analytical model for the shielding over the full range of cohesive strengths. A discrete dislocation (DD) simulation of a large (>10³) number of edge dislocations interacting with a cohesive crack described by a trapezoidal traction-separation law confirms the transition in shielding, showing that the cohesive crack does behave like a singular crack at very high cohesive strengths (∼7 GPa), but that significant deviations in shielding between singular and cohesive crack predictions arise at cohesive strengths around 1GPa, consistent with the analytic models. Both analytical and numerical studies indicate that an appropriate crack tip model is essential for accurately quantifying dislocation shielding for cohesive strengths in the GPa range.     

Elongational effects in the flow of viscoelastic fluid through a wavy channel

The numerical simulation of the viscoelastic flow through a wavy channel was carried out using the modified Giesekus model. It was found that the excess pressure loss relates to the stretch-thickening properties of elongational viscosity and the geometry of the wavy channel through a large elongational component of the flow at the winding part of the channel. The profiles of the axial component of the velocity become significantly asymmetric when the excess pressure loss occurs. Furthermore, the velocity profiles of a 0.1 wt% solution of polyacrylamide were measured using laser Doppler velocimetry. The results of these measurements are compared to the numerical results. Received: 30 June 1998 Accepted: 20 May 1999     

Comparison of approximate analytical and numerical solutions for heat fluxes in viscous supersonic flow past a body

The axisymmetric supersonic laminar flow of a viscous perfect gas past a blunt body is studied over a wide range of Reynolds and Mach numbers. An approximate analytical solution for the distribution of the heat flux, divided by its value at the stagnation point, over the body surface is derived. The expression for the relative heat flux depends only on the body geometry and the pressure distribution over the body surface. To evaluate the accuracy and the range of validity of the formulas obtained, a comparison with numerical solutions of the complete Navier-Stokes equations was carried out for bodies of various shapes, such as spheres, ellipsoids, and paraboloids.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 125–132, January–February, 1996.