This paper presents a numerical investigation to understand the transport and deposition of sprays emitted by an impinging-jet inhaler in the human respiratory tract under different inhalation flow rates. An injection model is used for the numerical simulations considering the spreading angles of the spray in the two directions, which are measured from experiments. The model parameter is adjusted to match the mean droplet size measured in the previous experiment. A time-varying sinusoidal inhalation flow rate is utilized as airflow conditions, which is closer to the actual situation when using an inhaler. The results demonstrate that the inhalation airflow rate significantly affects the spray\'s transport behavior and deposition results in the respiratory tract. Both excessively high and low inhalation flow rates lead to an increase in deposition in the mouth-throat. A moderate inhalation flow rate reduces throat deposition while maximizing lung deposition. Higher inhalation flow rates enable faster delivery of the droplets to the lungs, whereas lower inhalation flow rates achieve a more uniform deposition over time in the lungs. The amount of deposition in different parts of the lung lobes follows a fixed order. This study provides valuable insights for optimizing the inhalation flow rate conditions of the impinging-jet inhaler for clinical applications.

Under the condition that the working face was directly covered with hard roof, the abrupt breaking of hard roof release significant amount of energy, thus prone to triggering dynamic disasters such as roadway instability or rockburst. This paper based on the engineering background of the Xieqiao Coal Mine\'s 11,618 working face, a numerical simulation method was put forward to study the dynamic response of roadway under the disturbance of hard roof breaking and proposed an evaluation index IC for roadway stability. Research indicates that the elastic energy released during the periodic weighting of the hard roof is higher than that released during the first weighting. Under the dynamic disturbance caused by hard roof breaking, the peak stresses of the roadway was slight decreased, accompanied by a significant increase in the range of stress concentration and plastic zone expansion. Roadway deformation patterns are significantly influenced by hard roof breaking, with noticeable increases in deformation on the roof and right side. During the period of hard roof breaking, the possibility of instability of the roadway increase significantly due to the disturbance caused by the dynamic load. The research results reveal the instability mechanism of roadway under the condition of hard roof, and provide a more reliable basis for evaluating the stability of roadway.

This study investigates the deformation and damage characteristics of the surrounding rock along the top return mining roadway of an isolated island working face at different stages and reveals its damage mechanism and evolution law. Utilizing a mine in Yangquan City, Shanxi Province, China, as the engineering background, this research employs FLAC 3D numerical simulation and on-site measurements. The findings suggest that the evolution of the plastic zone along the top roadway of the 15,106 island face is largely similar during both the excavation and mining periods. The plastic zones on either side of the roadway are expanding asymmetrically and gradually merging into the plastic zone of the coal pillar. In the destructive stage, the sub-gangs of the roadway are penetrated, indicating the progression into the plastic zone. The investigation points to extensive damage on the larger side of the roadway, the development of fissures, and the significant depth of damage as primary causes of roadway deformation. Moreover, the extent of the plastic zones on both sides of the roadway correlates positively with their relative distance. Continuous monitoring reveals an ongoing increase in roadway displacement, consistent with general observations in coal mining. The results provide valuable insights for optimizing support structures in similar mining environments.

OBJECTIVE: This study focuses on evaluating the disruptions in key physiological parameters during microstroke events to assess their severity.
METHODS: A mathematical model was developed to simulate the changes in cerebral tissue pO2, glucose concentration, and temperature due to blood flow interruptions. The model considers variations in baseline cerebral blood flow (CBF), capillary density, and blood oxygen/glucose levels, as well as ambient temperature changes.
RESULTS: Simulations indicate that complete blood flow obstruction still allows for limited glucose availability, supporting nonoxidative metabolism and potentially exacerbating lactate buildup and acidosis. Partial obstructions decrease tissue pO2, with minimal impact on glucose level, which can remain almost unchanged or even slightly increase. Reduced CBF, capillary density, or blood oxygen due to aging or disease enhances hypoxia risk at lower obstruction levels, with capillary density having a significant effect on stroke severity by influencing both pO2 and glucose levels. Conditions could lead to co-occurrence of hypoxia/hypoglycemia or hypoxia/hyperglycemia, each worsening outcomes. Temperature effects were minimal in deep brain regions but varied near the skull by 0.2-0.8°C depending on ambient temperature.
CONCLUSIONS: The model provides insights into the conditions driving severe stroke outcomes based on estimated levels of hypoxia, hypoglycemia, hyperglycemia, and temperature changes.

A numerical simulation method, namely, SDNMR-WEBFIT, is reported for simulating proton spin diffusion NMR based on the Levenberg-Marquardt algorithm and a pseudo-2D diffusion model. This method is used for the precise quantification of dynamics heterogeneity of the interphase within multiphase polymer systems. The numerical simulation method provides measurements of spin-lattice relaxation time (T1), proton density (ρH), lamellar thickness (d), and spin diffusion coefficient (D) for each component. The pseudo-2D diffusion model is employed to simulate the proton spin diffusion build-up/decay curves, simultaneously calculating the lateral fraction of island-like structures (x-ratio). Such approach was successfully applied to various polymer systems, such as semi-crystalline polymer (Poly(ε-caprolactone), PCL), block copolymers (Styrene-butadiene-styrene triblock copolymer, SBS), and plasticized semi-polymers (Polvinyl alcohol, PVA).

Global warming caused by the use of fossil fuels is a common concern of the world today. It is of practical importance to conduct in-depth fundamental research and optimal design for modern engine combustors through high-fidelity computational fluid dynamics (CFD), so as to achieve energy conservation and emission reduction. However, complex hydrocarbon chemistry, an indispensable component for predictive modeling, is computationally demanding. Its application in simulation-based design optimization, although desirable, is quite limited. To address this challenge, we propose a methodology for representing complex chemistry with artificial neural networks (ANNs), which are trained with a comprehensive sample dataset generated by the Latin hypercube sampling (LHS) method. With a given chemical kinetic mechanism, the thermochemical sample data is able to cover the whole accessible pressure/temperature/species space in various turbulent flames. The ANN-based model consists of two different layers: the self-organizing map (SOM) and the back-propagation neural network (BPNN). The methodology is demonstrated to represent a 30-species methane chemical mechanism. The obtained ANN model is applied to simulate both a non-premixed turbulent flame (DLR_A) and a partially premixed turbulent flame (Flame D) to validate its applicability for different flames. Results show that the ANN-based chemical kinetics can reduce the computational cost by about two orders of magnitude without loss of accuracy. The proposed methodology can successfully construct an ANN-based chemical mechanism with significant efficiency gain and a broad scope of applicability, and thus holds a great potential for complex hydrocarbon fuels.

In this paper, the nonlinear mechanical response of elastic cable structures under mechanical load is studied based on the discrete catenary theory. A cable net is discretized into multiple nodes and edges in our numerical approach, which is followed by an analytical formulation of the elastic energy and the associated Hessian matrix to realize the dynamic simulation. A fully implicit framework is proposed based on the discrete differential geometry (DDG) theory. The equilibrium configuration of a target object is derived by adding damping force into the system, known as the dynamic relaxation method. The mechanical response of a single suspended cable is investigated and compared with the analytical solution for cross-validation. A more intricate scenario is further discussed in detail, where a structure consisting of multiple slender cables is connected through joints. Utilizing the robustness and efficiency of our discrete numerical framework, a systematic parameter sweep is performed to quantify the force displacement relationships of nets with the different number of cables and different directions of fibers. Finally, an empirical scaling law is provided to account for the rigidity of elastic cable net in terms of its geometric properties, material characteristics, component numbers, and cable orientations. Our results would provide new insight in revealing the connections between flexible structures and tensegrity structures, and could motivate innovative designs in both mechanical and civil engineered equipment.

Red blood cells are destroyed when the shear stress in the blood pump exceeds a threshold, which in turn triggers hemolysis in the patient. The impeller design of centrifugal blood pumps significantly influences the hydraulic characteristics and hemolytic properties of these devices. Based on this premise, the present study employs a multiphase flow approach to numerically simulate centrifugal blood pumps, investigating the performance of pumps with varying numbers of blades and blade deflection angles. This analysis encompassed the examination of flow field characteristics, hydraulic performance, and hemolytic potential. Numerical results indicated that the concentration of red blood cells and elevated shear stresses primarily occurred at the impeller and volute tongue, which drastically increased the risk of hemolysis in these areas. It was found that increasing the number of blades within a certain range enhanced the hydraulic performance of the pump but also raised the potential for hemolysis. Moreover, augmenting the blade deflection angle could improve the hemolytic performance, particularly in pumps with a higher number of blades. The findings from this study can provide valuable insights for the structural improvement and performance enhancement of centrifugal blood pumps.
血泵中剪切应力超过阈值时红细胞会被破坏，进而引发患者出现溶血。离心式血泵叶轮结构设计对血泵的水力特性及溶血特性有着显著影响。基于此，本文采用多相流方法对离心式血泵进行数值模拟，探究了具有不同叶片数量及偏转角叶轮形式血泵的性能，分析了血泵的流场特性、水力性能以及溶血性能。数值模拟结果表明：血泵主要在叶轮及隔舌处出现了红细胞集聚现象及较大的切应力，导致此处溶血急剧增加；在一定范围内增加叶片数会提升血泵水力性能，同时也会增加溶血风险；增加叶片偏转角有助于提升血泵溶血性能，在叶片数较多时更为明显。本文研究结果可为离心式血泵的结构改进及性能改善提供参考。.

It is difficult for the existing Burgers model to accurately depict the off-axis cyclic drawing process of woven coatings. In this paper, the mechanical deformation of woven PVC (polyvinyl chloride)-coated film at different temperatures is investigated. One-dimensional (1D) and two-dimensional (2D) constitutive models were established to characterize cyclic deformation processes. The 1D model is an improved Burgers model. The effects of the time dependence of the viscosity coefficient and the ratio of elastic to viscous deformation are considered simultaneously. The accuracy of the 1D model for predicting the cyclic nonlinear deformation at different temperatures and loading rates is improved. The 2D model is a nonlinear orthotropic model using polynomials. On the basis of the single-objective genetic algorithm, the inverse algorithm is used to obtain the shear polynomial coefficients in the tension phase and the shear modulus in the unloading phase, which circumvents performing the difficult shear test. UMAT subroutines of off-axis stretching and off-axis cyclic stretching are written separately. The intelligent inverse algorithm program consists of a single-objective genetic algorithm program, a finite element parametric modelling program, and a UMAT subroutine. The simulation results are compared with the off-axis cyclic tensile test data to validate the effectiveness and accuracy of the proposed 2D model for the analysis of the woven PVC-coated films in the tension-shear coupling state.

Unconventional reservoirs, such as shale and tight formations, have become increasingly vital contributors to oil and gas production. In these reservoirs, fractures serve as crucial spaces for fluid migration and storage, making their precise assessment essential. Array acoustic logging stands out as a pivotal method for evaluating fractures. To investigate the impact of fracture width, fracture-filling conditions, and acoustic frequency on compressional and shear waves, a three-dimensional variable mesh finite difference program was employed for acoustic logging numerical simulation. Firstly, numerical models representing fractured formations with varying fracture widths and distinct fluid-filling conditions were established, and array acoustic logging numerical simulations were conducted at different frequencies. Subsequently, the waveform data were processed to extract acoustic characteristic parameters, such as velocities and amplitude attenuations of compressional and shear waves. Finally, a quantitative analysis was conducted to examine the variation patterns of characteristic parameters of refracted compressional and shear waves in relation to fracture properties. The research results indicate that amplitude attenuation information derived from borehole wave modes is particularly sensitive to the changes in fracture properties. As fracture width increased, we observed a significant amplitude attenuation in both compressional and shear waves, proportional to the logarithm of the attenuation coefficients. Furthermore, when the fracture width was constant, gas-filled fractures exhibited more prominent amplitude attenuation than water-filled fractures, with shear wave attenuation being more sensitive to the filling material. Moreover, from a quantitative perspective, the analysis revealed that the attenuation coefficients of refracted compressional and shear waves exhibited an exponential variation with gas saturation. Notably, once fracture width and filling conditions were established, the amplitudes of compressional and shear waves at the dominant frequency of 40 kHz were significantly reduced compared to those at 8 kHz, accompanied by increased attenuation. Subsequent quantitative analysis revealed that, when the product of fracture width and dominant frequency remains constant, the corresponding attenuation coefficient ratios approach 1. This indicates that the attenuation process of acoustic propagation in fractured media follows the principle of acoustic similarity. The findings of this study provide reference for further research on fracture property evaluation methods based on array acoustic logging data.