Numerical Simulation of Stefan’s One-Dimensional Multi-Front Problems

Author(s):  
Boris G. Aksenov ◽  
Yuri E. Karyakin
2020 ◽  
pp. 014459872098361
Author(s):  
Zhongbao Wu ◽  
Qingjun Du ◽  
Bei Wei ◽  
Jian Hou

Foam flooding is an effective method for enhancing oil recovery in high water-cut reservoirs and unconventional reservoirs. It is a dynamic process that includes foam generation and coalescence when foam flows through porous media. In this study, a foam flooding simulation model was established based on the population balance model. The stabilizing effect of the polymer and the coalescence characteristics when foam encounters oil were considered. The numerical simulation model was fitted and verified through a one-dimensional displacement experiment. The pressure difference across the sand pack in single foam flooding and polymer-enhanced foam flooding both agree well with the simulation results. Based on the numerical simulation, the foam distribution characteristics in different cases were studied. The results show that there are three zones during foam flooding: the foam growth zone, stable zone, and decay zone. These characteristics are mainly influenced by the adsorption of surfactant, the gas–liquid ratio, the injection rate, and the injection scheme. The oil recovery of polymer-enhanced foam flooding is estimated to be 5.85% more than that of single foam flooding. Moreover, the growth zone and decay zone in three dimensions are considerably wider than in the one-dimensional model. In addition, the slug volume influences the oil recovery the most in the foam enhanced foam flooding, followed by the oil viscosity and gas-liquid ratio. The established model can describe the dynamic change process of foam, and can thus track the foam distribution underground and aid in optimization of the injection strategies during foam flooding.


Author(s):  
Qiao Luo ◽  
Xiaobing Zhang

Purpose – In engineering applications, gas-solid two-phase reaction flow with multi-moving boundaries is a common phenomenon. The launch process of multiple projectiles is a typical example. The flow of adjacent powder chambers is coupled by projectile’s motion. The purpose of this paper is to study this flow by numerical simulation. Design/methodology/approach – A one-dimensional two-phase reaction flow model and MacCormack difference scheme are implemented in a computational code, and the code is used to simulate the launch process of a system of multiple projectiles. For different launching rates and loading conditions, the simulated results of the launch process of three projectiles are obtained and discussed. Findings – At low launching rates, projectiles fired earlier in the series have little effect on the launch processes of projectiles fired later. However, at higher launching rates, the projectiles fired first have a great influence on the launch processes of projectiles fired later. As the launching rate increases, the maximum breech pressure for the later projectiles increases. Although the muzzle velocities increase initially, they reach a maximum at some launching rate, and then decrease rapidly. The muzzle velocities and maximum breech pressures of the three projectiles have an approximate linear relationship with the charge weight, propellant web size and chamber volume. Originality/value – This paper presents a prediction tool to understand the physical phenomenon of the gas-solid two-phase reaction flow with multi-moving boundaries, and can be used as a research tool for future interior ballistics studies of launch system of multiple projectiles.


Author(s):  
Doug Garrard ◽  
Milt Davis ◽  
Steve Wehofer ◽  
Gary Cole

The NASA Lewis Research Center (LeRC) and the Arnold Engineering Development Center (AEDC) have developed a closely coupled computer simulation system that provides a one dimensional, high frequency inlet / engine numerical simulation for aircraft propulsion systems. The simulation system, operating under the LeRC-developed Application Portable Parallel Library (APPL), closely coupled a supersonic inlet with a gas turbine engine. The supersonic inlet was modeled using the Large Perturbation Inlet (LAPIN) computer code, and the gas turbine engine was modeled using the Aerodynamic Turbine Engine Code (ATEC). Both LAPIN and ATEC provide a one dimensional, compressible, time dependent flow solution by solving the one dimensional Euler equations for the conservation of mass, momentum, and energy. Source terms are used to model features such as bleed flows, turbomachinery component characteristics, and inlet subsonic spillage while unstarted. High frequency events, such as compressor surge and inlet unstart, can be simulated with a high degree of fidelity. The simulation system was exercised using a supersonic inlet with sixty percent of the supersonic area contraction occurring internally, and a GE J85-13 turbojet engine.


2000 ◽  
Author(s):  
Zhilong He ◽  
Xueyuan Peng ◽  
Pengcheng Shu

Abstract This paper presents a numerical method for simulating the thermal and fluid-dynamic behavior of hermetic compressors in the whole compressor domain. The model of fluid flow is developed by integrating transient one-dimensional conservation equations of continuity, momentum and energy through all of the elements from suction line to discharge line. The model describing thermal behavior is based on heat balance in the components such as muffler, connecting tubes and orifices. The calculation of the thermodynamic and transport properties for different refrigerants at various conditions has been considered, and some numerical results for a hermetic compressor are presented. The present study has demonstrated that the numerical simulation is a fest and reliable tool for compressor design.


Author(s):  
Azka Umar ◽  
Chun Jiang

This paper focuses on manipulating thermal emission and radiation loss of heat energy in a heat waveguide. A One-Dimensional Photonic Crystal is used as a waveguide clad to prohibit the thermal emission from escaping. The model may reduce the radiation loss of heat energy in the waveguide core, and heat energy can be confined to propagate along the waveguide’s longitude axis. The waveguide clad comprises alternative layers of high and low refractive index materials containing sufficient electromagnetic stop bands to trap the thermal emission from escaping out of the waveguide. The numerical simulation of the model shows that the forbidden bandgap of photonic crystal structures with alternative layers of silica and silicon has width enough to make heat energy be confined within the waveguide core so that efficient heat energy transmission can be achieved along the longitude axis of the waveguide.


Sign in / Sign up

Export Citation Format

Share Document