Finite difference method

Abstract: to seismic migration.- Harmonic analysis, delta functions, and all that.- Equations of motion for the earth.- Elastic wave equations.- Ray theory.- Kirchhoff migration.- Kirchhoff migration/inversion.- The method of stationary phase.- Downward continuation of the seismic wavefield.- Plane wave decomposition of seismograms.- Numerical methods for tracing rays.- Finite difference methods for wave propagation and migration.- SU user's manual.- SUB user's guide.

Abstract: A numerical simulation of a pressure swing adsorption process is presented for a system in which a small concentration of an adsorbable component is separated from an inert carrier. Linear equilibrium and a linear rate expression are assumed. The model equations were solved by orthogonal collocation and by finite difference methods with consistent results. The theory is shown to provide a good representation of the experimental data of Mitchell and Shendalman (1973) for the system CO2-He-silica gel.

Abstract: The Kelvin–Helmholtz instability of multi-component (more than two) incompressible and immiscible fluids is studied numerically using a phase-field model. The instability is governed by the modified Navier–Stokes equations and the multi-component convective Cahn–Hilliard equations. A finite difference method is used to discretize the governing system. To solve the equations efficiently and accurately, we employ the Chorin’s projection method for the modified Navier–Stokes equations and the recently developed practically unconditionally stable method for the multi-component Cahn–Hilliard equations. Through our model and numerical solution, we investigate the effects of surface tension, density ratio, magnitude of velocity difference, and forcing on the Kelvin–Helmholtz instability of multi-component fluids. It is shown that increasing the surface tension or the density ratio reduces the growth of the Kelvin–Helmholtz instability. And it is also observed that as the initial horizontal velocity difference gets larger, the interface rolls up more. We also found that the billow height reaches its maximum more slowly as the initial wave amplitude gets smaller. And, for the linear growth rate for the Kelvin–Helmholtz instability of two-component fluids, the simulation results agree well with the analytical results. From comparison between the numerical growth rate of two- and three-component fluids, we observe that the inclusion of extra layers can alter the growth rate for the Kelvin–Helmholtz instability. Finally, we simulate the billowing cloud formation which is a classic example of the Kelvin–Helmholtz instability and cannot be seen in binary fluids. With our multi-component method, the details of the real flow (e.g., the asymmetry in the roll-up and the self-interaction of the shear layer) are well captured.