Simplified Model Research Articles

Ultraviolet Radiation Fields in Star-forming Disk Galaxies: Numerical Simulations with TIGRESS-NCR

With numerical simulations that employ adaptive ray-tracing (ART) for radiative transfer at the same time as evolving gas magnetohydrodynamics, thermodynamics, and photochemistry, it is possible to obtain a high-resolution view of ultraviolet (UV) fields and their effects in realistic models of the multiphase interstellar medium. Here, we analyze results from TIGRESS-NCR simulations, which follow both far-UV (FUV) wavelengths, important for photoelectric heating and polycyclic aromatic hydrocarbon excitation, and the Lyman continuum (LyC), which photoionizes hydrogen. Considering two models, representing solar neighborhood and inner-galaxy conditions, we characterize the spatial distribution and time variation of UV radiation fields, and quantify their correlations with gas. We compare four approximate models for the FUV to simulated values to evaluate alternatives when full ART is infeasible. By convolving FUV radiation with density, we produce mock maps of dust emission. We introduce a method to calibrate mid-IR observations, for example from JWST, to obtain high-resolution gas surface density maps. We then consider the LyC radiation field, finding most of the gas exposed to this radiation to be in ionization–recombination equilibrium and to have a low neutral fraction. Additionally, we characterize the ionization parameter as a function of the environment. Using a simplified model of the LyC radiation field, we produce synthetic maps of emission measure (EM). We show that the simplified model can be used to extract an estimate of the neutral fraction of the photoionized gas and mean free path of ionizing radiation from observed EM maps in galaxies.

Lattice Structures-Mechanical Description with Respect to Additive Manufacturing.

Lattice structures, characterized by their repetitive, interlocking patterns, provide an efficient balance of strength, flexibility, and reduced weight, making them essential in fields such as aerospace and automotive engineering. These structures use minimal material while effectively distributing stress, providing high resilience, energy absorption, and impact resistance. Composed of unit cells, lattice structures are highly customizable, from simple 2D honeycomb designs to complex 3D TPMS forms, and they adapt well to additive manufacturing, which minimizes material waste and production costs. In compression tests, lattice structures maintain stiffness even when filled with powder, suggesting minimal effect from the filler material. This paper shows the principles of creating finite element simulations with 3D-printed specimens and with usage of the lattice structure. The comparing of simulation and real testing is also shown in this research. The efficiency in material and energy use underscores the ecological and economic benefits of lattice-based designs, positioning them as a sustainable choice across multiple industries. This research analyzes three selected structures-solid material, pure latices structure, and boxed lattice structure with internal powder. The experimental findings reveal that the simulation error is less than 8% compared to the real measurement. This error is caused by the simplified material model, which is considering the isotropic behavior of the used material PA12GB (not the anisotropic model). The used and analyzed production method was multi jet fusion.

Modeling Method of Meshing Stiffness of High Speed Spur Gears

As one of the most important dynamic excitation sources, the gear time-varying mesh stiffness is the key parameter of gear system dynamic model. So how to calculate the gear mesh stiffness accurately is of great importance. Based on Euler beam theory, this paper proposes a primitive calculation algorithm (PCA) to calculate the time-varying mesh stiffness of spur gears. A simplified gear model is developed based on nonlinear strain-displacement relationships, and a finite element program is used to simulate the gear's meshing stiffness. The rationality of this method was verified through finite element analysis. The results indicate that mesh stiffness plays a crucial role in controlling mesh deformation. In view of this, this study lays a theoretical foundation for further analysis of the vibration and noise of gears.

Computational Methods for Verifying the Normative Requirements Regarding the Lateral Correction Force of a Powered Roof Support

The article discusses laboratory methods and the corresponding computational methods verifying compliance with the normative requirements regarding the lateral correction force of the powered roof support. The currently used flat model only allows for checking the normative requirement in relation to the sum of active forces of the correction cylinders installed in the roof support. Determining the required value of the active force of each cylinder is possible due to the simplified FEM model of a powered roof support, described in the research work, treated as a uniform weightless elastic body loaded with a concentrated force recreating the weight of the roof support located on an inclined longwall panel. The third analysed computational method involves determining the reaction in the four correction cylinders of the roof support, creating a spatial, statically indeterminate system of forces. It enables determining the range of variability of the response in the correction cylinders as a function of the distribution of floor pressure on the roof support base. The discussed computational methods were used to determine, for example, the lateral correction force of one of the types of powered supports used in a longwall panel inclined at an angle of 35°. The usefulness of the discussed calculation methods at various stages of the designing process of the powered support and its certification has been confirmed.

Simplified Modelling Techniques for Dynamic Wireless Power Transfer

Simplified Modelling Techniques for Dynamic Wireless Power Transfer

Electrokinetics-Driven Anodic Oxide Pore Formation: Linear and Weakly Nonlinear Analysis

Abstract Anodization of aluminium in an acidic medium facilitates the formation of well-ordered nanoporous anodic oxide films. The mechanism of pore formation is investigated as a morphological instability using a simplified model. The model accounts for the high field conduction law and field-assisted reactions (oxide formation/dissolution) only at the oxide-solution interface. The role of Butler–Volmer electrokinetics, electrolyte pH, anodic efficiency, and interface curvature on reaction kinetics are taken into account. Linear stability analysis suggests that the oxide film is unstable to well-defined wavelengths in specific ranges of parameters such as anodizing efficiency, applied voltage and electrolyte pH. Subsequently, a weakly nonlinear analysis is carried out to determine the nature of bifurcation beyond the stability threshold. Our findings indicate that the instability exhibits a subcritical nature, well in agreement with experimental observations.

Interaction of a pyrene derivative with cationic [60]fullerene in phospholipid membranes and its effects on photodynamic actions.

We have reported that upon visible light irradiation, ferrocene-porphyrin-[60]fullerene triad molecules yield long-lived charge-separated states, enabling the control of the plasma membrane potential (V m) in living cells. These previous studies indicated that the localization of the triad molecules in a specific intra-membrane orientation and the suppression of the photodynamic actions of the [60]fullerene (C60) moiety are likely important to achieve fast and safe control of V m, respectively. In this study, by mimicking our previous system of triad molecules and living cells, we report a simplified model system with a cationic C60 derivative (catC60) and a liposome with embedded 1-pyrenebutyric acid (PyBA) to demonstrate that the addition of PyBA was important to achieve fast and safer control of V m.

Power system frequency nadir prediction based on data-driven and power-frequency polynomial fitting

As the proportion of renewable energy and power electronics equipment continues to rise, the level of rotational inertia decreases considerably, resulting in severe frequency stability challenges to the power grid. It is of great significance to accurately predict the frequency nadir following a large disturbance. This paper proposes a novel data-model fusion-driven approach for the prediction of frequency nadir. As the physics-driven part, a Simplified Prediction Model (SPM) based on power-frequency polynomial fitting is developed to quickly produce the frequency nadir. As the data-driven part, Back Propagation Neural Network (BPNN) is deployed to correct the errors of the SPM to achieve more accurate results. This serial integration scheme not only obtains the final prediction result with higher accuracy, but also meets the computational efficiency requirements of online prediction. Compared with existing integration-driven methods, SPM only focuses on the active power-frequency characteristics of the system, which retains the most critical effects and greatly reduces the dependence of BPNN on sample data quality. Case studies on a modified IEEE 39-bus system verify the effectiveness of the proposed approach.

A simplified theoretical model to predict the shear behavior of concrete/rock joints with Gaussian irregularities

The inherent randomness in the roughness of concrete-rock joints poses challenges for analyzing their shear behaviour. Traditionally, theoretical models simplify these complex geometries by approximating them as regular triangular asperities. This study presents a novel micro-mechanical model developed to simulate the shear behaviour. The model accounts for the complex and random roughness of joint profiles by considering varying asperity heights as Gaussian distribution. Local stress distribution among individual asperities was determined. A failure envelope, mapping joint dilation against asperity height, was obtained, allowing us to identify the critical joint dilation for individual asperities at which individual asperity failure results. The parameter of shear area ratio, which captures the portion of failed asperities relative to all asperities, is introduced, and its semi-analytical solution was derived. Additionally, a wedge-shaped planar failure surface is determined based on laboratory observations and equations to estimate post-failure residual resistance is obtained. A series of direct shear tests on rock/concrete joints with a range of geometries was conducted to validate the proposed model and reasonably close agreement was obtained between predicted and measured performance. An additional case study analyzing the results of field tests on a drilled pile socketed in rock has been conducted using the proposed method. The results showed good agreement with the measured data, demonstrating the validity and practicality of the proposed approach.

TriPoD: Tri-Population size distributions for Dust evolution

Context. Dust coagulation and fragmentation impact the structure and evolution of protoplanetary disks and set the initial conditions for planet formation. Dust grains dominate the opacities, they determine the cooling times of the gas via thermal accommodation in collisions, they influence the ionization state of the gas, and the available grain surface area is an important parameter for the chemistry in protoplanetary disks. Therefore, dust evolution is an effect that should not be ignored in numerical studies of protoplanetary disks. Available dust coagulation models are, however, too computationally expensive to be implemented in large-scale hydrodynamic simulations. This limits detailed numerical studies of protoplanetary disks, including these effects, mostly to one-dimensional models. Aims. We aim to develop a simple – yet accurate – dust coagulation model that can be easily implemented in hydrodynamic simulations of protoplanetary disks. Our model shall not significantly increase the computational cost of simulations and provide information about the local grain size distribution. Methods. The local dust size distributions are assumed to be truncated power laws. Such distributions can be fully characterized by only two dust fluids (large and small grains) and a maximum particle size, truncating the power law. We compare our model to state- of-the-art dust coagulation simulations and calibrate it to achieve a good fit with these sophisticated numerical methods. Results. Running various parameter studies, we achieved a good fit between our simplified three-parameter model and DustPy, a state-of-the-art dust coagulation software. Conclusions. We present TriPoD, a sub-grid dust coagulation model for the PLUTO code. With TriPoD, we can perform twodimensional, vertically integrated dust coagulation simulations on top of a hydrodynamic simulation. Studying the dust distributions in two-dimensional vortices and planet-disk systems is thus made possible.

Correlation between Voltammetry of Immobilized Particles and Mott-Schottky analysis of metal corrosion patinas.

The conjoint application of the voltammetry of immobilized particles (VIMP) methodology and the Mott-Schottky analysis (MS) of impedance data to studying metal corrosion patinas is described. The study is applied to copper and bronze objects exploiting the semiconducting character of cuprite and other copper corrosion products. A simplified theoretical modeling of MS analysis at microparticulate deposits extracted from metal corrosion layers attached to graphite electrodes is provided. The proposed model compensates for the disturbing effect of the regions of the basal electrode directly exposed to the electrolyte. Alternative models accounting for the variation of the density of charge carriers with depth are tested as well as the correlation between VIMP and MS data with reasonably satisfactory results.

Observation of Space-Dependent Rotational Doppler Shifts with a Single Ion Probe.

We present an experiment investigating the rotational Doppler effect using a single trapped ion excited by two copropagating vortex laser beams. The setup isolates the azimuthal gradients of the fields, eliminating longitudinal and curvature effects. We provide a detailed characterization of the phenomenon by deterministically positioning a single ion across the beams and measuring fluorescence spectra with sharp "dark resonances" whose features depend on the angular velocity of the ion and the difference of optical orbital angular momentum between the two beams. The interpretation of the measurements is supported by numerical simulations and by a simplified analytical model. Our results reveal key properties of the rotational Doppler effect, showing that it increases approaching the center of the beam and that it is independent of the waist of the beam. This offers insights into the feasibility of superkicks or super-Doppler shifts for sensing and manipulating atomic motion transverse to the beams' propagation direction.

Explanation of the excesses in associated di-photon production at 152 GeV in 2HDM

Statistically significant excesses exist at around 152 GeV in associated di-photon production (γγ + X) in the sidebands of SM Higgs analyses of ATLAS (using the full run-2 dataset). They are most pronounced in the single-τ, missing-transverse-energy, four-jet and ⩾ 1ℓ+ ⩾ 1b-jet channels (≈ 3σ) and can be explained by the Drell-Yan production of new Higgs bosons, i.e. pp → W*→ H±H0. We first examine the excesses in a simplified model approach, considering that H± decays to τν, WZ or tb. Both the τν and tb decay modes individually lead to a significance of ⪅ 4σ while for WZ one can obtain at most 3.5σ. This is because the decays of WZ lead to multiple leptons contributing to the two-lepton channel which does not show an excess at 152 GeV. Next, we consider two-Higgs-doublet models where the charged Higgs does not decay to WZ at tree-level, finding a significance of ⪆ 4σ for a branching ratio of the new neutral Higgs to photons of ≈2%. Even though this branching fraction is quite sizable, it can be obtained in composite models or via the Lagrangian term \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\lambda }_{6}{H}_{1}^{†}{H}_{1}{H}_{2}^{†}{H}_{1}$$\\end{document} + h.c. breaking the commonly imposed Z2 symmetry.

Search for single vector-like B quark production in hadronic final states at the LHC

In this paper, we study the discovery potential of a Vector-Like B quark (VLB) via the process pp→B(→bZ)j→b(Z→νlνl¯)j\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$pp \\rightarrow B(\\rightarrow bZ)j\\rightarrow b(Z \\rightarrow \ u _l\\bar{\ u _l})j$$\\end{document} at the Large Hadron Collider (LHC) with s=14\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{s}=14$$\\end{document} TeV. In the framework of a simplified model, we perform a scan over its parameter space and test its viability following a Monte Carlo analysis developed to include all production and decay dynamics. We use cut-and-count combined with Extreme Gradient Boosting (XGBoost) methods to classify the signal and background events in order to improve the efficiency of signal identification and background rejection. We find that this approach can reduce background events significantly while the signal retention rate is much higher than that of traditional methods, thereby improving the VLB discovery potential. We then calculate the exclusion and discovery capabilities for VLBs and find that the advantages of the cut-and-count plus XGBoost method especially lie in the high-mass region, i.e., mB>1500GeV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_B > 1500 \ ext { GeV}$$\\end{document}. We finally obtain the following LHC results in terms of the coupling and chiral structure of a singlet heavy VLB interactions: (i) for g∗=0.2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g^{*}=0.2$$\\end{document} and RL=0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$R_L=0$$\\end{document} with 3000 fb-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1}$$\\end{document}, the B quark mass can be be excluded (discovered) up to 3000 GeV (2500 GeV); (ii) for g∗=0.2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g^{*}=0.2$$\\end{document} and RL=0.5\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$R_L=0.5$$\\end{document} with 3000 fb-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1}$$\\end{document}, the exclusion (discovery) region can reach up to 4750 GeV (4250 GeV).

Incorporating soil‐structure interaction into simplified numerical models for fragility analysis of RC structures

AbstractSimplified building models are a valuable option for seismic assessment at the regional scale. These models often use calibrated springs to model column behaviour, and recent advances have made them suitable for capturing torsional response in Reinforced‐Concrete Moment‐Resisting‐Frames. Nevertheless, their validation is typically achieved using fixed‐base models, which do not include the influence of soil‐structure interaction (SSI). This study introduces a novel approach to quantify the accuracy of a recently developed simplified model while accounting for dynamic SSI, using a newly implemented, refined 3D Finite Element non‐linear soil model in OpenSees. The accuracy of the simplified structural model is assessed by comparing the results of non‐linear dynamic analyses with those of a refined model in terms of (i) a peak structural demand parameter such as the interstorey‐drift ratio and (ii) fragility curves computed from cloud analysis and accounting for collapse cases. The study presents details of the proposed refined approach for 3D soil modelling in OpenSees, focusing on implementing free‐field boundary conditions and structure‐to‐soil connections. Results show that the accuracy of the simplified model is maintained, even in the presence of SSI, and it successfully captures the overall structural response measured at peak demand. For the proposed case study, the difference between the simplified and refined models’ fragility curves’ medians is 4% and 2% for fixed and SSI models, respectively. The simplified structural model, combined with the refined soil model for SSI effects, presents an innovative and conservative, yet computationally efficient, alternative for seismic risk analysis, even in the presence of structural irregularity.

A simplified bedside model for acute kidney injury prediction after percutaneous coronary intervention

Abstract Background Concern for acute kidney injury (AKI) after percutaneous coronary intervention (PCI) is a major determinant in the decision to pursue cardiac catheterization if otherwise indicated. A bedside AKI risk prediction model derived from the National Cardiovascular Data Registry (NCDR), most recently updated in 2023, employs relatively complex parameters for AKI prediction. Purpose To test the efficacy of a simplified pre-procedural 5-variable bedside score for post-PCI AKI prediction based on the NCDR model. Methods We analyzed our institution’s NCDR-reported catheterization laboratory registry between Q2 2022 and Q1 2023. Inclusion criteria, variable selection, and definitions were based on the NCDR model. A bedside prediction model was validated, including 5 pre-procedural variables: anemia (hemoglobin <10 mg/dl), history of heart failure, stage 4 or 5 chronic kidney disease (glomerular filtration rate <30 mg/dl), cardiac arrest, and shock status. Results We included 840 patients, of which 593 (70.6%) were men. Mean age was 69 (±12.5). Most procedures were performed on an urgent basis (444, 52.9%). AKI developed in 107 (12.7%) patients. Compared to the non-AKI group, patients who developed AKI were older (age 74±11.5 vs 68.3 ±12.5), and more frequently had anemia (23.6% vs 7.4%), diabetes (51.4% vs 35.2%), severe frailty (45.8% vs 23.7%), heart failure (62.6% vs 29.7%), cardiac arrest (13.1% vs 2.9%), CKD (30.8% vs 6%), shock (12.1% vs 1.5%), and mechanical circulatory support (7.5% vs 2.9%), p<0.001. Utilizing integer scores, the simplified bedside model employing 5 pre-procedural variables demonstrated discrimination similar to the bedside NCDR model, with a C-statistic of 0.77 (0.60-0.88) and excellent calibration (intercept of 0.31 and slope of 1.03) (Figure), indicating the robust performance of the model (compared to bedside NCDR model with a C-statistic 0.793, intercept of -0.197, and slope of 0.957). Conclusion A simplified NCDR-based bedside score using 5 pre-procedural variables is an effective tool for predicting AKI as a complication of PCI, with comparable performance to the 2023 updated bedside NCDR AKI risk model.

Adverse cardiovascular events in patients with cancer: a biomarker-enhanced clinical risk score

Abstract Background/Introduction Continuously improving cancer-specific survival puts a growing proportion of cancer patients at risk of major adverse cardiovascular events (MACE), but reliable tools for cardiovascular risk prediction remain unavailable (1). Purpose To assess a broad panel of cardiovascular biomarkers and clinical risk factors and to provide a novel risk score for the prediction of MACE in cancer patients. Methods Among 2’192 patients with newly diagnosed or recurrent cancer, a broad panel of cardiovascular biomarkers, including NT-proBNP, P-/E-/L-selectins, ICAM-1, VCAM-1, sLOX-1, CRP, and Lp(a) were measured by blinded study personnel. Patients were followed prospectively for the occurrence of 2-year MACE and 5-year cardiovascular death (2,3). Uni- and multivariable competing risk models were fit to assess independent associations of cardiovascular risk factors and biomarkers with MACE, and a risk score was developed. Results Beyond clinical risk factors, ICAM-1, P-/E-/L-selectins and NT-proBNP levels were independently linked to 2-year MACE risk, though shape of associations differed markedly. A clinical risk score was derived, assigning +1 point for male sex, smoking, and age ≥60 years, and +2 points for ASCVD history, yielding a bootstrapped C-statistic of 0.76 (95%CI:0.71-0.81). Implementation of cardiovascular biomarkers conferred improved performance (0.83, 95%CI:0.78-0.88). A simplified model, solely informed by clinical variables and readily available NT-proBNP (+1 point for NT-proBNP≥230pg/mL), achieved similar performance (0.80, 95%CI:0.74-0.86), with a cumulative 2-year MACE incidence of 0% in low- and 11.0% in high-risk patients. Conclusions Cancer patients are at substantial MACE risk. The herein presented first-of-its-kind risk score, integrating traditional cardiovascular risk factors and biomarkers, including NT-proBNP, effectively predicts 2-year MACE and 5-year cardiovascular death.

Seismic Performance of the Bridge Piers Reinforced with PP-ECC in the Plastic Hinge Region

Ductility deficiency in reinforced concrete (RC) piers is a substantial factor contributing to earthquake damage in bridge structures. In this study, the conventional concrete in the potential plastic hinge area was substituted with polypropylene fiber-reinforced engineered cementitious composite (PP-ECC), thereby leveraging its high ductility characteristics to explore the seismic performance of RC piers strengthened by PP-ECC. A finite element model was established to discuss the reinforced height, axial compression ratio, longitudinal reinforcement ratio and stirrup ratio on the hysteresis performance by referring to a published paper with a 1/5 scale test specimen. Based on the results of the parameter analysis, an ideal set of configuration parameters was proposed. Subsequently, a full-size simplified mechanical model was established with ideal reinforcement parameters. Time-history and pushover analysis were utilized to test the seismic response. The study indicates that the two analytical methods were largely consistent. The improvement in the displacement and shear force of the strengthened pier gradually decreased with the increase in acceleration amplitude. Time-history analysis reveals a decrease in the enhancement of displacement from 82.26% to 17.98% and a reduction at the bottom reaction from 12.2% to 1.5%. Under the pushover analysis method, the retrofitting level of the top displacement decreased from 77.22% to 11.53%, whereas the improvement level of the bottom shear decreased from 9.62% to 3.69%. These results provide a theoretical basis and reference standard for the application of PP-ECC.

On the Numerical Investigation of Natural-Convection Heat Sinks Across a Wide Range of Flow and Operating Conditions

Many designs for natural-convection heat sinks and semi-empirical correlations have been proposed in the recent years, but they are only valid in a limited range of Elenbaas numbers El and were mostly tested for laminar flows. To alleviate those limits, parametric studies with 2D and quasi-3D models were carried out, in ranges of Grashof numbers up to 1.55×1011 and Elenbaas numbers up to 3.42×107. Ansys Fluent’s laminar, transition-SST, SST k-ω and k-ϵ models were applied. In addition, when used in this valid range, i.e., mean Elenbaas numbers, with the simplified quasi-3D model, the transition-SST model could predict better results, overestimating the heat flux by 10 to 15% compared to semi-empirical correlations. The 2D model was not deemed satisfying, regarding turbulence models. Consequently, a quasi-3D model was developed: it appeared to be an efficient trade-off between computational time and prediction accuracy, in particular for turbulence models. New grouping factors were also found, to ensure proper dimensioning of natural-convection heat sinks. They corresponded to non-dimensional parameters that dictated the physical behaviour of the heat sink with respect to the semi-empirical correlations. Typically, the ratio of the spacing to the optimal spacing predicted by Bar-Cohen’s correlation turned out to be an appropriate grouping factor with a threshold of 1, above which the fins could safely be considered as isolated, thus greatly simplifying all further calculations.

Study and Analysis of the Thunder Source Location Error Based on Acoustic Ray-Tracing

Error analysis and estimation of thunder source location results is a prerequisite for obtaining accurate location results of thunder sources, which is of great significance for a deeper understanding of the physical process of lightning channel discharges. Most of the thunder source location algorithms are based on the simplified model of the straight-line propagation of acoustic waves to determine the location of the thunder source; however, the acoustic wave is affected by the inhomogeneity of the atmosphere medium in the propagation process and its acoustic ray will be bent. Temperature and humidity are the main factors affecting the vertical distribution of the velocity of sound in the atmosphere, therefore, it is necessary to study the changes in location errors under the models of uniform vertical distribution of temperature only and uniform vertical distribution of humidity only. This paper focuses on the theory of acoustic ray-tracing in neglecting the presence of the wind and the acoustic attenuation and the theoretical derivation of the location error of thunder source inversion for the three models is carried out by using MATLAB R2019b programming. Then, simulation analysis and comparative study on the variation law of thunder source location error with the height of the source, ground temperature, ground humidity, and array position under the three models are carried out. The results of the study show that the maximum location error can be obtained from the straight-line propagation model, the location error obtained from the model of uniform vertical distribution of temperature only is the second, and the location error obtained from the model of uniform vertical distribution of humidity only is the least and can be negligible compared to the first two models. In the trend of error variation, the variation of location error with temperature and humidity is relatively flat in the first two models; however, the variation of location error with the height of the thunder source is more drastic, which can be more than 80%. The location error obtained from the array inversion closer to the thunder source increases linearly with the height of the thunder source, the location error obtained from the more distant array inversion shows a fast-decreasing trend at the height of the thunder source from 500 to 3500 m, and a flat trend above 3500 m. The location error varies relatively smoothly with the height of the thunder source, the ground temperature, and the ground humidity in the model of uniform vertical distribution of humidity only. In addition, the position of the array also has an important effect on lightning location. The further the horizontal distance from the source, the greater the location error will be obtained in the first two models, and when the thunder source is at a low height and detected at a long distance, the location error will be very large, so relevant data should be modified in actual observation.