Numerical Simulation Method Research Articles

A Numerical Simulation Method for Investigating the Fluid–Structure–Ice Coupling Mechanism of a Wedge Breaking through Ice into Water

A Numerical Simulation Method for Investigating the Fluid–Structure–Ice Coupling Mechanism of a Wedge Breaking through Ice into Water

Shock buffet onset prediction with flow feature-informed neural network

Transonic shock buffet is a significant self-excited shock oscillations and aerodynamic instability phenomenon induced by shock-boundary layer interaction, which limits the flight envelope and even causes flight accidents. The aviation industry has a significant interest in accurately predicting the shock buffet onset boundary, defined by a specific combination of Mach number and angle of attack. While the current methods of steady and unsteady numerical simulation suffer from a contradiction of efficiency and accuracy. In the current paper, a flow feature-informed neural network (FINN) model is constructed to predict the buffet onset boundary over airfoils. Typical features associated with buffet onset are extracted from the steady base flow and subsequently integrated into the hidden layer of the neural network to impose physical constraints. With the test cases of the NACA0012 airfoil at various Mach numbers, the FINN model can accurately predict the damping representing the unsteady instability margin. Compared to the direct mapping input-output neural network (NN) model, the proposed method with shock wave feature-informed has enhanced accuracy, with an average relative error decreased by 70% at extrapolated Mach numbers. This research demonstrates the effectiveness of the FINN model in predicting the buffet onset, which leverages physics features derived from the more economical steady solution far from the onset boundary at a given predicted Mach number.

A method for numerical simulation of shock waves in rarefied gas mixtures based on direct solution of the Boltzmann kinetic equation

The paper proposes an approach to numerical simulation of shock waves in rarefied gas mixtures on the basis of direct solution of the Boltzmann kinetic equation. Software for simulating the gas flows was developed. The structure of a shock wave in a binary gas mixture was computed with an accuracy controlled by the computational parameters. The computations were performed for various molecular masses ratios and Mach numbers. The total accuracy of at least 1.4% for the local values of the molecular densities and temperatures of the mixture components was achieved. Numerical simulation of a shock wave propagation through a periodically perforated surface was performed. The distributions of the macroscopic characteristics of the mixture components at various points in time were obtained. Unsteady areas of strong separation of the gas mixture components were discovered.

Unsteady numerical simulations of aircraft icing based on the dynamic grid method

In this study, an unsteady numerical simulation method based on dynamic grids was established to solve unsteady aircraft icing processes. The dual-time step method was used to achieve the time advancement of the transient airflow and droplets fields. The droplets trajectory was calculated using the Euler method, which incorporates a supercooled large droplets model. The transient evolution of the water film and ice layer on the substrate surface was considered in the icing thermodynamic model. The accuracy of the unsteady method for ice-shape prediction was verified via comparison and analysis of numerical examples. Additionally, the rime and glaze ice formation processes were explored, and the coupling effects of ice horn growth and the air convective heat transfer coefficient, water collection coefficient, and water film flow on the substrate surface were studied. The ice shapes calculated by unsteady and quasi-steady methods were compared. The results indicated that the latter predicts a larger icing range than the former. Moreover, there may be more obvious distinctions between these two methods in predicting glaze ice.

Research on Chip Removal Performance of Rotary Percussive Ultrasonic Drill for Planetary Rock Sampling

The objective of the kinematic model of rock chips is to examine the impact of the driving parameters of drilling tools on the effective chip removal speed, with a view to analyze the chip removal performance of rotary percussion ultrasonic drilling. The discrete element numerical simulation method, EDEM-2022 simulation software, was employed to establish a simulation model, match the simulation parameters, and conduct a simulation study of chip removal under the action of rotary percussion. This was performed to analyze the influence of the driving parameters on the chip conveying process. In order to ascertain the order of parameters affecting the degree of conveyance rate as well as the optimal driving parameters, orthogonal experiments were designed based on Taguchi’s experimental analysis method.

Optimization Study on Nozzle Selection Based on the Influence of Nozzle Parameters on Jet Flow Field Structure

Currently, the primary method for controlling red tides in the ocean involves spraying water solutions with special chemicals as solutes. High-pressure spraying results in the formation of typical jet structures. In this study, numerical simulation methods are employed to investigate the velocity variations, turbulent characteristics, and gas content distribution of jet flow fields under different initial jet flow pressures, cone angles, and nozzle diameters. Based on practical application scenarios, cluster analysis is used to explore the similarities and differences in jet equivalent diameters under different parameter conditions. The research findings indicate the following. (1) The difference of jet velocity distribution at the far field exit will be enlarged with the increase in the nozzle cone angle. When the nozzle cone angle is 4 mm, the velocity uniformity at the outlet is the best. (2) The TKE of the flow field has no consistent change law along the central axis. At the jet exit, the TKE shows an obvious multi-peak structure. (3) The gas content demonstrates a typical “double-valley” feature at the jet outlet cross-section. Increasing the initial pressure leads to a decrease in the gas content within the jet due to reduced entrainment, while the entrainment range remains largely constant. (4) Cluster analysis reveals that the similarity of jet flow width when it reaches the water surface is minimal compared to other operating conditions when the initial pressure is 0.36 MPa, the cone angle is 115°, and the nozzle diameter is 2 mm. All conditions can be categorized into two or three groups to ensure jet effectiveness. The study results provide scientific guidance for selecting spray devices for controlling red tides in the ocean.

Study on the influence of slope morphology on seismic dynamic response based on numerical simulation

Geological disasters frequently occur in our country. Multiple earthquakes have caused landslides on slopes, resulting in serious casualties and property losses. Therefore, it is particularly important to study the dynamic response and instability mechanism of slopes under earthquakes of different frequencies. This paper uses numerical simulation methods to establish a slope model, loads sine waves of different frequencies and amplitudes, collects acceleration and displacement dynamic responses of different profiles, and studies the distribution law of slope acceleration and displacement under dynamic action. Finally, the phenomenon of acceleration amplification effect is discovered. The research results of this paper can provide a basis for slope disaster prevention and mitigation and seismic reinforcement measures.

A Review of the Degradation Research on the Single-Lap Bolted Joint

Bolted joints, with their advantages of simple structure and convenient disassembly and assembly, are widely used in complex equipment fields such as aerospace systems and weaponry. Subject to complex mechanical loads, the contact surfaces may undergo nonlinear behaviors such as contact–separation and viscous slip, leading to the nonlinear degradation of the connection stiffness, which severely threatens the safety and reliability. These have driven research on bolted joints to span multiple disciplines, from interfacial micro-friction to macro-structural dynamics. Therefore, from the field of micro-friction to macro-dynamics, this review summarizes and analyzes three major degradation models, outlines the experimental development in connected structures, and provides an overview of the numerical analysis methods for degradation simulation. This paper also looks forward to the development directions for future research on the degradation of connected structures.

Overview of Theory, Simulation, and Experiment of the Water Exit Problem

The water exit problem, which is ubiquitous in ocean engineering, is a significant research topics in the interaction between navigators and water. The study of the water exit problem can help to improve the structural design of marine ships and underwater weapons, allowing for better strength and movement status. However, the water exit problem involves complex processes such as three-phase gas–liquid–solid coupling, cavitation, water separation, liquid surface deformation, and fragmentation, making it challenging to study. Following work carried out by many researchers on this issue, we summarize recent developments from three aspects: theoretical research, numerical simulation, and experimental results. In theoretical research, the improved von Karman model and linearized water exit model are introduced. Several classical experimental devices, data acquisition means, and cavitation approaches are introduced in the context of experimental development. Three numerical simulation methods, namely, the BEM (Boundary Element Method), VOF (Volume of Fluid), and FVM (Finite Volume Method) with LES (Large Eddy Simulation) are presented, and the respective limitations and shortcomings of these three aspects are analyzed. Finally, an outlook on future research improvements and developments of the water exit problem is provided.

Research on layout optimization for propellers of a submarine operation vehicle

The submarine operation vehicle was an important equipment for exploring the marine environment and performing underwater tasks. In this paper, the calculation and analysis of hydrodynamic performance and the layout optimization of the propellers under different layouts were completed. Based on the left-right symmetrical layout of the four ducted propellers, the two ducted propellers were taken as the research object to optimize its vector layout, and an approximate model was constructed according to the numerical calculation results under different combinations of propeller spacing and deflection angle. The simulation calculation and verification are carried out, and the optimal layout of the four ducted propeller were determined considering the interference between the ducted propellers on both sides. Finally, an experimental platform for the thrust test of the two ducted propellers were built. It was found that the axial thrust ratio of the rear propeller and the front propeller in experimental results were in good agreement with the simulation results, which verified the feasibility of using the numerical simulation method to study the mutual interference between the ducted propellers.

Verification of frequency-dependent impedance boundaries for the Finite-Volume Time-Domain scheme

Numerical simulation methods are widely used to approximate solutions to the wave-equation, taking into account the initial condition ('initial sources') and boundary conditions ('wall materials'). Trustworthy simulations are essential; a primary step is to ensure the correctness of the implementation through a process known as code verification. In this work, we do code verification by studying the convergence of the known finite-volume time-domain (FVTD) method with a frequency-dependent boundary condition modeled as an inductance/resistance/capacitance (LRC) termination. Due to the difficulties of obtaining an analytical solution for frequency-dependent impedance boundaries, we thoroughly derive and implement The Method of Manufactured Solutions tailored to the wave-equation inside a rectangular room for uniform frequency- dependent wall impedance. Asymptotic estimates are found empirically based on a convergence study compared against the manufactured ground truth solution. We investigate the impact of the frequency-dependent boundaries on the convergence affected by the total sum of numerical errors and validate the theoretical properties of the FVTD scheme for various LRC model configurations. The method should also serve as a useful resource for others to verify the correctness of their code implementations.

Research on concrete-filled square CFRP steel tube under compressive-torsional hysteresis loading

Concrete-filled CFRP steel tube structures are presently employed on extensively in various engineering practices. This study examines the mechanical performance of concrete-filled square CFRP steel tubes under compressive-torsional hysteresis loads. Through the analysis of nine specimens, the research focuses on failure modes, torque–angle (T–θ) behavior, triaxial strain, and the synergistic effects of the steel tube and CFRP during cyclic loading. A numerical simulation method is initially introduced based on the test data to predict the actions of concrete-filled CFRP steel tube compressive-torsional specimens under hysteresis loading. This method is subsequently verified by contrasting it with the outcomes from representative experiments. In addition, the stress distribution of each material during the entire loading process of the sample is examined. Finally, using data from experiments and finite elements, combined with the trilinear model, a hysteretic model under compressive-torsional hysteresis stresses has been established for square CFRP steel tube filled with concrete. The calculated outcomes of the developed finite element model closely align with those of the hysteretic model.

Adaptive-coefficient finite difference frequency domain method for time fractional diffusive-viscous wave equation arising in geophysics

The diffusive-viscous wave (DVW) equation is a widely used model to describe frequency dependent attenuation of seismic wave in fluid-saturated porous medium. In this paper taking power law frequency dependent attenuation into account, we first introduce a modified DVW equation (time fractional DVW equation) in which the first order temporal derivative of viscous term is replaced with a fractional order temporal derivative. In consideration of that most of the existing numerical simulations for seismic wave equations are based on time domain methods and truncation with some specific boundary conditions, we incorporate the absorbing boundary condition as complex-frequency-shifted (CFS) perfectly matched layer (PML) into the time fractional DVW equation, and then develop an adaptive-coefficient (AC) finite difference frequency domain (FDFD) method for numerical simulation. The corresponding analytical solution for homogeneous time fractional DVW equation is provided for model validation, and the effectiveness of the developed AC FDFD method is verified by some numerical examples including the homogeneous model and the layered model. Numerical results show that AC FDFD method is more accurate than the traditional 2nd-order FDFD method for numerical modeling of time fractional DVW equation with CFS PML absorbing boundary condition, while requiring similar computational costs.

Exploration and Practice of Synergistic and Efficient Development Technology for Gas-Capped Condensate Gas Reservoirs and Oil Rings

Abstract The oil-ringed condensate gas reservoir is one of the more unique and complex types of oil and gas reservoirs. The stable oil and gas interface is crucial for the efficient development of such reservoirs, and the fundamental means of addressing the issue lies in maintaining the dynamic balance between the gas cap and the oil ring. In this paper, we address issues such as the rapid decline in oil ring production, gas cap breakthrough, and increased condensate oil loss exposed during the development of the R carbonate reservoir in the Caspian Sea Basin of Kazakhstan. Using reservoir engineering and numerical simulation methods, we optimized synergistic development ideas and the boundaries of exploitation policies for gas-capped condensate gas reservoirs and oil rings. We formulated reasonable development parameters such as extraction rate, injection-production ratio, barrier water injection ratio, and proposed dynamic control standards for zoning and classification. After implementing the plan, the gas-oil ratio of the oil ring significantly decreased, with a reduction rate of 12 percentage points. The gas production per unit pressure drop increased by 40%. The annual decline rate of condensate oil decreased from 23.6% to 13.8%. All development indicators exhibited characteristics of benign development, providing a theoretical foundation and technical methods for guiding the efficient development of similar oil and gas reservoirs.

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.1molL-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, Eo, 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.

Fracture performance of concrete-epoxy mortar interface subjected to bending and shear: Experimental and numerical study

To investigate the fracture performance of concrete-epoxy mortar interfaces, experimental and numerical simulation methods were employed to obtain the failure patterns and applied load P versus opening displacement curves of specimens. A special emphasis is placed on the discussions of crack initiation load (Pini), maximum load (Pmax), corresponding deformation (CTODini and CTODmax, CTSDini and CTSDmax), KR and GR curves of I-Type and stress intensity factors of I-II-Type cracks (KIini and KIIini) and facture energy (GF). Test results show that increasing both concrete compressive strength and interface roughness enhances the load carrying capacity and crack initiation toughness, as well as crack extension resistance. Higher concrete compressive strength improves the moment capacity of the interface and increases the proportion of KIini. An increase in the interface roughness enhances the mechanical interlocking interaction at the concrete-epoxy mortar interface, thereby decreasing the CTSD/CTOD ratio and increasing the proportion of KIIini for the I-II Type crack. As the mode ratio (KIIini/KIini) increases, the proportion of shear and bending moment acting on the interface changes continuously, with the shear effect becoming more dominant. Basically, the increase in both the crack-height ratio (a0/D) and initiated crack offset ratio (C1/C0) decreases the fracture performances of specimen. The observations provide a deeper understanding to the design and practice of concrete-epoxy mortar interface.

Study on the transient cavitation characteristics in an axial flow tandem pump during the rapid starting period

Abstract In axial flow pump systems constrained by limited internal space, the advantages of axial flow tandem pumps are high space utilization, strong power capacity, and great thrust. The onset of cavitation during the rapid starting period poses a classic issue for tandem axial flow pumps, resulting in issues like noise, vibration, cavitation breakdown, and material damage. This article aims to investigate the transient cavitation characteristics of each flow passage component of an axial flow tandem pump during the rapid starting period by numerical simulation method. It is compared that the evolution laws of the cavity of each flow passage component. The evolution of the cavity in all flow passage components is strongly correlated with the flow rate change rate in each stage of the rapid starting period. Among all flow passage components, the first impeller firstly experiences cavitation, and its cavitation severity is also the most severe. The Effective Net Positive Suction Head (NPSHa ) can reduce the influence of transient effects on cavity prediction and better characterize the details of the cavity evolution in different impellers than the cavitation number σ during the rapid starting period. Meanwhile, The NPSHa of the second impeller inlet is always higher than that of the first impeller inlet, indicating superior cavitation performance for the second impeller.

Dynamic mixing characteristics of a soft elastic tubular reactor (SETR): Experimental investigation and numerical simulation

In an innovative effort to mitigate material damage and enhance overall mixing performance, the use of a soft elastic tubular reactor (SETR) was explored. In this study, the mixing of tracer in SETR was observed with the help of non-invasive method. The effect of different parameters such as pusher squeezing frequency (F), distance between pushers (D) and pusher squeezing mode (M) on the flow field inside the SETR was also analyzed with the help of numerical simulations. The study also further optimized the operating parameters of the reactor. The experimental results proved the accuracy of the subsequent simulation model. Through the combination of response surface design and numerical simulation methods, it was found that all three factors have a significant effect on the mixing time. The results of this work add to our understanding of the mixing mechanism within SETR and present new opportunities for designing highly efficient industrial reactors.

Fault zone architecture and hydraulic properties characterization using earth tide analysis

Fault zone architecture and hydraulic properties, which determine subsurface fluid flow and storage, have received long-term research attention from hydrogeologists. Current field experimental techniques and numerical simulation methods are limited by the scale and complexity of the investigations and are inadequate for deep subsurface fault zones. Estimating the internal rock properties and structures of fracture zones requires a combination of methods from many different fields. In recent years, groundwater tidal response methods have been widely used to characterize subsurface hydrogeomechanical properties because of their advantages of being in situ and low cost. In this study, tidal analysis was employed for a well within the Blue Ridge Fault Zone, and we integrated the fracture zone properties revealed by a single well. Well W-03 revealed an aquifer with an overall permeability of ∼ 10–14 m2 and hydraulic diffusivity of 0.24 m2⋅s−1. The fault damage zone had a permeability of 10–14 ∼ 10–13 m2 while the host rock had a permeability of 10–16 ∼ 10–14 m2. In addition, we found that the heterogeneity of specific storage affects the assessment of poroelasticity parameters, we offer a method to estimate these parameters without relying on specific storage, yielding porosity estimates ranging from 0.05 to 0.06, aligning with field test results. Furthermore, the estimated dip and strike of the fault zone were also consistent with the outcrop and borehole observations.

Jet centrifugal pump internal flow numerical simulation and structural optimization

Abstract Jet centrifugal pump is essential equipment in multiple systems such as agricultural irrigation and energy production. This paper investigates the internal flow characteristics of jet centrifugal pumps at different flow rates and the cavitation behavior under high flow conditions by numerical simulation methods. Additionally, two performance optimization strategies are proposed based on simulation results. The findings indicate that high-velocity regions and low-pressure zones exist within the nozzle. As flow rate increases, the outlet velocity of the nozzle gradually decreases, and large areas of low pressure are observed near the inlet of the impeller. Cavitation occurs to varying degrees at flow rates of 4.056 m3/h and 3.840 m3/h. The turbulent kinetic energy contours for these conditions reveal uneven distribution within the jet nozzle, with higher values near the jet mixing layer, inner wall of the jet nozzle and impeller inlet. Structural optimization using hemispherical bolts significantly reduces the turbulent kinetic energy around the bolts, resulting in an increase in head ranging from 1.68% to 6.70% across various flow rates. Shortening the clearance between the impeller and the rear cover reduces energy losses, leading to an increase in head ranging from 2.72% to 5.86%.