Point Velocity Research Articles

Fluid-conveying pipes in the floating frame of reference formulation

Abstract This work presents a new formulation for flexible fluid-conveying pipe elements based on the widely used floating frame of reference formulation. The elements can be used as a tool for the analysis of flexible multibody systems that contain fluid-conveying pipes, as it is well known that the movement of the fluid can influence the behavior and stability of such systems. The pipe defines a control volume through which the fluid, which is considered to be a moving mass, axially flows. The velocity of a material point of the fluid is therefore a material derivative of its position, representing the large rigid movement and small elastic deformation of the pipe along with the velocity of the fluid with respect to the pipe. The equations of motion are derived through the principle of virtual work, spatial discretization by finite element interpolation functions, and model reduction. A simplification of the consistent equations of motion is proposed, which avoids the use of inertia shape integrals and reduces the effort required to implement the developed fluid-conveying pipe elements in existing multibody software. The developed elements are validated by simulation of a straight cantilevered pipe and a curved pipe constrained on one end by a hinge. A simulation of a concrete printing system illustrates the straightforward incorporation of the elements in larger multibody systems.

Flexible Model Predictive Control for Bounded Gait Generation in Humanoid Robots.

With advancements in bipedal locomotion for humanoid robots, a critical challenge lies in generating gaits that are bounded to ensure stable operation in complex environments. Traditional Model Predictive Control (MPC) methods based on Linear Inverted Pendulum (LIP) or Cart-Table (C-T) methods are straightforward and linear but inadequate for robots with flexible joints and linkages. To overcome this limitation, we propose a Flexible MPC (FMPC) framework that incorporates joint dynamics modeling and emphasizes bounded gait control to enable humanoid robots to achieve stable motion in various conditions. The FMPC is based on an enhanced flexible C-T model as the motion model, featuring an elastic layer and an auxiliary second center of mass (CoM) to simulate joint systems. The flexible C-T model's inversion derivation allows it to be effectively transformed into the predictive equation for the FMPC, therefore enriching its flexible dynamic behavior representation. We further use the Zero Moment Point (ZMP) velocity as a control variable and integrate multiple constraints that emphasize CoM constraint, embed explicit bounded constraint, and integrate ZMP constraint, therefore enabling the control of model flexibility and enhancement of stability. Since all the above constraints are shown to be linear in the control variables, a quadratic programming (QP) problem is established that guarantees that the CoM trajectory is bounded. Lastly, simulations validate the effectiveness of the proposed method, emphasizing its capacity to generate bounded CoM/ZMP trajectories across diverse conditions, underscoring its potential to enhance gait control. In addition, the validation of the simulation of real robot motion on the robots CASBOT and Openloong, in turn, demonstrates the effectiveness and robustness of our approach.

Stream Discharge Determinations Using Slug Additions and Specific Conductance

AbstractStream discharge is often determined by wading the stream and measuring the point velocity at fixed widths and depths. However, there are conditions when wading measurements are not safe or the measurements are poor because of high turbulence, rocky streambeds, non‐standard velocity distributions, shallow or sheet flow, aquatic plants, or inaccessibility due to ice. Under these conditions, it is often preferable to determine discharge using salt slug addition and downstream measurement of salt concentration with time. A new method for determining stream discharge using specific conductance as a surrogate for salt concentrations is presented. The method adapts an approach that accurately calculates the specific conductance by utilizing ionic molal conductivities to determine the concentration of salt. The method was applied at four mountainous stream sites where a total of twenty‐nine slug‐additions were performed. The discharge determined from the new method was compared to four alternative methods including discharge from continuous injection, slug addition with discrete sample calibration, wading measurements with velocity measurement, and a stream gage. The discharge ranged from 21.5 to 778 L/s and the median difference between the new method and the traditional methods was −0.01%. Additionally, the p‐value (0.75) determined from a paired t‐test indicates that there is no significant difference between the discharge determined from the new and alternative discharge methods. The primary advantage of the new method is that it obviates the need to collect and analyze discrete samples to accurately quantify the specific conductance‐salt surrogate relationship, allowing for rapid, low‐cost determination of discharge.

Prediction of peak particle velocity using hybrid random forest approach

Blasting excavation is widely used in mining, tunneling and construction industries, but it leads to produce ground vibration which can seriously damage the urban communities. The peak particle velocity (PPV) is one of main indicators for determining the extent of ground vibration. Owing to the complexity of blasting process, there is controversy over which parameters will be considered as the inputs for empirical equations and machine learning (ML) algorithms. According to current researches, the burden has controversial impact on the blast-induced ground vibration. To judge whether the burden affects blast-induced ground vibration, the data of ground vibration considering burden have been recorded at the Wujiata coal mine. Correlation coefficient is used to analyze the relationship between variables, the correlation between the distance from blasting center to monitored point (R) and peak particle velocity (PPV) is greatest and the value of correlation coefficient is − 0.67. This study firstly summarizes the most common empirical equations, and a new empirical equation is established by dimension analysis. The new equation shows better performance of predicting PPV than most other empirical equations by regression analysis. Secondly, the machine learning is confirmed the applicability of predicting PPV. Based on the performance assessments, regression error characteristic curve and Uncertainty analysis in the first round of predicting PPV, the random forest (RF) and K-Nearest Neighbors (KNN) show better performance than other four machine learning algorithms. Then, in the second round, based on the artithmetic optimization algorithm (AOA), the optimized random forest (AOA-RF) model as the most accurate model compared with the optimized K-Nearest Neighbors (AOA-KNN) presented in the literature. Finally, the points of predicted PPV which have been informed of danger are marked based on Chinese safety regulations for blasting.

Rotationary feedback control of the cylinder wake flow using a linear dynamic model

This study presents an active feedback control of the Kármán vortex shedding flow past a circular cylinder at low Reynolds numbers. The cylinder's rotational motion functions as the control actuator, while the transverse velocities of points along the wake axis serve as the feedback signals. First, using the autoregressive with exogenous input method, a linear reduced-order model (ROM) for the unstable flow is developed to capture the input–output behavior between the cylinder's rotational displacement and the feedback signals. This model is then utilized for controller design using the proportional and linear quadratic regulator (LQR) control methods, respectively, with their effectiveness analyzed and validated through high-fidelity numerical simulations. The results show that both methods can effectively suppress the unstable vortex shedding flow, while proportional control exhibits strong sensitivity to monitoring point locations and time delays. The ROM-based model can accurately predict the stability characteristics of the control system, providing valuable guidance for selecting optimal feedback signals. Moreover, we show that by appropriately adjusting the phase angle between the control input and feedback signals via time delays, the performance of proportional control can be significantly enhanced. Lastly, based on the ROM, an output-feedback suboptimal control law is designed using the LQR method. This suboptimal feedback control transforms unstable fluid modes into stable ones, resulting in complete suppression of the unsteady vortex shedding. It is further revealed that the inherent mechanism of suboptimal flow control is to construct an optimal phase shift through the linear superposition of multiple feedback signals. Overall, model-based analysis results agree well with those obtained from direct numerical simulations, confirming the validity of the proposed ROM-based feedback control procedure.

GNSS monitoring of recent crustal movements in the region of western Bulgaria. Results of 2024th

The region of Western Bulgaria and especially Southwestern Bulgaria is one of the most active geodynamic regions in Bulgaria with the presence of a large number of active faults. The study of geodynamic processes in this area using GNSS measurements has been ongoing for more than 27 years through the measurement of a large number of geodynamic stations. This article shows the obtained results of measurements in 2024 to monitor modern crustal movements in the area. Eight points were measured. Seven of them are newly added. GNSS data obtained from periodic measurement campaigns between 1997 and 2024 have been processed and analyzed to obtain the result of the geodynamic network of point movement and strain accumulation. As a result of many years of periodic measurements, much more accurate and reliable results for point velocities are obtained. The velocity rates are between 1–2 mm/y in the north and 3–4 mm/y in the southernmost part of the area. The new results clearly confirm them, as well as the broad southward movement of southwestern Bulgaria. The results obtained from the addition of 7 new points and the expansion of the network give better and more detailed information, and can be used for detailed geodynamic and geological studies of the active tectonics in the area.

Marine Trajectory Reconstruction Method Based on Navigation State Recognition and Bi-Directional Kinematic Interpolation

The trajectory data mining and analysis of maritime targets are of great significance in furthering the construction of maritime traffic facilities, improving the ability of marine supervision and maintaining national marine security. However, due to factors such as detection means and environmental interference, a large number of trajectory data have problems such as large space-time span, uneven sampling, and poor continuity, which seriously restrict the effect of trajectory mining. Therefore, this paper proposes a method of trajectory reconstruction based on navigation state recognition and bidirectional kinematic interpolation. The method mainly includes three steps: (1) data preprocessing, (2) navigation state recognition, and (3) trajectory interpolation. The method can recognize the navigation state of the targets in different segments, and then adaptively select the interpolation method to reconstruct the trajectories, that is, linear interpolation in the straight segments and bidirectional kinematic interpolation in the turning segments. Among them, bidirectional kinematic interpolation uses the cubic Hermite function to nonlinearly fit the acceleration of the interpolation section, and then calculates the velocity and coordinates of the interpolation points by time weighting from the positive and negative directions. The proposed method is verified and analyzed on the contest dataset of “Intelligent classification and recognition of XX trajectories”. Compared with the existing methods, the reconstruction results of the proposed method are closer to the real trajectories, and it can effectively reconstruct the target trajectories with better accuracy and stability. At the same time, the effect of trajectories classification based on the Long Short-Term Memory (LSTM) model, which uses trajectories before and after reconstruction, is compared and analyzed. The results show that the model has a higher classification accuracy for reconstructed trajectory, which proves the necessity of trajectory reconstruction.

A Study on the Attenuation Patterns of Underground Blasting Vibration and Their Impact on Nearby Tunnels

The natural caving method, as a new technique in underground mining, has been promoted and applied in several countries worldwide. The destruction of the bottom rock mass structure directly impacts the structural stability of underground engineering, resulting in damage and collapse of underground tunnels. Therefore, based on the principles of explosion theory and field monitoring data, a scaled three-dimensional numerical simulation model of underground blasting was constructed using LS-DYNA19.0 software to investigate the attenuation patterns of underground blasting vibrations and their impact on nearby tunnels. The results show that the relative error range between the simulated blasting vibration velocities based on the FEM-SPH (Finite Element Method–Smoothed Particle Hydrodynamics) algorithm and the measured values is between 7.75% and 9.85%, validating the feasibility of this method. Significant fluctuations in blasting vibration velocities occur when the blast center increases to within a range of 10–20 m. As the blast center distance exceeds 25 m, the vibration velocities are increasingly influenced by the surrounding stress. Additionally, greater stress results in higher blasting vibration velocities and stress wave intensities. Fitting the blasting vibration velocities of various measurement points using the Sadovsky formula yields fitting correlation coefficients ranging between 0.92 and 0.97, enabling the prediction of on-site blasting vibration velocities based on research findings. Changes in propagation paths lead to localized fluctuations in the numerical values of stress waves. These research findings are crucial for a deeper understanding of underground blasting vibration patterns and for enhancing blasting safety.

Location Detection and Numerical Simulation of Guided Wave Defects in Steel Pipes

At present, researchers in the field of pipeline inspection focus on pipe wall defects while neglecting pipeline defects in special situations such as welds. This poses a threat to the safe operation of projects. In this paper, a multi-node fusion and modal projection algorithm of steel pipes based on guided wave technology is proposed. Through an ANSYS numerical simulation, research is conducted to achieve the identification, localization, and quantification of axial cracks on the surface of straight pipelines and internal cracks in circumferential welds. The propagation characteristics and vibration law of ultrasonic guided waves are theoretically solved by the semi-analytical finite element method in the pipeline. The model section is discretized in one-dimensional polar coordinates to obtain the dispersion curve of the steel pipe. The T(0,1) mode, which is modulated by the Hanning window, is selected to simulate the axial crack of the pipeline and the L(0,2) mode to simulate the crack in the weld, and the correctness of the dispersion curve is verified. The results show that the T(0,1) and L(0,2) modes are successfully excited, and they are sensitive to axial and circumferential cracks. The time–frequency diagram of wavelet transform and the time domain diagram of the crack signal of Hilbert transform are used to identify the echo signal. The first wave packet peak point and group velocity are used to locate the crack. The pure signal of the crack is extracted from the simulation data, and the variation law between the reflection coefficient and the circumferential and radial dimensions of the defect is calculated to evaluate the size of the defect. This provides a new and feasible method for steel pipe defect detection.

Microstructure Evolution and Mechanical Properties of Ta-Cladded Ni-Cr-Mo Low Alloyed Steel via Explosive Welding

The effects of welding variables (stand-off distance and explosive thickness) on the interfacial microstructure evolution and mechanical properties of explosively welded Ta alloy to Ni-Cr-Mo low alloyed steel have been investigated. Regardless of welding conditions (a stand-off distance of 3-5 mm and explosive thickness of 40-80 mm), the Ta/steel interface consistently exhibited a wavy configuration. This wavy interface facilitated the formation of vortex, resulting in strong interlocking. The height of the vortex increased with a larger stand-off distance at a fixed explosive thickness of 60 mm. Similarly, increasing the explosive thickness at a stand-off distance of 3 mm also resulted in a greater vortex height. The explosive weldability window, plotting the collision angle (β) against the collision point velocity (vc), was successfully established for the dissimilar Ta and steel plates. The upper limit prediction with N=0.11, as proposed by Wittman, best matched the experimental results. This guided the determination of the optimal condition, which was a standoff distance of 3 mm and an explosive thickness of 40 mm. A vortex melted zone (VMZ) was identified, which resulted from the dynamic intermixing of Ta and steel, combined with localized melting caused by high-energy collisions and heat accumulation. The VMZ surrounded by a highly deformed Ta alloy, showed the highest hardness. Near the interface on the steel side, a fine recrystallized grain structure was observed. No significant inter-diffusion was detected at the wavy Ta/steel interface. The tension-shear properties of the wavy interface, which was subjected to loading parallel to interface, showed a good balance of strength and ductility, confirming the soundness of Ta/steel interface.

Research on the trajectory tracking method of aircraft towing systems based on nonlinear control

Abstract In this paper, for the towing problem of wheeled aircraft, an automatic control method is proposed based on the nonlinear control theory, which realizes the trajectory tracking of the aircraft by controlling the velocity of the towing point of the aircraft, so as to realize the aircraft’s outbound and inbound operations. Firstly, based on the nonholonomic constraint theory, the mathematical model of the aircraft towing system is established, and the kinematic relationship between the aircraft towing point and the aircraft is obtained. Then, a controller is designed based on the nonlinear control theory, and the stability analysis is carried out using the Lyapunov method to realize the trajectory tracking of the aircraft traction system. Finally, the effectiveness and accuracy of the control method are verified by the combination of MATLAB simulation analysis and scaled-down model experiment. The results show that the proposed control method can realize the trajectory control of the aircraft towing system.

Coating flow of a liquid film with colloidal particles on a vertical fiber

A thin liquid film of colloidal suspension is considered which is spreading down on a vertical cylinder under the effect of gravity. A precursor film model is employed at the three-phase contact line to relieve the stress singularity. The curvature pressure leads to the formation of a capillary ridge at the contact-line. Bulk and surface colloids are assumed to be present in the liquid film. The interfacial pattern is governed by the film evolution equation, obtained by simplifying mass and momentum balance equations within the lubrication assumption, while rapid vertical diffusion is assumed for the advection–diffusion equation of the bulk concentration. Fluid viscosity and diffusivity are considered to be functions of the particle volume fraction. The effects of the bulk and surface colloids on the spreading film dynamics are systematically studied. The presence of bulk colloids leads to the thinning of the capillary ridge for smaller inlet bulk colloid concentrations. On the other hand, a larger inlet bulk concentration leads to the formation of a secondary advancing front behind the capillary ridge in the film profile. A reduced contact point velocity is observed with an increase in the upstream bulk concentration. Surface concentration is found to result in a hump-like structure behind the capillary ridge in the upstream direction due to solutal Marangoni stress. The Marangoni stress also hinders the surface-tension-driven spreading, resulting in a smaller contact line velocity. Reducing the Bond number results in unstable film profiles, which lead to a wave-like structure in the region with no bulk concentration gradients.

Determination of Velocities of Thread Axis Points When Interacting with a Direct Large Curvature Taking into Account the Deformation in the Contact Zone

The theoretical study of the process of the movement of the thread, which is deformed in the contact zone with the guide, from the point of view of determining the speed, is of great importance for solving a number of specific applied problems. When determining the speed, it was assumed that the start of the Lagrangian and Euler coordinates coincided. Taking into account that the crumpling of the thread in the contact zone occurs at forces much lower than the forces required to stretch the thread, the thread was considered as non-extensible. To determine the law of the distribution of the speed of the points of the thread, an independent vector, which is always connected to the main triangle, was introduced into consideration. The relationship between the vector of total curvature and the vector of absolute angular velocity is obtained. The obtained values of the speed at an arbitrary point of the axis of the thread as a function of the deformation vector, which determines the degree of twisting of the thread in the contact zone. The obtained results can be used to study various technological processes in the sewing, knitting, and textile industries, where the movement of the thread takes place along the guide surface of great curvature.

High potency of application on an open direct-contact thermal storage using humid air

To reuse low-temperature wasted heat as a thermal resource for high temperature, a direct-contact adsorption thermal storage was focused using humid air and zeolite 13X particles as the working fluid and adsorbent, respectively. Only a few previous studies have chosen the working fluid in gaseous form because it is unavailable for latent heat in generating heat sources. Recovering waste heat in humid air to generate hotter steam is unique and becomes an originality of our present work. The time required to regenerate zeolite particles and the maximum temperature of the generated steam were investigated assuming a warm-up device for a vehicle. The time required to regenerate zeolite was investigated by changing the dew point, temperature, and superficial velocity of the inlet humid air. It was mainly affected by the temperature of the inlet air. The absorbent was regenerated within 30 min when the humid air preheated to 200 °C was supplied. On the other hand, the maximum steam temperature was investigated by changing the superficial velocity and temperature of saturated inlet humid air. As one of the significant and novel finding in this work, the steam of >200 °C was obtained as a high-temperature heat source even with saturated humid air unavailable latent heat. Moreover, as theoretical knowledge, it was revealed that the maximum temperature of the heat source can be estimated by the relationship between the heat balance on the packed bed and adsorption equilibrium.

Design and kinematics analysis of a parabolic cylinder deployment mechanism with modular over-constrained triangular pyramid

Space missions require novel mechanisms that can have compact form in complex space environments. This study proposes a modular triangular pyramid parabolic cylinder deployment mechanism. The modular deployable unit contains an over-constrained triangular pyramid deployable mechanism and a symmetrical trapezoidal Bricard linkage that drives the longitudinal and transverse motions of the mechanism. A cable net is added between two symmetrical linkages to form a parabolic cylindrical reflector and is analyzed by the force density method. The Denavit-Hartenberg (D-H) coordinate method is used to analyze the kinematic characteristics. Based on the analytical solution for the angular displacement, the degrees of freedom of the mechanism are derived from screw theory. Subsequently, the angular velocity and acceleration of the joint points are obtained. Finally, kinematic models of the modular triangular pyramid parabolic cylinder deployable mechanism are established, and the accuracy of the theoretical model is verified using a numerical method. This novel parabolic cylinder deployable mechanism will have a significant influence in aerospace domain. Crucially, the work has value to design deployment mechanism for parabolic cylinder antenna.

Symmetry‐based analysis of nonlinear mixed convection in 3D EMHD nano‐Carreau fluid flow with Riga stretched surface effects and multi‐physical interactions

AbstractThis study presents a comprehensive investigation into the dynamics of an electrically magneto‐hydrodynamic (EMHD) nano‐Carreau fluid under nonlinear mixed convection. We develop a 3D steady‐state framework that incorporates various influential factors such as nonuniform heat source‐sink terms, nonlinear thermal radiation, Joule heating, and chemical reactions, along with the effects of a Riga stretched surface. Through rigorous analysis, we explore the impact of thermophoretic and Brownian motions on flow patterns and stagnation point velocities. Our study encompasses scenarios involving a Riga stretching sheet, EMHD phenomena, porous media, suction‐injection processes, and diverse slip conditions (momentum, heat, volume fractions), in conjunction with chemical reactions. By employing symmetry transformations, we transform complex partial differential equations (PDEs) into more manageable ordinary differential equations (ODEs), facilitating effective numerical solutions using the Lobatto IIIa bvp4c method in Matlab. The findings are presented through detailed graphical representations and comparative tables. Key findings include the observation that elevated Hartmann numbers contribute to reduced velocity yet enhanced temperature profiles, influenced by factors such as nonuniform heat distribution, thermal radiation, and viscous dissipation. Additionally, concentration profiles exhibit a diminishing trend with increased Lewis numbers, chemical reactions, and specific slip parameters.

Dual-stage control method for continuous hopping on the surface of small celestial bodies

Surface movement exploration techniques are an important approach to further the understanding of small celestial bodies. Hopping movement is employed to take advantage of the weak gravitational field of small celestial bodies, facilitating long-distance mobility of the rover. The rover exhibits a strong attitude-orbit coupling characteristic during the hopping movement. Different contact attitudes contribute to the challenge of determining the rebound trajectory, rendering it difficult for the rover to maintain continuous movement towards the target position along the nominal trajectory. In this paper, a dual-stage control method for continuous surface hopping movement is developed, wherein the control is applied separately during the flight stage and the collision stage. During the flight stage, applying attitude control enables the rover to establish contact with the surface in a specific attitude, ensuring accuracy in movement direction and creating favorable conditions for the control process in the collision stage. The flywheel is controlled in the collision stage to correct the state change and energy loss due to the collision, enabling the rover to attain the desired hopping state. Considering the energy consumption of the rover, the minimum hopping velocity sequence and angular velocity sequence are designed based on the hopping angle and heading angle. The intricacies and uncertainties inherent in the surface environment of small celestial bodies pose challenges for the rover during hopping and flying. In response to the issue, a finite-time adaptive sliding-mode control law is devised. The control law facilitates the rapid adjustment of the rover to the target attitude in flight, demonstrating strong robustness. In addition, by integrating the changes in contact point velocity with the impulse contact model, the accuracy of the rover's rebound state updates is enhanced. Finally, a set of simulations is performed, and simulation results indicate that the dual-stage control method can achieve a final landing position error of less than 0.2 m for a continuously hopping. The finite-time adaptive sliding-mode control law can stabilize the rover's attitude within 10 s. Conducting 500 Monte Carlo simulations of the continuously hopping rover in three small celestial body surface simulation environments shows that the position endpoint error is within 0.3 m.

Kinematics analysis of planetary roller screw mechanism with misalignment

The misalignment is unavoidable for planetary roller screw mechanism (PRSM) in practical application due to many uncertainties. The misalignment has a strong effect on the kinematic characteristics of PRSM. However, the relevant research has been rarely reported in the past. In order to clearly present the motion characteristics, a novel kinematics model of PRSM with misalignment is established in this article. The rotational tensor theory is introduced to obtain the coordinate transformation equations. The position vectors of PRSM with misalignment are acquired to calculate the position of the contact point. The motion model of PRSM with misalignment is derived by analyzing the velocity of the contact point. Compared with some published literature, the kinematic model of PRSM is verified. Moreover, a relationship between the misalignment variates, which are the misalignment angle, the misalignment azimuth angle, and the offset distance of the nut's ending point, respectively, and the kinematic characteristics of PRSM are investigated. The results indicate that those misalignment variates have a great influence on the kinematics of PRSM and lead to fluctuations in the PRSM's transmission velocity and angular velocity ratio. The fluctuation amplitude of transmission velocity is positively correlated with misalignment angle.

Oscillatory instability of magneto-convection during the melting of non-Newtonian ferromagnetic nano-enhanced phase change materials at low Rayleigh number

This work complements mechanism of the oscillatory instability of magneto-convection in the process of melting heat transfer controlled by magnetic field at low Rayleigh number (Ra), aiming to provide a theoretical basis for improving the stability of energy storage system and melting manufacturing process. A latent heat storage cavity unit exposed in uniform magnetic field is designed to reveal the magnetic-controlled heat transfer in solid–liquid phase change process. The volume average method is used to describe porous media composite nano-enhanced phase change materials in latent heat unit. Special attentions are paid to the key parameters of Magnetic number (Mn) and Ra, and coupling flow heat transfer mechanism of magnetic force and natural convection effect is numerically studied. The evolution stages of different flow patterns are divided according to the Kelvin force, vorticity, velocity, flow field and isotherm distributions of characteristic points. Main results show that the magnetic field-modulated phase change process at low Ra can be divided into three stages, of which the second stage is characterized by periodic oscillations. Increasing Mn at fixed Ra accelerates the second stage of oscillatory convection heat transfer coming earlier and promotes melting. When Mn increase to 5 × 106, the full melting time is reduced to the greatest extent, resulting in 11.6 % reduction in heat storage but 17.4 % increase in heat storage efficiency. Three mechanisms of magnetic forced convection coupled with natural convection are found during the increase of Ra: synergistic (Ra < 2 × 103), antagonistic (2 × 103≤ Ra < 1 × 104), and natural convection-dominated (Ra = 1 × 104). Finally, the mode distribution map of Ra-Mn regulation is obtained, in which the unstable transformation mechanism of steady-state, periodic oscillation and chaotic oscillation heat transfer in the melting process under the control of magnetic field is described.

Numerical Study on Permeability of Reconstructed Porous Concrete Based on Lattice Boltzmann Method

The reconstruction of the porous media model is crucial for researching the mesoscopic seepage characteristics of porous concrete. Based on a self-compiled MATLAB program, a porous concrete model was modeled by controlling four parameters (distribution probability, growth probability, probability density, and porosity) with clear physical meanings using a quartet structure generation set (QSGS) along with the lattice Boltzmann method (LBM) to investigate permeability. The rationality of the numerical model was verified through Poiseuille flow theory. The results showed that the QSGS model exhibited varied pore shapes and disordered distributions, resembling real porous concrete. Seepage velocity distribution showed higher values in larger pores, with flow rates reaching up to 0.012 lattice point velocity. The permeability–porosity relationship demonstrated high linearity (the Pearson correlation coefficient is 0.92), consistent with real porous concrete behavior. The integration of QSGS-LBM represents a novel approach, and the research results can provide new ideas and new means for subsequent research on the permeability of porous concrete or similar porous medium materials.