The Citarum Watershed is indeed a critical water resource in Indonesia, playing a significant role in providing water to Jakarta and other areas in West Java. However, it faces severe environmental stress due to land use changes and climate changes. The Upper Citarum Watershed, considered to be a conservation area, has experienced rapid development due to human activities and economic growth. Climate change not only affects the rainfall value but also the rainfall patterns and sediment flow. The sedimentation process significantly impacts the soil characteristics around the river body. Several factors such as topography, flow velocity, and soil texture influence the sediment characteristics. Given the critical condition of climate and land use change, this study aims to analyse the impacts of the hydrological response of land use and climate change on the sediment rate in the Upper Citarum Watershed. The land use change analysis was conducted by comparing the land use data in 2000, 2010, and 2023 in the Upper Citarum Watershed. The deposition process of solid particles such as sand, silt, and gravel that are transported in the Upper Citarum River were examined in a soil investigation. The sediment rate and deposition by river flow were analysed using HEC-RAS quasi-unsteady flow. The impact of climate change in this study was assessed by simulating the discharge in three conditions, with the first simulation using the discharge from 2000 to 2010, the second simulation using the discharge from 2011 to 2023, and the last simulation using the discharge from 2000 to 2023. Due to the land use change, the developed area increased from 4% to 24% between 2000 until 2023. The magnitude of low flow during the simulation step for three discharge gauges (Majalaya, Dayeuhkolot, and Nanjung) decreased up to 48%, but, on other hand, the sediment rate increased by 20% in Dayeuhkolot.

Coastal defense structures play a crucial role in mitigating wave impacts; yet, existing breakwater designs often face challenges in balancing wave reflection, energy dissipation, and structural stability. This study leverages machine learning (ML) to predict the optimal 2D dimensions of rectangular breakwaters in two configurations: submerged at the bottom of a wave tank and positioned at the free surface. Further, the objective is to achieve controlled wave reflection allowing a specific wave run-up and optimized energy dissipation, while ensuring maritime stability. Thus, we used an analytical equation modeling the reflection coefficient versus relative water depth (KH), for different immersion ratios of obstacle (h/H), and relative length (l/H). Two datasets of 32,000 data points were generated for underwater and free-surface breakwaters, with an additional 10,000 data points for validation, totaling 42,000 data points per case. Five ML algorithms—Random Forest, Support Vector Regression, Artificial Neural Network, Decision Tree, and Gaussian Process—were applied and evaluated. Results demonstrated that Random Forest and Decision Tree balanced accuracy with computational efficiency, while the Gaussian Process closely matched analytical results but demanded higher computational resources. These findings support ML as a powerful tool to optimize breakwater design, complementing traditional methods and contributing to more sustainable and resilient coastal defense systems.

Although there do exist a few aeroacoustic studies on harmful artificial phenomena related to the usage of non-uniform Cartesian grids in lattice Boltzmann methods (LBM), a thorough quantitative comparison between different categories of grid arrangement is still missing in the literature. In this paper, several established schemes for hierarchical grid refinement in lattice Boltzmann simulations are analyzed with respect to spurious aeroacoustic emissions using a weakly compressible model based on a D3Q19 athermal velocity set. In order to distinguish between various sources of spurious phenomena, we deploy both the classical Bhatnagar–Gross–Krook and other more recent collision models like the hybrid recursive-regularization operator, the latter of which is able to filter out detrimental non-hydrodynamic mode contributions, inherently present in the LBM dynamics. We show by means of various benchmark simulations that a cell-centered approach, either with a linear or uniform explosion procedure, as well as a vertex-centered direct-coupling method, proves to be the most suitable with regards to aeroacoustics, as they produce the least amount of spurious noise. Furthermore, it is demonstrated how simple modifications in the selection of distribution functions to be reconstructed during the communication step between fine and coarse grids affect spurious aeroacoustic artifacts in vertex-centered schemes and can thus be leveraged to positively influence stability and accuracy.

Fine bubbles (FBs) are defined by the ISO/TC 281 as gas bubbles with a diameter of less than 100 μm, and they have interesting properties such as high surface-to-volume ratio, low buoyancy, long residence time, electric charge, and self-pressurization effect. Typically, FBs are characterized in terms of size distribution, concentration, and zeta potential through specialized microscopic and nanoscopic measuring devices. This work proposes a multi-objective optimization problem to find the optimal conditions to generate FBs from experimental macroscopic measurements in terms of dissolved oxygen (DO). Then, detailed microscopic measurements in terms of size distribution and zeta potential are conducted. Additionally, two venturi-type Fine Bubble Generators (FBGs) were 3D-printed in-house to evaluate the relationship between the internal structure and the generation of FBs. The best FBGs have an obstacle in the diverging section that promotes FB generation under the evaluated experimental conditions. Under the best operating conditions, FBs were stable over 7 days with a size distribution between 60 and 90 nm and with an average of −21 mV.

The stability of time-dependent compressible linear flows, which are characterized by periodic variations in either their shape or their shear, is investigated. Two novel parametric instabilities are found: an instability that occurs for periodically wobbling elliptic vortices at a number of discrete oscillation frequencies that are proportional to the Mach number and an instability that occurs for all linear flows at various frequencies of the shear oscillation that depend on the Mach number. In addition, the physical mechanism underlying the instabilities is explained in terms of the linear interaction of three waves with time-varying wavevectors that describe the evolution of perturbations: a vorticity wave representing the evolution of incompressible perturbations and two counter-propagating acoustic waves. Elliptical instability occurs because the scale of the acoustic waves decreases exponentially and their wave action is conserved, leading to an exponential increase in the acoustic waves’ energies. The instability in shear-varying flows is driven by the interaction between vorticity and the acoustic waves, which couple through the shear and for specific frequencies resonate parametrically, leading to exponential or linear growth.

The flow occurring at the bottom of a polycrystalline diamond compact (PDC) drill bit involves a complex process made up of drilling fluid and the drilled rock cuttings. A thorough understanding of the bottom-hole flow conditions is essential for accurately evaluating and optimizing the hydraulic structure design of the PDC drill bit. Based on a comprehensive understanding of the hydraulic structure and fluid flow characteristics of PDC drill bits, this study integrates computational fluid dynamics (CFD) with rock-breaking simulation methods to refine and enhance the numerical simulation approach for the liquid–solid two-phase flow field of PDC drill bits. This study further conducts a comparative analysis of simulation results between single-phase and liquid–solid two-phase flows, highlighting the influence of rock cuttings on flow dynamics. The results reveal substantial differences in flow behavior between single-phase and two-phase conditions, with rock cuttings altering the velocity distribution, flow patterns, and hydraulic performance near the bottom-hole region of the drill bit. The two-phase flow simulation results demonstrate higher accuracy and provide a more detailed depiction of the bottom-hole flow, facilitating the identification of previously unrecognized issues in the hydraulic structure design. These findings advance the methodology for multiphase flow simulation in PDC drill bit studies, providing significant academic and engineering value by offering actionable insights for optimizing hydraulic structures and extending bit life.

In this study, the influence of the presence of a Newtonian solvent on the flow of a Giesekus fluid in a plane channel or fracture is investigated with a focus on the determination of the flow rate for an assigned external pressure gradient. The pressure field is nonlinear due to the presence of the normal transverse stress component. As expected, the flow rate per unit width Q′ is larger than for a Newtonian fluid and decreases as the solvent increases. It is strongly dependent on the viscosity ratio ε (0≤ε≤1), the dimensionless mobility parameter β (0≤β≤1) and the Deborah number De, the dimensionless driving pressure gradient. The degree of dependency is notably strong in the low range of ε. Furthermore, Q′ increases with De and tends to a constant asymptotic value for large De, subject to the limitation of laminar flow. When the mobility factor β is in the range 0.5÷1, there is a minimum value of ε to obtain an assigned value of De. The ratio UN/U between Newtonian and actual mean velocity depends only on the product βDe, as for other non-Newtonian fluids.

The aquaculture industry encounters substantial obstacles, including organic pollution, oxygen insufficiency, and elevated levels of ammonia and carbon dioxide. Aeration systems are employed to enhance the process of oxygen transfer and promote circulation. The Lily impeller, a newly developed technology, has demonstrated reduced energy consumption in comparison to conventional impeller designs. The objective of this study is to examine how changes in the cut sweep ratio impact the distribution of dissolved oxygen in shrimp ponds, using computational fluid dynamics (CFD) simulation. A user-defined function (UDF) was utilized to incorporate a dissolved oxygen model into the pond. Five designs of Lily impellers were analyzed and compared with each other. This study demonstrated that alterations in the cut sweep ratio significantly affected the distribution of dissolved oxygen, dynamic pressure, and flow velocity in the pond. The “no cut” variant exhibited the highest average dissolved oxygen value of 0.00385 kg/m3, along with a maximum dynamic pressure of 11.5 Pa and a maximum flow velocity of 0.96 m/s, resulting in the most significant outcomes. This study determined that only the immediate area surrounding the aerator possesses dissolved oxygen levels that are sufficiently elevated to support the survival of shrimp. Consequently, the installation of additional aerators is necessary to guarantee the presence of adequate dissolved oxygen throughout the entire pond.

The Tesla valve is a fluidic diode that enables unidirectional flow while impeding the reverse flow without the assistance of any moving parts. Conventional Tesla valves share a distinctive feature of a bifurcated section that connects the inlet and outlet. This study uses computational fluid dynamic (CFD) simulations to analyze the importance of the bifurcated design to the efficiency of the Tesla valve, quantified by diodicity. Simulations over the range of the Reynolds number, Re = 50–2000, are performed for three designs: the T45-R, D-valve, and GMF valve, each with two versions with and without the bifurcated section. For the T45-R valve, removing the bifurcated section leads to a consistent increase in diodicity, particularly at high Re. In contrast, the diodicity of the GMF valve drops significantly when the bifurcated section is removed. The D-valve exhibits a mixed behavior. Without the bifurcated section, its diodicity is suppressed at low Re but begins to increase for Re > 1100, eventually matching the diodicity of the bifurcated version at Re = 2000. The results highlight the intricate relationship between valve geometry and efficiency of Tesla-type valves and the dependence of this relationship on the Reynolds number.

Quasi-direct current (Q-DC) filamentary electrical discharges are used to control the shock train in a back-pressured Mach 2 duct flow. The coupled interaction between the plasma filaments and the shock train leading edge (STLE) is studied for a variety of boundary conditions. Electrical parameters associated with the discharge are recorded during actuation, demonstrating a close correlation between the STLE position and dynamics. High-speed self-aligned focusing schlieren (SAFS) and high frame-rate color camera imaging are the primary optical diagnostics used to study the flowfield and plasma morphology. Shock tracking and plasma characterization algorithms are employed to extract time-resolved quantitative data during shock–plasma interactions. Four distinct shock–plasma interaction types are identified and outlined, revealing a strong dependence on the spacing between the uncontrolled STLE and discharge electrodes and a moderate dependence on flow parameters.