Fractional telegraph equation under moving time-harmonic impact.

The time-fractional telegraph equation with moving time-harmonic source is considered on a real line. We investigate two characteristic versions of this equation: the "wave-type" with the second and Caputo fractional time-derivatives as well as the "heat-type" with the first and Caputo fractional time-derivatives. In both cases the order of fractional derivative 1 < α < 2. For the time-fractional telegraph equation it is impossible to consider the quasi-steady-state corresponding to the solution being a product of a function of the spatial coordinate and the time-harmonic term. The considered problem is solved using the integral transforms technique. The solution to the "wave-type" equation contains wave fronts and describes the Doppler effect contrary to the solution for the "heat-type" equation. Numerical results are illustrated graphically for different values of nondimensional parameters.

Spontaneous Negative Entropy Increments in Granular Flows.

The entropy inequality, commonly taken as an axiom of continuum mechanics, is found to be spontaneously violated in macroscopic granular media undergoing collisional dynamics. The result falls within the fluctuation theorem of nonequilibrium thermodynamics, which is known to replace the Second Law for finite systems. This phenomenon amounts to the system stochastically displaying negative increments of entropy. The focus is on granular media in Couette flows, consisting of monosized circular disks (with 10 to 104 disks of diameters 0.01 m to 1 m) with frictional-Hookean contacts simulated by molecular dynamics accounting for micropolar effects. Overall, it is determined that the probability of negative entropy increments diminishes with the Eulerian velocity gradient increasing, while it tends to increase in a sigmoidal fashion with the Young modulus of disks increasing. This behavior is examined for a very wide range of known materials: from the softest polymers to the stiffest (i.e., carbyne). The disks' Poisson ratio is found to have a weak effect on the probability of occurrence of negative entropy increments.

Thermo-poromechanics of fractal media.

This article advances continuum-type mechanics of porous media having a generally anisotropic, product-like fractal geometry. Relying on a fractal derivative, the approach leads to global balance laws in terms of fractal integrals based on product measures and, then, converting them to integer-order integrals in conventional (Euclidean) space. Proposed is a new line transformation coefficient that is frame invariant, has no bias with respect to the coordinate origin and captures the differences between two fractal media having the same fractal dimension but different density distributions. A continuum localization procedure then allows the development of local balance laws of fractal media: conservation of mass, microinertia, linear momentum, angular momentum and energy, as well as the second law of thermodynamics. The product measure formulation, together with the angular momentum balance, directly leads to a generally asymmetric Cauchy stress and, hence, to a micropolar (rather than classical) mechanics of fractal media. The resulting micropolar model allowing for conservative and/or dissipative effects is applied to diffusion in fractal thermoelastic media. First, a mechanical formulation of Fick's Law in fractal media is given. Then, a complete system of equations governing displacement, microrotation, temperature and concentration fields is developed. As a special case, an isothermal model is worked out. This article is part of the theme issue 'Advanced materials modelling via fractional calculus: challenges and perspectives'.

Impact force and moment problems on random mass density fields with fractal and Hurst effects.

This paper reports the application of cellular automata to study the dynamic responses of Lamb-type problems for a tangential point load and a concentrated moment applied on the free surface of a half-plane. The medium is homogeneous, isotropic and linear elastic while having a random mass density field with fractal and Hurst characteristics. Both Cauchy and Dagum random field models are used to capture these effects. First, the cellular automata approach is tested on progressively finer meshes to verify the code against the continuum elastodynamic solution in a homogeneous continuum. Then, the sensitivity of wave propagation on random fields is assessed for a wide range of fractal and Hurst parameters. Overall, the mean response amplitude is lowered by the mass density field's randomness, while the Hurst parameter (especially, for ß < 0.2) is found to have a stronger influence than the fractal dimension on the response. The resulting Rayleigh wave is modified more than the pressure wave for the same random field parameters. Additionally, comparisons with previously studied Lamb-type problems under normal in-plane and anti-plane loadings are given. This article is part of the theme issue 'Advanced materials modelling via fractional calculus: challenges and perspectives'.

Elastic-plastic-brittle transitions and avalanches in disordered media.

A spring lattice model with the ability to simulate elastic-plastic-brittle transitions in a disordered medium is presented. The model is based on bilinear constitutive law defined at the spring level and power-law-type disorder introduced in the yield and failure limits of the springs. The key parameters of the proposed model effectively control the disorder distribution, significantly affecting the stress-strain response, the damage accumulation process, and the fracture surfaces. The model demonstrates a plastic strain avalanche behavior for perfectly plastic as well as hardening materials with a power-law distribution, in agreement with the experiments and related models. The strength of the model is in its generality and ability to interpolate between elastic-plastic hardening and elastic-brittle transitions.

Heat conduction in live tissue during radiofrequency electrosurgery.

In this paper, we propose a method to model radiofrequency electrosurgery to capture the phenomena at higher temperatures and present the methods for parameter estimation. Experimental data taken from our surgical trials performed on in vivo porcine liver show that a non-Fourier Maxwell-Cattaneo-type model can be suitable for this application when used in combination with an Arrhenius-type model that approximates the energy dissipation in physical and chemical reactions. The resulting model structure has the advantage of higher accuracy than existing ones, while reducing the computation time required.

Minimally Invasive Live Tissue High-fidelity Thermophysical Modeling using Real-time Thermography.

We present a novel thermodynamic parameter estimation framework for energy-based surgery on live tissue, with direct applications to tissue characterization during electrosurgery. This framework addresses the problem of estimating tissue-specific thermodynamics in real-time, which would enable accurate prediction of thermal damage impact to the tissue and damage-conscious planning of electrosurgical procedures. Our approach provides basic thermodynamic information such as thermal diffusivity, and also allows for obtaining the thermal relaxation time and a model of the heat source, yielding in real-time a controlled hyperbolic thermodynamics model. The latter accounts for the finite thermal propagation time necessary for modeling of the electrosurgical action, in which the probe motion speed often surpasses the speed of thermal propagation in the tissue operated on. Our approach relies solely on thermographer feedback and a knowledge of the power level and position of the electrosurgical pencil, imposing only very minor adjustments to normal electrosurgery to obtain a high-fidelity model of the tissue-probe interaction. Our method is minimally invasive and can be performed in situ. We apply our method first to simulated data based on porcine muscle tissue to verify its accuracy and then to in vivo liver tissue, and compare the results with those from the literature. This comparison shows that parameterizing the Maxwell-Cattaneo model through the framework proposed yields a noticeably higher fidelity real-time adaptable representation of the thermodynamic tissue response to the electrosurgical impact than currently available. A discussion on the differences between the live and the dead tissue thermodynamics is also provided.

Minimally Invasive Live Tissue High-Fidelity Thermophysical Modeling Using Real-Time Thermography.

We present a novel thermodynamic parameter estimation framework for energy-based surgery on live tissue, with direct applications to tissue characterization during electrosurgery. This framework addresses the problem of estimating tissue-specific thermodynamics in real-time, which would enable accurate prediction of thermal damage impact to the tissue and damage-conscious planning of electrosurgical procedures. Our approach provides basic thermodynamic information such as thermal diffusivity, and also allows for obtaining the thermal relaxation time and a model of the heat source, yielding in real-time a controlled hyperbolic thermodynamics model. The latter accounts for the finite thermal propagation time necessary for modeling of the electrosurgical action, in which the probe motion speed often surpasses the speed of thermal propagation in the tissue operated on. Our approach relies solely on thermographer feedback and a knowledge of the power level and position of the electrosurgical pencil, imposing only very minor adjustments to normal electrosurgery to obtain a high-fidelity model of the tissue-probe interaction. Our method is minimally invasive and can be performed in situ. We apply our method first to simulated data based on porcine muscle tissue to verify its accuracy and then to in vivo liver tissue, and compare the results with those from the literature. This comparison shows that parameterizing the Maxwell-Cattaneo model through the framework proposed yields a noticeably higher fidelity real-time adaptable representation of the thermodynamic tissue response to the electrosurgical impact than currently available. A discussion on the differences between the live and the dead tissue thermodynamics is also provided.

Blunt head impact causes a temperature rise in the brain.

A magnetic resonance imaging-based finite-element model is employed to assess the temperature in the human brain due to blunt head trauma. The model is based on a coupled thermoelasticity under small strain and Fourier or Maxwell-Cattaneo heat conduction assumptions, accompanied by a standard coupling of thermal fields to mechanics. It is found that mechanical impacts on the forehead cause a temperature rise of up to 0.3°C above the reference homogeneous temperature field.

Infinite-Dimensional Adaptive Boundary Observer for Inner-Domain Temperature Estimation of 3D Electrosurgical Processes using Surface Thermography Sensing.

We present a novel 3D adaptive observer framework for use in the determination of subsurface organic tissue temperatures in electrosurgery. The observer structure leverages pointwise 2D surface temperature readings obtained from a real-time infrared thermographer for both parameter estimation and temperature field observation. We introduce a novel approach to decoupled parameter adaptation and estimation, wherein the parameter estimation can run in real-time, while the observer loop runs on a slower time scale. To achieve this, we introduce a novel parameter estimation method known as attention-based noise-robust averaging, in which surface thermography time series are used to directly estimate the tissue's diffusivity. Our observer contains a real-time parameter adaptation component based on this diffusivity adaptation law, as well as a Luenberger-type corrector based on the sensed surface temperature. In this work, we also present a novel model structure adapted to the setting of robotic surgery, wherein we model the electrosurgical heat distribution as a compactly supported magnitude- and velocity-controlled heat source involving a new nonlinear input mapping. We demonstrate satisfactory performance of the adaptive observer in simulation, using real-life experimental ex vivo porcine tissue data.

Averaging of turbulent micropolar media: turbulent couple-stress, heat flux, and energy.

The equations governing turbulent flow of micropolar incompressible media are studied using the Reynolds decomposition. The average force-stress is augmented by the Reynolds stress of the same form as in classical continuum mechanics, while the average couple-stress is augmented by a new turbulent couple-stress. Additionally, the average heat flux and internal energy density are modified from the expressions known in classical turbulent fluids. On this basis, the entropy inequality is examined in both classical and micropolar continuum settings.

Doppler effect described by the solutions of the Cattaneo telegraph equation.

The Cattaneo telegraph equation for temperature with moving time-harmonic source is studied on the line and the half-line domain. The Laplace and Fourier transforms are used. Expressions which show the wave fronts and elucidate the Doppler effect are obtained. Several particular cases of the considered problem including the heat conduction equation and the wave equation are investigated. The quasi-steady-state solutions are also examined for the case of non-moving time-harmonic source and time-harmonic boundary condition for temperature.

Mach Fronts in Random Media with Fractal and Hurst Effects.

An investigation of transient second sound phenomena due to moving heat sources on planar random media is conducted. The spatial material randomness of the relaxation time is modeled by Cauchy or Dagum random fields allowing for decoupling of fractal and Hurst effects. The Maxwell-Cattaneo model is solved by a second-order central differencing. The resulting stochastic fluctuations of Mach wedges are examined and compared to unperturbed Mach wedges resulting from the heat source traveling in a homogeneous domain. All the examined cases are illustrated by simulation movies linked to this paper.

Violations of the Clausius-Duhem inequality in Couette flows of granular media.

Spontaneous violations of the Clausius-Duhem (CD) inequality in Couette-type collisional flows of model granular media are studied. Planar systems of monosized circular discs (with disc numbers from 10 to 204, and disc diameters from 0.001 m to 1 m) with frictional-Hookean contacts are simulated under periodic boundary conditions by a molecular dynamics. The scale-dependent homogenization of micropolar media is used to determine the energy balances and mechanical entropy production. The dissipation function exhibits spontaneous negative entropy increments described by the fluctuation theorem. The boundary between violations and non-violations of the CD inequality is mapped in the parameter space, where the probability of such events diminishes with the disc diameter, the disc number and the area fraction increasing. The dissipation function is a random process, tending to Gaussian as the number of discs increases, and possessing non-trivial fractal and anti-persistent Hurst properties.

Finite Element Methods in Human Head Impact Simulations: A Review.

Head impacts leading to traumatic brain injury (TBI) present a major health risk today, projected to become the third leading cause of death by 2020. While finite element (FE) models of the human brain are important tools to understand and mitigate TBI, many unresolved issues remain that need to be addressed to improve these models. This work aims to provide readers with background information regarding the current state of research in this field as well as to present recent advancements made possible by improvements to computational resources. Specifically, this has manifested as a drive to introduce more details in FE models in the form of increased spatial resolution and improved material models such as nonlinear and anisotropic constitutive models. The need to work with high-resolution FE meshes is underlined by the dominant wavelengths involved in transient pressure and shear wave propagation and the ability to model the brain surface. We also discuss improvements to experimental validation techniques which allow for better calibrated models. We review these recent developments in detail, highlighting their contributions to the field as well as identifying open issues where more research is needed.

Electrostatic and magnetostatic properties of random materials.

Scale dependence of electrostatic and magnetostatic properties is investigated in the setting of spatially random linear lossless materials with statistically homogeneous and spatially ergodic random microstructures. First, from the Hill-Mandel homogenization conditions adapted to electric and magnetic fields, uniform boundary conditions are formulated for a statistical volume element (SVE). From these conditions, there follow upper and lower mesoscale bounds on the macroscale (effective) electrical permittivity and magnetic permeability. Using computational electromagnetics methods, these bounds are obtained through numerical simulations for composites of two types: (i) two-dimensional (2D) random checkerboard (two-phase) microstructures and (ii) analogous 3D random (two-phase) media. The simulation results demonstrate a scale-dependent trend of these bounds toward the properties of a representative volume element (RVE). This transition from SVE to RVE is described using a scaling function dependent on the mesoscale δ, the volume fraction v_{f}, and the property contrast k between two phases. The scaling function is calibrated through fitting the data obtained from extensive simulations (∼10000) conducted over the aforementioned parameter space. The RVE size of a given microstructure can be estimated down to within any desired accuracy using this scaling function as parametrized by the contrast and the volume fraction of two phases.

Heat conduction in porcine muscle and blood: experiments and time-fractional telegraph equation model.

This paper presents experimental evidence for the damped-hyperbolic nature of transient heat conduction in porcine muscle tissue and blood. An examination of integer order and Maxwell-Cattaneo heat conduction models indicates that the latter, in effect resulting in a time-fractional telegraph (TFT) equation, provides the best fit to transient heat phenomena in such materials. The numerical method is verified on Dirichlet and Neumann initial boundary value problems using existing analytical results. Overall, the TFT equation captures the wave-like nature of heat conduction and temperature profiles obtained in experiments, while reducing the need for further tunable parameters.

Effect of cerebrospinal fluid modeling on spherically convergent shear waves during blunt head trauma.

The MRI-based computational model, previously validated by tagged MRI and harmonic phase imaging analysis technique on in vivo human brain deformation, is used to study transient wave dynamics during blunt head trauma. Three different constitutive models are used for the cerebrospinal fluid: incompressible solid elastic, viscoelastic, and fluid-like elastic using an equation of state model. Three impact cases are simulated, which indicate that the blunt impacts give rise not only to a fast pressure wave but also to a slow, and potentially much more damaging, shear (distortional) wave that converges spherically towards the brain center. The wave amplification due to spherical geometry is balanced by damping due to tissues' viscoelasticity and the heterogeneous brain structure, suggesting a stochastic competition of these 2 opposite effects. It is observed that this convergent shear wave is dependent on the constitutive property of the cerebrospinal fluid, whereas the peak pressure is not as significantly affected.

Scaling of slip avalanches in sheared amorphous materials based on large-scale atomistic simulations.

Atomistic simulations of binary amorphous systems with over 4 million atoms are performed. Systems of two interatomic potentials of the Lennard-Jones type, LJ12-6 and LJ9-6, are simulated. The athermal quasistatic shearing protocol is adopted, where the shear strain is applied in a stepwise fashion with each step followed by energy minimization. For each avalanche event, the shear stress drop (Δσ), the hydrostatic pressure drop (Δσ_{h}), and the potential energy drop (ΔE) are computed. It is found that, with the avalanche size increasing, the three become proportional to each other asymptotically. The probability distributions of avalanche sizes are obtained and values of scaling exponents fitted. In particular, the distributions follow a power law, P(ΔU)∼ΔU^{-τ}, where ΔU is a measure of avalanche sizes defined based on shear stress drops. The exponent τ is 1.25±0.1 for the LJ12-6 systems, and 1.15±0.1 for the LJ9-6 systems. The value of τ for the LJ12-6 systems is consistent with that from an earlier atomistic simulation study by Robbins et al. [Phys. Rev. Lett. 109, 105703 (2012)]PRLTAO0031-900710.1103/PhysRevLett.109.105703, but the fitted values of other scaling exponents differ, which may be because the shearing protocol used here differs from that in their study.

Frequency-dependent scaling from mesoscale to macroscale in viscoelastic random composites.

This paper investigates the scaling from a statistical volume element (SVE; i.e. mesoscale level) to representative volume element (RVE; i.e. macroscale level) of spatially random linear viscoelastic materials, focusing on the quasi-static properties in the frequency domain. Requiring the material statistics to be spatially homogeneous and ergodic, the mesoscale bounds on the RVE response are developed from the Hill-Mandel homogenization condition adapted to viscoelastic materials. The bounds are obtained from two stochastic initial-boundary value problems set up, respectively, under uniform kinematic and traction boundary conditions. The frequency and scale dependencies of mesoscale bounds are obtained through computational mechanics for composites with planar random chessboard microstructures. In general, the frequency-dependent scaling to RVE can be described through a complex-valued scaling function, which generalizes the concept originally developed for linear elastic random composites. This scaling function is shown to apply for all different phase combinations on random chessboards and, essentially, is only a function of the microstructure and mesoscale.