On Solutions of Some Non-Linear Differential Equations Arising in Newtonian and Non-Newtonian Fluids

Some steady as well as unsteady solutions of the equations of motion for an incompressible Newtonian and non-Newtonian (second-grade) fluids are obtained by applying different methods including the Lie symmetry group method. The flows considered are axially symmetric with the swirling motion, and the governing equations for second-grade fluid flow have been modeled. Expressions for streamlines, velocity and vorticity components are constructed explicitly in each case. Exact analytical solutions in second-grade fluid are obtained and compared with the corresponding viscous solutions.     

On the Investigation of Some Parameter Identification and Experimental Modal Filtering Issues for Nonlinear Reduced Order Models

This paper discusses modal filtering of experimental data and the corresponding identification of linear and nonlinear parameters in reduced order space. Specifically, several experimental configurations will be discussed in order to provide insight into such identification issues as spatial discretization, observability, and the linear independence of the assumed filter or basis. The two experiments considered herein represent different measurement configurations of the same clamped–clamped beam. First, asymmetric inertial loading via asymmetric sensor location was considered, while the second scenario presents a symmetric sensor configuration. Several important conclusions can be drawn from the two experimental scenarios. First, by asymmetrically loading the beam, a corresponding asymmetric beam mode was excited yet not observable. In the second scenario, the symmetric distribution of sensors minimized the impact of the respective asymmetric mode. The resulting spatial information allowed for the proper filtering of the remnants of the asymmetric mode. Nonlinear parameters in modal space as well as the underlying linear parameters were successfully identified simultaneously in both experimental scenarios, although the usefulness of the asymmetrically loaded beam was limited. Finally, successful comparisons were made between the identified reduced order model and experimental response at the beam quarter point using the symmetric case and the beam midpoint using both experimental scenarios.     

Analysis of the Applicability of Macroscopic Models of the Shock Wave Structure in a Binary Mixture of Inert Gases

The results of calculating the shock wave structure in Ne–Ar, He–Ar, He–Ne, and He–Xe mixtures by means of the relaxation method on the basis of the system of Navier-Stokes equations and complete and modified systems of Burnett equations are compared with the results of direct statistical simulation (Monte-Carlo method). The domain of applicability of these systems of equations for calculating gas dynamic variable profiles is analyzed as a function of both the molecular mass ratio and the initialconcentrations.     

Lie Symmetry Analysis of Differential Equations in Finance

Lie group theory is applied to differential equations occurring as mathematical models in financial problems. We begin with the complete symmetry analysis of the one-dimensional Black–Scholes model and show that this equation is included in Sophus Lie's classification of linear second-order partial differential equations with two independent variables. Consequently, the Black–Scholes transformation of this model into the heat transfer equation follows directly from Lie's equivalence transformation formulas. Then we carry out the classification of the two-dimensional Jacobs–Jones model equations according to their symmetry groups. The classification provides a theoretical background for constructing exact (invariant) solutions, examples of which are presented.     

Spatial Decay Estimates for a Coupled System of Second-Order Quasilinear Partial Differential Equations Arising in Thermoelastic Finite Anti-plane Shear

The spatial decay behavior of solutions of a coupled system of second-order quasilinear partial differential equations, in divergence form, defined on a two-dimensional semi-infinite strip, is investigated. Such equations arise in the theory of anti-plane shear deformations for isotropic nonlinearly thermoelastic solids. Differential inequality techniques are employed to obtain exponential decay estimates. The results are illustrated by several examples. The results are relevant to Saint-Venant principles for nonlinear thermoelasticity as well as to theorems of Phragmen-Lindelof type. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Connection Between Isometries and Symmetries of Geodesic Equations of the Underlying Spaces

A connection between the symmetries of manifolds and differential equations is sought through the geodesic equations of maximally symmetric spaces, which have zero, constant positive or constant negative curvature. It is proved that for a space admitting so(n+1) or so(n,1) as the maximal isometry algebra, the symmetry of the geodesic equations of the space is given by so( or (where d ₂ is the two-dimensional dilation algebra), while for those admitting (where represents semidirect sum) the algebra is sl(n+2). A corresponding result holds on replacing so(n) by so(p,q) with p+q = n. It is conjectured that if the isometry algebra of any underlying space of non-zero curvature is h, then the Lie symmetry algebra of the geodesic equations is given by , provided that there is no cross-section of zero curvature at the point under consideration. If there is a flat subspace of dimension m, then the symmetry group becomes ).     

The Conjugacy of Stochastic and Random Differential Equations and the Existence of Global Attractors

We consider stochastic differential equations in d-dimensional Euclidean space driven by an m-dimensional Wiener process, determined by the drift vector field f₀ and the diffusion vector fields f₁,...,f_m, and investigate the existence of global random attractors for the associated flows . For this purpose is decomposed into a stationary diffeomorphism given by the stochastic differential equation on the space of smooth flows on R^d driven by m independent stationary Ornstein Uhlenbeck processes z¹,...,z^m and the vector fields f₁,...,f_m, and a flow generated by the nonautonomous ordinary differential equation given by the vector field (_t/x)^–1[f₀(_t)+ _i=1 ¹ f_i(_t)z _t ⁱ ]. In this setting, attractors of are canonically related with attractors of . For , the problem of existence of attractors is then considered as a perturbation problem. Conditions on the vector fields are derived under which a Lyapunov function for the deterministic differential equation determined by the vector field f₀ is still a Lyapunov function for , yielding an attractor this way. The criterion is finally tested in various prominent examples.     

Comparison of Kinetic and Hybrid-Equilibrium Models Simulating Colloid-Facilitated Contaminant Transport in Porous Media

The presence of colloidal particles in groundwater can enhance contaminant transport by reducing retardation effects and carrying them to distances further than predicted by a conventional advective/dispersive equation with normal retardation values. When colloids exist in porous media and affect contaminant migration, the system can best be simulated as a three-phase medium. Mechanisms of mass transfer from one phase to another by colloids and contaminants can be kinetic or equilibrium-based, depending on the sorption–desorption reaction rate between the aqueous and solid phases. When the rate of sorption between the water phase and the solid phase(s) is not much greater than the rate of change in contaminant concentration in the water phase, kinetic sorption models may better describe the phenomenon. In some cases of modeling one or more mass transfer processes, a useful simplification may be to introduce the local equilibrium assumption. In this study, the local equilibrium assumption for sorption processes on colloidal surfaces (hybrid equilibrium model) was compared with kinetic-based models. Sensitivity analyses were conducted to deduce the effect of major parameters on contaminant transport. The results obtained from the hybrid equilibrium model in predicting the transport of colloid-facilitated groundwater contaminant are very similar to those of the kinetic model, when the point of interest is not at contaminant and colloid source vicinities and the time of interest is sufficiently long for imposed sources.     

钢筋混凝土剪力墙单元非线性分析模型及其应用

本文简要评述了几种常见的钢筋混凝土剪力墙非线性宏观单元模型后，着重对多垂直杆剪力墙非线性单元模型的几个重要问题如剪力变形的考虑方法，单元刚度矩阵的形式，垂直拉压杆及剪切弹簧的恢复力模型等进行了探讨与改进，最后给出了一个算例，并与试验结果比较，表明非线性宏观墙单元模型具有较好的计算精度。     

Analysis of one-dimensional horizontal two-phase flow in geothermal reservoirs

In the absence of capillarity the single-component two-phase porous medium equations have the structure of a nonlinear parabolic pressure (equivalently, temperature) diffusion equation, with derivative coupling to a nonlinear hyperbolic saturation wave equation. The mixed parabolic-hyperbolic system is capable of substaining saturation shock waves. The Rankine-Hugoniot equations show that the volume flux is continuous across such a shock. In this paper we focus on the horizontal one-dimensional flow of water and steam through a block of porous material within a geothermal reservoir. Starting from a state of steady flow we study the reaction of the system to simple changes in boundary conditions. Exact results are obtainable only numerically, but in some cases analytic approximations can be derived. When pressure diffusion occurs much faster than saturation convection, the numerical results can be described satisfactorily in terms of either saturation expansion fans, or isolated saturation shocks. At early times, pressure and saturation profiles are functionally related. At intermediate times, boundary effects become apparent. At late times, saturation convection dominates and eventually a steady-state is established. When both pressure diffusion and saturation convection occur on the same timescale, initial simple shock profiles evolve into multiple shocks, for which no theory is currently available. Finally, a parameter-free system of equations is obtained which satisfactorily represents a particular case of the exact equations.     

Performances of LES and RANS Models for Simulation of Complex Flows in a Coaxial Jet Mixer

Flow and mixing processes in a classical coaxial jet mixer have been investigated numerically. Calculations have been performed using three Large Eddy Simulation models and three unsteady RANS models. The time averaged mixture fraction and axial velocity, their rms values and energy spectra are compared with LIF and LDA measurements for both j- and r-modes of the jet mixer flow. A special attention is paid to the ability of different models to reproduce unsteady effects. The analysis demonstrates the superiority of the LES method with the dynamic mixed SGS model (DMM) with respect to other RANS and LES models.     

Some existence theorems of common and coincidence solutions for a class of systems of functional equations arising in dynamic programming

Some existence theorems of common and coincidence solutions for a class of moregeneral systems of functional equations arising in dynamic programming are shown.Theresults presented in this paper not only contain the corresponding results of[6,7]as specialcases,but also give an existence theorem of solutions for a class of functional equationssuggested by Wang recently.     

Assessment of Turbulence Models for Predicting the Interaction Region in a Wall Jet by Reference to LES Solution and Budgets

Highly-resolved LES and experimental data for a plane wall jet are used to study the characteristics of turbulence-closure proposals, mainly within the framework of second-moment-transport modelling. The study is motivated by the observed importance of diffusive Reynolds-stress transport in the interaction region between the outer shear layer and the near-wall layer of the wall jet, which gives the near-wall flow characteristics that are very different from those of a conventional boundary layer. Comparisons are presented for mean-flow quantities, second moments and budgets. Also included are a priori studies of approximations for the pressure-velocity interaction, pressure-fluctuation-driven transport and turbulent transport of the Reynolds stresses by triple correlations, the last observed to contribute significantly to the stress budgets. The study reveals major defects in the closure approximations for the pressure-velocity interaction terms, especially in the near-wall region. These defects result in a poor representation by the particular second-moment closures investigated of even the integral and mean-flow characteristics of the wall jet.     

Outflow boundary conditions for one-dimensional finite element modeling of blood flow and pressure waves in arteries

Flow and pressure waves, originating due to the contraction of the heart, propagate along the deformable vessels and reflect due to tapering, branching, and other discontinuities. The size and complexity of the cardiovascular system necessitate a “multiscale” approach, with “upstream” regions of interest (large arteries) coupled to reduced-order models of “downstream” vessels. Previous efforts to couple upstream and downstream domains have included specifying resistance and impedance outflow boundary conditions for the nonlinear one-dimensional wave propagation equations and iterative coupling between three-dimensional and one-dimensional numerical methods. We have developed a new approach to solve the one-dimensional nonlinear equations of blood flow in elastic vessels utilizing a space-time finite element method with GLS-stabilization for the upstream domain, and a boundary term to couple to the downstream domain. The outflow boundary conditions are derived following an approach analogous to the Dirichlet-to-Neumann (DtN) method. In the downstream domain, we solve simplified zero/one-dimensional equations to derive relationships between pressure and flow accommodating periodic and transient phenomena with a consistent formulation for different boundary condition types. In this paper, we also present a new boundary condition that accommodates transient phenomena based on a Green’s function solution of the linear, damped wave equation in the downstream domain.     

结构分支型初始后屈曲中有限元平衡方程组的解法

根据结构分支型初始后屈曲中有约束的有限元平衡方程组的特点，推导了一种具有工程实用价值的变换方法，将其化为无约束的线性方程组，并给出了数值算例。     

Macroscopic modelling of transport phenomena in porous media. 2: Applications to mass,momentum and energy transport

In this second paper, the averaging rules presented in Part 1 are employed in order to develop a general macroscopic balance equation and particular equations for mass, mass of a component, momentum and energy, all of a phase in a porous medium domain. These balance equations involve averaged fluxes. Then macroscopic equations are developed for advective, dispersive and diffusive fluxes, all in terms of averaged state variables of the system. These are combined with the macroscopic balance equations to yield field equations that serve as the core of the mathematical models that describe the transport of extensive quantities in a porous medium domain. It is shown that the methodology of averaging leads to a better understanding of the effective stress concept employed in dealing with transport phenomena in deformable porous media.     

The method of dynamic mode decomposition in shallow water and a swirling flow problem

In this paper, a new vector‐filtering criterion for dynamic modes selection is proposed that is able to extract dynamically relevant flow features from dynamic mode decomposition of time‐resolved experimental or numerical data. We employ a novel modes selection criterion in parallel with the classic selection based on modes amplitudes, in order to analyze which of these procedures better highlight the coherent structures of the flow dynamics. Numerical tests are performed on two distinct problems. The efficiency of the proposed criterion is proved in retaining the most influential modes and reducing the size of the dynamic mode decomposition model. By applying the proposed filtering mode technique, the flow reconstruction error is shown to be significantly reduced. Copyright © 2016 John Wiley & Sons, Ltd.     

The study of quadruple integral equations involving inverse Mellin transforms and its application to a contact problem for a wedge shaped elastic solid in antiplane shear distribution

Closed form solution of quadruple integral equations involving inverse Mellin transforms has been obtained. The solution of quadruple integral equations is used in solving a two dimensional four-part mixed boundary value contact problem for an elastic wedge-shaped region as an application. Closed form expression for shear stress has been obtained. Finally, numerical results for shear stress are obtained and shown graphically.     

The numerical study of roll‐waves in inclined open channels and solitary wave run‐up

In this paper, we introduce a finite‐volume kinetic BGK scheme and its applications to the study of roll and solitary waves. The current scheme is based on the numerical solution of the gas‐kinetic Bhatnagar–Gross–Krook model in the flux evaluation across each cell interface. An intrinsic connection between the BGK model and time‐dependent, non‐linear, non‐homogeneous shallow‐water equations enables us to solve shallow‐water equations automatically with our kinetic scheme. The analytical solution, experimental measurements, and numerical calculations for problems associated with roll‐waves down an inclined open channel and solitary waves incident on a sloped beach are also presented. Copyright © 2005 John Wiley & Sons, Ltd.