Search

This article deals with the numerical simulation of a two-dimensional instantaneous heat and mass transfer processes within a commonly used intermetallic compound (a.k.a. Mischmetal) packed in a unit disc of an annulus-disc reactors, during hydrogen gas desorption. Using the finite volumes technique bundled in the OpenFOAM® CFD code, temperature and amount of desorbed hydrogen and their time-averaged quantities inside the metal-hydride packed bed are obtained for various temperatures of fluids used in heat transfer, and several outlet pressure magnitudes. Using a set of numerical simulations, we have emphasized the impacts of both parameters on metal-hydride reactor performance related to discharging time. An excellent accord was recorded for the present simulations results compared against the literature-reported experimental data [20].

Abrasive flow jet machining (ASJM) is modern manufacturing technique which uses comparatively low-pressure abrasive flow jet for machining various machined surfaces like holes, channels and intricate shapes which is not possible from conventional machining processes. The effect of abrasive particle concentration on the impacting slurry velocity on the surface was examined in the present CFD simulation, resulting in target surface erosion. To predict the impact velocity of abrasive particles hitting the surface, a 2D CFD model was used. The depth of machined surface, surface irregularity and surface removal rate at normal incidence mainly depend on the KE of particle, impact angle, etc. in line with previously published research work on highly pressurised air and water driven abrasive jet. In the present work CFD simulation is performed to predict the effect of solid particle concentration on the impacting slurry jet velocity at the target surface and also turbulence kinetic energy near the surface is studied. As per the CFD results the simulation model predictions the velocity of impacting particle goes on decreasing due to internal frictional resistance between solid and liquid phase and it shows scope of further parametric analysis in this area.

Particle Impact Drilling (PID) is a new technology to effectively improve the rate of penetration (ROP) for oil and gas drilling in hard and strongly abrasive formations. In this paper, numerical simulation method is used to analyse the motion characteristics and the modulation method of particle swarm in high-pressure tank for the particle injection system based on differential pressure ejection in PID. The numerical simulation results show that: when there is no modulation elements, the motion of particle swarm in the high-pressure tank follows an asymmetric funnel flow with pulsating state, which could be divided into vertical flow domain, fast flow domain, slow flow domain and stagnation domain. The unstable dynamic arching effect of the funnel flow, the viscous effect of the liquid bridge force and the collapsing effect of the particle swarm could probably lead to the blockage of the discharge port of the high-pressure tank. When the semiapex angles of the high-pressure tank decreases, the volume flow rate of particles increases and the stagnation domain becomes smaller, but it becomes easier to form arching and blockage. The modelling results indicate that the pulsation of the funnel flow is minimum when the semiapex angle is 45° without the mutilation element, which means the funnel flow of the particle swarm is relatively stable. By introducing a conical modulating element above the discharge port, the unstable funnel flow of the particle swarm could be transformed to an overall uniform flow. The modelling results indicate that the installation height of the modulation element has the greatest influence on the pulsation degree. The optimized parameters for the conical modulation element based on numerical modelling tests are 70° for the vertex angle, 35 mm for the length of the flank and 70 mm for the installation heigh.

Based on thermoelectric analogy method, a thermal network model coupled with electromagnetic loss is presented. According to the actual structure of dry-type transformer windings, an accurate thermal network is established. Considering the distribution characteristics of heat source along cooling ducts inside the dry-type transformer, the thermal parameters such as heat flux and convective heat transfer coefficients are calculated by finite element method analysis (FEM) and computational fluid dynamics (CFD). Based on temperature rise experiment, the hottest spot located at the top of outmost high voltage winding is 151.8 °C, which verified the correctness of thermal network model.

The origami principle of folding planar sheets into functional three-dimensional devices promises a future with increased compactness and reconfigurability of biomedical robots. To inspire new origami-based biomedical robots, we highlight two categories of origami patterns that are applicable to build compact actuators, deployable stents, and minimally invasive surgery devices. Due to the requirements in system integration, biocompatibility, and payload capability, the conventional paper-based origamis are converted to physical robots by replacing the zero-thickness facets and zero-width folds with thick panels and flexible hinges, respectively. Therefore, the physical origami-based robots become inhomogeneous, making it more complex and challenging to analyze their mechanics. So far, no such model can analyze the mechanics of any physical origami robot. Herein, we propose a computational model that discretizes the continuous structures into truss and spring elements to fulfill this goal. Based on the model, we simulate the identified origami structures under various boundary conditions to investigate their application in biomedical scenarios. This article is expected to accelerate the design iteration of new functional biomedical origami robots.

Liquid–solid fluidized beds (LSFB) modeling validation is crucial for establishing design rules and monitoring tools. However, it generally relies on comparing global variables, which overlook dynamic features that influence reaction outputs. This work aims to implement time series analysis tools to compare Radioactive Particle Tracking data with a simulation consisting of Computational Fluid Dynamics coupled with Discrete-Element Method. Experiments have been performed in a pilot-scale LSFB of calcium alginate spheres fluidized with a calcium chloride solution. The Diks’ test indicates that the simulation can capture the LSFB behavior. It also allows diagnosing flow regime transitions from the simulation. Trends of solid dispersion coefficients and mixing times predicted by the simulation are in good agreement with the experiments.

We are developing a universal, advanced ventricular assist device (AVAD) with automatic pressure regulation suitable for both left and right ventricular support. The primary goal of this computational fluid dynamics (CFD) study was to analyze the biventricular performance of the AVAD across its wide range of operating conditions. An AVAD CFD model was created and validated using in vitro hydraulic performance measurements taken over conditions spanning both left ventricular assist device (LVAD) and right ventricular assist device (RVAD) operation. Static pressure taps, placed throughout the pump, were used to validate the CFD results. The CFD model was then used to assess the change in hydraulic performance with varying rotor axial positions and identify potential design improvements. The hydraulic performance was simulated and measured at rotor speeds from 2,300 to 3,600 revolutions/min and flow rates from 2.0 to 8.0 L/min. The CFD-predicted hydraulic pressure rise agreed well with the in vitro measured data, within 6.5% at 2300 rpm and within 3.5% for the higher rotor speeds. The CFD successfully predicted wall static pressures, matching experimental values within 7%. High degree of similarity and circumferential uniformity in the pump's flow fields were observed over the pump operation as an LVAD and an RVAD. A secondary impeller axial clearance reduction resulted in a 10% decrease in peak flow residence time and lower static pressures on the secondary impeller. These lower static pressures suggest a reduction in the upwards rotor forces from the secondary impeller and a desired increase in the pressure sensitivity of the pump. The CFD analyses supported the feasibility of the proposed AVAD's use as an LVAD or an RVAD, over a wide range of operating conditions. The CFD results demonstrated the operability of the pump in providing the desired circumferential flow similarity over the intended range of flow/speed conditions and the intended functionality of the AVAD's automated pressure regulation.

Due to current and future environmental and safety issues in space propulsion, typical propellants for upper stage or satellite rocket engines such as the toxic hydrazine are going to be replaced by green propellants like the combination of liquid oxygen and hydrogen or methane. The injection of that kind of cryogenic fluids into the vacuum atmosphere of space leads to a superheated state, which results in a sudden and eruptive atomization due to flash boiling. For a detailed experimental investigation of superheated cryogenic fluids, the new cryogenic test bench M3.3 with a temperature controlled injection system was built at DLR Lampoldshausen. After a first test campaign with high-speed shadowgraphy of flash boiling liquid nitrogen sprays, a laser-based Phase Doppler system was set-up to determine the spatial distributions of droplet velocities and diameters in highly superheated sprays. The spatial distributions revealed a core region with high mean velocities close to the injector orifice. With increasing distance from the injector orifice, the sprays develop a more and more monodisperse pattern. These distributions also showed that atomization due to flash boiling generates finer sprays with growing degrees of superheat. In certain spray regions, two droplet populations varying in their direction of motion, velocity and diameter due to possible recirculation zones were observed. The experimental data of flash boiling liquid nitrogen generated within this study provide a comprehensive data base for the validation of numerical models and further numerical investigations., Deutsche Forschungsgemeinschaft, Projekt DEAL

Enhancing damping characteristic is one of the effective methods to solve the instability problem of the rotor system. The three-dimensional numerical analysis model of scallop damper seal was established, and the effects of inlet pressures, preswirl ratios, rotational speeds, interlaced angles and seal cavity depths on the rotordynamic characteristics of scallop damper seal were studied based on dynamic mesh method and multi-frequencies elliptic whirling model. Results show that the direct stiffness of the scallop damper seal increases with decreasing inlet pressure and increasing rotational speed and cavity depth. When the seal cavity is interlaced by a certain angle, which shows positive direct stiffness. The effective damping of the scallop damper seal increases with the increasing inlet pressure, the decreasing preswirl ratio and the rotational speed and cavity depth. There exists an optimal interlaced angle to maximize the effective damping and the system stability. The leakage of the scallop damper seal is significantly reduced with decreasing inlet pressure. The preswirl will reduce the leakage flowrate, and the rotational speed has a slight effect on the leakage performance. The leakage of the scallop damper seal decreases with increasing seal cavity depth.

The three-dimensional computational fluid dynamics (3D-CFD) of a pulsating flow applied to the fluid catalytic cracking (FCC) reaction was investigated in the riser of a circulating fluidized bed reactor. The kinetic parameters of the FCC and coke burning reactions for predicting the reactant conversion and product yield percentages were applied. To increase the reactant conversion level and product yield, the effect of the pulsating flow operating parameters was considered using a 2k statistical experimental design with four factors (amplitude, frequency, types of the waveform, and amplitude ratio). The 3D-CFD simulation was successfully validated from the experimental literature data. The frequency and type of the waveform were found to be the significant operating parameters. The expression of the fitted regression model and response surface contour were derived and revealed that the pulsating flow provides a higher reactant conversion level and product yield percentages compared to a non-pulsating or steady flow.

In addition to photovoltaic cells, solar power plates, photovoltaic lights, and solar pumping water, solar energy is the primary source of heat from the sun. At the moment, researchers are looking at the use of nanotechnological and solar radiation to increase aeronautical efficiency. In this study, hybrid nanofluid flow linearly passes through a parabolic trough solor collector (PTSC) on the interior of solar aircraft wings to study heat transmission. Solar radiative flow was the term used to describe the heat source. The heat transfer efficiency of the wings is evaluated for several effects, such as a slanted magnetic field, viscous dissipation, play heating, and thermal radiative flow. Entropy generation study was performed on the Prandtl-Eyring hybrid nanofluid (P-EHNF). The Keller box technique was used to solve the predicted energy and momentum equations. As a typical fluid, EG (ethylene glycol) is used to disperse the nanosolid particles, which consist of copper (Cu) and cobalt ferrite (CoFe2O4). A variety of control factors, including velocity, shear stress, and temperature outlines as well as a frictional factor and Nusselt number, are addressed in detail. Thermal radiation amplification and variable thermal conduction parameters appear to increase the efficiency of aircraft wings in terms of thermal transfer. Hybrid nanofluid is superior to conventional nanofluid in terms of heat transmission. Cu-EG has a low thermal efficiency between 3.8% and 4.8% than CoFe2O4-Cu/EG nanofluid.

This work uses a computational fluid dynamics approach to evaluate the ability of the aft element of a slotted, natural-laminar-flow airfoil, designed for transonic applications, to function as a h...

In this study, based on a closed bomb test combined with computational fluid dynamics, a structural finite element method, and an acoustic boundary element method, a fluid-solid acoustic one-way coupling calculation model is established for the combustion process of energetic materials in a closed bomb, and the effectiveness of the model is verified by experiments. It is found that the maximum peak sound pressure increases exponentially with an increase in loading doses or gas pressure. However, a change in the combustion coefficient of the energetic materials has little effect on the noise generated during the combustion process in the closed bomb. When the combustion coefficient is reduced by a multiple of 16, the maximum transient sound pressure is reduced by 1.79 dB, and the sound pressure level in the frequency band is reduced by 1.75 dB. With an increase in shell thickness, the combustion noise of the energetic materials in the closed bomb decreases, and the reduction range of the combustion noise increases with the increase in shell thickness.

Thermoacoustics is about conversion between thermal and acoustical energies to provide alternative green technology for power cycle and cooling system. The oscillatory flow across porous structure inside the system is playing the role of energy conversion between the acoustic wave and the surface of the porous structure. Better understanding of fluid dynamics of oscillatory flow inside thermoacoustic system is therefore important for the thermoacoustic based energy conversion system. This paper presents the ‘entrance’ and ‘exit’ effects of the oscillatory flow within a parallel-plate structure that is placed inside a standing-wave thermoacoustic environment. Two-dimensional SST k-ω CFD models which were validated using experimental data and theoretical predictions were used for this investigation. Two different operating frequencies of 14.2 Hz and 23.6 Hz were studied for flow with five different amplitudes that were represented using drive ratios of 0.3%, 0.83%, 1.5%, 2.1% and 3%. These correspond to cases with Reynolds numbers between 5936 and 62926. Due to the cyclic nature of the flow, a region defined as an ‘exit’ region was observed in addition to the usual ‘entrance’ region and the fully developed region for flow inside a channel. The change of shape of velocity profiles from the ‘m’ shape profile, to the ‘slug-like’ profile and ‘parabolic-like’ profile was discussed in relation to the ‘entrance’ and ‘exit’ effects on flow inside the channel. The ‘entrance’ and ‘exit’ effects become bigger as drive ratio increases. The effect of ‘entrance’ and ‘exit’ are slightly reducing as frequency increases from 14.2 Hz to 23.6 Hz. This may be related to the shorter travel distance of fluid as the frequency increases. The results shown in this paper suggest that the flow within a 200 mm parallel-plate structure should be treated as developing flow especially for flow with low resonance frequency at drive ratio higher than 1%.

Mesoscale bubbles exist inherently in bubbling fluidized beds and hence should be considered in the constitutive modeling of the drag force. The energy minimization multiscale bubbling (EMMS/bubbling) drag model takes the effects of mesoscale structures (i.e., bubbles) into the modeling of drag coefficient and thus improves the coarse-grid simulation of bubbling and turbulent fluidized beds. However, its dependence on the bubble diameter correlation has not been thoroughly investigated. The hydrodynamic disparity between homogeneous and heterogeneous fluidization is accounted for by the heterogeneity index, Hd, which can be affected by choice of bubble diameter correlation. How this choice of bubble diameter correlation influences the model prediction calls for further fundamental research. This article incorporated seven different bubble diameter correlations into EMMS/bubbling drag model and studied their effects on Hd. The performance of these correlations has been compared with the correlation used previously by EMMS/bubbling drag model. We found that some of the correlations predicted lower Hd by order of a magnitude than the correlation used by the original EMMS/bubbling drag. Based on such analysis, we proposed a modification in the EMMS drag model for bubbling and turbulent fluidized beds. A computational fluid dynamics (CFD) simulation using two-fluid model with the modified EMMS/bubbling drag model was performed for two bubbling and one turbulent fluidized beds. Voidage distribution, time averaged solid concentration and axial solid concentration profiles were studied and compared with the previous version of the EMMS/bubbling drag model and experimental data. We found that the right choice of bubble diameter correlations can significantly improve the results for CFD simulations.

Modern gas turbine engines operate at high temperatures to improve thermal efficiency and power output. Increased rotor inlet temperatures increase the rate of heat transfer to the turbine blades, which requires sophisticated cooling schemes to keep the blade temperature at acceptable levels. This work is a numerical investigation of film cooling techniques as applied to gas turbines. The cooling performance of two differently shaped holes, namely, Ginkgo Forward and Ginkgo Reverse, were investigated in terms of centerline and local lateral cooling effectiveness, and a comprehensive comparison was made with the cooling performance of a cylindrical hole. The investigations were performed at a constant density ratio (DR = 2.0) and three different blowing ratios (BR = 1.0, 1.5, and 2.0). Under all of the operating conditions, the results demonstrated significant augmentation in centerline and lateral cooling effectiveness when the Ginkgo Reverse shaped hole was used, followed by the Ginkgo Forward and cylindrical cooling holes. For the shaped cooling holes, the low velocity gradient through the film alleviated the jet lift-off and turbulence, resulting in decreased entrainment of hot gas to the bottom surface. To conclude, the prominent lateral dispersal of the coolant due to the shaped cooling holes significantly enhanced thermal protection and the overall cooling performance.

Saturation flow is a fundamental performance measure parameter used for the design of intersections. However, in mixed traffic conditions, the evaluation of saturation flow is more sensitive due to heterogeneous traffic where vehicles are travelling in the same right of way without any lane separation. In such cases, the passenger car unit (PCU) plays a major role in analysing heterogeneity using an equivalent value for each vehicle. This paper presents a novel approach using the optimisation technique to determine PCU values and saturation flow for mixed traffic streams. Field data were collected at six signalised intersections in three Indian cities. The optimisation technique produced reasonable PCU values and an accurate estimation of saturation flow of 1906 PCUs/h. The saturation flow prediction model was developed using the multi-linear regression technique as per Pearson's correlation between saturation flow and vehicle compositional share. The developed saturation flow model was validated using field data from another site and the predicted value was in good agreement with the observed value. Finally, the predictive model was used to analyse the impact of traffic composition on saturation flow in mixed traffic streams.

The force-generation mechanism of a dovelike flapping-wing micro air vehicle was studied by numerical simulation and experiment. To obtain the real deformation pattern of the flapping wing, the digital image correlation technology was used to measure the dynamic deformation of the wing. The dynamic deformation data were subsequently interpolated and embedded into the CFD solver to account for the aeroelastic effects. The dynamic deformation data were further used to calculate the inertial forces by regarding the wing as a system of particles to take into account the wing flexibility. The temporal variation of the forces produced by the flapping wing was measured by a miniature load cell. The numerical results provide more flow details of the unsteady aerodynamics of the flapping wing in terms of vortex formation and evolution. The calculated results of the inertial forces are analyzed and compared with the CFD results which represent the aerodynamic forces. In addition, the total forces, i.e., the sum of the CFD result and inertial result, are compared with the experimental results, and an overall good agreement is obtained.

Cleveland Clinic is developing a continuous-flow total artificial heart (CFTAH). This novel design operates without valves and is suspended both axially and radially through the balancing of the magnetic and hydrodynamic forces. A series of long-term animal studies with no anticoagulation demonstrated good biocompatibility, without any thromboemboli or infarctions in the organs. However, we observed varying degrees of thrombus attached to the right impeller blades following device explant. No thrombus was found attached to the left impeller blades. The goals for this study were: (1) to use computational fluid dynamics (CFD) to gain insight into the differences in the flow fields surrounding both impellers, and (2) to leverage that knowledge in identifying an improved next-generation right impeller design that could reduce the potential for thrombus formation. Transient CFD simulations of the CFTAH at a blood flow rate and impeller rotational speed mimicking in vivo conditions revealed significant blade tip-induced flow separation and clustered regions of low wall shear stress near the right impeller that were not present for the left impeller. Numerous right impeller design variations were modeled, including changes to the impeller cone angle, number of blades, blade pattern, blade shape, and inlet housing design. The preferred, next-generation right impeller design incorporated a steeper cone angle, a primary/splitter blade design similar to the left impeller, and an increased blade curvature to better align the incoming flow with the impeller blade tips. The next-generation impeller design reduced both the extent of low shear regions near the right impeller surface and flow separation from the blade leading edges, while maintaining the desired hydraulic performance of the original CFTAH design.

The coaxial mixers enhance the suspension of concentrated slurries in an agitated reactor. In this research work, the complex slurry suspension and dissemination behavior in a coaxial slurry mixing system (comprised of a close clearance anchor rotating with a low speed and an inner axial impeller rotating with a high speed) was analyzed employing ERT (electrical resistance tomography, a non-intrusive flow visualization technique), and computational fluid dynamics (CFD). The numerical models were validated by comparing the axial solid concentration profiles generated using the ERT data and the CFD simulation results. The influences of various important parameters such as the diameter of the inner axial impeller, the inner impeller type, and the inner impeller spacing on the hydrodynamic characteristics of the slurry suspensions in a coaxial mixing vessel were thoroughly analyzed. The radial and axial velocity profiles of solid particles were generated using the validated mathematical models. The assessment of energy loss due to the solid–solid collisions, the particle–fluid frictions, and the particle–vessel wall collisions was conducted. The evaluation of optimum inner impeller clearance and inner impeller diameter is essential to attain a high degree of solids suspension and dissemination in a coaxial slurry mixing system.

The dynamic analysis of prestressed, bending-torsion coupled beams is revisited. The axially loaded beam is assumed to be slender, isotropic, homogeneous, and linearly elastic, exhibiting coupled flexural-torsional displacement caused by the end moment. Based on the Euler-Bernoulli bending and St. Venant torsion beam theories, the vibration and stability of such beams are explored. Using the closed-form solutions of the uncoupled portions of the governing equations as the basis functions of approximation space, the dynamic, frequency-dependent, interpolation functions are developed, which are then used in conjunction with the weighted residual method to develop the Dynamic Finite Element (DFE) of the system. Having implemented the DFE in a MATLAB-based code, the resulting nonlinear eigenvalue problem is then solved to determine the coupled natural frequencies of illustrative beam examples, subjected to various boundary and load conditions. The proposed method is validated against limited available experimental and analytical data, those obtained from an in-house conventional Finite Element Method (FEM) code and FEM-based commercial software (ANSYS). In comparison with FEM, the DFE exhibits higher convergence rates and in the absence of end moment it produces exact results. Buckling analysis is also carried out to determine the critical end moment and compressive force for various load combinations.

Solid-liquid mixing plays a significant role in crystallization, suspension polymerization, leaching, solid-catalyzed reaction and adsorption. In this study, a computational fluid dynamic (CFD) model was developed for solid-liquid mixing in a cylindrical tank equipped with a top-entering impeller. The multiple reference frame (MRF) technique, k-ε model and Eulerian-Eulerian approach were employed to simulate the impeller rotation, turbulent flow and multiphase flow, respectively. The effects of impeller speed, solid concentration, particle size, solid density and impeller clearance on the mixing performance of four different impellers (A310, marine propeller, pitched blad turbine and A320) were investigated. The CFD results were in good agreement with experimental data measured using electrical resistance tomography (ERT). In order to investigate the mixing quality in this study, the impeller speed required for maximum homogeneity, clouding height, and just-suspended impeller speed were investigated.

In this paper, we introduce techniques for animating explosions and their effects. The primary effect of an explosion is a disturbance that causes a shock wave to propagate through the surrounding medium. This disturbance determines the behavior of nearly all other secondary effects seen in explosions. We simulate the propagation of an explosion through the surrounding air using a computational fluid dynamics model based on the equations for compressible, viscous flow. To model the numerically stable formulation of shocks along blast wave fronts, we employ an integration method that can handle steep gradients without introducing inappropriate damping. The system includes two-way coupling between solid objects and surrounding fluid. Using this technique, we can generate a variety of effects including shaped explosive charges, a projectile propelled from a chamber by an explosion, and objects damaged by a blast. With appropriate rendering techniques, our explosion model can be used to create such visual effects such as fireballs, dust clouds, and the refraction of light caused by a blast wave., 7 pages, 12 figures, anc videos, 22 citations. Alternative online location: http://graphics.berkeley.edu/papers/Yngve-AEX-2000-07

Single-orifice plate is wildly utilized in the piping system of the nuclear power plant to throttle and depressurize the fluid of the pipeline. The cavitation induced by the single-orifice plate may cause some serious vibration of the pipeline. This study aims to find the optimal designs of the single-orifice plates that may have weak cavitation possibilities. For this purpose, a new single-orifice plate with a convergent-flat-divergent hole was modeled, a multi-objective optimization method was proposed to optimize the shape of a single-orifice plate, while computational fluid dynamics method was adopted to obtain the fluid physical quantities. The reciprocal cavitation number and the developmental integral were treated as cavitation indexes (e.g., objectives for the optimization algorithm). Two non-dominant designs ultimately achieved illustrated obvious reduction in the cavitation indexes at a Reynolds number Re = 1 × 105 defined based on fluid velocity. Besides, the sensitivity analysis and temperature effects were also performed. The results indicated that the convergent angle of the single-orifice plate dominants the cavitation behavior globally. The optimal designs of single-orifice plates result in lower downstream jet areas and lower upstream pressure. For a constant Reynolds number, the higher temperature of liquid water, the easier it is to undergo cavitation. Whereas there is a diametric phenomenon for a constant fluid velocity. Moreover, the regression models were carried out to establish the mathematical relation between temperature and cavitation indexes.

Multi-jet interactions are characterized by a number of shock and expansion waves, formation of turbulent shear layers and subsonic recirculation zones. The problem of the occurrence of unsteady, transient and oscillatory, processes arising in the nozzle devices of rocket engines at different flight modes is discussed. The main flow regimes, the sequence of their change and qualitative patterns of shock-wave structures that arise in each of the regimes are considered. The main emphasis is paid to the causes of occurrence and mechanisms of maintenance of low-frequency oscillations of the base pressure. The regularities of changes in the amplitude-frequency characteristics of oscillations are studied depending on the main parameters of the problem. Areas of existence of oscillatory processes are found, and a search is carried out for the ranges of design parameters, geometric characteristics of nozzle assemblies that provide a non oscillating flow of gas-dynamic processes. Hysteresis phenomena of shock-wave structure at increase or decrease in the reservoir total pressure is observed.

In this study, the effect of channel baffles and louver baffles on the flow pattern in the large-scale industrial fluidized beds was studied by computational fluid dynamics (CFD) methods. Then, the effect of flow pattern on the chemical reaction performance was studied for the first time. Simulation results showed that the gas velocity distributed more uniformly, solid particles dispersed more homogeneously and aggregation scarcely occurred in the fluidized bed with louver baffles than that with channel baffles. The residence time distribution indicated that louver baffles remarkably suppressed gas back-mixing in comparison with channel baffles. The reasonable agreements of pressure distribution and reaction results between the simulation in the bed with channel baffles and the data on a large-scale industrial apparatus demonstrated the accuracy of the CFD model. The predicted conversion of SiCl4 in the bed with louver baffles (27.44%) was higher than that with channel baffles (22.69%), indicating that louver baffles markedly improved the performance of the fluidized bed. This study could provide useful information for future structural improvements of baffles in large-scale fluidized beds.

A novel method for close-proximity formation flying under differential atmospheric drag using Coulomb forces is investigated for applications in Earth sensing, space-situational awareness (SSA), and aeronomy. Objects in LEO are supersonic with respect to the ambient environment, creating a thinned out wake region behind the craft as it travels through the ionosphere. Objects within this wake experience little drag acceleration and are able to attain voltages much greater than in the ambient ionospheric plasma, creating implications for the design and control of close-proximity leader–follower spacecraft pairs. The proposed system consists of a leader craft with a set of affixed, conducting spheres and a charged follower craft located in the wake of the leader. The differential drag acceleration between the leader and follower craft is countered by a controlled Coulomb repulsion to maintain precise separation. The charged structure on the rear of the leader craft is designed such that the charged follower craft sits in an electrostatic potential well which opposes off-axis perturbations. A conceptual method for controlling such a pair without the use of propellant using a set of charged spheres is investigated, with nonlinear models of the system’s relative motion derived and discussed. Linearized models are used to demonstrate the local controllability of the system to demonstrate the proposed system’s merit. This linear analysis is used to derive conditions on controllability and control performance under different charge geometries and environmental assumptions.

A wave rotor is suitable for compact and efficient pressure-exchange between gas flows. This work measured the circumferential pressure distribution of the rotor/stator interfaces and utilized a CFD method to simulate the unsteady pressure waves. The experimental and CFD results showed some slopes in the circumferential pressure distributions, and the slopes indicated the traces of specific unsteady pressure waves. Such traces varied regularly if the rotational speed varied within a range from −11% to +11% off the baseline value, but they were seriously disturbed if the rotational speed varied by −45% from the baseline value. It verified that a pressure wave in a wave rotor tended to keep its pressure ratio and propagation velocity unchanged if the rotational speed varied by a small extent, and that the pressure wave could not keep its propagation patterns if the rotational speed varied by a large extent. Because of the pressure wave behaviors, the wave rotor demonstrated specific regulations of the rotational speed effects on its operational states.

This paper studies a novel gas cyclone with a cylindrical filter face installed in the center from the vortex finder to the bottom hopper. The experimental results show that this composite cyclone has a higher collection efficiency and a lower pressure drop than the original cyclone. The mechanisms for the improvement are analyzed by both physical experiments and numerical simulations. By measuring dust samples collected at different places it is revealed that the center filter can prevent fine particles from entering the inner vortex and escaping, which accounts for the increase of the collection efficiency. In addition, the flow field of the composite cyclone is simulated by computational fluid dynamics and compared with that of the original cyclone. The analysis shows that with the filter layer installed, the swirling flow disappears in the vortex finder, which decreases the kinetic energy dissipation and hence lowers the pressure drop.

Radiative heat transfer in a concentric cylinder system with a participating medium is demonstrated to be inherently two-dimensional, even for the consideration of heat transfer to the total inner ...

The flow dynamics inside cavities associated with turbulent incoming flow and oscillations of the cavity shear layer is an important field of study that has attracted many researchers aiming to understand the complex flow features using both numerical and experimental techniques. Such a flow exists in many industries, such as in the automobile industry for flows over the sunroof of vehicles, and in the energy industry such as in finned heat-exchangers, combustion chambers, and turbines that utilizes wind energy harvesting in buildings. Most of those studies focused on the aero-acoustic mechanism and its control. Despite its wide industrial applications, heat transfer in cavity with oscillating shear layer remains less studied in the literature. The present paper presents a literature review of various studies that focuses on the heat transfer mechanism in those cavity configurations and how it is affected by the flow parameters.

Proppants transport is an advanced technique to improve the hydraulic fracture phenomenon, in order to promote the versatility of gas/oil reservoirs. A numerical simulation of proppants transport at both hydraulic fracture (HF) and natural fracture (NF) intersection is performed to provide a better understanding of key factors which cause, or contribute to proppants transport in HF–NF intersection. Computational fluid dynamics (CFD) in association with discrete element method (DEM) is used to model the complex interactions between proppant particles, host fluid medium and fractured walls. The effect of non-spherical geometry of particles is considered in this model, using the multi-sphere method. All interaction forces between fluid flow and particles are considered in the computational model. Moreover, the interactions of particle–particle and particle–wall are taken into account via Hertz–Mindlin model. The results of the CFD-DEM simulations are compared to the experimental data. It is found that the CFD-DEM simulation is capable of predicting proppant transport and deposition quality at intersections which are in agreement with experimental data. The results indicate that the HF–NF intersection type, fluid velocity and NF aperture affect the quality of blockage occurrence, presenting a new index, called the blockage coefficient which indicates the severity of the blockage.

This study was conducted to investigate the effect of absorption coefficient models on the P1 radiation model for a premixed hydrogen flame containing the water vapor. A CFD combustion simulation analysis was performed using XiFoam, one of the open-source CFD solvers in OpenFOAM. The solver using the flamelet combustion model has been modified to implement radiative heat transfer. The absorption coefficient models used in this study the grey-mean model and constant model, and for comparison, case without radiation was added. This CFD simulation study consisted of benchmarking the THAI HD-15 and HD-22 experiments. The difference between the two tests is the inclusion of water vapor in the condition before ignition. In the case of the HD-22 experiment containing water vapor in the initial condition, the simulation results show that the grey-mean absorption coefficient model has a strong influence on the temperature decrease of the flame and on the change in pressure inside the vessel.

The chemical looping hydrogen production (CLHP) process is a new approach for hydrogen ( H 2 ) fuel production. The proposed process consists of three reactors which are air reactor (AR), fuel reactor (FR) and steam reactor (SR). In this study, the gas solid flow behaviour in a cold model of the proposed CLHP process was simulated using computational fluid dynamics (CFD) with kinetic theory of granular flow. The effect of drag coefficient models on the pressure profile of each reactor was investigated and the results from the modified Syamlal–O’Brien drag model agreed well with the experimental data. The model was further used to investigate effects of the operating parameters on the hydrodynamics profiles of each reactor. The solid flux increased with the increasing of inlet velocity of AR and FR but decreased with the increasing of the inlet velocity of SR. At the same time, the solid flux increased with the increasing of the solid inventory.

An appropriate microenvironment for preserving cultural relics is essential, and the air-water direct contact technology is utilized to create the microenvironment recently. The influence of a deflector in a tank was numerically investigated based on uniform design method to improve the heat and mass transfer and pressure drop performance of the air-water direct tank. In this study, a simplified CFD-based model was established and validated between airstream and water surface within the tank, to analyze the heat and mass transfer and pressure drop processes. Meanwhile, regression models of the heat transfers rate, mass transfer rate and pressure drop were developed by uniform design method based on three parameters: installation position, tilt angle, and height of the deflector, in order to analyze the influences of these three parameters on the heat and mass transfer and pressure drop of the tank. Finally, all three optimal structural parameters of the deflector were obtained based on the proposed comprehensive evaluation index using a genetic algorithm. The results showed that the model established for air-water direct contact adopted well to predict the heat and mass transfer and pressure drop performance between airstream and still water surface. Furthermore, the results found that the flow field inside the water tank was affected by the deflector's structure, which affected the heat and mass transfer performance. The simulation results suggested that the deflector's optimal structural parameters are 8 mm of installation position, 88 ° of tilt angle and 19 mm of height, respectively, within a given extent in this study.

A 90° elbow equipped with guide vanes was developed with the intent of reducing elbow erosion. Numerical models were formed to predict the maximum erosion rate of elbow and Face-1, and the response surface methodology was used to study the relationship between the erosion rate and structural parameters of guide vane. A second-order response surface model was established to determine the relationship between R1 and variables, and a reduced cubic (RC) polynomial model was obtained to reveal the relationship between R2 and the three factors. The numerical results show that the low-speed region is expended and the maximum discrete particle matter (DPM) concentration is reduced after installing the guide vane. This internal component provides a shelter for the elbow from the direct impact of high-speed solids.

Cleveland Clinic’s continuous-flow total artificial heart (CFTAH) is a double-ended centrifugal blood pump that has a single rotating assembly with an embedded magnet, which is axially and radially suspended by a balance of magnetic and hydrodynamic forces. The key to the radial suspension is a radial offset between the stator-bearing bore and the magnet’s steel laminations. This offset applies a radial magnetic force, which is balanced by a hydrodynamic force as the rotating assembly moves to a “force-balanced” radial position. The journal-bearing blood passage is a narrow flow path between the left and right impellers. The intent of this study was to determine the impact of the stator-bearing bore radius on the journal-bearing hydraulic performance, while satisfying the geometric design constraints imposed by the pump and motor configuration. Electromagnetic forces on the journal-bearing were calculated using the ANSYS EMAG program, Version 18 (ANSYS, Canonsburg, PA). ANSYS CFX Version 19.2 was then used to model the journal-bearing flow paths of the most recent design of the CFTAH. A transient, moving mesh approach was used to locate the steady state, force-balanced position of the rotating assembly. The blood was modeled as a non-Newtonian fluid. The computational fluid dynamics simulations showed that, by increasing stator bore radius, rotor power, stator wall average shear stress, and blood residence time in journal-bearing decrease, while blood net flow rate through the bearing increases. The results were used to select a new bearing design that provides an improved performance compared with the baseline design. The performance of the new CFTAH-bearing design will be confirmed through upcoming in vitro and in vivo testing.

A theoretical treatment on the oxide-controlled dwell fatigue crack growth of a γ’ strengthened nickel-based superalloys is presented. In particular, this study investigates the influence of an externally applied load and variations in the γ’ dispersion on the grain boundary oxide growth kinetics. A dislocation-based viscoplastic constitutive description for high temperature deformation is used to simulate the stress state evolution in the vicinity of a crack at elevated temperature. The viscoplastic model explicitly accounts for multimodal γ’ particle size distributions. A multicomponent mass transport formulation is used to simulate the formation/evolution of an oxide wedge ahead of the crack tip, where stress-assisted vacancy diffusion is assumed to operate. The resulting set of constitutive and mass transport equations have been implemented within a finite element scheme. Comparison of predicted compositional fields across the matrix/oxide interface are compared with experiments and shown to be in good agreement. Simulations indicate that the presence of a fine γ’ size distribution has a strong influence on the predicted ow stress of the material and consequently on the relaxation in the vicinity of the crack-tip/oxide wedge. It is shown that a unimodal dispersion leads to reduced oxide growth rates (parabolic behavior) when compared to a bimodal one. Stability conditions for oxide formation are investigated and is associated with the prediction of compressive stresses within the oxide layer just ahead of the crack tip, which become progressively negative as the oxide wedge develops. However, mechanical equilibrium requirements induce tensile stresses at the tip of the oxide wedge, where failure of the oxide is predicted. The time taken to reach this critical stress for oxide failure has been calculated, from which dwell crack growth rates are computationally derived. The predicted rates are shown to be in good agreement with available experimental data.

In this study, an inverse-problem method was applied to estimate the solid concentration in a solid–liquid two-phase flow. An algebraic slip mixture model was introduced to solve the forward problem of solid–liquid convective heat transfer. The time-average conservation equations of mass, momentum, energy, as well as the volume fraction equation were computed in a computational fluid dynamics (CFD) simulation. The solid concentration in the CFD model was controlled using an external program that included the inversion iteration, and an optimal estimation was performed via experimental measurements. Experiments using a fly-ash–water mixture and sand–water mixture with different solid concentrations in a horizontal pipeline were conducted to verify the accuracy of the inverse-problem method. The estimated results were rectified using a method based on the relationship between the estimated results and estimation error; consequently, the accuracy of the corrected inversion results improved significantly. After a verification through experiments, the inverse-problem method was concluded to be feasible for predicting the solid concentration, as the estimation error of the corrected results was within 7% for all experimental samples for a solid concentration of less than 50%. The inverse-problem method is expected to provide accurate predictions of the solid concentration in solid–liquid two-phase flow systems.

Aiming at the main problems in the design and analysis of well wall strengthening for fractured formations, a three-dimensional (3D) numerical simulation scheme is proposed. The three-dimensional finite element software is used to analyze the mechanical behavior of fractures and the pressure sealing process, and evaluate the stress cage effect. The main features of the model are as follows: (1) The equivalent fractures in the analytical model represents the function sum of the mechanical behavior of all fractures on the well wall, which is a functionally equivalent crack.(2) When evaluating the stress cage effect, the shape of the crack wedge filled with the plugging agent particles is formed by simulating the fractures opening process under the injection pressure, not a given regular shape.(3) In the model of calculating bull heading of block agent, the liquid pressure on the well wall is the injection pressure, which is a variation increased with time. The fluid pressure on the well wall in the stress cage calculation model is generated by the initial pore pressure. (4) The numerical evaluation of the stress cage effect is achieved by calculating the increase amplitude FX of the minimum hoop stress on the well wall. Using this model, several sets of injection pressure design values can be used for pressure plugging numerical simulation and stress cage effect evaluation calculation, and then the optimal and accurate quantitative value of “injection pressure, block agent particle size, safe mud window upper bound” are found through comparison. Finally, through an engineering example of horizontal well drilling pressure sealing in a shale gas reservoir developed by a fracture, we use the above theoretical tools to introduce the process and results of the numerical analysis of the extended mud window drilled by the shale gas fissure reservoir.

s Feed-water valves (FWVs) are used to regulate the flow rate of water entering steam generators, which are very important devices in nuclear power plants. Due to the working environment of relatively high pressure and temperature, there is strength failure problem of valve body in some cases. Based on the thermo-fluid-solid coupling model, the valve body stress of the feed-water valve in the opening process is investigated. The flow field characteristics inside the valve and temperature change of the valve body with time are studied. The stress analysis of the valve body is carried out considering mechanical stress and thermal stress comprehensively. The results show that the area with relatively high-velocity area moves gradually from the bottom of the cross section to the top of the cross section with the increase of the opening degree. The whole valve body reaches the same temperature of 250 °C at the time of 1894 s. The maximum stress of the valve body meets the design requirements by stress assessment. This work can be referred for the design of FWVs and other similar valves.

Low fuel regression rate is one of the major drawbacks of hybrid rocket motors. In this study, a novel injector concept is proposed to provide a substantial enhancement in the fuel regression rate....

Ground effect dramatically improves the performance of the Wing-in-Ground (WIG) vehicles near the ground. However, making coordinate turnings in ground effect zone may raise complexity in flight control. The Angle of Attack (AOA) and height are believed to be important factors. To find the impact of these factors, three straight rectangular wings with different aspect ratios are simulated via CFD approach. The results show that in normal situations, the rolling moment tends to level the wing but less effective with the decreasing AOA and changes its direction when the AOA is less than a height dependent value. The yaw moment exhibits complex behavior because two different effects generate oppose moments. The rolling moment coefficient is more influenced by the ground if a wing has larger Aspect Ratio (AR) and works at low to medium AOA. The yawing coefficient is always more influenced if the wing has larger AR regardless AOAs. These findings can provide basic guidance for the design of control system and enable aircraft designers to reshape the crafts to avoid or utilize this effect.

In this research, a Graphical Processing Unit (GPU) accelerated Discrete Element Method (DEM) code was developed and coupled with the Computational Fluid Dynamic (CFD) software MFiX to simulate granular and multiphase flows with heat transfers and chemical reactions. The Fortran-based CFD solver was coupled with the CUDA/C++ based DEM solver through inter-process pipes. The speedup to the CPU version of MFiX-DEM is about 130–243 folds in the simulation of particle packings. In fluidized bed simulations, the DEM computation time is reduced from 91% to 17% with a speedup of 78 folds. The simulation of Geldart A particle fluidization revealed a similar level of importance of both fluid and particle coarse-graining. The filtered drag derived from the two-fluid model is suitable for Euler-Lagrangian simulations with both fluid and particle coarse-graining. It overcorrects the influence of sub-grid structures if used for simulations with only fluid coarse-graining.

A diaphragm can be installed inside a hybrid motor to enhance the fuel regression rate and combustion efficiency in hybrid rockets. However, significant pressure oscillations have been observed in ...

The micro-nozzle allows small spacecraft to realize the orbital maneuvering and precise station-keeping. Due to the rarefaction effect of the gas flow in micro-nozzles, the velocity-slip and temperature-jump exists near inner wall, and thus, affecting the operating characteristics of micro-nozzles. Various gas-solid theories and models have been developed to describe the velocity slip and temperature jump in computational fluid dynamics (CFD). In this study, numerical simulations of a micro-nozzle with five classical slip boundary models, which are based on the Maxwell scattering and Langmuir adsorption, are conducted to analyze the applicability of these models through comparing with the results obtained by direct simulated Monte Carlo (DSMC) method. The results show that the inlet pressure dramatically influences the flow regime of the micro-nozzle. The lower the inlet pressure, the higher the influence of the rarefaction effect. The results also indicate that the Langmuir slip model can accurately simulate the flow field in the transition regime. However, the modified Maxwell and Langmuir models proposed by Agrawal et al. (2008) and Le et al. (2012), respectively, are valid in the slip regime.

Liquid two-phase flow has been studied widely by many authors. Two phase flow exist in petroleum industries, nuclear power generation and chemical plants. In the present study a review on two phase flow such as oil/water was discussed, also the patterns of the flow and types of measurement are presented. Also a review on the new trends using the computational fluid dynamic (CFD) to predict the liquid-liquid two phase are viewed.

Hydraulic transient is a flow situation where the velocity of the flow and the pressure change rapidly with time in pipes networks filled with fluid. Water hammer is one of the consequences that may occur accidentally in the reactors cooling system. The pressure level can exceed the limits given by the manufacturer, and it leads to the failure of the pipes integrity or nuclear fuel elements fixation in the core. Cooling system involves many components, and when a flow control component changes abruptly its status (for example, a valve closing or pump trip), it causes rapid changes in velocity and propagates pressure pulses in short time. The cooling system is equipped with flapper valve and flywheel which overcomes the transient issues. A pump or check valve may be stopped or closed abruptly and may cause bad impacts on the system which is simulated with soft trip using flywheel and without flywheel in case of its failure for any accidental reasons. In addition the abrupt valve closure is deemed as a major issue in the system if occurred. AFT Impulse 5.0 software is used to simulate and study the transient of each case. The simulation appears that the reactor cooling system may suffer from surge that may occur during the sudden closure of check valve but in other hand it can be avoided by regular maintenance of flapper valve to open in timely manner.