New Approaches to the Propagation of Nonlinear Transients in Porous Media

Many problems in regional groundwater flow require the characterization and forecasting of variables, such as hydraulic heads, hydraulic gradients, and pore velocities. These variables describe hydraulic transients propagating in an aquifer, such as a river flood wave induced through an adjacent aquifer. The characterization of aquifer variables is usually accomplished via the solution of a transient differential equation subject to time-dependent boundary conditions. Modeling nonlinear wave propagation in porous media is traditionally approached via numerical solutions of governing differential equations. Temporal or spatial numerical discretization schemes permit a simplification of the equations. However, they may generate instability, and require a numerical linearization of true nonlinear problems. Traditional analytical solutions are continuous in space and time, and render a more stable solution, but they are usually applicable to linear problems and require regular domain shapes. The method of decomposition of Adomian is an approximate analytical series to solve linear or nonlinear differential equations. It has the advantages of both analytical and numerical procedures. An important limitation is that a decomposition expansion in a given coordinate explicitly uses the boundary conditions in such axis only, but not necessarily those on the others. In this article we present improvements of the method consisting of a combination of a partial decomposition expansion in each coordinate in conjunction with successive approximation that permits the consideration of boundary conditions imposed on all of the axes of a transient multidimensional problem; transient modeling of irregularly-shaped aquifer domains; and nonlinear transient analysis of groundwater flow equations. The method yields simple solutions of dependent variables that are continuous in space and time, which easily permit the derivation of heads, gradients, seepage velocities and fluxes, thus minimizing instability. It could be valuable in preliminary analysis prior to more elaborate numerical analysis. Verification was done by comparing decomposition solutions with exact analytical solutions when available, and with controlled experiments, with reasonable agreement. The effect of linearization of mildly nonlinear saturated groundwater equations is to underestimate the magnitude of the hydraulic heads in some portions of the aquifer. In some problems, such as unsaturated infiltration, linearization yields incorrect results.     

An analytical approximate technique for a class of strongly non-linear oscillators

An analytical approximate technique for large amplitude oscillations of a class of conservative single degree-of-freedom systems with odd non-linearity is proposed. The method incorporates salient features of both Newton's method and the harmonic balance method. Unlike the classical harmonic balance method, accurate analytical approximate solutions are possible because linearization of the governing differential equation by Newton's method is conducted prior to harmonic balancing. The approach yields simple linear algebraic equations instead of non-linear algebraic equations without analytical solution. With carefully constructed iterations, only a few iterations can provide very accurate analytical approximate solutions for the whole range of oscillation amplitude beyond the domain of possible solution by the conventional perturbation methods or harmonic balance method. Three examples including cubic-quintic Duffing oscillators are presented to illustrate the usefulness and effectiveness of the proposed technique.     

Approximate solutions for the problem of liquid film flow over an unsteady stretching sheet with thermal radiation and magnetic field

The proposed method is based on replacement of the unknown function by a truncated series of the shifted Legendre polynomial expansion. An approximate formula of the integer derivative is introduced. Special attention is given to study the convergence analysis and derive an upper bound of the error for the presented approximate formula. The introduced method converts the proposed equation by means of collocation points to a system of algebraic equations with shifted Legendre coefficients. Thus, after solving this system of equations, the shifted Legendre coefficients are obtained. This efficient numerical method is used to solve the system of ordinary differential equations which describe the thin film flow and heat transfer with the effects of the thermal radiation, magnetic field, and slip velocity.     

Thermoelastically coupled axisymmetric nonlinear vibration of shallow spherical and conical shells

The problem of axisymmetric nonlinear vibration for shallow thin spherical and conical shells when temperature and strain fields are coupled is studied. Based on the large deflection theories of yon Ktirrntin and the theory of thermoelusticity, the whole governing equations and their simplified type are derived. The time-spatial variables are separated by Galerkin ‘ s technique, thus reducing the governing equations to a system of time-dependent nonlinear ordinary differential equation. By means of regular perturbation method and multiple-scales method, the first-order approximate analytical solution for characteristic relation of frequency vs amplitude parameters along with the decay rate of amplitude are obtained, and the effects of different geometric parameters and coupling factors us well us boundary conditions on thermoelustically coupled nonlinear vibration behaviors are discussed.     

Falkner-Skan方程的近似解析解

研究了粘性流体绕流楔型物体的Falkner-Skan边界层方程求解问题.利用Adomian拆分方法,通过引入Crocco变量变换将无穷区间的边界值问题转为初值问题并利用Padé逼近技巧确定初值,给出了一种有效的解析分解方法.进一步,本文设计了一种数值解法,将本文得到的近似解析解及数值结果与早期研究者Hartree等人的结果进行了比较,证明了本文提出的解法的有效性和可靠性.     

PATH INTEGRAL SOLUTION OF NONLINEAR DYNAMIC BEHAVIOR OF STRUCTURE UNDER WIND EXCITATION

IntroductionThe dynamic behavior of the nonlinear structure under wind excitation has beenobserved very complicated.Taking guyed masts as an example,only a few collapsingaccidents occurred under extreme atmospheric conditions[1],many took place under mild…     

A linearization technique for multi-species transport problems

We study a one-dimensional multi-species system of dispersive-advective contaminant transport equations coupled by nonlinear biological (kinetic reactions) and physical (adsorption) processes. To deal with the nonlinearities and the coupling, and to avoid additional computational costs, we propose a linearization technique based on first-order Taylor’s series expansions. A stabilized finite element in space, combined with an Euler implicit finite difference discretization in time, is used to approximate the dispersive-advective transport problem. Three computational tests are performed with different boundary conditions, retardation factors and kinetic parameters for a nonlinear reactive multi-species transport model. The proposed methodology is shown to be accurate and decrease computational costs in the numerical implementation of nonlinear reactive transport problems.     

The influence of a magnetic field on free convection in a finite vertical channel

An approximate analytical solution is presented for developing free convection flows of electrically conducting fluids between finite vertical channels which are subjected to a uniformly applied transverse magnetic field. Specifically, the basic approximation lies in the linearization of the governing boundary layer type of equations. It is demonstrated that the application of a transverse magnetic field reduces the induced flow rate in the channel and the heat transfer to the fluid.     

Finite element solution of the incompressible Navier-Stokes equations by a Helmholtz velocity decomposition

Finite element solution methods for the incompressible Navier-Stokes equations in primitive variables form are presented. To provide the necessary coupling and enhance stability, a dissipation in the form of a pressure Laplacian is introduced into the continuity equation. The recasting of the problem in terms of pressure and an auxiliary velocity demonstrates how the error introduced by the pressure dissipation can be totally eliminated while retaining its stabilizing properties. The method can also be formally interpreted as a Helmholtz decomposition of the velocity vector. The governing equations are discretized by a Galerkin weighted residual method and, because of the modification to the continuity equation, equal interpolations for all the unknowns are permitted. Newton linearization is used and at each iteration the linear algebraic system is solved by a direct solver. Convergence of the algorithm is shown to be very rapid. Results are presented for two-dimensional flows in various geometries.     

Accurate analytical perturbation approach for large amplitude vibration of functionally graded beams

The present work derives the accurate analytical solutions for large amplitude vibration of thin functionally graded beams. In accordance with the Euler–Bernoulli beam theory and the von Kármán type geometric non-linearity, the second-order ordinary differential equation having odd and even non-linearities can be formulated through Hamilton's principle and Galerkin's procedure. This ordinary differential equation governs the non-linear vibration of functionally graded beams with different boundary constraints. Building on the original non-linear equation, two new non-linear equations with odd non-linearity are to be constructed. Employing a generalised Senator–Bapat perturbation technique as an ingenious tool, two newly formulated non-linear equations can be solved analytically. By selecting the appropriate piecewise approximate solutions from such two new non-linear equations, the analytical approximate solutions of the original non-linear problem are established. The present solutions are directly compared to the exact solutions and the available results in the open literature. Besides, some examples are selected to confirm the accuracy and correctness of the current approach. The effects of boundary conditions and vibration amplitudes on the non-linear frequencies are also discussed.     

Approximate dynamic boundary conditions for a thin piezoelectric layer

Approximate dynamic boundary conditions of different orders are derived for the case of a thin piezoelectric coating layer bonded to an elastic material. The approximate boundary conditions are derived using series expansions of the elastic displacements and the electric potential in the thickness coordinate of the layer. All the expansion functions are then eliminated with the aid of the equations of motion and boundary/interface conditions of the layer. This results in boundary conditions on the elastic material that may be truncated to different orders in the thickness of the layer to obtain approximate boundary conditions. The approximate boundary conditions may be used as a replacement for the piezoelectric layer and thus simplify the analysis significantly. Numerical examples show that the approximate boundary conditions give good results for low frequencies and/or thin piezoelectric layers.     

Transient responses of a magneto-electro-elastic hollow sphere for fully coupled spherically symmetric problem

By introducing a dependent variable and a special function satisfying the inhomogeneous mechanical boundary conditions, the governing equation for a new variable with homogeneous mechanical boundary conditions is derived. Then by means of the separation of variables technique and the electric and magnetic boundary conditions, the dynamic problem of a magneto-electro-elastic hollow sphere under spherically symmetric deformation is transformed to two Volterra integral equations of the second kind about two functions of time. Cubic Hermite polynomials are adopted to approximate the two undetermined functions at each time subinterval and the recursive formula is obtained to solve the integral equations successfully. The transient responses of displacements, stresses, electric and magnetic potentials are completely determined at the end. Numerical results are presented.     

Influence of magnetic field and thermal radiation by natural convection past vertical cone subjected to variable surface heat flux

A numerical study is performed to examine the heat transfer characteristics of natural convection past a vertical cone under the combined effects of magnetic field and thermal radiation.The surface of the cone is subjected to a variable surface heat flux.The fluid considered is a gray,absorbing-emitting radiation but a non-scattering medium.With approximate transformations,the boundary layer equations governing the flow are reduced to non-dimensional equations valid in the free convection regime.The dimensionless governing equations are solved by an implicit finite difference method of Crank-Nicolson type which is fast convergent,accurate,and unconditionally stable.Numerical results are obtained and presented for velocity,temperature,local and average wall shear stress,and local and average Nusselt number in air and water.The present results are compared with the previous published work and are found to be in excellent agreement.     

轴向功能梯度变截面梁的自由振动研究

摘 要：本文引入一种新的、简单易行的近似方法,求解轴向非均匀变截面梁的自由振动固有频率。将位移展开成切比雪夫多项式,从而变系数控制微分方程转化为含未知系数的齐次线性方程组。利用非零解的存在条件,进而得到含固有频率的特征方程。通过和特定梯度下已有的精确解进行比较,验证了该方法的精度和有效性,并分析了梯度参数、支承条件等对固有频率的影响。     

Self-similar thermal boundary layers on plane and axisymmetric bodies

This paper proposes an approximate solution procedure for the prediction of the forced convection heat transfer through self-similar laminar boundary layers. The differential equations governing the viscous and thermal boundary layers have been reduced to a pair of algebraic equations for the boundary layer shape factor and the boundary layer thickness ratio. The local Nusselt number predicted under various pressure gradients turns out to be in excellent agreement with that of the exact solution over a wide range of the Prandtl number.     

An inverse problem for a functionally graded elliptical plate with large deflection and slightly disturbed boundary

This paper deals with the inverse problem of a functionally graded material (FGM) elliptical plate with large deflection and disturbed boundary under uniform load. The properties of functionally graded material are assumed to vary continuously through the thickness of the plate, and obey a simple power law expression based on the volume fraction of the constituents. Based on the classical nonlinear von Karman plate theory, the governing equations of a thin plate with large deflection were derived. In order to solve this non-classical problem, a perturbation technique was employed on displacement terms in conjunction with Taylor series expansion of the disturbed boundary conditions. The displacements of in-plane and transverse are obtained in a non-dimensional series expansion form with respect to center deflection of the plate. The approximate solutions of displacements are solved for the first three terms, and the corresponding internal stresses can also be obtained.     

Nonlinear vibration of a two-mass system with nonlinear stiffnesses

An analytical approach is developed for nonlinear free vibration of a conservative, two-degree-of-freedom mass–spring system having linear and nonlinear stiffnesses. The main contribution of the proposed approach is twofold. First, it introduces the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation and, more significantly, the treatment a nonlinear differential system by linearization coupled with Newton’s method and harmonic balance method. New and accurate higher-order analytical approximate solutions for the nonlinear system are established. After solving the nonlinear differential equation, the displacement of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, this higher-order Newton–harmonic balance (NHB) method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation unlike the classical harmonic balance method which results in complicated algebraic equations requiring further numerical analysis. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and the classical harmonic balance methods. Two examples of nonlinear two-degree-of-freedom mass–spring system are analyzed and verified with published result, exact solutions and numerical integration data.     

Concurrent topology optimization for minimization of total mass considering load-carrying capabilities and thermal insulation simultaneously

The present work introduces a novel concurrent optimization formulation to meet the requirements of lightweight design and various constraints simultaneously. Nodal displacement of macrostructure and effective thermal conductivity of microstructure are regarded as the constraint functions, which means taking into account both the load-carrying capabilities and the thermal insulation properties. The effective properties of porous material derived from numerical homogenization are used for macrostructural analysis. Meanwhile, displacement vectors of macrostructures from original and adjoint load cases are used for sensitivity analysis of the microstructure. Design variables in the form of reciprocal functions of relative densities are introduced and used for linearization of the constraint function. The objective function of total mass is approximately expressed by the second order Taylor series expansion. Then, the proposed concurrent optimization problem is solved using a sequential quadratic programming algorithm, by splitting into a series of sub-problems in the form of the quadratic program. Finally, several numerical examples are presented to validate the effectiveness of the proposed optimization method. The various effects including initial designs, prescribed limits of nodal displacement, and effective thermal conductivity on optimized designs are also investigated. An amount of optimized macrostructures and their corresponding microstructures are achieved.     

Axisymmetric bending of functionally graded circular magneto-electro-elastic plates

Axisymmetric bending of functionally graded circular magneto-electro-elastic plates of transversely isotropic materials is analyzed based on linear three-dimensional theory of elasticity coupled with magnetic and electric fields. The transverse loads are expanded in Fourier-Bessel series and therefore can be arbitrarily distributed along the radial direction. The radial distributions of the displacements are assumed in combination of Fourier-Bessel series and polynomials as well as the electric potential and magnetic potential. If the material properties obey the exponential law along the thickness of the plate, two three-dimensional exact solutions for two unusual boundary conditions can be derived since they satisfy the governing equations and specified boundary conditions point by point. For simply supported or clamped boundary, the obtained solutions satisfy the governing equations exactly and the boundary conditions approximately. A layer wise model is also introduced to treat with the plates whose material property components vary independently and arbitrarily along the thickness of the plates. The numerical results are finally tabulated and plotted to demonstrate the presented method and agree well with those from finite element methods.     

A novel statistical linearization solution for randomly excited coupled bending-torsional beams resting on non-linear supports

A novel statistical linearization technique is developed for computing stationary response statistics of randomly excited coupled bending-torsional beams resting on non-linear elastic supports. The key point of the proposed technique consists in representing the non-linear coupled response in terms of constrained linear modes. The resulting set of non-linear equations governing the modal amplitudes is then replaced by an equivalent linear one via a classical statistical error minimization procedure, which provides algebraic non-linear equations for the second-order statistics of the beam response, readily solved by a simple iterative scheme. Data from Monte Carlo simulations, generated by a pertinent boundary integral method in conjunction with a Newmark numerical integration scheme, are used as benchmark solutions to check accuracy and reliability of the proposed statistical linearization technique.