Finite difference automatically generated non-uniform grids in the numerical solution of the Reynolds' compressible one-dimensional squeeze-film equation

2003 ◽  
Vol 217 (3) ◽  
pp. 243-249 ◽  
Author(s):  
E. A. M. Almas ◽  
F. A. P. De Silva
Keyword(s):  
Finite Difference ◽  
Finite Difference Methods ◽  
Computation Time ◽  
Squeeze Film ◽  
Uniform Grid ◽  
Efficient Computation ◽  
Uniform Grids ◽  
Uniform Model ◽  
Squeeze Number ◽  
Solution Accuracy

With the Reynolds equation, for compressible squeeze-film thrust bearings, the use of a finite difference discretization with a grid of nodes having a non-uniform spacing can result in a more efficient computation of the solution. Comparative tests of two variable spatial grid models against a uniform model were conducted with the classical finite difference methods for chosen combinations of squeeze number and excursion ratio values. For problems with a high value of σ, one of the non-uniform grid models has shown some advantages over the uniform model, requiring a smaller number of nodes and less computation time with the same solution accuracy.

Dynamics of Multilayer Microplates Considering Nonlinear Squeeze Film Damping

2006 ◽  
Author(s):  
M. T. Ahmadian ◽  
M. Moghimi Zand ◽  
H. Borhan
Keyword(s):  
Finite Element ◽  
Finite Difference ◽  
Reynolds Equation ◽  
Finite Difference Methods ◽  
Shear Deformation Theory ◽  
Single Layer ◽  
Squeeze Film ◽  
First Order ◽  
Squeeze Film Damping ◽  
Good Agreement

This paper presents a model to analyze pull-in phenomenon and dynamics of multi layer microplates using coupled finite element and finite difference methods. First-order shear deformation theory is used to model dynamical system using finite element method, while Finite difference method is applied to solve the nonlinear Reynolds equation of squeeze film damping. Using this model, Pull-in analysis of single layer and multi layer microplates are studied. The results of pull-in analysis are in good agreement with literature. Validating our model by pull-in results, an algorithm is presented to study dynamics of microplates. These simulations have many applications in designing multi layer microplates.

Compressible Squeeze Films and Squeeze Bearings

1964 ◽  
Vol 86 (2) ◽  
pp. 355-364 ◽  
Author(s):  
E. O. J. Salbu
Keyword(s):  
Steady State ◽  
Finite Difference ◽  
Journal Bearings ◽  
Squeeze Film ◽  
Thrust Bearings ◽  
Finite Difference Solution ◽  
Squeeze Number ◽  
Gas Bearing ◽  
Compressibility Effects ◽  
Piezoelectric Devices

Experimental agreement with a finite-difference solution of the isothermal squeeze film equation was obtained for steady-state sinusoidal squeeze motion of parallel, coaxial disks. At low squeeze number, the film force is in phase with the velocity; at high squeeze number, with the displacement. Compressibility effects at high squeeze number introduce a superambient mean film pressure, so that it is possible to operate a gas bearing on squeeze effects alone. Thrust bearings, spherical bearings, and journal bearings have been successfully operated as squeeze bearings, using both electromagnetic and piezoelectric devices to generate the squeeze motion.

COMPARISON OF FINITE VOLUME AND FINITE DIFFERENCE METHODS AND APPLICATION

2006 ◽  
Vol 04 (02) ◽  
pp. 163-208 ◽  
Author(s):  
S. FAURE ◽  
D. PHAM ◽  
R. TEMAM
Keyword(s):  
Finite Difference ◽  
Boundary Value Problems ◽  
Finite Volume ◽  
Finite Difference Methods ◽  
Finite Volume Methods ◽  
Uniform Grid ◽  
Parabolic Boundary ◽  
Convergence Results ◽  
Difference Methods ◽  
Parabolic Boundary Value Problems

In this article, we consider finite volume methods based on a non-uniform grid. Finite volume methods are compared to finite difference methods based on a related grid. As an application, various convergence results are proved for the finite volume function spaces and for some model elliptic and parabolic boundary value problems using these discretization spaces.

Cascadic Newton’s method for the elliptic Monge–Ampère equation

2020 ◽  
Vol 18 (03) ◽  
pp. 2050018
Author(s):  
Qin Li ◽  
Zhiyong Liu
Keyword(s):  
Finite Difference ◽  
Newton’S Method ◽  
Newton's Method ◽  
Finite Difference Methods ◽  
Computation Time ◽  
Type Interpolation ◽  
Local Type ◽  
Cascadic Multigrid ◽  
Monge Ampere ◽  
Prolongation Operator

In this paper, a cascadic Newton’s method is designed to solve the Monge–Ampère equation. In the process of implementing the cascadic multigrid, we use the Full-Local type interpolation as prolongation operator and Newton iteration as smoother. In order to obtain Full-Local type interpolation, we provide several finite difference stencils. Especially, the skewed finite difference methods are first applied by us for the elliptic Monge–Ampère equation. Based on Full-Local interpolation techniques and cascade principle, the new algorithm can save a large amount of computation time. Some numerical experiments are provided to confirm the efficiency of our proposed method.

Reversed time migration in spatial frequency domain

Geophysics ◽  
1983 ◽  
Vol 48 (5) ◽  
pp. 627-635 ◽  
Author(s):  
Dan Loewenthal ◽  
Irshad R. Mufti
Keyword(s):  
Finite Difference ◽  
Acoustic Waves ◽  
Seismic Waves ◽  
Finite Difference Methods ◽  
Computation Time ◽  
Real Data ◽  
Coordinate Frame ◽  
Seismic Section ◽  
Time Migration ◽  
Similar Scheme

During the past decade, finite‐difference methods have become important tools for direct modeling of seismic data as well as for certain interpretation processes. One of the earliest applications of these methods to seismics is the pioneering contribution of Alterman who, in a series of papers (Alterman and Karal, 1968; Alterman and Aboudi, 1968; Alterman and Rotenberg, 1969; Alterman and Loewenthal, 1972) demonstrated the usefulness of such numerical computations for the propagation of seismic waves in elastic media. A clear exposition of these techniques, as well as a comparison of results obtained from them with the corresponding analytical solutions, can be found in Alterman and Karal (1968). This subject was further developed and extended to more complicated models by Boore (1970), Ottaviani (1971), and Kelly et al (1976). Claerbout introduced a somewhat different finite‐difference approach (Claerbout, 1970; Claerbout and Johnson, 1971) for modeling the acoustic waves which often dominate the reflection seismogram. In his approach, the original wave equation, which governs the propagation of the acoustic waves, is modified in such a way so as to allow the propagation of either only upcoming or only downgoing waves. By moving the coordinate frame with the downgoing waves, Claerbout showed that one could greatly reduce computation time. Using the same concepts, he showed (Claerbout and Doherty, 1972) how to use a similar scheme for migrating a seismic section by downward continuation of the upcoming waves. This migration method is an interesting extension of the ideas of Hagedoorn (1954) and was found to be extremely useful with real data (Larner and Hatton, 1976; Loewenthal et al, 1976).

Third Order ENO Scheme on Non-Uniform Grid for Supersonic Flows II

2015 ◽  
Vol 798 ◽  
pp. 518-522
Author(s):  
Yekaterina Moisseyeva ◽  
Altynshash Naimanova ◽  
Asel Beketaeva
Keyword(s):  
Numerical Modeling ◽  
Finite Difference ◽  
Numerical Experiments ◽  
Mixing Layer ◽  
Shock Capturing ◽  
Uniform Grid ◽  
Third Order ◽  
Uniform Grids ◽  
Optimal Function ◽  
Eno Scheme

In the present paper we extend our earlier work on the finite-difference shock-capturing essentially non-oscillatory (ENO) scheme for a non-uniform grid. The design of the ENO scheme is based on the methodology for uniform grids. We provide the analysis of the different variations of the limiter functions for the developed algorithm to define the optimal function producing the smallest spread of the solution. The effect of the limiter choice on the mixing layer dynamics was studied for the non-uniform grid. The numerical experiments revealed that the nonoptimal choice of limiter can result in the overgrowth of the mixing layer, that is important for the numerical modeling of the combustion.

scholarly journals Modified Finite Difference Schemes on Uniform Grids for Simulations of the Helmholtz Equation at Any Wave Number

2014 ◽  
Vol 2014 ◽  
pp. 1-9 ◽  
Author(s):  
Hafiz Abdul Wajid ◽  
Naseer Ahmed ◽  
Hifza Iqbal ◽  
Muhammad Sarmad Arshad
Keyword(s):  
Finite Difference ◽  
Helmholtz Equation ◽  
Bloch Wave ◽  
Difference Schemes ◽  
Wave Number ◽  
Uniform Grid ◽  
Finite Difference Schemes ◽  
Radiation Boundary Conditions ◽  
Fine Mesh ◽  
Uniform Grids

We construct modified forward, backward, and central finite difference schemes, specifically for the Helmholtz equation, by using the Bloch wave property. All of these modified finite difference approximations provide exact solutions at the nodes of the uniform grid for the second derivative present in the Helmholtz equation and the first derivative in the radiation boundary conditions for wave propagation. The most important feature of the modified schemes is that they work for large as well as low wave numbers, without the common requirement of a very fine mesh size. The superiority of the modified finite difference schemes is illustrated with the help of numerical examples by making a comparison with standard finite difference schemes.

scholarly journals Staggered-Grid Finite Difference Method with Variable-Order Accuracy for Porous Media

2013 ◽  
Vol 2013 ◽  
pp. 1-10 ◽  
Author(s):  
Jinghuai Gao ◽  
Yijie Zhang
Keyword(s):  
Porous Media ◽  
Finite Difference ◽  
Finite Difference Methods ◽  
Computation Time ◽  
High Order ◽  
Staggered Grid ◽  
Difference Operators ◽  
Variable Order ◽  
Low Velocity ◽  
Order Accuracy

The numerical modeling of wave field in porous media generally requires more computation time than that of acoustic or elastic media. Usually used finite difference methods adopt finite difference operators with fixed-order accuracy to calculate space derivatives for a heterogeneous medium. A finite difference scheme with variable-order accuracy for acoustic wave equation has been proposed to reduce the computation time. In this paper, we develop this scheme for wave equations in porous media based on dispersion relation with high-order staggered-grid finite difference (SFD) method. High-order finite difference operators are adopted for low-velocity regions, and low-order finite difference operators are adopted for high-velocity regions. Dispersion analysis and modeling results demonstrate that the proposed SFD method can decrease computational costs without reducing accuracy.

Application of Backus average to simulate acoustic logs in thinly bedded formations

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. D93-D104
Author(s):  
Elsa Maalouf ◽  
Carlos Torres-Verdín ◽  
Jingxuan Li
Keyword(s):  
Finite Difference ◽  
Numerical Simulations ◽  
Finite Difference Methods ◽  
Homogeneous Medium ◽  
Computation Time ◽  
Transversely Isotropic ◽  
Simulation Method ◽  
Thin Layers ◽  
Homogeneous Layer ◽  
Spatial Averaging

Slowness logs acquired in layered formations are not only affected by spatial averaging associated with the borehole acoustic tool. Layers with thicknesses smaller than the acoustic wavelength can cause measurable effects on the associated wave propagation phenomena. While spatial averaging functions can be used to model tool averaging effects, computer-intensive numerical methods such as finite differences must be used to simulate slowness logs across formations with thin layers. We adopted Backus averaging as a faster alternative to model borehole slownesses when layer thicknesses are smaller than the acoustic wavelength (i.e., in the long-wavelength limit). Using synthetic models and numerical simulations via finite-element and finite-difference methods, we have determined that borehole slownesses of a stack of horizontal layers first approach the average slowness of the individual layers. However, as the layer thickness decreases, sonic slownesses approach the slowness of a homogeneous medium with elastic properties obtained from the Backus average. Therefore, to model acoustic logs acquired in layered formations, we first approximate thin layers as a single homogeneous layer with stiffness coefficients calculated using the Backus average. Next, we apply a spatial averaging function to reproduce the spatial averaging effect inherent to borehole acoustic tools. Results indicate that the latter method is accurate and efficient for fast modeling borehole slownesses of formations with thin layers that are isotropic and intrinsically vertical transversely isotropic. The fast simulation method decreases computation time by at least a factor of 10 and yields slowness logs with a relative error below 2% compared with finite-difference numerical simulations. We also determine that the moving Backus average that is typically applied to upscale acoustic logs for seismic applications is not accurate to model borehole acoustic logs acquired across thinly layered formations.

Sign in / Sign up

Export Citation Format

Share Document