A class of subgrid-scale models preserving the symmetry group of Navier–Stokes equations

Navier–Stokes equations (NS) admit transformations which transform a solution to another solution (galilean transformation, scaling transformation, …). They also admit viscosity dependent transformations which transform a solution to a solution of another NS with different viscosity. These particular transformations are called symmetries of NS. Each of them has a physical role (such as conservation laws, …). A consistent turbulence model should then remain invariant under these symmetry transformations. Unfortunately, this is not the case of several models.In this article, a class of subgrid-scale models preserving the symmetries of NS is built. This class is then refined such that the models respect the second law of thermodynamics. One of the simplest models of the class is tested to the flow in a ventilated room. Better results than those provided by Smagorinsky and dynamic models are obtained.     

The average path length of scale free networks

In this paper, the exact solution of average path length in Barabási–Albert model is given. The average path length is an important property of networks and attracts much attention in many areas. The Barabási–Albert model, also called scale free model, is a popular model used in modeling real systems. Hence it is valuable for us to examine the average path length of scale free model. There are two answers, regarding the exact solution for the average path length of scale free networks, already provided by Newman and Bollobas respectively. As Newman proposed, the average path length grows as log(n) with the network size n. However, Bollobas suggested that while it was true when m = 1, the answer changed to log(n)/log(log(n)) when m > 1. In this paper, as we propose, the exact solution of average path length of BA model should approach log(n)/log(log(n)) regardless the value of m. Finally, the simulation is presented to show the validity of our result.     

Magnetohydrodynamic stationary and oscillatory convective stability in a mushy layer during binary alloy solidification

For the case of solidification of a bottom cooled binary alloy, the magnetohydrodynamic stationary and oscillatory convective stability in the mushy layer is investigated analytically using normal mode linear stability analysis. In the limit of large Stefan number (S_t), a near–eutectic approximation with large far field temperature is considered in the present research. To ascertain the instability in the mushy layer, the strength of the superimposed magnetic field is so chosen that it corresponds to a given mush Hartmann number (Ha_m) of the problem. The results are presented for various values of mush Hartmann numbers in the range, 0 ≤ Ha_m ≤ 50. The critical Rayleigh number for stationary convection shows a linear relationship with increasing Ha_m. The magnetohydrodynamic effect imparts a stabilizing influence during stationary convection. In comparison to that of the stationary convective mode, the oscillatory mode appears to be critically susceptible at higher values of β (β = S_t/℘² ϒ², ℘ is the compositional ratio, ϒ = 1 + S_t/℘), and vice versa for lower β values. Analogous to the behavior for stationary convection, the magnetic field also offers a stabilizing effect in oscillatory convection and thus influences global stability of the mushy layer. Increasing magnetic strength shows reduction in the wavenumber and in the number of rolls formed in the mushy layer.     

Environmental flows for the Yangtze Estuary based on salinity objectives

Environmental flows based on the salinity objectives for China’s Yangtze Estuary were quantified to assess the impacts of changing freshwater inflow on the estuarine ecosystem. The salinity objectives for three benthos habitats in the Yangtze Estuary were studied to determine their requirements during the most critical period of the year. The temporal variation in the natural monthly river discharge represented the temporal variation of the ecological objectives used to maintain natural flow regimes. A numerical model, coupled a hydrodynamic processes model with a salinity model and validated using field data of tidal height, current velocity, and salinity at different stations, was developed to simulate the spatial distribution of salinity as a function of the variation in freshwater inflows. The prediction results of the proposed model agreed well with the field data. Considering the salinity objectives for the different habitats, the environmental flow requirements for the Yangtze Estuary were determined. Annual environmental flows should be 9.63 × 10¹¹, 6.32 × 10¹¹, and 4.70 × 10¹¹ m³ for the high, medium, and minimum objectives, which are equivalent to 104%, 68%, and 50% of the annual river discharge, respectively. The periods from September to November should be considered as the critical seasons to maintain the minimum environmental flows in the Yangtze Estuary.     

The new airfoil model NACA0015, modal analysis and flutter properties

The experimental airfoil model NACA0015 was used to study aeroelastic phenomena during self-excited profile vibration. It provides data for control of aeroelastic calculation programs at subsonic speeds of the stream. The model movability is two-dimensional with two-degree of freedom dynamic system, one in pitch and the second in plunge, and is proposed to be a dynamic system having two near corresponding eigenfrequencies. To quantitatively evaluate flow field using interferometry, a special test section design and profile was constructed. It utilized a large visual field for the optical system together with the option of changing support stiffness for both degrees of freedom. In this paper experimental results from the range of Reynolds numbers Re = (2.63–2.83) 10⁵ are published. The identified eigenvalues and eigenmodes for zero flow velocity are compared with measured flutter properties (frequency, modes and time evolutions) of the airfoil.     

Un modelo numérico para la simulación de disolución de precipitados en aleaciones de aluminio con endurecimiento utilizando redes neuronales

The motivation of this work is the modeling of the hardening precipitate and hardness evolutions of fully hardened heat treatable aluminium alloys during friction stir welding (FSW) and/or heat treatment processes. The model used is based on the kinetics of dissolution of precipitates model for hardened aluminium alloys given by Myhr and Grong (1991). This model contains a single independent variable, the time, and a single state variable, the volume fraction of hardening precipitates. A key point of this model is the identification of the effective activation energy for precipitates dissolution and the master curve defining the model, which was given by a look-up table. The goal of this work is to find an estimation of the effective activation energy and to model the dissolution rate of hardening precipitate in aluminium alloys using neural networks, avoiding the use of look-up tables. For this purpose a new and more convenient parametrization of the master curve is defined, a neural networks class is proposed, an objective functional is defined and a variational problem including independent parameters is solved. The novel methodology has been applied to different aluminium alloys, including the AA 6005A T6, AA 7449 T79 and AA 2198 T8. Experimental tests have been carried out for those aluminium alloys in order to get the HV1 hardness after isothermal heat treatments at different temperatures and for different treatment time durations. The effective activation energy for hardening precipitates dissolution and the master curve of the model have been obtained, using different network architectures, for the aluminium alloys considered in this work.     

Complex bifurcations in the oscillatory reaction model

The mixing of different types of bifurcations, i.e. supercritical Andronov–Hopf (SAH), double loop (DL) and saddle-loop (SL) bifurcations in the vicinity of their total annihilation, is examined on the highly nonlinear six-variable model for the Bray–Liebhafsky (BL) oscillatory reaction under continuously well-stirred tank reactor (CSTR) conditions. For this kind of the reaction system where the law of mass conservation is additional constraint that must be satisfied and where because of that, some simple bifurcations cannot be formed independently to the others, the considered transformations of the bifurcations are particularly important. That is why as the control parameters for bifurcation analysis, the specific flow rate (j₀), as well as the inflow hydrogen peroxide concentration (h = [H₂O₂]_in), were used. The complex bifurcations obtained from numerical simulations are compared with some experimental results. It was shown that these complex bifurcations cannot be easily recognized in experimental investigations without knowing their evolution.     

An approximate solution of nonlinear fractional reaction–diffusion equation

The article presents a mathematical model of nonlinear reaction diffusion equation with fractional time derivative α (0 < α ? 1) in the form of a rapidly convergent series with easily computable components. Fractional reaction diffusion equation is used for modeling of merging travel solutions in nonlinear system for popular dynamics. The fractional derivatives are described in the Caputo sense. The anomalous behaviors of the nonlinear problems in the form of sub- and super-diffusion due to the presence of reaction term are shown graphically for different particular cases.     

Effect of variable density on hydromagnetic mixed convection flow of a non-Newtonian fluid past a moving vertical plate

The effect of temperature-dependent density on MHD mixed convection flow of power-law fluid past a moving semi-infinite vertical plate for high temperature differences between the plate and the ambient fluid is studied. The fluid density is assumed to decrease exponentially with temperature. The usual Boussinesq approximations are not considered due to the large temperature differences. The surface temperature of the moving plate was assumed to vary according to a power-law form, that is, T_w(x) = T_∞ + Ax^γ. The fluid is permeated by a uniform magnetic field imposed perpendicularly to the plate on the assumption of small magnetic Reynolds number. A numerical shooting algorithm for two unknown initial conditions with fourth-order Runge–Kutta integration scheme has been used to solve the coupled non-linear boundary value problem. The effects of various parameters on the velocity and temperature profiles as well as the local skin-friction coefficient and the local Nusselt number are presented graphically and in the tabular form. The results show that application of Boussinesq approximations in a non-Newtonian fluid subjected to high temperature differences gives a significant error in the values of the skin-friction coefficient and the application of an external magnetic field reduces this error markedly in the case of shear-thickening fluid.     

Evaluation of thermal comfort conditions in a classroom equipped with radiant cooling systems and subjected to uniform convective environment

The aim of this work is to evaluate numerically the human thermal response that 24 students and 1 teacher feel in a classroom equipped with radiant cooling systems and subjected to uniform convective environments, in lightly warm conditions. The evolution of thermal comfort conditions, using the PMV index, is made by the multi-nodal human thermal comfort model.In this numerical model, that works in transient or steady-state conditions and simulates simultaneously a group of persons, the three-dimensional body is divided in 24 cylindrical and 1 spherical elements. Each element is divided in four parts (core, muscle, fat and skin), sub-divided in several layers, and protected by several clothing layers. This numerical model is divided in six parts: human body thermal system, clothing thermal system, integral equations resolution system, thermoregulatory system, heat exchange between the body and the environment and thermal comfort evaluation.Seven different radiant systems are combined to three convective environments. In the radiant systems (1) no radiant system without warmed curtain, (2) no radiant system with warmed curtain, (3) radiant floors cooling system with warmed curtain, (4) radiant panels cooling system with warmed curtain, (5) radiant ceiling cooling system with warmed curtain, (6) radiant floor and panels cooling system with warmed curtain and (7) radiant ceiling and panels cooling system with warmed curtain are analysed, while in the convective environments (1) without air velocity field and with uniform air velocity field of (2) 0.2 m/s and (3) 0.6 m/s are also analysed. The internal air temperature and internal surfaces temperature are 28 °C, the radiant cooling surfaces temperature are 19 °C and the warmed internal curtains surfaces temperatures, subjected to direct solar radiation, are 40 °C.The numerical model calculates the Mean Radiant Temperature field, the human bodies’ temperatures field and the thermal comfort level, for the 25 occupants, for the 21 analysed situations.Without uniform air velocity field, when only one individual radiant cooling system is used, the Predicted Percentage of Dissatisfied people is lowest when the radiant floor cooling system is applied and is highest when the radiant panel cooling system is applied. When are combined the radiant ceiling or the floor cooling systems with the radiant panel cooling system the Predicted Percentage of Dissatisfied people decreases.When the uniform air velocity increases the thermal comfort level, that the occupants are subjected, increases. When the radiant floor cooling system or the combination of radiant floor and panel cooling systems without uniform air velocity field is applied, the Category C is verified for some occupants. However, with a convective uniform air velocity field of 0.2 m/s the Category B is verified and with a convective uniform air velocity field of 0.6 m/s the Category A is verify for some occupants. In the last situation the Category C is verified, in general, for all occupants.     

The applications of a higher-dimensional Lie algebra and its decomposed subalgebras

With the help of invertible linear transformations and the known Lie algebras, a higher-dimensional 6 × 6 matrix Lie algebra sμ(6) is constructed. It follows a type of new loop algebra is presented. By using a (2 + 1)-dimensional partial-differential equation hierarchy we obtain the integrable coupling of the (2 + 1)-dimensional KN integrable hierarchy, then its corresponding Hamiltonian structure is worked out by employing the quadratic-form identity. Furthermore, a higher-dimensional Lie algebra denoted by E, is given by decomposing the Lie algebra sμ(6), then a discrete lattice integrable coupling system is produced. A remarkable feature of the Lie algebras sμ(6) and E is used to directly construct integrable couplings.     

The conservation laws and group properties of the equations of gas dynamics with zero velocity of sound

All the conservation laws of zero order are obtained by the method of A-operators for a system of n-dimensional (n ≥ 1) equations of gas dynamics with zero velocity of sound. A group subdivision is carried out of this system with respect to an infinite subgroup, which is a normal divider of its main Lie group of transformations; the main group of the resolving system is obtained. First-order non-local symmetries are obtained for the initial system. A special choice of the mass Lagrange variables enables this system to be converted to a reduced system equivalent to it, containing n - 1 spatial variables, which, for n = 2, is written in the form of a one-dimensional complex heat-conduction equation using complex dependent and independent variables.     

A geometrically nonlinear finite element model of nanomaterials with consideration of surface effects

In conventional continuum mechanics, the surface energy is usually small and negligible. But at nano-length scale, it becomes a significant part of the total elastic energy due to the high specific surface area of nanomaterials. A geometrically nonlinear finite element (FE) model of nanomaterials with considering surface effects is developed in this paper. The aim is to extend the conventional finite element method (FEM) to analyze the size-dependent mechanical properties of nanomaterials. A numerical example, analysis of InAs quantum dot (QD) on GaAs (0 0 1) substrate, is given in this paper to verify the validity of the method and demonstrate surface effects on the stress fields of QDs.     

Links between some (2 + 1)-dimensional nonlinear evolution equations and Painlevé-II equations

Using the idea of transformation, some links between (2 + 1)-dimensional nonlinear evolution equations and the ordinary differential equations Painlevé-II equations has been illustrated. The Kadomtsev–Petviashvili (KP) equation, generalized (2 + 1)-dimensional break soliton equation and (2 + 1)-dimensional Boussinesq equation are researched. As a result, some new interesting results about these (2 + 1)-dimensional PDEs have been obtained, such as the exact solutions with arbitrary functions, rich rational solutions and the nontrivial Bäcklund transformations have been derived.     

Heat transfer analysis of unsteady boundary layer flow by homotopy analysis method

This paper aims to present complete analytic solution to the unsteady heat transfer flow of an incompressible viscous fluid over a permeable plane wall. The flow is started due to an impulsively stretching porous plate. Homotopy analysis method (HAM) has been used to get accurate and complete analytic solution. The solution is uniformly valid for all time τ ∈ [0, ∞) throughout the spatial domain η ∈ [0, ∞). The accuracy of the present results is shown by giving a comparison between the present results and the results already present in the literature. This comparison proves the validity and accuracy of our present results. Finally, the effects of different parameters on temperature distribution are discussed through graphs.     

Simulación numérica del flujo turbulento de aire con gotas dispersas de agua a través de separadores de torres de refrigeración

This work presents a numerical study on the turbulent flow of air with dispersed water droplets in separators of mechanical cooling towers. The averaged Navier-Stokes equations are discretised through a finite volume method, using the Fluent and Phoenics codes, and alternatively employing the turbulence models k ? ?, k ? ω and the Reynolds stress model, with low-Re version and wall enhanced treatment refinements. The results obtained are compared with numerical and experimental results taken from the literature. The degree of accuracy obtained with each of the considered models of turbulence is stated. The influence of considering whether or not the simulation of the turbulent dispersion of droplets is analyzed, as well as the effects of other relevant parameters on the collection efficiency and the coefficient of pressure drop. Focusing on four specific eliminators (‘Belgian wave’, ‘H1-V’, ‘L-shaped’ and ‘Zig-zag’), the following ranges of parameters are outlined: 1 ≤ U_e ≤ 5 m/s for the entrance velocity, 2 ≤ D_p ≤ 50 μm for the droplet diameter, 650 ≤ Re ≤ 8.500 for Reynolds number, and 0.05 ≤ P_i ≤ 5 for the inertial parameter. Results reached alternately with Fluent and Phoenics codes are compared. The best results correspond to the simulations performed with Fluent, using the SST k ? ω turbulence model, with values of the dimensionless scaled distance to wall y⁺ in the range 0.2 to 0.5. Finally, correlations are presented to predict the conditions for maximum collection efficiency (100 %), depending on the geometric parameter of removal efficiency of each of the separators, which is introduced in this work.     

Invertible linear transformations and the Lie algebras

With the help of invertible linear transformations and the known Lie algebras, a way to generate new Lie algebras is given. These Lie algebras obtained have a common feature, i.e. integrable couplings of solitary hierarchies could be obtained by using them, specially, the Hamiltonian structures of them could be worked out. Some ways to construct the loop algebras of the Lie algebras are presented. It follows that some various loop algebras are given. In addition, a few new Lie algebras are explicitly constructed in terms of the classification of Lie algebras proposed by Ma Wen-Xiu, which are bases for obtaining new Lie algebras by using invertible linear transformations. Finally, some solutions of a (2 + 1)-dimensional partial-differential equation hierarchy are obtained, whose Hamiltonian form-expressions are manifested by using the quadratic-form identity.     

Construction of analytic solution to chaotic dynamical systems using the Homotopy analysis method

Based on a new kind of analytic method, namely the Homotopy analysis method, an analytic approach to solve non-linear, chaotic system of ordinary differential equations is presented. The method is applied to Lorenz system; this system depends on the three parameters: σ, b and the so-called bifurcation parameter R are real constants. Two cases are considered. The first case is when R = 20.5 which corresponds to the transition region and the second case corresponds to R = 23.5 which corresponds to the chaotic region.The validity of the method is verified by comparing the approximation series solution with the results obtained using the standard numerical techniques such as Runge-Kutta method.     

Modelling approach for biological control of insect pest by releasing infected pest

Models of biological control have a long history of theoretical development that have focused on the interactions between a predator and a prey. Here we have extended the classical epidemic model to include a continuous and impulsive pest control strategies by releasing the infected pests bred in laboratory. For the continuous model, the results imply that the susceptible pest goes to extinct if the threshold condition R₀ < 1. While R₀ > 1, the positive equilibrium of continuous model is globally asymptotically stable. Similarly, the threshold condition which guarantees the global stability of the susceptible pest-eradication periodic solution is obtained for the model with impulsive control strategy. Consequently, based on the results obtained in this paper, the control strategies which maintain the pests below an acceptably low level are discussed by controlling the release rate and impulsive period. Finally, the biological implications of the results and the efficiency of two control strategies are also discussed.     

On mathematical models and numerical simulation of the fluidization of polydisperse suspensions

The well-known Masliyah–Lockett–Bassoon (MLB) model for sedimentation of small particles is extended to fluidization of polydisperse suspensions. For N particle species that differ in size and density, this model leads to a first-order system of N conservation laws, which in general is of mixed (in the case N = 2, hyperbolic–elliptic) type. By a simple algebraic steady-state analysis, we derive necessary compatibility conditions on the size and density parameters that admit the formation of stationary fluidized beds. We then proceed to determine the composition of polydisperse fluidized beds of given compatible species by varying the fluidization velocity and the initial composition of the suspensions, and prove that, within the framework of the MLB model combined with the Richardson–Zaki formula, the constructed bidisperse beds always cause the equations to be hyperbolic. This means that these states are always predicted to be stable. The transient behaviour of the MLB model applied to fluidization is illustrated by three numerical examples, in which the system of conservation laws is solved for N = 2, N = 3 and N = 5, respectively. These examples illustrate the effects of bed expansion and layer inversion caused by successively increasing the applied fluidization velocity and show that the predicted fluidized states are indeed attained.