Dynamic Numerical Simulation Research Articles

Simulation of particle deposition on solar photovoltaic panels based on a new critical capture velocity criterion

Solar energy, as a clean energy source, is becoming increasingly important in the global energy mix. However, particle deposition on the surface of photovoltaic (PV) panels can significantly reduce their power generation efficiency. In this study, the collision-deposition behaviour between silica particles and the surface of PV modules is investigated. The impact process of 13 μm silica particles on the glass surface was recorded by using a high-speed digital camera at various incident velocities and angles. A particle dynamics model was developed to predict the critical capture velocity of particles at different incident angles. It was observed that the critical capture velocity of the particles decreases as the angle of incidence increases. Subsequently, a correlation equation was established between the incident angle and the critical capture velocity, serving as the deposition criterion. Computational Fluid Dynamics (CFD) numerical simulation was employed to simulate particle deposition on PV surfaces under different wind speeds and installation tilting angles. The simulation results demonstrate that the mass of 13 μm silica particles deposited on the surface of PV panels decreases with increasing wind speed. Moreover, under identical inlet wind speeds, the particle deposition mass exhibits an initial decrease followed by a subsequent increase as the installation tilt angle of the PV panel increases. The distribution pattern of particle deposition on PV panel surfaces is diverse; however, predominantly concentrated at the mid-bottom region.

Investigation on the formation mechanism of collapse defects, microstructural evolution and mechanical properties in full penetration laser welding of thick-section steel

In this paper, we described full penetration laser welding (FPLW), a process capable of achieving single-pass welding of thick-section steel. This work summarized a welding process window encompassing various typical formations, including those well-formed and collapse defects. Furthermore, collapse defects were categorized as “collapse formation with backside humping” and “collapse formation with underfill”. High-speed camera (HSC) images reveal that the formation mechanisms of these two defects are attributed to the accumulation of liquid metal at the rear edge of the bottom molten pool and the ejection of a large amount of liquid metal from the weld in a spatter form, respectively. A computational fluid dynamics (CFD) numerical simulation model well-aligned with the experimental results has been established. The analysis of the molten pool flow behavior and keyhole evolution revealed that during well-formed welding, the keyhole status in the welding process exhibits a periodic closure-collapse mode. This periodic mode can instantaneously release the pressure at the bottom of the keyhole, inhibiting the excessive accumulation of liquid metal at the bottom of the weld and thus inhibiting the formation of collapse defects. The cooling rate in the upper part of the weld exceeds that in the middle and lower parts, primarily resulting in the formation of fine acicular ferrite and martensite. The increased tensile strength exhibited by the welded joint relative to the base material predominantly emanates from the mechanisms of grain boundary reinforcement and dislocation strengthening, precisely caused by the refined microstructure.

Numerical simulation of gas dynamics of the inflow into a channel located behind a conical or plane shock wave

Numerical simulation of gas dynamics of the inflow into a channel located behind a conical or plane shock wave

Experimental and numerical study to optimize building integrated photovoltaic (BIPV) roof structure

The study addresses the issue of heat dissipation in Building Integrated Photovoltaic (BIPV) roofs by proposing a novel PV module and establishing its testbed in Beijing, China, which is extended to more structures by the validated numerical model. The validation reveals MBD (mean bias difference）values of 10.24% for electricity production and 9.32% for temperature, accompanied by corresponding RMSD (root mean square difference) values of 23.97% and 10.46%, respectively. Taguchi method and Computational Fluid Dynamics (CFD) numerical simulation were employed to design and simulate orthogonal schemes, optimizing BIPV roof structure for the first time by considering three factors: air gap thickness, PV panel spacing and heating power. The results indicate that the air gap thickness impacts most on PV panel temperature, followed by meteorological parameters and PV panel spacing. The optimal BIPV roof structure, featured an air gap of 68 mm and a PV panel spacing of 30 mm, exhibits a 25.35% reduction in PV panel temperature, an 8.78% increase in signal-to-noise (S/N) ratio and a 2.49% growth in electricity production compared to the performance of the intermediate level scheme.

Study on the Influence of Deep Soil Liquefaction on the Seismic Response of Subway Stations

Subway systems are a crucial component of urban public transportation, especially in terms of safety during seismic events. Soil liquefaction triggered by earthquakes is one of the key factors that can lead to underground structural damage. This study investigates the impact of deep soil liquefaction on the response of subway station structures during seismic activity, aiming to provide evidence and suggestions for earthquake-resistant measures in underground constructions. The advanced finite element software PLAXIS was utilized for dynamic numerical simulations. Non-linear dynamic analysis methods were employed to construct models of subway stations and the surrounding soil layers, including soil–structure interactions. The UBC3D-PLM liquefaction constitutive model was applied to describe the liquefaction behavior of soil layers, while the HS constitutive model was used to depict the dynamic characteristics of non-liquefied soil layers. The study examined the influence of deep soil liquefaction on the dynamic response of subway station structures under different seismic waves. The findings indicate that deep soil liquefaction significantly increases the vertical displacement and acceleration responses of subway stations compared to non-liquefied conditions. The liquefaction behavior of deep soil layers leads to increased horizontal effective stress on both sides of the structure, thereby increasing the horizontal deformation of the structure and posing a potential threat to the safety and functionality of subway stations. This research employed detailed numerical simulation methods, incorporating the non-linear characteristics of deep soil layer liquefaction, providing an analytical framework based on regulatory standards for evaluating the impact of deep soil liquefaction on the seismic responses of subway stations. Compared to traditional studies, this paper significantly enhances simulation precision and practical applicability. Results from this research indicate that deep soil layer liquefaction poses a non-negligible risk to the structural safety of subway stations during earthquakes. Therefore, the issue of deep soil liquefaction should receive increased attention in engineering design and construction, with effective prevention and mitigation measures being implemented.

Evaluating hydraulic dissipation in a reversible mixed-flow pump for micro-pumped hydro storage based on entropy production theory

The promotion of renewables improves energy security but degrades the electricity quality and grid stability. This paradox can be resolved by enlarging the capacity of pumped hydro storage. Although high-head hydropower resources have been exhausted, untapped low-head resources are ideal alternatives and are utilized for micro-pumped hydro storage. Instead of the miniaturized pump-turbine or pump as turbine, we apply a low-head reversible mixed-flow pump (RMFP), which can efficiently fulfil both storage and generation requirements in microgrids without complicated guide vanes. To enhance the efficiency of the RMFP in two reverse directions simultaneously, high-hydraulic dissipation zones associated with deformed flows and unreasonable structural parameters must be investigated. This study introduces entropy production theory into computational fluid dynamics numerical simulation to evaluate hydraulic dissipation under the rated discharge in both pump and turbine modes of the RMFP. Refined mesh and validation from the energy characteristics experiment demonstrate the accuracy of the numerical scheme. The head loss assessed by the entropy Production method is compared with the predictions of the pressure drop method, and its precision exceeds 7%. In pump mode, entropy production in the volute constitutes the largest proportion (43%), followed by that in the runner (35%). However, in turbine mode, the largest proportion of the total entropy is generated in the runner (45%), followed by that in the draft tube (33%). The dissipation in the runner is bound up with the secondary flows over the blades. Backflow occurs on the suction surface near the leading edge; however, across different spans, flow separation occurs on the suction surface near the trailing edge. The requests for blade modification are opposite in pump and turbine modes. Significantly, hydraulic dissipation of the draft tube concentrates within its 0.5D2 segment upstream. The dissipation there is bound up with the transport of two types of vortex ropes. The originality lies in evaluating hydraulic dissipation and explaining the high-dissipation zone with the specific unsteady secondary flows in pump and turbine modes of an RMFP.

Dynamical chaos in the integrable Toda chain induced by time discretization.

We use the Toda chain model to demonstrate that numerical simulation of integrable Hamiltonian dynamics using time discretization destroys integrability and induces dynamical chaos. Specifically, we integrate this model with various symplectic integrators parametrized by the time step τ and measure the Lyapunov time TΛ (inverse of the largest Lyapunov exponent Λ). A key observation is that TΛ is finite whenever τ is finite but diverges when τ→0. We compare the Toda chain results with the nonintegrable Fermi-Pasta-Ulam-Tsingou chain dynamics. In addition, we observe a breakdown of the simulations at times TB≫TΛ due to certain positions and momenta becoming extremely large ("Not a Number"). This phenomenon originates from the periodic driving introduced by symplectic integrators and we also identify the concrete mechanism of the breakdown in the case of the Toda chain.

Microclimate Investigation in a Conference Room with Thermal Stratification: An Investigation of Different Air Conditioning Systems

This paper investigates the microclimate in a conference room with thermal stratification, taking as a case study the chapel of Villa San Saverio, now the seat of the “Scuola Superiore” of the University of Catania (Italy). Surveys of the former chapel were conducted to monitor air temperature and relative humidity. Subsequently, the investigation relied on numerical simulations of a simplified computational fluid dynamics (CFD) model built with the DesignBuilder v7.0 software and validated by comparison with measured values. Simulations were then carried out considering three different scenarios: the current state without any HVAC system and two possible HVAC system configurations providing both air conditioning and ventilation. The results show that, from a comfort perspective, a lightweight radiant floor heating system, assisted by an appropriate ventilation system for air renewal placed at the floor level near the occupants, is preferable to floor-level fan coils and high ventilation channels. Furthermore, this was also confirmed by a preliminary energy analysis of the two HVAC options, where the ventilation effectiveness of the winter period, the temperature of the water the emitters are fed, the consequent COP value of the heat pump, and the electricity consumption were taken into consideration.

Influence of runner blade number on hydraulic performance and flow control in draft tube of Francis turbine

As a widely used core component of hydropower generation, the stable and safe operation of the Francis turbine plays a very important role in engineering applications. Therefore, the effects of different rotor blade numbers on the hydraulic performance and flow control of the draft tube of the Francis turbine are of great research value. This paper focuses on the rotor blade number 13, 14, 15, 16, and 17 in the prohibited and stabilized operating zones, each of which is taken as a working condition point. (Condition αref - HM and Condition 2 αref - HM). Computational Fluid Dynamics (CFD) numerical simulations were performed and analyzed. In terms of hydraulic performance, the increase in the number of blades increases the efficiency and power of the turbine under Condition αref - HM by 2.26% and 0.0021 MW respectively. The efficiency and power of the turbine under Condition αref - HM is also elevated, with the efficiency reaching a maximum of 91.77% at blade number 15, and the power of the same turbine reaching a maximum of 0.182 MW at runner blade number 13. From the 3D flow analysis of the turbine, the increase in the number of blades does not significantly change the flow state of the turbine’s draft tube. By defining the turbulence energy E, the fast Fourier transform of the turbulence energy signals on the two monitoring surfaces of the draft tube is used to obtain its main frequency, and the amplitude and phase at the typical main frequency are visualized, through which we can analyze the vortex band motion state on the two monitoring surfaces of the draft tube, and speculate the change of pressure pulsation of the draft tube. Such an analysis method has a guiding significance for us to improve the performance and stability of the turbine, which can be achieved by increasing or decreasing the number of blades.

Vibrational spectrum of Granular packings with random matrices

The vibrational spectrum of granular packings can be used as a signature of the jamming transition, with the density of states at zero frequency becoming nonzero at the transition. It has been proposed previously that the vibrational spectrum of granular packings can be approximately obtained from random matrix theory. Here, we show that the autocorrelation function of the density of states shows good agreement between dynamical numerical simulations of frictionless bead packs near the jamming point and the analytic predictions of the Laguerre orthogonal ensemble of random matrices; there is clear disagreement with the Gaussian orthogonal ensemble, establishing that the Laguerre ensemble correctly reproduces the universal statistical properties of jammed granular matter and excluding the Gaussian orthogonal ensemble. We also present a random lattice model which is a physically motivated variant of the random matrix ensemble. Numerical calculations reveal that this model reproduces the known features of the vibrational density of states of frictionless granular matter, while also retaining the correlation structure seen in the Laguerre random matrix theory.Graphic abstract

A 5-D memristive hyperchaotic system with extreme multistability and its application in image encryption

By coupling memristors to nonlinear circuits, more complex dynamical behaviors can be induced. However, to date, there has been insufficient attention given to high-dimensional chaotic systems based on memristors. In this paper, a magnetic-controlled memristor is combined with a three-dimensional chaotic system, resulting in a five-dimensional memristive chaotic system. Through dynamic analysis and numerical simulations, the chaotic nature of the system is elucidated based on fundamental system behaviors, including Lyapunov dimension, dissipativity, stability of equilibrium points, 0–1 test, and Poincaré mapping. During the complex dynamical analysis of this system, unique dynamical behaviors are discovered, including intermittent chaos, transient chaos, extreme multistability, and offset-boosting. Moreover, the consistency between numerical calculations and the physical implementation of the actual system is verified through equivalent circuit design. Finally, this system is applied to image encryption, leading to the design of an efficient and secure hyper-chaotic image encryption algorithm, whose effectiveness is confirmed through several security tests.

Seismic response of deep circular tunnels subjected to S-waves: Axial bending

Ovaling deformation of circular tunnels has received great interest from the tunneling community because this mode of seismic-induced deformation is considered the most critical. However, there is growing evidence that other deformation modes can also be important and thus need to be considered in design. This study presents a new analytical solution to estimate axial bending (snaking), a mode of deformation caused by S-waves impinging on a tunnel parallel to the tunnel axis. The solution is developed using the soil-structure interaction approach with the assumption that the interface between the ground and the tunnel lining is frictionless (full-slip). Full dynamic numerical simulations are conducted to verify the new full-slip solution, together with the existing no-slip solution. Effects of dynamic amplification are also explored for both full-slip and no-slip interface conditions by changing the wavelength (or frequency) of the seismic input motions.

Dynamic Numerical Simulation of Curved Surface Coating Trajectory Based on STL Slicing Algorithm

The thickness of the coating on the surface of a workpiece is an important factor in determining the quality of spraying. However, it is challenging to estimate the distribution of film thickness accurately before the actual spraying process. This lack of estimation hinders the optimization of spraying process parameters and trajectory. To overcome this, a numerical simulation of surface spray coating thickness was conducted to provide guidance for the actual coating process. The research consists of three main parts. Firstly, the spray trajectory of the spray gun is determined using the proposed Stereo Lithography (STL) model slicing algorithm. Secondly, a two-phase flow spray model and collision adhesion model are established to construct the spray film model. The surface mesh is determined, and the spraying process parameters are set. Finally, numerical simulation is conducted to analyze the dynamic spraying trajectory and the distribution of coating thickness. The results show that the coating thickness distribution on an arc surface is thicker in the middle and thinner on the edges. The distribution is symmetric with respect to both the transverse and longitudinal directions of the arc surface. The coating thickness distribution at both ends is not as uniform as in the middle section. The concave part of the free surface has the largest coating thickness, while the coating thickness distribution on the convex part is not as uniform as on the relatively flat part. This method of simulating the coating thickness distribution on complex surfaces provides a solid foundation for further optimization of spraying process parameters and trajectory, ultimately improving the qualification rate of workpiece spraying processing.

Propulsion curve analysis and optimisation of biomimetic pectoral fin with three degree of freedom based on multi-layer perception

This paper aims to improve the swimming efficiency of biomimetic robotic fish by optimising its propulsion curve with three degree of freedom (DOF), in which the key is to minimise the resistance during swing backstroke of the pectoral fin. For this purpose, a reference point on the pectoral fin is first selected, and then an 8-shaped propulsion curve is constructed to coordinate rowing and flapping motions. Meanwhile, a propulsion curve of feathering motion is also given correspondingly. On this basis, the overall curvature and curve offset are proposed as measurement parameters for the deformation of the propulsion curve and the proportion of abdominal/dorsal circles in the propulsion curve, respectively. As the motion law of feathering is given in a sinusoidal pattern, a data set of hydrodynamic parameters is obtained for pectoral fin propulsion under different motion parameter conditions by computational fluid dynamics(CFD) numerical simulation. Once again, a propulsion curve optimisation model is established based on multi-layer perception(MLP), which uses the hydrodynamic parameters data set as input. The results show that the resistance of swing backstroke reached the minimum when overall curvature is π/10 and curve offset is 0.4, which is enhanced by 30.13% and 33.33% at thrust and lift direction compared to the existing design, respectively. Furthermore, the flow field structure and hydrodynamic parameters around the robotic fish are calculated by CFD. The results show that two sets of vortex rings are generated around the fish body during a propulsion cycle. The vortex rings detached towards the lower and upper sides of the fish body such that the lift is minimal, the thrust direction is relatively stable, and the resistance is minimal during the swing backstroke. The experimental results indicate that it can raise propulsion efficiency during the pectoral fin propulsion process, verifying the effectiveness of this method.

Numerical simulation of long-range open quantum many-body dynamics with tree tensor networks

Numerical simulation of long-range open quantum many-body dynamics with tree tensor networks

Research on wind pressure characteristics of traditional timber buildings: a case study of the main hall of Shisi Temple

Traditional timber buildings are sensitive to wind action. Studying the wind pressure characteristics is the premise for the preventive conservation of traditional timber buildings. To investigate the computational fluid dynamics (CFD) numerical simulation method for wind pressure on traditional timber buildings, a typical traditional timber building, the main hall of Shisi Temple, is chosen as a case to carry out the study. A comparative analysis is conducted to examine the effects of curve simplification of the roof slope, as well as the Dougong (bracket sets) and roof tile components, on the numerical simulation results of wind pressure on the building surface. Additionally, simplification schemes of geometric modeling are provided for the efficient and accurate simulation. The results indicate that moderate simplification of the roof curve has a relatively minor impact on the overall calculation of wind pressure, and the difference between the drag coefficients of the simplified model and the accurate model is no more than 3%. However, excessive simplification can lead to distorted simulation results, and a three-segment curve simplification method is recommended for roof cornices. The influence of Dougong on the wind pressure calculation results is negligible (within 5%), whereas roof tiles significantly reduce the drag coefficient, with an impact of over 30% at various wind directions. The impact of roof tiles on wind pressure distribution in traditional timber buildings lies in their alteration of the building aerodynamic shape rather than an increase in roof thickness. The findings can provide a basis for assessing the wind resistance of traditional timber buildings and helpful insights for improving the efficiency of wind pressure analyses of traditional timber structures.

Computational Fluid Dynamics Simulation Study on Aerodynamic Characteristics under Unfavorable Conditions during Flight Phase in Ski Jumping

Objective: The stability of the flight phase in ski jumping is crucial for athletes’ performance and safety. This study aims to investigate the influence of unfavorable conditions on aerodynamic characteristics and flight stability through computational fluid dynamics (CFD) numerical simulations. Methods: The ski jumper and the skis are considered a multi-body system. A detailed three-dimensional (3D) model of this multi-body system under a commonly observed posture during flight is established. The Partially Averaged Navier–Stokes (PANS) turbulence model is employed, and CFD simulations are conducted to predict the aerodynamic characteristics of the multi-body system under lateral environmental wind and asymmetric postures during the flight phase. The conditions of asymmetric postures include yaw rotation and roll rotation. Results: (1) Lateral environmental wind generated a yaw force, yaw moment, and roll moment, which influenced the lift, drag, and pitch moment of the athlete. These forces and moments were relatively small at lower wind speeds (less than 3 m/s) and became more significant at higher wind speeds (greater than 4.5 m/s). (2) Under the influence of yaw rotation or roll rotation, the multi-body system exhibited a noticeable yaw force, yaw moment, and roll moment, all showing a monotonic increasing trend. Moreover, they had a significant impact on the lift, drag, and pitch moment of the multi-body system. Conclusion: (1) The influence of unfavorable conditions was complex, resulting in a significant yaw force, yaw moment, and roll moment on the multi-body system. The adverse effects of roll rotation were generally greater than those of yaw rotation. (2) The multi-body system exhibited self-stabilizing tendencies in yaw and roll. This phenomenon can provide a solution to maintain flight stability by employing appropriate yaw or (and) roll rotation angles, effectively compensating for or even eliminating the adverse effects of lateral environmental wind. (3) Understanding the mechanisms of how unfavorable conditions affect aerodynamic characteristics and stability during flight in ski jumping can provide valuable assistance for real-time prediction and decision making during competitions, as well as scientific guidance for training athletes’ stable flight control and techniques for improving their sport performance.

A comprehensive fluid–solid coupling dynamic simulation for spatiotemporal distribution of regression rate in hybrid rocket motors

The spatiotemporal distribution characteristics of the regression rate are crucial aspects of the research on Hybrid Rocket Motor (HRM). This study presents a pioneering effort in achieving a comprehensive numerical simulation of fluid dynamics and heat transfer in both the fluid and solid regions throughout the entire operation of an HRM. To accomplish this, a dynamic grid technique that incorporates fluid–solid coupling is utilized. To validate the precision of the numerical simulations, a firing test is conducted, with embedded thermocouple probes being used to measure the inner temperature of the fuel grain. The temperature variations in the solid fuel obtained from both experiment and simulations show good agreement. The maximum combustion temperature and average thrust obtained from the simulations are found to deviate from the experimental results by only 3.3% and 2.4%, respectively. Thus, it can be demonstrated that transient numerical simulations accurately capture the fluid–solid coupling characteristics and transient regression rate. The dynamic simulation results of inner flow field and solid region throughout the entire working stage reveal that the presence of vortices enhances the blending of combustion gases and improves the regression rate at both the front and rear ends of the fuel grain. In addition, oscillations of the regression rate obtained in the simulation can also be well corresponded with the corrugated surface observed in the experiment. Furthermore, the zero-dimension regression rate formula and the formula describing the axial location dependence of the regression rate are fitted from the simulation results, with the corresponding coefficients of determination (R2) of 0.9765 and 0.9298, respectively. This research serves as a reference for predicting the performance of HRM with gas oxygen and polyethylene, and presents a credible way for investigating the spatiotemporal distribution of the regression rate.

Numerical Simulation Research on Aerodynamic Characteristics during Take-Off Phase in Ski Jumping

In view of the inability to directly and accurately obtain an athlete’s aerodynamic force during the take-off phase through the wind tunnel test, the athlete’s aerodynamic force and surrounding flow field form under different take-off postures are obtained through numerical simulation research, and the effects of different take-off modes on the aerodynamic characteristics during take-off in ski jumping are discussed. The multi-body system composed of the athlete and skis was selected as the research object. By using a partially averaged Navier–Stokes (PANS) turbulence model and a 3D numerical simulation of computational fluid dynamics (CFD), the aerodynamic characteristics of the athlete under different take-off postures were predicted. The take-off modes include the knee-push-hip (KPH) mode and hip-drive-knee (HDK) mode, and the hip joint angle of the HDK mode is significantly greater than that of the KPH mode. First, the aerodynamic force ratio of the athlete’s torso and legs is obviously large. Although the aerodynamic forces of arms themselves are not obvious, they have a great impact on the overall aerodynamic characteristics of the athlete, so the posture of the arms cannot be ignored. The total drag and moment of the HDK mode are significantly higher than that of the KPH mode, and the lift-to-drag ratio of the HDK mode is significantly lower than that of the KPH mode. At first, the total lift of the HDK mode is higher than that of the KPH mode, but in the last attitude, the total lift of the HDK mode does not rise but fall, and finally, the total lift of the HDK mode is lower than that of the KPH mode. The aerodynamic characteristics change dramatically during the take-off phase, and the aerodynamic characteristics of the two take-off modes are quite different, and these changes and differences are difficult to observe during real training and at the competition site. The KPH mode has an obvious aerodynamic advantage over the HDK mode. During the take-off process, the athlete should increase the force generated by the knee joint extension and appropriately reduce the speed of the hip joint extension, control the using force order of the lower limb joints, and push the hip joint extension by the knee joint extension in order to avoid issues, such as the hip joint angle being too large, the hip joint extension angle being too fast, the center of gravity being too far back, and other problems. Studying the aerodynamic characteristics during the take-off phase provides valuable insights for athletes to achieve favorable flight postures after take-off, offering scientific guidance to improve their training strategies and enhance their competitive performance.

Seismic Response Analysis of Multi Span Hyperbolic Brick Arch Thin Shell Structure

The research objective is to investigate the multi-span hyperbolic brick arch thin-shell structure located in Changleyuan, Baoji City, Shaanxi Province. The study utilizes on-site dynamic testing and finite element numerical simulation methods to analyze the dynamic characteristics of the thin-shell structure and its seismic response under various seismic waves and directions. The simulation results indicate the following findings: Vertical earthquakes cause more severe damage to structures compared to horizontal earthquakes. Bidirectional seismic excitation has a greater impact on the structure than unidirectional and three-dimensional seismic excitation. The obtained results provide a scientific basis for the seismic protection, reinforcement, and repair of this multi-span double-curved brick arch thin-shell structure.