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.

Influence of the In-Stream Structure on Solute Transport in the Hyporheic Zone

The hyporheic zone (HZ) plays an important role in the river ecosystem, and hyporheic exchange and solute transport in the HZ are important ecological functions. However, the relationship between the design parameters of river structure and solute transport is still poorly understood. In this study, we combined flume experiments and numerical simulations to systematically evaluate how in-stream structures impact the solute transport depth (D_P), hyporheic vertical exchange flux (Q), and solute flux (Qs). The results showed that the in-stream structure had a significant influence on solute transport in the HZ and could obviously increase the intensity of hyporheic exchange and promote solute transport. Model results indicated that D_P, Q, and Qs increased with the ratio of ground height to underground height of structure (H/D) and structure number (N), while Q, D_P, and Qs increased with the structural spacing (S) to begin with; then, Q remained constant, and D_P and Qs decreased as S continued to increase. This study deepened our understanding of the influence of in-stream structural design parameters on HZ solute transport, which is helpful to provide a theoretical basis for ecological restoration projects in the river HZ.

Implications of Coastal Conditions and Sea-Level Rise on Mangrove Vulnerability: A Bio-Morphodynamic Modeling Study

Mangrove forests are valuable coastal ecosystems that have been shown to persist on muddy intertidal flats through bio-morphodynamic feedbacks. However, the role of coastal conditions on mangrove behavior remains uncertain. This study conducts numerical experiments to systematically explore the effects of tidal range, small wind waves, sediment supply and coastal slope on mangrove development under sea-level rise (SLR). Our results show that mangroves in micro-tidal conditions are more vulnerable because of the gentler coastal equilibrium slope and the limited ability to capture sediment, which leads to substantial mangrove landward displacement even under slow SLR. Macro-tidal conditions with large sediment supply promote accretion along the profile and platform formation, reducing mangrove vulnerability for slow and medium SLR, but still cause rapid mangrove retreat under fast SLR. Small wind waves promote sediment accretion, and exert an extra bed shear stress that confines the mangrove forest to higher elevations with more favorable inundation regimes, offsetting SLR impacts. These processes also have important implications for the development of new landward habitats under SLR. In particular, our experiments show that landward habitat can be created even with limited sediment supply and thus without complete infilling of the available accommodation space. Nevertheless, new accommodation space may be filled over time with sediment originating from erosion of the lower coastal profile. Consistent with field data, model simulations indicate that sediment accretion within the forest can accelerate under SLR, but the timing and magnitude of accretion depend non-linearly on coastal conditions and distance from the mangrove seaward edge.

Patient Specific Numerical Modeling for Renal Blood Monitoring Using Electrical Bio-Impedance

Knowledge of renal blood circulation is considered as an important physiological value, particularly for fast detection of acute allograft rejection as well as the management of critically ill patients with acute renal failure. The electrical impedance signal obtained from kidney with an appropriate electrode system and optimal electrode system position regarding to the kidney projection on skin surface reflects the nature of renal blood circulation and tone of renal blood vessels. This paper proposes a specific numerical modelling based on prior information from MRI-data. The numerical modelling was conducted for electrical impedance change estimation due to renal blood distribution. The proposed model takes into the account the geometrical and electrophysiological parameters of tissues around the kidney as well as the actual blood distribution within the kidney. The numerical modelling had shown that it is possible to register the electrical impedance signal caused by renal blood circulation with an electrode system commensurate with the size of kidney, which makes it possible to reduce the influence of surrounding tissues and organs. Experimental studies were obtained to prove the numerical modelling and the effectiveness of developed electrode systems based on the obtained simulation results. The obtained electrical impedance signal with the appropriate electrode system shows very good agreement with the renal blood change estimated using Doppler ultrasound. For the measured electrical impedance signal, it is possible to obtain the amplitude-time parameters, which reflect the hemodynamic characteristics of the kidneys and used in diagnostics, which is the subject of further research.

Issues with Modeling a Tunnel Communication Channel through a Plasma Sheath

We consider two of the most relevant problems that arise when modeling the properties of a tunnel radio communication channel through a plasma layer. First, we studied the case of the oblique incidence of electromagnetic waves on a layer of ionized gas for two wave polarizations. The resonator parameters that provide signal reception at a wide solid angle were found. We also took into account the unavoidable presence of a protective layer between the plasma and the resonator, as well as the conducting elements of the antenna system in the dielectric itself. This provides the first complete simulation for a tunnel communication channel. Noise immunity and communication range studies were conducted for a prospective spacecraft radio line.

Rainfall Triggered Slope Instability Analysis with Changing Climate

Slope instability triggered by rainfall is common in tropical countries like India, where the soil of hill slopes are generally in unsaturated condition, having higher strength than in the saturated condition due to matric suction. This study aims to investigate the change in pore water pressure (PWP) with rainfall infiltration into the slope and the corresponding effect of the slope stability on different hill slopes from different regions of India. By applying rainfall on the surface of slope as a flux boundary condition, a transient seepage analysis is performed using SEEP/W, and the corresponding results in the form of PWP distribution within the slope are used as input for SLOPE/W program of Geostudio 2018 R2, to check the stability of slope. This study involves the use of CMIP5 (Coupled model intercomparision project phase 5) simulations for rainfall data of different regions for the duration of 2015–2050, 2051–2075 and 2076–2100; which is completely based on non-stationary approach (i.e., the statistics of extremes changes with time) rather than the conventional methods of extracting rainfall data that are based on stationary approach. The stability of different slopes is checked for 95 percentile rainfall intensity as well as for maximum rainfall intensity for a period of 2015–2100, and the analysis is performed for 5 days for which the duration of rainfall is considered to be 24 h. It is concluded that due to climate change, slopes which may be marginally stable under historical and present rainfall events can fail under projected precipitation estimates, and this can be understood by the results of Haflong slope stability analysis as the 30° slope, which is marginally stable under maximum rainfall intensity for the period 2015–2050 having factor of safety (FOS) value 1.0013, will fail under the maximum rainfall intensity for the period 2051–2075 and 2076–2100 having FOS value 0.9935 and 0.9989, respectively. A continuous decrease in FOS with increasing rainfall intensity is observed for all the study areas under consideration.

Ensemble-Based Assimilation of Satellite All-Sky Microwave Radiances Improves Intensity and Rainfall Predictions for Hurricane Harvey (2017)

Ensemble-based data assimilation of radar observations across inner-core regions of tropical cyclones (TCs) in tandem with satellite all-sky infrared (IR) radiances across the TC domain improves TC track and intensity forecasts. This study further investigates potential enhancements in TC track, intensity, and rainfall forecasts via assimilation of all-sky microwave (MW) radiances using Hurricane Harvey (2017) as an example. Assimilating Global Precipitation Measurement constellation all-sky MW radiances in addition to GOES-16 all-sky IR radiances reduces the forecast errors in the TC track, rapid intensification (RI), and peak intensity compared to assimilating all-sky IR radiances alone, including a 24-hr increase in forecast lead-time for RI. Assimilating all-sky MW radiances also improves Harvey's hydrometeor fields, which leads to improved forecasts of rainfall after Harvey's landfall. This study indicates that avenues exist for producing more accurate forecasts for TCs using available yet underutilized data, leading to better warnings of and preparedness for TC-associated hazards in the future.

Numerical Modeling and Evaluation of PEM Used for Fuel Cell Vehicles

The present study sought to analyze a novel type of polymer membrane fuel cell to be used in vehicles. The performance of the fuel cell was evaluated by modeling the types of production-consumption heat in the anode and cathode (including half-reaction heat, activation heat, and absorption/desorption heat) and waterflood conditions. The meshing of flow channels was carried out by square cells and the governing equations were numerically discretized in the steady mode using the finite difference method followed by solving in MATLAB software. Based on the simulation results, the anodic absorption/desorption heat, anodic half-reaction heat, and cathodic activation heat are positive while the cathodic absorption/desorption heat and cathodic half-reaction heat show negative values. All heat values exhibit a decremental trend over the flow channel. Considering the effect of relative humidity, the relative humidity of the cathode showed no significant change while the anode relative humidity decreased along the flow channel. The velocity at the membrane layer was considerably lower, due to the smaller permeability coefficient of this layer compared to the gas diffusion and reactants (cathode) layers.

Biomechanics of Lumbar Spine Injury in Road Barrier Collision-Finite Element Study

Literature and field data from CIREN database have shown that lumbar spine injuries occur during car crashes. There are multiple hypotheses regarding how they occur; however, there is no biomechanical explanation for these injuries during collisions with road safety barriers (RSBs). Therefore, the objective of this study was to investigate the mechanics of vertebral fractures during car collisions with concrete RSBs. The finite element method was used for the numerical simulations. The global model of the car collision with the concrete RSB was created. The lumbar spine kinematics were extracted from the global simulation and then applied as boundary conditions to the detailed lumbar spine model. The results showed that during the collision, the occupant was elevated, and then dropped during the vehicle landing. This resulted in axial compression forces 2.6 kN with flexion bending moments 34.7 and 37.8 Nm in the L2 and L3 vertebrae. It was shown that the bending moment is the result of the longitudinal force on the eccentricity. The lumbar spine index for the L1-L5 section was 2.80, thus indicating a lumbar spine fracture. The minimum principal strain criterion of 7.4% and damage variable indicated L2 and L3 vertebrae and the inferior part of L1, as those potentially prone to fracture. This study found that lumbar spine fractures could occur as a consequence of vehicle landing during a collision with a concrete RSB mostly affecting the L1-L3 lumbar spine section. The fracture was caused by a combination of axial forces and flexion bending moments.

Numerical analysis of hemodynamic effect under different enhanced external counterpulsation (EECP) frequency for cerebrovascular disease: a simulation study

Based on the changes in phase characteristics of blood flow and pressure, enhanced external counterpulsation (EECP) reduces cardiac load and improves cerebral perfusion in patients with cerebrovascular diseases. However, increased cerebral blood flow (CBF) is associated with the rise in blood pressure and its complications. Increased EECP frequency is a valuable solution when combined with the electrical equivalent impedance characteristics of the lumped parameter model (LPM) of the human blood circulation system. Herein, to investigate the effect of different EECP frequencies on CBF perfusion, an LPM was established with cardiopulmonary circulation and eight systemic blood flow units with cerebral autoregulation module of ischemic stroke patients. Then, using differential equations, we analyzed those parameters through hemodynamic simulations in four EECP modes. With related influencing parameters remaining constant, we adjusted the pressure frequency of EECP and found that when compared to the traditional sequential EECP mode, the relative increase rate of CBF was 16.68%, 18.95%, and 21.21% from 1 to 3 Hz, respectively. This study validates the effect of improving blood prefusion with increasing EECP frequency through numerical analysis.

A numerical study on drug delivery via multiscale synergy of cellular hitchhiking onto red blood cells

Red blood cell (RBC)-hitchhiking, in which different nanocarriers (NCs) shuttle on the erythrocyte membrane and disassociate from RBCs to the first organ downstream of the intravenous injection spot, has recently been introduced as a solution to enhance target site uptake. Several experimental studies have already approved that cellular hitchhiking onto the RBC membrane can improve the delivery of a wide range of NCs in mice, pigs, and ex vivo human lungs. In these studies, the impact of NC size, NC surface chemistry, and shear rate on the delivery process and biodistribution has been widely explored. To shed light on the underlying physics in this type of drug delivery system, we present a computational platform in the context of the lattice Boltzmann method, spring connected network, and frictional immersed boundary method. The proposed algorithm simulates nanoparticle (NP) dislodgment from the RBC surface in shear flow and biomimetic microfluidic channels. The numerical simulations are performed for various NP sizes and RBC-NP adhesion strengths. In shear flow, NP detachment increases upon increasing the shear rate. RBC-RBC interaction can also significantly boost shear-induced particle detachment. Larger NPs have a higher propensity to be disconnected from the RBC surface. The results illustrate that changing the interaction between the NPs and RBCs can control the desorption process. All the findings agree with in vivo and in vitro experimental observations. We believe that the proposed setup can be exploited as a predictive tool to estimate optimum parameters in NP-bound RBCs for better targeting procedures in tissue microvasculature.

On the Dynamic Electro-Mechanical Failure Behavior of Automotive High-Voltage Busbars Using a Split Hopkinson Pressure Bar

High-voltage busbars are important electrical components in today's electric vehicle battery systems. Mechanical deformations in the event of a vehicle crash could lead to electrical busbar failure and hazardous situations that pose a threat to people and surroundings. In order to ensure a safe application of busbars, this study investigated their mechanical behavior under high strain rate loading using a split Hopkinson pressure bar. Two different types of high-voltage busbars, consisting of a polyamide 12 and a glass-fiber-reinforced (30%) polyamide 6 insulation layer, were tested. Additionally, the test setup included a 1000 V electrical short circuit measurement to link the electrical with the mechanical failure. It was found that the polyamide 12 insulated busbars' safety regarding insulation failure increases at high loading speed compared to quasi-static measurements. On the contrary, the fiber-reinforced polyamide 6 insulated busbar revealed highly brittle material behavior leading to reduced bearable loads and intrusions. Finally, the split Hopkinson pressure bar tests were simulated. Existing material models for the thermoplastics were complemented with an optimized generalized incremental stress state-dependent model (GISSMO) with strain rate dependency. A good agreement with the experimental behavior was achieved, although the absence of viscoelasticity in the underlying material models was notable.

The Influence of Welding Heat Source Inclination on the Melted Zone Shape, Deformations and Stress State of Laser Welded T-Joints

The paper concerns the numerical analysis of the influence for three different of welding heat source inclinations on the weld pool shape and mechanical properties of the resulting joint. Numerical analysis is based on the experimental tests of single-side welding of two sheets made of X5CrNi18-10 stainless steel. The joint is made using a laser welding heat source. Experimental test was performed for one heating source inclination. As a part of the work metallographic tests are performed on which the quality of obtained joints are determined. Numerical calculations are executed in Abaqus FEA. The same geometrical model is assumed as in the experiment. Material model takes into account changing with temperature thermophysical properties of austenitic steel. Modeling of the motion of heating source is performed in additional subroutine. The welding source parameters are assumed in accordance with the welding process parameters. Numerical calculations were performed for three different inclinations of the source. One inclination is consistent with experimental studies. The performed numerical calculations allowed to determine the temperature field, shape of welding pool as well as deformations and stress state in welded joint. The obtained results are compared to results of the experiment.

Mind the Exposure Gaps-Modeling Chemical Transport in Sediment Toxicity Tests

Chemical exposure in flow-through sediment toxicity tests can vary in time, between pore and overlying water, and amid free and bound states, complicating the link between toxicity and observable concentrations such as free pore (C_free,pore), free overlying (C_free,over), or the corresponding dissolved concentrations (C_diss, free + bound to dissolved organic carbon, DOC). We introduce a numerical model that describes the desorption from sediments to pore water, diffusion through pores and the sediment-water boundary, DOC-mediated transport, and mixing in and outflow from overlying water. The model explained both the experimentally measured gap between C_free,over and C_free,pore and the continuous decrease in overlying C_diss. Spatially resolved modeling suggested a steep concentration gradient present in the upper millimeter of the sediment due to slow chemical diffusion in sediment pores and fast outflux from the overlying water. In contrast to continuous decrease in overlying C_diss expected for any chemical, C_free,over of highly hydrophobic chemicals was kept relatively constant following desorption from DOC, a mechanism comparable to passive dosing. Our mechanistic analyses emphasize that exposure will depend on the chemical's hydrophobicity, the test organism habitat and uptake of bound chemicals, and the properties of sediment components, including DOC. The model can help to re-evaluate existing toxicity data, optimize experimental setups, and extrapolate laboratory toxicity data to field exposure.

Investigating Optimal Confinement Behaviour of Low-Strength Concrete through Quantitative and Analytical Approaches

Reinforced concrete is used worldwide in the construction industry. In past eras, extensive research has been conducted and has clearly shown the performance of stress–strain behaviour and ductility design for high-, standard-, and normal-strength concrete (NSC) in axial compression. Limited research has been conducted on the experimental and analytical investigation of low-strength concrete (LSC) confinement behaviour under axial compression and relative ductility. Meanwhile, analytical equations are not investigated experimentally for the confinement behaviour of LSC by transverse reinforcement. The current study experimentally investigates the concrete confinement behaviour under axial compression and relative ductility of NSC and LSC using volumetric transverse reinforcement (VTR), and comparison with several analytical models such as Mander, Kent, and Park, and Saatcioglu. In this study, a total of 44 reinforced-column specimens at a length of 18 in with a cross-section of 7 in × 7 in were used for uniaxial monotonic loading of NSC and LSC. Three columns of each set were confined with 2 in, 4 in, 6 in, and 8 in c/c lateral ties spacing. The experimental results show that the central concrete stresses are significantly affected by decreasing the spacing between the transverse steel. In the case of the LSC, the core stresses are double the central stress of NSC. However, increasing the VTR, the capacity and the ductility of NSC and LSC increases. Reducing the spacing between the ties from 8 in to 2 in center to center can affect the concrete column’s strength by 60% in LSC, but 25% in the NSC. The VTR and the spacing between the ties greatly affected the LSC compared to NSC. It was found that the relative ductility of the confined column samples was almost twice that of the unrestrained column samples. Regarding different models, the Manders model best represents the performance before the ultimate strength, whereas Kent and Park represents post-peak behaviour.

Boundary element fast multipole method for modeling electrical brain stimulation with voltage and current electrodes

Objective. To formulate, validate, and apply an alternative to the finite element method (FEM) high-resolution modeling technique for electrical brain stimulation-the boundary element fast multipole method (BEM-FMM). To include practical electrode models for both surface and embedded electrodes.Approach. Integral equations of the boundary element method in terms of surface charge density are combined with a general-purpose fast multipole method and are expanded for voltage, shunt, current, and floating electrodes. The solution of coupled and properly weighted/preconditioned integral equations is accompanied by enforcing global conservation laws: charge conservation law and Kirchhoff's current law.Main results.A sub-percent accuracy is reported as compared to the analytical solutions and simple validation geometries. Comparison to FEM considering realistic head models resulted in relative differences of the electric field magnitude in the range of 3%-6% or less. Quantities that contain higher order spatial derivatives, such as the activating function, are determined with a higher accuracy and a faster speed as compared to the FEM. The method can be easily combined with existing head modeling pipelines such as headreco or mri2mesh.Significance.The BEM-FMM does not rely on a volumetric mesh and is therefore particularly suitable for modeling some mesoscale problems with submillimeter (and possibly finer) resolution with high accuracy at moderate computational cost. Utilizing Helmholtz reciprocity principle makes it possible to expand the method to a solution of EEG forward problems with a very large number of cortical dipoles.

Effect of Dispersion-Enhanced Sensitivity in a Two-Mode Optical Waveguide with an Asymmetric Diffraction Grating

Analysis of trends in the development of silicon photonics shows the high efficiency regarding the creation of optical sensors. The concept of bimodal sensors, which suggests moving away from the usual paradigm based only on single-mode waveguides and using the inter-mode interaction of guided optical waves in a two-mode optical waveguide, is developed in the present paper. In this case, the interaction occurs in the presence of an asymmetric periodic perturbation of the refractive index above the waveguide surface. Such a system has unique dispersion properties that lead to the implementation of collinear Bragg diffraction with the mode number transformation, in which there is an extremely high dependence of the Bragg wavelength on the change in the refractive index of the environment. This is called the "effect of dispersion-enhanced sensitivity". In this paper, it is shown by numerical calculation methods that the effect can be used to create optical sensors with the homogeneous sensitivity higher than 3000 nm/RIU, which is many times better than that of sensors in single-mode waveguide structures.

Long-Term Fertilization Alters Nitrous Oxide Cycling Dynamics in Salt Marsh Sediments

Salt marsh sediments are known hotspots for nitrogen cycling, including the production and consumption of nitrous oxide (N₂O), a potent greenhouse gas and ozone-depleting agent. Coastal eutrophication, particularly elevated nitrogen loading from the application of fertilizers, is accelerating nitrogen cycling processes in salt marsh sediments. Here, we examine the impact of long-term fertilization on nitrogen cycling processes with a focus on N₂O dynamics in a New England salt marsh. By combining ¹⁵N-tracer experiments with numerical modeling, we found that both nitrification and denitrification contribute to net N₂O production in fertilized sediments. Long-term fertilization increased the relative importance of nitrification to N₂O production, likely a result of increased oxygen penetration from nutrient-induced increases in marsh elevation. Substrate utilization rates of key nitrogen cycling processes revealed links between functions and the corresponding microbial communities. Higher specific substrate utilization rates leading to N₂O production from nitrification in fertilized sediments indicate a shift in the community composition of ammonia oxidizers, whereas the lack of change in specific substrate utilization of N₂O production from denitrification under long-term fertilization suggests resilience of the denitrifying communities. Both are consistent with previous studies on the functional gene community composition in these experimental plots.

A Computational Study on the Role of Parameters for Identification of Thyroid Nodules by Infrared Images (and Comparison with Real Data)

According to experts and medical literature, healthy thyroids and thyroids containing benign nodules tend to be less inflamed and less active than those with malignant nodules. It seems to be a consensus that malignant nodules have more blood veins and more blood circulation. This may be related to the maintenance of the nodule's heat at a higher level compared with neighboring tissues. If the internal heat modifies the skin radiation, then it could be detected by infrared sensors. The goal of this work is the investigation of the factors that allow this detection, and the possible relation with any pattern referent to nodule malignancy. We aim to consider a wide range of factors, so a great number of numerical simulations of the heat transfer in the region under analysis, based on the Finite Element method, are performed to study the influence of each nodule and patient characteristics on the infrared sensor acquisition. To do so, the protocol for infrared thyroid examination used in our university's hospital is simulated in the numerical study. This protocol presents two phases. In the first one, the body under observation is in steady state. In the second one, it is submitted to thermal stress (transient state). Both are simulated in order to verify if it is possible (by infrared sensors) to identify different behavior referent to malignant nodules. Moreover, when the simulation indicates possible important aspects, patients with and without similar characteristics are examined to confirm such influences. The results show that the tissues between skin and thyroid, as well as the nodule size, have an influence on superficial temperatures. Other thermal parameters of thyroid nodules show little influence on surface infrared emissions, for instance, those related to the vascularization of the nodule. All details of the physical parameters used in the simulations, characteristics of the real nodules and thermal examinations are publicly available, allowing these simulations to be compared with other types of heat transfer solutions and infrared examination protocols. Among the main contributions of this work, we highlight the simulation of the possible range of parameters, and definition of the simulation approach for mapping the used infrared protocol, promoting the investigation of a possible relation between the heat transfer process and the data obtained by infrared acquisitions.

Vein wall shrinkage induced by thermal coagulation with high-intensity-focused ultrasound: numerical modeling and in vivo experiments in sheep

BACKGROUND

Varicose veins are a common disease that may significantly affect quality of life. Different approaches are currently used in clinical practice to treat this pathology.

MATERIALS AND METHODS

In thermal therapy (radiofrequency or laser therapy), the vein is directly heated to a high temperature to induce vein wall coagulation, and the heat induces denaturation of the intramural collagen, which results macroscopically in vein shrinkage. Thermal vein shrinkage is a physical indicator of the efficiency of endovenous treatment. High-intensity focused ultrasound (HIFU) is a noninvasive technique that can thermally coagulate vein walls and induce vein shrinkage. In this study, we evaluated the vein shrinkage induced in vivo by extracorporeal HIFU ablation of sheep veins: six lateral saphenous veins (3.4mm mean diameter) were sonicated for 8 s with 3MHz continuous waves. Ultrasound imaging was performed before and immediately post-HIFU to quantify the HIFU-induced shrinkage.

RESULTS

Luminal constriction was observed in 100% (6/6) of the treated veins. The immediate findings showed a mean diameter constriction of 53%. The experimental HIFU-induced shrinkage data were used to validate a numerical model developed to predict the thermally induced vein contraction during HIFU treatment.

CONCLUSIONS

This model is based on the use of the k-wave library and published contraction rates of vessels immersed in hot water baths. The simulation results agreed well with those of in vivo experiments, showing a mean percent difference of 5%. The numerical model could thus be a valuable tool for optimizing ultrasound parameters as functions of the vein diameter, and future clinical trials are anticipated.

Determination of Optimal Flat-End Head Geometries for Pressure Vessels Based on Numerical and Experimental Approaches

The experimental and numerical analyses of the pressure vessels with different flat ends are presented and discussed in the paper. The main aim of the study is to propose the optimal flat head end geometry. The analyses are focused on the comparison of standardized geometries and with the proposed elliptical cut-out. The experimental tests with the application of strain-gauge measurements and numerical modeling of the pressure vessel are conducted. The behavior under low and high pressures and the influence of the residual welding stresses, material properties, and geometrical tolerances on the level of the plastic deformation in the flat end is discussed. It is presented that the rules given in the recent standard are not sufficient for optimal selection of the optimal geometry. It is observed that in certain geometries the deviations of the pipe thickness may lead to a significant increase of the equivalent stresses. The residual welding stresses have a significant influence on the stress and strain level—particularly in the stress relief groove (SRG). The performed study and comparison of the different geometries allow for the proposal of the optimal shape of the flat end. It appeared that the pressure vessels with SRG are the most optimal choice, particularly when elliptic shapes are in use. In some cases (i.e., pipe with wall-thickness equal to 40 mm and the flat end with circular SRG), the optimal configuration is reached for dimensions beyond the admissible by code range.

Physical and Numerical Modeling of the Slag Splashing Process

The influence of technological factors on the process of slag splashing was analyzed in the paper. The problems were solved in several stages using our own and commercial calculation programs and laboratory tests. Based on the performed calculations and simulations, factors affecting the slag splashing were determined. It was observed that the high efficiency of the process can be achieved by optimizing numerous technological parameters, e.g., flow parameters, pressure, and temperature of the nitrogen stream, height and angle of the lance position, as well as slag height into which the gas stream enters and MgO consumption. In addition, the chemical and mineralogical composition of the slag and its physicochemical parameters should be also considered. The obtained results of numerical simulations of slag splashing in the oxygen converter coincide with the results of experiments carried out using the physical model of oxygen converter. This means that the simulations well represent the real course of the slag splashing process for the studied variants.

Evolution of the Numerical Model Describing the Distribution of Non-Metallic Inclusions in the Tundish

Continuous casting is one of the steel production stages, during which the improvement in the metallurgical purity of steel can be additionally affected by removing nonmetallic inclusions (NMIs). This can be achieved by means of various types of flow controllers, installed in the working space of the tundish. The change in the steel flow structure, caused by those flow controllers, should lead to an intensification of NMIs removal from the liquid metal to the slag. Therefore, it is crucial to understand the behavior of nonmetallic inclusions during the flow of liquid steel through the tundish, and particularly during their distribution. The presented paper reports the results of the modeling studies of NMI distribution in liquid steel, flowing through the tundish. CFD modeling methods—using different models and computation variants—were employed in the study. The obtained CFD results were compared with the results of laboratory tests (using a tundish water model). The results of the performed investigations allow us to compare both methods of modeling; the investigated phenomena were microparticle distribution and mass microparticle concentration in the model fluid. The validation of the CFD results verified the analyzed computation variants. The aim of the research was to determine which numerical model is the best for describing the studied phenomenon. This will be used as the first phase of a larger research program which will provide for a comprehensive study of the distribution of NMIs flowing through tundish steel.

An Illustrative Analysis of Atypical Gas Production Profiles Obtained from In Vitro Digestibility Studies Using Fecal Inoculum

Gas production profiles typically show a monotonically increasing monophasic pattern. However, atypical gas production profiles exist whereby at least two consecutive phases of gas production or additional extraneous features that distort the typical profile are present. Such profiles are more likely to occur with the use of a fecal inoculum and are much less well described. The presence of multiple phases or non-descript extraneous features makes it difficult to apply directly recommended modeling approaches such as standard response functions or classical growth functions. To overcome such difficulties, extensions of the Mitscherlich equation and a numerical modeling option also based on the Mitscherlich are explored. The numerical modeling option uses an estimate of relative rate obtained from the smoothed data profile and an estimate of maximum gas produced together with any lag time information drawn from the raw data to construct a simple Mitscherlich equation. In summary, this article illustrates the analysis of atypical gas production profiles obtained using a fecal inoculum and explores the methodology of numerical modeling to reconstruct equivalent typical growth-like trends.