Model-order reduction of nonlinear models of electromagnetic phased-array hyperthermia

A method based on the Karhunen-Loeve (KL) transform is proposed for the reduction of large-scale, nonlinear ordinary differential equations such as those arising from the finite difference modeling of biological heat transfer. The method of snapshots is used to expedite computation of the required quantities in the KL procedure. Guidelines are presented and validated for snapshot selection and resultant basis series truncation, emphasizing the special physical features of the electromagnetic phased-array heat transfer physics. Applications to fast temperature prediction are presented.

[1]  P B Hoffer,et al.  Computerized three-dimensional segmented human anatomy. , 1994, Medical physics.

[2]  W. Dewey,et al.  Thermal dose determination in cancer therapy. , 1984, International journal of radiation oncology, biology, physics.

[3]  P Wust,et al.  ESHO quality assurance guidelines for regional hyperthermia. , 1998, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[4]  Michel Loève,et al.  Probability Theory I , 1977 .

[5]  F. Bardati,et al.  SAR optimization in a phased array radiofrequency hyperthermia system , 1995, IEEE Transactions on Biomedical Engineering.

[6]  H. H. Pennes.  Analysis of tissue and arterial blood temperatures in the resting human forearm. 1948. , 1948, Journal of applied physiology.

[7]  R B Roemer,et al.  Engineering aspects of hyperthermia therapy. , 1999, Annual review of biomedical engineering.

[8]  M. Seebass,et al.  Impact of nonlinear heat transfer on temperature control in regional hyperthermia , 1999, IEEE Transactions on Biomedical Engineering.

[9]  J. Overgaard,et al.  Randomised trial of hyperthermia as adjuvant to radiotherapy for recurrent or metastatic malignant melanoma , 1995, The Lancet.

[10]  R B Roemer,et al.  Optimal actuator placement for large scale systems: a reduced-order modelling approach. , 1998, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[11]  J.K. Potocki,et al.  Reduced-order modeling for hyperthermia control , 1992, IEEE Transactions on Biomedical Engineering.

[12]  Om P. Gandhi,et al.  Numerical simulation of annular phased arrays for anatomically based models using the FDTD method , 1989 .

[13]  Jian-Ming Jin,et al.  Determination of electromagnetic phased-array driving signals for hyperthermia based on a steady-state temperature criterion , 2000 .

[14]  A W Dutton,et al.  A comparison of reduced-order modelling techniques for application in hyperthermia control and estimation. , 1998, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[15]  C Gabriel,et al.  The dielectric properties of biological tissues: I. Literature survey. , 1996, Physics in medicine and biology.

[16]  Paul F. Turner,et al.  Regional Hyperthermia with an Annular Phased Array , 1984, IEEE Transactions on Biomedical Engineering.

[17]  Dennis M. Sullivan,et al.  Mathematical methods for treatment planning in deep regional hyperthermia , 1991 .

[18]  Carey M. Rappaport,et al.  FDTD verification of deep-set brain tumor hyperthermia using a spherical microwave source distribution , 1996 .

[19]  K D Paulsen,et al.  Finite element simulation of Sigma 60 heating in the Utah phantom: computed and measured data compared. , 1994, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[20]  J. Marsden,et al.  A subspace approach to balanced truncation for model reduction of nonlinear control systems , 2002 .

[21]  Ricardo C. H. del Rosario,et al.  Reduced-order model feedback control design: numerical implementation in a thin shell model , 2000, IEEE Trans. Autom. Control..

[22]  Dennis M. Sullivan,et al.  Three-dimensional computer simulation in deep regional hyperthermia using the finite-difference time-domain method , 1990 .

[23]  M. Fahl,et al.  Computation of POD basis functions for fluid flows with lanczos methods , 2001 .

[24]  Dennis M. Sullivan,et al.  Comparison of measured and simulated data in an annular phased array using an inhomogeneous phantom , 1992 .

[25]  S K Das,et al.  Computational techniques for fast hyperthermia temperature optimization. , 1999, Medical physics.

[26]  P Wust,et al.  Simulation studies promote technological development of radiofrequency phased array hyperthermia. , 1996, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[27]  Peter Deuflhard,et al.  A high-resolution interpolation at arbitrary interfaces for the FDTD method , 1998 .

[28]  Panagiotis D. Christofides,et al.  Output feedback control of parabolic PDE systems with nonlinear spatial differential operators , 1999 .

[29]  Magdy F. Iskander,et al.  FDTD analysis of power deposition patterns of an array of interstitial antennas for use in microwave hyperthermia , 1992 .

[30]  O. Gandhi,et al.  Numerical simulation of annular-phased arrays of dipoles for hyperthermia of deep-seated tumors , 1992, IEEE Transactions on Biomedical Engineering.

[31]  Joel R. Phillips.  Automated extraction of nonlinear circuit macromodels , 2000, Proceedings of the IEEE 2000 Custom Integrated Circuits Conference (Cat. No.00CH37044).

[32]  R. W. Lau,et al.  The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues. , 1996, Physics in medicine and biology.

[33]  K. Paulsen,et al.  Evaluation of perfectly matched layer mesh terminations in finite-element bioelectromagnetic scattering computations , 1999 .

[34]  R Vanderby,et al.  Temperature-dependent versus constant-rate blood perfusion modelling in ferromagnetic thermoseed hyperthermia: results with a model of the human prostate. , 1994, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[35]  Francis A. Duck,et al.  Physical properties of tissue : a comprehensive reference book , 1990 .

[36]  Kenneth R. Diller,et al.  Heat Transfer in Living Systems: Current Opportunities , 1998 .

[37]  I. Jolliffe.  Principal Component Analysis , 2002 .

[38]  J K Potocki,et al.  Concurrent hyperthermia estimation schemes based on extended Kalman filtering and reduced-order modelling. , 1993, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[39]  P. Deuflhard,et al.  Strategies for optimized application of annular-phased-array systems in clinical hyperthermia. , 1991, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[40]  Andrew G. Webb,et al.  Optimization of electromagnetic phased-arrays for hyperthermia via magnetic resonance temperature estimation , 2002, IEEE Transactions on Biomedical Engineering.

[41]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[42]  H. Tran,et al.  Modeling and control of physical processes using proper orthogonal decomposition , 2001 .

[43]  H. Arkin,et al.  Recent developments in modeling heat transfer in blood perfused tissues , 1994, IEEE Transactions on Biomedical Engineering.

[44]  K. Paulsen,et al.  Finite element computations of specific absorption rates in anatomically conforming full-body models for hyperthermia treatment analysis , 1993, IEEE Transactions on Biomedical Engineering.

[45]  M. Mattingly,et al.  Reduced-order modeling for hyperthermia: an extended balanced-realization-based approach , 1998, IEEE Transactions on Biomedical Engineering.

[46]  W. Wulff,et al.  The Energy Conservation Equation for Living Tissue , 1974 .

[47]  N. Maratos,et al.  Optimization of the deposited power distribution inside a layered lossy medium irradiated by a coupled system of concentrically placed waveguide applicators , 1998, IEEE Transactions on Biomedical Engineering.

[48]  Michel Nakhla,et al.  Nonlinear circuit-reduction of high-speed interconnect networks using congruent transformation techniques , 2001 .

[49]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .