Influence of Residual Stress Field on the Fatigue Crack Propagation in Prestressing Steel Wires

This paper deals with the effect of several residual stress profiles on the fatigue crack propagation in prestressing steel wires subjected to tension loading or bending moment. To this end, a computer program was developed to evaluate the crack front evolution on the basis of the Walker law. Results demonstrate that the absence of residual stresses makes the crack propagate towards a preferential crack path. When surface residual stresses are tensile and, correspondingly, core residual stresses are compressive, the fatigue crack fronts rapidly converge towards a quasi-straightshape. When surface residual stresses are compressive, with their corresponding tensile stresses in the core area, a preferential crack path also appears.

Transcranial Assessment and Visualization of Acoustic Cavitation: Modeling and Experimental Validation

The interaction of ultrasonically-controlled microbubble oscillations with tissues and biological media has been shown to induce a wide range of bioeffects that may have significant impact on therapy and diagnosis of brain diseases and disorders. However, the inherently non-linear microbubble oscillations combined with the micrometer and microsecond scales involved in these interactions and the limited methods to assess and visualize them transcranially hinder both their optimal use and translation to the clinics. To overcome these challenges, we present a framework that combines numerical simulations with multimodality imaging to assess and visualize the microbubble oscillations transcranially. In the present work, microbubble oscillations were studied with an integrated US and MR imaging guided clinical FUS system. A high-resolution brain CT scan was also co-registered to the US and MR images and the derived acoustic properties were used as inputs to two- and three-dimensional Finite Difference Time Domain simulations that matched the experimental conditions and geometry. Synthetic point sources by either a Gaussian function or the output of a microbubble dynamics model were numerically excited and propagated through the skull towards a virtual US imaging array. Using passive acoustic mapping (PAM) that was refined to incorporate variable speed of sound, we were able to correct the aberrations introduced by the skull and substantially improve the PAM resolution. The good agreement between the simulations incorporating microbubble emissions and experimentally-determined PAMs suggest that this integrated approach can provide a clinically-relevant framework and more control over this nonlinear and dynamic process.

Biomechanical rationale of coronary artery bypass grafting of multivessel disease

The biomechanical model of human coronary arteries was modified for improving the quality of diagnosis and surgical treatment for coronary heart disease. The problem of hemodynamics in the left coronary artery with multivessel bed disease - 45% stenosis of the anterior descending branch and 75% stenosis of the circumflex branch - was particularly considered. Numerical simulation of the coronary arterial bypass of the main trunk was carried out to estimate the functional condition of the coronary arteries after restoring myocardial blood supply by surgery.

Evaluating simulant materials for understanding cranial backspatter from a ballistic projectile

In cranial wounds resulting from a gunshot, the study of backspatter patterns can provide information about the actual incidents by linking material to surrounding objects. This study investigates the physics of backspatter from a high-speed projectile impact and evaluates a range of simulant materials using impact tests. Next, we evaluate a mesh-free method called smoothed particle hydrodynamics (SPH) to model the splashing mechanism during backspatter. The study has shown that a projectile impact causes fragmentation at the impact site, while transferring momentum to fragmented particles. The particles travel along the path of least resistance, leading to partial material movement in the reverse direction of the projectile motion causing backspatter. Medium-density fiberboard is a better simulant for a human skull than polycarbonate, and lorica leather is a better simulant for a human skin than natural rubber. SPH is an effective numerical method for modeling the high-speed impact fracture and fragmentations.

Modeling of time dependent localized flow shear stress and its impact on cellular growth within additive manufactured titanium implants

Bone augmentation implants are porous to allow cellular growth, bone formation and fixation. However, the design of the pores is currently based on simple empirical rules, such as minimum pore and interconnects sizes. We present a three-dimensional (3D) transient model of cellular growth based on the Navier-Stokes equations that simulates the body fluid flow and stimulation of bone precursor cellular growth, attachment, and proliferation as a function of local flow shear stress. The model's effectiveness is demonstrated for two additive manufactured (AM) titanium scaffold architectures. The results demonstrate that there is a complex interaction of flow rate and strut architecture, resulting in partially randomized structures having a preferential impact on stimulating cell migration in 3D porous structures for higher flow rates. This novel result demonstrates the potential new insights that can be gained via the modeling tool developed, and how the model can be used to perform what-if simulations to design AM structures to specific functional requirements.

Biogeochemical control of marine productivity in the Mediterranean Sea during the last 50 years

The temporal dynamics of biogeochemical variables derived from a coupled 3-D model of the Mediterranean Sea are evaluated for the last 50 years (1960-2010) against independent data on fisheries catch per unit effort (CPUE) for the same time period. Concordant patterns are found in the time series of all of the biological variables (from the model and from fisheries statistics), with low values at the beginning of the series, a later increase, with maximum levels reached at the end of the 1990s, and a posterior stabilization. Spectral analysis of the annual biological time series reveals coincident low-frequency signals in all of them. The first, more energetic signal peaks around the year 2000, while the second, less energetic signal peaks near 1982. Almost identical low-frequency signals are found in the nutrient loads of the rivers and in the integrated nutrient levels in the surface marine ecosystem. Nitrate concentration shows a maximum level in 1998, with a later stabilization to present-day values, coincident with the first low-frequency signal found in the biological series. Phosphate shows maximum concentrations around 1982 and a posterior sharp decline, in concordance with the second low-frequency signal observed in the biological series. That result seems to indicate that the control of marine productivity (plankton to fish) in the Mediterranean is principally mediated through bottom-up processes that could be traced back to the characteristics of riverine discharges. The high sensitivity of CPUE time series to environmental conditions might be another indicator of the overexploitation of this marine ecosystem.

KEY POINTS

Biogeochemical evolution of the Mediterranean over the past 50 yearsRiver nutrient loads drive primary and secondary productionsStrong link between low trophic levels and fisheries.

Evaluation of the biodegradation of Alaska North Slope oil in microcosms using the biodegradation model BIOB

We present the details of a numerical model, BIOB that is capable of simulating the biodegradation of oil entrapped in the sediment. The model uses Monod kinetics to simulate the growth of bacteria in the presence of nutrients and the subsequent consumption of hydrocarbons. The model was used to simulate experimental results of Exxon Valdez oil biodegradation in laboratory columns (Venosa et al., 2010). In that study, samples were collected from three different islands: Eleanor Island (EL107), Knight Island (KN114A), and Smith Island (SM006B), and placed in laboratory microcosms for a duration of 168 days to investigate oil bioremediation through natural attenuation and nutrient amendment. The kinetic parameters of the BIOB model were estimated by fitting to the experimental data using a parameter estimation tool based on Genetic Algorithms (GA). The parameter values of EL107 and KN114A were similar whereas those of SM006B were different from the two other sites; in particular biomass growth at SM006B was four times slower than at the other two islands. Grain size analysis from each site revealed that the specific surface area per unit mass of sediment was considerably lower at SM006B, which suggest that the surface area of sediments is a key control parameter for microbial growth in sediments. Comparison of the BIOB results with exponential decay curves fitted to the data indicated that BIOB provided better fit for KN114A and SM006B in nutrient amended treatments, and for EL107 and KN114A in natural attenuation. In particular, BIOB was able to capture the initial slow biodegradation due to the lag phase in microbial growth. Sensitivity analyses revealed that oil biodegradation at all three locations were sensitive to nutrient concentration whereas SM006B was sensitive to initial biomass concentration due to its slow growth rate. Analyses were also performed to compare the half-lives of individual compounds with that of the overall polycyclic aromatic hydrocarbons (PAHs).

Numerical Modeling Calcium and CaMKII Effects in the SA Node

Sinoatrial node (SAN) is the primary heart pacemaker which initiates each heartbeat under normal conditions. Numerous experimental data have demonstrated that Ca(2+-) and CaMKII-dependent processes are crucially important for regulation of SAN cells. However, specific mechanisms of this regulation and their relative contribution to pacemaker function remain mainly unknown. Our review summarizes available data and existing numerical modeling approaches to understand Ca(2+) and CaMKII effects on the SAN. Data interpretation and future directions to address the problem are given within the coupled-clock theory, i.e., a modern view on the cardiac pacemaker cell function generated by a system of sarcolemmal and intracellular proteins.

Modeling of electric field distribution in tissues during electroporation

BACKGROUND

Electroporation based therapies and treatments (e.g. electrochemotherapy, gene electrotransfer for gene therapy and DNA vaccination, tissue ablation with irreversible electroporation and transdermal drug delivery) require a precise prediction of the therapy or treatment outcome by a personalized treatment planning procedure. Numerical modeling of local electric field distribution within electroporated tissues has become an important tool in treatment planning procedure in both clinical and experimental settings. Recent studies have reported that the uncertainties in electrical properties (i.e. electric conductivity of the treated tissues and the rate of increase in electric conductivity due to electroporation) predefined in numerical models have large effect on electroporation based therapy and treatment effectiveness. The aim of our study was to investigate whether the increase in electric conductivity of tissues needs to be taken into account when modeling tissue response to the electroporation pulses and how it affects the local electric distribution within electroporated tissues.

METHODS

We built 3D numerical models for single tissue (one type of tissue, e.g. liver) and composite tissue (several types of tissues, e.g. subcutaneous tumor). Our computer simulations were performed by using three different modeling approaches that are based on finite element method: inverse analysis, nonlinear parametric and sequential analysis. We compared linear (i.e. tissue conductivity is constant) model and non-linear (i.e. tissue conductivity is electric field dependent) model. By calculating goodness of fit measure we compared the results of our numerical simulations to the results of in vivo measurements.

RESULTS

The results of our study show that the nonlinear models (i.e. tissue conductivity is electric field dependent: σ(E)) fit experimental data better than linear models (i.e. tissue conductivity is constant). This was found for both single tissue and composite tissue. Our results of electric field distribution modeling in linear model of composite tissue (i.e. in the subcutaneous tumor model that do not take into account the relationship σ(E)) showed that a very high electric field (above irreversible threshold value) was concentrated only in the stratum corneum while the target tumor tissue was not successfully treated. Furthermore, the calculated volume of the target tumor tissue exposed to the electric field above reversible threshold in the subcutaneous model was zero assuming constant conductivities of each tissue.Our results also show that the inverse analysis allows for identification of both baseline tissue conductivity (i.e. conductivity of non-electroporated tissue) and tissue conductivity vs. electric field (σ(E)) of electroporated tissue.

CONCLUSION

Our results of modeling of electric field distribution in tissues during electroporation show that the changes in electrical conductivity due to electroporation need to be taken into account when an electroporation based treatment is planned or investigated. We concluded that the model of electric field distribution that takes into account the increase in electric conductivity due to electroporation yields more precise prediction of successfully electroporated target tissue volume. The findings of our study can significantly contribute to the current development of individualized patient-specific electroporation based treatment planning.

A numerical investigation of vapor intrusion--the dynamic response of contaminant vapors to rainfall events

The U.S. government and various agencies have published guidelines for field investigation of vapor intrusion, most of which suggest soil gas sampling as an integral part of the investigation. Contaminant soil gas data are often relatively more stable than indoor air vapor concentration measurements, but meteorological conditions might influence soil gas values. Although a few field and numerical studies have considered some temporal effects on soil gas vapor transport, a full explanation of the contaminant vapor concentration response to rainfall events is not available. This manuscript seeks to demonstrate the effects on soil vapor transport during and after different rainfall events, by applying a coupled numerical model of fluid flow and vapor transport. Both a single rainfall event and seasonal rainfall events were modeled. For the single rainfall event models, the vapor response process could be divided into three steps: namely, infiltration, water redistribution, and establishment of a water lens atop the groundwater source. In the infiltration step, rainfall intensity was found to determine the speed of the wetting front and wash-out effect on the vapor. The passage of the wetting front led to an increase of the vapor concentration in both the infiltration and water redistribution steps and this effect is noted at soil probes located 1m below the ground surface. When the mixing of groundwater with infiltrated water was not allowed, a clean water lens accumulated above the groundwater source and led to a capping effect which can reduce diffusion rates of contaminant from the source. Seasonal rainfall with short time intervals involved superposition of the individual rainfall events. This modeling results indicated that for relatively deeper soil that the infiltration wetting front could not flood, the effects were damped out in less than a month after rain; while in the long term (years), possible formation of a water lens played a larger role in determining the vapor intrusion risk. In addition, soil organic carbon retarded the transport process, and damped the contaminant concentration fluctuations.

Characterization of nonlinear ultrasound fields of 2D therapeutic arrays

A current trend in high intensity focused ultrasound (HIFU) technologies is to use 2D focused phased arrays that enable electronic steering of the focus, beamforming to avoid overheating of obstacles (such as ribs), and better focusing through inhomogeneities of soft tissue using time reversal methods. In many HIFU applications, the acoustic intensity in situ can reach thousands of W/cm2 leading to nonlinear propagation effects. At high power outputs, shock fronts develop in the focal region and significantly alter the bioeffects induced. Clinical applications of HIFU are relatively new and challenges remain for ensuring their safety and efficacy. A key component of these challenges is the lack of standard procedures for characterizing nonlinear HIFU fields under operating conditions. Methods that combine low-amplitude pressure measurements and nonlinear modeling of the pressure field have been proposed for axially symmetric single element transducers but have not yet been validated for the much more complex 3D fields generated by therapeutic arrays. Here, the method was tested for a clinical HIFU source comprising a 256-element transducer array. A numerical algorithm based on the Westervelt equation was used to enable 3D full-diffraction nonlinear modeling. With the acoustic holography method, the magnitude and phase of the acoustic field were measured at a low power output and used to determine the pattern of vibrations at the surface of the array. This pattern was then scaled to simulate a range of intensity levels near the elements up to 10 W/cm2. The accuracy of modeling was validated by comparison with direct measurements of the focal waveforms using a fiber-optic hydrophone. Simulation results and measurements show that shock fronts with amplitudes up to 100 MPa were present in focal waveforms at clinically relevant outputs, indicating the importance of strong nonlinear effects in ultrasound fields generated by HIFU arrays.

Impact of climate change on ecological quality indicators and biogeochemical fluxes in the Baltic sea: a multi-model ensemble study

Multi-model ensemble simulations using three coupled physical-biogeochemical models were performed to calculate the combined impact of projected future climate change and plausible nutrient load changes on biogeochemical cycles in the Baltic Sea. Climate projections for 1961-2099 were combined with four nutrient load scenarios ranging from a pessimistic business-as-usual to a more optimistic case following the Helsinki Commission's (HELCOM) Baltic Sea Action Plan (BSAP). The model results suggest that in a future climate, water quality, characterized by ecological quality indicators like winter nutrient, summer bottom oxygen, and annual mean phytoplankton concentrations as well as annual mean Secchi depth (water transparency), will be deteriorated compared to present conditions. In case of nutrient load reductions required by the BSAP, water quality is only slightly improved. Based on the analysis of biogeochemical fluxes, we find that in warmer and more anoxic waters, internal feedbacks could be reinforced. Increased phosphorus fluxes out of the sediments, reduced denitrification efficiency and increased nitrogen fixation may partly counteract nutrient load abatement strategies.

Modeling the dispersion of Ambrosia artemisiifolia L. pollen with the model system COSMO-ART

Common ragweed (Ambrosia artemisiifolia L.) is a highly allergenic plant that is spreading throughout Europe. Ragweed pollen can be transported over large distances by the wind. Even low pollen concentrations of less than 10 pollen m(-3) can lead to health problems in sensitive persons. Therefore, forecasting the airborne concentrations of ragweed pollen is becoming more and more important for public health. The question remains whether distant pollen sources need to be considered in reliable forecasts. We used the extended numerical weather prediction system COSMO-ART to simulate the release and transport of ragweed pollen in central Europe. A pollen episode (September 12-16, 2006) in north-eastern Germany was modeled in order to find out where the pollen originated. For this purpose, several different source regions were taken into account and their individual impact on the daily mean pollen concentration and the performance of the forecast were studied with the means of a 2 × 2 contingency table and skill scores. It was found that the majority of the pollen originated in local areas, but up to 20% of the total pollen load came from distant sources in Hungary. It is concluded that long-distance transport should not be neglected when predicting pollen concentrations.

Evaluation of Common RF Coil Setups for MR Imaging at Ultrahigh Magnetic Field: A Numerical Study

This study is an evaluation of the ratio of electric field to magnetic field (E/B₁), specific absorption rate (SAR) and signal-to-noise ratio (SNR) generated by three different RF transceiver coil setups: surface coil, surface coil with shielding, and microstrip using a finite discrete time domain (FDTD) simulation in the presence of a head phantom. One of our main focuses in this study is to better understand coil designs that would improve patient safety at high fields by studying a coil type that may potentially minimize SAR while examining potential changes in SNR. In the presence of a human head load, the microstrip's E/B₁ ratio was on average smallest while its SAR was also on average smallest of the three setups, suggesting the microstrip may be a better RF coil choice for MRI concerning patient safety and parallel excitation applications than the other two coils. In addition, the study suggests that the microstrip also has a higher SNR compared with the other two coils demonstrating the possibility that the microstrip could lead to higher quality MRI images.

Three-dimensional modeling of the dynamics of therapeutic ultrasound contrast agents

A 3-D thick-shell contrast agent dynamics model was developed by coupling a finite volume Navier-Stokes solver and a potential boundary element method flow solver to simulate the dynamics of thick-shelled contrast agents subjected to pressure waves. The 3-D model was validated using a spherical thick-shell model validated by experimental observations. We then used this model to study shell break-up during nonspherical deformations resulting from multiple contrast agent interaction or the presence of a nearby solid wall. Our simulations indicate that the thick viscous shell resists the contrast agent from forming a re-entrant jet, as normally observed for an air bubble oscillating near a solid wall. Instead, the shell thickness varies significantly from location to location during the dynamics, and this could lead to shell break-up caused by local shell thinning and stretching.

Numerical Modeling of Silicon Photodiodes for High-Accuracy Applications Part II. Interpreting Oxide-Bias Experiments

The semiconductor device modeling program PC-1D and the programs that support its use in high-accuracy modeling of photodiodes, all of which were described in Part I of this series of papers, are used to simulate oxide-bias self-calibration experiments on three different types of silicon photodiodes. It is shown that these simulations can be used to determine photodiode characteristics, including the internal quantum efficiency for the different types of photodiodes. In the latter case, the simulations provide more accurate values than can be determined by using the conventional data reduction procedure, and an uncertainty estimate can be derived. Finally, it is shown that 0.9997 ± 0.0003 is a nominal value for the internal quantum efficiency of one type of photodiode over the 440 to 460 nm spectral region.