Modified Reynolds Equation for Squeeze-Film Air Damping of Slotted Plates in MEMS Devices

A differential equation for calculating squeeze-film air damping in slotted plates is developed by modifying the Reynolds equation. A term is added to account for the effect of airflow through the slots on the air damping of the plate. The end effect of the airflow in the slots is also treated by substituting an effective channel length for the geometric channel length （i. e. the thickness of the plate）..The damping pressure distribution, damping force, and damping force coefficient of the slotted plates can be found by solving the equation under appropriate boundary conditions. With restrictions on the thickness and the lateral dimensions of the slotted plate removed,the equation provides a useful tool for analysing the squeeze-film air damping effect of slotted plates with finite thickness and finite lateral dimensions. For a typical slotted plate structure, the damping force coefficient obtained by this equation agrees well with that generated by ANSYS.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.

Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

Derivation of a Revised Tsiolkovsky Rocket Equation That Predicts Combustion Oscillations

Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx1, and slow-jet, vx2. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.

High-Order Spatial FDTD Solver of Maxwell’s Equations for Terahertz Radiation Production

We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.

MODELING THE INTERACTION OF SOLITARY WAVES AND SEMI-CIRCULAR BREAKWATERS BY USING UNSTEADY REYNOLDS EQUATIONS

A vertical 2-D numerical wave model was developed based on unsteady Reynolds equations. In this model, the k-epsilon models were used to close the Reynolds equations, and volume of fluid(VOF) method was used to reconstruct the free surface. The model was verified by experimental data. Then the model was used to simulate solitary wave interaction with submerged, alternative submerged and emerged semi-circular breakwaters. The process of velocity field, pressure field and the wave surface near the breakwaters was obtained. It is found that when the semi-circular breakwater is submerged, a large vortex will be generated at the bottom of the lee side wall of the breakwater; when the still water depth is equal to the radius of the semi-circular breakwater, a pair of large vortices will be generated near the shoreward wall of the semi-circular breakwater due to wave impacting, but the velocity near the bottom of the lee side wall of the breakwater is always relatively small. When the semi-circular breakwater is emerged, and solitary wave cannot overtop it, the solitary wave surface will run up and down secondarily during reflecting from the breakwater. It can be further used to estate the diffusing and transportation of the contamination and transportation of suspended sediment.

A GENERALIZED REYNOLDS EQUATION BASED ON NON-NEWTONIAN FLOW IN LUBRICATION MECHANICS

A generalized Reynolds equation based on non-Newtonian flow is derived in this paper.This equation is suitable for a number of non-Newtonian flow models and can be solved numerically to obtain pressure fields in thermalhydrodynamically or elastohydrodynamically lubricated fluid films.A mathematical ap- proach is given for solving simultaneously the shearing stress,shearing rate,velocity and equivalent viscosity.To show the application of this equation,two rheological models which have been widely used in lubrication mechnaics are incorporated into this equation to obtain numerical solutions to the line contact thermal elastohydrodynamic lubrication problem.

Travelling Waves: Interplay of Low to High Reynolds Number and Tan-Cot Function Method to Solve Burger’s Equations

We study the nonlinear parabolic equations for travelling wave solutions of Burger’s equations. The purpose of the present work is to study various types of Burger’s equations describing waves and those are based on nonlinear equations. We focus on to describe the analytic solution in the special pattern of travelling wave solutions using tan-cot function method. We discuss about inviscid and viscous version of Burger’s equation for fluid flow and investigate the effects of internal friction of a fluid via Reynolds number. By changing the velocity amplitude, the nature of flows with shock wave and disturbance are observed. For numerical solutions, the Crank-Nicolson scheme is introduced to establish the wave solutions.

Water Profile over an Ogee-Shaped Crest and the Reynolds Fractional Equations

Weir crest must have the correct shape in the concave side of an ogee-shaped crest to diminish erosion. This shape can be obtained using an approximation of the fractional Reynolds equations when the water interacts with the surface. A model is introduced for the Reynolds stresses complemented with a closure relation of fractional origin. A power type solution is obtained for the main velocity and stress. Velocity profile is found based on the assumption of a steady flow and the energy conservation equation. A Froude number and the cubic equation of the weir are built. The dimensionless upstream velocity head is also determined which allow graphically showing the exponent and coefficient of the water-profile over an ogee-shaped crest. It is possible to estimate the occupied-space index trough an exponents＇ ratio of profile over the velocity head.

Review of Kelvin's Equation and Its Modification in Characterization of Mesoporous Materials

Physical and mathematical models as well as calculation methods of nitrogen bed on porous media have been introduced to evaluate the structural parameters of mesoporous materials. Kelvin＇s equation is a link between the relative adsorbate pressure, the mean pore radius, and pore capillarity on the basis of macroscopic capillary condensation. However, Kelvin＇s equation has been identified that it underestimates the calculated pore size of a material especially in the boundary of pore size which is between 2 and 4 nm.Various modifications on Kelvin＇s equation were mentioned in order to develop a new model to improve the accuracy of pore size calculation. The problems on conventional mathematical models were analyzed and discussed. A number of calculation methods on physisorption and pore size, especially fundamental theories of physisorption, basis of models and their deficiencies are reviewed. It can provide guidance on developing a modified Kelvin＇s equation for pore size calculation.

Efficient Cell-Centered Multigrid Scheme for the Three-Dimensional Navier-Stokes Equations

A cell centered scheme for three dimensional Navier Stokes equations, which is based on central difference approximations and Runge Kutta time stepping, is described. By using local time stepping, implicit residual smoothing, a multigrid method, and carefully controlled artificial dissipative terms, good convergence rates are obtained for two and three dimensional flows. The emphases are on the implicit smoothing and artificial dissipative terms with locally variable coefficients which depend on cel...

PROPER SPLITTINGS FOR RESTRICTED LINEAR EQUATIONS AND THE GENERALIZED INVERSE A_(T,S)~2

This paper presents a proper splitting iterative method for comparing the general restricted linear euqations Ax=b, x ∈T （where, b ∈AT, and T is an arbitrary but fixed subspace of Cm） and the generalized in AT,S For the special case when b ∈AT and dim（T）=dim（AT）, this splitting iterative methverse AT,S hod converges to AT,Sb （the unique solution of the general restricted system Ax=bx ∈T）.

Exact invariants and adiabatic invariants of Raitzin＇s canonical equations of motion for a nonlinear nonholonomic mechanical system

The exact invariants and the adiabatic invariants of Raitzin＇s canonical equations of motion for a nonlinear nonholonomic mechanical system are studied. The relations between the invariants and the symmetries of the system are established. Based on the concept of higher-order adiabatic invariant of a mechanical system under the action of a small perturbation, the forms of the exact invariants and adiabatic invariants and the conditions for their existence are proved. Finally, the inverse problem of the perturbation to symmetries of the system is studied and an example is also given to illustrate the application of the results.

Local projection stabilized finite element method for Navier-Stokes equations

This paper extends the results of Matthies, Skrzypacz, and Tubiska for the Oseen problem to the Navier-Stokes problem. For the stationary incompressible Navier- Stokes equations, a local projection stabilized finite element scheme is proposed. The scheme overcomes convection domination and improves the restrictive inf-sup condition. It not only is a two-level approach but also is adaptive for pairs of spaces defined on the same mesh. Using the approximation and projection spaces defined on the same mesh, the scheme leads to much more compact stencils than other two-level approaches. On the same mesh, besides the class of local projection stabilization by enriching the approximation spaces, two new classes of local projection stabilization of the approximation spaces are derived, which do not need to be enriched by bubble functions. Based on a special interpolation, the stability and optimal prior error estimates are shown. Numerical results agree with some benchmark solutions and theoretical analysis very well.

A two-grid algorithm based on Newton iteration for the stream function form of the Navier-Stokes equations

In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems （with same stiffness matrix but different right-hand side） on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the Mgorithm produces a numerical solution with the optimal asymptotic H^2-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.

A new full discrete stabilized viscosity method for transient Navier-Stokes equations

A new full discrete stabilized viscosity method for the transient Navier-Stokes equations with the high Reynolds number （small viscosity coefficient） is proposed based on the pressure projection and the extrapolated trapezoidal rule. The transient Navier-Stokes equations are fully discretized by the continuous equal-order finite elements in space and the reduced Crank-Nicolson scheme in time. The new stabilized method is stable and has many attractive properties. First, the system is stable for the equal-order combination of discrete continuous velocity and pressure spaces because of adding a pres- sure projection term. Second, the artifical viscosity parameter is added to the viscosity coefficient as a stability factor, so the system is antidiffusive. Finally, the method requires only the solution to a linear system at every time step. Stability and convergence of the method is proved. The error estimation results show that the method has a second-order accuracy, and the constant in the estimation is independent of the viscosity coefficient. The numerical results are given, which demonstrate the advantages of the method presented.

Finite difference streamline diffusion method using nonconforming space for incompressible time-dependent Navier-Stokes equations

This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and the Crouzeix-Raviart （CR） element combined with the P0 element in space are used. The result shows that this scheme has good stabilities and error estimates independent of the viscosity coefficient.

EXISTENCE OF SOLUTION FOR BOUNDARY VALUE PROBLEM OF NONLINEAR FRACTIONAL DIFFERENTIAL EQUATION

This paper is concerned with the boundary value problem of a nonlinear fractional differential equation. By means of Schauder fixed-point theorem, an existence result of solution is obtained.

Generalized Extended tanh-function Method for Traveling Wave Solutions of Nonlinear Physical Equations

In this paper, the generalized extended tanh-function method is used for constructing the traveling wave solutions of nonlinear evolution equations. We choose Fisher＇s equation, the nonlinear schr＆#168;odinger equation to illustrate the validity and ad-vantages of the method. Many new and more general traveling wave solutions are obtained. Furthermore, this method can also be applied to other nonlinear equations in physics.