[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (12,694)

Search Parameters:
Keywords = variational parameter

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
25 pages, 4475 KiB  
Article
Characterizing the Temporal Variation of Airborne Particulate Matter in an Urban Area Using Variograms
by Gokul Balagopal, Lakitha Wijeratne, John Waczak, Prabuddha Hathurusinghe, Mazhar Iqbal, Rittik Patra, Adam Aker, Seth Lee, Vardhan Agnihotri, Christopher Simmons and David J. Lary
Air 2025, 3(1), 7; https://doi.org/10.3390/air3010007 (registering DOI) - 3 Mar 2025
Abstract
This study aims to determine the optimal frequency for monitoring airborne pollutants in densely populated urban areas to effectively capture their temporal variations. While environmental organizations worldwide typically update air quality data hourly, there is no global consensus on the ideal monitoring frequency [...] Read more.
This study aims to determine the optimal frequency for monitoring airborne pollutants in densely populated urban areas to effectively capture their temporal variations. While environmental organizations worldwide typically update air quality data hourly, there is no global consensus on the ideal monitoring frequency to adequately resolve pollutant (particulate matter) time series. By applying temporal variogram analysis to particulate matter (PM) data over time, we identified specific measurement intervals that accurately reflect fluctuations in pollution levels. Using January 2023 air quality data from the Joppa neighborhood of Dallas, Texas, USA, temporal variogram analysis was conducted on three distinct days with varying PM2.5 (particulate matter of size ≤ 2.5 μm in diameter) pollution levels. For the most polluted day, the optimal sampling interval for PM2.5 was determined to be 12.25 s. This analysis shows that highly polluted days are associated with shorter sampling intervals, highlighting the need for highly granular observations to accurately capture variations in PM levels. Using the variogram analysis results from the most polluted day, we trained machine learning models that can predict the sampling time using meteorological parameters. Feature importance analysis revealed that humidity, temperature, and wind speed could significantly impact the measurement time for PM2.5. The study also extends to the other size fractions measured by the air quality monitor. Our findings highlight how local conditions influence the frequency required to reliably track changes in air quality. Full article
Show Figures

Figure 1

Figure 1
<p>Sensors used in the MINTS Central Node air module.</p>
Full article ">Figure 2
<p>Variation in the daily average <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>2.5</mn> </mrow> </msub> </mrow> </semantics></math> concentration recorded in the Joppa neighborhood of Dallas, Texas, during January 2023. The calendar format displays the daily average concentrations using a gradient of color tones, where darker shades indicate higher pollution levels, and lighter shades represent lower levels. The color bar below the calendar provides a reference scale for the concentrations in <math display="inline"><semantics> <mi mathvariant="sans-serif">μ</mi> </semantics></math><math display="inline"><semantics> <mrow> <msup> <mi>g/m</mi> <mn>3</mn> </msup> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Time series plots of <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>2.5</mn> </mrow> </msub> </mrow> </semantics></math> concentrations for three days in January 2023: the most polluted day (2 January 2023), the intermediate polluted day (29 January 2023), and the least polluted day (24 January 2023) in Joppa, Dallas, Texas. (<b>a</b>,<b>c</b>,<b>e</b>) Hourly <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>2.5</mn> </mrow> </msub> </mrow> </semantics></math> concentrations for these days. (<b>b</b>,<b>d</b>,<b>f</b>) The <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>2.5</mn> </mrow> </msub> </mrow> </semantics></math> concentrations at one-second intervals. The horizontal dashed lines indicate air quality standards set by the EPA (<math display="inline"><semantics> <mrow> <mn>35</mn> <mo> </mo> <mi mathvariant="sans-serif">μ</mi> <msup> <mi>g/m</mi> <mn>3</mn> </msup> </mrow> </semantics></math>, indicated by a red dashed line), the EU (<math display="inline"><semantics> <mrow> <mn>25</mn> <mo> </mo> <mi mathvariant="sans-serif">μ</mi> <msup> <mi>g/m</mi> <mn>3</mn> </msup> </mrow> </semantics></math>, indicated by a blue dashed line), and the WHO (<math display="inline"><semantics> <mrow> <mn>15</mn> <mo> </mo> <mi mathvariant="sans-serif">μ</mi> <msup> <mi>g/m</mi> <mn>3</mn> </msup> </mrow> </semantics></math>, indicated by a green dashed line).</p>
Full article ">Figure 4
<p>Variogram plots of <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>2.5</mn> </mrow> </msub> </mrow> </semantics></math> concentration in decreasing order of pollution (2 January 2023, 29 January 2023, and 24 January 2023). (<b>a</b>,<b>c</b>,<b>e</b>) Variograms for Window Number 10 on the respective days. (<b>b</b>,<b>d</b>,<b>f</b>) Window Number 50,000 on the same days. Window Numbers 10 and 50,000 were randomly selected from a total of 85,500 windows for each day. The theoretical variogram model is represented by yellow dashed lines, while the empirical model is depicted in blue. The cross indicates the coordinates of the sill and range, as shown in the legend, and the white dot represents the nugget. The <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math> value for each plot is also included in the legend.</p>
Full article ">Figure 5
<p>Probability Density Plot of <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>2.5</mn> </mrow> </msub> </mrow> </semantics></math> measurement time across 3 distinctly polluted days of January 2023 for Joppa, Dallas, Texas. The peak measurement times are marked at 12.25 s, 15.05 s, and 18.92 s for the highest, intermediate, and least polluted days, respectively.</p>
Full article ">Figure 6
<p>Scatter plots of the hyperparameter-optimized Random Forest Regression models for PM measurement time (minutes) on 2 January 2023, in Joppa, Dallas, Texas. Panels (<b>a</b>–<b>g</b>) correspond to the measurement times of <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>0.1</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>0.3</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>0.5</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>1.0</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>2.5</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>5.0</mn> </mrow> </msub> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>10.0</mn> </mrow> </msub> </mrow> </semantics></math>, respectively. The blue dots represent the training data, while the orange dots represent the testing data. The <math display="inline"><semantics> <msup> <mi>R</mi> <mn>2</mn> </msup> </semantics></math> values and the marginal distributions for both training and testing data are provided, along with the 1:1 reference line for evaluating prediction accuracy.</p>
Full article ">Figure 7
<p>Permutation importance rankings of features influencing estimated PM measurement time (in minutes) for various particulate matter (PM) size fractions on 2 January 2023, in Joppa, Texas. Panels (<b>a</b>–<b>g</b>) correspond to the measurement times of <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>0.1</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>0.3</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>0.5</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>1.0</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>2.5</mn> </mrow> </msub> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>5.0</mn> </mrow> </msub> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msub> <mi>PM</mi> <mrow> <mn>10.0</mn> </mrow> </msub> </mrow> </semantics></math>, respectively. Each bar chart displays the relative importance of meteorological parameters—humidity, temperature, pressure, and wind speed—in influencing the measurement time for each PM size fraction.</p>
Full article ">
23 pages, 17900 KiB  
Article
Unveiling the Impact of Microfractures on Longitudinal Dispersion Coefficients in Porous Media
by Muyuan Wang, Keliu Wu, Qingyuan Zhu and Jiawei Ye
Processes 2025, 13(3), 722; https://doi.org/10.3390/pr13030722 (registering DOI) - 2 Mar 2025
Abstract
Longitudinal dispersion coefficient is a key parameter governing solute transport in porous media, with significant implications for various industrial processes. However, the impact of microfractures on the longitudinal dispersion coefficient remains insufficiently understood. In this study, pore-scale direct numerical simulations are performed to [...] Read more.
Longitudinal dispersion coefficient is a key parameter governing solute transport in porous media, with significant implications for various industrial processes. However, the impact of microfractures on the longitudinal dispersion coefficient remains insufficiently understood. In this study, pore-scale direct numerical simulations are performed to analyze solute transport in microfractured porous media during unstable miscible displacement. Spatiotemporal concentration profiles were fitted to the analytical solution of the convection–dispersion equation to quantify the longitudinal dispersion coefficient across different microfracture configurations. The results indicate that the longitudinal dispersion coefficient is highly sensitive to microfracture characteristics. Specifically, an increased projection length of microfractures in the flow direction and a reduced lateral projection length enhance longitudinal dispersion at the outlet. When Peclet number ≥1, the longitudinal dispersion coefficient follows a three-stage variation pattern along the flow direction, with microfracture connectivity and orientation dominating its scale sensitivity. Furthermore, both diffusion-dominated and mixed advective-diffusion regimes are observed. In diffusion-dominated regimes, significant channeling alters the applicability of traditional scaling laws, with the relationship between longitudinal dispersion coefficient and porosity holding only when the Peclet number is below 0.07. These results provide a comprehensive scale-up framework for CO2 miscible flooding in unconventional reservoirs and CO2 storage in saline aquifers, offering valuable insights for the numerical modeling of heterogeneous reservoir development. Full article
(This article belongs to the Section Energy Systems)
Show Figures

Figure 1

Figure 1
<p>Granular porous medium before microfracture incorporation. (Note: The computational domain (<span class="html-italic">L</span> × <span class="html-italic">W</span>) dimensions are 1000 μm × 50 μm; the white regions represent solid particles and the blue regions represent the pore space).</p>
Full article ">Figure 2
<p>SSD and PSD in porous medium.</p>
Full article ">Figure 3
<p>Porous media domains with embedded microfractures. Microfractures in M-PM1 to M-PM4 are parallel to the flow direction but vary in length and density. (<b>a</b>) M-PM1. (<b>b</b>) M-PM2. (<b>c</b>) M-PM3. (<b>d</b>) M-PM4. (<b>e</b>) M-PM5. (<b>f</b>) M-PM6. (<b>g</b>) M-PM7. (<b>h</b>) M-PM8. (Note: In M-PM1 to M-PM4, microfractures are aligned parallel to the flow direction but differ in length and density. M-PM5 and M-PM6 contain microfractures oriented at 90° and 45° to the flow direction, respectively, while M-PM8 exhibits enhanced microfracture connectivity compared to M-PM7).</p>
Full article ">Figure 4
<p>Outlet concentration curves for mesh resolutions during CO<sub>2</sub> miscible displacements in M-PM5 at <span class="html-italic">Pe</span> = 100 (<b>a</b>) and local mesh details in Mesh 3 (<b>b</b>).</p>
Full article ">Figure 5
<p>Comparison of simulated and analytical outlet concentration data in M-PM1 at <span class="html-italic">Pe</span> = 0.01.</p>
Full article ">Figure 6
<p>Comparison of numerical simulation and analytical solutions (Berkowitz and Zhou [<a href="#B48-processes-13-00722" class="html-bibr">48</a>]; Wang et al. [<a href="#B49-processes-13-00722" class="html-bibr">49</a>]) for <span class="html-italic">D<sub>L</sub></span> and <span class="html-italic">Pe</span> in stable miscible flow between two parallel flat plates.</p>
Full article ">Figure 7
<p>Variation in dimensionless <span class="html-italic">D<sub>L</sub></span> (<span class="html-italic">D<sub>L</sub></span>/<span class="html-italic">D</span>) along flow direction from M-PM1 to M-PM8 under different <span class="html-italic">Pe</span>. (<b>a</b>) <span class="html-italic">Pe</span> = 0.01. (<b>b</b>) <span class="html-italic">Pe</span> = 0.1. (<b>c</b>) <span class="html-italic">Pe</span> = 1. (<b>d</b>) <span class="html-italic">Pe</span> = 10. (<b>e</b>) <span class="html-italic">Pe</span> = 50. (<b>f</b>) <span class="html-italic">Pe</span> = 100.</p>
Full article ">Figure 7 Cont.
<p>Variation in dimensionless <span class="html-italic">D<sub>L</sub></span> (<span class="html-italic">D<sub>L</sub></span>/<span class="html-italic">D</span>) along flow direction from M-PM1 to M-PM8 under different <span class="html-italic">Pe</span>. (<b>a</b>) <span class="html-italic">Pe</span> = 0.01. (<b>b</b>) <span class="html-italic">Pe</span> = 0.1. (<b>c</b>) <span class="html-italic">Pe</span> = 1. (<b>d</b>) <span class="html-italic">Pe</span> = 10. (<b>e</b>) <span class="html-italic">Pe</span> = 50. (<b>f</b>) <span class="html-italic">Pe</span> = 100.</p>
Full article ">Figure 8
<p>Impact of strong flow direction heterogeneity on CO<sub>2</sub> concentration front propagation. (<b>a</b>) M-PM1. (<b>b</b>) M-PM2. (<b>c</b>) M-PM8. (Note: <span class="html-italic">c<sub>D</sub></span> distributions for M-PM1, M-PM2, and M-PM3 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 0.1).</p>
Full article ">Figure 9
<p>Impact of single horizontal microfracture length on CO<sub>2</sub> concentration front propagation. (<b>a</b>) M-PM1, <span class="html-italic">Pe</span> = 1. (<b>b</b>) M-PM2, <span class="html-italic">Pe</span> = 1. (<b>c</b>) M-PM1, Pe = 10. (<b>d</b>) M-PM2, <span class="html-italic">Pe</span> = 10. (<b>e</b>) M-PM1, <span class="html-italic">Pe</span> = 100. (<b>f</b>) M-PM2, <span class="html-italic">Pe</span> = 100. (Note: <span class="html-italic">c<sub>D</sub></span> distributions for M-PM1 and M-PM2 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 1, 10, 100).</p>
Full article ">Figure 10
<p>Impact of microfracture parallel density on CO<sub>2</sub> concentration front propagation. (<b>a</b>) M-PM3, <span class="html-italic">Pe</span> = 1. (<b>b</b>) M-PM4, <span class="html-italic">Pe</span> = 1. (<b>c</b>) M-PM3, <span class="html-italic">Pe</span> = 10. (<b>d</b>) M-PM4, <span class="html-italic">Pe</span> = 10. (<b>e</b>) M-PM3, <span class="html-italic">Pe</span> = 100. (<b>f</b>) M-PM4, <span class="html-italic">Pe</span> = 100. (Note: <span class="html-italic">c<sub>D</sub></span> distributions for M-PM3 and M-PM4 at <span class="html-italic">t<sub>D</sub></span> = 0.8 and <span class="html-italic">Pe</span> = 1, 10, 100).</p>
Full article ">Figure 11
<p>Dimensionless CO<sub>2</sub> concentration profiles at <span class="html-italic">x<sub>D</sub></span> = 1 for M-PM1 to M-PM4 as a function of <span class="html-italic">t<sub>D</sub></span> at <span class="html-italic">Pe</span> = 1.</p>
Full article ">Figure 12
<p>Impact of microfracture orientation on CO<sub>2</sub> concentration front propagation at high <span class="html-italic">Pe</span>. (<b>a</b>) M-PM3, <span class="html-italic">Pe</span> = 50. (<b>b</b>) M-PM3, <span class="html-italic">Pe</span> = 100. (<b>c</b>) M-PM5, <span class="html-italic">Pe</span> = 50. (<b>d</b>) M-PM5, <span class="html-italic">Pe</span> = 100. (<b>e</b>) M-PM6, <span class="html-italic">Pe</span> = 50. (<b>f</b>) M-PM6, <span class="html-italic">Pe</span> = 100. (Note: <span class="html-italic">c<sub>D</sub></span> distributions for M-PM3, M-PM5, and M-PM6 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 50, 100).</p>
Full article ">Figure 13
<p>Impact of microfracture connectivity on CO<sub>2</sub> concentration front propagation at high <span class="html-italic">Pe</span>. (<b>a</b>) M-PM7, <span class="html-italic">Pe</span> = 100. (<b>b</b>) M-PM7, <span class="html-italic">Pe</span> = 50. (<b>c</b>) M-PM8, <span class="html-italic">Pe</span> = 100. (Note: <span class="html-italic">c<sub>D</sub></span> for M-PM7 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 50, 100; <span class="html-italic">c<sub>D</sub></span> distribution for M-PM8 at <span class="html-italic">t<sub>D</sub></span> = 0.5 and <span class="html-italic">Pe</span> = 100).</p>
Full article ">Figure 14
<p>Perturbation of CO<sub>2</sub> concentration front in Stage 1 (<b>a</b>–<b>d</b>), Stage 2 (<b>e</b>), and Stage 3 (<b>f</b>).</p>
Full article ">Figure 14 Cont.
<p>Perturbation of CO<sub>2</sub> concentration front in Stage 1 (<b>a</b>–<b>d</b>), Stage 2 (<b>e</b>), and Stage 3 (<b>f</b>).</p>
Full article ">Figure 15
<p>Impact of microfracture orientation and connectivity on CO<sub>2</sub> concentration distribution along domain length. (<b>a</b>) M-PM1. (<b>b</b>) M-PM6. (<b>c</b>) M-PM7. (<b>d</b>) M-PM8.</p>
Full article ">Figure 16
<p>Identification of dispersion regimes in unstable displacement process.</p>
Full article ">
20 pages, 6337 KiB  
Article
Vehicle–Bridge Coupling of Road–Rail Dual-Use Network Arch Bridge Based on a Noniterative Approach: Parametric Analysis and Case Study
by Haocheng Chang, Rujin Ma, Baixue Ge and Qiuying Zhu
Buildings 2025, 15(5), 801; https://doi.org/10.3390/buildings15050801 (registering DOI) - 2 Mar 2025
Viewed by 41
Abstract
In the realm of road–rail dual-use bridges, conducting accurate vehicle–bridge coupling analysis is crucial, as the combined effects of road traffic and rail transit induce complex dynamic challenges. This study investigates a road–rail dual-use network arch bridge, highlighting the dynamic effects induced by [...] Read more.
In the realm of road–rail dual-use bridges, conducting accurate vehicle–bridge coupling analysis is crucial, as the combined effects of road traffic and rail transit induce complex dynamic challenges. This study investigates a road–rail dual-use network arch bridge, highlighting the dynamic effects induced by light rail loadings. By employing a noniterative vehicle–bridge coupling analysis method, the dynamic responses of hangers caused by vehicular and light rail loads are effectively captured. Additionally, this study explores the influence of various parameters, including vehicle types, driving lanes, and road surface roughness on the responses of hangers positioned at different locations along the bridge. The findings reveal that light rail induces significantly larger dynamic effects compared to motor vehicles. When the light rail operates closer to the hanger, the responses of hangers are more pronounced. Furthermore, different road surface roughness level notably affects the amplitude of axial stress and bending moment fluctuations. Poorer road conditions amplify these dynamic effects, leading to increased stress variations. These insights underscore the necessity of integrating considerations for both road and rail traffic in the structural analysis and design of network arch bridges to ensure their reliability and serviceability. Full article
Show Figures

Figure 1

Figure 1
<p>Illustration of the light-rail-vehicle model.</p>
Full article ">Figure 2
<p>(<b>a</b>) Elevation layout of the Qilu Yellow River Bridge; (<b>b</b>) lane layout of the main span. Note: Reprinted with permission from ref. [<a href="#B27-buildings-15-00801" class="html-bibr">27</a>], 2024, Elsevier.</p>
Full article ">Figure 3
<p>Finite element model and boundary conditions of the Qilu Yellow River Bridge.</p>
Full article ">Figure 4
<p>The axial forces in the hanger rods under dead load.</p>
Full article ">Figure 5
<p>Flowchart of the noniterative method for vehicle–bridge coupling analysis.</p>
Full article ">Figure 6
<p>Time-dependent axial stress of hangers considering different vehicle types: (<b>a</b>): hanger A; (<b>b</b>) hanger B; (<b>c</b>): hanger C.</p>
Full article ">Figure 6 Cont.
<p>Time-dependent axial stress of hangers considering different vehicle types: (<b>a</b>): hanger A; (<b>b</b>) hanger B; (<b>c</b>): hanger C.</p>
Full article ">Figure 7
<p>(<b>a</b>) Vertical displacement of the upper and lower vertices of hanger C; (<b>b</b>) vertical displacement difference between upper and lower vertices of hanger C.</p>
Full article ">Figure 8
<p>Time-dependent axial stress of the hangers under the light rail loading: (<b>a</b>): hanger A; (<b>b</b>) hanger B; (<b>c</b>) hanger C.</p>
Full article ">Figure 8 Cont.
<p>Time-dependent axial stress of the hangers under the light rail loading: (<b>a</b>): hanger A; (<b>b</b>) hanger B; (<b>c</b>) hanger C.</p>
Full article ">Figure 9
<p>Time-dependent bending moment of hanger C, considering different vehicle types: (<b>a</b>) longitudinal bending moment; (<b>b</b>) transverse bending moment.</p>
Full article ">Figure 10
<p>Time-dependent bending moment of hanger C under the light rail loading: (<b>a</b>) longitudinal bending moment; (<b>b</b>) transverse bending moment.</p>
Full article ">Figure 11
<p>Time-dependent responses of hanger C, considering different lanes: (<b>a</b>) axial stress; (<b>b</b>) longitudinal bending moment; (<b>c</b>) transverse bending moment.</p>
Full article ">Figure 12
<p>Time−dependent responses of hanger C, considering light rail traveling in different lanes: (<b>a</b>) axial stress; (<b>b</b>) longitudinal bending moment; (<b>c</b>) transverse bending moment.</p>
Full article ">Figure 12 Cont.
<p>Time−dependent responses of hanger C, considering light rail traveling in different lanes: (<b>a</b>) axial stress; (<b>b</b>) longitudinal bending moment; (<b>c</b>) transverse bending moment.</p>
Full article ">Figure 13
<p>Time-dependent responses of hanger C, considering different road surface roughness levels: (<b>a</b>) axial stress; (<b>b</b>) longitudinal bending moment; (<b>c</b>) transverse bending moment.</p>
Full article ">
15 pages, 4314 KiB  
Article
Research on Predictive Control Strategy for Phase Shift Full Bridge Transform Based on Improved Nonlinear Disturbance Observer
by Yunbo Wu, Cheng Liu, Qing Zhao and Lixin Liu
Electronics 2025, 14(5), 1002; https://doi.org/10.3390/electronics14051002 - 2 Mar 2025
Viewed by 144
Abstract
To enhance the power output performance of high-power electrolytic plating equipment and improve the dynamic response capability and disturbance rejection ability of the DC-DC converter—a core component in electrolytic plating systems—this study proposes a predictive control strategy for phase-shifted full-bridge converters based on [...] Read more.
To enhance the power output performance of high-power electrolytic plating equipment and improve the dynamic response capability and disturbance rejection ability of the DC-DC converter—a core component in electrolytic plating systems—this study proposes a predictive control strategy for phase-shifted full-bridge converters based on an improved nonlinear disturbance observer. The implementation framework comprises three key technical components: Firstly, a model predictive control (MPC)-based inner-loop controller architecture is constructed to optimize the dynamic response characteristics of the current inner-loop system. Subsequently, an enhanced nonlinear disturbance observer is designed to accurately estimate parameter variations in converter electronic components and external disturbances. Finally, a feedforward compensation module is developed to mitigate the inherent duty cycle loss phenomenon in phase-shifted full-bridge converters. Simulation results demonstrate that the proposed control strategy significantly improves system dynamic performance, achieving a 25% reduction in settling time compared with conventional methods while maintaining robust disturbance rejection capabilities under ±10% voltage fluctuations. This integrated approach effectively addresses the conflicting requirements between dynamic response speed and anti-interference performance in high-power electrochemical process control systems. Full article
(This article belongs to the Special Issue Control and Optimization of Power Converters and Drives)
Show Figures

Figure 1

Figure 1
<p>Topological structure diagram of phase-shifted full bridge DC-DC converter.</p>
Full article ">Figure 2
<p>Schematic diagram of drive logic.</p>
Full article ">Figure 3
<p>Disturbance observer predictive current control block diagram.</p>
Full article ">Figure 4
<p>Model predictive control flowchart.</p>
Full article ">Figure 5
<p>12.5 KW phase-shifted full bridge converter simulation.</p>
Full article ">Figure 6
<p>Voltage response curve of different PI parameters.</p>
Full article ">Figure 7
<p>Model predictive control voltage response curve.</p>
Full article ">Figure 8
<p>Current response curve under PI control and model predictive control.</p>
Full article ">Figure 9
<p>Current response diagram under sudden load changes.</p>
Full article ">Figure 10
<p>Current response diagram of output inductance parameter mismatch under working conditions.</p>
Full article ">Figure 11
<p>Current response diagram of random current disturbance under working conditions.</p>
Full article ">
26 pages, 13225 KiB  
Article
A New and Tested Ionospheric TEC Prediction Method Based on SegED-ConvLSTM
by Yuanhang Liu, Yingkui Gong, Hao Zhang, Ziyue Hu, Guang Yang and Hong Yuan
Remote Sens. 2025, 17(5), 885; https://doi.org/10.3390/rs17050885 (registering DOI) - 2 Mar 2025
Viewed by 196
Abstract
Total electron content (TEC) serves as a key parameter characterizing ionospheric conditions. Accurate prediction of TEC plays a crucial role in improving the precision of Global Navigation Satellite Systems (GNSS). However, existing research have predominantly emphasized spatial variations in the ionosphere, neglecting the [...] Read more.
Total electron content (TEC) serves as a key parameter characterizing ionospheric conditions. Accurate prediction of TEC plays a crucial role in improving the precision of Global Navigation Satellite Systems (GNSS). However, existing research have predominantly emphasized spatial variations in the ionosphere, neglecting the periodic changes of the ionosphere with the diurnal cycle. In this paper, we propose a TEC prediction model, which simultaneously considers both spatial and temporal characteristics to extract spatiotemporal features of ionospheric distribution. Additionally, we integrate several space weather element datasets into the prediction model framework, allowing the generation of multiple space weather feature values that represent the influence of space weather on the ionosphere at different latitudes and longitudes. Moreover, we apply Gaussian process regression (GPR) interpolation to geomagnetic data to characterize impact on the ionosphere, thereby enhancing the prediction accuracy. We compared our model with traditional image-based models such as convolutional neural networks (CNNs), convolutional long short-term memory networks (ConvLSTMs), a self-attention mechanism-integrated ConvLSTM (SAM-ConvLSTM) model, and one-day predicted ionospheric products (C1PG) provided by the Center for Orbit Determination in Europe (CODE). We also examined the effect of using different numbers of space weather feature values in these models. Our model outperforms the comparison models in terms of prediction error metrics, including mean absolute error (MAE), root mean square error (RMSE), correlation coefficient (CC), and the structural similarity index (SSIM). Furthermore, we analyzed the influence of different batch sizes on model training accuracy to find the best performance of each model. In addition, we investigated the model performance during geomagnetic quiet periods, where our model provided the most accurate predictions and demonstrates higher prediction accuracy in the equatorial anomaly region. We also analyzed the prediction performance of all models during space weather events. The results indicate that the proposed model is the least affected during geomagnetic storms and demonstrates superior prediction performance compared to other models. This study presents a more stable and high-performance spatiotemporal prediction model for TEC. Full article
Show Figures

Figure 1

Figure 1
<p>The time series of different space weather event elements and total TEC unit per hour during 2012–2024.</p>
Full article ">Figure 2
<p>The sliding window method used in this study divides continuous time segments into 60-day periods, with each period allocated as follows: 40 days for the training set and 20 days for the test set.</p>
Full article ">Figure 3
<p>SegED-ConvLSTM architecture.</p>
Full article ">Figure 4
<p>Global TEC distribution at the same time points for two consecutive days.</p>
Full article ">Figure 5
<p>ConvLSTM Structure.</p>
Full article ">Figure 6
<p>ED-ConvLSTM structure.</p>
Full article ">Figure 7
<p>Model Prediction Under Space Weather Anomaly Events.</p>
Full article ">Figure 8
<p>2014-05-12 Comparison between GroundTruth, TS-ConvLSTM-6, SegED-ConvLSTM-1, SegED-ConvLSTM-6.</p>
Full article ">Figure 9
<p>Absolute error distribution of TEC predictions for 2014-05-12 across different models.</p>
Full article ">Figure 10
<p>Comparison of ionospheric TEC predictions: ground truth vs. COPG and SegED-ConvLSTM-6 models.</p>
Full article ">Figure 11
<p>Comparison of SegED-ConvLSTM-6 and COPG at 12 May 2014.</p>
Full article ">
19 pages, 7258 KiB  
Article
Exploring Post-Machining Alternatives Under Dry Conditions for Thin-Walled Additive Manufacturing Components Aided by Infrared Thermography
by Eduard Garcia-Llamas, Giselle Ramirez, Miguel Fuentes, Eduard Vidales and Jaume Pujante
Processes 2025, 13(3), 717; https://doi.org/10.3390/pr13030717 (registering DOI) - 1 Mar 2025
Viewed by 277
Abstract
Additive manufacturing (AM) techniques have transformed the production of parts and components with intricate geometries and customized designs, driving innovation in sustainable manufacturing practices. The additive manufacturing technology used in this work was selective laser melting (SLM), a process that uses laser energy [...] Read more.
Additive manufacturing (AM) techniques have transformed the production of parts and components with intricate geometries and customized designs, driving innovation in sustainable manufacturing practices. The additive manufacturing technology used in this work was selective laser melting (SLM), a process that uses laser energy to sinter powdered metals into solid structures. Among the various materials utilized in AM, Ti6Al4V titanium alloys are of particular interest due to their favorable mechanical properties, corrosion resistance, biocompatibility, and potential for reducing material waste. However, the machining of additively manufactured titanium parts presents challenges due to the material’s low conductivity, elastic modulus, and chemical affinity with cutting tools, which impact tool wear and surface finish quality. Milling, a commonly employed process for finishing titanium parts, often involves significant energy use and tool wear, highlighting the need for optimized and eco-conscious machining strategies. This study aims to establish correlations among four key aspects: (1) surface finish of machined Ti6Al4V AM parts, (2) cutting tool damage, (3) dry milling parameters including different cutting tools, and (4) variation of temperature at the contact surface of AM parts and tools using infrared thermography. By examining parameters such as feed per tooth (Fz), axial depth of cut (Ap), spindle trajectories (trochoidal, helicoidal, and linear), and cutting tool diameters, this work identifies conditions that enhance process efficiency while reducing environmental impact. Infrared thermography provides insights into temperature variations during milling, correlating these changes to surface roughness and critical machining parameters, thus contributing to the development of sustainable and high-performance manufacturing practices. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Schematic of the SLM component machined by milling process in this study, along its dimensions.</p>
Full article ">Figure 2
<p>(<b>a</b>) Image of cutting tool XDLX 09T308ER-F40 CTC5240 and the tool holder GA SD090 C 032 Z3; (<b>b</b>) Image of cutting tool T290 LNMT 100405TR IC808 and the tool holder T290 ELN D20-03-W20-10.</p>
Full article ">Figure 3
<p>Scheme of the three different trajectories used in the tests.</p>
Full article ">Figure 4
<p>From left to right, the figure shows the machined component with the analyzed surface highlighted by a black rectangle, followed by the extracted 3D image and the graph illustrating the surface roughness variation.</p>
Full article ">Figure 5
<p>(<b>a</b>) Experimental setup with the five-axis CNC machining center, the HAAS UMC-750. Inside, the SLM component is visible, and the black square indicates the location of the infrared camera; (<b>b</b>) FLIR SC645 high-resolution LWIR infrared camera; (<b>c</b>) Example of an infrared image acquired by the infrared camera. The two registered temperature areas are represented by two rectangles. The bottom rectangle focuses on the component temperature, while the top rectangle focuses on the tool temperature. The red triangles indicate the maximum temperature in each area.</p>
Full article ">Figure 6
<p>Variation of Ra and RSM roughness parameters while keeping Ap fixed at 0.1 and increasing Fz with values of 0.1 and 0.5 (<b>left chart</b>), and while keeping Fz fixed at 0.1 and increasing Ap with values of 0.1 and 0.25 (<b>right chart</b>).</p>
Full article ">Figure 7
<p>(<b>a</b>) The chart shows the maximum temperatures recorded during the milling process for different trajectories. (<b>b</b>) Surface roughness parameters measured on the component for different milling trajectories. (<b>c</b>) Image of the heat-affected zone (HAZ) on the cutting tool edge (highlighted in a yellow square) for the various trajectories, along with a chart (<b>d</b>) illustrating the extent of the HAZ.</p>
Full article ">Figure 8
<p>(<b>a</b>) Correlation between three levels of cutting tool wear and the corresponding maximum temperature changes measured by the infrared camera. (<b>b</b>) Evolution of Ra and RSM roughness of the AM component obtained with the three levels of tool wear.</p>
Full article ">Figure 8 Cont.
<p>(<b>a</b>) Correlation between three levels of cutting tool wear and the corresponding maximum temperature changes measured by the infrared camera. (<b>b</b>) Evolution of Ra and RSM roughness of the AM component obtained with the three levels of tool wear.</p>
Full article ">Figure 9
<p>Temperature evolution recorded by the infrared camera for Setup 1 and Setup 2 during the milling process with a trochoidal trajectory.</p>
Full article ">Figure 10
<p>Evolution of Ra and RSM roughness parameters with the three levels of tool wear for Setup 1 and Setup 2.</p>
Full article ">Figure 11
<p>At (<b>left</b>), 3D images of the worn cutting tools for Setup 1 and for Setup 2 are shown. At the (<b>right</b>), charts depicting tool damage, including adhesion and wear, are presented.</p>
Full article ">Figure 12
<p>On the (<b>left</b>), images of the cutting tools for the compound fixed with a vise and the compound with 3D-printed thin-wall support, both for Setup 2, are shown. On the (<b>right</b>), a chart depicting tool damage, including adhesion and wear, is presented.</p>
Full article ">Figure 13
<p>Three-dimensional surface images of the milled, 3D-printed Ti-6Al-4V surface fixed with the (<b>a</b>) vise and (<b>b</b>) thin-walled supported for Setup 2 and the Ra and RSM roughness values.</p>
Full article ">
15 pages, 3561 KiB  
Article
Classification and Recognition of Soybean Quality Based on Hyperspectral Imaging and Random Forest Methods
by Man Chen, Zhichang Chang, Chengqian Jin, Gong Cheng, Shiguo Wang and Youliang Ni
Sensors 2025, 25(5), 1539; https://doi.org/10.3390/s25051539 - 1 Mar 2025
Viewed by 313
Abstract
To achieve the rapid and accurate classification and identification of soybean components, this study selected soybeans harvested by the 4LZ-1.5 soybean combine harvester as the research subject. Hyperspectral images of soybean samples were collected using the Pika L spectrometer, and spectral information was [...] Read more.
To achieve the rapid and accurate classification and identification of soybean components, this study selected soybeans harvested by the 4LZ-1.5 soybean combine harvester as the research subject. Hyperspectral images of soybean samples were collected using the Pika L spectrometer, and spectral information was extracted from the regions of interest (ROI) in the images. Eight preprocessing methods, including baseline correction (BC), moving average (MA), Savitzky–Golay derivative (SGD), normalization, standard normal variate transformation (SNV), multiplicative scatter correction (MSC), first derivative (DS), and Savitzky–Golay smoothing (SGS), were applied to the raw spectral data to eliminate irrelevant information. Feature wavelengths were selected using the successive projections algorithm (SPA) and the competitive adaptive reweighted sampling (CARS) algorithm to reduce spectral redundancy and enhance model detection performance, retaining eight and ten feature wavelengths, respectively. Subsequently, a random forest (RF) model was developed for soybean component classification. The model parameters were optimized using particle swarm optimization (PSO) and differential evolution (DE) algorithms to improve performance. Experimental results showed that the RF classification model based on SPA-BC preprocessed spectra and DE-tuned parameters achieved an optimal prediction accuracy of 1.0000 during training. This study demonstrates the feasibility of using hyperspectral imaging technology for the rapid and accurate detection of soybean components, providing technical support for the assessment of breakage and impurity levels during soybean harvesting and storage processes. It also offers a reference for the development of future machine-harvested soybean breakage and impurity detection systems. Full article
(This article belongs to the Section Smart Agriculture)
Show Figures

Figure 1

Figure 1
<p>Mechanically harvested soybean components.</p>
Full article ">Figure 2
<p>Spectral imaging system.</p>
Full article ">Figure 3
<p>Extraction of spectral information of soybean components.</p>
Full article ">Figure 4
<p>Spectral information of different components.</p>
Full article ">Figure 5
<p>The importance of different preprocessed spectral modeling bands. (<b>a</b>) Raw data; (<b>b</b>) BC; (<b>c</b>) normalization.</p>
Full article ">Figure 6
<p>Recognition accuracy of models with different characteristic wavelength numbers (SPA).</p>
Full article ">Figure 7
<p>The optimal feature wavelength obtained through SPA.</p>
Full article ">Figure 8
<p>Recognition accuracy of models with different characteristic wavelength numbers (CARS).</p>
Full article ">Figure 9
<p>The optimal feature wavelength obtained through CARS.</p>
Full article ">Figure 10
<p>The contributions of characteristic wavelength modeling. (<b>a</b>) DE-SPA-BC; (<b>b</b>) DE-CARS.</p>
Full article ">
19 pages, 9441 KiB  
Article
Modeling and Design Aspects of Shallow Geothermal Energy Piles—A Case Study on Large Commercial Building Complex in Zagreb, Croatia
by Marija Macenić and Tomislav Kurevija
Geosciences 2025, 15(3), 90; https://doi.org/10.3390/geosciences15030090 (registering DOI) - 1 Mar 2025
Viewed by 188
Abstract
With ambitious targets set by the EU for the reduction of emissions from the energy sector by 2030, there is a need to design and develop more building projects using renewable energy sources. Even though in Europe, heating and cooling share from renewable [...] Read more.
With ambitious targets set by the EU for the reduction of emissions from the energy sector by 2030, there is a need to design and develop more building projects using renewable energy sources. Even though in Europe, heating and cooling share from renewable resources is increasing, and in 2021, the total share in this sector in Croatia was at 38%, the share of heat production by heat pumps is rather low. One possibility to increase this share is to install energy piles when constructing a building, which is becoming an increasingly common practice. This case study focuses on such a system designed for a large, non-residential building in Zagreb, Croatia. The complex was designed as 13 separate dilatations, with central heating and cooling of all facilities, covered by 260 energy piles (130 pairs in serial connection), with a length of the polyethylene pipe of 20 m in a double loop inserted within the pile. The thermo-technical system was designed as a bivalent parallel system, with natural gas covering peak heating loads and a dry cooler covering cooling peak loads when the loads cannot be covered only by ground-source heat pumps. In the parallel bivalent system, the geothermal source will work with a much higher number of working hours at full load than is the case for geothermal systems that are dimensioned to peak consumption. Therefore, the thermal response test was conducted on two energy piles, connected in series, to obtain thermogeological parameters and determine the heat extraction and rejection rates. The established steady-state heat rate defines the long-term ability to extract heat energy during constant thermal load, with the inlet water temperature from the pile completely stabilized, i.e., no significant further sub-cooling is achieved in the function of the geothermal field operation time. Considering the heating and cooling loads of the building, modeling of the system was performed in such a manner that it utilized renewable energy as much as possible by finding a bivalent point where the geothermal system works efficiently. It was concluded that the optimal use of the geothermal field covers total heating needs and 70% for cooling, with dry coolers covering the remaining 30%. Additionally, based on the measured thermogeological parameters, simulations of the thermal response test were conducted to determine heat extraction and rejection rates for energy piles with various geometrical parameters of the heat exchanger pipe and fluid flow variations. Full article
Show Figures

Figure 1

Figure 1
<p>Principle of steady-state thermal response step test evolution as a function of heat power and time.</p>
Full article ">Figure 2
<p>Legend for geology map: <b>a</b>–alluvium: gravels, sands, and clays; <b>a<sub>2</sub></b>–middle terrace: gravels and sands; <b>ap</b>–flood facies: clayey sand, clayey silt; <b>b</b>–marsh sediments: clays, clayey silt; <b>dpr</b>–diluvium–proluvium: fine-grained clays, sands, and gravels; <b>l</b>–clayey silt; <b>lb</b>–marshy loes: silty clays. General geological map of the wider Zagreb area and location of the non-residential case study [<a href="#B17-geosciences-15-00090" class="html-bibr">17</a>].</p>
Full article ">Figure 3
<p>Coring samples were taken from geotechnical research of the location prior to the drilling of the energy piles.</p>
Full article ">Figure 4
<p>(<b>a</b>) Pilot setup with heat exchangers, (<b>b</b>) energy piles of 20 m in length ready to be installed with U-loop pipe, (<b>c</b>) installation of TRT equipment on two serial coupled piles at the construction field.</p>
Full article ">Figure 5
<p>First heat step and determination of effective thermal conductivity of the heterogeneous ground.</p>
Full article ">Figure 6
<p>Entire thermal response test of the ground with four heat steps and one extended condition (dotted lines).</p>
Full article ">Figure 7
<p>Steady-state heat rejection/extraction rates during long-term operation in the function of the fluid temperature for the installed energy piles (with corrected static temperature during the winter period).</p>
Full article ">Figure 8
<p>Influence of thermal interference between energy pile system for different radii and the annual full-load hours.</p>
Full article ">Figure 9
<p>(<b>a</b>) Pile grid arrangement of 256 × 20 m energy piles in different concrete dilatations; (<b>b</b>) current view of commercial buildings after system completion.</p>
Full article ">Figure 10
<p>Temperature evolution of brine in energy piles over 20 years for Scenario 1.</p>
Full article ">Figure 11
<p>Temperature evolution of brine in energy piles over 20 years for Scenario 2.</p>
Full article ">
21 pages, 1271 KiB  
Article
Human and Machine Reliability Estimation in Discrete Simulations and Machine Learning for Industry 4.0 and 5.0
by Wojciech M. Kempa, Iwona Paprocka, Bożena Skołud and Grzegorz Ćwikła
Symmetry 2025, 17(3), 377; https://doi.org/10.3390/sym17030377 (registering DOI) - 1 Mar 2025
Viewed by 178
Abstract
Currently, Industry 4.0 creates new opportunities for analyzing data on production processes and extracting knowledge from them. With the Internet of Things, data is continuously collected from machine sensors to analyze machine health. Thanks to artificial intelligence methods and discrete simulation, it is [...] Read more.
Currently, Industry 4.0 creates new opportunities for analyzing data on production processes and extracting knowledge from them. With the Internet of Things, data is continuously collected from machine sensors to analyze machine health. Thanks to artificial intelligence methods and discrete simulation, it is possible to process data and dynamically adjust the operating conditions of the production line to the expected time of failure-free operation of the machine or reliable work of an employee. Recently, machine learning techniques have been used to automatically adapt the production line to changes in a given production environment. The paper presents various methods of modeling actions, i.e., forecasting the failure-free operation time of a machine or the error-free working time of an employee. The possible actions the agent can perform, the possible prediction techniques that can be selected are presented. The time between failures is described by a log-normal distribution. The asymmetric lognormal distribution is much more flexible for practical modeling compared to the “perfectly” symmetric normal distribution. In practice, the asymmetric lognormal distribution, strongly shifted to the left, can be used to describe the decreasing time between failures due to human error, as well as the time between failures of a machine in the third phase of its life cycle, which decreases as the machine ages and its components wear out. The parameters of the distribution are estimated using the maximum-likelihood approach, theempirical moments approach, the renewal-theory approach, the empirical distribution function and the method based on coefficient of variation. Numerical examples of predicting failure-free operation times described by the log-normal distribution are presented. The results are compared assuming that failure-free times are described by exponential, normal and Weibull distributions. The results are also compared with an example of the simplest learning method. Full article
(This article belongs to the Section Engineering and Materials)
Show Figures

Figure 1

Figure 1
<p>Data on (<b>a</b>) failure-free, repair, and operation times of pumps in November and (<b>b</b>) failure-free times for 5 months.</p>
Full article ">Figure 2
<p>Failure -free, repair, and operation times of the first pump for November–March.</p>
Full article ">Figure 3
<p>Parameter estimation of log-normal distributions: (<b>a</b>) maximum likelihood approach and (<b>b</b>) empirical moments approach.</p>
Full article ">Figure 4
<p>Parameter estimation of distributions: (<b>a</b>) Weibull, (<b>b</b>) exponential, and (<b>c</b>) truncated normal.</p>
Full article ">Figure 5
<p>SCADA file from 2 February containing the operating conditions of the first pump.</p>
Full article ">Figure 6
<p>Ex post errors of failure-free time predictions using predictive maintenance, condition-based maintenance, and the Monte Carlo-based approach.</p>
Full article ">Figure 7
<p>Parameter estimation of the log-normal distribution using the Monte Carlo-based approach.</p>
Full article ">Figure 8
<p>The agent–environment interaction in a Markov decision process for predictive maintenance of wastewater treatment.</p>
Full article ">
21 pages, 12426 KiB  
Article
Scientific Molding and Adaptive Process Quality Control with External Sensors for Injection Molding Process
by Chen-Hsiang Chang, Chien-Hung Wen, Ren-Ho Tseng, Chieh-Hsun Tsai, Yu-Hao Chen, Sheng-Jye Hwang and Hsin-Shu Peng
Technologies 2025, 13(3), 97; https://doi.org/10.3390/technologies13030097 (registering DOI) - 1 Mar 2025
Viewed by 228
Abstract
This study established a real-time measurement system to monitor the melt quality in an injection molding process using a pressure sensor installed on the nozzle and a strain gauge installed on the tie bar. Based on the sensing curves from these two external [...] Read more.
This study established a real-time measurement system to monitor the melt quality in an injection molding process using a pressure sensor installed on the nozzle and a strain gauge installed on the tie bar. Based on the sensing curves from these two external sensors, the characteristic values of nozzle pressure and clamping force were used to optimize parameters. This study defined product weight as a quality indicator and developed a scientific molding parameter setup process. The optimization sequence of parameters is injection speed, V/P switchover point, packing pressure, packing time, and clamping force. Finally, an adaptive process control system was established based on the online quality characteristic values to maintain product quality consistency. Continuous production experiments were conducted at two sites to verify the system’s effectiveness. The results revealed that the optimized process parameters can ensure product weight stability during long-term production. Furthermore, using the adaptive process control system further enhanced product weight stability at both sites, reducing the standard deviation of product weight to 0.0289 g and 0.0148 g, and the coefficient of variation to 0.065% and 0.035%, respectively. Full article
(This article belongs to the Section Manufacturing Technology)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Changes in melt viscosity at various shear rates (Courtesy: Suhas Kulkarni, FIMMTECH INC).</p>
Full article ">Figure 2
<p>Nozzle pressure characteristic values.</p>
Full article ">Figure 3
<p>Clamping force difference value.</p>
Full article ">Figure 4
<p>Adaptive process control system flowchart.</p>
Full article ">Figure 5
<p>(<b>a</b>) The dimensions of the phone stand sample and (<b>b</b>) mold of the part.</p>
Full article ">Figure 6
<p>(<b>a</b>) The nozzle pressure sensor and (<b>b</b>) the strain sensor.</p>
Full article ">Figure 7
<p>Experimental measurement system.</p>
Full article ">Figure 8
<p>Nozzle pressure profile with various injection speeds.</p>
Full article ">Figure 9
<p>The relationship between nozzle peak pressure and timing of peak pressure.</p>
Full article ">Figure 10
<p>Relative viscosity with various injection speeds.</p>
Full article ">Figure 11
<p>Nozzle pressure profile with various <span class="html-italic">V</span>/<span class="html-italic">P</span> switchover points.</p>
Full article ">Figure 12
<p>The relationship between nozzle peak pressure and product weight.</p>
Full article ">Figure 13
<p>Nozzle pressure profiles with various packing pressures.</p>
Full article ">Figure 14
<p>Nozzle pressure profiles and screw positions for packing pressures of (<b>a</b>) 25 bar, (<b>b</b>) 125 bar, and (<b>c</b>) 225 bar.</p>
Full article ">Figure 15
<p>Product weights and weight difference with various packing pressures.</p>
Full article ">Figure 16
<p>Nozzle pressure profiles with various packing times.</p>
Full article ">Figure 17
<p>Product weights and weight difference with various packing times.</p>
Full article ">Figure 18
<p>Clamping force profiles with various clamping force settings.</p>
Full article ">Figure 19
<p>Clamping force difference value and product weights with various clamping force settings.</p>
Full article ">Figure 20
<p>Product weight control results of the adaptive process control system.</p>
Full article ">Figure 21
<p>Results of parameter adjustment.</p>
Full article ">Figure 22
<p>The results of nozzle peak pressure and viscosity index (<b>a</b>) without system, and (<b>b</b>) with system.</p>
Full article ">Figure 23
<p>The results of clamping force difference value (<b>a</b>) without system, and (<b>b</b>) with system.</p>
Full article ">Figure 24
<p>Product weight control results of the adaptive process control system at the Beta site.</p>
Full article ">Figure 25
<p>Results of parameter adjustment at the Beta site.</p>
Full article ">Figure 26
<p>Results of nozzle peak pressure and viscosity index at the Beta site (<b>a</b>) without system, and (<b>b</b>) with system.</p>
Full article ">Figure 27
<p>Results of clamping force difference value at the Beta site (<b>a</b>) without system, and (<b>b</b>) with system.</p>
Full article ">
24 pages, 21314 KiB  
Article
RELAY: Robotic EyeLink AnalYsis of the EyeLink 1000 Using an Artificial Eye
by Anna-Maria Felßberg and Dominykas Strazdas
Vision 2025, 9(1), 18; https://doi.org/10.3390/vision9010018 - 1 Mar 2025
Viewed by 171
Abstract
The impact of ambient brightness surroundings on the peak velocities of visually guided saccades remains a topic of debate in the field of eye-tracking research. While some studies suggest that saccades in darkness are slower than in light, others question this finding, citing [...] Read more.
The impact of ambient brightness surroundings on the peak velocities of visually guided saccades remains a topic of debate in the field of eye-tracking research. While some studies suggest that saccades in darkness are slower than in light, others question this finding, citing inconsistencies influenced by factors such as pupil deformation during saccades, gaze position, or the measurement technique itself. To investigate these, we developed RELAY (Robotic EyeLink AnalYsis), a low-cost, stepper motor-driven artificial eye capable of simulating human saccades with controlled pupil, gaze directions, and brightness. Using the EyeLink 1000, a widely employed eye tracker, we assessed accuracy and precision across three illumination settings. Our results confirm the reliability of the EyeLink 1000, demonstrating no artifacts in pupil-based eye tracking related to brightness variations. This suggests that previously observed changes in peak velocities with varying brightness are likely due to human factors, warranting further investigation. However, we observed systematic deviations in measured pupil size depending on gaze direction. These findings emphasize the importance of reporting illumination conditions and gaze parameters in eye-tracking experiments to ensure data consistency and comparability. Our novel artificial eye provides a robust and reproducible platform for evaluating eye tracking systems and deepening our understanding of the human visual system. Full article
(This article belongs to the Section Visual Neuroscience)
Show Figures

Figure 1

Figure 1
<p>Mean peak velocities across participants for a 10° amplitude from our previous study [<a href="#B2-vision-09-00018" class="html-bibr">2</a>].</p>
Full article ">Figure 2
<p>Lateral distance changes between corneal reflection (CR) and pupil based on gaze direction.</p>
Full article ">Figure 3
<p>Pupil foreshortening effect: the observable shape of the pupil changes with gaze directions.</p>
Full article ">Figure 4
<p>Different visible pupil sizes depending on the thickness of the iris (A: with thick iris, B: with thin iris).</p>
Full article ">Figure 5
<p>Double pivot hinge that allows unrestricted movement in two axes, <math display="inline"><semantics> <mi>α</mi> </semantics></math> and <math display="inline"><semantics> <mi>β</mi> </semantics></math>.</p>
Full article ">Figure 6
<p>The final RELAY assembly: two stepper motors connected to an analog stick with an artificial prop eye. An Arduino UNO with a CNC expansion shield and two stepper drivers powers the motors (12v). The lower motor is attached to the base using a dual tilt adjustment adapter.</p>
Full article ">Figure 7
<p>Schematic overview of the test setup.</p>
Full article ">Figure 8
<p>Setup in the soundproof cabin. The operator sat on the other side of the glass window, which was occluded for the experiment.</p>
Full article ">Figure 9
<p>Laser beam reflected by the mirror on the eye.</p>
Full article ">Figure 10
<p>Laser experiment setup.</p>
Full article ">Figure 11
<p>Visualization of the trial design for all three brightness conditions.</p>
Full article ">Figure 12
<p>Planned patterns vs. long exposure laser patterns.</p>
Full article ">Figure 13
<p>Validation of the 13-point calibration pattern as used by the EyeLink 1000 Eye Tracker.</p>
Full article ">Figure 14
<p>Eye tracker data for 13-point calibration pattern, including theoretical, laser, raw, and calculated saccade values.</p>
Full article ">Figure 15
<p>Temporal Graph View (EyeLink Data Viewer 4.4.1) of a random trial from the Artificial Saccades Experiment. The X-axis represents time, while the Y-axis displays multiple gaze-related parameters: X (red) and Y (green) gaze locations in screen pixels, pupil size (blue, arbitrary units), and velocity data in degrees per second for the X-axis (violet) and Y-axis (dark gray).</p>
Full article ">Figure 16
<p>Saccade waveforms for position and velocity signals.</p>
Full article ">Figure 17
<p>Main sequence for random “participant 15” overall conditions. All three linear regressions lie close to each other.</p>
Full article ">Figure 18
<p>Trials from 13-point calibration grid split into two groups evenly. Half trials for each horizontal rm ANOVA (red) with three levels (left, middle, and right) and vertical rm ANOVA (blue) with levels (up, middle and down).</p>
Full article ">Figure 19
<p>Mean pupil sizes from the rm ANOVAs for the gaze direction groups. Error bars indicate standard error. Double asterisks indicate highly significant mean differences.</p>
Full article ">
14 pages, 5366 KiB  
Article
Investigation of Mn2+-Doped Stearic-Acid Through XRD, Raman, and FT-IR, and Thermal Studies
by Rodrigo M. Rocha, Marinaldo V. de Souza Junior, Luiz F. L. Silva, Paulo T. C. Freire, Gardênia S. Pinheiro, Waldomiro Paschoal, Francisco F. de Sousa and Sanclayton G. C. Moreira
Quantum Beam Sci. 2025, 9(1), 8; https://doi.org/10.3390/qubs9010008 (registering DOI) - 1 Mar 2025
Viewed by 128
Abstract
In this research, we investigated the influence of Mn2+ ions on the packing in stearic acid (SA) crystals through the use of Raman spectroscopy, X-ray diffraction (XRD), and Fourier transform infrared (FT-IR) spectroscopy. The crystals investigated were obtained utilizing the slow evaporation [...] Read more.
In this research, we investigated the influence of Mn2+ ions on the packing in stearic acid (SA) crystals through the use of Raman spectroscopy, X-ray diffraction (XRD), and Fourier transform infrared (FT-IR) spectroscopy. The crystals investigated were obtained utilizing the slow evaporation methodology in a hexane solution under varying manganese (Mn) concentrations sourced from MnSO4 5H2O (0.5, 1.0, 1.5, 2.0, 4.0, and 6.0%). XRD studies indicated that all SA crystals were grown in the Bm form (monoclinic), favoring the gauche conformation in molecular packing. Additionally, crystalline lattice modifications were observed through Raman spectral changes in the low-vibrational energy region. Variations in the intensities and Raman shifts in two lattice vibrational modes, centered at approximately 59 and 70 cm−1, revealed that two types of hydrogen bonds are distinctly affected within the crystalline lattice. Furthermore, the unit cell parameters (a, b, c, and β) were determined via Rietveld refinement, and their behavior was analyzed as a function of Mn concentration. The results indicated that Mn2+ ions exert a strain and deformation effect on the unit cell. Lastly, differential scanning calorimetry (DSC) was employed to evaluate the thermal stability of the Bm form of SA crystals. Full article
(This article belongs to the Section Engineering and Structural Materials)
Show Figures

Figure 1

Figure 1
<p>Organizational structure of dimers in the (<b>a</b>) all-<span class="html-italic">trans</span> and (<b>b</b>) <span class="html-italic">gauche</span> conformations for the SA crystal lattice.</p>
Full article ">Figure 2
<p>XRD patterns of Mn<sup>2+</sup>-doped SA crystals compared to the theoretical patterns of both the B<sub>m</sub> and C forms and the experimental pattern of the MnSO<sub>4</sub> 5H<sub>2</sub>O phase in the 2θ = 5–50° angular range.</p>
Full article ">Figure 3
<p>Cell parameters vs. Mn concentration plot for the SA crystals.</p>
Full article ">Figure 4
<p>Raman spectra of Mn<sup>2+</sup>-doped SA crystals compared to the spectra of MnSO<sub>4</sub> 5H<sub>2</sub>O and pure SA crystals (C form) in the ranges of (<b>a</b>) 20–300 cm<sup>−1</sup> and (<b>b</b>) 300–1000 cm<sup>−1</sup>.</p>
Full article ">Figure 5
<p>Raman spectra of Mn<sup>2+</sup>-doped SA crystals compared to the spectra of MnSO<sub>4</sub> 5H<sub>2</sub>O and pure SA crystals (C form) in the ranges of (<b>a</b>) 1000–1700 cm<sup>−1</sup> and (<b>b</b>) 2800–3000 cm<sup>−1</sup>.</p>
Full article ">Figure 6
<p>Wavenumber vs. Mn concentration plot for the lattice modes initially observed at 59 and 70 cm<sup>−1</sup>.</p>
Full article ">Figure 7
<p>IR spectra of Mn<sup>2+</sup>-doped SA crystals compared to the spectra of MnSO<sub>4</sub> 5H<sub>2</sub>O and pure SA crystals (C form) in the ranges of (<b>a</b>) 100–500 cm<sup>−1</sup> and (<b>b</b>) 500–1000 cm<sup>−1</sup>.</p>
Full article ">Figure 8
<p>IR spectra of Mn<sup>2+</sup>-doped SA crystals compared to the spectra of MnSO<sub>4</sub> 5H<sub>2</sub>O and pure SA crystals (C form) in the ranges of (<b>a</b>) 1000–1800 cm<sup>−1</sup> and (<b>b</b>) 2700–3000 cm<sup>−1</sup>.</p>
Full article ">Figure 9
<p>DSC analyses of the Mn<sup>2+</sup>-doped SA samples measured in the 300–390 K temperature interval.</p>
Full article ">
21 pages, 33924 KiB  
Article
Multiparameter Inversion of Seismic Pre-Stack Amplitude Variation with Angle Based on a New Propagation Matrix Method
by Qianlong Ding, Shuangquan Chen, Jinsong Shen, Zuzhi Hu and Guoquan Wang
Appl. Sci. 2025, 15(5), 2636; https://doi.org/10.3390/app15052636 - 28 Feb 2025
Viewed by 254
Abstract
The classical pre-stack seismic inversion technique uses the Zoeppritz equation and its simplified versions to calculate the PP and PS reflection coefficients at different incidence angles, aiding in inverting the subsurface velocity and density parameters. Despite its widespread application, the amplitude variation with [...] Read more.
The classical pre-stack seismic inversion technique uses the Zoeppritz equation and its simplified versions to calculate the PP and PS reflection coefficients at different incidence angles, aiding in inverting the subsurface velocity and density parameters. Despite its widespread application, the amplitude variation with angle (AVA) inversion based on the Zoeppritz equation has limitations regarding the accuracy. The AVA neglects transmission losses and the effects of multiple reflections during seismic wave propagation, resulting in reduced resolution. In contrast, the propagation matrix theory offers a comprehensive range of reflection coefficients for P- and S-waves in multilayered media at arbitrary incidence angles, thereby theoretically enhancing the inversion accuracy. However, the seismic responses obtained using this method exist in the slowness–frequency domain and require constant slowness for consistency along a profile. This assumption is violated when variations in the P-wave velocity occur within the subsurface, affecting the incidence angle of propagating seismic waves. This study modifies the propagation matrix theory to compute AVA seismic responses and applies it to pre-stack multiparameter inversion. The effectiveness of the modified method was validated by deriving theoretical AVA seismic responses and comparing them to solutions from a typical layered media model. The modified theory was also employed for seismic pre-stack inversion. Numerical simulations and field data tests demonstrated that the new propagation matrix method offers a high accuracy and stability. Full article
Show Figures

Figure 1

Figure 1
<p>Three-layer model with increasing velocity and density. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">P</mi> </msub> </mrow> </semantics></math>; (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">S</mi> </msub> </mrow> </semantics></math>; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>density</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 2
<p>Three-layer model with a low-velocity body. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">P</mi> </msub> </mrow> </semantics></math>; (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">S</mi> </msub> </mrow> </semantics></math>; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>density</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Three-layer model with a high-velocity body. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">P</mi> </msub> </mrow> </semantics></math>; (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">S</mi> </msub> </mrow> </semantics></math>; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>density</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Synthetic PP-wave seismogram record of <a href="#applsci-15-02636-f001" class="html-fig">Figure 1</a>. (<b>a</b>) Synthetic AVA PP records. (<b>b</b>) Synthetically modified PM PP records. (<b>c</b>) Synthetic PM PP records. (<b>d</b>) Synthetic Zoeppritz PP records.</p>
Full article ">Figure 5
<p>Synthetic PP-wave seismogram record of <a href="#applsci-15-02636-f002" class="html-fig">Figure 2</a>. (<b>a</b>) Synthetic AVA PP records. (<b>b</b>) Synthetically modified PM PP records. (<b>c</b>) Synthetic PM PP records. (<b>d</b>) Synthetic Zoeppritz PP records.</p>
Full article ">Figure 6
<p>Synthetic PP-wave seismogram record of <a href="#applsci-15-02636-f003" class="html-fig">Figure 3</a>. (<b>a</b>) Synthetic AVA PP records. (<b>b</b>) Synthetically modified PM PP records. (<b>c</b>) Synthetic PM PP records. (<b>d</b>) Synthetic Zoeppritz PP records.</p>
Full article ">Figure 7
<p>Synthetic PS-wave seismogram record of <a href="#applsci-15-02636-f001" class="html-fig">Figure 1</a>. (<b>a</b>) Synthetic AVA PS records. (<b>b</b>) Synthetically modified PM PS records. (<b>c</b>) Synthetic PM PS records. (<b>d</b>) Synthetic Zoeppritz PS records.</p>
Full article ">Figure 8
<p>Synthetic PS-wave seismogram record of <a href="#applsci-15-02636-f002" class="html-fig">Figure 2</a>. (<b>a</b>) Synthetic AVA PS records. (<b>b</b>) Synthetically modified PM PS records. (<b>c</b>) Synthetic PM PS records. (<b>d</b>) Synthetic Zoeppritz PS records.</p>
Full article ">Figure 9
<p>Synthetic PS-wave seismogram record of <a href="#applsci-15-02636-f003" class="html-fig">Figure 3</a>. (<b>a</b>) Synthetic AVA PS records. (<b>b</b>) Synthetically modified PM PS records. (<b>c</b>) Synthetic PM PS records. (<b>d</b>) Synthetic Zoeppritz PS records.</p>
Full article ">Figure 10
<p>The reflection coefficient of the first reflection interface of <a href="#applsci-15-02636-f001" class="html-fig">Figure 1</a>. (<b>a</b>) PP-wave amplitude. (<b>b</b>) PS-wave amplitude.</p>
Full article ">Figure 11
<p>The reflection coefficient of the second reflecting interface of <a href="#applsci-15-02636-f001" class="html-fig">Figure 1</a>. (<b>a</b>) PP-wave amplitude. (<b>b</b>) PS-wave amplitude.</p>
Full article ">Figure 12
<p>The reflection coefficient of the second reflecting interface of <a href="#applsci-15-02636-f002" class="html-fig">Figure 2</a>. (<b>a</b>) PP-wave amplitude. (<b>b</b>) PS-wave amplitude.</p>
Full article ">Figure 13
<p>The reflection coefficient of the second reflecting interface of <a href="#applsci-15-02636-f003" class="html-fig">Figure 3</a>. (<b>a</b>) PP-wave amplitude. (<b>b</b>) PS-wave amplitude.</p>
Full article ">Figure 14
<p>Multilayered models with complex velocity and density profiles. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">P</mi> </msub> </mrow> </semantics></math>; (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">S</mi> </msub> </mrow> </semantics></math>; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>density</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 15
<p>Synthetic PP-wave seismogram record of <a href="#applsci-15-02636-f014" class="html-fig">Figure 14</a>. (<b>a</b>) Synthetic AVA PP records. (<b>b</b>) Synthetically modified PM PP records. (<b>c</b>) Synthetic PM PP records. (<b>d</b>) Synthetic Zoeppritz PP records.</p>
Full article ">Figure 16
<p>Data error. (<b>a</b>) shows the new propagation matrix to calculate the error of the seismic response. (<b>b</b>) shows the error of the seismic response calculated by the traditional propagation matrix. (<b>c</b>) shows the error value of the seismic response calculated by the Zoeppritz equation.</p>
Full article ">Figure 17
<p>Two-dimensional real model. (<b>a</b>) The P-wave velocity model, (<b>b</b>) S-wave velocity model, and (<b>c</b>) density model.</p>
Full article ">Figure 18
<p>Initial model. Two-dimensional real model. (<b>a</b>) The P-wave velocity model, (<b>b</b>) S-wave velocity model, and (<b>c</b>) density model.</p>
Full article ">Figure 19
<p>Two-dimensional inversion results (SNR = 2). (<b>a</b>) The P-wave velocity model, (<b>b</b>) S-wave velocity model, and (<b>c</b>) density model.</p>
Full article ">Figure 20
<p>One-dimensional inversion results. (<b>a</b>–<b>c</b>) show the model data. (<b>d</b>–<b>f</b>) show the absolute error value between the different methods and the real logging data. The black curve represents the logging data. The green curve represents the initial model. The blue curve denotes the inversion result obtained using the Zoeppritz equation. The red curve is the inversion result obtained using the new propagation matrix method.</p>
Full article ">Figure 21
<p>Field data. (<b>a</b>–<b>c</b>) represent the field seismic data for incidence angles of 7°, 13°, and 19°, respectively.</p>
Full article ">Figure 22
<p>Inversion profile of the field data. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">P</mi> </msub> </mrow> </semantics></math>; (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">S</mi> </msub> </mrow> </semantics></math>; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>density</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 23
<p>Field data inversion results. (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">P</mi> </msub> </mrow> </semantics></math>; (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi>V</mi> <mi mathvariant="normal">S</mi> </msub> </mrow> </semantics></math>; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>density</mi> </mrow> </semantics></math>. The black curve represents the logging data and the green curve represents the initial model. The blue curve is the inversion result obtained using the Zoeppritz equation. The red curve is the inversion result obtained using the new propagation matrix method. (<b>d</b>) illustrates the lithology map, in which blue denotes sandstone, yellow denotes mudstone, and green indicates a mixture of sandstone and mudstone.</p>
Full article ">
22 pages, 5673 KiB  
Article
Effects of Sensor Speed and Height on Proximal Canopy Reflectance Data Variation for Rice Vegetation Monitoring
by Md Rejaul Karim, Md Asrakul Haque, Shahriar Ahmed, Md Nasim Reza, Kyung-Do Lee, Yeong Ho Kang and Sun-Ok Chung
Agronomy 2025, 15(3), 618; https://doi.org/10.3390/agronomy15030618 (registering DOI) - 28 Feb 2025
Viewed by 126
Abstract
Sensing distance and speed have crucial effects on the data of active and passive sensors, providing valuable information relevant to crop growth monitoring and environmental conditions. The objective of this study was to evaluate the effects of sensing speed and sensor height on [...] Read more.
Sensing distance and speed have crucial effects on the data of active and passive sensors, providing valuable information relevant to crop growth monitoring and environmental conditions. The objective of this study was to evaluate the effects of sensing speed and sensor height on the variation in proximal canopy reflectance data to improve rice vegetation monitoring. Data were collected from a rice field using active and passive sensors with calibration procedures including downwelling light sensor (DLS) calibration, field of view (FOV) alignment, and radiometric calibration, which were conducted per official guidelines. The data were collected at six sensor heights (30–130 cm) and speeds (0–0.5 ms–1). Analyses, including peak signal-to-noise ratio (PSNR) and normalized difference vegetation index (NDVI) calculations and statistical assessments, were conducted to explore the impacts of these parameters on reflectance data variation. PSNR analysis was performed on passive sensor image data to evaluate image data variation under varying data collection conditions. Statistical analysis was conducted to assess the effects of sensor speed and height on the NDVI derived from active and passive sensor data. The PSNR analysis confirmed that there were significant impacts on data variation for passive sensors, with the NIR and G bands showing higher noise sensitivity at increased speeds. The NDVI analysis showed consistent patterns at sensor heights of 70–110 cm and sensing speeds of 0–0.3 ms–1. Increased sensing speeds (0.4–0.5 ms–1) introduced motion-related variability, while lower heights (30–50 cm) heightened ground interference. An analysis of variance (ANOVA) indicated significant individual effects of speed and height on four spectral bands, red (R), green (G), blue (B), and near-infrared (NIR), in the passive sensor images, with non-significant interaction effects observed on the red edge (RE) band. The analysis revealed that sensing speed and sensor height influence NDVI reliability, with the configurations of 70–110 cm height and 0.1–0.3 ms–1 speed ensuring the stability of NDVI measurements. This study notes the importance of optimizing sensor height and sensing speed for precise vegetation index calculations during field data acquisition for agricultural crop monitoring. Full article
(This article belongs to the Section Precision and Digital Agriculture)
Show Figures

Figure 1

Figure 1
<p>The experimental rice field for data acquisition: (<b>a</b>) the location of the experimental rice field, red color indicates the experimental rice field and the pink sections are the plot used in this experiment, and (<b>b</b>) the rice field condition during the data collection period.</p>
Full article ">Figure 2
<p>The custom portable aluminum structure for data acquisition along with the active sensor, passive sensor, and GPS.</p>
Full article ">Figure 3
<p>DLS calibration procedure: (<b>a</b>) radiometric calibration process using reference panel and DLS, and (<b>b</b>) magnetometer calibration process for passive sensor.</p>
Full article ">Figure 4
<p>Overall workflow of data acquisition and processing for data from GPS, active sensor, and passive sensor.</p>
Full article ">Figure 5
<p>The FOVs of the active and passive sensors for the sensor height. The FOV calculations were based on the horizontal angular coverage with both sensors mounted at varying heights above the crop canopy.</p>
Full article ">Figure 6
<p>The NDVI patterns estimated using an active sensor at varying sensor heights (30–130 cm above the canopy) and data collection speeds (0–0.5 ms<sup>−1</sup>): (<b>a</b>) 30 cm, (<b>b</b>) 50 cm, (<b>c</b>) 70 cm, (<b>d</b>) 90 cm, (<b>e</b>) 110 cm, and (<b>f</b>) 130 cm. Each data point represents the mean NDVI for individual grid plots, with different colors corresponding to different speeds.</p>
Full article ">Figure 7
<p>The effects of speed and height on NDVI measurements using an active sensor: (<b>a</b>) NDVI variation across different heights (30, 50, 70, 90, 110, and 130 cm above the canopy) at varying speeds, and (<b>b</b>) NDVI variation across different speeds (0–0.5 ms<sup>−1</sup>) at varying heights. Each data point represents the mean NDVI for individual grid plots, with colors indicating different (<b>a</b>) sensor heights and (<b>b</b>) sensing speeds.</p>
Full article ">Figure 8
<p>The NDVI patterns estimated using a passive sensor at varying sensor heights (30–130 cm above the canopy) and data collection speeds (0–0.5 ms<sup>−1</sup>): (<b>a</b>) 30 cm, (<b>b</b>) 50 cm, (<b>c</b>) 70 cm, (<b>d</b>) 90 cm, (<b>e</b>) 110 cm, and (<b>f</b>) 130 cm. Each data point represents the mean NDVI for individual grid plots, with different colors corresponding to different speeds.</p>
Full article ">Figure 9
<p>The effects of speed and height on NDVI measurements using a passive sensor: (<b>a</b>) NDVI variation across different heights (30, 50, 70, 90, 110, and 130 cm above the canopy) at varying speeds, and (<b>b</b>) NDVI variation across different speeds (0–0.5 ms<sup>−1</sup>) at varying heights. Each data point represents the mean NDVI for individual grid plots, with colors indicating different heights (<b>a</b>) and speeds (<b>b</b>).</p>
Full article ">
18 pages, 3194 KiB  
Article
Ritz Solution of Vibration Analysis of Functionally Graded Porous Elliptic Shells and Panels Under Various Arbitrary Boundary Types
by Qingtao Gong, Tao Liu, Yao Teng, Binjie Ma and Xin Li
Materials 2025, 18(5), 1101; https://doi.org/10.3390/ma18051101 - 28 Feb 2025
Viewed by 107
Abstract
This paper seeks to establish a generalized numerical model to examine the free vibration behavior of functionally graded porous (FGP) elliptical shells and panels with various boundary types. The model is built on first-order shear deformation theory (FSDT) to express structural displacements. A [...] Read more.
This paper seeks to establish a generalized numerical model to examine the free vibration behavior of functionally graded porous (FGP) elliptical shells and panels with various boundary types. The model is built on first-order shear deformation theory (FSDT) to express structural displacements. A segmentation technique is used to maintain continuity between shell elements, and virtual spring boundary techniques are employed to simulate arbitrary boundaries. Variable-coefficient Jacobi polynomials are introduced as admissible functions for displacement. Finally, the Ritz variational method, combined with the least-squares weighted residual method (LSWRM), is used for constructing the energy functional and solving the energy equations. Validation of the numerical model against finite element and literature results confirms its reliability and convergence properties. This study also explores the effects of geometric parameters and boundary conditions on FG elliptical shells and panels, providing a theoretical basis for future research. Full article
(This article belongs to the Special Issue Numerical Analysis of Sandwich and Laminated Composites)
Show Figures

Figure 1

Figure 1
<p>Structural schematic diagram and material cross-section diagram of FGP elliptical shells and panels.</p>
Full article ">Figure 2
<p>Variation in FGP structure frequencies with different segment numbers: (<b>a</b>) elliptic shells; (<b>b</b>) elliptic panels.</p>
Full article ">Figure 3
<p>Variation in FGP structure frequencies with different truncation terms: (<b>a</b>) elliptic shells; (<b>b</b>) elliptic panels.</p>
Full article ">Figure 4
<p>Vibration modes of FGP elliptic shells under various boundary types.</p>
Full article ">Figure 5
<p>Vibration modes of FGP elliptic panels under various boundary types.</p>
Full article ">Figure 6
<p>Variation in ΔΩ versus the porosity coefficients <span class="html-italic">e</span><sub>0</sub> for functionally graded porous elliptic shells: (<b>a</b>) C-C; (<b>b</b>) C-F; (<b>c</b>) E1-E1; (<b>d</b>) E4-E4.</p>
Full article ">Figure 6 Cont.
<p>Variation in ΔΩ versus the porosity coefficients <span class="html-italic">e</span><sub>0</sub> for functionally graded porous elliptic shells: (<b>a</b>) C-C; (<b>b</b>) C-F; (<b>c</b>) E1-E1; (<b>d</b>) E4-E4.</p>
Full article ">Figure 7
<p>Variation in ΔΩ versus the porosity coefficients <span class="html-italic">e</span><sub>0</sub> for functionally graded porous elliptic panels: (<b>a</b>) CCCC; (<b>b</b>) CFCF; (<b>c</b>) E1E1E1E1; (<b>d</b>) E4E4E4E4.</p>
Full article ">Figure 8
<p>Variation in the frequency parameters ΔΩ versus the circumferential angle <span class="html-italic">θ</span><sub>0</sub> for functionally graded porous elliptic panels: (<b>a</b>) CCCC boundary; (<b>b</b>) CFCF boundary; (<b>c</b>) E1E1E1E1 boundary; (<b>d</b>) E4E4E4E4 boundary.</p>
Full article ">
Back to TopTop