Equilibria of axially moving beams in the supercritical regime

Equilibria of axially moving beams are computationally investigated in the supercritical transport speed ranges. In the supercritical regime, the pattern of equilibria consists of the straight configuration and of non-trivial solutions that bifurcate with transport speed. The governing equations of coupled planar is reduced to a partial-differential equation and an integro-partial-differential equation of transverse vibration. The numerical schemes are respectively presented for the governing equations and the corresponding static equilibrium equation of coupled planar and the two governing equations of transverse motion for non-trivial equilibrium solutions via the finite difference method and differential quadrature method under the simple support boundary. A steel beam is treated as example to demonstrate the non-trivial equilibrium solutions of three nonlinear equations. Numerical results indicate that the three models predict qualitatively the same tendencies of the equilibrium with the changing parameters and the integro-partial-differential equation yields results quantitatively closer to those of the coupled equations.     

Transient Heat and Moisture Flow Around Heat Source Buried in an Unsaturated Half Space

This article presents solutions for the transient heat and moisture transport due to both disk heat source and cylindrical heat source buried in an unsaturated half space. The solutions are presented in Hankel–Laplace transform domain and in dimensionless style. Coupled effect of thermally driven moisture transport is especially investigated because of its importance to alter the flow field in low-permeability medium. Parametric study has been performed to assess the effects of five independent dimensionless parameters on flow field. The stability and accuracy of the present solutions are demonstrated from the comparison between the results obtained from these solutions and those by using a well-established finite element code CODE_BRIGHT. Despite the simplified assumptions required in order to obtain analytical solutions in Hankel–Laplace transform domain, the results incorporate the main mechanisms involved in the coupled thermo-hydraulic (T-H) problem, and they may be eventually used for validation purposes.     

Analytical Solutions for Partitioned Diffusion in Laminates: I. Initial Value Problem with Steady Cauchy Conditions

Studies of the transport of contaminants and nutrients in industrial and environmental systems are complicated by the heterogeneous nature of the supporting porous or permeable media, and by the numerical problems associated with high Peclet number advection and sharp interface models. In order to provide independent theoretical checks of numerical transport theories, this set of papers presents analytical solutions to diffusive transport equations in simplified (one-dimensional) laminate systems subject to partitioning interactions. Here, in Part I, a standard separation of variables technique is used to develop analytical eigenfunction expansions of the concentration solution in an N-laminate system subject to steady Cauchy (third-type) nonhomogeneous boundary conditions. Both Cartesian and radial (axisymmetric) coordinate systems are considered. The solutions are developed for two different interface partitioning formulations, allowing the partitioning processes to be described by instantaneous equilibration mechanisms, or in terms of gradual equilibration mediated by mass transfer coefficients. Worked examples are presented and limitations of the approach discussed.     

Some Stratigraphic Control Problems

In this paper, we are interested in lithology diffusion models applied in the field of stratigraphic basin simulations for large-scale depositional transport processes of sediments. Such models describe erosion-sedimentation processes and take into account limited weathering via nonstandard unilateral problems. Various theoretical results, illustrations, and numerical solutions are presented for the monolithologic column case. A new conservation law involved in modeling is formulated, and mathematical tools for solving the problem are described.     

Delay Model for a Cycling Transport Through Porous Medium

A solute transport through a porous medium is examined provided that the fluid leaving the porous sample returns back in a continuous way. The porous medium is thus included into a closed hydrodynamic circuit. This cycling process is suggested as an experimental tool to determine porous medium parameters describing transport. In the present paper the mathematical theory of this method is developed. For the advective type of transport with solute retention and degradation in porous medium, the system of transport equations in a closed circuit is transformed to a delay differential equation. The exact analytical solution to this equation is obtained. The solute concentration manifests both the oscillatory and monotonous behaviors depending on system parameters. The number of oscillation splashes is shown to be always finite. The maximum/minimum points are determined as solutions of a polynomial equation whose degree depends on the unknown solution itself. The cyclic methods to determine porous medium parameters as porosity and retention rate are developed.     

Normal flow boundary conditions in 3D circulation models

The problem of enforcing normal transport conditions on 3D velocity fields is considered in the context of ‘wave equation’ finite element models. A procedure for strong enforcement of the transport constraint is given. The procedure is identical for Neumann (transport known) and Dirichlet (pressure known) problems, which are treated reversibly. All local mass and force balance relations are retained in the FEM system. A global mass conservation property is proven for the general 3D, discrete-time case. Examples demonstrate the quality of the solutions and the practicality of the approach. © 1997 John Wiley & Sons, Ltd.     

The Permeability of a Representative Carbonate Volume with a Large Vug

There are fundamental challenges in characterizing transport properties of a porous medium with large cavities, or vugs, when the characteristic size of the cavity is larger than or equal to 1 cm. Neither existing flow models nor common laboratory measurements are well suited to deal with such challenges. The present study determines the effective permeability of a representative carbonate volume with an arbitrary connected-through large vug. It also determines the minimum matrix permeability that impacts the flow in vuggy carbonates. The Stokes and Darcy models are used to capture the flows in the vug and in the matrix, respectively. At the interface, the flow models interact through slippage based on the Beavers–Joseph–Saffman model. The developed model is tested with analytical solutions available in the literature and using independent laboratory measurements. The results have major implications for characterizing the effective transport properties of carbonates, which constitute more than half of the world’s hydrocarbon reserves.     

Modelling of Advection-Dominated Transport in a Porous Column with a Time-Decaying Flow Field

This paper examines the problem of the advective-dispersive movement of a non-decaying, inert chemical dye solution through the pore space of a fluid saturated porous column. The objective of the paper is to present a complete study of the one-dimensional advective-dispersive transport problem by considering certain analytical solutions, experimental results and their comparisons with specific computational simulations. Dye concentrations obtained by means of an image processing method are used in conjunction with an analytical solution to identify the hydrodynamic dispersion coefficient that governs the advective-dispersive transport problem. The experimental results and identified parameters are also used to assess the computational estimates derived from several stabilized computational schemes available in the literature, for examining advection-dominated transport processes in porous media.     

A Decomposition Method for Solving Coupled Multi–Species Reactive Transport Problems

Concerns over the problems associated with mixed waste groundwater contamination have created a need for more complex models that can represent reactive contaminant fate and transport in the subsurface. In the literature, partial differential equations describing solute transport in porous media are solved either for a single reactive species in one, two or three dimensions, or for a limited number of reactive species in one dimension. Those solutions are constrained by many simplifying assumptions. Often, it is desirable to simulate transport in two or three dimensions for a more practical system that might have multiple reactive species. This paper presents a decomposition method to solve the partial differential equations of multi–dimensional, multi–species transport problems that are coupled by linear reactions. A matrix method is suggested as a tool for describing the reaction network. In this way, the level of complexity required to solve the multi–species reactive transport problem is significantly reduced.     

Effective Molecular Dynamics Model of Ionic Solutions for Large-Scale Calculations

A model of ionic solutions is proposed which can be used to calculate aqueous salt solutions in different nanostructures. The interaction potential of the model includes the Lennard-Jones potential and angularly averaged dipole–dipole and ion–dipole interactions. Lennard-Jones potential parameters for different ions are obtained. Characteristics of aqueous solutions at different salt concentrations are calculated using the molecular dynamics method. It is shown that the calculated values of the hydration shells of ions parameters are in good agreement with the theoretical and experimental data at a salt concentration of 1 mol/kg. The computational scheme used in the calculations is described. It is shown that calculations using the proposed model require less computing resources compared with the standard models of ionic solutions.     

Asymptotic Nonlinear Behaviors in Transverse Vibration of an Axially Accelerating Viscoelastic String

This paper investigates longtime dynamical behaviors of an axially accelerating viscoelastic string with geometric nonlinearity. Application of Newton's second law leads to a nonlinear partial-differential equation governing transverse motion of the string. The Galerkin method is applied to truncate the partial-differential equation into a set of ordinary differential equations. By use of the Poincare maps, the dynamical behaviors are presented based on the numerical solutions of the ordinary differential equations. The bifurcation diagrams are presented for varying one of the following parameter: the mean transport speed, the amplitude and the frequency of transport speed fluctuation, the string stiffness or the string dynamic viscosity, while other parameters are fixed.     

Interfacial dynamics‐based modelling of turbulent cavitating flows,Part‐1: Model development and steady‐state computations

The merits of transport equation‐based models are investigated by adopting an enhanced pressure‐based method for turbulent cavitating flows. An analysis of the mass and normal‐momentum conservation at a liquid–vapour interface is conducted in the context of homogeneous equilibrium flow theory, resulting in a new interfacial dynamics‐based cavitation model. The model offers direct interpretation of the empirical parameters in the existing transport‐equation‐based models adopted in the literature. This and three existing cavitation models are evaluated for flows around an axisymmetric cylindrical body and a planar hydrofoil, and through a convergent–divergent nozzle. Although all models considered provide qualitatively comparable wall pressure distributions in agreement with the experimental data, quantitative differences are observed in the closure region of the cavity, due to different compressibility characteristics of each cavitation model. In particular, the baroclinic effect of the vorticity transport equation plays a noticeable role in the closure region of the cavity, and contributes to the highest level of turbulent kinetic energy there. Copyright © 2004 John Wiley & Sons, Ltd.     

Supercritical vibration of nonlinear coupled moving beams based on discrete Fourier transform

Natural frequencies of nonlinear coupled planar vibration are investigated for axially moving beams in the supercritical transport speed ranges. The straight equilibrium configuration bifurcates in multiple equilibrium positions in the supercritical regime. The finite difference scheme is developed to calculate the non-trivial static equilibrium. The equations are cast in the standard form of continuous gyroscopic systems via introducing a coordinate transform for non-trivial equilibrium configuration. Under fixed boundary conditions, time series are calculated via the finite difference method. Based on the time series, the natural frequencies of nonlinear planar vibration, which are determined via discrete Fourier transform (DFT), are compared with the results of the Galerkin method for the corresponding governing equations without nonlinear parts. The effects of material parameters and vibration amplitude on the natural frequencies are investigated through parametric studies. The model of coupled planar vibration can reduce to two nonlinear models of transverse vibration. For the transverse integro-partial-differential equation, the equilibrium solutions are performed analytically under the fixed boundary conditions. Numerical examples indicate that the integro-partial-differential equation yields natural frequencies closer to those of the coupled planar equation.     

Hall and ion slip effects on peristaltic flow and heat transfer analysis with Ohmic heating

The peristaltic transport of a magnetohydrodynamic （MHD） fluid is exam- ined for both symmetric and asymmetric channels. Hall and ion slip effects are taken into account. The heat transfer is analyzed by considering the effects of viscous and Ohmic dissipations. The relevant flow problems are first modeled, and then the closed form solutions are constructed under the assumptions of long wavelength and low Reynolds number. The solutions are analyzed through graphical illustration. It is noted that the velocity increases but the temperature decreases with the increases in the Hall and ion slip parameters. The axial pressure gradient is less in magnitude in the presence of Hall and ion slip currents. The Hall and ion slip effects are to decrease the maximum pres- sure against which peristalsis works as a pump. The free pumping flux decreases with the increases in the Hall and ion slip parameters. The increases in the Hall and ion slip parameters result in an increase in the size of the trapped bolus.     

Theoretical modeling of microstructured liquids: a simple thermodynamic approach

A new method is presented for accounting for microstructural features of flowing complex fluids at the level of mesoscopic, or coarse-grained, models by ensuring compatibility with macroscopic and continuum thermodynamics and classical transport phenomena. In this method, the microscopic state of the liquid is described by variables that are local expectation values of microscopic features. The hypothesis of local thermodynamic equilibrium is extended to include information on the microscopic state, i.e., the energy of the liquid is assumed to depend on the entropy, specific volume, and microscopic variables. For compatibility with classical transport phenomena, the microscopic variables are taken to be extensive variables (per unit mass or volume), which obey convection-diffusion-generation equations. Restrictions on the constitutive laws of the diffusive fluxes and generation terms are derived by separating dissipation by transport (caused by gradients in the derivatives of the energy with respect to the state variables) and by relaxation (caused by non-equilibrated microscopic processes like polymer chain stretching and orientation), and by applying isotropy. When applied to unentangled, isothermal, non-diffusing polymer solutions, the equations developed according to the new method recover those developed by the Generalized Bracket [J. Non-Newtonian Fluid Mech. 23 (1987) 271; A.N. Beris, B.J. Edwards, Thermodynamics of Flowing Systems with Internal Microstructure, first ed., Oxford University Press, Oxford, 1994] and by the Matrix Model [J. Rheol. 38 (1994) 769]. Minor differences with published results obtained by the Generalized Bracket are found in the equations describing flow coupled to heat and mass transfer in polymer solutions. The new method is applied to entangled polymer solutions and melts in the general case where the rate of generation of entanglements depends nonlinearly on the rate of strain. A link is drawn between the mesoscopic transport equations of entanglements and conformation and the microscopic equation describing the configurational distribution of polymer segment stretch and orientation. Constraints are derived on the generation terms in the transport equations of entanglements and conformation, and the formula for the elastic stress is generalized to account for reversible formation and destruction of entanglements. A simplified version of the transport equation of conformation is presented which includes many previously published constitutive models, separates flow-induced polymer stretching and orientation, yet is simple enough to be useful for developing large-scale computer codes for modeling coupled fluid flow and transport phenomena in two- and three-dimensional domains with complex shapes and free surfaces.     

Analytical and numerical solution for viscous dissipation effect on the heat transfer in a deformable channel

In this paper, shooting method and homotopy perturbation technique are applied for the flow analysis of temporal energy transport in a deformation channel with isothermal walls. An incompressible viscous fluid fills the space inside the channel. Analytical and numerical solutions are developed for the momentum and energy equations. The viscous dissipation effects are taken into account. Graphs for pertinent flow parameters are sketched and discussed. Comparison between the analytical and numerical solutions indicates an excellent agreement. It is noticed that behaviors of Prandtl and Eckert numbers on the temperature are qualitatively similar. Copyright © 2011 John Wiley & Sons, Ltd.     

Models of the Interaction between a Laminar Boundary Layer and a Transonic Flow

Possible regimes of viscous-inviscid interaction at transonic external flow velocities are investigated. It is shown that different flow regimes can exist depending on the relation between such parameters as the disturbance amplitude and the Mach and Reynolds numbers. Corresponding mathematical models are formulated and the solutions of some problems describing linear regimes of disturbance development are obtained. The models developed make it possible to describe all the possible interaction regimes.     

Principal interval decomposition framework for POD reduced‐order modeling of convective Boussinesq flows

A principal interval decomposition (PID) approach is presented for the reduced‐order modeling of unsteady Boussinesq equations. The PID method optimizes the lengths of the time windows over which proper orthogonal decomposition (POD) is performed and can be highly effective in building reduced‐order models for convective problems. The performance of these POD models with and without using the PID approach is investigated by applying these methods to the unsteady lock‐exchange flow problem. This benchmark problem exhibits a strong shear flow induced by a temperature jump and results in the Kelvin–Helmholtz instability. This problem is considered a challenging benchmark problem for the development of reduced‐order models. The reference solutions are obtained by direct numerical simulations of the vorticity and temperature transport equations using a compact fourth‐order‐accurate scheme. We compare the accuracy of reduced‐order models developed with different numbers of POD basis functions and different numbers of principal intervals. A linear interpolation model is constructed to obtain basis functions when varying physical parameters. The predictive performance of our models is then analyzed over a wide range of Reynolds numbers. It is shown that the PID approach provides a significant improvement in accuracy over the standard Galerkin POD reduced‐order model. This numerical assessment of the PID shows that it may represent a reliable model reduction tool for convection‐dominated, unsteady‐flow problems. Copyright © 2015 John Wiley & Sons, Ltd.     

A nonlinear ALE‐FCT scheme for non‐equilibrium reactive solute transport in moving domains

In this paper, we present a conservative, positivity‐preserving, high‐resolution nonlinear ALE‐flux‐corrected transport (FCT) scheme for reactive transport models in moving domains. The mathematical model is a convection–diffusion equation with a nonlinear flux equation on the moving channel wall. The reactive transport is assumed to have dominant Péclet and Damköhler numbers, a phenomenon that often results in non‐physical negative solutions. The scheme presented here is proven to be mass conservative in time and positive at all times for a small enough Δt. Reactive transport examples are simulated using this scheme for its validation, to show its convergence, and to compare it against the linear ALE‐FCT scheme. The nonlinear ALE‐FCT is shown to perform better than the linear ALE‐FCT schemes for large time steps. Copyright © 2014 John Wiley & Sons, Ltd.