A comprehensive modeling framework for the thermoforming of polymer matrix woven laminate composite was developed. Two numerical indicators, the slip path length and traction magnitude, have been identified to be positively correlated to matrix smearing and wrinkling defects. The material model has been calibrated with picture-frame experimental results, and the prediction accuracy for intra-ply shear and thickness distribution was examined with measurements of the physically formed parts. Specifically, thickness prediction for most locations on the formed parts was accurate within an 11.6% error margin. However, at two points with significant intra-ply shear, the prediction errors increased to around 20%. Finally, a parametric study was conducted to determine the relationship between various process parameters and the quality of the formed part. For the trapezoidal part, orienting the laminate at 45 degrees to the mold axis reduces the likelihood of matrix smear and wrinkling defects. Although this laminate orientation yielded a greater spatial variation in part thickness, the thickness deviation is lower than that for the 0-degree orientation case. Two forming analyses were conducted with ramp rates of 25 mm/s and 80 mm/s to match the equipment\'s operational limits. It was observed that higher forming rates led to a greater likelihood of defects, as evidenced by a 15% and 10% increase in the formed part areas with longer slip paths and higher traction magnitudes, respectively. It was discovered that shallower molds benefit from faster ramp rates, while deeper molds require slower rates to manage extensive shearing, stretching and bending. Faster forming rates lead to smaller thickness increases at high intra-ply shear regions, indicating a shift from intra-ply shear to out-of-plane bending due to the visco-plastic effect of the molten laminate and can negatively impact part quality. Lastly, it was shown that a well-conceived strategy using darts could improve the part quality by reducing the magnitude of the defect indicators.

Cracks in rock and concrete have a great adverse effect on the stability of engineering structures; however, there are few studies on X-shaped fissures which widely exist in rock and concrete structures. Based on this background, three-point bending fracture tests of SCB specimens containing X-shaped fissures are carried out. The momentum equations in the SPH method are improved, and the crack propagations of SCB specimens under three-point bending are simulated. The results show that cracks grow simply along the vertical direction in the sample with no X-shaped fissures, and the existence of an X-shaped fissure changes the crack growth path and final failure modes of the SCB samples. The crack propagation simulation results are consistent with the experimental results, which verifies the rationality of the improved SPH method. The load-displacement curves mainly present three typical stages: the initial compaction stage, linear elastic deformation stage, and failure stage. The peak load decreases first then increases with an increase in eccentricity. With an increase in X-shaped fissure length and decrease in X-shaped fissure angle, the peak load decreases. The damage counts remain at 0 at the initial loading stage, corresponding to the initial compaction stage and the linear elastic deformation stage, and increase sharply at the later loading stage, corresponding to the failure stage, which is consistent with the experimental results. The influence mechanisms of X-shaped fissures on the crack propagation paths are discussed; the existence of different X-shaped fissure morphologies aggravate the tensile stress concentration at specific positions, leading to different crack propagation modes in the experiments. The research results can provide a certain reference for understanding the failure mechanisms of engineering structures containing X-shaped fissures and promote the applications of the SPH method into the simulations of cross-fissure crack propagations.

This study introduces a fractional order model to investigate the dynamics of polio disease spread, focusing on its significance, unique results, and conclusions. We emphasize the importance of understanding polio transmission dynamics and propose a novel approach using a fractional order model with an exponential decay kernel. Through rigorous analysis, including existence and stability assessment applying the Caputo Fabrizio fractional operator, we derive key insights into the disease dynamics. Our findings reveal distinct disease-free equilibrium (DFE) and endemic equilibrium (EE) points, shedding light on the disease\'s stability. Furthermore, graphical representations and numerical simulations demonstrate the behavior of the disease under various parameter values, enhancing our understanding of polio transmission dynamics. In conclusion, this study offers valuable insights into the spread of polio and contributes to the broader understanding of infectious disease dynamics.

The site-specific signal enhancement provided by parahydrogen induced polarization (PHIP) may be combined with magnetic resonance imaging (MRI) to study chemical and biomolecular processes. However, imaging of hydrogen nuclei (1H) is hampered by background signals arising from the presence of thermally polarized nuclei. Additionally, fast imaging sequences are commonly based on multiple radio-frequency pulses, where the signals resulting from PHIP oscillate due to the evolution with a J-coupling Hamiltonian. In this article, an innovative imaging scheme for single-scan MRI is presented that effectively detects hyperpolarized components while simultaneously canceling out thermal contributions. This method is based on the quenching of inherent oscillations of PHIP-originated signals due to J-couplings during the multipulse sequence and the suppression of thermal signals by spin dynamics and a tailored restructuring of the k-space. A series of numerical simulations on specific two- and three-spin systems serve to support the feasibility of the approach. Furthermore, this theoretical study demonstrates the potential of combining hyperpolarization and long-lived states (PHIP and LLS) in the selected molecules, which could be seen as a preliminary step towards the development of fast imaging techniques, for example in the field of biomolecular research.

The practical applications of underwater optical devices, such as cameras or sensors, often suffer from widespread surface biofouling. Current antifouling techniques are primarily hindered by low efficiency, poor compatibility, as well as environmental pollution issues. This paper presents a transparent electrode coating as antifouling system of underwater optics as potential substitute for alternating current electrokinetic (ACEK)-based systems. A strong-coupling model is established to predict the Joule heating induced fluid flows and the negative dielectrophoretic (nDEP) effect for mobilizing organisms or deposited sediments on optic surfaces. The performance of the proposed antifouling system is numerically evaluated through simulations of electrostatic, fluid and temperature fields as well as trajectories of submicron particles, which is then experimentally verified and found to be in good agreement. A parametric study revealed that the degree of electrodes asymmetry is the key factor affecting the flow pattern and therefore the overall performance of the system. This ACEK-based universal strategy is expected to shed light on designing high performance and non-toxic platforms toward energy-efficient surface antifouling applications of underwater optics.

Monitoring the microstructural change in cementitious materials during hydration is an essential but challenging task. Therefore, a non-invasive and sophisticated technique is warranted to understand the microscopic behaviour of the multiphase cementitious materials (where the length scale of the constituents varies from centimeters to micrometers) in different stages of hydration. Due to exothermic hydration reactions, different hydration products start to evolve with individual mechanical properties. In concrete, an interface transition zone (ITZ) appears between the aggregate surface and paste matrix, which influences the overall properties of concrete material. In the present research, 1) several wave characteristics, such as wave velocity, energy distribution, and signal phase are found out using Ultrasonic Pulse Velocity (UPV), Wavelet Packet Energy (WPE) and Hilbert Transform (HT) methods, to monitor the hydration mechanism (1d-28d) in cement-based materials with two levels of heterogeneities (cement paste and concrete, representing microscale and mesoscale, respectively). Also, the unique nonlinear behaviour is studied in the frequency domain using the promising Sideband Energy Ratio (SER) and Sideband Peak Count Index (SPC-I) methods. 2) Numerical simulations are carried out to understand the wave interaction in the developing microstructure. A discretized microstructure of cement shows microscopic details of each phase at any instant of hydration (e.g., formation stage and after complete maturity level). The experimental and numerical investigations on the characteristics of the nonlinear ultrasonic wave propagation show the impact of microstructural development of multi-scale cementitious materials during hydration.

In this paper, piezoceramic-based excitation of shear horizontal waves is investigated. A thickness-shear d15 piezoceramic transducer is modeled using the finite-element method. The major focus is on the directivity and excitability of the shear horizontal fundamental mode with respect to the maximization of excited shear and minimization of Lamb wave modes. The results show that the geometry of the transducer has more effect on the directivity than on the excitability of the analyzed actuator. Numerically simulated results are validated experimentally. The experimental results show that transducer bonding significantly affects the directivity and amplitude of the excited modes. In conclusion, when the selected actuator is used for shear excitation, the best solution is to tailor the transducer in such a way that at the resonant frequency the desired directivity is achieved.

The rate of mass transfer of lower molecular weight hydrocarbons (naphtha) from bitumen drops in mature fine tailings of oil sand tailings ponds (OSTPs) may control their bioavailability and the associated rate of GHG production. Experiments were conducted using bitumen drops spiked with o-xylene and 1-methylnaphthalene to determine the mass transfer rate of these naphtha components from bitumen drops. The results were compared to simulations using a multi-component numerical model that accounted for transport in the drop and across the oil-water interface. The results demonstrate rate-limited mass transfer, with aqueous concentrations after 60 days of dissolution that were different than those in equilibrium with the initial drop composition (less for o-xylene and greater for 1-methylnaphthalene). The simulations suggest that mole fractions were unchanged at the center of the drop, resulting in concentration gradients out to the oil-water interface. Numerical simulations conducted using different drop sizes and bitumen viscosities also suggest the potential for persistent naphtha dissolution, where the time required to deplete 80% of the o-xylene and 1-methylnaphthalene mass from an oil drop was estimated to be on the order of months to years for mm-sized drops, and years to decades for cm-sized drops assuming instantaneous biodegradation in the aqueous phase surrounding the bitumen.

Leakage accidents of buried pipelines have become increasingly common due to the prolonged service of some pipelines which have been in use for more than 150 years. Therefore, there is an urgent need for accurate prediction of pollution scope to aid in the development of emergency remediation strategies. This study investigated the distribution of a light non-aqueous phase liquid in soils containing gas and water through numerical simulations and laboratory experiments. Firstly, a three-dimensional porous medium model was established using ANSYS FLUENT, and for the first time, the distribution of gas and groundwater in soil environments was simulated in the model. Subsequently, the distribution of the three phases of diesel, gas, and water in soil was studied with different leakage velocities and it was found that the leakage velocity played a significant role in the distribution. The areas of diesel in soils at 60 min were 0.112 m2, 0.194 m2, 0.217 m2, and 0.252 m2, with corresponding volumes of 0.028 m3, 0.070 m3, 0.086 m3, and 0.106 m3, respectively, for leakage velocities of 1.3 m/s, 3.4 m/s, 4.6 m/s, and 4.9 m/s. Calculation formulas for distribution areas and volumes were also developed to aid in future prevention and control strategies under different leakage velocities. The study also compared the distribution areas and volumes of diesel in soils with and without groundwater, and it was found that distribution scopes were larger in soils containing groundwater due to capillary force. In order to validate the accuracy of the numerical simulation, laboratory experiments were conducted to study the diffusion of oil, gas, and water under different leakage velocities. The results showed good agreement between the experiments and the simulations. The research findings are of great significance for preventing soil pollution and provide a theoretical basis for developing scientifically sound soil remediation strategies.

In recent years, blending mechanistic knowledge with machine learning has had a major impact in digital healthcare. In this work, we introduce a computational pipeline to build certified digital replicas of cardiac electrophysiology in paediatric patients with congenital heart disease. We construct the patient-specific geometry by means of semi-automatic segmentation and meshing tools. We generate a dataset of electrophysiology simulations covering cell-to-organ level model parameters and using rigorous mathematical models based on differential equations. We previously proposed Branched Latent Neural Maps (BLNMs) as an accurate and efficient means to recapitulate complex physical processes in a neural network. Here, we employ BLNMs to encode the parametrized temporal dynamics of in silico 12-lead electrocardiograms (ECGs). BLNMs act as a geometry-specific surrogate model of cardiac function for fast and robust parameter estimation to match clinical ECGs in paediatric patients. Identifiability and trustworthiness of calibrated model parameters are assessed by sensitivity analysis and uncertainty quantification.