Geometric Model Research Articles

On tiling spherical triangles into quadratic subpatches

Various interpolation and approximation methods arising in several practical applications in geometric modeling deal, at a particular step, with the problem of computing suitable rational patches (of low degree) on the unit sphere. Therefore, we are concerned with the construction of a system of spherical triangular patches with prescribed vertices that globally meet along common boundaries. In particular, we investigate various possibilities for tiling a given spherical triangular patch into quadratically parametrizable subpatches. We revisit the condition that the existence of a quadratic parameterization of a spherical triangle is equivalent to the sum of the interior angles of the triangle being π, and then circumvent this limitation by studying alternative scenarios and present constructions of spherical macro-elements of the lowest possible degree. Applications of our method include algorithms relying on the construction of (interpolation) surfaces from prescribed rational normal vector fields.

UAV vision-based crack quantification and visualization of bridges: system design and engineering application

Accurately measuring visible cracks in bridges is crucial for their structural health diagnosis, damage detection, performance evaluation, and maintenance planning. The primary means of visual crack detection still relies heavily on manual visual inspection, an inefficient process that can pose significant safety risks. This article develops a unmanned aerial vehicle (UAV) vision-based surface crack measurement methodology and visualization scheme for the bridges that can detect and measure cracks automatically with improved efficiency. The surface crack measurement methodology is achieved by designing a three-stage crack sensing system including the You Only Look Once-based crack recognition, U-shaped network-based crack segmentation, and deep-vision-based crack width calculation. This workflow is integrated into a comprehensive UAV inspection system, which is intended for operation at the field. The surface crack visualization scheme is accomplished by taking advantage of time-series image fusion, GPS information migration, and three-dimensional (3D) point cloud technique to reconstruct the 3D geometrical model of the tested bridge, which is convenient for unveiling the crack information in the bridge. The proposed methodology was successfully validated by a case study on an arch bridge. The achievement of this article promotes the UAV vision-based bridge’s surface crack inspection technology to a new status that no preparation for pasting calibration marker is needed, and crack identification, segmentation, and width calculation are realized promptly during the UAV flying on-site, as well as damage evaluation for bridges is visually fulfilled based on the reconstructed digital-graphical 3D model. The working environments and influencing factors to the developed system are sufficiently discussed. Certain limitations in the current application are pointed out for future improvements.

Utilizing Artificial Neural Networks for Geometric Bone Model Reconstruction in Mandibular Prognathism Patients

Patient-specific 3D models of the human mandible are finding increasing utility in medical fields such as oral and maxillofacial surgery, orthodontics, dentistry, and forensic sciences. The efficient creation of personalized 3D bone models poses a key challenge in these applications. Existing solutions often rely on 3D statistical models of human bone, offering advantages in rapid bone geometry adaptation and flexibility by capturing a range of anatomical variations, but also a disadvantage in terms of reduced precision in representing specific shapes. Considering this, the proposed parametric model allows for precise manipulation using morphometric parameters acquired from medical images. This paper highlights the significance of employing the parametric model in the creation of a personalized bone model, exemplified through a case study targeting mandibular prognathism reconstruction. A personalized model is described as 3D point cloud determined through the utilization of series of parametric functions, determined by the application of geometrical morphometrics, morphology properties, and artificial neural networks in the input dataset of human mandible samples. With 95.05% of the personalized model’s surface area displaying deviations within −1.00–1.00 mm relative to the input polygonal model, and a maximum deviation of 2.52 mm, this research accentuates the benefits of the parametric approach, particularly in the preoperative planning of mandibular deformity surgeries.

Deep Blind Fault Activity—A Fault Model of Strong Mw 5.5 Earthquake Seismogenic Structures in North China

North China is one of the high-risk areas for destructive and strong earthquakes in mainland China and has experienced numerous strong historical earthquakes. An earthquake of magnitude MW 5.5 struck Pingyuan County, Dezhou city, in Shandong Province, China, on 6 August 2023. This earthquake was the strongest in the eastern North China Craton since the 1976 Tangshan earthquake. Since the earthquake did not produce surface ruptures, the seismogenic structure for fault responsible for the Pingyuan MW 5.5 earthquake is still unclear. To reveal the subsurface geological structure near the earthquake epicenter, this study used high-resolution two-dimensional (2D) seismic reflection profiles and constructed a three-dimensional (3D) geometric model of the Tuqiao Fault by interpreting the faults in the seismic reflection profiles. This study further combined focal mechanism solutions, aftershock clusters, and other seismological data to discuss the seismogenic fault of the Pingyuan MW 5.5 earthquake. The results show that the Tuqiao Fault is not the seismogenic fault of the MW 5.5 earthquake. The actual seismogenic structure may be related to the NE-oriented high-angle strike-slip blind fault developed in the basement. We further propose three possible fault models for the strong seismogenic structure in North China to discuss the potential seismotectonics in this region.

Geometric Error Identification of Gantry-Type CNC Machine Tool Based on Multi-Station Synchronization Laser Tracers

Laser tracers are a three-dimensional coordinate measurement system that are widely used in industrial measurement. We propose a geometric error identification method based on multi-station synchronization laser tracers to enable the rapid and high-precision measurement of geometric errors for gantry-type computer numerical control (CNC) machine tools. This method also improves on the existing measurement efficiency issues in the single-base station measurement method and multi-base station time-sharing measurement method. We consider a three-axis gantry-type CNC machine tool, and the geometric error mathematical model is derived and established based on the combination of screw theory and a topological analysis of the machine kinematic chain. The four-station laser tracers position and measurement points are realized based on the multi-point positioning principle. A self-calibration algorithm is proposed for the coordinate calibration process of a laser tracer using the Levenberg–Marquardt nonlinear least squares method, and the geometric error is solved using Taylor’s first-order linearization iteration. The experimental results show that the geometric error calculated based on this modeling method is comparable to the results from the Etalon laser tracer. For a volume of 800 mm × 1000 mm × 350 mm, the maximum differences of the linear, angular, and spatial position errors were 2.0 μm, 2.7 μrad, and 12.0 μm, respectively, which verifies the accuracy of the proposed algorithm. This research proposes a modeling method for the precise measurement of errors in machine tools, and the applied nature of this study also makes it relevant both to researchers and those in the industrial sector.

A Calculation Method for the Hyperspectral Imaging of Targets Utilizing a Ray-Tracing Algorithm

This paper proposes a hyperspectral imaging simulation method based on a ray-tracing algorithm. The algorithm combines calculations based on solar and atmospheric visible light radiation as well as the spectral bidirectional reflection distribution function (BRDF) of the target surface material and can create its own scenarios for simulation calculations on demand. Considering the presence of multiple scattering between the target and background, using the ray-tracing algorithm enables the precise computation of results involving multiple scattering. To validate the accuracy of the algorithm, we compared the simulated results with the theoretical values of the visible light scattering intensity from a Lambertian sphere. The relative error obtained was 0.8%. Subsequently, a complex scene of engineering vehicles and grass was established. The results of different observation angles and different coating materials were calculated and analyzed. In summary, the algorithm presented in this paper has the following advantages. Firstly, it is applicable to geometric models composed of any triangular mesh elements and accurately computes the effects of multiple scattering. Secondly, the algorithm combines the spectral BRDF information of materials and improves the efficiency of multiple scattering calculations using nonuniform sampling. The computed hyperspectral scattering data can be applied to simulate airborne or space-borne remote sensing data.

Optimizing Muscle-Tendon Lengths in Reverse Total Shoulder Arthroplasty: Evaluation of Surgical and Implant-Design-Related Parameters.

Optimizing the function of muscles that cross the glenohumeral articulation in reverse total shoulder arthroplasty (RTSA) is controversial. The current study used a geometric model of the shoulder to systematically examine surgical placement and implant-design parameters to determine which RTSA configuration most closely reproduces native muscle-tendon lengths of the deltoid and rotator cuff. A geometric model of the glenohumeral joint was developed and adjusted to represent small, medium, and large shoulders. Muscle-tendon lengths were assessed for the anterior deltoid, middle deltoid, posterior deltoid, and supraspinatus from 0 to 90° of scaption; for the subscapularis from 0° to 60° of internal rotation (IR) and 0° to 60° of scaption; for the infraspinatus from 0° to 60° of external rotation (ER) and 0° to 60° of scaption; and for the teres minor from 0° to 60° of ER at 90° of scaption. RTSA designs were virtually implanted using the following parameters: (1) surgical placement with a centered or inferior glenosphere position and a humeral offset of 0, 5, or 10 mm relative to the anatomic neck plane, (2) implant design involving a glenosphere size of 30, 36, or 42 mm, glenosphere lateralization of 0, 5, or 10 mm, and humeral neck-shaft angle of 135°, 145°, or 155°. Thus, 486 RTSA-shoulder size combinations were analyzed. Linear regression assessed the strength of association between parameters and the change in each muscle-tendon length from the native length. The configuration that most closely restored anatomic muscle-tendon lengths in a small shoulder was a 30-mm glenosphere with a centered position, 5 mm of glenoid lateralization, 0 mm of humeral offset, and a 135° neck-shaft angle. For a medium shoulder, the corresponding combination was 36 mm, centered, 5 mm, 0 mm, and 135°. For a large shoulder, it was 30 mm, centered, 10 mm, 0 mm, and 135°. The most important implant-design parameter associated with restoration of native muscle-tendon lengths was the neck-shaft angle, with a 135° neck-shaft angle being favored (β = 0.568 to 0.657, p < 0.001). The most important surgical parameter associated with restoration of native muscle-tendon lengths was humeral offset, with a humeral socket placed at the anatomic neck plane being favored (β = 0.441 to 0.535, p < 0.001). A combination of a smaller, lateralized glenosphere, a humeral socket placed at the anatomic neck plane, and an anatomic 135° neck-shaft angle best restored native deltoid and rotator cuff muscle-tendon lengths in RTSA. This study of surgical and implant factors in RTSA highlighted optimal configurations for restoration of native muscle-tendon lengths of the deltoid and rotator cuff, which has direct implications for surgical technique and implant selection. Additionally, it demonstrated the most influential surgical and implant factors with respect to muscle-tendon lengths, which can be used to aid intraoperative decision-making.

Inverse optical scatterometry using sketch-guided deep learning

Optical scatterometry, also referred to as optical critical dimension (OCD) metrology, is a widely used technique for characterizing nanostructures in semiconductor industry. As a model-based optical metrology, the measurement in optical scatterometry is not straightforward but involves solving a complicated inverse problem. So far, the methods for solving the inverse scattering problem, whether traditional or deep-learning-based, necessitate a predefined geometric model, but they are also constrained by this model with poor applicability. Here, we demonstrate a sketch-guided neural network (SGNN) for nanostructure reconstruction in optical scatterometry. By learning from training data based on the designed generic profile model, the neural network acquires not only scattering knowledge but also sketching techniques, that allows it to draw the profiles corresponding to the input optical signature, regardless of whether the sample structure is the same as the generic profile model or not. The accuracy and strong generalizability of proposed approach is validated by using a series of one-dimensional gratings. Experiments have also demonstrated that it is comparable to nonlinear regression methods and outperforms traditional deep learning methods. To our best knowledge, this is the first time that the concept of sketching has been introduced into deep learning for solving the inverse scattering problem. We believe that our method will provide a novel solution for semiconductor metrology, enabling fast and accurate reconstruction of nanostructures.

A fast scene geometric modeling approach for digital twins combining neural rendering and model retrieval

ABSTRACT Digital twin technology is a crucial driving technology for the realization of Industry 4.0, which enables simulation, analysis, and prediction by constructing geometric scenes that correspond to physical scenes. Current automation in scene modeling is limited, relying heavily on costly and inefficient manual operations. This hinders digital twin technology progress. This paper proposes a framework for a fast scene geometric modeling method that combines neural rendering and model retrieval for digital twins. Neural rendering techniques first train the collected image data, and the point cloud data of physical entities are rendered. Then based on the semantic mapping algorithm between the point cloud data of physical entities and the multi-view image data of 3D CAD models proposed in this study, the corresponding geometric models are retrieved from the geometric asset library by inputting the point cloud data, and all the retrieved geometric models are embedded in the geometric scene to complete the scene geometric modeling. Finally, in a case study of the digital twin-based scene construction for decommissioned lithium battery dismantling, the effectiveness of this method is demonstrated, which can improve the speed and reduce the cost of scene geometric modeling.

A distributed geometric rewiring model

We propose a distributed rewiring model which starts with a planar graph embedded into the Euclidean space and then behaves as a distributed system, where each node is provided with a set of dynamic links. The proposed rewiring evolves through cycles, where nodes explore the network to identify possible shortcuts and rewire their dynamic links. The rewiring decisions are subject to Euclidean and geodesic distance constrains. The emerging networks were assessed through topological and robustness analyses. We found that the networks display a variety of characteristics observed in complex networks encompassing phenomena such as preferential attachment, the distinctive traits of small-world networks, the presence of community structures, and robustness against degradation process. We consider that our proposal can be applied in the design of those self-managed systems in which there is a limitation on communication resources that can be represented by the Euclidean distance and, however, the components themselves can deploy strategies to optimize the transport of information and develop tolerance before contingencies.

Optimizing the Thickness Uniformity of Magnetron Sputtering Deposited Films on a Large-Scale Curved Workpiece Surface

This study prepared magnetron-sputtering deposition film on a large-scale curved workpiece with a favorable thickness uniformity using a rectangular target cathode. For a substrate rotation structure with eccentric rotation/revolution composite motion, the geometrical model of the thickness distribution of the coating film on a large substrate was established in combination with the principle of rectangular cathode magnetron sputtering. The integral formula of the film thickness was derived, and various film thickness distribution patterns under different parameters were simulated. According to the present research results, the film thickness distribution can be optimized by adjusting the revolution-to-rotation radius and eccentric ratios. At an eccentric distance of 400 mm and a radius ratio of 2.78, the most favorable uniformity degree of 0.87 was obtained, which was controlled by the motion path of the fixed point on the substrate. The proposed approach helps to obtain large-scale films with favorable uniformity on the substrate.

Impact of the place of lightwave in secondary school curricula on conceptual understanding of their basic concepts: a moroccan case study

We aim to highlight the importance of using Huygens’ geometric model when instructing high school students about light. To achieve this, we conducted a research study on teaching lightwaves in Morocco’s secondary schools. We investigated how the place of the lightwave in this teaching affects understanding its basic concepts. We carried out a two-dimensional analysis to examine the place of the lightwave in the high school physics curriculum. We administered a test to assess the conceptual understanding of secondary school and university students. The analysis of the curricula revealed several shortcomings in the teaching of geometrical and wave optics. Furthermore, the assessment of conceptual understanding revealed that a significant proportion of students have misconceptions about the wave concepts of lightwaves and that these misconceptions persist throughout their higher education.

Economic evaluation of mineral deposit: A bottom-up approach

Proper economic evaluation of a mineral deposit is critical to effective investment decision-making in a mineral project. However, this often requires detailed mine production scheduling to produce a schedule that reflects the actual cash flow when the project comes on stream. Because of the dexterity required for this task many mine planners and explorationist attempt to use statistical formulas that approximate the mine scheduling and value of the project. The mathematical model developed for production scheduling often produces a constant production rate schedule over an approximate life span of the project. In this paper, we have attempted to apply a bottom-up approach that begins with geometrical modelling and equipment deployment pattern to define the number of equipment required for each sequence of operation based on available workfront in the development of each bench. Then, based on the number of equipment and the production rate at each sequence of operations, a production schedule is developed. This production schedule therefore will reflect annual cash flow since it is based on the sequence of operation. Application of this method on the west pit of Itakpe mine shows a considerable net present value and internal rate of return of the deposit compared with the evaluation made using statistical models. The NPV of the west pit was found to be USD621 million as against USD122.41 million using Nwosu’s formula and USD123.85 million Taylor’s formula. The above value of NPV using the proposed method shows the maximum expected NPV of the mineral project-based technical restrictions. An understanding of this value can guide the mineral property owner in decision-making.

Geometric Error-Based Multi-Source Error Identification and Compensation Strategy for Five-Axis Side Milling

Based on a multi-source error model, this paper discusses the principle of error element identification and uses the mirror bias method to compensate the geometric errors of a process system. Firstly, a nine-line measurement method to determine the geometric error of the linear feed axes of machine tools is introduced, and the geometric error identification model based on the “nine-line method” is established. Then, using a ballbar mounted in the axial, tangential, and radial directions of the machine, the geometric error elements of the rotation axis are identified by three simple measurements in each direction. Subsequently, for the more common flat vise clamping workpiece in actual production, the workpiece position error is identified by using the traditional process of dimensional chain, and the workpiece attitude error is identified by fitting the angle between the positioning plane and the horizontal plane by the least squares method. Finally, based on the tool position points and tool axis vectors obtained from the multi-source error model, the error compensation value is solved using inverse machine tool kinematics to offset the machining error by mirroring the error value of the same size, and based on the “S-shaped specimen” to compensate the processing experiments, after compensation, the processing error is reduced by 30~45%, verifying the effectiveness of the compensation method.

Analytical Solution for Seismic Stability of 3D Rock Slope Reinforced with Prestressed Anchor Cables

Currently, the study of analytical solutions for the seismic stability of slopes under anchorage conditions is one of the hottest subjects in engineering. In this paper, an analytical solution for the seismic safety of the three-dimensional (3D) two-stage rock slope reinforced with prestressed anchor cables governed by the nonlinear Hoek–Brown criterion was deduced, in which the analyses of seismicity were performed by the latest modified pseudo-dynamic method. This method supplements the consideration of the damping effect of the rock medium on seismic waves, which is more in line with the real seismic situation. A mathematical geometric model was developed for calculating the external forces work and internal energy dissipation acting on 3D rotating rock masses reinforced with prestressed anchor cables, in which the seismic work rate was calculated using a new layer-by-layer superposition summation method. The analytical solution of the safety factor could be collated as an explicit function of several variables, and then, the optimal value was obtained by the Particle Swarm Optimization (PSO) algorithm. To corroborate the accuracy of new analytical solutions, the results were contrasted with those of the pertinent literature. The results of the two comparisons were very close. Ultimately, the sensitivity analyses and coupling effects of seismic pseudo-dynamic factors and prestressing anchorage factors were carried out. It was found that even small seismic intensities had a large effect on the stability of rock slopes with developed joints. Increasing the number of steps and prestressing anchors can effectively improve the stability of rock slopes under seismic effects. The conclusions have significant implications for the anchorage design of the 3D two-stage rock slope as seismic events occur.

Tritium management in ITER test blanket systems port cell for maintenance operations

Four Test Blanket Systems (TBS) will be tested in the International Thermonuclear Experimental Reactor equatorial ports #16 and #18 to verify tritium breeding and heat extraction technology. A significant quantity of tritium would be produced in TBM, and partly released into the port cell from the pipework of TBS or other high-temperature components due to its strong mobility and high permeation. The port cell should be accessible during equipment maintenance and human intervention. This work built a multi-dimensional geometric model to characterize HTO transport in the port cell, absorption/desorption, and diffusion in walls and discussed the effect of paint thickness, ventilation rate, source term, and epoxy properties on detritiation efficiency. The results suggest that a 0.1–0.16 mm paint with the lowest HTO solubility is optimal from the compromise between quick cleanup and tritiated waste decommission. A higher ventilation rate could accelerate detritiation while minimizing the radioactive source by a tritium-resisting layer is the most direct method. The optimized design options for reducing the time required to reach 1 DAC in 12 h still need further discussion because of the delayed HTO source from epoxy paint and dead zone of the flow field.

Multiscalar 3D temporal structural characterisation of Smøla island, mid-Norwegian passive margin: an analogue for unravelling the tectonic history of offshore basement highs

Abstract. Smøla island, situated within the mid-Norwegian passive margin, contains crystalline-basement-hosted intricate fracture and fault arrays formed during a polyphase brittle tectonic evolution. Its detailed study may strengthen correlation attempts between the well-exposed onshore domain and the inaccessible offshore domain, further the understanding of the passive margin evolution, and provide useful constraints on petrophysical properties of fractured basement blocks. A combination of geophysical and remote sensing lineament analysis, field mapping, high-resolution drill hole logging, 3D modelling, petrographic and microstructural studies, and fault gouge K–Ar geochronology made it possible to define five deformation episodes (D1 to D5). These episodes occurred between the post-Caledonian evolution of the regional-scale Møre–Trøndelag Fault Complex (MTFC) and the Late Cretaceous and younger crustal extension preceding the final stages of Greenland–Norway break-up. Each reconstructed deformation stage is associated with different structural features, fault and fracture geometries, and kinematic patterns. Synkinematic mineralisations evolved progressively from epidote–prehnite, sericite–chlorite–calcite, chlorite–hematite, hematite–zeolite–calcite, to quartz–calcite. K–Ar geochronology constrains brittle deformation to discrete localisation events spanning from the Carboniferous to the Late Cretaceous. Multiscalar geometrical modelling at scales of 100, 10, and 1 m helps constrain the extent and size of the deformation zones of each deformation episode, with D2 structures exhibiting the greatest strike continuity and D1 features the most localised. Overall, the approach highlighted here is of great utility for unravelling complex brittle tectonic histories within basement volumes. It is also a prerequisite to constrain the dynamic evolution of the petrophysical properties of basement blocks.

A Novel UAV Air-to-Air Channel Model Incorporating the Effect of UAV Vibrations and Diffuse Scattering

In this paper, we propose a geometric channel model for air-to-air (A2A) unmanned aerial vehicle (UAV) communication scenarios. The model is established by incorporating line-of-sight, specular reflection, and diffuse scattering components, and it can capture the impacts of UAV vibrations induced by the propeller’s rotation. Based on UAV heights and ground scatterer density, a closed-form expression is derived to jointly capture the zenith and azimuth angular distributions of diffuse rays. The power of diffuse rays is modeled according to the grazing angle of the rays and the electrical properties and roughness of the ground materials. Key statistics, including the temporal autocorrelation function, spatial cross-correlation function, Doppler power spectrum density, and coherence time are derived, providing an in-depth understanding of the time-variant characteristics of the channel. The results indicate that the presented model is capable of capturing certain A2A channel characteristics, which align with the corresponding theoretical analysis. The findings suggest that the scattering effect of the A2A channel is significantly influenced by the altitude of the UAV. Additionally, it is shown that UAV vibrations can introduce extra Doppler frequencies, notably decreasing the temporal correlation and coherence time of the channel. This effect is more prominent when the system operates at high-frequency bands. The effectiveness of the presented model is confirmed through a comparison of its statistics with those of an existing model and with available measurement data.

Evaluation of the Truncated Cone-Rhomboid Pyramid Formula for Simplified Right Ventricular Quantification: A Cardiac Magnetic Resonance Study.

Background/Objective: Cardiac magnetic resonance (CMR) is the reference method for right ventricular (RV) volume and function analysis, but time-consuming manual segmentation and corrections of imperfect automatic segmentations are needed. This study sought to evaluate the applicability of an echocardiographically established truncated cone-rhomboid pyramid formula (CPF) for simplified RV quantification using CMR. Methods: A total of 70 consecutive patients assigned to RV analysis using CMR were included. As standard method, the manual contouring of RV-short axis planes was performed for the measurement of end-diastolic volume (EDV) and end-systolic volume (ESV). Additionally, two linear measurements in four-chamber views were obtained in systole and diastole: basal diameters at the level of tricuspid valve (Dd and Ds) and baso-apical lengths from the center of tricuspid valve to the RV apex (Ld and Ls) were measured for the calculation of RV-EDV = 1.21 × Dd2 × Ld and RV-ESV = 1.21 × Ds 2 × Ls using CPF. Results: RV volumes using CPF were slightly higher than those using standard CMR analysis (RV-EDV index: 86.2 ± 29.4 mL/m2 and RV-ESV index: 51.5 ± 22.5 mL/m2 vs. RV-EDV index: 81.7 ± 24.1 mL/m2 and RV-ESV index: 44.5 ± 23.2 mL/m2) and RV-EF was lower (RV-EF: 41.1 ± 13.5% vs. 48.4 ± 13.7%). Both methods had a strong correlation of RV volumes (ΔRV-EDV index = -4.5 ± 19.0 mL/m2; r = 0.765, p < 0.0001; ΔRV-ESV index = -7.0 ± 14.4 mL/m2; r = 0.801, p < 0.0001). Conclusions: Calculations of RV volumes and function using CPF assuming the geometrical model of a truncated cone-rhomboid pyramid anatomy of RV is feasible, with a strong correlation to measurements using standard CMR analysis, and only two systolic and diastolic linear measurements in four-chamber views are needed.

Simulation and Experimental Verification of Dispersion and Explosion of Hydrogen–Methane Mixture in a Domestic Kitchen

Hydrogen is a carbon-free energy source that can be obtained from various sources. The blending of hydrogen represents a transitional phase in the shift from natural gas systems to hydrogen-based systems. However, concerns about the safety implications of introducing hydrogen have led to extensive discussions. This paper utilizes Fluent 17.0 numerical simulation software to simulate the leakage of hydrogen-blended natural gas in a closed domestic kitchen and analyze the concentration distribution and its variation pattern after a leakage. An experimental platform is set up, and a mixture of nitrogen and helium gas is used as a substitute for hydrogen-blended natural gas for the simulations and experiments. The simulation results demonstrate that the leaked gas spreads and accumulates towards the top of the space, gradually filling the entire area as the leak persists. As the hydrogen content in the leaked gas increases, the dispersion capacity of the gas in the confined space also increases. Furthermore, as the flow rate of the leaked gas increases, the average concentration of the leaked gas rises, and the gas stratification in the confined kitchen diminishes. The concentration distribution observed in the experiments aligns with the simulation results. After establishing the feasibility conditions of the model, the dispersion of the hydrogen-blended natural gas in the kitchen is further simulated. The results suggest that blending hydrogen into the gas enhances the dispersion of the gas after a leak, leading to a wider distribution within the kitchen and an increased risk in the event of a leak. Additionally, this paper employs the CASD module of FLACS 11.0 software to construct a three-dimensional geometric model of the domestic kitchen for simulation studies on the explosion of hydrogen-blended natural gas in a confined space. By adjusting the hydrogen ratio in the combustible gases present in the space and examining the variations in hydrogen concentration and external conditions, such as opening or closing the door, the influence on parameters including the peak explosion pressure, explosion overpressure, explosion flame temperature, and explosion response time are examined. Furthermore, the extent of the explosion area is determined, and the effect of hydrogen on the blast is clarified.