Modeling transcranial magnetic stimulation coil with magnetic cores

Objective.Accurate modeling of transcranial magnetic stimulation (TMS) coils with the magnetic core is largely an open problem since commercial (quasi) magnetostatic solvers do not output specific field characteristics (e.g. induced electric field) and have difficulties when incorporating realistic head models. Many open-source TMS softwares do not include magnetic cores into consideration. This present study reports an algorithm for modeling TMS coils with a (nonlinear) magnetic core and validates the algorithm through comparison with finite-element method simulations and experiments.Approach.The algorithm uses the boundary element fast multipole method applied to all facets of a tetrahedral core mesh for a single-state solution and the successive substitution method for nonlinear convergence of the subsequent core states. The algorithm also outputs coil inductances, with or without magnetic cores. The coil-core combination is solved only once i.e. before incorporating the head model. The resulting primary TMS electric field is proportional to the total vector potential in the quasistatic approximation; it therefore also employs the precomputed core magnetization.Main results.The solver demonstrates excellent convergence for typical TMS field strengths and for analyticalB-Happroximations of experimental magnetization curves such as Froelich's equation or an arctangent equation. Typical execution times are 1-3 min on a common multicore workstation. For a simple test case of a cylindrical core within a one-turn coil, our solver computed the small-signal inductance nearly identical to that from ANSYS Maxwell. For a multiturn rodent TMS coil with a core, the modeled inductance matched the experimental measured value to within 5%.Significance.Incorporating magnetic core in TMS coil design has advantages of field shaping and energy efficiency. Our software package can facilitate model-informed design of more efficiency TMS systems and guide selection of core material. These models can also inform dosing with existing clinical TMS systems that use magnetic cores.

Developments on Constitutive Material Model for Architectural Soda-Lime Silicate (SLS) Glass and Evaluation of Key Modelling Parameters

Architectural soda-lime silicate glass (SLS) is increasingly taking on complex shapes that require more detailed numerical analysis. Glass modeling is a thoroughly described topic with validated constitutive models. However, these models require a number of precise material parameters for SLS glass, and these are very sensitive to changes in glass composition. The currently available information is based on SLS glass tested in the late 1990s. As a result, most current publications are based on the above data. The object of this work was to analyze the available sources and update the information on selected key parameters for modeling. Using the currently utilized SLS glass in construction, the coefficient of thermal expansion (CTE), glass transition temperature, and the Young's modulus have been experimentally investigated. The updated material parameters will allow for more accurate modeling of the SLS glass currently used in construction, and in consequence will make the prototyping process for glass with complex geometries possible to be transferred from the production stage to the design stage, resulting in shorter production times.

Total ice content and lipid saturation determine adipose tissue cryolipolysis by injection of ice-slurry

OBJECTIVES

Cryolipolysis uses tissue cooling to solidify lipids, preferentially damaging lipid-rich cells. Topical cooling is popular for the reduction of local subcutaneous fat. Injection of biocompatible ice-slurry is a recently introduced alternative. We developed and verified a quantitative model that simulates the heat exchange and phase changes involved, offering insights into ice-slurry injection for treating subcutaneous fat.

METHODS

Finite element method was used to model the spatial and temporal progression of heat transfer between adipose tissue and injected ice-slurry, estimating dose-response relationships between properties of the slurry and size of tissue affected by cryolipolysis. Phase changes of both slurry and adipose tissue lipids were considered. An in vivo swine model was used to validate the numerical solutions. Oils with different lipid compositions were exposed to ice-slurry in vitro to evaluate the effects of lipid freezing temperature. Microscopy and nuclear magnetic resonance (NMR) were performed to detect lipid phase changes.

RESULTS

A ball of granular ice was deposited at the injection site in subcutaneous fat. Total injected ice content determines both the effective cooling region of tissue, and the duration of tissue cooling. Water's high latent heat of fusion enables tissue cooling long after slurry injection. Slurry temperature affects the rate of tissue cooling. In swine, when 30 ml slurry injection at -3.5°C was compared to 15 ml slurry injection at -4.8°C (both with the same total ice content), the latter led to almost twice faster tissue cooling. NMR showed a large decrease in diffusion upon lipid crystallization; saturated lipids with higher freezing temperatures were more susceptible to solidification after ice-slurry injection.

CONCLUSIONS

Total injected ice content determines both the volume of tissue treated by cryolipolysis and the cooling duration after slurry injection, while slurry temperature affects the cooling rate. Lipid saturation, which varies with diet and anatomic location, also has an important influence.

Microstructure-Based Computational Analysis of Deformation Stages of Rock-like Sandy-Cement Samples in Uniaxial Compression

This work presents a new finite-difference continuum damage mechanics approach for assessment of threshold stresses based on the mechanical response of a representative volume element of a sandy-cement rock-like material. An original experimental study allows validating the mathematical model. A new modification of the damage accumulation kinetic equation is proposed. Several approaches based on acoustic emission, instantaneous Poisson's ratio and reversal point method are employed to determine the threshold stresses. Relying on the numerical modeling of deformation and failure of model samples, the threshold stresses and the deformation stages are determined. The model predicts the crack initiation stress threshold with less than 10% error. The model prediction of the crack damage stress threshold corresponds to the upper boundary of the experimental range. The model predicts the peak stress threshold with less than 0.2% error in comparison with the average experimental peak stress. The results of numerical modeling are shown to correlate well with the available experimental and literature data and sufficiently complement them.

Validation of a cerebral hemodynamic model with personalized calibration in patients with aneurysmal subarachnoid hemorrhage

This study aims to validate a numerical model developed for assessing personalized circle of Willis (CoW) hemodynamics under pathological conditions. Based on 66 computed tomography angiography images, investigations were obtained from 43 acute aneurysmal subarachnoid hemorrhage (aSAH) patients from a local neurovascular center. The mean flow velocity of each artery in the CoW measured using transcranial Doppler (TCD) and simulated by the numerical model was obtained for comparison. The intraclass correlation coefficient (ICC) over all cerebral arteries for TCD and the numerical model was 0.88 (N = 561; 95% CI 0.84-0.90). In a subgroup of patients who had developed delayed cerebral ischemia (DCI), the ICC had decreased to 0.72 but remained constant with respect to changes in blood pressure, Fisher grade, and location of ruptured aneurysm. Our numerical model showed good agreement with TCD in assessing the flow velocity in the CoW of patients with aSAH. In conclusion, the proposed model can satisfactorily reproduce the cerebral hemodynamics under aSAH conditions by personalizing the numerical model with TCD measurements. Clinical trial registration: [http://www.trialregister.nl/], identifier [NL8114].

Isolation Improvement in Reflectarray Antenna-Based FMCW Radar Systems

This paper presents an optimization of reflectarray-based RF sensors for detecting UAV and human presence. Our previous human detection radar system adapted a center-fed reflectarray antenna to a commercially available radar system, successfully increasing the gains of the transmit (TX) and receive (RX) antennas by 21.18 dB and the range for detecting human targets 3.4 times. However, because the TX and RX antennas were placed in the focal point of the reflectarray, the TX signal reflected by the reflectarray was directly propagated into the RX antenna, causing desensitization or damage to the receiving circuit if high powers were used. To reduce this direct reflection, we propose a novel radar antenna configuration in which the TX and RX antennas are placed back-to-back with each other. In this configuration, the RX antenna does not directly face the reflectarray, thus direct path between the TX to RX through the reflectarray is removed. The results demonstrate that this approach achieves the optimum isolation level of 51.3 dB. With the reflectarray, the TX antenna gain increases to 30.6 dBi, but the RX antenna gain remains at 16 dBi since the RX antenna does not utilize the reflectarray. The TX and RX gain difference (14.6 dB) is a trade-off for good isolation and may be reduced by utilizing a high-gain receiver amplifier.

Results of Numerical Modeling of Blood Flow in the Internal Jugular Vein Exhibiting Different Types of Strictures

The clinical relevance of nozzle-like strictures in upper parts of the internal jugular veins remains unclear. This study was aimed at understanding flow disturbances caused by such stenoses. Computational fluid dynamics software, COMSOL Multiphysics, was used. Two-dimensional computational domain involved stenosis at the beginning of modeled veins, and a flexible valve downstream. The material of the venous valve was considered to be hyperelastic. In the vein models with symmetric 2-leaflets valve without upstream stenosis or with minor 30% stenosis, the flow was undisturbed. In the case of major 60% and 75% upstream stenosis, centerline velocity was positioned asymmetrically, and areas of reverse flow and flow separation developed. In the 2-leaflet models with major stenosis, vortices evoking flow asymmetry were present for the entire course of the model, while the valve leaflets were distorted by asymmetric flow. Our computational fluid dynamics modeling suggests that an impaired outflow from the brain through the internal jugular veins is likely to be primarily caused by pathological strictures in their upper parts. In addition, the jugular valve pathology can be exacerbated by strictures located in the upper segments of these veins.

Plasmon-Enhanced Ultraviolet Luminescence in Colloid Solutions and Nanostructures Based on Aluminum and ZnO Nanoparticles

Aluminum nanoparticles attract scientific interest as a promising low-cost material with strong plasmon resonance in the ultraviolet region, which can be used in various fields of photonics. In this paper, for the first time, ultraviolet luminescence of zinc oxide nanoparticles in colloid solutions and nanostructure films in the presence of plasmonic aluminum nanoparticles 60 nm in size with a metal core and an aluminum oxide shell were studied. Mixture colloids of ZnO and Al nanoparticles in isopropyl alcohol solution with concentrations from 0.022 to 0.44 g/L and 0.057 to 0.00285 g/L, correspondingly, were investigated. The enhancement of up to 300% of ZnO emission at 377 nm in colloids mixtures with metal nanoparticles due to formation of Al-ZnO complex agglomerates was achieved. Plasmon nanostructures with different configurations of layers, such as Al on the surface of ZnO, ZnO on Al, sandwich-like structure and samples prepared from a colloidal mixture of ZnO and Al nanoparticles, were fabricated by microplotter printing. We demonstrated that photoluminescence can be boosted 2.4-fold in nanostructures prepared from a colloidal mixture of ZnO and Al nanoparticles, whereas the sandwich-like structure gave only 1.1 times the amplification of luminescence. Calculated theoretical models of photoluminescence enhancement of ideal and weak emitters near aluminum nanoparticles of different sizes showed comparable results with the obtained experimental data.

Analysis of Stress and Displacement Fields in Prosthetic Crowns Made of Zirconium Dioxide Using Numerical Approach of Homogenization Hypothesis

The main goal of this paper is to analyze the stress and displacement fields in prosthetic crowns made of zirconium dioxide using the numerical approach of homogenization hypothesis. The simple engineering model is developed and applied in case of vertical forces. The model is a three-dimensional simulation of molars subjected to crushing, mastication, and clenching. Two basic approaches are considered: the single prosthetic crown on a single molar, and the prosthetic bridge on two molars. The distributions of material parameters are determined for the rigid support and the elastic gum structure of the homogenized properties. The crown on a single molar is analyzed in respect of caries, which are represented by weak material parameters. Irrespective of the problem, the maximal stresses are always insignificant compared to the compressive strength for enamel, dentin, periodontium, and zirconium dioxide. In case of caries, the maximal stresses are located at the contact surface caries/crown, whereas the displacement was higher than the same parameter without caries. The stresses inside the prosthetic bridge on two molars were comparable for elastic and rigid support, and located at the same areas. The molar displacement for elastic gum was higher than for the rigid base, and additionally supplemented by the displacement of the supporting structure.

Optimization of Transpedicular Electrode Insertion for Electroporation-Based Treatments of Vertebral Tumors

Electroporation-based treatments such as electrochemotherapy and irreversible electroporation ablation have sparked interest with respect to their use in medicine. Treatment planning involves determining the best possible electrode positions and voltage amplitudes to ensure treatment of the entire clinical target volume (CTV). This process is mainly performed manually or with computationally intensive genetic algorithms. In this study, an algorithm was developed to optimize electrode positions for the electrochemotherapy of vertebral tumors without using computationally intensive methods. The algorithm considers the electric field distribution in the CTV, identifies undertreated areas, and uses this information to iteratively shift the electrodes from their initial positions to cover the entire CTV. The algorithm performs successfully for different spinal segments, tumor sizes, and positions within the vertebra. The average optimization time was 71 s with an average of 4.9 iterations performed. The algorithm significantly reduces the time and expertise required to create a treatment plan for vertebral tumors. This study serves as a proof of concept that electrode positions can be determined (semi-)automatically based on the spatial information of the electric field distribution in the target tissue. The algorithm is currently designed for the electrochemotherapy of vertebral tumors via a transpedicular approach but could be adapted for other anatomic sites in the future.

Experimental and Numerical Determination of Strength Characteristics Related to Paraglider Wing with Fourier Transform Infrared Spectroscopy of Applied Materials

The aim of paper is to determine experimentally and numerically the strength characteristics related to the paraglider wing with Fourier transform infrared spectroscopy of applied materials. The applied method consists in theoretical modeling supplemented by the tests of material parameters. First, the set of 10 lightweight fabrics was selected for the tests; the samples are representative for these structures. The materials were tested using the spectroscopy to determine the FTIR spectra. The samples differ in the content of certain characteristic groups. Air permeability change of the materials was determined for the different pressure drops. The air permeability of almost all the analyzed samples was close to zero with the exception of only one material. The tensile strength and elongation at the break of samples were determined on the testing machine. The paraglider samples were characterized by slightly decreased mechanical properties compared to the parachute fabrics. The material characteristics determined during the tests are the input data for the theoretical analysis. The numerical model of the paraglider wing is based on a 3D geometry from previous research, but the stress, strain, and deformation were determined using the ANSYS Structural program and the finite elements method. To determine the strength correctly, we introduce two basic values: the absolute maximal and the representative values that are the biggest repetitive values of stress, strain, and deformation. The stress value was determined by the main factors: (i) the thinner the material, the bigger the stresses that were accumulated; (ii) the stronger the material, the bigger the stresses that were accumulated. The results are similar for all materials and differ mainly by the values. The biggest stresses were observed inside the material contacting the ribs, whereas the biggest deformation and strain were in the regions between ribs, and the smallest were in the contact areas with the fixed supports. Their highest intensity was observed on the leading edge of the paraglider. We conclude that the obtained stresses were far from the breaking level for the wing.

The effect of plaque morphology, material composition and microcalcifications on the risk of cap rupture: A structural analysis of vulnerable atherosclerotic plaques

Background

The mechanical rupture of an atheroma cap may initiate a thrombus formation, followed by an acute coronary event and death. Several morphology and tissue composition factors have been identified to play a role on the mechanical stability of an atheroma, including cap thickness, lipid core stiffness, remodeling index, and blood pressure. More recently, the presence of microcalcifications (μCalcs) in the atheroma cap has been demonstrated, but their combined effect with other vulnerability factors has not been fully investigated.

Materials and methods

We performed numerical simulations on 3D idealized lesions and a microCT-derived human coronary atheroma, to quantitatively analyze the atheroma cap rupture. From the predicted cap stresses, we defined a biomechanics-based vulnerability index (VI) to classify the impact of each risk factor on plaque stability, and developed a predictive model based on their synergistic effect.

Results

Plaques with low remodeling index and soft lipid cores exhibit higher VI and can shift the location of maximal wall stresses. The VI exponentially rises as the cap becomes thinner, while the presence of a μCalc causes an additional 2.5-fold increase in vulnerability for a spherical inclusion. The human coronary atheroma model had a stable phenotype, but it was transformed into a vulnerable plaque after introducing a single spherical μCalc in its cap. Overall, cap thickness and μCalcs are the two most influential factors of mechanical rupture risk.

Conclusions

Our findings provide supporting evidence that high risk lesions are non-obstructive plaques with softer (lipid-rich) cores and a thin cap with μCalcs. However, stable plaques may still rupture in the presence of μCalcs.

A New Mechanism for Early-Time Plasmaspheric Refilling: The Role of Charge Exchange Reactions in the Transport of Energy and Mass Throughout the Ring Current-Plasmasphere System

Cold H⁺ produced via charge exchange reactions between ring current ions and exospheric neutral hydrogen constitutes an additional source of cold plasma that further contributes to the plasmasphere and affects the plasma dynamics in the Earth's magnetosphere system; however, its production and associated effects on the plasmasphere dynamics have not been fully assessed and quantified. In this study, we perform numerical simulations mimicking an idealized three-phase geomagnetic storm to investigate the role of heavy ion composition in the ring current (O⁺ vs. N⁺) and exospheric neutral hydrogen density in the production of cold H⁺ via charge exchange reactions. It is found that ring current heavy ions produce more than 50% of the total cold H⁺ via charge exchange reactions, and energetic N⁺ is more efficient in producing cold H⁺ via charge exchange reactions than O⁺. Furthermore, the density structure of the cold H⁺ is highly dependent on the mass of the parent ion; that is, cold H⁺ deriving from charge exchange reactions involving energetic O⁺ with neutral hydrogen, populates the lower L-shells, while cold H⁺ deriving from charge exchange reactions involving energetic N⁺ with neutral hydrogen populates the higher L-shells. In addition, the density of cold H⁺ produced via charge exchange reactions involving N⁺ can be peak at values up to one order of magnitude larger than the local plasmaspheric density, suggesting that solely considering the supply of cold plasma from the ionosphere to the plasmasphere can lead to a significant underestimation of plasmasphere density.

Utilizing low oxygen to mitigate resistance of stored product insects to phosphine

BACKGROUND

Data are provided on the utilization of modified atmospheres, at a commercial scale, against stored product insect populations that are resistant to phosphine. The method is evaluated on different populations of two major stored-product beetle species, Rhyzopertha dominica and Oryzaephilus surinamensis. The trials were carried out in commercial facilities, in which nitrogen was introduced through an embedded nitrogen generator. Each chamber contained three or four pallets of either currants or herbs. A computational model was developed to evaluate the nitrogen concentration.

RESULTS

In most trials, 100% mortality was recorded for both beetle species and all populations, regardless of the temperature and exposure intervals tested. Control progeny production ranged between 20 and 45 adults per vial for R. dominica, and 29 and 27 adults per vial for O. surinamensis. Simulation results reveal that nitrogen can easily penetrate the currants, and its concentration is uniform (differences are below 1.5%) across the pallet. Additionally, the simulation model revealed that lower temperatures do not have an impact on the nitrogen concentration profiles.

CONCLUSIONS

The modified atmosphere applications evaluated here were proved to be effective for all populations, regardless of the level of resistance to phosphine, and any survival could be attributed to the short exposure intervals. Modified atmosphere applications can be effective at a considerably short exposure interval, even at 2.5 days, which is an incontestable advantage for the use of this method against insects, at exposures comparable with those of commercial fumigations. © 2022 Society of Chemical Industry.

Gene Electrotransfer into Mammalian Cells Using Commercial Cell Culture Inserts with Porous Substrate

Gene electrotransfer is one of the main non-viral methods for intracellular delivery of plasmid DNA, wherein pulsed electric fields are used to transiently permeabilize the cell membrane, allowing enhanced transmembrane transport. By localizing the electric field over small portions of the cell membrane using nanostructured substrates, it is possible to increase considerably the gene electrotransfer efficiency while preserving cell viability. In this study, we expand the frontier of localized electroporation by designing an electrotransfer approach based on commercially available cell culture inserts with polyethylene-terephthalate (PET) porous substrate. We first use multiscale numerical modeling to determine the pulse parameters, substrate pore size, and other factors that are expected to result in successful gene electrotransfer. Based on the numerical results, we design a simple device combining an insert with substrate containing pores with 0.4 µm or 1.0 µm diameter, a multiwell plate, and a pair of wire electrodes. We test the device in three mammalian cell lines and obtain transfection efficiencies similar to those achieved with conventional bulk electroporation, but at better cell viability and with low-voltage pulses that do not require the use of expensive electroporators. Our combined theoretical and experimental analysis calls for further systematic studies that will investigate the influence of substrate pore size and porosity on gene electrotransfer efficiency and cell viability.

Bond Behavior of a Bio-Aggregate Embedded in Cement-Based Matrix

This paper investigates the bond behavior between a bio-aggregate and a cement-based matrix. The experimental evaluation comprised physical, chemical, image, and mechanical characterization of the bio-aggregate. The image analyses about the bio-aggregate's outer structure provided first insights to understand the particularities of this newly proposed bio-aggregate for use in cementitious materials. A mineral aggregate (granitic rock), largely used as coarse aggregate in the Brazilian civil construction industry, was used as reference. The bond behavior of both aggregates was evaluated via pull-out tests. The results indicated that both aggregates presented a similar linear elastic branch up to each respective peak loads. The peak load magnitude of the mineral aggregate indicated a better chemical adhesion when compared to the bio-aggregate's. The post-peak behavior, however, indicated a smoother softening branch for the bio-aggregate, corroborated by the microscopy image analyses. Although further investigation is required, the macaúba crushed endocarp was found to be a thriving bio-material to be used as bio-aggregate.

Geometric and Scaling Effects in the Speed of Catalytic Enzyme Micropumps

Self-powered, biocompatible pumps in the nanometer to micron length scale have the potential to enable technology in several fields, including chemical analysis and medical diagnostics. Chemically powered, catalytic micropumps have been developed but are not able to function well in biocompatible environments due to their intolerance of salt solutions and the use of toxic fuels. In contrast, enzymatically powered catalytic pumps offer good biocompatibility, selectivity, and scalability, but their performance at length scales below a few millimeters, which is important to many of their possible applications, has not been well tested. Here, urease-based enzyme pumps of millimeter and micrometer dimensions were fabricated and studied. The scaling of the pumping velocity was measured experimentally and simulated by numerical modeling. Pumping speeds were analyzed accurately by eliminating Brownian noise from the data using enzyme patches between 5 mm and 350 μm in size. Pumping speeds of microns per second could be achieved with urease pumps and were fastest when the channel height exceeded the width of the catalytic pump patch. In all cases, pumping was weak when the dimensions of the patch were 100 μm or less. Experimental and simulation results were consistent with a density-driven pumping mechanism at all sizes studied and served as a framework for the in silico study of more complex two-dimensional (2D) and three-dimensional (3D) geometries. Attempts to create directional flow by juxtaposing inward and outward pumps were unsuccessful because of the symmetry of convection rolls produced by millimeter-size pump patches and the slow speeds of smaller pumps. However, simulations of a corrugated ratchet structure showed that directional pumping could be achieved with pump patches in the millimeter size range.

On the Numerical Modeling of Flax/PLA Bumper Beams

Significant progress has been made in green composites developing fully biodegradable composites made of microbially degradable polymers reinforced with natural fibers. However, an improvement in the development of numerical models to predict the damage of green composites is necessary to extend their use in industrial applications of structural responsibility. This paper is focused on developing a numerical model that can predict the failure modes of four types of bumper beams made of flax/PLA green composites with different cross sections. The predictions regarding energy absorption, contact force history, and extension of delamination were compared with experimental results to validate the FEM model, and both results revealed a good agreement. Finally, the FEM model was used to analyze the failure modes of the bumper beams as a function of the impact energy and cross-section roundness. The impact energy threshold defined as the maximum absorbed-energy capability of the beam match with the impact energy that produces delaminations extended through all the cross sections. Experimental and numerical results revealed that the threshold energy, where the maximum energy-absorption capability is reached, for Type A is over 60 J; for Type B and C is around 60 J; and for Type D is at 50 J. Since delamination is concentrated at the cross-section corners, the threshold energy decreases with the cross-section roundness because the higher the roundness ratio, the wider the delamination extension.

Characteristics of Earthquake Cycles: A Cross-Dimensional Comparison of 0D to 3D Numerical Models

High-resolution computer simulations of earthquake sequences in three or even two dimensions pose great demands on time and energy, making lower-cost simplifications a competitive alternative. We systematically study the advantages and limitations of simplifications that eliminate spatial dimensions in quasi-dynamic earthquake sequence models, from 3D models with a 2D fault plane down to 0D or 1D models with a 0D fault point. We demonstrate that, when 2D or 3D models produce quasi-periodic characteristic earthquakes, their behavior is qualitatively similar to lower-dimension models. Certain coseismic characteristics like stress drop and fracture energy are largely controlled by frictional parameters and are thus largely comparable. However, other observations are quantitatively clearly affected by dimension reduction. We find corresponding increases in recurrence interval, coseismic slip, peak slip velocity, and rupture speed. These changes are to a large extent explained by the elimination of velocity-strengthening patches that transmit tectonic loading onto the velocity-weakening fault patch, thereby reducing the interseismic stress rate and enhancing the slip deficit. This explanation is supported by a concise theoretical framework, which explains some of these findings quantitatively and effectively estimates recurrence interval and slip. Through accounting for an equivalent stressing rate at the nucleation size h* into 2D and 3D models, 0D or 1D models can also effectively simulate these earthquake cycle parameters. Given the computational efficiency of lower-dimensional models that run more than a million times faster, this paper aims to provide qualitative and quantitative guidance on economical model design and interpretation of modeling studies.

Numerical Parametric Study of Coda Wave Interferometry Sensitivity to Microcrack Change in a Multiple Scattering Medium

The expansion of cracks in 3D printing concrete materials may lead to structural failure, so it is essential to monitor crack propagation development. Coda wave interferometry (CWI) has been proven to be sensitive to microcracks, however, the evolution pattern of ultrasonic coda waves during crack growth is still not clear. This paper reports a numerical study of the sensitivity and feasibility of CWI for monitoring microcrack growth in heterogeneous materials. A two-phase concrete model, which contains microcracks with different angles and lengths, was developed using the finite element analysis software ABAQUS. The relative velocity change (Δv/v) and the decorrelation coefficient (Kd) at different crack increments were quantitatively analyzed. The numerical simulation results show that coda waves are sensitive to microcrack length as well as the crack angle. The Δv/v increases linearly with the increase of the length of a single microcrack, and the Kd could be linked to the crack length quadratically. Furthermore, a quantitative functional relationship between the CWI observations (Kd, Δv/v) and the angle of the crack to the source/receiver and the relative length growth of the crack are established. In addition, the nonlinear relationship between slope and angle can be fitted with a sinusoidal function. The reported results quantitatively assess the coda wave variation pattern during crack propagation, which is important for the promotion and application of CWI technology.

Investigation of Adhesion Properties of Tire-Asphalt Pavement Interface Considering Hydrodynamic Lubrication Action of Water Film on Road Surface

To obtain the tire−pavement peak adhesion coefficient under different road states, a field measurement and FE simulation were combined to analyze the tire−pavement adhesion characteristics in this study. According to the identified texture information, the power spectral distribution of the road surface was obtained using the MATLAB Program, and a novel tire hydroplaning FE model coupled with a textured pavement model was established in ABAQUS. Experimental results show that here exists an “anti-skid noncontribution area” for the insulation and lubrication of the water film. Driving at the limit speed of 120 km/h, the critical water film thickness for the three typical asphalt pavements during hydroplaning was as follows: AC pavement, 0.56 mm; SMA pavement, 0.76 mm; OGFC pavement, 1.5 mm. The road state could be divided into four parts dry state, wet sate, lubricated state, and ponding state. Under the dry road state, when the slip rate was around 15%, the adhesion coefficient reached the peak value, i.e., around 11.5% for the wet road state. The peak adhesion coefficient for the different asphalt pavements was in the order OGFC > SMA > AC. This study can provide a theoretical reference for explaining the tire−pavement interactions and improving vehicle brake system performance.

Elevated Temperature Baseplate Effect on Microstructure, Mechanical Properties, and Thermal Stress Evaluation by Numerical Simulation for Austenite Stainless Steel 316L Fabricated by Directed Energy Deposition

In the present study, the effect of material deposition at the elevated temperature baseplate on the microstructure and mechanical properties was investigated and correlated to the unique thermal history by using numerical simulation. Numerical results agreed well with the experimental results of microstructure and mechanical properties. Numerical results revealed a significant decrease in temperature gradient and a 40% decrease in thermal stress due to material deposition on the elevated temperature baseplate. The reduced thermal stress and temperature gradient resulted in coarser grain features, which in turn led to a decrease in hardness and tensile strength, especially for the bottom region near the baseplate. Meanwhile, no significant effect could be found for ductility. In addition, an elevated temperature baseplate promoted less heterogeneity in hardness and tensile properties along the building direction. The current work demonstrates a collective and direct understanding of the baseplate preheating effect on thermal stress, microstructure and mechanical properties and their correlations, which is believed beneficial for the better utilization of baseplate preheating positive effects.

Conventional Meso-Scale and Time-Efficient Sub-Track-Scale Thermomechanical Model for Directed Energy Deposition

Thermally-induced distortion and residual stresses in parts fabricated by the additive manufacturing (AM) process can lead to part rejection and failure. Still, the understanding of thermo-mechanical behavior induced due to the process physics in AM process is a complex task that depends upon process and material parameters. In this work, a 3D thermo-elasto-plastic model is proposed to predict the thermo-mechanical behavior (thermal and distortion field) in the laser-directed energy deposition (LDED) process using the finite element method (FEM). The predicted thermo-mechanical responses are compared to stainless steel 316L (SS 316L) deposition, with single and double bead 42-layer wall samples subject to different inter-layer dwell times, which govern the thermal response of deposited parts in LDED. In this work, the inter-layer dwell times used in experiments vary from 0 to 10 s. Based on past research into the LDED process, it is assumed that fusion and thermal cycle-induced annealing leads to stress relaxation in the material, and is accounted for in the model by instantaneously removing stresses beyond an inversely calibrated relaxation temperature. The model predicts that, for SS 316L, an increase in dwell time leads to a decrease in in situ and post-process distortion values. Moreover, increasing the number of beads leads to an increase in in situ and post-process distortion values. The calibrated numerical model's predictions are accurate when compared with in situ and post-process experimental measurements. Finally, an elongated ellipsoid heat source model is proposed to speed up the simulation.

Concentration, Propagation and Dilution of Toxic Gases in Underground Excavations under Different Ventilation Modes

The drill-and-blast method is widely used for the excavation of hard rock tunnels. Toxic gases such as carbon monoxide and nitrogen oxides are released immediately after blasting by the detonation of explosives. To provide a safe working environment, the concentration of noxious gases must be reduced below the threshold limit value according to health and safety regulations. In this paper, one-dimensional mathematical models and three-dimensional CFD numerical simulations were conducted to analyze the concentration, propagation and dilution of the blasting fumes under different operating conditions. Forced, exhaust and mixed ventilation modes were compared to determine the safe re-entry times after blasting in a 200 m-long tunnel excavated using the top-heading-and-benching method. Based on the numerical simulations, carbon monoxide was the most critical gas, as it required a longer ventilation time to reduce its concentration below the threshold limit value. The safe re-entry time reached 480 s under the typical forced ventilation mode, but was reduced to 155 s when a mixed ventilation system was used after blasting, reducing the operating costs. The reduction of the re-entry time represents a significant improvement in the excavation cycle. In addition, the results obtained show that 1D models can be used to preliminary analyze the migration of toxic gases. However, to reliably determine the safe re-entry times, 3D numerical models should be developed. Finally, to verify the accuracy of the CFD results, field measurements were carried out in a railway tunnel using gas sensors. In general, good agreements were obtained between the 3D numerical simulations and the measured values.

An Evaluation of 3D-Printed Materials' Structural Properties Using Active Infrared Thermography and Deep Neural Networks Trained on the Numerical Data

This article describes an approach to evaluating the structural properties of samples manufactured through 3D printing via active infrared thermography. The mentioned technique was used to test the PETG sample, using halogen lamps as an excitation source. First, a simplified, general numerical model of the phenomenon was prepared; then, the obtained data were used in a process of the deep neural network training. Finally, the network trained in this manner was used for the material evaluation on the basis of the original experimental data. The described methodology allows for the automated assessment of the structural state of 3D-printed materials. The usage of a generalized model is an innovative method that allows for greater product assessment flexibility.