Finite ramp time correction method for non-linear viscoelastic material model

In this paper, a simple formula for the prediction of stress (strain) in the relaxation (creep) period is derived for a non-linear viscoelastic material model which takes into account the finite ramp time. Usually, it is assumed that the ramp time is small and, therefore, loading can be described via Heaviside step function. This assumption, when applied to the material parameters identification process, can lead to a large errors in the values of the approximated material parameters. Especially, for the materials which undergo significant stress decay in the beginning of relaxation the assumption of infinite small ramp time can induce severe errors. With the help of the derived formula more reliable material parameter identification can be accomplished. The proposed method is tested with numerical simulations and compared with analytical results, Heaviside step loading case and method described by Nordin and Varna. Simulations show good agreement with analytical results.     

Excitation symmetry analysis method (ESAM) for calculation of higher order non-linearities

A non-linear system with 3rd order non-linearity is fully characterized using symmetry analysis (SA) applied to the excitation, as it is done in 2nd order non-linear systems using the pulse inverted method. Symmetry analysis is performed using irreductible representations and the character table of C₃ rotation point group, which leads to the construction of three eigen-excitations allowing extraction of the 3rd order non-linearity parameter without the perturbation of fundamental and 2nd order terms. Validation of this concept is based on excitation symmetry analysis method (ESAM) which was tested on simulated noisy signals and compared with classical spectral analysis.     

Three-dimensional discrete element models for the granular statics and dynamics of powders in cavity filling

Rapid granular flow from a moving container and angle of repose formation were investigated by numerical simulations using the discrete element method (DEM) and experiments. Grain models of various geometrical complexity were studied and their ability to reproduce the experiments in those regimes was explored. The predictive power of the most realistic model for gravity driven cavity filling was assessed. Good agreement between computed and measured density distributions within the filled cavities provides a basis for numerical process variations aiming at homogenized density distributions. The effect of numerical coarse graining was found to be negligible for all properties of interest provided that force laws are scaled properly and corrections for boundary effects are taken into account. The proposed scaling was tested for a certain set of force laws but could be applied to different DEM forces as well. An analytic mass flow law for powder discharge from a moving container was derived and verified by our DEM simulations.     

Magnetic torque attitude control of a satellite using the state-dependent Riccati equation technique

A non-linear attitude control method for a satellite with magnetic torque rods using the state-dependent Riccati equation (SDRE) technique has been developed. The magnetic torque caused by the interaction with the Earth's magnetic field and the magnetic moment of torque rods plays a role of the control torque. The detailed equations of motion for this system are presented using angular velocity and quaternions. The SDRE controller is developed for the non-linear systems which can be formed in pseudo-linear representations using the state-dependent coefficient (SDC) method without linearization procedure. The aim of this control system is to achieve a stable attitude within 5°, and minimize the control effort. The stability regions for the SDRE controlled satellite system are estimated through the investigation of the stability conditions developed for pseudo-linear systems and the application of Lyapunov's theorem. For comparisons, the Linear Quadratic Regulator (LQR) method using the solution of the algebraic Riccati equation (ARE) is also applied to this non-linear system. The performance of the SDRE non-linear control system demonstrates more robustness and stability than the LQR control system when subjected to a wide range of initial conditions.     

Modeling and simulation of martensitic phase transitions with a triple point

A framework for modeling complex global energy landscapes in a piecewise manner is presented. Specifically, a class of strain-dependent energy functions is derived for the triple point of Zirconia (ZrO₂), where tetragonal, orthorhombic (orthoI) and monoclinic phases are stable. A simple two-dimensional framework is presented to deal with this symmetry breaking. An explicit energy is then fitted to the available elastic moduli of Zirconia in this two-dimensional setting. First, we use the orbit space method to deal with symmetry constraints in an easy way. Second, we introduce a modular (piecewise) approach to reproduce or model elastic moduli, energy barriers and other characteristics independently of each other in a sequence of local steps. This allows for more general results than the classical Landau theory (understood in the sense that the energy is a polynomial of invariant polynomials). The class of functions considered here is strictly larger. Finite-element simulations for the energy constructed here demonstrate the pattern formation in Zirconia at the triple point.     

Non-linear principal resonance of an orthotropic and magnetoelastic rectangular plate

Based on the von Karman plate theory of large deflection, we have derived a non-linear partial differential equation for the vibration of a thin orthotropic plate under the combined action of a transverse magnetic field and a transverse harmonic mechanical load. The influence of the magnetic field is due to the magnetic Lorentz force induced by the eddy current. By employing the Bubnov-Galerkin method, the non-linear partial differential equation is transformed into a third-order non-linear ordinary differential equation. The amplitude-frequency equations are further derived by means of the multiple-scale method. As numerical examples for an orthotropic plate made of silver, the influence of the magnetic field, orthotropic material property, plate thickness, and the mechanical load on the principal resonance behavior is investigated. The higher-order effect and stability of the solution are also discussed.     

Time reversal and non-linear elastic wave spectroscopy (TR NEWS) techniques

Non-linear elastic wave spectroscopy (NEWS) has been shown to exhibit a high degree of sensitivity to both distributed and isolated non-linear scatterers in solids. In the case of an isolated non-linear scatterer such as a crack, by combining the elastic energy localization of time reversal (TR) with NEWS, it is shown that one can isolate non-linear scatterers in solids. The experiments reviewed here present two distinct methods of combining TR and NEWS for this purpose. The techniques each have there own advantages and disadvantages, with respect to each other and other non-linear methods, which are discussed.     

Convection in a Porous Medium with Variable Gravity Field and Magnetic Field Effects

The problem of convection in a variable gravity field with magnetic field effect is studied using methods of linear instability theory and non-linear energy theory. Then, the accuracies of both the linear instability and global non-linear energy stability thresholds are tested using a three-dimensional simulation. The strong stabilizing effect of gravity field and magnetic field is shown. Moreover, the results support the assertion that the linear theory is very accurate in predicting the onset of convective motion, and thus regions of stability.     

Non-linear vibration analysis of an elastic plate subjected to heavy fluid loading in magnetic field

In the present study, the non-linear vibration of an elastic plate subjected to heavy fluid loading in an inclined magnetic field is investigated. The structural non-linearity, fluid non-linearity, and the effects of magnetic field are all incorporated in the formulations to derive the governing equation of the plate. The method of multiple scales is adopted to determine the eigenvalues and mode shapes of the linear vibration, and then the amplitude of the non-linear vibration response of the plate is calculated. Based on the assumptions of ordering and formulations of multiple scales, it can be concluded that the linear dynamic behavior of the plate under heavy fluid loading but weak near-resonant loading is influenced by the effects of the fluid loading, linear structural rigidity and linear magnetic field, furthermore, the non-linear dynamic behavior of the plate under heavy fluid loading but weak near-resonant loading is dominated and controlled by the effects of the fluid loading, non-linear structural rigidity and non-linear magnetic field. Both thick and thin plates are investigated; the contributions due to the structural non-linearity and acoustic linear radiation damping are of the same order for a rather thick plate. For a thin plate, the structural non-linearity completely controls the behavior of the plate, which implies that in this case the effect of fluid loading is considerably negligible. In general, it can be concluded that both the effects of magnetic field and structural non-linearity play important roles only on the first few modes of the plate.     

Chaotic diffusion in steady wavy vortex flow—Dependence on wave state and correlation with Eulerian symmetry measures

The dimensionless effective axial diffusion coefficient, D_z, calculated from particle trajectories in steady wavy vortex flow in a narrow gap Taylor-Couette system, has been determined as a function of Reynolds number (R=Re/Re_c), axial wavelength (λ_z), and the number of azimuthal waves (m). Two regimes of Reynolds number were found: (i) when R<3.5, D_z has a complex and sometimes multi-modal dependence on Reynolds number; (ii) when R>3.5, D_z decreases monotonically.Eulerian quantities measuring the departure from rotational symmetry, ?_θ, and flexion-free flow, ?_ν, were calculated. The space-averaged quantities and were found to have, unlike D_z, a simple unimodal dependence on R. In the low R regime the correlation between D_z and ?_θ?_ν was complicated and was attributed to variations in the spatial distribution of the wavy disturbance occurring in this range of R. In the large R regime, however, the correlation simplified to for all wave states, and this was attributed to the growth of an integrable vortex core and the concentration of the wavy disturbance into narrow regions near the outflow and inflow jets.A reservoir model of a wavy vortex was used to determine the rate of escape across the outflow and inflow boundaries, the size of the ‘escape basins’ (associated with escape across the outflow and inflow boundaries), and the size of the trapping region in the vortex core. In the low R regime after the breakup of all KAM tori, the outflow basin (γ_O) is larger than the inflow basin (γ_I), and both γ_O and γ_I are (approximately) independent of R. In the large R regime, with increasing Reynolds number the trapping region grows, the outflow basin decreases, and the inflow basin shows a slight increase. This implies that the growth of the integrable core occurs at the expense of the outflow escape basin. Finally, it is shown that the variation of the weighted escape rates (γ_Or_O,γ_Ir_I) with Reynolds number was in excellent qualitative agreement with the variation of .     

Non-linear dynamics of a soft magneto-elastic Cartesian manipulator

This paper studies the non-linear dynamics of a soft magneto-elastic Cartesian manipulator with large transverse deflection. The system has been subjected to a time varying magnetic field and a harmonic base excitation at the roller-supported end. Unlike elastic and viscoelastic manipulators, here the governing temporal equation of motion contains additional two frequency forced, and linear and non-linear parametric excitation terms. Method of multiple scales has been used to solve the temporal equation of motion. The influences of various system parameters such as amplitude and frequency of magnetic field strength, amplitude and frequency of support motion, and the payload on the frequency response curves have been investigated for three different resonance conditions. With the help of numerical results, it has been shown that by using suitable amplitude and frequency of magnetic field, the vibration of the manipulator can be significantly controlled. The developed results and expressions can find extensive applications in the feed-forward vibration control of the flexible Cartesian manipulator using magnetic field.     

Wilton ripples between two uniform streaming magnetic fluids

An analysis is made of the small-amplitude capillary-gravity waves which occur on the interface of two incompressible inviscid magnetic fluids of different densities. The waves arise as a result of second harmonic resonance. The fluids moving with uniform velocities parallel to their interface are stressed by an oblique magnetic field. The linear relations between the oblique magnetic field and the instability criteria of the linear waves are analyzed. At the stability region (away from the neutral curve) of the linear theory, a pair of coupled non-linear partial differential equations are presented. On the neutral curve, a pair of coupled non-linear partial differential equations are introduced. The last pair of equations may be regarded as the counterparts of the single Klein-Gordon equation which occurs in the non-resonant case. In all cases, the wave profile and its stability conditions are obtained. These conditions are discussed analytically and graphically.     

Stability of conducting liquid flowing down an inclined plane at moderate Reynolds number in the presence of constant electromagnetic field

The stability of a conducting viscous film flowing down an inclined plane at moderate Reynolds number in the presence of electromagnetic field is investigated under induction-free approximation. Using momentum integral method a non-linear evolution equation for the development of the free surface is derived. The linear stability analysis of the evolution equation shows that the magnetic field stabilizes the flow whereas the electric field stabilizes or destabilizes the flow depending on its orientation with the flow. The weakly non-linear study reveals that both the supercritical stability and subcritical instability are possible for this type of thin film flow. The influence of magnetic field on the different zones is very significant, while the impact of electric field is very feeble in comparison.     

Electromagnetically conducting elastic plates in a magnetic field: modeling and dynamic implications

The basic field equations and boundary conditions necessary for the dynamic approach of electromagnetically conducting flat plates subjected to an external magnetic field are derived.Whereas the structural equations include the geometrical non-linearities of elastic isotropic plates, the electromagnetic equations are used in a linearized form that is obtained from their non-linear counterpart by applying the small disturbance concept. In this context, it was shown that the governing equations involve the bending-stretching coupling arising from both the geometrical non-linear approach of the problem and the inclusion of the Lorentz ponderomotive forces that are reduced to a 2-D plate counterpart.A number of special cases are considered, implications of the external magnetic field on non-linear/linear eigenfrequencies are highlighted, and pertinent conclusions are outlined.     

On non-linear magnetohydrodynamic problems of an Oldroyd 6-constant fluid

In this work we develop a mathematical model to predict the velocity profile for an unidirectional, incompressible and steady flow of an Oldroyd 6-constant fluid. The fluid is electrically conducting by a transverse magnetic field. The developed governing equation is non-linear. This equation is solved analytically to obtain the general solution. The governing non-linear equation is also solved numerically subject to appropriate boundary conditions (three cases of typical plane shearing flows) by an iterative technique with the finite-difference discretizations. A parametric study of the physical parameters involved in the problems such as the applied magnetic field and the material constants is conducted. The obtained results are illustrated graphically to show salient features of the solutions. Numerical results show that the applied magnetic field tends to reduce the flow velocity. Depending on the choice of the material parameters, the fluid exhibits shear-thickening or shear-thinning behaviours.     

Effects of magnetohydrodynamic flow past a vertical plate with variable surface temperature

An analysis is performed to study the magnetohydrodynamic flow of an electrically conducting, viscous incompressible fluid past a semi-infinite vertical plate with variable surface temperature under the action of transversely applied magnetic field. The heat due to viscous dissipation and the induced magnetic field are assumed to be negligible. The dimensionless governing equations are unsteady, two-dimensional, coupled and non-linear governing equations. It is found that the magnetic field parameter has a retarding effect on the velocities of air and water.     

Similarity analysis in magnetohydrodynamics: Hall effects on free convection flow and mass transfer past a semi-infinite vertical flat plate

Heat and mass transfer along a semi-infinite vertical flat plate under the combined buoyancy force effects of thermal and species diffusion is investigated in the presence of a strong non-uniform magnetic field and the Hall currents are taken into account. The induced magnetic field due to the motion of the electrically conducting fluid is negligible. This assumption is valid for a small magnetic Reynolds number. The similarity solutions are obtained using the scale group of transformations. These are the only symmetry transformations admitted by the field equations. The non-linear boundary layer equations with the boundary conditions are transferred to a system of non-linear ordinary differential equations with the appropriate boundary conditions. Furthermore, the similarity equations are solved numerically by using a fourth order Runge-Kutta scheme with the shooting method. Numerical results for the velocity profiles, the temperature profiles and the concentration profiles are presented graphically for various values of the magnetic parameter M in the range of 0-1 with the Hall parameter m taking the values 0.5, 1, 2, and 3.     

MHD free convective flow and mass transfer over a stretching sheet with chemical reaction

The effect of chemical reaction on free convective flow and mass transfer of a viscous, incompressible and electrically conducting fluid over a stretching surface is investigated in the presence of a constant transverse magnetic field. The non-linear boundary layer equations with the boundary conditions are transferred by a similarity transformation into a system of non-linear ordinary differential equations with the appropriate boundary conditions. Furthermore, the similarity equations are solved numerically by using a fourth order Runge-Kutta scheme with the shooting method. Numerical results of the skin friction coefficient, the local Nusselt number Nu, the local Sherwood number Sh, as will as the velocity, temperature and concentration profiles are presented for gases with a Prandtl number of 0.71 for various values of chemical reaction parameter, order of reaction, magnetic parameter and Schmidt number.     

Non-linear based time reversal acoustic applied to crack detection: Simulations and experiments

In non-destructive evaluation, non-linear elastic wave spectroscopy (NEWS) methods represent powerful tools to explore damaged zones in a sample. The combination of these methods with time reversal (TR) process can be used to either increase the stress on the retrofocusing position or to retrofocuse elastic waves on the defect, when only the non-linear components of the received signal are time reversed. In this paper, we propose two different approaches to detect cracks in metallic samples in coupling NEWS methods and TR process. The first one, NEWS-TR, is defined by sending back only the non-linear components which are preliminary time reversed. Two different techniques to filter non-linear components have been numerically studied: classical harmonic filter and pulse inversion. In these two cases, performances of retrofocusing on different defect positions are analyzed and compared. The second approach, TR-NEWS, consists in measuring the non-linear signature on the focal spot. An experimental study of parametric interaction between two reversed signals (f₁ and f₂) is led. Measurements of components at f₂-f₁ and f₂+f₁ around a crack are performed and discussed.