Strain hardening in 2D discrete dislocation dynamics simulations: A new ‘2.5D’ algorithm

The two-dimensional discrete dislocation dynamics (2D DD) method, consisting of parallel straight edge dislocations gliding on independent slip systems in a plane strain model of a crystal, is often used to study complicated boundary value problems in crystal plasticity. However, the absence of truly three dimensional mechanisms such as junction formation means that forest hardening cannot be modeled, unless additional so-called ‘2.5D’ constitutive rules are prescribed for short-range dislocation interactions. Here, results from three dimensional dislocation dynamics (3D DD) simulations in an FCC material are used to define new constitutive rules for short-range interactions and junction formation between dislocations on intersecting slip systems in 2D. The mutual strengthening effect of junctions on preexisting obstacles, such as precipitates or grain boundaries, is also accounted for in the model. The new ‘2.5D’ DD model, with no arbitrary adjustable parameters beyond those obtained from lower scale simulation methods, is shown to predict athermal hardening rates, differences in flow behavior for single and multiple slip, and latent hardening ratios. All these phenomena are well-established in the plasticity of crystals and quantitative results predicted by the model are in good agreement with experimental observations.     

Elastic cables–rigid body coupled dynamics: asymptotic modeling and analysis

An elastic cables–rigid body coupled model is proposed for investigating dynamic interactions between cables’ nonlinear transversal vibrations and boundary tower’s torsional dynamics, arising in large transmission line–tower systems and suspended cable–bridge tower systems. By introducing a weak torsion assumption and a large moment of inertia for the tower, an asymptotic expansion of cables–tower coupled dynamics is conducted in a weakly nonlinear framework, and a cables–tower reduced coupled model is eventually established. After model’s validations using direct numerical simulations, two distinct kinds of coupled dynamics are fully investigated. The first is that an external torque is applied to the tower and the two cables would both be indirectly excited, asymmetrically, by the torsional/oscillating tower. The two cables’ responses are the same in this case. The second is that only one of the two cables, i.e., the leader cable, is directly excited, and the other cable, i.e., the follower one, is only indirectly excited through cables–tower dynamic interactions. In such kind of leader–follower dynamics, the leader cable is quite different from the follower one. Nonlinear coupled frequency response diagrams for both systems are constructed using numerical continuation algorithms, mainly focused on the coupled steady solutions’ stabilities and bifurcations. Furthermore, the dynamic effects of tower’s moment of inertia, wing span and damping are thoroughly investigated.

     

Advanced parametric space-frequency separated representations in structural dynamics: A harmonic–modal hybrid approach

This paper is concerned with the solution to structural dynamics equations. The technique here presented is closely related to Harmonic Analysis, and therefore it is only concerned with the long-term forced response. Proper Generalized Decomposition (PGD) is used to compute space-frequency separated representations by considering the frequency as an extra coordinate. This formulation constitutes an alternative to classical methods such as Modal Analysis and it is especially advantageous when parametrized structural dynamics equations are of interest. In such case, there is no need to solve the parametrized eigenvalue problem and the space-time solution can be recovered with a Fourier inverse transform. The PGD solution is valid for any forcing term that can be written as a combination of the considered frequencies. Finally, the solution is available for any value of the parameter. When the problem involves frequency-dependent parameters the proposed technique provides a specially suitable method that becomes computationally more efficient when it is combined with a modal representation.     

A physico-mechanical approach to modeling of metal forming processes—Part II: damage analysis in processes with plastic flow of metals

A damage analysis is presented for the extrusion of a case-shaped cylindrical part by using a physico-mechanical approach for modeling metal forming processes. Two integral measures related to the hydrostatic and deviatoric parts of the damage tensor are used for the calculation of strain damage. The combined use of two damage measures in contrast to only one allows us to assess not only a risk of macro-fracture of the deformed material but also the stage of formation of large cavernous defects due to coalescence of ellipsoidal voids. Such a refined prediction of the actual quality of the material’s micro-structure is important when producing metalware that is supposed to operate under intense loading and thermal conditions. In case study of this paper the kinetic equations of damage are solved by using mutually consistent fields of stresses, flow velocities, and strains. It is shown that the predicted damage is less than its permissible value since a high hydrostatic pressure in the plastic zone heals the micro-defects, prevents their growth, and, thereby, increases the processing ductility of deformed metals during extrusion.      

Wavelet analysis on particle dynamics in a horizontal air–solid two-phase pipe flow at low air velocity

High-speed particle image velocimetry (PIV) is first used to measure two components of the fluctuating particle velocities for different particle sizes and solid mass flow rates at low air velocity in a horizontal pipe. Then, the continuous wavelet transform and orthogonal wavelet multi-resolution techniques are employed to analyze and decompose the fluctuating particle velocities to provide both quantitative and qualitative information on the particle fluctuation velocity of various frequencies. It is revealed that the fluctuating energy of axial particle velocity is mainly contributed from the wavelet components of low frequency, accounting for about 84%, near the bottom part of the pipe cross-section. However, the contribution to the fluctuating energy of vertical particle velocity accounts for about 82% from the wavelet components of high frequency. The auto-correlation analysis suggests a quasi-periodical large-scale axial particle fluctuating velocity. On the other hand, the spatial correlation analysis indicates that the low-frequency components of the axial particle velocity exhibit a large correlation near the bottom part of the pipe cross-section. From the probability density function (PDF) distribution, it is found that the low-frequency components of the axial particle velocity exhibit larger fluctuation, and this fluctuation reduces as the frequencies increase near the bottom part of the pipe cross-section. Near the top part of the pipe cross-section, however, a larger fluctuating axial particle velocity appears in the high-frequency range.     

Diffusion-based coarse graining in hybrid continuum–discrete solvers: Theoretical formulation and a priori tests

Coarse graining is an important ingredient in many multi-scale continuum–discrete solvers such as CFD–DEM (computational fluid dynamics–discrete element method) solvers for dense particle-laden flows. Although CFD–DEM solvers have become a mature technique that is widely used in multiphase flow research and industrial flow simulations, a flexible and easy-to-implement coarse graining algorithm that can work with CFD solvers of arbitrary meshes is still lacking. In this work, we proposed a new coarse graining algorithm for continuum–discrete solvers for dense particle-laden flows based on solving a transient diffusion equation. Via theoretical analysis we demonstrated that the proposed method is equivalent to the statistical kernel method with a Gaussian kernel, but the current method is much more straightforward to implement in CFD–DEM solvers. A priori numerical tests were performed to obtain the solid volume fraction fields based on given particle distributions, the results obtained by using the proposed algorithm were compared with those from other coarse graining methods in the literature (e.g., the particle centroid method, the divided particle volume method, and the two-grid formulation). The numerical tests demonstrated that the proposed coarse graining procedure based on solving diffusion equations is theoretically sound, easy to implement and parallelize in general CFD solvers, and has improved mesh-convergence characteristics compared with existing coarse graining methods. The diffusion-based coarse graining method has been implemented into a CFD–DEM solver, the results of which are presented in a separate work.     

Neimark–Sacker bifurcation analysis and complex nonlinear dynamics in a heterogeneous quadropoly game with an isoelastic demand function

In this paper, a nonlinear quadropoly game based on Cournot model with fully heterogeneous players is established. This game extends the model introduced by Tramontana and Elsadany (Nonlinear Dyn 68:187–193, 2012) who considered a heterogeneous triopoly game with an isoelastic demand function. Here, four different types of players and potentially different marginal costs are considered. Moreover, the assumption of an isoelastic demand function increases the nonlinearity of the final four-dimensional map. The stability of the resulting discrete-time dynamical system is analyzed. The existence of Neimark–Sacker bifurcation near the Nash equilibrium point of the game is shown. Also, based on the Kuznetsov’s normal form technique for discrete-time system, the stability of the Neimark–Sacker bifurcation is also discussed which indicates that the bifurcation is supercritical. Moreover, it is shown that the Nash equilibrium point of the game undergoes period-doubling (flip) bifurcation. Furthermore, the double route to chaotic dynamics in this model, via flip bifurcations and via Neimark–Sacker bifurcation of the Nash equilibrium point, is illustrated. Coexistence of multi-chaotic attractors of the model is illustrated. Simulation tools like bifurcation diagrams, stability regions of parameters, Lyapunov exponent spectrum, phase plots and basins of attraction are used to verify the complex dynamics of the game.     

Analysis of gas–solid flow and shaft-injected gas distribution in an oxygen blast furnace using a discrete element method and computational fluid dynamics coupled model

Ironmaking using an oxygen blast furnace is an attractive approach for reducing energy consumption in the iron and steel industry. This paper presents a numerical study of gas–solid flow in an oxygen blast furnace by coupling the discrete element method with computational fluid dynamics. The model reliability was verified by previous experimental results. The influences of particle diameter, shaft tuyere size, and specific ratio (X) of shaft-injected gas (SIG) flowrate to total gas flowrate on the SIG penetration behavior and pressure field in the furnace were investigated. The results showed that gas penetration capacity in the furnace gradually decreased as the particle diameter decreased from 100 to 40 mm. Decreasing particle diameter and increasing shaft tuyere size both slightly increased the SIG concentration near the furnace wall but decreased it at the furnace center. The value of X has a significant impact on the SIG distribution. According to the pressure fields obtained under different conditions, the key factor affecting SIG penetration depth is the pressure difference between the upper and lower levels of the shaft tuyere. If the pressure difference is small, the SIG can easily penetrate to the furnace center.     

The Interface Between Fresh and Salt Groundwater in Horizontal Aquifers: The Dupuit–Forchheimer Approximation Revisited

We analyze the motion of a sharp interface between fresh and salt groundwater in horizontal, confined aquifers of infinite extend. The analysis is based on earlier results of De Josselin de Jong (Proc Euromech 143:75–82, 1981). Parameterizing the height of the interface along the horizontal base of the aquifer and assuming the validity of the Dupuit–Forchheimer approximation in both the fresh and saltwater, he derived an approximate interface motion equation. This equation is a nonlinear doubly degenerate diffusion equation in terms of the height of the interface. In that paper, he also developed a stream function-based formulation for the dynamics of a two-fluid interface. By replacing the two fluids by one hypothetical fluid, with a distribution of vortices along the interface, the exact discharge field throughout the flow domain can be determined. Starting point for our analysis is the stream function formulation. We derive an exact integro-differential equation for the movement of the interface. We show that the pointwise differential terms are identical to the approximate Dupuit–Forchheimer interface motion equation as derived by De Josselin de Jong. We analyze (mathematical) properties of the additional integral term in the exact interface motion formulation to validate the approximate Dupuit–Forchheimer interface motion equation. We also consider the case of flat interfaces, and we study the behavior of the toe of the interface. In particular, we give a criterion for finite or infinite speed of propagation.     

Aerosol–cloud–precipitation interactions: A challenging problem in regional environment and climate research

Aerosols affect clouds in two broad ways: (i) presence of more number of aerosols leads to formation of more smaller droplets, and reduces coalescence, resulting in brighter clouds that reflect more solar energy back to space, hence they contribute to cooling of the Earth's surface and (ii) numerous smaller cloud droplets tend to reduce precipitation and change the extent of cloud cover and increase cloud lifetime and albedo. One of our recent studies on aerosols over the Indo-Gangetic Plains (IGP) relative to the pristine oceans to the south of Indian Ocean showed that highly absorbing aerosols could potentially lead to the revival of active condition preceded by long break. The absorption of solar radiation by aerosols such as black carbon and desert dust produces surface cooling and local stabilization of lower atmosphere. This stability effect is overcome by the enhanced moisture convergence due to the meridional gradient of aerosol-induced heating. In some other studies, we showed association between cloud thickness and cloud to sub-cloud ratio (SCR), aerosol variability (in terms of aerosol optical depth and aerosol index) and monsoon precipitation and climate over regional scale. This paper provides an overview of some salient results that have been obtained from the studies conducted, using the ground- and space-based active and passive remote sensing techniques, at the Indian Institute of Tropical Meteorology (IITM), Pune, India in the recent decade.     

The axisymmetric contraction–expansion: the role of extensional rheology on vortex growth dynamics and the enhanced pressure drop

The flow of a polystyrene Boger fluid through axisymmetric contraction–expansions having various contraction ratios (2≤β≤8) and varying degrees of re-entrant corner curvatures are studied experimentally over a large range of Deborah numbers. The ideal elastic fluid is dilute, monodisperse and well characterized in both shear and transient uniaxial extension. A large enhanced pressure drop above that of a Newtonian fluid is observed independent of contraction ratio and re-entrant corner curvature. Streak images, laser Doppler velocimetry (LDV) and digital particle image velocimetry (DPIV) are used to investigate the flow kinematics upstream of the contraction plane. LDV is used to measure velocity fluctuation in the mean flow field and to characterize a global elastic flow instability which occurs at large Deborah numbers. For a contraction ratio of β=2, a steady elastic lip vortex is observed while for contraction ratios of 4≤β≤8, no lip vortex is observed and a corner vortex is seen. Rounding the re-entrant corner leads to shifts in the onset of the flow transitions at larger Deborah numbers, but does not qualitatively change the overall structure of the flow field. We describe a simple rescaling of the deformation rate which incorporates the effects of lip curvature and allows measurements of vortex size, enhanced pressure drop and critical Deborah number for the onset of elastic instability to be collapsed onto master curves. Transient extensional rheology measurements are utilized to explain the significant differences in vortex growth pathways (i.e. elastic corner vortex versus lip vortex growth) observed between the polystyrene Boger fluids used in this research and polyisobutylene and polyacrylamide Boger fluids used in previous contraction flow experiments. We show that the role of contraction ratio on vortex growth dynamics can be rationalized by considering the dimensionless ratio of the elastic normal stress difference in steady shear flow to those in transient uniaxial extension. It appears that the differences in this normal stress ratio for different fluids at a given Deborah number arise from variations in solvent quality or excluded volume effects.     

The Beavers–Joseph Boundary Condition and Related Matters: A Historical and Critical Note

This note is concerned with the historical background to the Beavers–Joseph boundary condition at the interface of a porous medium and a clear fluid. Relevant papers published prior to 1975 are discussed. The merits of the alternative methodology utilizing the Brinkman equation are also discussed.     

The numerical modeling of rotor–stator rubbing in rotating machinery: a comprehensive review

Nonlinear Dynamics - The rotor–stator rubbing in rotating machinery generated as a consequence of rotor imbalance, shaft misalignment, and casing deformation is a potential threat to the...     

Particle statistics in a gas–solid coaxial strongly swirling flow: A direct numerical simulation

A direct numerical simulation of a strongly coaxial swirling particle-laden flow is conducted with reference to a previous experiment. The carrier phase is simulated as a coaxial swirling flow through a short nozzle injecting into a large container. The particle phase is carried by the primary jet, and simulated in the Lagrangian approach. The drag force, slip-shear force and slip-rotation force experienced by particles are calculated. A partial validation of the results is followed. The results are analyzed in Eulerian approach focusing on the statistical behavior of particle motion. The relative importance of the drag, slip-shear and slip-rotation forces under different Stokes numbers is indicated quantitatively. The particle velocity profiles, fluctuations, Reynolds stress, and turbulence intensity are demonstrated and analyzed respectively. An important “choke” behavior for large particles within the mainstream is found and interpreted. Additionally, the patterns of particle distribution and the helical structures of particle motion under different Stokes numbers are demonstrated qualitatively and analyzed quantitatively.     

A three-fluid approach in bipolar semiconductors with generation–recombination: constitutive laws and Onsager symmetry

If electrons (e) and holes (h) in metals or semiconductors are heated to the temperatures \(T_{e}\) and \(T_{h}\) greater than the lattice temperature \(T_{p}\), the electron–phonon interaction causes energy relaxation. In the nonuniform case, a momentum relaxation occurs as well. In view of such an application, a new model based on an asymptotic procedure for solving the kinetic equations of carriers and phonons is proposed, with generation–recombination of electrons and holes, which gives naturally the displaced Maxwellian at the leading order. After that, balance equations for the electron number, hole number, energy densities, and momentum densities are constructed, which constitute now a system of eight equations for the chemical potentials (carriers), the temperatures (carriers and phonons), and the drift velocities (carriers and phonons). In the drift-diffusion approximation the constitutive laws are derived and the Onsager relations recovered.     

The linearized Boltzmann equation with Cercignani–Lampis boundary conditions: Basic flow problems in a plane channel

A polynomial expansion procedure and the ADO (analytical discrete-ordinates) method are used to solve a collection of basic flow problems based on the linearized Boltzmann equation for rigid-sphere interactions and the Cercignani–Lampis boundary conditions with a free choice of the accommodation coefficients at each boundary. In particular, three classical problems defined by flow in a plane-parallel channel (Poiseuille, thermal-creep, and Couette flow) are solved (essentially) analytically and evaluated to a very high numerical standard. Some comparisons with known kinetic models are also reported.     

A fictitious domain finite element method for simulations of fluid–structure interactions: The Navier–Stokes equations coupled with a moving solid

The paper extends a stabilized fictitious domain finite element method initially developed for the Stokes problem to the incompressible Navier–Stokes equations coupled with a moving solid. This method presents the advantage to predict an optimal approximation of the normal stress tensor at the interface. The dynamics of the solid is governed by Newton׳s laws and the interface between the fluid and the structure is materialized by a level-set which cuts the elements of the mesh. An algorithm is proposed in order to treat the time evolution of the geometry and numerical results are presented on a classical benchmark of the motion of a disk falling in a channel.