Equations Of Motion Research Articles

Shear and normal stresses of electroosmotic magnetized physiological nanofluid on curved artery with moderate Reynolds number: application on electroshock therapy

PurposeStudying the shear stress and pressure resulting on the walls of blood vessels, especially during high-pressure cases, which may lead to the explosion or rupture of these vessels, can also lead to the death of many patients. Therefore, it was necessary to try to control the shear and normal stresses on these veins through nanoparticles in the presence of some external forces, such as exposure to some electromagnetic shocks, to reduce the risk of high pressure and stress on those blood vessels. This study aims to examines the shear and normal stresses of electroosmotic-magnetized Sutterby Buongiorno’s nanofluid in a symmetric peristaltic channel with a moderate Reynolds number and curvature. The production of thermal radiation is also considered. Sutterby nanofluids equations of motion, energy equation, nanoparticles concentration, induced magnetic field and electric potential are calculated without approximation using small and long wavelengths with moderate Reynolds numbers.Design/methodology/approachThe Adomian decomposition method solves the nonlinear partial differential equations with related boundary conditions. Graphs and tables show flow features and biophysical factors like shear and normal stresses.FindingsThis study found that when curvature and a moderate Reynolds number are present, the non-Newtonian Sutterby fluid raises shear stress across all domains due to velocity decay, resulting in high shear stress. Additionally, modest mobility increases shear stress across all channel domains. The Sutterby parameter causes fluid motion resistance, which results in low energy generation and a decrease in the temperature distribution.Originality/valueEquations of motion, energy equation, nanoparticle concentration, induced magnetic field and electric potential for Sutterby nano-fluids are obtained without any approximation i.e. the authors take small and long wavelengths and also moderate Reynolds numbers.

Analyzing the dynamics of a charged rotating rigid body under constant torques

This study explores the dynamical rotary motion of a charged axisymmetric spinning rigid body (RB) under the effect of a gyrostatic moment (GM). The influence of transverse and invariable body fixed torques (IBFTs), and an electromagnetic force field, is also considered. Euler’s equations of motion (EOM) are utilized to derive the regulating system of motion for the problem in a suitable formulation. Due to the lack of torque exerted along the spin axis and the nearly symmetrical nature of the RB, the spin rate is nearly unchanged. Assuming slight angular deviations of the spin axis relative to a fixed direction in space, it is possible to derive approximate analytical solutions (AS) in closed form for the attitude, translational, and rotational movements. These concise solutions that are expressed in complex form are highly effective in analyzing the maneuvers performed by spinning RBs. The study focuses on deriving the AS for various variables including angular velocities, Euler’s angles, angular momentum, transverse displacements, transverse velocities, axial displacement, and axial velocity. The graphical simulation of the subsequently obtained solutions is presented to show their precision. Furthermore, the positive impacts that alterations in the body’s parameters have on the motion’s behavior are presented graphically. The corresponding phase plane curves, highlighting the influence of different values in relation to the electromagnetic force field, the GM, and the IBFTs are drawn to analyze the stability of the body’s motion. This study has a significant role in various scientific and engineering disciplines. Its importance lies in its ability to optimize mechanical systems, explain celestial motion, and enhance spacecraft performance.

Shape stability of a microbubble in a power–law liquid

AbstractThe onset of non-spherical oscillations of a microbubble in an unbounded power–law liquid, important for biomedical ultrasound applications, is studied. Two sets of evolution equations are obtained from the equation of motion: a Rayleigh Plesset-type equation for the spherical oscillations and an equation for the non-spherical oscillations. The non-spherical oscillations are modeled using the perturbation method via the Legendre polynomials. Two kinds of instabilities, namely parametric and Rayleigh-Taylor instabilities, are investigated. A higher power–law index causes the damping of the oscillations for both spherical and non-spherical oscillations. The power–law index damping effect depends on the ultrasonic drive frequency. At natural frequency, the amplitude of the perturbations is high compared to the non-resonant cases. At a low consistency index, the damping effect of the power–law index decreases. Unlike Newtonian liquids, the viscosity of power–law liquids is affected by the frequency of the acoustic field, thereby affecting Rayleigh-Taylor instability.

Development of equation of motion deciphering locomotion including omega turns of Caenorhabditis elegans.

Locomotion is a fundamental behavior of Caenorhabditis elegans (C. elegans). Previous works on kinetic simulations of animals helped researchers understand the physical mechanisms of locomotion and the muscle-controlling principles of neuronal circuits as an actuator part. It has yet to be understood how C. elegans utilizes the frictional forces caused by the tension of its muscles to perform sequenced locomotive behaviors. Here, we present a two-dimensional rigid body chain model for the locomotion of C. elegans by developing Newtonian equations of motion for each body segment of C. elegans. Having accounted for friction-coefficients of the surrounding environment, elastic constants of C. elegans, and its kymogram from experiments, our kinetic model (ElegansBot) reproduced various locomotion of C. elegans such as, but not limited to, forward-backward-(omega turn)-forward locomotion constituting escaping behavior and delta-turn navigation. Additionally, ElegansBot precisely quantified the forces acting on each body segment of C. elegans to allow investigation of the force distribution. This model will facilitate our understanding of the detailed mechanism of various locomotive behaviors at any given friction-coefficients of the surrounding environment. Furthermore, as the model ensures the performance of realistic behavior, it can be used to research actuator-controller interaction between muscles and neuronal circuits.

Development of equation of motion deciphering locomotion including omega turns of Caenorhabditis elegans

Locomotion is a fundamental behavior of Caenorhabditis elegans (C. elegans). Previous works on kinetic simulations of animals helped researchers understand the physical mechanisms of locomotion and the muscle-controlling principles of neuronal circuits as an actuator part. It has yet to be understood how C. elegans utilizes the frictional forces caused by the tension of its muscles to perform sequenced locomotive behaviors. Here, we present a two-dimensional rigid body chain model for the locomotion of C. elegans by developing Newtonian equations of motion for each body segment of C. elegans. Having accounted for friction-coefficients of the surrounding environment, elastic constants of C. elegans, and its kymogram from experiments, our kinetic model (ElegansBot) reproduced various locomotion of C. elegans such as, but not limited to, forward-backward-(omega turn)-forward locomotion constituting escaping behavior and delta-turn navigation. Additionally, ElegansBot precisely quantified the forces acting on each body segment of C. elegans to allow investigation of the force distribution. This model will facilitate our understanding of the detailed mechanism of various locomotive behaviors at any given friction-coefficients of the surrounding environment. Furthermore, as the model ensures the performance of realistic behavior, it can be used to research actuator-controller interaction between muscles and neuronal circuits.

Period-doubling cascade route to chaos in an initially curved microbeam resonator exposed to fringing-field electrostatic actuation

This paper explores the chaotic dynamics of a piezoelectrically laminated initially curved microbeam resonator subjected to fringing-field electrostatic actuation, for the first time. The resonator is fully clamped at both ends and is coated with two piezoelectric layers, encompassing both the top and bottom surfaces. The nonlinear motion equation which is obtained by considering the nonlinear fringing-field electrostatic force, includes geometric nonlinearities due to the mid-plane stretching and initial curvature. The motion equation is discretized using Galerkin method and the reduced order system is numerically integrated over the time for the time response. The variation of the first three natural frequencies with respect to the applied electrostatic voltage is determined and the frequency response curve is obtained using the combination of shooting and continuation methods. The bifurcation points have been examined and their types have been clarified based on the loci of the Floquet exponents on the complex plane. The period-doubled branches of the frequency response curves originating from the period doubling (PD) bifurcation points are stablished. It's demonstrated that the succession PD cascades leads to chaotic behavior. The chaotic behavior is identified qualitatively by constructing the corresponding Poincaré section and analyzing the response's associated frequency components. The bifurcation diagram is obtained for a wide range of excitation frequency and thus the exact range in which chaotic behavior occurs for the system is determined. The chaotic response of the system is regularized and controlled by applying an appropriate piezoelectric voltage which shifts the frequency response curve along the frequency axis.

A new higher-order shear deformation theory for flutter instability of a hybrid sandwich panel resting on a viscoelastic foundation

Magnetorheological elastomers (MREs) are smart materials characterized by their ability to adjust mechanical properties through the manipulation of applied magnetic fields, facilitating rapid and reversible tunability. The incorporation of reinforcing fibers introduces distinctive materials termed Magnetorheological Elastomer Composites (MRECs), combining tunability with enhanced rigidity. This study delves into the flutter analysis of a sandwich panel with a mid-layer composed of laminated MRECs and face sheets composed of functionally graded materials (FGMs) featuring porosity. The sandwich panel is situated on a viscoelastic foundation and exposed to a supersonic air current. Employing a higher-order (hyperbolic) shear deformation theory (HSDT), the displacement field of the sandwich plate is estimated. This selection eliminates the requirement for shear correction factors. The mechanical characteristics of the three sheets are determined using the Halpin-Tsai micromechanical method and the rule of mixture. The First-order Piston theory is employed to evaluate the effect of aerodynamic pressure and airflow damping on the sandwich panel. The derivation of the motion equation is achieved through Hamilton's principle. In order to solve the motion equation of the structure, the Finite Element Method (FEM) is utilized. The investigation systematically explores the influence of diverse parameters on the flutter boundary of the sandwich panel.

Numerical modeling of turbulent air-water mixture flows in spilling breaking waves

This paper presents a numerical study of turbulent air-water mixture flows under spilling breaking waves by using a mixture-flow model. The model is extended from the earlier 2D numerical framework NEWFLUME, which solves the Reynolds-averaged Navier-Stokes (RANS) equations for mean flow motions. The turbulence closure is accomplished by a modified k−ε two-equations model, in which the effect of entrained bubbles on turbulence production and dissipation is considered by incorporating additional buoyancy terms. A transient equation based on two-fluid formulation is solved for air bubble transport. The present model is validated by comparing the simulated free surface displacement, void fraction, mean velocity and turbulence kinetic energy against existing experimental data and the previous numerical results. The effects of the additional convective term in transient equation of void fraction and the bubble-induced turbulence term in turbulence transport equation are further explored. The results reveal that the additional convective term in the transient equation of void fraction plays an essential role of suppressing unphysical turbulence growth and bringing correct amount of air entrainment into water during wave breaking processes. The additional bubble-induced turbulence production and dissipation terms become effective after wave breaks when the void fraction has been appropriately modeled, and they can further reduce the simulated turbulence intensities in spilling breaking wave, rendering a better agreement with the experimental measurements.

Design of a sliding mode controller for suppressing coupled axial & torsional vibrations in horizontal drill strings using Extended Kalman Filter

Nowadays, drill string technologies are among the leading fields in the industries of oil and gas extraction and are widely used in the construction of subway tunnels. Due to their significant applications, horizontal drill strings have attracted the attention of researchers. Various vibrations that occur during the drilling process are also sources of system uncertainty. To solve the mentioned problems, in this paper, a coupled axial-torsional model is considered for a horizontal drill string, and all the forces and torques that are involved in the drilling process, are modeled. All the primary modes of the axial and torsional vibrations are calculated using a numerical solution. For solving the PDE equations of motion, the mode summation method is used, and in order to determine the minimum necessary modes, mode convergence analysis has been accomplished for this coupled axial-torsional plant for the first time. In the next step, a novel nonlinear sliding mode controller (SMC) is designed for this plant to eliminate the coupled axial and torsional vibrations, using only two control inputs. As the proposed controller requires the state variables, and the sensors are incapable of recording all the state variables, an Extended Kalman Filter (EKF) is designed to estimate all the 12 state variables, which is one of the novelties of this research. Eventually, a comparison analysis with a conventional SMC and a Pole-Placement controller is provided to ensure the performance of the proposed controller. Moreover, unmodelled dynamics and parametric uncertainties are considered to ensure the robustness of the proposed controller.

Higher Time-Derivative Theories from Space–Time Interchanged Integrable Field Theories

We compare a relativistic and a nonrelativistic version of Ostrogradsky’s method for higher-time derivative theories extended to scalar field theories and consider as an alternative a multi-field variant. We apply the schemes to space–time rotated modified Korteweg–de Vries systems and, exploiting their integrability, to Hamiltonian systems built from space–time rotated inverse Legendre transformed higher-order charges of these systems. We derive the equal-time Poisson bracket structures of these theories, establish the integrability of the latter theories by means of the Painlevé test and construct exact analytical period benign solutions in terms of Jacobi elliptic functions to the classical equations of motion. The classical energies of these partially complex solutions are real when they respect a certain modified CPT-symmetry and complex when this symmetry is broken. The higher-order Cauchy and initial-boundary value problem are addressed analytically and numerically. Finally, we provide the explicit quantization of the simplest mKdV system, exhibiting the usual conundrum of having the choice between having to deal with either a theory that includes non-normalizable states or spectra that are unbounded from below. In our non-Hermitian system, the choice is dictated by the correct sign in the decay width.

Revealing the inherent running mechanism of quadruped mammals based on a novel bionic stiffness model

In this paper, the limb of a goat is chosen as the research object, and according to mammalian anatomy, a bionic model called the quasi inverted pendulum with “J” curve spring (QIPJCS) model with nonlinear stiffness is built, and the equations of motion are derived. Based on these equations, the advantages of the QIPJCS model are illustrated from the aspect of the stable motion region by the SFA (step-to-fall analysis) numerical simulation method. These results are compared with the traditional SLIP model. Furthermore, the ARM (Apex-Return-Map) of this model is built, and the fixed points are analyzed. Finally, according to the locomotion law of goats running with gallop gaits and the analysis of the dead-point support effect, the dynamic motion mechanism of goat limbs is elucidated, and the equivalent mechanism model is built. Based on the mechanism, the dynamic mechanical analysis indicates that the joint driving torque can be minimized to conserve energy by optimizing the landing angle. The running mechanism research of quadruped mammals, which is based on the novel bionic stiffness model, provides theoretical support for the design of high-performance mechanical legs and the motion control of bionic robots.

A variational approach to deriving Hamiltonian functions under non-holonomic constraints — the case of the Gaussian thermostats

The Hamiltonian function (H), if it exists, provides a powerful platform for studying the dynamics of constrained many-particle systems. For systems under holonomic constraints, H can be derived from the principle of least action (PLA) formulated in terms of a suitable set of generalized coordinates. In molecular dynamics (MD) simulations, where the system temperature often prescribes a non-holonomic constraint, the challenge is to find the system’s Hamiltonian. The Gaussian Isokinetic (GIK) thermostat, the first deterministic physics-preserving thermostat in MD, was derived from the principle of least constraint (PLC) in the early 1980s: its H, interestingly, was discovered almost 15 years later. A mathematically rigorous connection of the GIK Hamiltonian with PLA, however, has been missing until now. In this article, we adopt a variational approach to establish the equivalence of PLC with PLA for a class of non-holonomically constrained many-particle systems, and present the general formulation of the Euler–Lagrange equations of motion through time rescaling. Making use of the particular forms of the Lagrange multiplier obtained from the principle of least constraint, we derive H for the GIK and the Gaussian Isoenergetic thermostats.

Geometrically nonlinear vibration of toroidal sandwich shells with auxetic honeycomb core under periodic/impulsive pressure

The main objective of this paper is to evaluate the dynamic behaviour of sandwich toroidal shells composed of a core with negative Poisson ratio (auxetic honeycomb material) and nanocomposite enriched coating layers under impulsive/periodic pressures considering the geometrical nonlinearity. The set of nonlinear equilibrium equations based on the first-order shear deformation theory are and the resultant nonlinear differential motion equations are approximately solved by the Galerkin method and 4th-order Runge–Kutta solution, respectively. The achieved responses show that the amplitude of outward deformations is intensively increased when the proposed sandwich shells are subjected to rectangular periodic and impulsive loadings which shows the destructive effects of this type dynamic loading. The performed parametric studies shows that these destructive effects can be controlled by changing the geometrical properties of the honeycomb cells and the volume fractions of the nanocomposites in the coating layers. It is notable that in the previous studies the nonlinear dynamic response of toroidal shells subjected to various periodic and impulsive pressures was not investigated. However, in the present work, the mentioned dynamic responses were achieved and investigated without spending time-consuming computational. In fact, the novelty of this work is to develop such a robust and efficient tool.

Noether’s currents for conformable fractional scalar field theories

The construction of fractional derivatives with the right properties for use in field theory is reputed to be a difficult task, essentially because of the absence of a unique definition and uniform properties. The conformable fractional derivative introduced in 2014 by Khalil et al. in their seminal paper is a novel and well-behaved definition of fractional derivative for a function that is derivable in the usual sense. In this paper, we investigate the consistency of the Euler–Lagrange formalism for a field theory defined on such a fractional space–time. We especially focus on the relation between symmetries and conservation laws (Noether’s currents), about the symmetry group introduced to construct the Lagrangian of the field. In particular, we show that the use of the conformable derivative induces additional terms in the calculation of the action variation. We also investigate the conservation of the Noether current and show that this property only takes place on condition that the equations of motion are verified with a new definition of the conserved law.

Vibroacoustic analysis of submerged fluid-filled cylindrical shell

This study presents a comprehensive Chebyshev-Fourier spectral method for analyzing the free and forced vibroacoustic characteristics of submerged fluid-filled cylindrical shells in heavy fluids. The vibroacoustic interaction system of submerged cylindrical shells consists of three discrete regions: The Sanders shell theory is applied for the shell structure, whereas the wave equation is used to determine the fluid inside. The external fluid pressure on the shell is described using the Helmholtz boundary integral equation. The governing equations for the shell's motion, the formulation of the external fluid pressure on the fluid-structure-interaction (FSI) surface using first-class Chebyshev polynomial and Fourier series expansions, and the internal sound pressure employing two-dimensional Chebyshev polynomial and Fourier series expansions, are seamlessly integrated with Gauss-Chebyshev-Lobatto sampling techniques. Artificial springs are implemented to replicate the boundary conditions of the cylindrical shell. The precision of the proposed methodology is corroborated by comparisons with existing literature and finite element method (FEM) results. Furthermore, the research examines how the properties of the acoustical medium on the vibroacoustic responses of the cylindrical shell. Results indicate that the presence of air inside the shell slightly influence on its acoustic radiation. However, replacing the internal medium with water significantly alters the shell's original radiation characteristics.

Numerical Modeling of Non-Isothermal Laminar Flow and Heat Transfer of Paraffinic Oil with Yield Stress in a Pipe

This paper presents the results of a study on the non-isothermal laminar flow and heat transfer of oil with Newtonian and viscoplastic rheologies. Heat exchange with the surrounding environment leads to the formation of a near-wall zone of viscoplastic fluid. As the flow proceeds, the transformation of a Newtonian fluid to a viscoplastic state occurs. The rheology of the Shvedoff–Bingham fluid as a function of temperature is represented by the effective molecular viscosity apparatus. A numerical solution to the system of equations of motion and heat transfer was obtained using the Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. The calculated data are obtained at Reynolds number Re from 523 to 1046, Bingham number Bn from 8.51 to 411.16, and Prandl number Pr = 45. The calculations’ novelty lies in the appearance of a “stagnation zone” in the near-wall zone and the pipe cross-section narrowing. The near-wall “stagnation zone” is along the pipe’s radius from r/R = 0.475 to r/R = 1 at Re = 523, Bn = 411.16, Pr = 45, u1 = 0.10 m/s, t1 = 25 °C, and tw = 0 °C. The influence of the heat of phase transition of paraffinic oil on the development of flow and heat transfer characteristics along the pipe length is demonstrated.

Frequency Regulation in Nanocomposite Sandwich Joined Conical-Conical Shells with Open Cell Foam Partially Supported by Pasternak Elastic Foundation Using GDQ

This study investigated the natural vibration of a sandwich conical-conical shell that is partially supported by Winkler–Pasternak foundations. The sandwich system comprises an open-cell foam (OCF) core and laminated composite face sheets reinforced with graphene platelets using functionally graded models (FG-GPLRC). To account for through-the-thickness shear deformations and rotary inertias, the first-order theory of shells is combined with Donnell-type kinematic assumptions. Hamilton’s principle is utilized to establish the general motion equations and related boundary and continuity conditions. The resulting system of equations is discretized using the semi-analytical generalized differential quadrature (GDQ) approach. An eigenvalue problem is formulated to analyze the corresponding mode shapes and vibration frequencies for the shell ends and continuity criteria under various boundary conditions. Parametric experiments are conducted for sandwich conical-conical shells partially supported by Winkler–Pasternak foundations. This study examined the impact of multiple factors on the frequency of a joined shell structure, including the total thickness, FG-GPL face sheet thickness, various boundary conditions, the length of each cone segment, and the Winkler and shear elastic foundation. Notably, this research also analyzed the effects of a partial elastic foundation on the structure’s frequency for the first time.

TSDT vs CPT and FSDT for Free Vibration Analysis of Functionally Graded Incompressible Plates

This article studies vibrational behavior of incompressible functionally graded plates through utilizing classical, first-order shear, and third-order shear deformation plate theories. Plate material properties are presumed to differ continuously in the thickness direction concurringto a power law function. The motion and continuity equations are produced using three different plate theories, and then they are analytically solved for rectangular plates with simple supports utilizingNavier's method. The results are verified by obtaining the plate natural frequency for the power law index value equals zero and comparing them to those reported in previous works. It is shown that transverse vibrational analysis of the classical and first-order shear deformationplate theories for incompressible plates are not as precise as thecompressible ones. It is shown that this issue is due to the fact that according to those theories, unlike the higher order theories, the hydrostatic pressure cannot participate in carrying the bending loads. Consequently, the equivalent flexural stiffness and as a result flexural frequencies decrease. So, to analyze the vibrational behavior of functionally graded incompressibleplates, whether the plate is either thin or thick, higher order theories should be adopted. Also, it is demonstrated that TSDT is the simplest shear deformation plate theory for which the hydrostatic pressure can contribute to withstandig the bending moments. So, it can be a practical theory for free vibration analysis of functionally graded incompressible plates. Finally, this theory has been taken into consideration to analyze the vibrational behavior of rectangular plates made of functionally graded incompressible materials Also detailed parametric studies have been carried out.

Studying the nonlinear response of incompressible hyperelastic thin circular cylindrical shells with geometric imperfections

This study presents a comprehensive analysis of hyperelastic thin cylindrical shells exhibiting initial geometrical imperfections. The nonlinear equations of motion are derived using an improved formulation of Donnell’s nonlinear shallow-shell theory and Lagrange’s equations, incorporating the small strain hypothesis. Mooney–Rivlin constitutive model is employed to capture the hyperelastic behavior of the material. The coupled nonlinear equations of motion are analytically solved using Multiple-Scale method, which effectively accounts for the inherent nonlinearity of the system. To ensure the model’s accuracy, the linear model is verified by comparing the results with those obtained through hybrid finite element method. Subsequently, the model with only geometrical nonlinearity is evaluated against other research works existing in the open literature to ensure its reliability and precision. Finally, the results of the model, considering both geometrical and physical nonlinearity, are verified against the results obtained from Abaqus software. The main objective of this research is to provide a detailed understanding of the response of hyperelastic thin cylindrical shells in the presence of initial geometric imperfections. In this order, the impact of three distinct geometric imperfections – axisymmetric, asymmetric, and a combination of driven and companion modes – on the natural frequency is examined. The behavior of each of these geometric imperfections is investigated by varying their respective coefficients. The numerical results indicate that geometric imperfections enhance the natural frequency, and employing different models for imperfections leads to a variation in this trend. In the amplitude response of hyperelastic cylindrical shells, two peaks coexist, reflecting the softening and hardening responses of the system. Distinct initial geometric imperfections influence these two peaks.

Nonlinear thermoelastic wave propagation in general FGM sandwich rectangular plates

The present work is dedicated to investigating the thermoelastic wave propagation behavior of sandwich rectangular plates (SRP) made of functionally graded material (FGM). The main contribution lies in the partial modification of basic theoretical expressions and solution methods to improve the accuracy of practical system models. An analytical model with three types of general configurations is established. The porosity distribution in FGM layers depends on the degree of mixture of the constituent materials, with the FGM layers without porosity taken as a reference model. The effect of porosity within FGMs is addressed through a refined analytical formulation of material properties, and the temperature-dependent material properties of FGM sandwich structures (FGMSS) maintain continuity through the thickness. This improved framework introduces a porosity function encompassing three distinct structural and geometrical functions: the core-to-thickness ratio (CTR), porosity volume fraction (PVF), and porosity distribution function (PDF). It is worth mentioning that the theoretical expressions maintain good continuity and reliability under the influence of thermal conditions and system parameters of the proposed structures. Furthermore, considering the generation of thermal strain energy (TSE) caused by thermal expansion of the structure in the normal direction, an improved analytical approach for determining TSE in rectangular plate structures is then investigated by introducing the Green's nonlinear strain (GNS). Hamilton's principle is applied to derive the wave motion equations and analytical solutions for the wave dispersion relations are derived. Furthermore, accurate numerical simulation is performed and the solution is verified with data available in published resources. In addition, we present a systematic parametric analysis to examine the effects of porosity, configuration, power-law exponent (PLE), PVF, CTR, temperature, and wave number on the thermoelastic wave propagation behavior of FGMSRP.