Search results

The purpose of this paper is to present a numerical and an analytical study of the thermohaline convection with Soret effect in a square enclosure filled with a binary fluid mixture.

The horizontal boundaries of the enclosure are impermeable and heated from below while its vertical walls are assumed to be adiabatic and impermeable. The Navier‐Stokes equations under the Boussinesq‐Oberbeck approximation are solved numerically. The results are given for different values of the separation ratio. The critical Rayleigh number at the onset of convection is determined analytically and numerically. The Hopf frequency at the onset of convection is obtained.

The existence of two stable stationary bifurcation branches is illustrated. Furthermore, it is shown that the existence of stable traveling waves in the transition from one branch to the other depends on the value of the separation ratio. For some values of Rayleigh number, asymmetric flows are observed. A good agreement is found between the numerical solution and analytical analysis.

The present work is the first to consider thermosolutal convection with Soret effect in a square enclosure.

The purpose of this paper is to deal with the study of plane waves and fundamental solution in a modified couple stress generalized thermoelastic solid with three-phase-lag (TPL) model of thermoelasticity.

It is found that for two-dimensional model, there exists two longitudinal waves, namely, longitudinal wave (P-wave), thermal wave (T-wave), and a set of coupled transverse waves (SV1 and SV2 waves). In addition, the fundamental solution for the system of differential equations for steady oscillations in terms of elementary functions has been constructed. Some properties of fundamental solution are also established. Various particular cases of interest are also deduced from the present investigations and compared with the known results.

The phase velocity, attenuation coefficient, specific loss and penetration depth are computed numerically and presented graphically to see the effect of TPL model, dual-phase-lag (DPL) model and GN-III model in the presence of couple stress parameter.

The results are compared with couple stress TPL model, couple stress DPL model and GN-III model.

In the recent work, a new blast-wave impact-mitigation concept involving the use of a protective structure consisting of bimolecular reactants (polyvinyl pyridine+cyclohexyl chloride), capable of undergoing a chemical reaction (to form polyvinyl pyridinium ionic salt) under shockwave loading conditions, was investigated using all-atom reactive equilibrium and non-equilibrium molecular-dynamics analyses. The purpose of this paper is to reveal the beneficial shockwave dispersion/attenuation effects offered by the chemical reaction, direct simulations of a fully supported single planar shockwave propagating through the reactive mixture were carried out, and the structure of the shock front examined as a function of the extent of the chemical reaction (i.e. as a function of the strength of the incident shockwave). The results obtained clearly revealed that chemical reactions give rise to considerable broadening of the shockwave front. In the present work, the effect of chemical reactions and the structure of the shockwaves are investigated at the continuum level.

Specifically, the problem of the (conserved) linear-momentum accompanying the interaction of an incident shockwave with the protective-structure/protected-structure material interface has been investigated, within the steady-wave/structured-shock computational framework, in order to demonstrate and quantify an increase in the time period over which the momentum is transferred and a reduction in the peak loading experienced by the protected structure, both brought about by the occurrence of the chemical reaction (within the protective structure).

The results obtained clearly revealed the beneficial shock-mitigation effects offered by a protective structure capable of undergoing a chemical reaction under shock-loading conditions.

To the authors’ knowledge, the present manuscript is the first report dealing with a continuum-level analysis of the blast-mitigation potential of chemical reactions.

The purpose of this paper is to determine both analytically and numerically the existence of smooth, cusped and sharp shock wave solutions to a one-dimensional model of microfluidic droplet ensembles, water flow in unsaturated flows, infiltration, etc., as functions of the powers of the convection and diffusion fluxes and upstream boundary condition; to study numerically the evolution of the wave for two different initial conditions; and to assess the accuracy of several finite difference methods for the solution of the degenerate, nonlinear, advection--diffusion equation that governs the model.

The theory of ordinary differential equations and several explicit, finite difference methods that use first- and second-order, accurate upwind, central and compact discretizations for the convection terms are used to determine the analytical solution for steadily propagating waves and the evolution of the wave fronts from hyperbolic tangent and piecewise linear initial conditions to steadily propagating waves, respectively. The amplitude and phase errors of the semi-discrete schemes are determined analytically and the accuracy of the discrete methods is assessed.

For non-zero upstream boundary conditions, it has been found both analytically and numerically that the shock wave is smooth and its steepness increases as the power of the diffusion term is increased and as the upstream boundary value is decreased. For zero upstream boundary conditions, smooth, cusped and sharp shock waves may be encountered depending on the powers of the convection and diffusion terms. For a linear diffusion flux, the shock wave is smooth, whereas, for a quadratic diffusion flux, the wave exhibits a cusped front whose left spatial derivative decreases as the power of the convection term is increased. For higher nonlinear diffusion fluxes, a sharp shock wave is observed. The wave speed decreases as the powers of both the convection and the diffusion terms are increased. The evolution of the solution from hyperbolic tangent and piecewise linear initial conditions shows that the wave back adapts rapidly to its final steady value, whereas the wave front takes much longer, especially for piecewise linear initial conditions, but the steady wave profile and speed are independent of the initial conditions. It is also shown that discretization of the nonlinear diffusion flux plays a more important role in the accuracy of first- and second-order upwind discretizations of the convection term than either a conservative or a non-conservative discretization of the latter. Second-order upwind and compact discretizations of the convection terms are shown to exhibit oscillations at the foot of the wave’s front where the solution is nil but its left spatial derivative is largest. The results obtained with a conservative, centered second--order accurate finite difference method are found to be in good agreement with those of the second-order accurate, central-upwind Kurganov--Tadmor method which is a non-oscillatory high-resolution shock-capturing procedure, but differ greatly from those obtained with a non-conservative, centered, second-order accurate scheme, where the gradients are largest.

A new, one-dimensional model for microfluidic droplet transport, water flow in unsaturated flows, infiltration, etc., that includes high-order convection fluxes and degenerate diffusion, is proposed and studied both analytically and numerically. Its smooth, cusped and sharp shock wave solutions have been determined analytically as functions of the powers of the nonlinear convection and diffusion fluxes and the boundary conditions. These solutions are used to assess the accuracy of several finite difference methods that use different orders of accuracy in space, and different discretizations of the convection and diffusion fluxes, and can be used to assess the accuracy of other numerical procedures for one-dimensional, degenerate, convection--diffusion equations.

The choke flutter is a fluid-structure interaction that can lead to the failure of fan or compressor blade in turbojet engines. In ultra high bypass ratio (UHBR) fans, the choke flutter appears at part-speed regimes and at low or negative incidence when a strong shock-wave chokes the blade to blade channel. The purpose of this study is to locate the main excitation sources and improving the understanding of the different work exchange mechanisms. This work contributes to avoiding deficient and dangerous fan design.

In this paper, an UHBR fan is analyzed using a time-linearized Reynolds-averaged Navier–Stokes equation solver to investigate the choke flutter. The steady-state and the imposed vibration (inter blade phase angle, reduced frequency and mode shape) are selected to be in choke flutter situation. Superposition principle induced by the linearization allow to decompose the blade in numerous small subsections to track the contribution of each local vibration to the global damping. All simulations have been performed on a two-dimensional blade to blade extraction.

Result analysis points to a restricted number of excitation sources at the trailing edge which induce a large part of the work exchange in a limited region of the airfoil. Main phenomena suspected are the shock-wave motion and the shock-wave/boundary layer interaction.

An original excitation source tracking methodology allowed by the linearized calculation is addressed and applied to a UHBR fan test case.

The purpose of this paper is to study of propagation of plane wave and the fundamental solution of the system of differential equations in the theory of a microstretch thermoelastic diffusion medium in phase-lag models for the case of steady oscillations in terms of elementary functions.

Wave propagation technique along with the numerical methods for computation using MATLAB software has been applied to investigate the problem.

Characteristics of waves like phase velocity and attenuation coefficient are computed numerically and depicted graphically. It is found that due to the presence of diffusion effect, these characteristics get influenced significantly. However, due to decoupling of CD-I and CD-II waves from rest of other, no effect on these characteristics can be perceived.

Basic properties of the fundamental solution are established by introducing the dual-phase-lag diffusion (DPLD) and dual-phase-lag heat transfer (DPLT) models.

The purpose of this paper is to introduce and analyze a new blast-wave impact-mitigation concept using advanced computational methods and tools. The concept involves the use of a protective structure consisting of bimolecular reactants displaying a number of critical characteristics, including: a high level of thermodynamic stability under ambient conditions (to ensure a long shelf-life of the protective structure); the capability to undergo fast/large-yield chemical reactions under blast-impact induced shock-loading conditions; large negative activation and reaction volumes to provide effective attenuation of the pressure-dominated shockwave stress field through the volumetric-energy storing effects; and a large activation energy for efficient energy dissipation. The case of a particular bimolecular chemical reaction involving polyvinyl pyridine and cyclohexyl chloride as reactants and polyvinyl pyridinium ionic salt as the reaction product is analyzed.

Direct simulations of single planar shockwave propagations through the reactive mixture are carried out, and the structure of the shock front examined, as a function of the occurrence of the chemical reaction. To properly capture the shockwave-induced initiation of the chemical reactions during an impact event, all the calculations carried out in the present work involved the use of all-atom molecular-level equilibrium and non-equilibrium reactive molecular-dynamics simulations. In other words, atomic bonding is not pre-assigned, but is rather determined dynamically and adaptively using the concepts of the bond order and atomic valence.

The results obtained clearly reveal that when the chemical reactions are allowed to take place at the shock front and in the shockwave, the resulting shock front undergoes a considerable level of dispersion. Consequently, the (conserved) linear momentum is transferred (during the interaction of the protective-structure borne shockwaves with the protected structure) to the protected structure over a longer time period, while the peak loading experienced by the protected structure is substantially reduced.

To the authors’ knowledge, the present work is the first attempt to simulate shock-induced chemical reactions at the molecular level, for purposes of blast-mitigation.

The purpose of this paper is to study the wave propagation in a generalized piezothermoelastic rotating bar of circular cross-section using three-dimensional linear theory of elasticity.

A mathematical model is developed to study the wave propagation in a generalized piezothermelastic rotating bar of circular cross-section by using Lord-Shulman (LS) and Green-Lindsay (GL) theory of thermoelasticity. After developing the formal solution of the mathematical model consisting of partial differential equations, the frequency equations have been derived by using the thermally insulated/isothermal and electrically shorted/charge free boundary conditions prevailing at the surface of the circular cross-sectional bar. The roots of the frequency equation are obtained by using the secant method, applicable for complex roots.

In order to include the time requirement for the acceleration of the heat flow and the coupling between the temperature and strain fields, the analytical terms have been derived for the non-classical thermo-elastic theories, LS and GL theory. The computed physical quantities such as thermo-mechanical coupling, electro-mechanical coupling, frequency shift, specific loss and frequency have been presented in the form of dispersion curves. From the graphical patterns of the structure, the effect of thermal relaxation times and the rotational speed as well as the anisotropy of the of the material on the various considered wave characteristics is more significant and dominant in the flexural modes of vibration. The effect of such physical quantities provides the foundation for the construction of temperature sensors, acoustic sensor and rotating gyroscope.

In this paper, the influence of thermal relaxation times and rotational speed on the wave number with thermo-mechanical coupling, electro-mechanical coupling, frequency shift, specific loss and frequency has been observed and are presented as dispersion curves. The effect of thermal relaxation time and rotational speed on wave number for the case of generalized piezothermoelastic material of circular cross-section was never reported in the literature. These results are new and original.

The purpose of this paper is to mainly review the state-of-the-art developments in the field of hydrodynamics of offshore pipelines, identifying the key tools for analysis of pipeline free spans, their applications, their qualifying characteristics and capabilities and limitations.

These different analytical, numerical and semi-empirical tools available for predicting such hydrodynamic loads and their effects include VIVANA, PIPESIN, VIVSIM, SIMULATOR, FATFREE, amongst others. Inherent in these models are current effects, wave effects and/ or pipe–soil interactions.

Amongst these models, the most attention was given to the new VIVANA model because this model take into account the vortex-induced effects with respect to free-spanning pipelines (which have dominant effect in the span analysis in deep water) better than other semi-empirical models (such as Shear 7). Recent improvements in VIVANA include its ability to have arbitrary variation in speed and direction of current, as well as the ability for calculation of pure IL and combined IL-CF response. Improvements in fatigue assessments at free spans, i.e. pipe–soil interaction have been achieved through the combined frequency domain and non-linear time domain analysis methodology adopted. Semi-empirical models are still the de facto currently used in the design of free-spanning pipelines. However, there is need for further research on free-span hydrodynamic coefficients and on how in-line and cross-flow vibrations interact. Again, there is still the challenge due to VIV complexity in fully understanding the fluid structure interaction problem, as there is no consolidated procedure for its analysis. It has been observed that there is large scatter between the different codes adopted in the prediction of fatigue damage, as there lacks full-scale test data devoted to determination/validation of the coefficients used in the semi-empirical models. A case study of the preliminary design of a typical 48 in. pipeline has been presented in this study to demonstrate the use of the free-span analysis tool, DNV RP F105. Excel spreadsheet has been applied in the execution of formulas.

This review paper is the first of its kind to study the state-of-the-art development in pipeline free-span analysis models and demonstrate the use of analysis tool, DNV for MAFSL calculation. Hence, information obtained from this paper would be invaluable in assisting designers both in the industry and academia.

A novel numerical formulation of the two‐phase macroscopic balance equations governing the flow field in incompressible porous media is presented. The numerical model makes use of the weighted average flux method and total variation diminishing flux limiting techniques, and results in a second‐order accurate scheme. A shock tube study was carried out to examine the interaction of a normal shock wave with a thin layer of porous, incompressible cellular ceramic foam. Particular attention was paid to the transmitted and reflected flow fields. The numerical model was used to simulate the experimental test cases, and their results compared with a view to validate the numerical model. A phenomenological model is proposed to explain the behaviour of the transmitted flow field.