
According to the traditional method of high intensity focused ultrasound (HIFU) treatment, the acoustic and thermal characteristic parameters of constant temperature (room temperature or body temperature) are used to predict thermal lesion. Based on the nonlinear spherical beam equation (SBE) and Pennes bio-heat transfer equation, and a new acoustic-thermal coupled model is proposed. The constant and temperature-dependent acoustic and thermal characteristic parameters are used to predict thermal lesion, and the predicted lesion area are compared with each other. Moreover, the relationship between harmonic amplitude ratio (P2/P1) and thermal lesion is studied. Combined with the known experimental data of acoustic and thermal characteristic parameters of biological tissue and data fitting method, the relationship between acoustic and thermal characteristic parameters and temperature is obtained; and the thermal lesion simulation calculation is carried out by using the acoustic and thermal characteristic parameters under constant temperature and temperature- dependent acoustic and thermal characteristic parameters, respectively. The simulation results show that under the same irradiation condition, the thermal lesion predicted by temperature-dependent acoustic and thermal characteristic parameters is larger than that predicted by traditional method, and the thermal lesion increases with the decrease of harmonic amplitude ratio.
Citation: Hu Dong, Gang Liu, Xin Tong. Influence of temperature-dependent acoustic and thermal parameters and nonlinear harmonics on the prediction of thermal lesion under HIFU ablation[J]. Mathematical Biosciences and Engineering, 2021, 18(2): 1340-1351. doi: 10.3934/mbe.2021070
[1] | Yan Ling Fu, Wei Zhang . Some results on frames by pre-frame operators in Q-Hilbert spaces. AIMS Mathematics, 2023, 8(12): 28878-28896. doi: 10.3934/math.20231480 |
[2] | Gang Wang . Some properties of weaving K-frames in n-Hilbert space. AIMS Mathematics, 2024, 9(9): 25438-25456. doi: 10.3934/math.20241242 |
[3] | Sergio Verdú . Relative information spectra with applications to statistical inference. AIMS Mathematics, 2024, 9(12): 35038-35090. doi: 10.3934/math.20241668 |
[4] | Cure Arenas Jaffeth, Ferrer Sotelo Kandy, Ferrer Villar Osmin . Functions of bounded (2,k)-variation in 2-normed spaces. AIMS Mathematics, 2024, 9(9): 24166-24183. doi: 10.3934/math.20241175 |
[5] | Chibueze C. Okeke, Abubakar Adamu, Ratthaprom Promkam, Pongsakorn Sunthrayuth . Two-step inertial method for solving split common null point problem with multiple output sets in Hilbert spaces. AIMS Mathematics, 2023, 8(9): 20201-20222. doi: 10.3934/math.20231030 |
[6] | Osmin Ferrer Villar, Jesús Domínguez Acosta, Edilberto Arroyo Ortiz . Frames associated with an operator in spaces with an indefinite metric. AIMS Mathematics, 2023, 8(7): 15712-15722. doi: 10.3934/math.2023802 |
[7] | Abdullah Ali H. Ahmadini, Amal S. Hassan, Ahmed N. Zaky, Shokrya S. Alshqaq . Bayesian inference of dynamic cumulative residual entropy from Pareto Ⅱ distribution with application to COVID-19. AIMS Mathematics, 2021, 6(3): 2196-2216. doi: 10.3934/math.2021133 |
[8] | Messaoud Bounkhel . V-Moreau envelope of nonconvex functions on smooth Banach spaces. AIMS Mathematics, 2024, 9(10): 28589-28610. doi: 10.3934/math.20241387 |
[9] | Jamilu Adamu, Kanikar Muangchoo, Abbas Ja'afaru Badakaya, Jewaidu Rilwan . On pursuit-evasion differential game problem in a Hilbert space. AIMS Mathematics, 2020, 5(6): 7467-7479. doi: 10.3934/math.2020478 |
[10] | Yasir Arfat, Muhammad Aqeel Ahmad Khan, Poom Kumam, Wiyada Kumam, Kanokwan Sitthithakerngkiet . Iterative solutions via some variants of extragradient approximants in Hilbert spaces. AIMS Mathematics, 2022, 7(8): 13910-13926. doi: 10.3934/math.2022768 |
According to the traditional method of high intensity focused ultrasound (HIFU) treatment, the acoustic and thermal characteristic parameters of constant temperature (room temperature or body temperature) are used to predict thermal lesion. Based on the nonlinear spherical beam equation (SBE) and Pennes bio-heat transfer equation, and a new acoustic-thermal coupled model is proposed. The constant and temperature-dependent acoustic and thermal characteristic parameters are used to predict thermal lesion, and the predicted lesion area are compared with each other. Moreover, the relationship between harmonic amplitude ratio (P2/P1) and thermal lesion is studied. Combined with the known experimental data of acoustic and thermal characteristic parameters of biological tissue and data fitting method, the relationship between acoustic and thermal characteristic parameters and temperature is obtained; and the thermal lesion simulation calculation is carried out by using the acoustic and thermal characteristic parameters under constant temperature and temperature- dependent acoustic and thermal characteristic parameters, respectively. The simulation results show that under the same irradiation condition, the thermal lesion predicted by temperature-dependent acoustic and thermal characteristic parameters is larger than that predicted by traditional method, and the thermal lesion increases with the decrease of harmonic amplitude ratio.
Quantitative, informative and applicable models have long been proven beneficial to assist law enforcement to curb urban crimes. Following the seminal work (see reference [18]) on the mathematics of agent-based models for residential burglary, many works have been done on mathematical crime modeling and prediction, see e.g. references [1], [3], [7], [9], [12], [13], [14], [15], [16], [18], [19] and [20], and the references cited therein, and in [18].
There are roughly two classes of crime models. One is agent-based aiming to describe the individual activities, and the other is event-based without specifying the individual agents, aiming to predict the patterns of observed events (see e.g. references [15] and [16]). Here we address the first class of models.
In reference [18], "hotspots", clusters of residential burglary well documented in real life, are quantitatively studied for the first time. A discrete lattice grid is imposed, upon which residents are distributed each labeled with its level of attractiveness, and burglary agents walk over sites seeking for targets. Burglary dynamics is coupled with the environment variable based on the repeat and near-repeat victimization and the broken-windows effects. These notions maintain that a previously burgled house and its neighbors are more attractive to potential offenders, as an environment encouraging further illegal activities is more likely to be created by past crimes with visible signs (see the relevant references in references [18], [19] and [20]).
In reference [18] the time steps are set as deterministic and all types of events occur at fixed regular intervals. Random arrivals are incorporated later in references [19] and [20], governed by stochastic clocks generated from Poisson processes. In reference [19], there is one uniform Poisson-clock, and in reference [20] independent Poisson clocks are assigned to each individual agent. To summarize, these agent-based models are referred to as DTS (deterministic-timestep) Model (see reference [18]), SSRB (stochastic-statistical residential burglary) Model (see reference [20]), and SSRB-IPC (independent Poisson-clock) Model (see reference [20]), respectively. Continuum analogues parallel to these models are derived as a mean-field limit in reference [18] or a potential hydrodynamic limit in references [19] and [20], both of which lead to the same coupled deterministic reaction-diffusion equation system.
The models mentioned above serve as the first-generation agent-based criminal behavior models, whose prototype will be referred to as First-Generation Model. They are stochastic models with a continuum limit that can be viewed as a statistical average. The continuum limit is simpler to analyze and simulate by computer, and when total population is sufficiently large, the discrete and continuum simulations exhibit similar hotspot dynamics. However, as agent number decreases, the discrete simulations exhibit more transient hotspot regimes, generating the finite size effects. Quantitative analysis in references [19] and [20] show that these effects arise due to stochasticity of the system, and criminals act randomly in reality. It is called for that we combine the strengths of First-Generation Model and its continuum limit to build a second-generation criminal behavior model.
In this work, we construct a stochastic multi-scale criminal behavior model with hybrid dynamics, where agent actions and environment variables are set on different spatio-temporal scales. This model will be referred to as M-IPC (Multiscale-Independent-Poisson-Clock) Model. In contrast, the clocks in First-Generation Model (deterministic and random) all have rates at the same order of magnitude. Here we assume instead that the environment variables evolve as the fast component continuously in time and in space, while criminal agents move along discrete lattice grids following independent Poisson clocks on slow time scales. The evolution of the environment variables can be sociologically interpreted as the spreading of information, which we assume to be in a continuous mode. That is, the rate of change of information at each spatial point is "fast" while the amount of change is infinitesimal during each infinitesimally short time period. And we adopt suitable temporal-spatial scaling to make the change of information at each spatial point over a time period to have the same order of magnitude. The result is a partial differential equation whose parameter depends upon agent distributions, which switches values according to the Poisson jumps governing agent activities. The simulation cost of agent-based M-IPC Model is vastly reduced compared to the SSRB-IPC model. Moreover, transient hotspot regimes as well as stationary hotspot regimes arise in the simulations. And our M-IPC Model with the multiscale setup can also be applied to types of crimes other than residential burglary, e.g. drug smugglings.
The separation of scales not only facilitates simulations, but also brings in the mathematical framework of piecewise deterministic Markov processes (PDMP), while First-Generation Model evolves according to (discrete or continuous time) Markov processes. A martingale approach similar to that developed in references [19] and [20] is applicable. The martingale formulation can be used to represent the model as the sum of a deterministic and a stochastic part, which will be useful for further study of hotspot dynamics. PDMP, also referred to as the iterated random functions (see e.g. reference [6]), hybrid stochastic processes, or jumping Markov processes (see e.g. reference [11]), has been extensively studied (see e.g. references [2], [5], [8], [10] and [17], and the references cited therein). Also PDMP has been widely applied in physics, chemistry, biology and other fields in natural science (see e.g. references [2], [8] and [17], and the references cited therein). However, as far as we know, this is the first time that PDMP is applied to crime modeling, and our methodology can be useful for quantitative study in social and behavioral science in general.
The paper is organized as follows. We first introduce the set-up of agent-based M-IPC Model (Section 2), where the slow and the fast components are described in Sections 2.2 and 2.3, respectively. The computer simulations are shown in Section 3, where both transient and stationary hotspot regimes are displayed. We review the general characterizations of PDMP and apply them to our model (Section 4.1). Finally the martingale representation of our model is established in Section 4.2.
We start with a domain
1With minor changes we can also consider e.g. the Dirichlet boundary conditions, which is more realistic.
(ⅰ) Natural districts corresponding to those in an urban setting. For example, in Los Angeles, we can consider the natural counties as what those districts correspond to.
(ⅱ) Artificial districts to facilitate the collecting of data etc.
Attached to each
ˉAK=∫DKA(x)dx, | (1) |
We also assume that
A(x,t)=B(x,t)+A0(x), | (2) |
where
(NK(t)),B(x,0))=(N0K,B0(x)). | (3) |
Following SSRB-IPC Model in reference [20], three types of independent Poisson clocks govern the time increments of arrivals of events, and all the Poisson clocks are independent with each other. Unlike the setup in First-Generation Model (see references [18], [19] and [20]), here we assume that the attractiveness evolves over
(Ⅰ) Committing crimes.
A Type (Ⅰ) Poisson clock is assigned to each agent to govern his action of committing crimes like burgling. Suppose that a Type (Ⅰ) clock of an agent in district
For agent activity Type (Ⅰ) described as above, there is an alternative description. Suppose that each criminal has a Poisson clock with rate
(Ⅱ) Moving along districts.
A Type (Ⅱ) Poisson clock is assigned to each agent to govern his movement. Suppose that a Type (Ⅱ) clock of an agent advances at time
qK→J(t)=ˉAJ(t−)TK(t−), | (4) |
where
δt≅1c2L2. | (5) |
(Ⅲ) Replacement.
A Type (Ⅲ) Poisson clock is assigned to each district
Remark 1. The Poisson clock is particularly suitable to model random arrivals and has been used to build Poisson processes which are widely applied in biology, economics, and physics (for references see those cited in references [19] and [20]).
Following the setup of First-Generation Model (see references [18], [19] and [20]), we assume that the attractiveness field react to the agent activities described above according to the repeat and near-repeat victimization, and the broken-windows effects (see the related references in references [7], [18], [19] and [20]).
To model the repeat victimization and broken windows effects, we assume that the attractiveness increases continuously in time, with an increasing rate depending upon the average number of criminal events at the current location. However, this increase has a finite lifetime and decays at a certain speed. And we model the near-repeat victimization effect by letting the attractiveness to spread diffusively in space. Summing up, between the Poisson jumps of the criminal activities, the following partial different equation describes the spatio-temporal evolution of
dB(x,t)dt=c1θNK(t)A(x,t)−ωB(x,t)−ηΔB(x,t),x∈DK,K∈D. | (6) |
Here
For
N(x,t)=∑K∈DNK1K(t)(x), | (7) |
where
dA(t,x)dt=θc1N(x,t)A(x,t)−ωB(x,t)+ηΔB(x,t). | (8) |
We note that to completely specify the coupled dynamics of criminal activities and attractiveness field evolution, it is sufficient to keep track of the total number of criminals in each district. To conclude, the scales of the criminal activity and the evolution of the attractiveness field separate automatically.
We simulate M-IPC Model described above under various combinations of the parameters, in order to give insight into the behavior regime of the model, and compare it with the DTS, SSRB, and SSRB-IPC Models. The resulting attractiveness plots are displayed in Figs. 1 and 2.
Behavioral regimes similar to those in references [18], [19] and [20] for the attractiveness field
Case 1 Stationary hotspots. In this regime, the system tends towards an almost steady state in which stationary spots of high attractiveness are found, surrounded by areas of extremely low attractiveness.
Case 2 Dynamic hotspots. In this regime, loclized spots of increased attractiveness form and are transient. These spots remain for varying lengths of time, and may appear and disappear at seemingly random locations, and move about in space and deform over time. There are two sub categories, depending on whether the hotspots interact between each other or not.
Case 2.1. Independent hotspots. Hotspots each change in size and shape but any two hotspots rarely have interactions.
Case 2.2. Interactive hotspots. Hotspots repeatedly and constantly merge with each other and split. This regime reveals the intrinsic randomness of the model.
All the simulations were run with
In Fig. 1 (a), we set
The parameter regimes are robust for all three cases of hotspot dynamics mentioned above, and the same hotspot dynamics is also observed over other random paths with the same parameters used to create the plots in Fig. 1. In Fig. 2, distinct-stochastic-basis realizations for Fig. 1(c) are displayed.
Remark 2. The numerical scheme for (6) is similar to that in reference [18]. We use a semi-implicit scheme in which
[1+ωδt−ηδtΔ]Am+1=Am+θc1δtNmAm. | (9) |
Here
The coupled dynamics of
First introduced as in reference [4], a PDMP, roughly speaking, is a dynamic system occasionally interrupted by a pure jump process whose states then switch according to certain probability distributions (see e.g. references [2], [5], [8] and [10]).
We start with a standard stochastic basis
dU(t)dt=L(U(t)). | (10) |
Alternatively, starting from
P(τ1>t|U(0)=U0)=exp(−∫t0Λ(U(s))ds). | (11) |
Then
The infinitesimal generator of the stochastic process can be derived following e.g. references [4], [5], [8], [10] and [17] as
Tf(u)=Λ(u)[Qf(u)−∫Ef(u)Q(u;dv)]+L(u)∘∇f(u), | (12) |
for every
f(U(t))=f(U(0))+∫t0Tf(U(s))ds+M(t), | (13) |
where
For every
⟨N(t),ϕ⟩=∫MN(x,t)ϕ(x)dx=L−2∑K∈DNK(t)ˉϕK, | (14) |
where
ˉϕK:=∫DKϕ(x)dx. |
We define
⟨(N(t),B(t)),ϕ⟩:=(⟨N(t),ϕ⟩,⟨B(t),ϕ⟩). | (15) |
Then
Theorem 4.1. Starting with the initial data
{⟨N(t),ϕ⟩=⟨N0,ϕ⟩+∫t0Λ(N(s),B(s))(Q1−I1)⟨(N(s),B(s)),ϕ⟩ds+M1(⟨(N(t),B(t)),ϕ⟩),⟨B(t),ϕ⟩=⟨B0,ϕ⟩+∫t0L2(⟨(N(s),B(s)),ϕ⟩)ds+M2(⟨(N(t),B(t)),ϕ⟩), | (16) |
where
Λ(N(t),B(t))=c2L2∑K∈DNK(t)+L2c1∑K∈DNK(t)ˉAK(t)+L2Γ, | (17) |
(Q1−I1)⟨(N(t),B(t)),ϕ⟩=c2∑K∈DNK(t)(∑K′K′∼KˉAK′(t)TK(t)(ˉϕK′−ˉϕK))c2L2∑K∈DNK(t)+L2c1∑K∈DNK(t)ˉAK(t)+L2Γ−c1∑K∈DˉAK(t)NK(t)ˉϕKc2L2∑K∈DNK(t)+L2c1∑K∈DNK(t)ˉAK(t)+L2Γ+ΓL−2∑K∈DˉϕKc2L2∑K∈DNK(t)+L2c1∑K∈DNK(t)ˉAK(t)+L2Γ, | (18) |
L2(⟨(N(t),B(t)),ϕ⟩)=∫M[θc1N(x,t)A(x,t)−ωB(x,t)+ηΔB(x,t)]ϕ(x)dx. | (19) |
In other words, the deterministic flow of the PDMP driven by
Proof of Theorem 4.1. Substituting
Tu=Λ(u)[Qu−u∫EQ(u;dv)]+L(u)=Λ(u)(Q−I)u+L(u). | (20) |
Thus we have for the vector-valued stochastic process
T⟨(N(t),B(t)),ϕ⟩=(L1(⟨(N(t),B(t)),ϕ⟩),L2(⟨(N(t),B(t)),ϕ⟩))+Λ(N(t),B(t))((Q1−I1)⟨(N(t),B(t)),ϕ⟩,(Q2−I2)⟨(N(t),B(t)),ϕ⟩), | (21) |
where
In between advancements of Poisson clocks, the value of
L1(⟨(N(t),B(t)),ϕ⟩)=d⟨N(t),ϕ⟩dt=⟨dN(t)dt,ϕ⟩=0. | (22) |
Similarly as
L2(⟨(N(t),B(t)),ϕ⟩)=d⟨B(t),ϕ⟩dt=⟨dB(t)dt,ϕ⟩, | (23) |
which together with (6) implies (19).
Since all the Poisson clocks in the system are independent,
Suppose that a Poisson clock advances at time
Λ(N(t),B(t))(Q1−I1)⟨(N(t),B(t)),ϕ⟩=limδt→01δtE[⟨B(δt+t−),ϕ⟩−⟨B(t−),ϕ⟩|(N(t−),B(t−))]=L2c1∑K∈DˉAK(t)NK(t)(−L−2ˉϕK)+c2L2∑K∈D∑K′K′∼KˉAK′(t)TK(t)(L−2NK(t)(ˉϕK′−ˉϕK))+Γ∑K∈DL−2ˉϕK=−c1∑K∈DˉAK(t)NK(t)ˉϕK+c2∑K∈D∑K′K′∼KˉAK′(t)TK(t)NK(t)(ˉϕK′−ˉϕK)+L−2Γ∑K∈DˉϕK. | (24) |
This together with (17) implies (18). From (6), we infer that the flow
limδt→0[⟨B(δt+t−),ϕ⟩−⟨B(t−),ϕ⟩]=0, | (25) |
which implies that
Λ(N(t−),B(t−))(Q2−I2)⟨(N(t−),B(t−)),ϕ⟩=limδt→01δtE[⟨B(δt+t−),ϕ⟩−⟨B(t−),ϕ⟩|(N(t−),B(t−))]=Λ(N(t−),B(t−))×0=0. | (26) |
With every term on the right-hand-side of (21) derived, from (12), (13), (20), and (21) we obtain (16).
The proof of Theorem 4.1 is completed.
Remark 3. Unlike First-Generation Model, the finite size effects are circumvented here. Nevertheless, our model still generates interesting pattern formation of hotspot dynamics. Pattern formation in complex system has been studied extensively recently in natural and social sciences (see e.g. the references cited in the introduction of [20]). And the mathematical framework that we develop here based on PDMP theory and martingale formulation can be useful for quantitative study in this area, which has been lacking so far.
Ever since the pioneering agent-based residential burglary model of criminal behavior reference [18] is published, many works have been done aiming to study and improve the original model e.g. by incorporating random arrivals as in references [19] and [20]. In these models, residents are assumed to be located on the grids of a lattice field. Two components, the agent activities over the lattice and the level of attractiveness of target sites, are assumed to interact with each other on the same spatio-temporal scale. The continuum limit of First-Generation Model turns out to be reaction-advection-diffusion equations. However in the continuum simulations, transient hotspot regime is missing, which is not realistic as transient hotspots are well-document in empirical observations. A new generation of agent-based models of criminal behavior are needed.
In this work, we propose a multiscale model, where the attractiveness evolves continuously based on a diffusive partial different equation, whose parameter changes values whenever a Poisson clock associated with an agent action advances. Compared with First-Generation Model with two slow components, our model has a fast component and easier to simulate and analyze. Moreover, the stochastic nature of the model is successfully maintained and transient hotspot regimes with intrinsic randomness are picked up in the simulations as desired.
Furthermore, our multi-scale hybrid model can be characterized as a PDMP, which is applied in such a circumstance for the first time to the best of our knowledge. A martingale formula is derived, which consists of the deterministic component corresponding to the diffusive partial differential equation of attractiveness evolution, and the stochastic component corresponding to the Poisson arrivals. Results presented here will be transformative for both elements of application and analysis of agent-based models for criminal behavior.
We would like to thank the helpful discussions with Professors A. Debussche, A. L. Bertozzi, M. B. Short, T. Liggett, C. Mueller, J. Quastel, Da Kuang, Yifan Chen, Fangbo Zhang, Yatin Chow, Hangjie Ji, Yuanzhen Shao, and Kenneth Van. Yuan Zhang has been partially supported by Beijing Academy of Artificial Intelligence (BAAI).
[1] | M. Marinova, M. Rauch, H. H. Schild, H. M. Strunk, Novel non-invasive treatment with high-intensity focused ultrasound (HIFU), Ultraschall Med., 37 (2016), 46-55. |
[2] | R. Kovatcheva, J. Vlahov, J. Stoinov, F. Lacoste, C. Ortuno, K. Zaletel, US-guided high-intensity focused ultrasound as a promising non-invasive method for treatment of primary hyperparathyroidism, Eur. Radiol., 24 (2014), 2052-2058. |
[3] | S. M. Huang, H. L. Liu, D. W. Li, M. L. Li, Ultrasonic Nakagami imaging of high-intensity focused ultrasound-induced thermal lesions in porcine livers: ex vivo study, Ultrason. imaging, 40 (2018), 310-324. |
[4] | C. Liu, Y. Zhou, Detection of gaps between high-intensity focused ultrasound (HIFU)-induced lesions using transient axial shear strain elastograms, Med. Phys., 45 (2018), 3831-3847. |
[5] | D. Zhang, S. Zhang, M. Wan, S. Wang, A fast tissue stiffness-dependent elastography for HIFU-induced lesions inspection, Ultrasonics, 51 (2011), 857-869. |
[6] | J. Zhang, S. Chauhan, Neural network methodology for real-time modelling of bio-heat transfer during thermo-therapeutic applications, Artif. Intell. Med., 101 (2019), 101728. |
[7] | J. Zhang, N. D. Bui, W. Cheung, S. K. Roberts, S. Chauhan, Fast computation of desired thermal dose: Application to focused ultrasound-induced lesion planning, Numer. Heat Transfer, Part A, 77 (2020), 666-682. |
[8] | Q. Tan, X. Zou, H. Dong, Y. Ding, X. Zhao, Influence of blood vessels on temperature during high-intensity focused ultrasound hyperthermia based on the thermal wave model of bioheat transfer, Adv. Condens. Matter Phys., 2018 (2018), 1-10. |
[9] | C. W. Huang, M. K. Sun, B. T. Chen, J. Shieh, C. S. Chen, W. S. Chen, Simulation of thermal ablation by high-intensity focused ultrasound with temperature-dependent properties, Ultrason. Sonochem., 27 (2015), 456-465. |
[10] | T. Kamakura, T. Ishiwata, K. Matsuda, A new theoretical approach to the analysis of nonlinear sound beams using the oblate spheroidal coordinate system, J. Acoust. Soc. Am., 105 (1999), 3083-3086. |
[11] | S. H. Chang, R. Cao, Y. B. Zhang, P. G. Wang, S. J. Wu, Y. H. Qian, et al., Treatable focal region modulated by double excitation signal superimposition to realize platform temperature distribution during transcranial brain tumor therapy with high-intensity focused ultrasound, Chinese Phys. B, 27 (2018), 078701. |
[12] | S. A. Sapareto, W. C. Dewey, Thermal dose determination in cancer therapy, Int. J. Radiat. Oncol., 10 (1984), 787-800. |
[13] | M. J. Choi, S. R. Guntur, J. M. Lee, D. G. Paeng, K. I. Lee, A. Coleman, Changes in ultrasonic properties of liver tissue in vitro during heating-cooling cycle concomitant with thermal coagulation, Ultrasound Med. Biol., 37 (2011), 2000-2012. |
[14] | S. R. Guntur, K. I. Lee, D. G. Paeng, A. J. Coleman, M. J. Choi, Temperature-dependent thermal properties of ex vivo liver undergoing thermal ablation, Ultrasound Med. Biol., 39 (2013), 1771-1784. |
[15] | A. Ciubara, D. Dorohoi, F. Severcan, D. Creanga, Quantitative model of ultrasound propagation in biological media, Univ. Politeh. Bucharest Sci. Bull., Series A, 76 (2014), 221-226. |
[16] | S. R. Guntur, M. J. Choi, Influence of temperature-dependent thermal parameters on temperature elevation of tissue exposed to high-intensity focused ultrasound: numerical simulation, Ultrasound Med. Biol., 41 (2015), 806-813. |
[17] | A. Bhattacharya, R. L. Mahajan, Temperature dependence of thermal conductivity of biological tissues, Physiol. Meas., 24 (2003), 769-783. |
[18] | C. W. Connor, K. Hynynen, Bio-acoustic thermal lensing and nonlinear propagation in focused ultrasound surgery using large focal spots: a parametric study, Phys. Med. Biol., 47 (2002), 1911-1928. |
[19] | X. Liu, X. Gong, C. Yin, J. Li, D. Zhang, Noninvasive estimation of temperature elevations in biological tissues using acoustic nonlinearity parameter imaging, Ultrasound Med. Biol., 34 (2008), 414-424. |
[20] | R. T. Beyer, Parameter of nonlinearity in fluids, J. Acoust. Soc. Am., 32 (1960), 719-721. |
[21] | J. N Tjo/tta, S. Tjo/tta, E. H. Vefring, Effects of focusing on the nonlinear interaction between two collinear finite amplitude sound beams, J. Acoust. Soc. Am., 89 (1991), 1017-1027. |
[22] | I. M. Hallaj, R. O. Cleveland, FDTD simulation of finite-amplitude pressure and temperature fields for biomedical ultrasound, J. Acoust. Soc. Am., 105 (1999), L7-L12. |
[23] | T. Kamakura, T. Ishiwata, K. Matsuda, Model equation for strongly focused finite-amplitude sound beams, J. Acoust. Soc. Am., 107 (2000), 3035-3046. |
[24] | D. Yang, Z. Ni, Y. Yang, G. Xu, J. Tu, X. Guo, P. Huang, D. Zhang, The enhanced HIFU-induced thermal effect via magnetic ultrasound contrast agent microbubbles, Ultrason. Sonochem., 49 (2018), 111-117. |
[25] | E. J. Jackson, C. C. Coussios, R. O. Cleveland, Nonlinear acoustic properties of ex vivo bovine liver and the effects of temperature and denaturation, Phys. Med. Biol., 59 (2014), 3223-3238. |
[26] | J. Zhang, J. Hills, Y. Zhong, B. Shirinzadeh, J. Smith, G. Gu, Temperature-dependent thermomechanical modeling of soft tissue deformation, J. Mech. Med. Biol., 18 (2019), 1840021. |
[27] | J. Zhang, R. J. Lay, S. K. Roberts, S. Chauhan, Towards real-time finite-strain anisotropic thermo-visco-elastodynamic analysis of soft tissues for thermal ablative therapy, Comput. Methods Programs Biomed., 198 (2020), 105789. |
[28] | J. Zhang, S. Chauhan, Fast computation of soft tissue thermal response under deformation based on fast explicit dynamics finite element algorithm for surgical simulation, Comput. Methods Programs Biomed., 187 (2020), 105244. |