Injury risk prediction from computational simulations of ocular blast loading

A predictive Lagrangian–Eulerian finite element eye model was used to analyze 2.27 and 0.45 kg trinitrotoluene equivalent blasts detonated from 24 different locations. Free air and ground level blasts were simulated directly in front of the eye and at lateral offset locations with box, average, less protective, and more protective orbital anthropometries, resulting in 96 simulations. Injury risk curves were developed for hyphema, lens dislocation, retinal damage, and globe rupture from experimental and computational data to compute risk from corneoscleral stress and intra-ocular pressure computational outputs. Corneoscleral stress, intra-ocular pressure, and injury risks increased when the blast size was larger and located nearer to the eye. Risks ranged from 20–100 % for hyphema, 1–100 % for lens dislocation, 2–100 % for retinal damage, and 0–98 % for globe rupture depending on the blast condition. Orbital geometry affected the stresses, pressures, and associated ocular injury risks of the blast conditions simulated. Orbital geometries that more fully surrounded the eye such as the more protective orbit tended to produce higher corneoscleral stresses and compression of the eye against the surrounding rigid orbit contributing to high stresses as the blast wave propagated. However, the more protective orbit tended to produce lower intra-ocular pressures in comparison with the other three orbital geometries which may indicate that the more protective orbit inhibits propagation of the blast wave and reduces ocular loading. Results of this parametric computational study of ocular blast loading are valuable to the design of eye protection equipment and the mitigation of blast-related eye injuries.     

A computational model of blast loading on the human eye

Ocular injuries from blast have increased in recent wars, but the injury mechanism associated with the primary blast wave is unknown. We employ a three-dimensional fluid–structure interaction computational model to understand the stresses and deformations incurred by the globe due to blast overpressure. Our numerical results demonstrate that the blast wave reflections off the facial features around the eye increase the pressure loading on and around the eye. The blast wave produces asymmetric loading on the eye, which causes globe distortion. The deformation response of the globe under blast loading was evaluated, and regions of high stresses and strains inside the globe were identified. Our numerical results show that the blast loading results in globe distortion and large deviatoric stresses in the sclera. These large deviatoric stresses may be indicator for the risk of interfacial failure between the tissues of the sclera and the orbit.     

Modelling human eye under blast loading

Primary blast injury (PBI) is the general term that refers to injuries resulting from the mere interaction of a blast wave with the body. Although few instances of primary ocular blast injury, without a concomitant secondary blast injury from debris, are documented, some experimental studies demonstrate its occurrence. In order to investigate PBI to the eye, a finite element model of the human eye using simple constitutive models was developed. The material parameters were calibrated by a multi-objective optimisation performed on available eye impact test data. The behaviour of the human eye and the dynamics of mechanisms occurring under PBI loading conditions were modelled. For the generation of the blast waves, different combinations of explosive (trinitrotoluene) mass charge and distance from the eye were analysed. An interpretation of the resulting pressure, based on the propagation and reflection of the waves inside the eye bulb and orbit, is proposed. The peculiar geometry of the bony orbit (similar to a frustum cone) can induce a resonance cavity effect and generate a pressure standing wave potentially hurtful for eye tissues.     

Numerical simulation of wave propagation in a realistic model of the human external ear

In this study, a numerical investigation is performed to evaluate the effects of high-pressure sinusoidal and blast wave's propagation around and inside of a human external ear. A series of computed tomography images are used to reconstruct a realistic three-dimensional (3D) model of a human ear canal and the auricle. The airflow field is then computed by solving the governing differential equations in the time domain using a computational fluid dynamics software. An unsteady algorithm is used to obtain the high-pressure wave propagation throughout the ear canal which is validated against the available analytical and numerical data in literature. The effects of frequency, wave shape, and the auricle on pressure distribution are then evaluated and discussed. The results clearly indicate that the frequency plays a key role on pressure distribution within the ear canal. At 4 kHz frequency, the pressure magnitude is much more amplified within the ear canal than the frequencies of 2 and 6 kHz, for the incident wave angle of 90° investigated in this study, attributable to the ‘4-kHz notch’ in patients with noise-induced hearing loss. According to the results, the pressure distribution patterns at the ear canal are very similar for both sinusoidal pressure waveform with the frequency of 2 kHz and blast wave. The ratio of the peak pressure value at the eardrum to that at the canal entrance increases from about 8% to 30% as the peak pressure value of the blast wave increases from 5 to 100 kPa for the incident wave angle of 90° investigated in this study. Furthermore, incorporation of the auricle to the ear canal model is associated with centerline pressure magnitudes of about 50% and 7% more than those of the ear canal model without the auricle throughout the ear canal for sinusoidal and blast waves, respectively, without any significant effect on pressure distribution pattern along the ear canal for the incident wave angle of 90° investigated in this study.     

Mesh considerations for finite element blast modelling in biomechanics

Finite element (FE) modelling is a popular tool for studying human body response to blast exposure. However, blast modelling is a complex problem owing to more numerous fluid–structure interactions (FSIs) and the high–frequency loading that accompanies blast exposures. This study investigates FE mesh design for blast modelling using a sphere in a closed-ended shock tube meshed with varying element sizes using both tetrahedral and hexahedral elements. FSI was consistent for sphere-to-fluid element ratios between 0.25 and 4, and acceleration response was similar for both element types (R² = 0.997). Tetrahedral elements were found to become increasingly volatile following shock loading, causing higher pressures and stresses than predicted with the hexahedral elements. Deviatoric stress response was dependent on the sphere mesh size (p < 0.001), while the pressure response was dependent on the shock tube mesh size (p < 0.001). The results of this study highlight the necessity for mesh sensitivity analysis in blast models.     

Dynamic material properties of the human sclera

As a result of trauma, approximately 30,000 people become blind in one eye every year in the United States. A common injury prediction tool is computational modeling, which requires accurate material properties to produce reliable results. Therefore, the purpose of this study was to determine the dynamic material properties of the human sclera. A high-rate pressurization system was used to create dynamic pressure to the point of rupture in 12 human eyes. Measurements were obtained for the internal pressure, the diameter of the globe, the thickness of the sclera, and the changing coordinates of the optical markers using high-rate video. A relationship between true stress and true strain was determined for the sclera tissue in two directions. It was found that the average maximum true stress was 13.89±4.81 MPa for both the equatorial and meridional directions, the average maximum true strain along the equator was 0.041±0.014, and the average maximum true strain along the meridian was 0.058±0.018. Results show a significant difference in the maximum strain in the equatorial and meridional directions (p=0.02). In comparing these data with previous studies, it is concluded that the human sclera is both anisotropic and viscoelastic. The dynamic material properties presented in this study can be used for advanced models of the human eye to help prevent eye injuries in the future.     

Blast related neurotrauma: a review of cellular injury

Historically, blast overpressure is known to affect primarily gas-containing organs such as the lung and ear. More recent interests focus on its ability to cause damage to solid organs such as the brain, resulting in neurological disorders. Returning veterans exposed to blast but without external injuries are being diagnosed with mild traumatic brain injury (Warden 2006) and with cortical dysfunction (Cernak et al 1999). Decades of studies have been conducted to elucidate the effects of primary blast wave on the central nervous system. These studies were mostly concerned with systemic effects (Saljo et al 2000-2003; Kaur et al 1995-1997, 1999; Cernak et al 1996, 2001). The molecular mechanism of blast-induced neurotrauma is still poorly understood. This paper reviews studies related to primary blast injury to the nervous system, particularly at the cellular level. It starts with a general discussion of primary blast injury and blast wave physics, followed by a review of the literature related to 1) the blast wave/body interaction, 2) injuries to the peripheral nervous system, 3) injuries to the central nervous system, and 4) injury criteria. Finally, some of our preliminary data on cellular injury from in vitro and in vivo studies are presented. Specifically, we report on the effects of overpressure on astrocytes. In the discussion, possible mechanisms of blast-related brain injury are discussed, as well as the concerns and limitations of the published studies. A clearer understanding of the injury mechanisms at both the molecular and macroscopic (organ) level will lead to the development of new treatment, diagnosis and preventive measures.     

Effect of exercise on blood flow through the aortic valve: a combined clinical and numerical study

The aim of this study was to measure the cardiac output and stroke volume for a healthy subject by coupling an echocardiogram Doppler (echo-Doppler) method with a fluid–structure interaction (FSI) simulation at rest and during exercise. Blood flow through aortic valve was measured by Doppler flow echocardiography. Aortic valve geometry was calculated by echocardiographic imaging. An FSI simulation was performed, using an arbitrary Lagrangian–Eulerian mesh. Boundary conditions were defined by pressure loads on ventricular and aortic sides. Pressure loads applied brachial pressures with (stage 1) and without (stage 2) differences between brachial, central and left ventricular pressures. FSI results for cardiac output were 15.4% lower than Doppler results for stage 1 (r = 0.999). This difference increased to 22.3% for stage 2. FSI results for stroke volume were undervalued by 15.3% when compared to Doppler results at stage 1 and 26.2% at stage 2 (r = 0.94). The predicted mean backflow of blood was 4.6%. Our results show that numerical methods can be combined with clinical measurements to provide good estimates of patient-specific cardiac output and stroke volume at different heart rates.     

Epidemiology of Pediatric Ocular Trauma in the Chaoshan Region,China, 2001–2010

Background

Ocular trauma is the leading cause of monocular visual disability and noncongenital unilateral blindness in children. This study describes the epidemiology and medical care associated with nonfatal pediatric (≤17 years of age) eye injury-related hospitalization in the largest industrial base for plastic toy production in China.

Methods

A population-based retrospective study of patients hospitalized for ocular and orbital trauma in the ophthalmology departments of 3 major tertiary hospitals from 1st January 2001 to 31st December 2010 was performed.

Results

The study included 1035 injured eyes from 1018 patients over a 10-year period: 560 (54.1%) eyes exhibited open globe injuries, 402 (38.8%) eyes suffered closed globe injuries, 10 (1.0%) eyes suffered chemical injuries and 8 (0.8%) eyes exhibited thermal injuries, representing an average annual hospitalization rate of 0.37 per 10,000 (95% confidence interval [CI], 0.36–0.38) due to pediatric eye injury in the Chaoshan region. The mean patient age was 9.2±4.4 years with a male-to-female ratio of 3.3∶1 (P = 0.007). Children aged 6 to 11 years accounted for the highest percentage (40.8%, 416/1018) of hospitalization, 56.7% (236/416) of whom were hospitalized for open globe wounds. Injury occurred most frequently at home (73.1%). Open globe wounds cost the single most expensive financial burden (60.8%) of total charges with $998±702 mean charges per hospitalization.

Conclusions

Open globe wounds occurred at home are earmarked for the priorities to prevention strategies. Higher public awareness of protecting primary schoolchildren from home-related eye injuries should be strengthened urgently by legislation or regulation since the traditional industrial mode seems to remain the pattern for the foreseeable future. Further research that provide detailed information on the specific inciting agents of pediatric eye injuries are recommended for facilitating the development and targeting of appropriate injury prevention initiatives.     

New head models extracted from thermal infrared (IR) images for dosimetry computations

In electromagnetic dosimetry, anatomical human models are commonly obtained by segmentation of magnetic resonance imaging or computed tomography scans. In this paper, a human head model extracted from thermal infrared images is examined in terms of its applicability to specific absorption rate (SAR) calculations. Since thermal scans are two-dimensional (2D) representation of surface temperature, this allows researchers to overcome the extensive computational demand associated with 3D simulation. The numerical calculations are performed using the finite-difference time-domain method with mesh sizes of 2 mm at 900 MHz plane wave irradiation. The power density of the incident plane wave is assumed to be 10 W/m². Computations were compared with a realistic anatomical head model. The results show that although there were marked differences in the local SAR distribution in the various tissues in the two models, the 1 g peak SAR values are approximately similar in the two models.     

Multiphysics simulation of left ventricular filling dynamics using fluid-structure interaction finite element method

To relate the subcellular molecular events to organ level physiology in heart, we have developed a three-dimensional finite-element-based simulation program incorporating the cellular mechanisms of excitation-contraction coupling and its propagation, and simulated the fluid-structure interaction involved in the contraction and relaxation of the human left ventricle. The FitzHugh-Nagumo model and four-state model representing the cross-bridge kinetics were adopted for cellular model. Both ventricular wall and blood in the cavity were modeled by finite element mesh. An arbitrary Lagrangian Eulerian finite element method with automatic mesh updating has been formulated for large domain changes, and a strong coupling strategy has been taken. Using electrical analog of pulmonary circulation and left atrium as a preload and the windkessel model as an afterload, dynamics of ventricular filling as well as ejection was simulated. We successfully reproduced the biphasic filling flow consisting of early rapid filling and atrial contraction similar to that reported in clinical observation. Furthermore, fluid-structure analysis enabled us to analyze the wave propagation velocity of filling flow. This simulator can be a powerful tool for establishing a link between molecular abnormality and the clinical disorder at the macroscopic level.     

African clawed toads (Xenopus laevis) sense the distance of lateral line stimuli

Sighted African clawed toads use their lateral lines to detect stimulus distance, although accuracy and precision are poorer than for stimulus direction. Single surface wave trains elicited discrete turns and/or swims towards the wave origin. Most responses were brief, ending with the toad stationary (70 % overall; 54–86 % individual toads) or pausing before turning away (11 %; 1–24 %). Lunges or capturing movements with the arms (13 %; 10–22 %) also indicated where toads expected to find prey. Overall, 94 % (88–100 %) of oriented responses had well-defined endpoints. Swim distance—measured as means, medians, and upper and lower quartiles—and the number of bilateral leg kicks increased with stimulus distance. Swim distance also depended upon stimulus angle due to features of turning. Most responses (81 %; 62–92 %) ended short of the wave origin. Regression slopes were 0.45 ± 0.04 mm/mm for stimulus distances up to 85 mm (ca. 2–3x body lengths), 0.16 ± 0.07 mm/mm for distances of 85–130 mm, and non-significant for larger distances to 220 mm. Slopes were steeper for responses that included lunges or capture movements. In only 15 % (3–26 %) of responses were both turn direction and swim distance sufficiently accurate for the toad to sweep through the wave origin.     

Reconstruction of internal orbital fractures with Vitallium mesh

Trauma to the face frequently results in internal orbital fractures that may produce large orbital defects involving multiple walls. Accurate anatomic reconstruction of the bony orbit is essential to maintain normal appearance and function of the eye following such injuries. Autogenous bone grafts do not always produce predictable long-term support of the globe. Displacement and varying amounts of bone-graft resorption can lead to enophthalmos. This study examines the use of Vitallium mesh in the acute reconstruction of internal orbital defects. Fifty-four patients with 66 orbits underwent reconstruction of internal orbital defects with Vitallium mesh. Associated fractures were anatomically reduced and rigidly fixed. Forty-six patients and 57 orbits had adequate follow-up for analysis of results. The average follow-up was 9 months, with 85 percent of the patients followed 6 months or longer. There were no postoperative orbital infections, and none of the Vitallium mesh required removal. Large internal orbital defects can be reconstructed using Vitallium mesh with good results and little risk of infection. Vitallium mesh appears to be well tolerated in spite of free communication with the sinuses. Stable reconstruction of the internal orbit can be achieved and predictable eye position maintained without donor-site morbidity.     

Three-dimensional Particle Tracking Velocimetry for Turbulence Applications: Case of a Jet Flow

3D-PTV is a quantitative flow measurement technique that aims to track the Lagrangian paths of a set of particles in three dimensions using stereoscopic recording of image sequences. The basic components, features, constraints and optimization tips of a 3D-PTV topology consisting of a high-speed camera with a four-view splitter are described and discussed in this article. The technique is applied to the intermediate flow field (5 <x/d <25) of a circular jet at Re ≈ 7,000. Lagrangian flow features and turbulence quantities in an Eulerian frame are estimated around ten diameters downstream of the jet origin and at various radial distances from the jet core. Lagrangian properties include trajectory, velocity and acceleration of selected particles as well as curvature of the flow path, which are obtained from the Frenet-Serret equation. Estimation of the 3D velocity and turbulence fields around the jet core axis at a cross-plane located at ten diameters downstream of the jet is compared with literature, and the power spectrum of the large-scale streamwise velocity motions is obtained at various radial distances from the jet core.     

Explosive Blast Loading on Human 3D Aggregate Minibrains

The effects of primary explosive blast on brain tissue still remain mostly unknown. There are few in vitro models that use real explosives to probe the mechanisms of injury at the cellular level. In this work, 3D aggregates of human brain cells or brain microphysiological system were exposed to military explosives at two different pressures (50 and 100 psi). Results indicate that membrane damage and oxidative stress increased with blast pressure, but cell death remained minimal.     

Characterisation of a phenomenological model for commercial pneumatic muscle actuators

This study focuses on the parameter characterisation of a three-element phenomenological model for commercially available pneumatic muscle actuators (PMAs). This model consists of a spring, damping and contractile element arranged in parallel. Data collected from static loading, contraction and relaxation experiments were fitted to theoretical solutions of the governing equation for the three-element model resulting in prediction profiles for the spring, damping and contractile force coefficient. For the spring coefficient, K N/mm, the following relationships were found: K = 32.7 ? 0.0321P for 150 ≤ P ≤ 314 kPa and K = 17 + 0.0179P for 314 ≤ P ≤ 550 kPa. For the damping coefficient, B Ns/mm, the following relationship was found during contraction: B = 2.90 for 150 ≤ P ≤ 550 kPa. During relaxation, B = 1.57 for 150 ≤ P ≤ 372 kPa and B = 0.311 + 0.00338P for 372 ≤ P ≤ 550. The following relationship for the contractile force coefficient, F _ce N, was also determined: F _ce = 2.91P+44.6 for 150 ≤ P ≤ 550 kPa. The model was then validated by reasonably predicting the response of the PMA to a triangular wave input in pressure under a constant load on a dynamic test station.