Results Of Numerical Modeling Research Articles

Linear transverse flux generator for wave energy conversion: design optimization and analysis

Abstract The electric power industry impacts each state’s economy significantly, driven by increasing electricity consumption that necessitates expanding power plants and finding alternative energy sources. Among alternative energy sources, ocean and sea wave energy converters can be distinguished as a separate class. Wave energy converters transform wave energy into mechanical and then electrical energy. The purpose of the study is to analyze and optimize the magnetic system of a transverse flux machine (TFM) linear generator and to determine the influence of the distance between the stator cores on the efficiency of the generator. This research included conducting 3D modeling and analysis to identify this rational distance. The methods for investigating the magnetic system and calculating the magnetic field pattern are divided into analytical and numerical. Thanks to advanced software for solving such tasks, numerical calculation methods based on the finite element method play a decisive role. Meanwhile, analytical calculations of the magnetic circuit are performed using Kirchhoff’s second law for preliminary analysis. The article discusses a two-phase linear TFM generator with a U-shaped core and permanent magnets. The results of numerical modeling show that the distance between the stator cores should have a specific size and requires detailed selection when designing the magnetic system in each particular case. In the design studied, it was calculated that 6 mm between the stator cores increases the machine’s performance. 3D modeling is necessary for accurate analysis, considering the axial magnetic flux to minimize stray fields and their mutual demagnetization. Future research will explore an E-shaped core TFM design.

Peculiarities of flexural wave propagation in a notched bar

The results of numerical modeling and experimental studies of the propagation of flexural elastic waves in a metal notched bar approximates the effect of an acoustic black hole are presented. The sample is a bar with notches, the depth of which increases according to the power law with an exponent equal to (4/3). It has been confirmed experimentally and with the simulation results, that such bars slow down the propagation of an elastic wave to the end of the bar. It is shown in such structures flexural waves have dispersion and their localization at the end of the bar is higher for some natural frequencies than that of a solid rod. The natural oscillations of the whole and notched bars are compared, i.e. the shape of the amplitude of the flexural wave along the rods. The dependence of the flexural wave length in a notched bar on the frequency is investigated as a wave propagates to the end of the bar.

Phase synthesis method in multi-aperture optical systems based on iterative image processing algorithms

Background. Modern technologies in the field of optics and photonics place high demands on the quality and output power of the radiation source. The fulfillment of these requirements can currently be achieved by creating multi-aperture laser systems with coherent beam addition. The main problem of creating such systems is the development of phase synthesis methods in a multichannel optical system. Aim. In this paper, we consider a phase synchronization method without a reference beam for a coherent addition system with active feedback, which is based on the Gerchberg–Saxton algorithm. This algorithm makes it possible to reconstruct complex field amplitudes in the aperture and focal planes from intensity distributions in these fields. The application of this algorithm for multichannel systems is analyzed and its features, such as the occurrence of stagnation and ambiguity of the solution, are revealed. Methods. Solutions are proposed to eliminate the problem of convergence of the algorithm to side solutions and getting the iterative procedure into the local extremum using global optimization algorithms, methods of reducing the dimension of the problem and the introduction of antisymmetric amplitude modulation. Results. The paper demonstrates the results of phase field reconstruction for a large number of optical sources. For a seven-aperture system, a physical simulation was performed to restore phase information, confirming the results of numerical modeling. Conclusion. The proposed approach is a reasonable alternative to currently existing methods and can be used in the problem of coherent addition in multichannel optical systems.

Effect of the Effective Metal Surface Area of Two Different Flow Diverter Stents on the Stagnation Region Formation Inside the Aneurysm Sac

Abstract Objective Flow diversion (FD) is a relatively new technique for treating large, wide-necked, or fusiform aneurysms. Although FD is a more preferred option than coiling or clipping techniques in neurosurgery and neuroradiology clinics, the blood flow mechanism inside the aneurysm sac is not fully understood after the treatment. Besides, effective metal surface area (EMSA), a property of an FD related to porosity, shows variation at the patient's aneurysm neck by providing more or less blood flow inside an aneurysm sac than planned, causing nonstagnant or stagnant fluid region formation in the sac, respectively. Thus, the change in FD's EMSA can significantly affect the treatment's effectiveness, making even operation unsuccessful when variation in FD's EMSA at the aneurysm neck is overlooked. Materials and Methods In this study, a large aneurysm of a 52-year-old female patient was numerically investigated by virtually placing two commercially available FDs with different EMSA values one by one into the aneurysm-carrying artery. Results While FD stents at the aneurysm site substantially reduced the blood flow into the aneurysm, an FD with a 15.6% EMSA caused blood to flow in the aneurysm sac to have six times more kinetic energy than that of FD with a 29.5% EMSA. Conclusion Although FD's EMSA value demonstrated nearly up to 20% reduction at the patient's aneurysm neck based on a product catalog value, numerical model results revealed that the stagnated region's formation inside the aneurysm sac could be determined within a 9% difference based on digital subtraction angiography reformat image.

Stationary regimes and parametrization of ekman friction in the Karman model of flow induced by external vortical body force

The detailed study of stationary regimes of Karman axisimmetric flow induced by external vortical body force is done. It is extracted two stationary regimes – with small (Batchelor regime) and with substantial (Stewartson regime) secondary circulation. The diagram of regimes existence is plotted in the space of flow parameters – Rossby and small Ekman numbers. For the flow decaying to the stationary flow in the Batchelor regime a theoretical model is proposed with which it was possible to derive a parametrization of linear friction coefficient, Ekman pumping velocity, stationary pressure from mean flow characterictics (vorticity and divergence). In the Stewartson regime a parameterization of the stationary flow is proposed and also numerically studied a decay rate. It is shown a good agreement between theoretical and numerical model results.

Dynamics of the plasma ejection structure in laboratory modeling of young star jets at plasma focus facilities

The use of “plasma focus” type facilities, such as PF-3 (Kurchatov Institute), makes it possible to carry out well-controlled and diagnosable laboratory experiments to study laboratory jets with scale parameters close to the jets of young stars. In this paper, we present the results of numerical modeling of plasma outburst propagation in PF-3. A self-consistent configuration was chosen as the initial conditions, which correctly takes into account the internal structure of the jet. This allowed us to obtain a detailed structure of the interaction between the magnetized emission and the ambient gas. Due to the scalability of such a structure, one should expect such a structure from the head shock waves of jets of young stars.

Assessment of Tsunami Hazard for the Coast of Oktyabrsky Spit, Western Kamchatka Based on Results of Numerical Modeling

Assessment of Tsunami Hazard for the Coast of Oktyabrsky Spit, Western Kamchatka Based on Results of Numerical Modeling

Shape of fragments cloud behind heterogeneous screen by a space debris particle impact

Technogenic pollution of near-Earth space poses a threat to the functioning of spacecraft. Collisions between space debris (SD) and spacecraft (SC) structures can have catastrophic consequences or cause localized damage, leading to the loss of SC operability or the failure of certain functions. The SC body must effectively protect the internal equipment from various external impacts, be technologically feasible to manufacture, and have as little mass as possible. As a result, the task of designing spacecraft bodies and protective screens for low-orbit SC is particularly relevant due to the large concentration of SD in low Earth orbits.A comparative numerical study was conducted to evaluate the effectiveness of various thin shields in protecting against impacts from space debris particles. The study examined homogeneous shields made of A356 aluminum alloy and 316L stainless steel, as well as volumetrically reinforced composite shields produced using additive manufacturing with steel inclusions, and shields with a gradient distribution of steel throughout the thickness of an aluminum matrix, all with the same areal density. In all the heterogeneous plates considered, the volumetric concentration of steel was 36 %. The study covered an interaction velocity range of 2–9 km/s. Numerical modeling results indicated that the structure of the thin heterogeneous plate does not affect the shape of the debris cloud formed behind the protective shield. The findings of this study can serve as a basis for selecting materials for the development of more effective protection for spacecraft against high-velocity impacts.

A New Methodology for Quantitative Risk Assessment of CO2 Leakage in CCS Projects

Summary Assessing risks during carbon dioxide (CO2) sequestration involves a systematic process to identify them (risk assessment), analyze, prioritize, score them (risk analysis), develop mitigation strategies, and plan for measurement, monitoring, and verification (MMV) activities. In this context, risk assessment and analysis are often conducted involving qualitative approaches and/or subjective evaluations based on experts judgments. In this paper, we propose to use numerical results of dynamic subsurface modeling to assess risks in a more precise, less subjective, and quantitative manner, enabling a comprehensive definition of effective MMV plans. In this context, MMV evolves into modeling, measurements, monitoring, and verification [(M)MMV]. The conceptual model used for this study integrates flow and geomechanical-coupled behavior to assess caprock integrity, faults reactivation, and ultimately the risk of CO2 leakage. Uncertainty analysis and a 4D probabilistic workflow are used to quantify the likelihood and consequences of CO2 leakage and to score risks. The proposed approach offers a general framework to enable subsurface-modeling-based quantitative risk assessment (QRA) for the definition of successful (M)MMV plans. In this study, we have focused on the QRA of leakage risks through the fault and intact caprock. Other applications are possible, such as: (1) leakage risks through the caprock resulting from caprock failure; (2) lateral CO2 migration leading to permit limits violation or inducing interferences with neighboring activities; (3) leakage risks resulting from loss of well integrity; (4) risks of faults reactivation; (5) risks of surface heave; and (6) risks of induced seismicity. All these results can be used to formulate effective (M)MMV strategies for an optimized and cost-effective risk mitigation.

Multiprobe Planar Near-field Range Antenna Measurement System with Improved Performance

This article presents a novel multiprobe planar near-field range (PNFR) measurement system. The said system simplifies the overall mechanical design, making it simpler than the existing scanning probe PNFR measurements, and also significantly reduces testing time. A dielectric-based probe is introduced to reduce the antenna size, thereby improving resolution. The probe under consideration is an antipodal Vivaldi antenna offering broadband support and ensuring wideband characteristics of aerials. Numerical results for representative X-band antenna models, presented in the Matlab environment, demonstrate robust performance of the developed measurement system.

Measurement of pulse current during high voltage tests using Rogowski coil with non-inverting integrator

The work is devoted to the development and modelling of the system of pulse current measurement during high-voltage tests and experiments. The proposed measurement system is based on the application of Rogowski coil with non-inverting integrator. The results of numerical modelling of the developed current measurement system are compared with the results of application of a low-inductance shunt resistor, Rogowski coil with passive RC-integrator and Rogowski coil with inverting active integrator in similar measurement systems. Various standard pulse currents including steep edge pulse (1/20 µs), large amplitude current pulse (4/10 µs), lightning pulse current (8/20 µs), switching pulse current (30/80 µs), and 1-3 ms pulse currents are modelled for the comparative test. From the simulation results, it follows that the developed Rogowski coil with the proposed non-inverting integrator has sufficiently high accuracy in a wide frequency band overlapping the frequency band of standard pulse currents in high-voltage tests.

Three-dimensional numerical simulation and experimental validation of airborne ground-penetrating radar based on CNCS-FDTD method under undulating terrain conditions

Traditional Finite Difference Time Domain (FDTD) approaches face challenges with increased computational demands and errors as terrain complexity and flight altitude rising. This study introduces the Crank-Nicolson-Cycle-Sweep-FDTD (CNCS-FDTD) method, enhancing airborne ground penetrating radar (GPR) simulations over undulating terrains. CNCS-FDTD, as an unconditionally stable implicit algorithm, overcomes these by allowing larger time steps without the constraints of the Courant-Friedrichs-Lewy (CFL) condition. Our research aims to assess how CNCS-FDTD can improve computational efficiency and accuracy in modeling airborne GPR responses across varied terrains. Initial simulations using a three-dimensional graben model indicate that CNCS-FDTD can maintain calculation stability with significantly larger time steps, reducing computational time by 42%. To further verify the reliability of the numerical simulation results, the paper also presents experimental tests of the undulating terrain model. By comparing the numerical and physical simulation results of models under different flight altitudes and terrain conditions, the accuracy of the simulation results is validated.

Dynamic Characteristics of a 1950s Heritage Building: A Comparison of Original Design Methods and Modern Techniques

Research on design rules and methods for architectural heritage is an important aspect of conservation practice. Nevertheless, efforts to recover and divulge design methods for Modern Heritage remain limited. This paper is related to the recent structural assessment of a 15-storey heritage building built in 1950, during which a document describing the original seismic analysis of this structure was identified. The methodology employed is of particular interest, as it involves the application of pioneer concepts of dynamic analysis in the design of the first tall buildings in Mexico. The primary aim of this paper is to review the seismic design criteria for the case under study in order to contribute to the state of the art in Modern Heritage. The review includes a comparison between the dynamic characteristics estimated during the design and the results of recent ambient vibration tests and numerical modeling. Several sources of error among the design criteria were identified. Notably, the fundamental period estimated during the design was 38% larger than the experimental value due to an underestimation in stiffness, which introduces significant uncertainty into the design. Overall, the review shows the evolution of seismic analysis over time and provide valuable insights for the study of similar buildings.

The impact of the destruction of the Kakhovka reservoir dam on the oceanographic conditions in the north-western part of the Black Sea according to the results of modeling

Based on the results of numerical modelling, characteristic features of distribution of transformed and polluted waters of the Dnieper River in the water area of the north-western part of the Black Sea (NWBS) are determined and analyzed. This pollution of the marine environment resulted from the extremely large man-made flood caused by the destruction of the Kakhovka Reservoir dam by the Russian occupation forces in June 2023. Various types of pollutants, which contained in the Kakhovka Reservoir, in the lower Dnieper River and in the bottom sediments, were carried out into the sea by the flood discharge. A significant amount of pollutants was washed into the sea from the flooded areas of the lower Dnieper basin. Special attention is given to changes in oceanographic conditions in the Dnieper-Bug estuary region (DBR) of the Black Sea, which determine the extent of marine pollution. Sea water salinity was used as an indicator of the degree of infiltration of polluted seawater from the estuary into the sea. It was identified that the least transformed river waters with low salinity correspond to high pollution levels. Delft3D-Flow Flexible Mesh hydrodynamic model, developed by Deltares, was used in this study. The simulation was performed for the time period 01 to 30 June 2023. During the first days after the dam break the desalinated and polluted transitional waters, flowing from the Dnieper-Bug estuary, spread over the entire water column down to 15-20 m depth in the DBR area. In this case, the spreading of transitional waters was atypical, as it usually occurs only within the thin near-surface layer of the sea. By mid-June the flow of transitional water with a salinity of 4-6 ‰ covered the entire water column from surface to bottom near the northern coast of the NWBS and over the Odessa Bank. This was facilitated by the supermassive volumes of transformed Dnieper river waters flowing out over a short period of time due to dam break, as well as by the resulting intense gradient currents and mixing caused by spatial gradients in the velocity of the generated currents. The width of the desalinated water plume and its area in the bottom layer were smaller than in the surface layer. The largest vertical salinity gradients were formed to the south of the Odessa Bank, where salinity was 4-6 ‰ in the surface layer and 14-17 ‰ in the bottom. Dilution of the desalinated water plume in the bottom water layer was much more intensive than in the surface water. The character of water circulation pattern in the northern part of the NWBS in the first days after the inflow of transient waters from the Dnieper-Bug estuary was determined by the runoff gradient currents, which later were replaced by wind-driven currents in the Odessa area and over the Odessa Bank, and by the meandering flow of gradient density currents further south. The findings of this study contribute to a better understanding of the impact of the inflow of transformed river waters from the Dnieper-Bug estuary on the variability of oceanographic characteristics and the ecological state of marine waters in the Dnieper-Bug region of the northwestern Black Sea.

2D discrete element analysis of the footing above excavated circle in soil

The bearing capacity and settlements of surface foundations located on a soil slope are the important issues that have to be considered by geotechnical engineers for the design. The presence of an underground void beneath the footing can affect the foundation stability and can lead to serious structure damages. In this study, the results of two-dimensional (2D) discrete element (DE) and finite element (FE) analyses of a surface footing on a soil slope above a void are presented. To validate the numerical model results, the DE results obtained have been compared with experimental test presented in the previous study. After validation of the DE numerical model, parametric studies were carried out to evaluate the effect of important factors on the surface footing performance. The studied parameters include the horizontal spacing of the void axis relative to the slope edge (SH), the vertical spacing of the void crown relative to the footing base (SV), the horizontal spacing of the footing edge relative to the slope edge (De) and the void diameter (Dv). The effects of these parameters on the pressure-settlement curves and the contact force distributions in the soil slope are presented and discussed. The results showed that the footing bearing pressure increases with an increase of SH, SV and De but decreases when Dv increases. The behavior of a surface footing on a soil slope above a void significantly depends on the SV value.

Simulation of intra-chamber processes in a low-thrust rocket engine with a hydrogen-air mixture and counterflow cooling

The principles and methods of improvement of thrust, efficiency and cooling of rocket engines for space flight safety are of great interest for aerospace industry. The development of reliable and durable low-thrust rocket engines for space missions seems to be one of the important areas of development of the rocket and space industry. Issues related to the implementation of a new direction related to the design of propulsion systems for orientation systems of upper stages of launch vehicles using environmentally friendly fuel components are considered. A mathematical model is developed and numerical modelling of processes in the combustion chamber of a low-thrust rocket engine with gaseous fuel components (hydrogen and oxygen) is carried out. The model takes into account the presence of a cooling jacket. The mathematical model includes the effects of turbulence, combustion of a gaseous fuel mixture, kinetics of hydrogen combustion, and radiation. As a result of calculations, distributions of flow quantities in the combustion chamber and thermal loads on its walls are obtained. The influence of hydrogen mass flow rate on the efficiency of the working process and the dependence of thrust on hydrogen mass flow rate are discussed. The results of numerical modelling are compared with the data of a physical experiment obtained through fire tests of an engine model made of stainless steel on a three-dimensional printer.

Numerical and Theoretical Study on Flexural Performance and Reasonable Structural Parameters of New Steel Grating–UHPFRC Composite Bridge Deck in Negative Moment Zone

As the bridge’s structural component is directly subjected to vehicle loads, the stress performance of the bridge deck has a significant impact on the safety, durability, and driving comfort of the bridge. In order to improve the bending performance of the bridge deck in the negative moment zone, a new type of steel grating–UHPFRC composite bridge deck was proposed in this paper. Firstly, structural details and advantages of the new steel grating-UHPFRC composite bridge deck were introduced. Secondly, the finite element program ABAQUS was used to establish a refined solid finite element model of the new bridge deck. The mathematical program MATLAB (PYTHON) was also used to analyze the effects of the structural parameters on bending bearing capacity and put forward reasonable structural parameters of the new bridge deck, considering the technical and economic indexes. Thirdly, the simplified plasticity theory was applied to analyze the bending bearing capacity of the new bridge deck, and the corresponding formula for bending bearing capacity calculation was derived and verified by numerical model results. In addition, the cost–benefit analysis and environmental impact assessment of the new bridge deck were also conducted. The results show that the bending bearing capacity of the new bridge deck in the negative moment zone increases with the increase of the width of the bridge deck, the thickness of the wing plate, and the height of the web plate, with a trend of increasing and then decreasing when the horizontal inclination of the web plate decreases. The bridge deck width does not have a significant effect on improving the bearing capacity. The bearing capacity calculated by theoretical formulas is close to that calculated by numerical models and the maximum relative deviation is 9.1%. The new steel grating-UHPFRC composite bridge deck proposed in this paper is superior to conventional steel-UHPC composite bridge deck in terms of cost-benefit and environmental impact.

Sensitivity of Frequency Domain Near Infrared Spectroscopy for Neurovascular Structure Detection in Biotissue Volume: Numerical Modeling Results.

Through numerical modeling, it has been determined that near infrared spectroscopy with a frequency domain approach can detect neurovascular structures with diameters from 0.5 mm at source-detector distances of 5-8 mm, depending on optical parameters and technical implementation of the method. Among the five classical machine learning methods considered, quadratic discriminant analysis is the most effective for detection. Furthermore, it has been demonstrated that the use of a photomultiplier tube and the registration of both amplitude and phase signal components exhibit the highest sensitivity. Spectroscopy can rival modern ultrasound for detecting arterial vessels. A cross-shaped probe configuration improves sensitivity, and the ratio of reduced scattering coefficient values at different wavelengths is informative for blood-filled vessel detection. These findings are consistent with and significantly extend previous experimental in vivo and in situ studies and could be valuable for intraoperative diagnostic tasks, particularly in neurosurgery.

Comparative analysis of the impedances of key radio engineering devices with PWM based on a linear model and a closed loop method

The problem of comparative analysis of the impedances of switching mode devices with negative feedback by two methods is considered. The use of a linear injection model is justified, on the basis of which analytical expressions of impedances in Laplace form are obtained. The results of numerical modeling and spectral analysis of processes in a closed dynamic system with PWM are presented, confirming the correctness of the analytical model in the frequency range not exceeding half the clock frequency.

Impact Craters on Earth with a Diameter of More than 200 km: Numerical Modeling

The three largest impact craters, the remains of which have been found on Earth to date, had diameters of about 200 km immediately after formation. The search for traces of larger impact structures continues. This paper presents the results of numerical modeling of the formation of terrestrial impact craters larger than those already found. It is shown that the inferred geothermal gradient significantly influences the initial geometry of the impact melt region, which may facilitate the search for the remains of deeply eroded ancient impact structures.