Control of muscle relaxation during anesthesia: a novel approach for clinical routine

During general anesthesia drugs are administered to provide hypnosis, ensure analgesia, and skeletal muscle relaxation. In this paper, the main components of a newly developed controller for skeletal muscle relaxation are described. Muscle relaxation is controlled by administration of neuromuscular blocking agents. The degree of relaxation is assessed by supramaximal train-of-four stimulation of the ulnar nerve and measuring the electromyogram response of the adductor pollicis muscle. For closed-loop control purposes, a physiologically based pharmacokinetic and pharmacodynamic model of the neuromuscular blocking agent mivacurium is derived. The model is used to design an observer-based state feedback controller. Contrary to similar automatic systems described in the literature this controller makes use of two different measures obtained in the train-of-four measurement to maintain the desired level of relaxation. The controller is validated in a clinical study comparing the performance of the controller to the performance of the anesthesiologist. As presented, the controller was able to maintain a preselected degree of muscle relaxation with excellent precision while minimizing drug administration. The controller performed at least equally well as the anesthesiologist.

Clinical validation of electromyography and acceleromyography as sensors for muscle relaxation

BACKGROUND AND OBJECTIVE: The aim of this study was to determine which of two clinically applied methods, electromyography or acceleromyography, was less affected by external disturbances, had a higher sensitivity and which would provide the better input signal for closed loop control of muscle relaxation. METHODS: In 14 adult patients, anaesthesia was induced with intravenous opioids and propofol. The response of the thumb to ulnar nerve stimulation was recorded on the same arm. Mivacurium was used for neuromuscular blockade. Under stable conditions of relaxation, the infusion-rate was decreased and the effects of turning the hand were investigated. RESULTS: Electromyography and acceleromyography both reflected the change of the infusion rate (P = 0.015 and P < 0.001, respectively). Electromyography was significantly less affected by the hand-turn (P = 0.008) than acceleromyography. While zero counts were detected with acceleromyography, electromyography could still detect at least one count in 51.1%. CONCLUSIONS: Electromyography is more reliable for use in daily practice as it is less influenced by external disturbances than acceleromyography.

Two dimensional Lattice Boltzmann simulation of fluid flow through an idealized micro-structure of disordered media

This technical report discusses the application of Lattice Boltzmann Method (LBM) in the fluid flow simulation through porous filter-wall of disordered media. The diesel particulate filter (DPF) is an example of disordered media. DPF is developed as a cutting edge technology to reduce harmful particulate matter in the engine exhaust. Porous filter-wall of DPF traps these soot particles in the after-treatment of the exhaust gas. To examine the phenomena inside the DPF, researchers are looking forward to use the Lattice Boltzmann Method as a promising alternative simulation tool. The lattice Boltzmann method is comparatively a newer numerical scheme and can be used to simulate fluid flow for single-component single-phase, single-component multi-phase. It is also an excellent method for modelling flow through disordered media. The current work focuses on a single-phase fluid flow simulation inside the porous micro-structure using LBM. Firstly, the theory concerning the development of LBM is discussed. LBM evolution is always related to Lattice gas Cellular Automata (LGCA), but it is also shown that this method is a special discretized form of the continuous Boltzmann equation. Since all the simulations are conducted in two-dimensions, the equations developed are in reference with D2Q9 (two-dimensional 9-velocity) model. The artificially created porous micro-structure is used in this study. The flow simulations are conducted by considering air and CO2 gas as fluids. The numerical model used in this study is explained with a flowchart and the coding steps. The numerical code is constructed in MATLAB. Different types of boundary conditions and their importance is discussed separately. Also the equations specific to boundary conditions are derived. The pressure and velocity contours over the porous domain are studied and recorded. The results are compared with the published work. The permeability values obtained in this study can be fitted to the relation proposed by Nabovati [8], and the results are in excellent agreement within porosity range of 0.4 to 0.8.

LATTICE BOLTZMANN METHOD AND CELLULAR AUTOMATA SIMULATION OF PARTICLE MOTION AND DEPOSITION IN 2-D CASE

This technical report discusses the application of the Lattice Boltzmann Method (LBM) and Cellular Automata (CA) simulation in fluid flow and particle deposition. The current work focuses on incompressible flow simulation passing cylinders, in which we incorporate the LBM D2Q9 and CA techniques to simulate the fluid flow and particle loading respectively. For the LBM part, the theories of boundary conditions are studied and verified using the Poiseuille flow test. For the CA part, several models regarding simulation of particles are explained. And a new Digital Differential Analyzer (DDA) algorithm is introduced to simulate particle motion in the Boolean model. The numerical results are compared with a previous probability velocity model by Masselot [Masselot 2000], which shows a satisfactory result.

T1(Gd) gives comparable information as Delta T1 relaxation rate in dGEMRIC evaluation of cartilage repair tissue

OBJECTIVES: To evaluate the relationship between T1 after intravenous contrast administration (T1Gd) and Delta relaxation rate (DeltaR1) = (1/T1(Gd) - 1/T1o) in the delayed Gadolinium-Enhanced MRI of cartilage (dGEMRIC) evaluation of cartilage repair tissue. MATERIALS AND METHODS: Thirty single MR examinations from 30 patients after matrix-associated autologous chondrocyte transplantations of the knee joint with different postoperative intervals were examined using an 8-channel knee-coil at 3T. T1 mapping using a 3D GRE sequence with a 35/10 degrees flip angle excitation pulse combination was performed before and after contrast administration (dGEMRIC technique). T1 postcontrast (T1(Gd)) and the DeltaR1 (relative index of pre- and postcontrast R1 value) were calculated for repair tissue and the weight-bearing normal appearing control cartilage. For evaluation of the different postoperative intervals, MR exams were subdivided into 3 groups (up to 12 months, 12-24 months, more than 24 months). For statistical analysis Spearman correlation coefficients were calculated. RESULTS: The mean value for T1 postcontrast was 427 +/- 159 ms, for DeltaR1 1.85 +/- 1.0; in reference cartilage 636 +/- 181 ms for T1 postcontrast and 0.83 +/- 0.5 for DeltaR1.The correlation coefficients were highly significant between T1 (Gd) and DeltaR1 for repair tissue (0.969) as well as normal reference cartilage (0.928) in total, and for the reparative cartilage in the early, middle postoperative, and late postoperative interval after surgery (R values: -0.986, -0.970, and -0.978, respectively). Using either T1(Gd) or DeltaR1, the 2 metrics resulted in similar conclusions regarding the time course of change of repair tissue and control tissue, namely that highly significant (P > 0.01) differences between cartilage repair tissue and reference cartilage were found for all follow-up groups. Additionally, for both metrics highly significant differences (P < 0.01) between early follow up and the 2 later postoperative groups for cartilage repair tissue were found. No statistical differences were found between the 2 later follow-up groups of reparative cartilage either for T1 (Gd) or DeltaR1. CONCLUSION: The high correlation between T1 (Gd) and DeltaR1 and the comparable conclusions reached utilizing metric implies that T1 mapping before intravenous administration of MR contrast agent is not necessary for the evaluation of repair tissue. This will help to reduce costs, inconvenience for the patients, simplifies the examination procedure, and makes dGEMRIC more attractive for follow-up of patients after cartilage repair surgeries.

Integrated data acquisition and processing to determine metabolite contents, relaxation times, and macromolecule baseline in single examinations of individual subjects

Absolute quantitation of clinical (1)H-MR spectra is virtually always incomplete for single subjects because the separate determination of spectrum, baseline, and transverse and longitudinal relaxation times in single subjects is prohibitively long. Integrated Processing and Acquisition of Data (IPAD) based on a combined 2-dimensional experimental and fitting strategy is suggested to substantially improve the information content from a given measurement time. A series of localized saturation-recovery spectra was recorded and combined with 2-dimensional prior-knowledge fitting to simultaneously determine metabolite T(1) (from analysis of the saturation-recovery time course), metabolite T(2) (from lineshape analysis based on metabolite and water peak shapes), macromolecular baseline (based on T(1) differences and analysis of the saturation-recovery time course), and metabolite concentrations (using prior knowledge fitting and conventional procedures of absolute standardization). The procedure was tested on metabolite solutions and applied in 25 subjects (15-78 years old). Metabolite content was comparable to previously found values. Interindividual variation was larger than intraindividual variation in repeated spectra for metabolite content as well as for some relaxation times. Relaxation times were different for various metabolite groups. Parts of the interindividual variation could be explained by significant age dependence of relaxation times.

High-field 1H T1 and T2 NMR Relaxation Time Measurements of H2O in Homeopathic Preparations of Quartz, Sulfur and Copper Sulfate

Quantitative meta-analyses of randomized clinical trials investigating the specific therapeutic efficacy of homeopathic remedies yielded statistically significant differences compared to placebo. Since the remedies used contained mostly only very low concentrations of pharmacologically active compounds, these effects cannot be accounted for within the framework of current pharmacology. Theories to explain clinical effects of homeopathic remedies are partially based upon changes in diluent structure. To investigate the latter, we measured for the first time high-field (600/500 MHz) 1H T1 and T2 nuclear magnetic resonance relaxation times of H2O in homeopathic preparations with concurrent contamination control by inductively coupled plasma mass spectrometry (ICP-MS). Homeopathic preparations of quartz (10c–30c, n = 21, corresponding to iterative dilutions of 100−10–100−30), sulfur (13x–30x, n = 18, 10−13–10−30), and copper sulfate (11c–30c, n = 20, 100−11–100−30) were compared to n = 10 independent controls each (analogously agitated dilution medium) in randomized and blinded experiments. In none of the samples, the concentration of any element analyzed by ICP-MS exceeded 10 ppb. In the first measurement series (600 MHz), there was a significant increase in T1 for all samples as a function of time, and there were no significant differences between homeopathic potencies and controls. In the second measurement series (500 MHz) 1 year after preparation, we observed statistically significant increased T1 relaxation times for homeopathic sulfur preparations compared to controls. Fifteen out of 18 correlations between sample triplicates were higher for controls than for homeopathic preparations. No conclusive explanation for these phenomena can be given at present. Possible hypotheses involve differential leaching from the measurement vessel walls or a change in water molecule dynamics, i.e., in rotational correlation time and/or diffusion. Homeopathic preparations thus may exhibit specific physicochemical properties that need to be determined in detail in future investigations.