Numerical Simulation Method Research Articles

Numerical simulation of fuselage sliding on calm water for an amphibious aircraft based on global moving mesh method

Abstract Aiming at the large motion fluctuation of fuselage sliding on calm water, a numerical simulation method based on the global moving mesh method is introduced by using the computational fluid dynamics software of STAR-CCM+, the two equations SST k-w turbulence model to solve the RANS equations, the VOF method to capture free surface, the DFBI module to simulate the movement of the fuselage. We established a simulation scheme on the sliding movement of the fuselage on calm water, which effectively solves the problem of calculation divergence caused by the excessive stretching deformation of the mesh. Comparison between calculation results and model test results shows that the trends of resistance and attitude angles changing with velocity are basically the same, and the relative errors are around 10%, which indicates that the global moving mesh method is feasible in the simulation of fuselage sliding on calm water.

Simultaneous effects of capillary number, viscosity ratio, and contraction ratio on droplet dynamics in contraction microchannel

Abstract In droplet-based microfluidic systems, fluid properties, flow conditions, and channel geometry play crucial roles in the movement and manipulation of droplets. This study employs a three-dimensional numerical simulation method to investigate droplet dynamics in a contraction microchannel under the simultaneous influence of different factors, such as capillary number (Ca), viscosity ratio (λ), and contraction ratio (C). Specifically, a theoretical model is introduced for predicting the transition of droplets from trap to squeeze, and it is observed to be 〖Ca〗_I = f(C,λ). Meanwhile, droplet dynamics are investigated for the transition from squeeze to breakup for finding critical capillary number value (CaII). In addition, the critical viscosity ratio for the process from squeeze to breakup has been figured out by numerical results for a wide range of contraction ratios. Eventually, droplet dynamics including deformation, velocity ratio, and breakup during contraction microchannel are also quantitatively studied under the factors above.

Numerical simulation of influence of jet flow parameters on cone-cylinder-flare configuration

Abstract Reaction control system technology is used to realize attitude control or trajectory change. The hot jet effect should be considered to get more accurate aerodynamic characteristics as the flight vehicle develops. The simulation for hot jet flow needs more calculation resources, so it is important to research the suitable similarity parameters of jet flow. Then, we can use cold jet flow to simulate hot jet flow. In this paper, the influence of jet flow parameters is numerically studied on cone-cylinder-flare configuration. The numerical simulation method for hot jet flow is established, containing finite volume method, multi-component three-dimensional Reynolds Navier-Stokes equation with source terms, multi-block patched grid technique, and fully implicit processing method. The numerical method is verified by using experimental reference results, and the grid independence is assured. The numerical results reveal that the cold jet flowfield under the same mass flow ratio is more similar to the hot jet flowfield than the cold jet flowfield under the total temperature of 293 K. The maximum amplification factor difference between cold and hot jet flow is reduced from 14.1% to 1.2%.

Spatio-temporal neural networks for monitoring and prediction of CO[formula omitted] plume migration from measurable field data

Carbon capture and storage (CCS) technologies are crucial for mitigating greenhouse gas emissions. The success of CCS projects hinges on accurately predicting and monitoring the CO2 plume migration during and after injection. To address the computational burden of traditional numerical simulation methods, previous studies have used neural networks as proxy models to expedite the prediction of the CO2 plume migration. However, these models rely on uncertain inputs, such as the distribution of heterogeneous rock permeability and porosity maps, which can lead to erroneous predictions and limit their adoption in real-world applications. To address this issue, this study introduces a framework for reconstruction and short-term prediction of CO2 plume migration based solely on field measurements that directly provide information on the migration of the CO2 plumes, thus eliminating the dependency on uncertain geological information as model input. The framework trains a spatio-temporal neural network model using simulated data under various geologic settings to capture the plume evolution dynamics without constraining it to a specific geological scenario. Once trained, the model integrates global and local field measurements from multiple sources to reconstruct the CO2 plume and predict its spatio-temporal evolution. The effectiveness of the proposed framework is tested using two case studies: one using a synthetic dataset and another with data generated from a model of a real field in the Southern San Joaquin Basin. The results show that the proposed framework can accurately reconstruct and perform short-term predictions of the CO2 plume migration by integrating various forms of field measurements.

Prediction and Analysis of Borosilicate Glass Surface Deformation Induced by Flame Jet

To address the issues of low processing efficiency, poor forming accuracy, and internal damage in glass material processing, this study proposes the use of flame jet forming. However, the mechanism of flame jet processing requires further elucidation. This research investigates the relationship between the indentation morphology on the glass surface and the inlet velocity of the flame jet. A theoretical model was established through mathematical analysis to reflect this relationship. The model’s accuracy was validated using numerical simulation methods. By comparing experimental data with theoretical model results, surface tension was incorporated, and the model was iteratively optimized using MATLAB R2024a. The final optimized model demonstrated an absolute error range of 0.009 to 0.069 mm. This study confirms the feasibility of flame jet processing and enriches the understanding of its mechanism, providing a novel, efficient, and precise method for processing glass materials.

Research on the Mechanism of Thermal Power Enhancement in an Interior Permanent Magnet Eddy Current Heater Driven by Wind

This article uses numerical simulation methods to study the comprehensive influences of the stator structure and materials on the thermal power of an interior permanent magnet eddy current heater (IPMECH). By analyzing the air-gap magnetic flux density (MFD), stator MFD, thermal power, and torque at different rotational speeds, the mechanism of thermal power enhancement has been revealed in depth. The results indicate that the armature magnetic field (MF) generated by the eddy current produces a magnetization effect on the side of its rotation direction, but the MF in the stator will be weakened in general, and this effect becomes more significant with the increase in the rotational speed. The stator material of the IPMECH has higher permeability, which has higher thermal power and torque, and a lower proportion of high-order harmonics, which is beneficial for reducing the radial vibration of the IPMECH. A permanent magnet with high remanence can increase the thermal power and torque of the IPMECH. Reducing the length of the air gap is beneficial for improving the thermal power, but it also increases the harmonic MFD. The rotational speed is 200 rpm, the air gap is 0.1 mm and 2 mm, and the thermal power is 1.12 kW and 0.35 kW, respectively. The fundamental amplitudes of the 0.1 mm and 2 mm air-gap lengths are 0.94 T and 0.64 T, respectively, and the 3rd harmonic Bi* values are 0.24 and 0.18, respectively. At rotational speeds of 200 rpm, 800 rpm, and 1600 rpm, the δPmax values are 17 mm, 11 mm, and 8 mm, respectively. When designing a heater, the higher the rotational speed, the smaller the stator wall thickness should be.

Research on the gravity sedimentation law of watery crude oil in large oil tanks

ABSTRACT Based on field survey data, a numerical simulation method that combines an indoor settling experiment with the Mixture multiphase flow model and the Population Balance Model (PBM) was adopted in this article to study and analyze the sedimentation law of water-containing crude oil in a tank under different initial particle sizes and different continuous-phase viscosities. The results indicate that the maximum error in the comparison of the water content changes at different positions in the sedimentation cylinder after 1 h with experimental values is approximately 8.34%, confirming the precision of the developed settling model. When the initial particle size of liquid droplets in the tank increases from 8 μm to 64 μm, the maximum water content in the tank increases by 1.115%, and the magnitude of the increase in the average droplet size decreases by 90%. As the viscosity of the oil increases (from 35.26 mPa · s to 77.77 mPa · s), the maximum water content in the tank decreases by 1.331%, while the height of the stable water content region increases by 13.2%.

Research on aerodynamic performance of centrifugal compressors for hydrogen-mixed natural gas.

Natural gas-doped hydrogen transportation has been widely used in industrial engineering. The change of the physical parameters of the conveying medium affects the working performance of the centrifugal compressor. This study aims to explore the aerodynamic performance of centrifugal compressors for hydrogen mixed natural gas, including the effects of hydrogen blending ratio (HBR) and inlet temperature on the pressure ratio and isentropic efficiency of compressors. The numerical simulation method is used and the results are compared with the experimental results to determine the reliability of the numerical simulation method. The results show that the pressure ratio decreases with the increase of HBR where the inlet temperature and rotation speed are constant. The pressure ratio and efficiency of the compressor are the highest when the conveying medium is pure natural gas (HBR = 0). The maximum pressure ratio reduction is 4.8% when the HBR is 5% and the inlet temperature is increased by 20 K. In summary, when the conveying medium of the compressor changes from pure natural gas to hydrogen-mixed natural gas, the working range and efficiency of the compressor will be reduced. Therefore, it is necessary to increase the rotation speed of the compressor or redesign the centrifugal compressor for hydrogen-mixed natural gas in order to achieve a constant flow rate at the outlet of the compressor.

Experimental and numerical study on cavitation pulsating pressure of water-jet propulsion axial-flow pump.

Cavitation may occur in water-jet pump during operation of water-jet propulsion vessel, and once cavitation occurs, the tip clearance pulsating pressure of the impeller may be intensified, resulting in increased vibration of the water-jet propulsion unit. In this paper, the cavitation pulsating pressure characteristics at different positions in the pump are studied by experiment and numerical simulation, and the pulsating pressure characteristics in tip clearance are mainly researched. Based on Star-CCM+ commercial software, unsteady Reynolds-averaged Navier-Stokes equations(RANS) numerical simulation is carried out, and the feasibility of the numerical simulation method is verified by uncertainty analysis. The results show that the cavitation pulsating pressure near the leading edge of the impeller in the tip clearance is the largest. The variation of the tip clearance pulsating pressure with the intensification of cavitation is studied by numerical simulation, and its mechanism is revealed. A dimensionless coefficient of net positive suction head (CNPSH) is proposed, and the study shows that the cavitation pulsation pressure coefficients of pumps of different scales are equal when the working conditions are similar and the CNPSH are equal, which indicates that the cavitation pulsating pressure performance of full scale pump can be predicted by model scale. It is of great significance to evaluate the vibration performance of the full scale water-jet propulsion.

A novel model for low-permeability reservoir numerical simulation considering dynamic capillary force

ABSTRACT Capillary force is not only affected by water saturation but also affected by fluid flow velocity, which is necessary to modify the prior numerical simulation methods. First, a dynamic capillary force testing device was developed, then, an experimental procedure was adopted to measure the dynamic capillary force curve. Three dynamic capillary force curves were obtained under different displacement velocity conditions. Finally, the dynamic capillary force was introduced into the flow equation of low-permeability reservoir numerical simulation. The oil-phase equation and water-phase equation were discretized in three-dimensional space, and the IMPES method was used to solve the water saturation and pressure difference. The calculation results between the traditional numerical simulation method and the novel method under heterogeneous and inhomogeneous conditions were analyzed. The results reveal that the difference in water cut is mainly reflected in the rapid rise of water cut. In the early stage of rapid increase, the water cut calculated by the new method is slightly higher than traditional method. When the water cut is more than 65%, the water cut calculated by the new method is slightly lower than the traditional method.

A Systematic Review of the Current State of Numerical Groundwater Modeling in American Countries: Challenges and Future Research

In arid and semi-arid regions, groundwater is often the only available water source. However, overexploitation and pollution have led to a decrease in groundwater quantity and quality. Therefore, the proper management of groundwater resources is essential to promote sustainable development. Numerical simulation models (NSMs) have emerged as a valuable tool to address these challenges due to their ability to accurately and efficiently model groundwater systems. This study provides a comprehensive systematic review to evaluate the current knowledge on using numerical groundwater flow models for planning and water resource management in countries in the American region. A total of 166 research articles were published between the years 2000 and 2024. We analyzed and summarized details such as the study regions, numerical simulation methods and applied software, performance metrics, modeling units, modeling limitations, and prediction scenarios. In addition, we discuss alternatives to address the constraints and difficulties and suggest recommendations for future research. The continued research, improvement, and development of numerical groundwater models are essential to ensure the sustainability of groundwater resources.

Optimizing Wide Band Gap Cu(In,Ga)Se2 Solar Cell Performance: Investigating the Impact of "Cliff" and "Spike" Heterostructures.

In recent years, the efficiency of high-efficiency Cu(In,Ga)Se2 (CIGS) solar cells has been significantly improved, particularly for narrow-gap types. One of the key reasons for the enhancement of narrow-gap device performance is the formation of the "Spike" structure at the CdS/CIGS heterojunction interface. Wide-gap CIGS solar cells excel in modular production but lag behind in efficiency compared to narrow-gap cells. Some studies suggest that the "Cliff" structure at the heterojunction of wide-gap CIGS solar cells may be one of the factors contributing to this decreased efficiency. This paper utilizes the SCAPS software, grounded in the theories of semiconductor physics and photovoltaic effects, to conduct an in-depth analysis of the impact of "Cliff" and "Spike" heterojunction structures on the performance of wide band gap CIGS solar cells through numerical simulation methods. The aim is to verify whether the "Spike" structure is also advantageous for enhancing wide-gap CIGS device performance. The simulation results show that the "Spike" structure is beneficial for reducing interfacial recombination, thereby enhancing the VOC of wide-gap cells. However, an electronic transport barrier may form at the heterojunction interface, resulting in a decrease in JSC and FF, which subsequently reduces device efficiency. The optimal heterojunction structure should exhibit a reduced "Cliff" degree, which can facilitate the reduction of interfacial recombination while simultaneously preventing the formation of an electronic barrier, ultimately enhancing both VOC and device performance.

The analysis of Gurney flap applications for autonomous underwater gliders

ABSTRACT In this paper, the Gurney flap structure (GF) from aeronautical research was applied to underwater gliders for improving the lift-to-drag ratio. Firstly, the numerical simulation method was used to verify the lift enhancement principle of the GF structure and the feasibility of its application to autonomous underwater gliders (AUG). Secondly, the performance of different GF was simulated and analysed, and a variable GF structure was proposed with the requirements of the operating modes. The GF parameters were optimised by a genetic algorithm. The results show that the performance of GF in water is similar to that in air, so it is feasible to apply GF to AUG. When the height of GF is close to the turbulent boundary layer of the trailing edge of the wing, the lift-to-drag ratio is significantly improved, especially in the range of small angle of attack and the maximum increase is up to 30%.

Study on the Dynamic Crushing Behaviors of Hourglass Honeycomb Sandwich Panels

In response to the problem of enclosed internal spaces in existing honeycomb sandwich panels, the concept of an hourglass honeycomb sandwich panel model is proposed for the first time, which provides a breakthrough approach for achieving the multifunctional integration of honeycomb sandwich panels. Numerical simulation methods are employed to investigate the dynamic performance of the hourglass honeycomb sandwich panels. The focus is on discussing the influences of the geometric parameters on the deformation mode, dynamic response, load uniformity, and energy absorption capacity of the hourglass honeycomb sandwich panel under different impact velocity conditions. The research results indicate that under low-velocity-impact conditions, the influence of the geometric parameters is predominant. In contrast, under high-velocity-impact conditions, the influence of the impact velocity conditions is predominant. Hourglass honeycomb sandwich panels with low density, a large inclination angle of the honeycomb wall, and small contact distances between the hourglass honeycomb cell and the panel have excellent load uniformity, and the distances between the contact points of the hourglass honeycomb cell and the panel have a great influence on the energy absorption capacity of the sandwich panels. This study provides a theoretical basis for the application of honeycombs in aerospace and other engineering areas.

Numerical Simulation and Experimental Verification of Rotor Airflow Field Based on Finite Volume Method and Lattice Boltzmann Method

The primary focus of research in agricultural unmanned aerial vehicle (UAV) pesticide application technology is the investigation of droplet drift and deposition. The influence of the rotor airflow on droplets is particularly significant, making numerical simulations a crucial tool for airflow field analysis. Among existing numerical simulation methods, the Finite Volume Method (FVM) and the Lattice Boltzmann Method (LBM) are commonly used, but there is limited research that compares the two approaches. Therefore, this paper conducts numerical simulations of the rotor airflow of an agricultural UAV using Fluent, representing the FVM, and XFlow, representing the LBM. This research aims to reveal the distribution patterns of airflow field numerical simulations under different theoretical methods, validate them through practical experiments, and select the optimal method for simulating rotor airflow. The ultimate goal is to establish an effective airflow field model to enhance the precision of pesticide application by an agricultural UAV. The results indicate that the lift error calculated by XFlow in this paper is 2.57% smaller than that by Fluent. The wind field of Fluent entered the “stable state” earlier than XFlow, and the speed value of Fluent was smaller than that of XFlow. The difference between the two speed values became larger and larger as the distance from the rotor was longer. Compared with XFlow, Fluent changes more obviously in the core region, and the center region gradually disappears with the distance from the rotor. However, in the velocity field calculated by XFlow, there are still more turbulent flows outside the core region, indicating that the transient calculation method based on the LBM can better show the details of fluid flow than the steady-state calculation method based on the FVM. Through comparison with the actual test data, it is found that the relative error of the velocity value of XFlow at 0.2 m and 0.4 m is small, while that of Fluent at 0 and 0.2 m is small. It shows that XFlow simulation has higher accuracy for external turbulent flow, while Fluent simulation has higher accuracy for steady laminar flow. The research results provide data comparison and a basis for further analysis of the wind field model of the rotor wing of the plant protection UAV, and they lay a foundation for further research on precision application technology.

Structural response and damage assessment method for subsea pipe-in-pipe subjected to anchor impact

The pipe-in-pipe system is widely used in the development and transportation of offshore oil and gas resources. The damage caused by anchor impacts on the pipe-in-pipe system has long been a focal point and challenge in both engineering and academic fields. Numerical simulation method is employed to investigate the impact process of anchors on pipe-in-pipe systems, revealing the variation patterns of pipeline deformation response over time. The structural response of the pipe-in-pipe system to dropped anchor impacts is studied from multiple perspectives, including deformation, energy, and load transfer. The influence of pipe-in-pipe characteristic parameters, anchor characteristic parameters, centralizers, and bulkheads on the extent of impact damage is analyzed, establishing a relationship between the absorbed energy and the damage to the pipe-in-pipe. The results indicate that the impact energy primarily depends on the mass and velocity of the anchor. Energy is conserved during the impact process. It is found that even with different anchor masses and velocities, the extent of damage to the pipe-in-pipe system remains the same if the impact energy is identical. High-strength steel, large diameter, and thick-walled pipelines exhibit better resistance to impact loads. Additionally, centralizers and bulkheads mitigate the impact damage to the pipe-in-pipe system. Finally, an evaluation method for pipe-in-pipe impact damage due to dropped anchors is proposed. It is discovered that the assessment method for pipelines impacted by dropped anchors in the DNV-RP-F107 standard tends to be conservative when the impact energy is high. The findings of this study provide technical support for the integrity management of pipe-in-pipe systems.

Research on the Arrangement Scheme of Spirally Twisted Tape Inserts in a Grate Furnace

To eliminate the flow dead zone and homogenize the asymmetric flow field of a grate furnace, spirally twisted tape inserts (STTIs) with a pitch ratio of 1.5 were installed in the vertical flues of an SCL1000-13.5/450 grate boiler. The arrangement schemes found to be present inside the chosen 1000 t/d grate furnace, determined using the orthogonal experimental method, included separate installation in chamber II, separate placement in chamber III, and simultaneous arrangement in both chamber II and chamber III. The effects of row spacing H, column spacing W, and mounting angle φ were investigated by means of the practicable and feasible numerical simulation method. With a focus on the uniformity degree of the flue gas, the results showed that temperature distribution is directly correlated with the velocity field. When it comes to the uniformity of the flow field, the exclusive use of STTIs in chamber II was better than that in chamber III. Under the optimal combination scheme of STTIs in both chamber II and chamber III (scheme N323), the exhaust gas temperature reached the minimum value and the uniformity index of temperature increased to the range of 0.994~0.997. The findings in this work could provide a reference for the optimization of the flow field in a grate furnace.

Study on Surrounding Rock Control of Withdrawal Space in Fully Mechanized Caving Mining of a 19 m Extra-Thick Coal Seam

The section span of the withdrawal space of fully mechanized top coal caving in an extra-thick coal seam is large, and with the gradual withdrawal of the hydraulic support, a series of strong dynamic pressure disasters occur in the withdrawal space, and the difficulty of surrounding rock support control increases sharply. In order to study the control mechanism of surrounding rock in the final mining withdrawal space in detail and put forward a reasonable support technology scheme, taking the large-section withdrawal space of an 8309 fully mechanized caving face in an extra-thick coal seam of a mine as the research object—through the theoretical investigation of whether the key blocks of the main roof are stably hinged under varied stopping coal caving distances and fracture locations of the main roof—the reasonable and optimal stopping coal caving distances and roadway formation time are determined. Using numerical simulation and similar simulation methods, the vertical stress and the maximum shear stress research indicators were introduced to verify the accuracy of the theoretical analysis results. The results show the following: (1) The reasonable stopping coal caving span is 1~2 times the cycle weighting interval, the best stopping coal caving distance in this geological condition is 30 m, and the best fracture position of the main roof is located above the goaf. (2) The migration of overlying strata in the withdrawal space has obvious zoning characteristics, and the zoning is as follows: a stopping coal caving area, support area of the hydraulic support, withdrawal channel area, and stopping coal pillar area. (3) According to the zoning characteristics of overlying strata movement, the asymmetric zoning support control scheme of the withdrawal space is proposed. The field monitoring results show that the maximum roof subsidence in the withdrawal space is 151 mm, the maximum internal squeezing amount of the stopping coal pillar is 82 mm, and the supporting and anchoring effect of each partition in the withdrawal space is good. The set of partition asymmetric support control schemes has been successfully applied to field practice.

The electrochemical characteristics of derivatives of anthracen-9-ylmethylene at a glassy carbon electrode in methylene chloride

This study reports the electrochemical investigation of the Anthracen-9-ylmethylene derivatives (anthracen-9-ylmethylene-(3,4-dime-thyl-isoxazol-5-yl)-amine) compound on a glassy carbon electrode (GCE), employing chronoamperometry, chronocoulometry, convolutive voltammetry transforms, and numerical simulation techniques. The experiments were performed in a solution containing 0.1 mol L−1 tetraethylammonium perchlorate (TEAP) in methylene chloride solvent (CH2Cl2). The oxidation process of the compound involved two sequential electron transfers. The first redox process involves the loss of a single electron transfer, leading to the formation of a radical cation followed by an irreversible chemical step, and the subsequent transformation of electron transfer produces a di-cation through an additional irreversible chemical process. The experimental kinetic parameters (α, ks, E°, D, and kc) were determined via the electrochemical techniques used. The numerical simulation method was used for confirmation and verification of the determined chemical and electrochemical parameters. The electrode pathway was found to follow an ECirrECirr mechanism, providing valuable insights into the examination electrochemical performance of Anthracen-9-ylmethylene derivatives' on a GC electrode. Examination of the anthracen-9-ylmethylene derivatives via chronocoulometry experiments confirmed that the mass transport of the investigated species is controlled by the diffusion process.

Cooling characteristics of nanofluids in a snowflake and squid fin composite bionic microchannel heat sink

To improve the heat dissipation performance of electronic components, this study proposed a novel biomimetic microchannel design inspired by snowflakes and squid fins. A numerical simulation method was employed to investigate the cooling characteristics of nanofluids at different concentrations within this biomimetic microchannel heat sink. The microchannel is shaped like a snowflake, while the base plate features a biomimetic squid fin structure. The study focused on analyzing the effects of waveform height (amplitude A = 0.1–0.4 mm) and waveform frequency (B = 0.5–2π) on heat transfer performance. The results indicated that when nanofluids flow over a base plate with a wavelength frequency of B = 0.5π, there is a significant enhancement in heat exchange while simultaneously reducing flow resistance. In the microchannel heat sink (A = 0.3 mm, B = π), the thermal conductivity of the nanofluids can be increased by up to 32.56 %, with a maximum comprehensive evaluation index of 1.68. This research provides important guidance for the application of nanofluids and biomimetic microchannels in the thermal control of electronic devices.