A quasi 1D model of two-phase flow for a urea-selective catalytic reduction (SCR) system is developed which can calculate not only the generation of reducing agent but also the formation of deposits in the exhaust pipe. The gas phase flow is solved through Euler method, variables are stored on staggered grids, and the semi-implicit method for pressure-linked equation (SIMPLE) algorithm is applied to decouple the pressure and velocity. The liquid phase is treated in a Lagrangian way, which solves the equations of droplet motion, evaporation, thermolysis, and spray wall interaction. A combination of a direct decomposition model and a kinetic model is implemented to describe the different decomposition behaviors of urea in the droplet phase and wall film, respectively. A new 1D wall film model is proposed, and the equations of wall film motion, evaporation, thermolysis, and species transport are solved. The position, weight, and components of deposits can be simulated following implementation of the semi-detailed kinetic model. The simulation results show that a decrease in the exhaust temperature will increase the wall film region and the weight of deposits. Deposit components are highly dependent on temperature. The urea-water-solution (UWS) injection rate can affect the total mass of wall film and expand the film region, but it has little influence on deposit components. An increase in exhaust mass flow can decrease the total weight of deposits on the pipe wall because of the promotion of the mass and heat transfer process both in the droplets and wall film.

In this paper, a dynamic model is established for a two-stage rotor system connected by a gear coupling and supported on ball bearings with squeeze film dampers (SFDs). The nonlinear dynamic behavior of the rotor system is studied under misalignment fault condition. The meshing force of the gear coupling is calculated considering the deformation of the tooth caused by torque transmission and dynamic vibration. The contact force between the ball and race is computed based on the Hertzian elastic contact deformation theory and the elastohydrodynamic lubrication theory. The supported force of SFD is simulated by integrating the pressure distribution derived from Reynolds’s equation. The equations of motion are rewritten in non-dimensional differential form, and the fourth-order Runge–Kutta method is employed to solve the nonlinear dynamic equilibrium equations iteratively. To verify the validity of the dynamic model and the correctness of the numerical solution method, the experimental power spectra of the rotor system under various misalignment degrees are compared with the analytical results. The effects of several important parameters, such as the lubrication of the ball bearing, the centralizing spring stiffness, the radial clearance of SFD, and the misalignment of gear coupling, on the dynamic characteristics of the rotor system are investigated and discussed mainly focusing on the system stability. The response spectra, bifurcation diagrams, and Pointcaré maps are analyzed accordingly. These parametric analyses are very helpful in the development of a high-speed rotor system and provide a theoretical reference for the vibration control and optimal design of rotating machinery.

An excellent airfoil with a high lift-to-drag ratio may decrease oil consumption and enhance the voyage. Based on NACA 0012, an improved airfoil is explored in this paper. The class/shape function transformation has been proved to be a good method for airfoil parameterization, and in this paper it is modified to improve imitation accuracy. The computational fluid dynamics method is applied to obtain numerically the aerodynamic parameters of the parameterized airfoil, and the result is proved credible by comparison with available experimental data in the open literature. A polynomial-based response surface model and the uniform Latin hypercube sampling method are employed to decrease computational cost. Finally, the nonlinear programming by quadratic Lagrangian method is utilized to modify the multi-island genetic algorithm, which has an improved optimization effect than the method used on its own. The obtained result shows that the modified class/shape function transformation method produces a better imitation of an airfoil in the nose and tail regions than the original method, and that it will satisfy the tolerance zone of the model in a wind tunnel. The response surface model based on the uniform Latin hypercube sampling method gives an accurate prediction of the lift-to-drag ratio with changes in the design variables. The numerical result of the flow around the airfoil shows reasonable agreement with the experimental data graphically and quantitatively. Ultimately, an airfoil with better capacity than the original one is acquired using the multi-island genetic algorithm based nonlinear programming by quadratic Lagrangian optimization method. The pressure contours and lift-to-drag ratio along with the attack angle have been compared with those of the original airfoil, and the results demonstrate the strength of the optimized airfoil. The process for exploring an improved airfoil through parameterization to optimization is worth referencing in future work.

A modified 3D finite element (3D-FE) model is developed under the FE software environment of LS-DYNA based on characteristics of stagger spinning process and actual production conditions. Several important characteristics of the model are proposed, including full model, hexahedral element, speed boundary mode, full simulation, double-precision mode, and no-interference. Modeling procedures and key technologies are compared and summarized: speed mode is superior to displacement mode in simulation accuracy and stability; time truncation is an undesirable option for analysis of the distribution trend of time-history parameters to guarantee that the data has reached the stable state; double-precision mode is more suitable for stagger spinning simulation, as truncation error has obvious effects on the accuracy of results; interference phenomenon can lead to obvious oscillation and mutation simulation results and influence the reliability of simulation significantly. Then, based on the modified model, some improvements of current reported results of roller intervals have been made, which lead to higher accuracy and reliability in the simulation.

Cover systems are used to prevent water infiltration into a waste body. They also play an important role in controlling landfill gas transport from the waste body to the atmosphere. It is important to assess the flux of landfill gas at the surface of a cover system by considering the coupled effects of rainwater infiltration and gas transport in the cover soils. We have developed a 1D mathematical model for coupled transient gas and water transport in unsaturated cover soils. The coupled model was solved by the finite element method. Results obtained by the proposed model agreed well with experimental data. Based on the proposed solution, the influences of gas pressure, gas permeability, and the thickness of the cover soils on soil gas concentration profiles were investigated. The difference in soil gas concentration reached up to 31% as the thickness of cover increased from 1 to 2 m. Gas concentration at a depth of 0.2 m decreased by 6% as the amplitude of atmospheric gas pressure fluctuation increased from 20 to 100 Pa. The gas concentration increased by only 3% when gas permeability increased by a factor of 2 for a relatively long period of gas migration (e.g., 60 h) under the given conditions. Results suggest that both diffusion and advection should be considered when estimating gas transport in unsaturated cover soils. The numerical model can be used in the design of cover systems in relation to gas breakthrough time, breakthrough concentration, and flux.