Muscle Swelling and Neuromuscular Responses Following Blood Flow Restricted Exercise in Untrained Women.

Purpose: There is conflicting evidence related to the prevalence and magnitude of exercise-induced muscle damage (EIMD) following four sets to volitional failure with BFR (BFR-F) or 75 total repetitions with BFR (1 × 30, 3 × 15, BFR-75). The purpose of this investigation was to examine muscle swelling, peak torque, and neuromuscular responses following BFR-75 and BFR-F. Methods: Thirteen untrained women completed unilateral isokinetic (120°s-1) leg extensions concentric-eccentric at 30% of their maximal voluntary isometric contraction (MVIC) using BFR-75 and BFR-F protocols, separated by 15 minutes. Ultrasound was used to assess muscle thickness, cross sectional area, and echo intensity of the rectus femoris and vastus lateralis before, 0-, 24-, 48-, 72-, and 96-hours post-exercise. Peak torque and surface electromyography (sEMG) were recorded during MVICs before, 24-, 48-, 72-, and 96-hours post-exercise to determine sEMG amplitude, frequency, and neuromuscular efficiency. Results: There were no differences between conditions. Collapsed across conditions, muscle thickness and cross-sectional area increased at 0-hours for the rectus femoris (2.5 ± 0.4, 2.8 ± 0.4 cm, 10.6 ± 1.8, 12.1 ± 1.8 cm2, respectively) and vastus lateralis (2.1 ± 0.5, 2.5 ± 0.7 cm; 22.2 ± 3.9, 25.1 ± 4.5 cm2, respectively), but returned to baseline at 24-hours. There were no changes in echo intensity, sEMG amplitude, sEMG frequency, or neuromuscular efficiency. MVIC peak torque increased relative to pre-exercise at 24-, 48-, 72-, and 96-hours (159.9 ± 34.9, 171.4 ± 30.1-179.1 ± 35.6 Nm). Conclusion: These results suggest that BFR-75 and BFR-F did not cause EIMD but caused an acute increase in muscle swelling that returned to baseline 24-hours post-exercise.

SCG variability and spectral energy distribution during normal breathing and breath hold at different lung volumes and airway pressures.

Seismocardiographic (SCG) signals are chest wall vibrations induced by cardiac activity and are potentially useful for cardiac monitoring and diagnosis. SCG waveform is observed to vary with respiration, but the mechanism of these changes is poorly understood as alterations in autonomic tone, lung volume, heart location and intrathoracic pressure are all varying during the respiratory cycle. Understanding SCG variability and its sources may help reduce variability and increase SCG clinical utility. This study investigated SCG variability during breath holding (BH) at two different lung volumes (i.e., end inspiration and end expiration) and five airway pressures (i.e., 0, ± 2-4, and ± 15-20 cm H2O). Variability during normal breathing was also studied with and without grouping SCG beats into two clusters of similar waveform morphologies (performed using the K-medoid algorithm in an unsupervised machine learning fashion). The study included 15 healthy subjects (11 Females and 4 males, Age: 21 ± 2 y) where SCG, ECG, and spirometry were simultaneously acquired. SCG waveform variability was calculated at each experimental state (i.e., lung volume and airway pressure). Results showed that breath holding was more effective in reducing the intra-state variability of SCG than clustering normal breathing data. For the BH states, the intra-state variability increased as the airway pressure deviated from zero. The subaudible-to-audible energy ratio of the BH states increased as the airway pressure decreased below zero which may be related to the effect of the intrathoracic pressure on cardiac afterload and blood ejection. When combining the BH waveforms at end inspiration and end expiration states (at the same airway pressures) into one group, the intra-state variability increased, which suggests that the lung volume and associated change in heart location were a significant source of variability. The linear trend between airway pressure and waveform changes was found to be statistically significant for BH at end expiration. To confirm these findings, more studies are needed with a larger number of airway pressure levels and larger number of subjects.

Blood flow restriction attenuates surface mechanomyography lateral and longitudinal, but not transverse oscillations during fatiguing exercise.

Objective. Surface mechanomyography (sMMG) can measure oscillations of the activated muscle fibers in three axes (i.e.X,Y, andZ-axes) and has been used to describe motor unit activation patterns (X-axis). The application of blood flow restriction (BFR) is common in exercise studies, but the cuff may restrict muscle fiber oscillations. Therefore, the purpose of this investigation was to examine the acute effects of submaximal, fatiguing exercise with and without BFR on sMMG amplitude in theX,Y, andZ-axes among female participants.Approach. Sixteen females (21 ± 1 years) performed two separate exercise bouts to volitional exhaustion that consisted of unilateral, submaximal (50% maximal voluntary isometric contraction [MVIC]) intermittent, isometric, leg extensions with and without BFR. sMMG was recorded and examined across percent time to exhaustion (%TTE) in 20% increments. Separate 2-way repeated measures ANOVA models were constructed: (condition [BFR, non-BFR]) × (time [20, 40, 60, 80, and 100% TTE]) to examine absolute (m·s-2) and normalized (% of pretest MVIC) sMMG amplitude in theX-(sMMG-X),Y-(sMMG-Y), andZ-(sMMG-Z) axes.Main results. The absolute sMMG-X amplitude responses were attenuated with the application of BFR (mean ± SD = 0.236 ± 0.138 m·s-2) relative to non-BFR (0.366 ± 0.199 m·s-2, collapsed across time) and for sMMG-Y amplitude at 60%-100% of TTE (BFR range = 0.213-0.232 m·s-2versus non-BFR = 0.313-0.445 m·s-2). Normalizing sMMG to pretest MVIC removed most, but not all the attenuation which was still evident for sMMG-Y amplitude at 100% of TTE between BFR (72.9 ± 47.2%) and non-BFR (98.9 ± 53.1%). Interestingly, sMMG-Z amplitude was not affected by the application of BFR and progressively decreased across %TTE (0.332 ± 0.167 m·s-2to 0.219 ± 0.104 m·s-2, collapsed across condition.)Significance. The application of BFR attenuated sMMG-X and sMMG-Y amplitude, although normalizing sMMG removed most of this attenuation. Unlike theXandY-axes, sMMG-Z amplitude was not affected by BFR and progressively decreased across each exercise bout potentially tracking the development of muscle fatigue.

Postural and longitudinal variability in seismocardiographic signals.

Objective. Low frequency cardiovascular vibrations detectable on the chest surface (termed seismocardiography or SCG) may be useful for non-invasive diagnosis and monitoring of various cardiovascular conditions. A potential limitation of using SCG for longitudinal patient monitoring is the existence of intra-subject variability, which can contribute to errors in calculating SCG features. Improved understanding of the contribution of intra-subject variability sources may lead to improved SCG utility. This study aims to quantify postural and longitudinal SCG variability in healthy resting subjects during normal breathing.Approach. SCG and ECG signals were longitudinally acquired in 19 healthy subjects at different postures (supine, 45° head up, and sitting) during five recording sessions over five months. SCG cycles were segmented using the ECG R wave. Unsupervised machine learning was used to reduce SCG variability due to respiration by grouping the SCG signals into two clusters with minimized intra-cluster waveform heterogeneity. Several SCG features were assessed at different postures and longitudinally.Main results. SCG waveform morphological variability was calculated within each cluster (intra-cluster) and between two clusters (inter-cluster) at each posture and data collection session. The variabilities were significantly different between the supine and sitting but not between supine and 45° postures. For the 45° and sitting postures, the intra-cluster variability was not significantly different, while the inter-cluster variability difference was significant. The energy ratio between different frequency bands to total spectral energy in 0.5-50 Hz were calculated and were comparable for all postures. The combined cardiac timing intervals from the two clusters showed significant variation with postural changes. There was significant heart rate difference between the clusters and between postural positions. The SCG features were compared between longitudinal sessions and all features were not significantly different,Significance. Several SCG features significantly varied with posture suggesting that posture needs to be specified when comparing SCG changes over time. Longitudinally comparable SCG feature values suggests that significant longitudinal differences, if observed, may reflect true alternations in the cardiac functioning over time.

Noninvasive detection of elevated ICP using spontaneous tympanic membrane pulsation.

Neurological conditions such as traumatic brain injury (TBI) and hydrocephalus may lead to intracranial pressure (ICP) elevation. Current diagnosis methods rely on direct pressure measurement, while CT, MRI and other expensive imaging may be used. However, these invasive or expensive testing methods are often delayed because symptoms of elevated ICP are non-specific. Invasive methods, such as intraventricular catheter, subdural screw, epidural sensor, lumbar puncture, are associated with an increased risk of infection and hemorrhage. On the other hand, noninvasive, low-cost, accurate methods of ICP monitoring can help avoid risks and reduce costs while expediting diagnosis and treatment. The current study proposes and evaluates a novel method for noninvasive ICP monitoring using tympanic membrane pulsation (TMp). These signals are believed to be transmitted from ICP to the auditory system through the cochlear aqueduct. Fifteen healthy subjects were recruited and TMp signals were acquired noninvasively while the subjects performed maneuvers that are known to change ICP. A custom made system utilizing a stethoscope headset and a pressure transducer was used to perform these measurements. Maneuvers included head-up-tilt, head-down-tilt and hyperventilation. When elevated ICP was induced, significant TMp waveform morphological changes were observed in each subject (p < 0.01). These changes include certain waveform slopes and high frequency wave features. The observed changes were reversed by the maneuvers that decreased ICP (p < .01). The study results suggest that TMp waveform measurement and analysis may offer an inexpensive, noninvasive, accurate tool for detection and monitoring of ICP elevations. Further studies are warranted to validate this technique in patients with pathologically elevated ICP.

The "Brain Stethoscope": A Non-Invasive Method for Detecting Elevated Intracranial Pressure.

Introduction Minimally invasive intracranial pressure (ICP) screening has long been desired by neurosurgeons. A novel approach deriving ICP from tympanic membrane (TM) pulsation may offer the solution. The ICP waveform appears to be transmitted to the TM by the cochlear aqueduct. The resulting TM infrasonic pulsations can be measured by certain sensors. Elevated ICP alters brain compliance, which appears to yield slower rise times of the TM pulsation waveform. Measurement of this change may be useful in screening for elevated ICP. This paper investigates one such technique. Methods A stethoscope was modified for airtight external ear canal fit; the dome was exchanged for a magnetic reluctance pressure sensor, allowing measurement of TM pulsations. Analog TM pulsations were analyzed by measuring the pulsation's slope ratio between the waveform's downslope and upslope. Seventeen normal subjects (ages 18-32 years) underwent hyperventilation and tilt table testing to induce ICP changes. An algorithm processed this data and predicted the subject's ICP status. Results The slope ratio method showed consistent and stable changes with the expected alterations in ICP from the tilt test and hyperventilation maneuvers. The classification algorithm correctly identified subjects with elevated ICP in 60 of 60 independent recordings on 17 subjects. Conclusion This paper has four conclusions. First, the "brain stethoscope" can detect increased ICP from the TM pulsation waveform in healthy subjects. Second, analysis of the TM waveform using slope ratio calculations is capable of distinguishing normal versus elevated ICP. Third, the tilt and hyperventilation maneuvers showed the expected physiologic trends. Last, further studies are needed on patients with pathological ICP before the brain stethoscope can be implemented into clinical practice.

Sound transmission in human thorax through airway insonification: an experimental and computational study with diagnostic applications.

Pulmonary diseases and injury lead to structural and functional changes in the lung parenchyma and airways, often resulting in measurable sound transmission changes on the chest wall surface. Additionally, noninvasive imaging of externally driven mechanical wave motion in the chest (e.g., using magnetic resonance elastography) can provide information about lung stiffness and other structural property changes which may be of diagnostic value. In the present study, a comprehensive computational simulation (in silico) model was developed to simulate sound wave propagation in the airways, parenchyma, and chest wall under normal and pathological conditions that create distributed structural (e.g., pneumothoraces) and diffuse material (e.g., fibrosis) changes, as well as a localized structural and material changes as may be seen with a neoplasm. Experiments were carried out in normal subjects to validate the baseline model. Sound waves with frequency content from 50 to 600 Hz were introduced into the airways of three healthy human subjects through the mouth, and transthoracic transmitted waves were measured by scanning laser Doppler vibrometry at the chest wall surface. The computational model predictions of a frequency-dependent decreased sound transmission due to pneumothorax were consistent with experimental measurements reported in previous work. Predictions for the case of fibrosis show that while shear wave motion is altered, changes to compression wave propagation are negligible, and thus, insonification, which primarily drives compression waves, is not ideal to detect the presence of fibrosis. Results from the numerical simulation of a tumor show an increase in the wavelength of propagating waves in the immediate vicinity of the tumor region. Graphical abstract.

Extraction of Peak Velocity Profiles from Doppler Echocardiography Using Image Processing.

The objective of this study is to extract positive and negative peak velocity profiles from Doppler echocardiographic images. These profiles are currently estimated using tedious manual approaches. Profiles can be used to establish realistic boundary conditions for computational hemodynamic studies and to estimate cardiac time intervals, which are of clinical utility. In the current study, digital image processing algorithms that rely on intensity calculations and two different thresholding methods were proposed and tested. Image intensity histograms were used to guide threshold choices, which were selected such that the resulting velocity profiles appropriately represent Doppler shift envelopes. The resulting peak velocity profiles contained artifacts in the form of sudden velocity changes and possible outliers. To reduce these artifacts, image smoothing using a moving average process was then implemented. Bland-Altman analysis suggested good agreement between the two thresholding methods. Artifacts decreased when image smoothing was performed. Results also suggested that one thresholding method tended to provide the lower limit (i.e., underestimate) of velocities, while the second tended to provide the velocity upper limit (i.e., overestimate). Combining estimates from both methods appeared to provide a smoother peak velocity profile estimate. The proposed automated approach may be useful for objective estimation of peak velocity profiles, which may be helpful for computational hemodynamic studies and may increase the efficiency of current clinical diagnostic tools.

Recent Advances in Seismocardiography.

Cardiovascular disease is a major cause of death worldwide. New diagnostic tools are needed to provide early detection and intervention to reduce mortality and increase both the duration and quality of life for patients with heart disease. Seismocardiography (SCG) is a technique for noninvasive evaluation of cardiac activity. However, the complexity of SCG signals introduced challenges in SCG studies. Renewed interest in investigating the utility of SCG accelerated in recent years and benefited from new advances in low-cost lightweight sensors, and signal processing and machine learning methods. Recent studies demonstrated the potential clinical utility of SCG signals for the detection and monitoring of certain cardiovascular conditions. While some studies focused on investigating the genesis of SCG signals and their clinical applications, others focused on developing proper signal processing algorithms for noise reduction, and SCG signal feature extraction and classification. This paper reviews the recent advances in the field of SCG.

Adverse Hemodynamic Conditions Associated with Mechanical Heart Valve Leaflet Immobility.

Artificial heart valves may dysfunction, leading to thrombus and/or pannus formations. Computational fluid dynamics is a promising tool for improved understanding of heart valve hemodynamics that quantify detailed flow velocities and turbulent stresses to complement Doppler measurements. This combined information can assist in choosing optimal prosthesis for individual patients, aiding in the development of improved valve designs, and illuminating subtle changes to help guide more timely early intervention of valve dysfunction. In this computational study, flow characteristics around a bileaflet mechanical heart valve were investigated. The study focused on the hemodynamic effects of leaflet immobility, specifically, where one leaflet does not fully open. Results showed that leaflet immobility increased the principal turbulent stresses (up to 400%), and increased forces and moments on both leaflets (up to 600% and 4000%, respectively). These unfavorable conditions elevate the risk of blood cell damage and platelet activation, which are known to cascade to more severe leaflet dysfunction. Leaflet immobility appeared to cause maximal velocity within the lateral orifices. This points to the possible importance of measuring maximal velocity at the lateral orifices by Doppler ultrasound (in addition to the central orifice, which is current practice) to determine accurate pressure gradients as markers of valve dysfunction.

Modeling Inspiratory Flow in a Porcine Lung Airway.

Inspiratory flow in a multigeneration pig lung airways was numerically studied at a steady inlet flow rate of 3.2 × 10-4 m3/s corresponding to a Reynolds number of 1150 in the trachea. The model was validated by comparing velocity distributions with previous measurements and simulations in simplified airway geometries. Simulation results provided detailed maps of the axial and secondary flow patterns at different cross sections of the airway tree. The vortex core regions in the airways were visualized using absolute helicity values and suggested the presence of secondary flow vortices where two counter-rotating vortices were observed at the main bifurcation and in many other bifurcations. Both laminar and turbulent flows were considered. Results showed that axial and secondary flows were comparable in the laminar and turbulent cases. Turbulent kinetic energy (TKE) vanished in the more distal airways, which indicates that the flow in these airways approaches laminar flow conditions. The simulation results suggested viscous pressure drop values comparable to earlier studies. The monopodial asymmetric nature of airway branching in pigs resulted in airflow patterns that are different from the less asymmetric human airways. The major daughters of the pig airways tended to have high airflow ratios, which may lead to different particle distribution and sound generation patterns. These differences need to be taken into consideration when interpreting the results of animal studies involving pigs before generalizing these results to humans.

Time-Frequency Distribution of Seismocardiographic Signals: A Comparative Study.

Accurate estimation of seismocardiographic (SCG) signal features can help successful signal characterization and classification in health and disease. This may lead to new methods for diagnosing and monitoring heart function. Time-frequency distributions (TFD) were often used to estimate the spectrotemporal signal features. In this study, the performance of different TFDs (e.g., short-time Fourier transform (STFT), polynomial chirplet transform (PCT), and continuous wavelet transform (CWT) with different mother functions) was assessed using simulated signals, and then utilized to analyze actual SCGs. The instantaneous frequency (IF) was determined from TFD and the error in estimating IF was calculated for simulated signals. Results suggested that the lowest IF error depended on the TFD and the test signal. STFT had lower error than CWT methods for most test signals. For a simulated SCG, Morlet CWT more accurately estimated IF than other CWTs, but Morlet did not provide noticeable advantages over STFT or PCT. PCT had the most consistently accurate IF estimations and appeared more suited for estimating IF of actual SCG signals. PCT analysis showed that actual SCGs from eight healthy subjects had multiple spectral peaks at 9.20 ± 0.48, 25.84 ± 0.77, 50.71 ± 1.83 Hz (mean ± SEM). These may prove useful features for SCG characterization and classification.

Geometric features of pig airways using computed tomography.

Accurate knowledge of the airway geometry is needed when constructing physical models of the airway tree and for numerical modeling of flow or sound propagation in the airways. Human and animal experiments are conducted to validate these models. Many studies documented the geometric details of the human airways. However, information about the geometry of pig airways is scarcer. Earlier studies suggested that the morphology of animal airways can be significantly different from that of humans. The objective of this study is to measure the airway diameter, length and bifurcation angles in domestic pigs using computed tomography. In this study, lungs of six pigs were imaged, then segmentation software tools were used to extract the geometry of the airway lumen. The airway dimensions were measured from the resulting 3-D models for the first 24 airway generations. Results showed that the size and morphology of the airways of the six pigs were similar. The trachea diameters were found to be comparable to the typical human adult, but the diameter, length and branching angles of other airways were noticeably different from that of humans. For example, pig airways consistently had an early branching from the trachea that feeds the top right lung lobe and precedes the main carina. This branch is absent in the human airways. The results suggested that the pig airways geometry may not be accurately approximated by human airways and this approximation may contribute to increasing the errors in computational models of the pig chest.

Generation of Pig Airways using Rules Developed from the Measurements of Physical Airways.

BACKGROUND: A method for generating bronchial tree would be helpful when constructing models of the tree for benchtop experiments as well as for numerical modeling of flow or sound propagation in the airways. Early studies documented the geometric details of the human airways that were used to develop methods for generating human airway tree. However, methods for generating animal airway tree are scarcer. Earlier studies suggested that the morphology of animal airways can be significantly different from that of humans. Hence, using algorithms for the human airways may not be accurate in generating models of animal airway geometry. OBJECTIVE: The objective of this study is to develop an algorithm for generating pig airway tree based on the geometric details extracted from the physical measurements. METHODS: In the current study, measured values of branch diameters, lengths and bifurcation angles and rotation of bifurcating planes were used to develop an algorithm that is capable of generating a realistic pig airway tree. RESULTS: The generation relations between parent and daughter branches were found to follow certain trends. The diameters and the length of different branches were dependent on airway generations while the bifurcation angles were primarily dependent on bifurcation plane rotations. These relations were sufficient to develop rules for generating a model of the pig large airways. CONCLUSION: The results suggested that the airway tree generated from the algorithm can provide an approximate geometric model of pig airways for computational and benchtop studies.

Sound transmission in porcine thorax through airway insonification.

Many pulmonary injuries and pathologies may lead to structural and functional changes in the lungs resulting in measurable sound transmission changes on the chest surface. Additionally, noninvasive imaging of externally driven mechanical wave motion in the chest (e.g., using magnetic resonance elastography) can provide information about lung structural property changes and, hence, may be of diagnostic value. In the present study, a comprehensive computational simulation (in silico) model was developed to simulate sound wave propagation in the airways, lung, and chest wall under normal and pneumothorax conditions. Experiments were carried out to validate the model. Here, sound waves with frequency content from 50 to 700 Hz were introduced into airways of five porcine subjects via an endotracheal tube, and transmitted waves were measured by scanning laser Doppler vibrometry at the chest wall surface. The computational model predictions of decreased sound transmission with pneumothorax were consistent with experimental measurements. The in silico model can also be used to visualize wave propagation inside and on the chest wall surface for other pulmonary pathologies, which may help in developing and interpreting diagnostic procedures that utilize sound and vibration.

Experimental and Computational Studies of Sound Transmission in a Branching Airway Network Embedded in a Compliant Viscoelastic Medium.

Breath sounds are often used to aid in the diagnosis of pulmonary disease. Mechanical and numerical models could be used to enhance our understanding of relevant sound transmission phenomena. Sound transmission in an airway mimicking phantom was investigated using a mechanical model with a branching airway network embedded in a compliant viscoelastic medium. The Horsfield self-consistent model for the bronchial tree was adopted to topologically couple the individual airway segments into the branching airway network. The acoustics of the bifurcating airway segments were measured by microphones and calculated analytically. Airway phantom surface motion was measured using scanning laser Doppler vibrometry. Finite element simulations of sound transmission in the airway phantom were performed. Good agreement was achieved between experiments and simulations. The validated computational approach can provide insight into sound transmission simulations in real lungs.

A model of lung parenchyma stress relaxation using fractional viscoelasticity.

Some pulmonary diseases and injuries are believed to correlate with lung viscoelasticity changes. Hence, a better understanding of lung viscoelastic models could provide new perspectives on the progression of lung pathology and trauma. In the presented study, stress relaxation measurements were performed to quantify relaxation behavior of pig lungs. Results have uncovered certain trends, including an initial steep decay followed by a slow asymptotic relaxation, which would be better described by a power law than exponential decay. The fractional standard linear solid (FSLS) and two integer order viscoelastic models - standard linear solid (SLS) and generalized Maxwell (GM) - were used to fit the stress relaxation curves; the FSLS was found to be a better fit. It is suggested that fractional order viscoelastic models, which have nonlocal, multi-scale attributes and exhibit power law behavior, better capture the lung parenchyma viscoelastic behavior.

Pneumothorax effects on pulmonary acoustic transmission.

Pneumothorax (PTX) is an abnormal accumulation of air between the lung and the chest wall. It is a relatively common and potentially life-threatening condition encountered in patients who are critically ill or have experienced trauma. Auscultatory signs of PTX include decreased breath sounds during the physical examination. The objective of this exploratory study was to investigate the changes in sound transmission in the thorax due to PTX in humans. Nineteen human subjects who underwent video-assisted thoracic surgery, during which lung collapse is a normal part of the surgery, participated in the study. After subjects were intubated and mechanically ventilated, sounds were introduced into their airways via an endotracheal tube. Sounds were then measured over the chest surface before and after lung collapse. PTX caused small changes in acoustic transmission for frequencies below 400 Hz. A larger decrease in sound transmission was observed from 400 to 600 Hz, possibly due to the stronger acoustic transmission blocking of the pleural air. At frequencies above 1 kHz, the sound waves became weaker and so did their changes with PTX. The study elucidated some of the possible mechanisms of sound propagation changes with PTX. Sound transmission measurement was able to distinguish between baseline and PTX states in this small patient group. Future studies are needed to evaluate this technique in a wider population.

Comparison of Poroviscoelastic Models for Sound and Vibration in the Lungs.

Noninvasive measurement of mechanical wave motion (sound and vibration) in the lungs may be of diagnostic value, as it can provide information about the mechanical properties of the lungs, which in turn are affected by disease and injury. In this study, two previously derived theoretical models of the vibroacoustic behavior of the lung parenchyma are compared: (1) a Biot theory of poroviscoelasticity and (2) an effective medium theory for compression wave behavior (also known as a "bubble swarm" model). A fractional derivative formulation of shear viscoelasticity is integrated into both models. A measurable "fast" compression wave speed predicted by the Biot theory formulation has a significant frequency dependence that is not predicted by the effective medium theory. Biot theory also predicts a slow compression wave. The experimentally measured fast compression wave speed and attenuation in a pig lung ex vivo model agreed well with the Biot theory. To obtain the parameters for the Biot theory prediction, the following experiments were undertaken: quasistatic mechanical indentation measurements were performed to estimate the lung static shear modulus; surface wave measurements were performed to estimate lung tissue shear viscoelasticity; and flow permeability was measured on dried lung specimens. This study suggests that the Biot theory may provide a more robust and accurate model than the effective medium theory for wave propagation in the lungs over a wider frequency range.

A comprehensive computational model of sound transmission through the porcine lung.

A comprehensive computational simulation model of sound transmission through the porcine lung is introduced and experimentally evaluated. This "subject-specific" model utilizes parenchymal and major airway geometry derived from x-ray CT images. The lung parenchyma is modeled as a poroviscoelastic material using Biot theory. A finite element (FE) mesh of the lung that includes airway detail is created and used in comsol FE software to simulate the vibroacoustic response of the lung to sound input at the trachea. The FE simulation model is validated by comparing simulation results to experimental measurements using scanning laser Doppler vibrometry on the surface of an excised, preserved lung. The FE model can also be used to calculate and visualize vibroacoustic pressure and motion inside the lung and its airways caused by the acoustic input. The effect of diffuse lung fibrosis and of a local tumor on the lung acoustic response is simulated and visualized using the FE model. In the future, this type of visualization can be compared and matched with experimentally obtained elastographic images to better quantify regional lung material properties to noninvasively diagnose and stage disease and response to treatment.