Effect of regularized functions on the dynamic response of a clutch system using a high-order algorithm

The main objective of this work is to propose some regularization techniques for modeling contact actions in a clutch system and to solve the obtained nonlinear dynamic problem by a high-order algorithm. This device is modeled by a discrete mechanical system with eleven degrees of freedom. In several works, the discontinuous models of the contact actions are replaced by the smoothed functions using the hyperbolic tangent. We propose, in this work, to replace the discontinuous model by a regularized model with new continuous functions that permit us to search the solution under Taylor series expansion. This regularized model approaches better the discontinuous model than the model based on the smoothing functions, especially in the vicinity of the zone of singularities. To solve the equations of motion of discrete mechanical systems, we propose to use a high-order algorithm combining a time discretization, a change of variable based on the previous time, a homotopy transformation and Taylor series expansion in the continuation process. The results obtained by this modeling are compared with those computed by the Newton–Raphson algorithm.     

Asymptotic Analysis of Three-Scale Model of pH-Dependent Flows in 1:1 Clays with Danckwerts’ Boundary Conditions

The electroremediation process is an efficient method for removing pollutants from clayey soils. We model this process in kaolinite clays considering three scales—nano, micro and macro—under the assumption of a stratified geometry in conjunction with more realistic Danckwerts’ boundary conditions imposed at the electrodes. The resulting multiscale model is a coupled system of nonlinear partial differential equations. We derive analytical solutions of the macroscopic equations considering the asymptotic behavior of strongly convective and diffusive regimes. We perform numerical simulations of different scenarios for the electroremediation using the Galerkin finite element method together with a staggered algorithm and the Newton–Raphson method. We validate the accuracy of the proposed algorithm by comparing the discrete solutions to the analytical ones. Finally, we explore and discuss optimal scenarios for the electroremediation process depending on the input values of pH, electrical current, and mass inflow using dimensionless numbers defined from the analytical solutions.     

Analysis of elastic layers with dilative eigenstrains varying through the thickness

IntroductionBrittlematerials,suchasceramicsandinorganicglasses ,oftenexhibittensilestrengthfarbelowtheircompressivestrength .Moreover,thesematerialsusuallyfailinanunstableandcatastrophicmannerwhensubjectedtoappliedloads.Toovercomethesedrawbacks,greatefforthasbeenmadeinthepastdecades.Astrategyistointroduceproperresidualcompressioninthematherialstosuperposethetensilestressandenhancethefractureresistance .Forexample,Greenetal.[1] haverecentlyproposedameansfortougheningglassplatesthroughtheuseofd…     

A multigrid finite difference approach to steady flow between eccentric rotating cylinders

A simple finite difference scheme over a non‐uniform grid is proposed to solve the two‐dimensional, steady Navier–Stokes equations. Instead of the Newton method, a more straightforward line search algorithm is used to solve the resultant system of non‐linear equations. By adopting the multigrid methodology, a fast convergence is achieved, at least for low‐Reynolds number flow. This scheme is applied, in particular, to flow between eccentric rotating cylinders. The computed results are shown in good agreement with some analytic findings. Copyright © 2000 John Wiley & Sons, Ltd.     

On the boundary layer methods

In this paper, the defect of the traditionary boundary layer methods (including the method of matched asymptotic expansions and the method of Visik-Lyusternik) is noted, from those methods we can not construct the asymptotic expansion of boundary layer term substantially. So the method of multiple scales is proposed for constructing the asymptotic expansion of boundary layer term, the reasonable result is obtained. Furthermore, we compare this method with the method used by Levinson, and find that both methods give the same asymptotic expansion of boundary layer term, but our method is simpler.Again, we apply this method to study some known works on singular perturbations. The limitations of those works have been noted, and the asymptotic expansion of solution is constructed in general condition.     

非稳态流动的隐式最小二乘等几何方法

对非稳态粘性不可压流动问题提出了一种隐式最小二乘等几何计算方法。该方法先用隐式的向后多步差分格式对Navier-Stokes方程进行时间离散,再用Newton法线性化对流项,最后在每个时间步上用最小二乘等几何方法进行求解。根据该算法编制了计算程序,通过构造解析解的方法验证了程序的正确性,用该程序求解了雷诺数为5000时的非稳态二维顶盖驱动流问题,计算结果捕捉到了流动过程中涡的演化过程,表明本文方法可用于非稳态流动的求解。     

牛顿迭代一致性算法及其在板弹塑性有限元分析中的应用

本文简略讨论了有限载荷增量弹塑性有限元分析中传统切线刚度法丧失精度和牛顿迭代平方收敛速度的原因,并提出保持牛顿迭代平方收敛速度、保证一阶精度和无条件稳定性的一致性算法.一致性算法具备以下两个特征:1)采用路径无关计算格式;2)采用一致弹塑性切线模量。根据一致性算法构造出以弯矩和曲率为基本变量的弹塑性板弯曲有限元NIDKQ元。数值结果表明NIDKQ元具有令人满意的精度,同时验证了有限载荷增量下牛顿迭代一致性算法的平方收敛率特性,而传统切线刚度法随着塑性区的扩展将大大降低收敛速度。     

Parallel domain decomposition method for finite element approximation of 3D steady state non‐Newtonian fluids

We introduce a stabilized finite element method for the 3D non‐Newtonian Navier–Stokes equations and a parallel domain decomposition method for solving the sparse system of nonlinear equations arising from the discretization. Non‐Newtonian flow problems are, generally speaking, more challenging than Newtonian flows because the nonlinearities are not only in the convection term but also in the viscosity term, which depends on the shear rate. Many good iterative methods and preconditioning techniques that work well for the Newtonian flows do not work well for the non‐Newtonian flows. We employ a Galerkin/least squares finite element method, with stabilization parameters adjusted to count the non‐Newtonian effect, to discretize the equations, and the resulting highly nonlinear system of equations is solved by a Newton–Krylov–Schwarz algorithm. In this study, we apply the proposed method to some inelastic power‐law fluid flows through the eccentric annuli with inner cylinder rotation and investigate the robustness of the method with respect to some physical parameters, including the power‐law index and the Reynolds number ratios. We then report the superlinear speedup achieved by the domain decomposition algorithm on a computer with up to 512 processors. Copyright © 2015 John Wiley & Sons, Ltd.     

A conservative semi‐implicit method for coupled surface–subsurface flows in regional scale

A semi‐implicit method for coupled surface–subsurface flows in regional scale is proposed and analyzed. The flow domain is assumed to have a small vertical scale as compared with the horizontal extents. Thus, after hydrostatic approximation, the simplified governing equations are derived from the Reynolds averaged Navier–Stokes equations for the surface flow and from the Darcy's law for the subsurface flow. A conservative free‐surface equation is derived from a vertical integral of the incompressibility condition and extends to the whole water column including both, the surface and the subsurface, wet domains. Numerically, the horizontal domain is covered by an unstructured orthogonal grid that may include subgrid specifications. Along the vertical direction a simple z‐layer discretization is adopted. Semi‐implicit finite difference equations for velocities and a finite volume approximation for the free‐surface equation are derived in such a fashion that, after simple manipulation, the resulting discrete free‐surface equation yields a single, well‐posed, mildly nonlinear system. This system is efficiently solved by a nested Newton‐type iterative method that yields simultaneously the pressure and a non‐negative fluid volume throughout the computational grid. The time‐step size is not restricted by stability conditions dictated by friction or surface wave speed. The resulting algorithm is simple, extremely efficient, and very accurate. Exact mass conservation is assured also in presence of wetting and drying dynamics, in pressurized flow conditions, and during free‐surface transition through the interface. A few examples illustrate the model applicability and demonstrate the effectiveness of the proposed algorithm. Copyright © 2015 John Wiley & Sons, Ltd.     

Discrete and non-local elastica

In this paper, the buckling and post-buckling behavior of an elastic lattice system referred to as the discrete elastica problem is investigated using an equivalent non-local continuum approach. The geometrically exact post-buckling analysis of the elastic chain, also called Hencky system, is first numerically solved using the shooting method. This discrete physical model is also mathematically equivalent to a finite difference formulation of the continuum elastica. Starting from the exact difference equations of the discrete problem, a continualization method is applied for approximating the difference operators by differential ones, in order to better characterize the discrete system by an enriched continuous one. It is shown that the new continuum associated with the discrete system exactly fits the discrete elastica post-buckling problem, where the non-locality is of Eringen׳s type (also called stress gradient non-local model). An asymptotic expansion is performed for both the discrete and the non-local continuum models, in order to approximate the post-buckling branches of the discrete system. Some numerical investigations show the efficiency of the non-local approach, especially for capturing the scale effects inherent to the cell size of the lattice model.     

Permanence and asymptotic properties of nonlinear delay difference equations

The asymptotic behavior of a class of nonlinear delay difference equation wax studied . Some sufficient conditions are obtained for permanence and global attractivity . The results can be applied to a claxs of nonlinear delay difference equations and to the delay discrete Logistic model and some known results are included.     

Low-rank solution of an optimal control problem constrained by random Navier-Stokes equations

We develop a low-rank tensor decomposition algorithm for the numerical solution of a distributed optimal control problem constrained by two-dimensional time-dependent Navier-Stokes equations with a stochastic inflow. The goal of optimization is to minimize the flow vorticity. The inflow boundary condition is assumed to be an infinite-dimensional random field, which is parametrized using a finite- (but high-) dimensional Fourier expansion and discretized using the stochastic Galerkin finite element method. This leads to a prohibitively large number of degrees of freedom in the discrete solution. Moreover, the optimality conditions in a time-dependent problem require solving a coupled saddle-point system of nonlinear equations on all time steps at once. For the resulting discrete problem, we approximate the solution by the tensor-train (TT) decomposition and propose a numerically efficient algorithm to solve the optimality equations directly in the TT representation. This algorithm is based on the alternating linear scheme (ALS), but in contrast to the basic ALS method, the new algorithm exploits and preserves the block structure of the optimality equations. We prove that this structure preservation renders the proposed block ALS method well posed, in the sense that each step requires the solution of a nonsingular reduced linear system, which might not be the case for the basic ALS. Finally, we present numerical experiments based on two benchmark problems of simulation of a flow around a von Kármán vortex and a backward step, each of which has uncertain inflow. The experiments demonstrate a significant complexity reduction achieved using the TT representation and the block ALS algorithm. Specifically, we observe that the high-dimensional stochastic time-dependent problem can be solved with the asymptotic complexity of the corresponding deterministic problem.     

A variational multiscale Newton–Schur approach for the incompressible Navier–Stokes equations

In the following paper, we present a consistent Newton–Schur (NS) solution approach for variational multiscale formulations of the time‐dependent Navier–Stokes equations in three dimensions. The main contributions of this work are a systematic study of the variational multiscale method for three‐dimensional problems and an implementation of a consistent formulation suitable for large problems with high nonlinearity, unstructured meshes, and non‐symmetric matrices. In addition to the quadratic convergence characteristics of a Newton–Raphson‐based scheme, the NS approach increases computational efficiency and parallel scalability by implementing the tangent stiffness matrix in Schur complement form. As a result, more computations are performed at the element level. Using a variational multiscale framework, we construct a two‐level approach to stabilizing the incompressible Navier–Stokes equations based on a coarse and fine‐scale subproblem. We then derive the Schur complement form of the consistent tangent matrix. We demonstrate the performance of the method for a number of three‐dimensional problems for Reynolds number up to 1000 including steady and time‐dependent flows. Copyright © 2009 John Wiley & Sons, Ltd.     

Thermally induced dynamic instability of laminated composite conical shells

Thermally induced dynamic instability of laminated composite conical shells is investigated by means of a perturbation method. The laminated composite conical shells are subjected to static and periodic thermal loads. The linear instability approach is adopted in the present study. A set of initial membrane stresses due to the elevated temperature field is assumed to exist just before the instability occurs. The formulation begins with three-dimensional equations of motion in terms of incremental stresses perturbed from the state of neutral equilibrium. After proper nondimensionalization, asymptotic expansion and successive integration, we obtain recursive sets of differential equations at various levels. The method of multiple scales is used to eliminate the secular terms and make an asymptotic expansion feasible. Using the method of differential quadrature and Bolotin's method, and imposing the orthonormality and solvability conditions on the present asymptotic formulation, we determine the boundary frequencies of dynamic instability regions for various orders in a consistent and hierarchical manner. The principal instability regions of cross-ply conical shells with simply supported–simply supported boundary conditions are studied to demonstrate the performance of the present asymptotic theory.     

Newton linearization of the incompressible Navier–Stokes equations

The present study aims to accelerate the convergence to incompressible Navier–Stokes solution. For the sake of computational efficiency, Newton linearization of equations is invoked on non‐staggered grids to shorten the sequence to the final solution of the non‐linear differential system of equations. For the sake of accuracy, the resulting convection–diffusion–reaction finite‐difference equation is solved line‐by‐line using the proposed nodally exact one‐dimensional scheme. The matrix size is reduced and, at the same time, the CPU time is considerably saved due to the decrease of stencil points. The effectiveness of the implemented Newton linearization is demonstrated through computational exercises. Copyright © 2004 John Wiley & Sons, Ltd.     

Wavelet multiresolution interpolation Galerkin method for nonlinear boundary value problems with localized steep gradients

The wavelet multiresolution interpolation for continuous functions defined on a finite interval is developed in this study by using a simple alternative of transformation matrix. The wavelet multiresolution interpolation Galerkin method that applies this interpolation to represent the unknown function and nonlinear terms independently is proposed to solve the boundary value problems with the mixed Dirichlet-Robin boundary conditions and various nonlinearities, including transcendental ones, in which the discretization process is as simple as that in solving linear problems, and only common two-term connection coefficients are needed. All matrices are independent of unknown node values and lead to high efficiency in the calculation of the residual and Jacobian matrices needed in Newton’s method, which does not require numerical integration in the resulting nonlinear discrete system. The validity of the proposed method is examined through several nonlinear problems with interior or boundary layers. The results demonstrate that the proposed wavelet method shows excellent accuracy and stability against nonuniform grids, and high resolution of localized steep gradients can be achieved by using local refined multiresolution grids. In addition, Newton’s method converges rapidly in solving the nonlinear discrete system created by the proposed wavelet method, including the initial guess far from real solutions.

     

A continuum model for size-dependent deformation of elastic films of nano-scale thickness

Ultra-thin elastic films of nano-scale thickness with an arbitrary geometry and edge boundary conditions are analyzed. An analytical model is proposed to study the size-dependent mechanical response of the film based on continuum surface elasticity. By using the transfer-matrix method along with an asymptotic expansion technique of small parameter, closed-form solutions for the mechanical field in the film is presented in terms of the displacements on the mid-plane. The asymptotic expansion terminates after a few terms and exact solutions are obtained. The mid-plane displacements are governed by three two-dimensional equations, and the associated edge boundary conditions can be prescribed on average. Solving the two-dimensional boundary value problem yields the three-dimensional response of the film. The solution is exact throughout the interior of the film with the exception of a thin boundary layer having an order of thickness as the film in accordance with the Saint-Venant’s principle.     

Residual a posteriori error estimate two-grid methods for the steady Navier-Stokes equation with stream function form

Residual based on a posteriori error estimates for conforming element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.     

Efficient computation of steady solitary gravity waves

An efficient numerical method to compute solitary wave solutions to the free surface Euler equations is reported. It is based on the conformal mapping technique combined with an efficient Fourier pseudo-spectral method. The resulting nonlinear equation is solved via the Petviashvili iterative scheme. The computational results are compared to some existing approaches, such as Tanaka’s method and Fenton’s high-order asymptotic expansion. Several important integral quantities are computed for a large range of amplitudes. The integral representation of the velocity and acceleration fields in the bulk of the fluid is also provided.     

Analysis of a two-grid method for semiconductor device problem

The mathematical model of a semiconductor device is governed by a system of quasi-linear partial differential equations.The electric potential equation is approximated by a mixed finite element method,and the concentration equations are approximated by a standard Galerkin method.We estimate the error of the numerical solutions in the sense of the Lqnorm.To linearize the full discrete scheme of the problem,we present an efficient two-grid method based on the idea of Newton iteration.The main procedures are to solve the small scaled nonlinear equations on the coarse grid and then deal with the linear equations on the fine grid.Error estimation for the two-grid solutions is analyzed in detail.It is shown that this method still achieves asymptotically optimal approximations as long as a mesh size satisfies H=O(h^1/2).Numerical experiments are given to illustrate the efficiency of the two-grid method.