Well-posedness of the problem of fiber suspension flows

The well-posedness of the equations governing the flow of fiber suspensions is studied. The fluid is assumed to be Newtonian and incompressible, and the presence of fibers is accounted for through the use of second- and fourth-order orientation tensors, which model the effects of the orientation of fibers in an averaged sense. The fourth-order orientation tensor is expressed in terms of the second-order tensor through various closure relations. It is shown that the linear closure relation leads to anomalous behavior, in that the rest state of the fluid is unstable, in the sense of Liapounov, for certain ranges of the fiber particle number. No such anomalies arise in the case of quadratic and hybrid closure relations. For the quadratic closure relation, it is shown that a unique solution exists locally in time for small data.     

An efficient three-dimensional semi-implicit finite element scheme for simulation of free surface flows

An efficient semi-implicit finite element model is proposed for the simulation of three-dimensional flows in stratified seas. The body of water is divided into a number of layers and the two horizontal momentum equations for each layer of water are first integrated vertically. Nine-node Lagrangian quadratic isoparametric elements are employed for spatial discretization in the horizontal domain. The time derivatives are approximated using a second-order-accurate semi-implicit time-stepping scheme. The distinguishing feature of the proposed numerical scheme is that only nodal values on the same vertical line are coupled. Two test cases for which analytic solutions are available are employed to test the proposed scheme. The test results show that the scheme is efficient and stable. A numerical experiment is also included to compare the proposed scheme with a finite difference scheme.     

Stochastic non-linear flutter of a panel subjected to random in-plane forces

—An analysis of non-linear flutter of a simply-supported panel exposed to supersonic gas flow and random in-plane forces is presented for two- and three-mode interactions. A first order quasi-steady state aerodynamic piston theory is used to model the aerodynamic loading. The Fokker-Planck equation is used to derive a general moment equation for two- and three-mode interactions. For stability analysis the moment equation is consistent and the mean square stability boundaries of the equilibrium are obtained in terms of the system parameters. The stability boundaries reveal common features to those predicted by the deterministic theory of panel nutter. For the non-linear response the moment equation is found inconsistent and a cumulant-neglect closure is used by setting cumulants of fifth and sixth orders to zero. This first order non-Gaussian closure is carried out to solve for the response statistics in terms of the air-to-plate mass ratio, aerodynamic pressure, modal damping, and in-plane random force spectral density. It is found that the non-Gaussian solution yields higher levels for the response statistics than those obtained by the Gaussian solution. The inclusion of more modes results in a reduction of the response levels and expands the stability region.     

A multi‐moment transport model on cubed‐sphere grid

A conservative, single‐cell‐based semi‐Lagrangian transport model is proposed in this paper. Using multi‐moment concept, an additional moment, i.e. volume‐integrated average (VIA), is treated as the model variable besides the point value (PV) updated in the traditional semi‐Lagrangian schemes. A quadratic interpolation function is constructed based on local degrees of freedom defined within each single cell. The PV moment is advanced by the semi‐Lagrangian formulation, whereas the VIA moment is updated by a finite volume formulation to rigorously ensure the numerical conservation. The numerical fluxes are computed from the PV moments defined along the boundary edges of the control volume. The scheme is extended to the spherical geometry through the application of the cubed‐sphere grid that eliminates the polar singularity in the conventional longitude/latitude coordinates by using the quasi‐uniform grid spacing covering the whole sphere. The single‐cell‐based scheme is well suited for the treatment of the connections between different patches. A simple quasi‐monotone limiter to the PV moment is applied to suppress non‐physical oscillations. The proposed scheme has been validated via representative benchmark tests and the performance is competitive to other existing transport schemes. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite Element Analysis of Dynamically Loaded Homogeneous Viscous Beams

The paper is concerned with impulsively loaded beams in which the material is treated as homogeneous viscous as an approximation of a rigid-viscoplastic constitutive relation. As opposed to the standard displacement method finite element formulation, where interpolation functions describing the velocity field across elements is given, a mixed formulation is used in which nodal velocities and nodal moments are carried as parameters. At each instant the accelerations (by the Tamuzh principle) and the rates of change of moment (by a virtual velocities formulation) are found, and velocities and moments are integrated forward independently. The properties of the mode solution are also introduced, and the forward integration is carried through only for the difference between the mode solution and the actual solution. This leads to a very efficient scheme for the numerical solution of a cantilever beam problem shown as an illustration.     

RCM‐TVD hybrid scheme for hyperbolic conservation laws

We describe a hybrid method for the solution of hyperbolic conservation laws. A third‐order total variation diminishing (TVD) finite difference scheme is conjugated with a random choice method (RCM) in a grid‐based adaptive way. An efficient multi‐resolution technique is used to detect the high gradient regions of the numerical solution in order to capture the shock with RCM while the smooth regions are computed with the more efficient TVD scheme. The hybrid scheme captures correctly the discontinuities of the solution and saves CPU time. Numerical experiments with one‐ and two‐dimensional problems are presented. Copyright © 2007 John Wiley & Sons, Ltd.     

多自由度强非线性耦合参激系统随机响应计算方法

本文采用非高斯矩闭合方法计算了一个含参数激振的二自由度强非线性耦合系统受随机激励的响应方差。用逐次消除法消除方程中的惯性耦合后推导出直到6阶的矩方程,采用中心累积量截断技术计算了白噪音激励下的响应方差,并与时域直接积分和数字模拟方法的结果进行比较,取得了一些有意义的结论。     

A stable algorithm for bed friction in three-dimensional shallow sea modal models

The mathematical formulation of a three-dimensional shallow sea model using a modal expansion in the veitical is briefly described. The importance of the time discretization of the vertical diffusion term and bottom friction term is discussed in some detail. Both stability theory and numerical calculations show the importance of time centring or evaluating the modal form of the viscosity term at the higher time step in order to develop a numerically efficient algorithm. Similar analysis and calculations show that in shallow water it is essential to time centre or evaluate bottom friction at the higher time step. In the case of linear bottom friction it is shown that this condition can be readily accomplished. However, using a quadratic friction formulation (a more physically realistic form), this cannot be readily achieved. A new algorithm is presented whereby a stable solution can be obtained even in shallow water using quadratic bottom friction.     

Orientation tensors in simple flows of dilute suspensions of non-Brownian rigid ellipsoids,comparison of analytical and approximate solutions

General analytical solutions are obtained for the planar orientation structure of rigid ellipsoid of revolutions subjected to an arbitrary homogeneous flow in a Newtonian fluid. Both finite and infinite aspect ratio particles are considered. The orientation structure is described in terms of two-dimensional, time-dependent tensors that are commonly employed in constitutive equations for anisotropic fluids such as fiber suspensions. The effect of particle aspect ratio on the evolution of orientation structure is studied in simple shear and planar elongational flows. With the availability of analytical solutions, accuracies of quadratic closure approximations used for nonhomogeneous flows are analyzed, avoiding numerical integration of orientation distribution function. In general, fourth-order orientation evolution equations with sixth-order quadratic closure approximations yield more accurate representations compared to the commonly used second-order evolution equations with fourth-order quadratic closure approximations. However, quadratic closure approximations of any order are found to give correct maximum orientation angle (i.e., preferred direction) results for all particle aspect ratios and flow cases.     

On fallacies in the decision between the Caputo and Riemann–Liouville fractional derivatives for the analysis of the dynamic response of a nonlinear viscoelastic oscillator

Recently Dal [Dal, F., 2011. Multiple time scale solution of an equation with quadratic and cubic nonlinearities having fractional-order derivative. Mathematical and Computational Applications 16 (1), 301–308] presented ‘a new analytical scheme’ to calculate the dynamic response of a fractionally damped nonlinear oscillator possessing both quadratic and cubic nonlinearities via the method of multiple time scales. It has been claimed that damping features are modeled via the Caputo fractional derivative. In the present paper, it is shown that both the scheme and the object of investigation are not new, and moreover, the above mentioned author's statement is inconsistent, since under the assumptions made in the paper under consideration these two fractional-order derivatives coincide. Besides, the utilized procedure was inconsequential. It has been proved that the investigation of the dynamic response of a nonlinear viscoelastic oscillator presents the case that, with some minimal restrictions, the Riemann–Liouville and Caputo definitions produce completely equivalent mathematical models of the nonlinear viscoelastic phenomenon.     

Presumed Mapping Functions for Eulerian Modelling of Turbulent Mixing

The mapping closure of Chen et al. [Phys. Rev. Lett., 63, 1989] is a transported probability density function (PDF) method that has proven very efficient for modelling of turbulent mixing in homogeneous turbulence. By utilizing a Gaussian reference field, the solution to the mapping function (in homogeneous turbulence) can be found analytically for a range of initial conditions common for turbulent combustion applications, e.g. for binary or trinary mixing. The purpose of this paper is to investigate the possibility of making this solution a presumed mapping function (PMF) for inhomogeneous flows. The PMF in turn will imply a presumed mixture fraction PDF that can be used for a wide range of models in turbulent combustion, e.g. flamelet models, the conditional moment closure (CMC) or large eddy simulations. The true novelty of the paper, though, is in the derivation of highly efficient, closed algebraic expressions for several existing models of conditional statistics, e.g. for the conditional scalar dissipation/diffusion rate or the conditional mean velocity. The closed form expressions nearly eliminates the overhead computational cost that usually is associated with nonlinear models for conditional statistics. In this respect it is argued that the PMF is particularly well suited for CMC that relies heavily on manipulations of the PDF for consistency. The accuracy of the PMF approach is shown with comparison to DNS of a single scalar mixing layer to be better than for the β-PDF. Not only in the shape of the PDF itself, but also for all conditional statistics models computed from the PDF.     

泊松白噪声激励下强非线性系统的半解析瞬态解

自然界与工程中都普遍存在着随机扰动,且大多数呈现出固有的非高斯性质,若采用高斯激励建模可能会导致巨大的误差.泊松白噪声作为一种典型且重要的非高斯激励模型,已引起了广泛的关注.目前,泊松白噪声激励下系统的动态特性分析主要集中于稳态响应的研究,而针对瞬态响应的求解难度仍较大,需进一步发展.本文引入径向基神经网络,提出了一种泊松白噪声激励下单自由度强非线性系统瞬态响应预测的高效半解析方法.首先将广义Fokker-Plank-Kolmogorov (FPK)方程的瞬态解表示为一组含时变待定权值系数的高斯径向基神经网络;然后采用有限差分法离散时间导数项,并结合随机取样技术构造含时间递推式的损失函数;最后通过拉格朗日乘子法使得损失函数最小化获得时变最优权值系数.作为算例,探究了两个经典强非线性系统,并采用蒙特卡罗模拟方法对解析结果加以验证.结果表明:本文方法所获得的瞬时概率密度函数与蒙特卡罗模拟数据吻合地较好,并且算法具备较高的计算效率.在系统响应的整个演化过程中,本文所提方法能够非常有效地捕捉到系统响应在各个时刻下的复杂非线性特征.此外,本文方法所获得的高精度半解析瞬态解,不仅可作为基准解检验其...     

A direct probabilistic approach to solve state equations for nonlinear systems under random excitation

In this paper, a direct probabilistic approach (DPA) is presented to formulate and solve moment equations for nonlinear systems excited by environmental loads that can be either a stationary or nonstationary random process. The proposed method has the advantage of obtaining the response’s moments directly from the initial conditions and statistical characteristics of the corresponding external exci-tations. First, the response’s moment equations are directly derived based on a DPA, which is completely independent of the It?/filtering approach since no specific assumptions regarding the correlation structure of excitation are made. By solving them under Gaussian closure, the response’s moments can be obtained. Subsequently, a multiscale algo-rithm for the numerical solution of moment equations is exploited to improve computational efficiency and avoid much wall-clock time. Finally, a comparison of the results with Monte Carlo (MC) simulation gives good agreement. Furthermore, the advantage of the multiscale algorithm in terms of efficiency is also demonstrated by an engineering example.     

窄带随机噪声作用下单自由度非线性碰撞系统的响应

研究了单自由度非线性单边碰撞系统在窄带随机噪声激励下的次共振响应问题。用Zhuravlev变换将碰撞系统转化为速度连续的非碰撞系统,然后用随机平均法得到了关于慢变量的随机微分方程。在没有随机扰动情形,得到了系统响应幅值满足的代数方程;在有随机扰动的情形下,给出了系统响应稳态矩计算的迭代公式。讨论了系统阻尼项、非线性项、随机扰动项和碰撞恢复系数等参数对于系统响应的影响。理论计算和数值模拟表明,系统响应幅值将在激励频率接近于次共振频率时达到最大。而当激励频率逐渐偏离次共振频率时,系统响应迅速衰减。     

Implicit time accurate simulation of unsteady flow

Implicit time integration was studied in the context of unsteady shock‐boundary layer interaction flow. With an explicit second‐order Runge–Kutta scheme, a reference solution to compare with the implicit second‐order Crank–Nicolson scheme was determined. The time step in the explicit scheme is restricted by both temporal accuracy as well as stability requirements, whereas in the A‐stable implicit scheme, the time step has to obey temporal resolution requirements and numerical convergence conditions. The non‐linear discrete equations for each time step are solved iteratively by adding a pseudo‐time derivative. The quasi‐Newton approach is adopted and the linear systems that arise are approximately solved with a symmetric block Gauss–Seidel solver. As a guiding principle for properly setting numerical time integration parameters that yield an efficient time accurate capturing of the solution, the global error caused by the temporal integration is compared with the error resulting from the spatial discretization. Focus is on the sensitivity of properties of the solution in relation to the time step. Numerical simulations show that the time step needed for acceptable accuracy can be considerably larger than the explicit stability time step; typical ratios range from 20 to 80. At large time steps, convergence problems that are closely related to a highly complex structure of the basins of attraction of the iterative method may occur. Copyright © 2001 John Wiley & Sons, Ltd.     

Some considerations for high‐order ‘incremental remap’‐based transport schemes: edges,reconstructions, and area integration

The problem of two‐dimensional tracer advection on the sphere is extremely important in modeling of geophysical fluids and has been tackled using a variety of approaches. A class of popular approaches for tracer advection include ‘incremental remap’ or cell‐integrated semi‐Lagrangian‐type schemes. These schemes achieve high‐order accuracy without the need for multistage integration in time, are capable of large time steps, and tend to be more efficient than other high‐order transport schemes when applied to a large number of tracers over a single velocity field. In this paper, the simplified flux‐form implementation of the Conservative Semi‐LAgrangian Multi‐tracer scheme (CSLAM) is reformulated using quadratic curves to approximate the upstream flux volumes and Gaussian quadrature for integrating the edge flux. The high‐order treatment of edge fluxes is motivated because of poor accuracy of the CSLAM scheme in the presence of strong nonlinear shear, such as one might observe in the midlatitudes near an atmospheric jet. Without the quadratic treatment of upstream edges, we observe at most second‐order accuracy under convergence of grid resolution, which is returned to third‐order accuracy under the improved treatment. A shallow‐water barotropic instability also reveals clear evidence of grid imprinting without the quadratic correction. Consequently, these tests reveal a problem that might arise in tracer transport near nonlinearly sheared regions of the real atmosphere, particularly near cubed‐sphere panel edges. Although CSLAM is used as the foundation for this analysis, the conclusions of this paper are applicable to the general class of incremental remap schemes. Copyright © 2012 John Wiley & Sons, Ltd.     

Upstream monotonic interpolation for scalar transport with application to complex turbulent flows

A new monotonic scheme for the approximation of steady scalar transport is formulated and implemented within a collocated finite-volume/pressure-correction algorithm for general turbulent flows in complex geometries. The scheme is essentially a monotonic implementation of the quadratic QUICK interpolation and uses a continuous and compact limiter to secure monotonicity. The principal purpose is to allow an accurate and fully bounded, hence stable, approximation of turbulence convection in the context of two-equation eddy viscosity and Reynolds stress transport modelling of two- and three-dimensional flows, both subsonic and transonic. Among other benefits, this capability permits an assessment to be made of the adequacy of approximating turbulence convection with first-order upwind schemes in conjunction with higher-order formulations for mean-flow properties—a widespread practice. The performance characteristics of the bounded scheme are illustrated by reference to computations for scalar transport, for a transonic flow in a Laval nozzle, for one separated laminar flow and for two separated turbulent flows computed with a non-linear RNG model and full Reynolds stress closure.     

Employment of the second‐moment turbulence closure on arbitrary unstructured grids

The paper presents a finite‐volume calculation procedure using a second‐moment turbulence closure. The proposed method is based on a collocated variable arrangement and especially adopted for unstructured grids consisting of ‘polyhedral’ calculation volumes. An inclusion of 23k in the pressure is analysed and the impact of such an approach on the employment of the constant static pressure boundary is addressed. It is shown that this approach allows a removal of a standard but cumbersome velocity–pressure –Reynolds stress coupling procedure known as an extension of Rhie‐Chow method (AIAA J. 1983; 21 : 1525–1532) for the Reynolds stresses. A novel wall treatment for the Reynolds‐stress equations and ‘polyhedral’ calculation volumes is presented. Important issues related to treatments of diffusion terms in momentum and Reynolds‐stress equations are also discussed and a new approach is proposed. Special interpolation practices implemented in a deferred‐correction fashion and related to all equations, are explained in detail. Computational results are compared with available experimental data for four very different applications: the flow in a two‐dimensional 180o turned U‐bend, the vortex shedding flow around a square cylinder, the flow around Ahmed Body and in‐cylinder engine flow. Additionally, the performance of the methodology is assessed by applying it to different computational grids. For all test cases, predictions with the second‐moment closure are compared to those of the k–εmodel. The second‐moment turbulence closure always achieves closer agreement with the measurements. A moderate increase in computing time is required for the calculations with the second‐moment closure. Copyright © 2004 John Wiley & Sons, Ltd.     

Improved quadratic isogeometric element simulation of one-dimensional elastic wave propagation with central difference method

Two improved isogeometric quadratic elements and the central difference scheme are used to formulate the solution procedures of transient wave propagation problems. In the proposed procedures, the lumped matrices corresponding to the isogeometric elements are obtained. The stability conditions of the solution procedures are also acquired. The dispersion analysis is conducted to obtain the optimal Courant-Friedrichs-Lewy (CFL) number or time-step sizes corresponding to the spatial isogeometric elements. The dispersion analysis shows that the isogeometric quadratic element of the fourth-order dispersion error (called the isogeometric analysis (IGA)-f quadratic element) provides far more desirable numerical dissipation/dispersion than the element of the second-order dispersion error (called the IGA-s quadratic element) when appropriate time-step sizes are selected. The numerical simulations of one-dimensional (1D) transient wave propagation problems demonstrate the effectiveness of the proposed solution procedures.