Three-dimensional Numerical Model Research Articles

Cyclic performance of RC frame joints with drilled openings

The construction industry commonly employs reinforced concrete (RC) beams with utility pipe openings to facilitate the passage of pipes through structures. These openings can be made through pre-construction or post-construction processes. Pre-construction openings in RC members are included in the design, and in most cases, designers include reinforcements around the opening to deal with the stress concentration and failure mechanism caused by the presence of the opening. However, post-construction or drilled openings are added after the structure is built, potentially disrupting established load paths, and developing stress concentrations that are not anticipated by the designer. In this study, three-dimensional nonlinear numerical finite element models (FEM) models were developed and validated according to four sets of experimental data. These models formed the basis for a parametric study exploring different design parameters. Parameters such as the width and depth of the specimen and the size of the opening were investigated. The results suggest that the presence of the opening leads to a brittle behavior and a sudden drop in the strength and stiffness of the specimen after the ultimate load. Increasing the depth of the beam results in more brittle behavior in specimens with openings and decreased the drift ratio where the specimen failed by up to 75%. However, increasing the width of the specimen increases the ductility of these specimens. Finally, opening size was found to lead the specimen to behave more brittle in all cases, and decrease the maximum drift ratio from 5% to 1.5%, and decrease the load capacity in 1.5% drift ratio by up to 61%. This study will help engineers and designers to understand the effect of beam-column joints with post-construction openings under cyclic loading, whether these openings are intentional or by human error.

Experimental and Numerical Investigation of Single-Phase Liquid Cooling for Heterogeneous Integration Multichip Module

Abstract The variety of new electronic packaging technologies has grown significantly over the last 20 years as a result of market demands for higher device performance at lower costs and in less space. Those demands have pushed for heterogeneous packaging, where computer chips with different stack heights are closely packed, creating nonuniform heat flux and temperature and additional challenges for thermal management. Without implementing an appropriate thermal management strategy for heterogeneous packages, large temperature gradients can be observed within the package, which would increase the thermal stresses on the chip and raise reliability issues. To mimic this real-life scenario of such packaging, an experimental setup was designed and built. The design of the new experimental setup consists of four identical 1.2 cm × 1.2 cm ceramic heaters, each of which is connected to a separate power supply and can reach a heat flux of 140 W/cm2. Accordingly, this mock package is capable of delivering different power levels to mimic different multicore microprocessor conditions. To give the heater the ability to move precisely in the x-, y-, and z-directions, each heater is mounted to an XYZ linear stage. Deionized water (DI) was used as the working fluid, and a pin-fin heat sink was used to run the initial steady-state tests on the experimental rig. The tests showed how different flow rates at a constant fluid temperature and input power affect the temperatures of the heaters and the thermohydraulic performance of the heat sink. In addition, a three-dimensional numerical model has been developed and validated with experimental data in terms of heat sink pressure drop and the temperatures of the heaters.

Research on the Correction of Soil Pressure Distribution Pattern in Circular Working Wells in Soft Soil Areas

This paper investigates the distribution characteristics of earth pressure on circular caisson foundations in soft soil areas based on practical engineering considerations. Utilizing a three-dimensional numerical analysis model and a displacement-based method for earth pressure calculation, the lateral earth pressure on the retaining structure is examined. Adjustments are made to the distribution pattern of earth pressure on circular caisson foundations in soft soil areas. The study demonstrates that the radius of the caisson significantly influences the spatial distribution of soil pressure. Additionally, empirical formulas for the depth and radius of the reduction in soil pressure around circular caisson foundations are obtained through the research. Comparative analysis reveals that the modified earth pressure planar foundation beam method, which accounts for the arching effect, is simpler, features more mature parameter selection, and is more cost-effective than the three-dimensional finite element method. In comparison to the standard planar elastic foundation beam method for earth pressure, it aligns more closely with the actual direction of displacement.

Earthquake source impacts on the generation and propagation of seismic infrasound to the upper atmosphere

Summary Earthquakes with moment magnitude (Mw) ranging from 6.5 to 7.0 have been observed to generate sufficiently strong acoustic waves (AWs) in the upper atmosphere. These AWs are detectable in Global Navigation Satellite System satellite signals-based total electron content (TEC) observations in the ionosphere at altitudes ∼250-300 km. However, the specific earthquake source parameters that influence the detectability and characteristics of AWs are not comprehensively understood. Here, we extend our approach of coupled earthquake-atmosphere dynamics modeling by combing dynamic rupture and seismic wave propagation simulations with two- and three-dimensional atmospheric numerical models, to investigate how the characteristics of earthquakes impact the generation and propagation of AWs. We developed a set of idealized dynamic rupture models varying faulting types and fault sizes, hypocentral depths and stress drops. We focus on earthquakes of Mw 6.0-6.5, which are considered the smallest detectable with TEC, and find that the resulting AWs undergo nonlinear evolution and form acoustic shock N-waves reaching thermosphere at ∼90-140 km. The results reveal that the magnitude of the earthquakes is not the sole or primary factor determining the amplitudes of AWs in the upper atmosphere. Instead, various earthquake source characteristics, including the direction of rupture propagation, the polarity of seismic wave imprints on the surface, earthquake mechanism, stress drop, and radiated energy, significantly influence the amplitudes and periods of AWs. The simulation results are also compared with observed TEC fluctuations from AWs generated by the 2023 Mw 6.2 Suzu (Japan) earthquake, finding preliminary agreement in terms of model-predicted signal periods and amplitudes. Understanding these nuanced relationships between earthquake source parameters and AW characteristics is essential for refining our ability to detect and interpret AW signals in the ionosphere.

Effects of environmental humidity on the gap of mortise-tenon joints: experimental and numerical investigations

Mortise-tenon (MT) joints play a vital role in the lateral deformability and seismic performance of traditional timber structures. However, moisture exposure can cause shrinkage-induced gaps and then affect the durability of MT joints. Limited knowledge exists on the humidity response of MT joints. To investigate the effect of environmental humidity on the gaps between mortise and tenon, multiple measurements of moisture content (MC) and shrinkage deformation were conducted on straight MT joints that were subjected to humidity variations. Different ambient conditions and exposure levels were taken into account. Three-dimensional numerical models were established to simulate the humidity response of MT joints. Model validation was then performed by comparing the simulation and experimental results in terms of moisture distribution and gaps. On this basis, quantitative analysis was made to reveal the development mechanism and main influence factors of moisture-induced gaps. Results show that straight MT joints are vulnerable to environmental humidity changes. The observed gap is estimated to cause a decrease of up to 23.1% in the initial elastic stiffness. Although only contributed by the shrinkage deformation of tenons, the vertical gap is theoretically larger than the horizontal one which is attributed to a combination of both mortise and tenon. A clear linear relationship is found between the vertical gaps and the MC decrease, with the product of the radial shrinkage coefficient and tenon height as the theoretical ratio. Controlling the moisture exposure is considered the most effective way to reduce moisture-induced gaps in MT joints.

Research on the kinetics characteristics and in-cylinder combustion of opposed rotary piston engine for different power output modes

As a new type of internal combustion engine, opposed rotary piston (ORP) engine has the advantages of simple structure and high working frequency, and can achieve high power density. Consequently, it stands as an ideal power source for hybrid power systems and unmanned aerial vehicles. The ORP engine has different power output modes, corresponding to different piston rotation and volume evolutions of combustion chamber. However, the effect of power output modes on kinetics characteristics and in-cylinder combustion of ORP engines are still uncovered. In this paper, the kinematic model of an ORP engine is established for scenarios involving power output by different shafts (Shaft 1 and Shaft 3). Meanwhile, the piston rotation patterns and cylinder volume variations are investigated correspondingly. Further, three-dimensional numerical model of the ORP engine is established, enabling analysis of in-cylinder combustion characteristics, engine performance, as well as nitrogen oxides (NOx) emissions. The results indicate that the volume variations of all cylinders follow a consistent pattern under the power output by Shaft 3. A pair of opposed cylinders features longer durations of intake and expansion strokes, and shorter durations of compression and exhaust strokes, the other pair of combustion chambers is the opposite under the power output by Shaft 1. The engine demonstrates a charging efficiency of 95.36 %, with corresponding indicated thermal efficiency and indicated power output of 38.08 % and 40.9 kW, respectively. The two adjacent cylinders present significantly different operation processes in the case of the power output by Shaft 1, achieving charging efficiencies of 94.29 % and 89.11 %, respectively, with the indicated thermal efficiency of 34.49 % and 31.60 %. The corresponding minimum indicated specific fuel consumption under the two power output modes are 211.7 g/kW∙h and 243.7 g/kW∙h, respectively, accompanied by NOx emission factors of 25.3 g/kW∙h and 17.6 g/kW∙h. With the given range of ignition timing, a trade-off relationship exists between the engine performance and NOx emissions. The power density of the case of power output by Shaft 3 is superior to that of Shaft 1.

Large-scale engineered meta-barriers for attenuation of seismic surface waves

This study evaluates the performance of large-scale meta-barriers in attenuating seismic surface waves and transportation-induced vibrations in the most disastrous frequency range of 0–30 Hz. By developing analytical and three-dimensional numerical models, it evaluates the properties of the proposed meta-barrier and presents their effects on the attenuation range. Moreover, it discusses a simplified mass-spring-mass model of the meta-barrier, demonstrating that the attenuation range of a unit cell can be approximated using this model. The proposed design consists of three-dimensional two-layers (3D-2C) unit cells arranged periodically with different configurations to be placed between the wave source and structure. Specifically, the meta-barriers can be arranged gradually (periodic spacing along the length), widening the attenuation range to block nearly 85% of the most devastating waves using typically arranged stacked unit cells and triangular-like arrangement (double wedge). The results show that a typical stacked unit cells reflect the incoming waves towards the source, while using the double wedge design, the surface wave can be converted to bulk wave, thereby reducing the risk to the surroundings. Moreover, the study validates the attenuation range using three-dimensional time history analysis subjected to an artificial wavelet (Ormsby wavelet) and three actual seismic records (1975 Oroville, 1984 Bishop, and 2001 Anza time histories).

Numerical Analysis on Performance Improvement of a Vertical Plate Indirect Evaporative Cooler with Baffles

The performance of the Plate Indirect Evaporative Cooler (PIEC) can be effectively improved by incorporating baffles in the dry channel. However, in the dimensional influence of the baffles on PIEC performance there remains a research gap. In order to investigate the impact of baffle dimensions on the wet bulb efficiency, namely the average heat transfer coefficient and the cooling capacity of the PIEC, this paper proposed and verified a three-dimensional numerical model and method based on the species transport model and the Euler wall film model. At the same time, in order to obtain the equilibrium point between the enhanced heat transfer performance and the additional resistance induced by baffles, a comprehensive performance evaluation index is introduced. The results indicate that, under the same conditions, (1) the baffle effect on PIEC performance is significant at a lower inlet air velocity, and the wet bulb efficiency of the PIEC with baffles can be improved by 22.8%; (2) the baffle effect on PIEC performance is negative if its relative length exceeds 60% or the primary air inlet velocity surpasses 4 m/s under the conditions specified in this paper; and (3) the baffle effect on PIEC performance is significant when its channel height is lower and its channel width is larger, and the wet bulb efficiency of the PIEC with baffles can be improved by 29.3%.

Three-Dimensional Finite Element Modeling of Ultrasonic Vibration-Assisted Milling of the Nomex Honeycomb Structure

Machining of Nomex honeycomb composite (NHC) structures is of critical importance in manufacturing parts to the specifications required in the aerospace industry. However, the special characteristics of the Nomex honeycomb structure, including its composite nature and complex geometry, require a specific machining approach to avoid cutting defects and ensure optimal surface quality. To overcome this problem, this research suggests the adoption of RUM technology, which involves the application of ultrasonic vibrations following the axis of revolution of the UCK cutting tool. To achieve this objective, a three-dimensional finite element numerical model of Nomex honeycomb structure machining is developed with the Abaqus/Explicit software, 2017 version. Based on this model, this research examines the impact of vibration amplitude on the machinability of this kind of structure, including cutting force components, stress and strain distribution, and surface quality as well as the size of the chips. In conclusion, the results highlight that the use of ultrasonic vibrations results in an important reduction in the components of the cutting force by up to 42%, improves the quality of the surface, and decreases the size of the chips.

The effects of damming and dam regulation on a river–lake-aquifer system: 3D groundwater flow modeling of Poyang Lake (China)

Large numbers of dams are planned or constructed for hydroelectric power production, water supply, and flood control in many rivers around the world. While the social and environmental impacts of dams are now well-understood, there is little knowledge about the effects of dams on groundwater-surface water interactions in a complex lake–river-aquifer system. This study adopts a three-dimensional groundwater numerical model to quantify the potential impacts of dam construction, on regional groundwater flow system, exemplified by the proposed Poyang Lake Hydraulic Dam (PLHD) in China. Once implemented, the PLHD will regulate lake levels, artificially controlling the lake inundation area during the recession season (September-October) and the dry season (November-December). Modeling results indicated that the increase in lake infiltration induced by PLHD regulation would alter the regional water exchange pattern. The PLHD regulation will likely enhance the groundwater storage by 20.44 % in the general dry year (i.e., 2018), alleviating 78.1 % of the groundwater deficit in the extreme dry year (i.e., 2019). Additionally, the response of groundwater dynamics to the PLHD regulation exhibits significant spatial differences between the upstream and downstream areas of the dam, primarily depending on topographical changes. Overall, the PLHD regulation is most likely to impact groundwater levels across an area of approximately 5,700 km2, and the associated groundwater storage is anticipated to reach a stable state, under future long-term regulation. Our findings may guide policy-makers wishing to develop more suitable PLHD scheduling plans and floodplain protection strategies by providing reliable information regarding groundwater storage change and ecosystem succession.

A Novel Numerical Model Considering Unsaturated Soil Properties and Computational Study on Heat and Moisture Transfer Characteristics of Helix Ground Heat Exchanger

Abstract A novel three-dimensional numerical simulation model for the helix ground heat exchanger was proposed in this paper, which takes into account unsaturated soil properties. This model is more suitable for real working conditions. To validate its accuracy, a miniature model heat transfer experimental platform was constructed. Additionally, the study conducted simulation research using three types of soil with significantly different thermal and moisture characteristics. Moreover, the comprehensive thermal conductivity and water diffusion coefficient of these three soil types were determined by relevant literature and experimental tests. The aim was to comprehensively explore the impact of different soil types on the heat and mass transfer of the helix ground heat exchanger. The results indicate that the numerical model developed in this paper accurately captures the heat and mass transfer characteristics of the helix ground heat exchanger to a certain extent. Increasing the comprehensive thermal conductivity and water diffusion coefficient of the soil can significantly enhance the heat exchange capacity of the exchanger. For instance, under sandy loam conditions, the heat exchange capacity is approximately 20.73% higher compared to clay loam conditions. The study also identifies two distinct areas around the helix ground heat exchanger: the severe change region and the soft change region. In the severe change region, there is a notable decrease in soil water content near the exchanger, which inevitably weakens the thermal conductivity of the soil. It is advised to minimize this effect through measures like active water spraying.

Triaxial behavior of unreinforced and geogrid-reinforced calcareous gravelly sand: Experiment and discrete element modelling

Geogrid reinforcement is a feasible technique for improving the mechanical and deformation behaviors of calcareous gravelly sand that is predominantly distributed in reclaimed foundations as construction materials. However, the shear behavior of geogrid-reinforced calcareous gravelly sand still remains unclear. In this paper, a series of large-scale drained consolidated triaxial compression tests (the specimen 300 mm in diameter by 600 mm in height) were conducted to investigate the effects of confining pressure, moisture content and geogrid inclusion on the mechanical and deformation behaviors of calcareous gravelly sand. The dilation effect, strain-softening behavior, reinforcing effect and particle breakage characteristics of unreinforced and geogrid-reinforced calcareous gravelly sand were comprehensively analyzed and discussed. Additional analysis was also performed from a microscopic perspective using the three-dimensional discrete numerical modelling. The experimental results indicate that calcareous gravelly sand exhibited a greater resistance to strain-softening tendency, and produced a higher susceptibility to crushing under the same input energy than calcareous coarse sand. Geogrid inclusion can effectively mitigate the susceptibility of calcareous gravelly sand to volumetric dilation through constraining the lateral expansion within the reinforced zone, and thus substantially increased the peak shear strength. Moreover, particle breakage appeared to be insensitive to geogrid reinforcement and water content, but increased significantly with increasing confining pressures. The numerical results further indicate that the mechanical response of geogrids was regulated by axial strain and confining pressure. As axial strain was greater than 10%, increased confining pressure resulted in a wider area of stress concentration within the geogrid meshes, and the upper geogrid sustained a greater contact force and tensile strain. The findings in this study are of practical significance for the design and construction of structures in calcareous gravelly sand reinforced with geogrids.

A three-dimensional numerical model for computing the actual voltage induced in a Barkhausen surface coil

This investigation presents a novel three-dimensional simulation model designed for the analysis of Magnetic Barkhausen Noise (MBN). The model’s formulation is directly deduced from Maxwell’s equations, establishing a Boundary Value Problem (BVP) with boundary conditions tailored to a specific configuration involving a lengthy steel rod encased within a solenoid, subjected to sinusoidal current excitation, all contained within an air enclosure. This chosen configuration holds physical equivalence to the majority of setups commonly utilized in Barkhausen experiments. Utilizing the proposed three-dimensional model, it became feasible, for the first time, to computationally assess the voltage induced in the coil by the flux emanating from all spatial directions. The findings indicate that the MBN signal is primarily generated by the flux rate aligned with the direction of the applied field. Nevertheless, simulations confirm additional contributions to the signal originating from magnetic flux in directions significantly deviating from the applied field direction. It is demonstrated that these contributions from lateral flux significantly influence the angular dependence of the Barkhausen signal energy, showcasing a notable alignment between the proposed model and experimental results. The results underscore the capability of the proposed model to facilitate the simulation of Barkhausen Noise.

Local entropy generation optimization and thermodynamic irreversibility analysis of oval-shaped coaxial ground heat exchangers: A detailed numerical investigation

Optimizing the heat transfer and thermodynamic efficiency of Ground Heat Exchangers (GHE) is crucial for maximizing the energy efficiency of Ground-coupled Heat Pump (GCHP) systems utilized in building heating and cooling applications. This study pioneers an effective approach of both the first and second laws of thermodynamics of the newly oval-shaped coaxial GHEs (oval-CGHEs), assessing heat transfer, thermohydraulic performance, and local entropy generation rate to explore the optimum design and operational conditions for heat exchange maximization and entropy generation minimization. Utilizing a three-dimensional numerical model validated by experimental results, this study examines the trade-offs between heat transfer improvements and entropy generation rate intensification of the oval-CGHEs in comparison to the typical circular-shaped CGHE (circle-CGHE) further than exploring the effect of key operational and design parameters and performing sensitivity analysis. The results reveal that oval-CGHEs exhibit higher irreversibility as a penalty for their superior heat transfer capabilities compared to traditional circular-shaped CGHEs (circle-CGHE), which can be mitigated by optimizing the position of the inner tube. Moreover, the study underscores the criticality of maintaining a balance between enhanced heat transfer and the associated rise in irreversibility, particularly concerning the factors of mass flow rates and installation depths.

A simplified mathematical modeling approach for thermal runaway of Li(Ni0.8Co0.1Mn0.1)O2 pouch cells based on thermal runaway experimental data

Currently, the numerical modeling methods used to characterize the thermal runaway and propagation mechanisms in NCM811 pouch batteries still face numerous challenges. This is primarily due to the uncertainties in thermal runaway reactions and the immense computational demands encountered when modeling begins at the granular level of the cell. To address this issue, this study designed a series of experiments for an NCM811 pouch battery. By employing experimental data, reasonable assumptions were made regarding the combined effects of solid-state thermal conduction, electrochemical heating, and fluid-solid heat transfer through gas jet streams. A simplified thermal runaway mathematical modeling method based on experimental data was proposed, effectively tackling the problem. This framework was then applied to thermal runaway experiments on a single NCM811 battery used in this study, providing essential input parameters for three-dimensional numerical modeling of modules composed of multiple batteries in multiphysics simulation software. The simulation results derived from the proposed method were in good agreement with experimental data. This study can characterize the thermal runaway and thermal propagation characteristics of battery modules using single-cell thermal runaway experiments. Moreover, the proposed modeling approach not only predicted the experimental outcomes of the NCM811 pouch battery with reasonable accuracy but also provided smoother data, elucidating the highly fluctuating trends observed in experimental data. This method also offers guidance for designing fire prevention and smoke exhaust systems for battery packs. This study presents a novel approach for the mathematical modeling of thermal runaway and thermal propagation in battery modules using NCM811 pouch batteries, contributing to the ongoing development of lithium-ion battery module design.

Numerical simulation of jet boiling characteristics of high-pressure water injected into high-temperature liquid lead–bismuth in a confined space

The steam generator heat transfer tube rupture (SGTR) accident is potentially devastating and is one of the most serious design basis accidents for lead–bismuth cooled fast reactor (LFR). It is important to study the mechanism of water jet injecting into liquid lead–bismuth eutectic (LBE) for the safety analysis of SGTR accidents. In this paper, based on the volume of fluid (VOF) model, the large eddy simulation (LES) turbulence model and the Lee phase transition model, a three-dimensional numerical model is developed to investigate the typical characteristics of a high-pressure sub-cooled water jet injecting into high-temperature LBE in a confined space. The study focuses on examining the impact of water temperature, inlet pressure, and nozzle diameter on jet development and its surroundings, including the cover gas layer and LBE phase. The results indicate that the central region of the jet is primarily subject to flash boiling, particularly at the nozzle position, whereas the jet interface undergoes heat-transfer boiling. The residual liquid phase migrates upward along the interface and gradually boils. Based on the phase distribution characteristics, a typical jet can be categorized into four zones: the water-steam transition zone, the multiphase flow zone, the end water phase zone, and the steam block zone. Under the conditions analyzed in this paper (water temperature 140–260 °C, injection pressure 2–10 MPa, nozzle diameter 2–10 mm, LBE temperature 400 °C), the jet volume expands approximately exponentially, correlating positively with the three factors studied. It is also found that at high injection pressures and large nozzle diameters, the steam at the top of the jet is prone to destabilization and fragmentation, with steam volume and pressure oscillations being the dominant factors. Steam migration triggers pressure oscillations in the LBE region, and the intensity of these oscillations is closely linked to the migration rate. Under the conditions studied, the maximum pressure oscillation reached 5.5 times the initial pressure. The expansion of the jet steam drives the pressurization of the cover gas layer, which correlates positively with all three factors studied. The nozzle diameter has the greatest effect, increasing the maximum pressure by 168 % from the initial pressure.

Coupled optimization of auxiliary channels and porosity gradient of GDL for PEMFC

In this study, we proposed a novel under-rib auxiliary channel and developed a systematic methodology combining computational fluid dynamics (CFD), an artificial neural network (ANN), and a particle swarm optimization (PSO) algorithm for the coupled optimization of the auxiliary-channel structure and gas diffusion layer (GDL) porosity gradient in a proton exchange membrane fuel cell (PEMFC). First, we developed a three-dimensional multiphase numerical model of a non-isothermal parallel flow-field PEMFC. ANN was then applied to obtain a mapping relationship between auxiliary channels structure, GDL porosity gradient and output performance of PEMFC with finite CFD data. Finally, PSO was employed to perform an optimization search to determine the optimal parameter pairing, resulting in a 16.7% increase in maximum power density and a 20.6% increase in the current operating range. In addition, the optimized PEMFC exhibited superior mass transfer and water removal capabilities, where the average oxygen concentration in the intermediate section of cathode GDL was increased by 33.3% and the maximum liquid water saturation was reduced by 23.5%. Moreover, the coupled optimization of the auxiliary channels and the GDL porosity gradient can increase power density by at least 8.7% and current density by 9.1% compared to single component optimization.

Microparticle focusing and micromixing with two-dimensional acoustic waves

Acoustofluidics has recently been popularized as a crucial element of lab-on-a-chip (LoC) platforms to efficiently manipulate microparticles and continuous matter alike. In this study, a three-dimensional (3D) numerical model is proposed to simulate the focusing of polystyrene microparticles with three diameters and micromixing of dilute species using two orthogonally oriented standing waves, contrasting them with one-dimensional (1D) waves. The limiting velocity method is modified to explore the 3D acoustic streaming in a symmetric microchannel. In contrast to 1D standing acoustic waves, the simultaneous excitation of two orthogonal waves generates an acoustic streaming velocity field that does not counteract the radiation force. The obtained results show that the focusing efficiency of 5-μm particles reaches 97% with two dimensional (2D) standing acoustic waves, which was unachievable using 1D waves. Moreover, by reducing the flow rate to 1 μL min−1, the focusing of critical microparticle diameter peaked at 94%, indicating an approximately 9% improvement over a flow rate of 2.5 μL min−1. Increasing the viscosity of the background fluid resulted in 16% better 2D focusing with a single vortex compared to other cases, and higher amplitudes did not change focusing efficiency with a single vortex, while reducing efficiency in other cases. Finally, using 2D acoustic waves remarkably improved the mixing efficiency of dilute species, underscoring the advantage of 2D acoustic waves over their 1D counterpart. The proposed numerical model can play a meaningful role in cutting fabrication costs of next-generation LoC devices by identifying the most crucial parameters influencing acoustofluidic matter transport.

Numerical simulation of local scour around an inclined monopile in steady current

Local scour around the underwater supports of constructions would reduce the embedment depth and bearing capacity of support structures, seriously threatening the safety of the structures. Most existing investigations focused on the local scour of a vertical pile. However, in practical engineering, inclined piles are widely applied in underwater foundation structures to withstand lateral loads. The pile inclination can result in flow field alteration around the pile, thereby affecting the development of local scour. Consequently, a threedimensional numerical model of local scour around an inclined pile under steady current is established using the open-source CFD package OpenFOAM, in which the hydrodynamics are solved based on the Reynolds-averaged Navier-Stokes (RANS) equations together with k – ω SST turbulence closure in this model. It is then coupled with both bed load and suspended load transport descriptions, which drive the resultant morphology of the bed. The numerical model is verified by comparing it with the existing laboratory experiment. Discussions are carried out on the effects of inclination angles on the flow field alteration and the local scour depth variation.

A novel design of inlet/outlet header structure for improving flow uniformity and performance in planar solid oxide fuel cell

The poor flow uniformity of parallel flow field in Solid oxide fuel cell limits its performance. Based on computational fluid dynamics, a three-dimensional numerical model of a solid oxide fuel cell was established to improve performance by improving flow uniformity. The influence of different number, shape and arrangement of obstacles on fuel cell is studied, and optimization is also made according to flow distribution. Results show that flow uniformity can be improved by using rectangular obstacles and arranging them at the rear of inlet header. Compared with structure without obstacles, the structure with four rectangular obstacles in inlet header can reduce flow non-uniformity coefficient by 55.8%. Optimized structure uses obstacles with different heights, and two obstacles are also set in outlet header. Compared with basic structure, the flow non-uniformity coefficient of optimized structure is reduced by 87%, while current density and net power density are increased by 10.1% and 24.2% respectively.