Vortex Funnel Formation by a Mass Sink within the Shallow-Water Model

The problem of the evolution of an axisymmetric vortex in the presence of given mass sinks is considered within the shallow-water model. Using the cyclostrophic balance approximation, i.e. the balance between pressure gradient and centrifugal force, the nonlinear process of development of a steady vortex flow regime (vortex funnel) is described.     

Viscoelastic flow of polymer solutions around a periodic,linear array of cylinders: comparisons of predictions for microstructure and flow fields

Numerical simulation is used to investigate the flow of polymer solutions around a periodic, linear array of cylinders by using three constitutive equations derived from kinetic theory of dilute polymer solutions: the Giesekus model; the finitely extensible, nonlinear elastic dumbbell model with Peterlin's approximation (FENE-P); and the FENE dumbbell model of Chilcott–Rallison (CR). In the Giesekus model, intramolecular forces are described by a Hookean spring, whereas a finitely extensible spring whose modulus is given by the Warner approximation is used in both the FENE-P and CR models. Hydro dynamic drag on the beads is taken to be anisotropic for the Giesekus model and isotropic for the other two models. The CR and FENE-P models differ subtly in their approximate treatment of the nonlinear force law. The three models exhibit very similar rheological behavior in viscometric flow and steady elongational flow, with the notable exception that the viscosity for the CR model is shear-rate independent. Finite element simulations are performed by using two different formulations: the elastic-viscous split-stress gradient (EVSS-G) method and a new variant of this formulation, the discrete EVSS-G (DEVSS-G) formulation, in which the elliptic stabilization term is added only to the discrete version of the momentum equation, and the constitutive equation is solved directly in terms of the polymer contribution to the stress tensor. Calculations are performed for all models up to a Weissenberg number We, where the configuration tensor 〈QQ〉 loses positive definiteness. However, by locally refining the mesh in the gap region, the positive definiteness of 〈QQ〉 is recovered. The flow and stress fields predicted by the three constitutive equations are qualitatively similar. A `birefringent strand' of highly stretched polymer molecules, which appears to emanate from the rear stagnation point in the cylinder, strengthens as We is increased. Not surprisingly, the molecular extension computed for the Giesekus model is considerably larger than that of the two FENE spring models. The drag force on the cylinders differs for the FENE-P and CR models, because of the difference in the shear-thinning viscosity resulting from the different approximations used in these models.     

Law of nonlinear flow in saturated clays and radial consolidation

It was derived that micro-scale amount level of average pore radius of clay changed from 0.01 to 0.1 micron by an equivalent concept of flow in porous media.There is good agreement between the derived results and test ones.Results of experiments show that flow in micro-scale pore of saturated clays follows law of nonlinear flow.Theoretical analyses demonstrate that an interaction of solid-liquid interfaces varies inversely with permeability or porous radius.The interaction is an important reason why nonlinear flow in saturated clays occurs.An exact mathematical model was presented for nonlinear flow in micro-scaie pore of saturated clays.Dimension and physical meanings of parameters of it are definite.A new law of nonlinear flow in saturated clays was established.It can describe characteristics of flow curve of the whole process of the nonlinear flow from low hydraulic gradient to high one.Darcy law is a special case of the new law.A math- ematical model was presented for consolidation of nonlinear flow in radius direction in saturated clays with constant rate based on the new law of nonlinear flow.Equations of average mass conservation and moving boundary,and formula of excess pore pressure distribution and average degree of consolidation for nonlinear flow in saturated clay were derived by using an idea of viscous boundary layer,a method of steady state in stead of transient state and a method of integral of an equation.Laws of excess pore pressure distribution and changes of average degree of consolidation with time were obtained.Re- suits show that velocity of moving boundary decreases because of the nonlinear flow in saturated clay.The results can provide geology engineering and geotechnical engineering of saturated clay with new scientific bases.Calculations of average degree of consolidation of the Darcy flow are a special case of that of the nonlinear flow.     

Wavelet method applied to specific adhesion of elastic solids via molecular bonds

We propose a wavelet method to analyze the stochastic-elastic problem of specific adhesion between two elastic solids via ligand-receptor bond clusters, which is governed by a nonlinear integro-differential equation with a singular Cauchy kernel to describe the mean-field coupling between deformation of elastic materials and stochastic behavior of the molecular bonds. To solve this problem, Galerkin method based on a wavelet approximation scheme is adopted, and special treatment which transforms the singular Cauchy kernel into a smooth one has been proposed to avoid the cumbersome calculation of singular integrals. Numerical results demonstrate that the method is fully capable of solving the specific adhesion problems with complex nonlinear and singular equations. Based on the proposed method, investigations are performed to reveal the relation between steady-state pulling force and mean surface separation under different stress concentration indexes, which is crucial for assembling the overall constitutive relations for multicellular tumor spheroids and polymer-matrix microcomposites.     

From super-Bloch oscillations to sudden self-trapping in Bose–Einstein condensates with inter-atomic interactions

Nonlinear Dynamics - We study super-Bloch oscillations of ultracold atoms loaded in a lattice potential. We consider nonlinear contribution emerging from a mean-field approximation for the...     

Sharp approximation to a filtration problem by maximum principle

Consider an incompressible fluid, filtrating through a saturated cylindrical porous layer with rectangular cross-section. A steady pressure gradient, parallel to the axis of the layer, drives a one-directional stationary non recirculating flow when the Darcy law has to include inertial and viscous corrections. This is the case, for instance, when the porosity of the medium or the seeping flow rate are not very small. The resulting nonlinear problem belongs to a class of equations which was proved to have positive solutions. It also satisfies a comparison principle from which approximations from above and from below are derived for the steady flow. The estimate from above is the flat profile which solves the Darcy-Forchheimer equation, which does not take account of viscous effects, and the approximation is excellent when the modified Darcy number is small, under the additional condition that the Forchheimer coefficient be small also. The flow still solves the problem when gradient forces, orthogonal to the axis of the layer, are also present.     

Upscaling Forchheimer law

We investigate the high velocity flow in heterogeneous porous media. The model is obtained by upscaling the flow at the heterogeneity scale where the Forchheimer law is assumed to be valid. We use the method of multiple scale expansions, which gives rigorously the macroscopic behaviour without any prerequisite on the form of the macroscopic equations. We show that Forchheimer law does not generally survive upscaling. The macroscopic flow law is strongly non-linear and anisotropic. A 2-point Padé approximation of the flow law in the form of a Forchheimer law is given. However, this approximation is generally poor. These results are illustrated in two particular cases: a layered composite porous media and a composite constituted by a square array of circular porous inclusions embedded in a porous matrix. We show that non-linearities are sources of anisotropy.     

非线性周期波及破碎波对杆件的作用力

本文采用十五阶 Stokes 波的 Pade 逼近,获得了与实验较为一致的流场,并且利用已有的破碎波的速度、加速度场,计算了非线性波和破碎波对各种杆件作用力,比较了它们的主要特征,为海洋工程设计提供依据。     

A detailed comparison of various FENE dumbbell models

The rheological behaviour of dilute solutions of finitely extensible non-linear elastic (FENE) dumbbells in both steady state and transient shear and simple elongational flow is investigated. Three dumbbell models are compared: the original FENE model with the Warner spring force, which is treated by brownian dynamics simulations, and the FENE-P model based on the Peterlin approximation and the FENE-CR model as suggested by Chilcott and Rallison, which are treated by standard numerical techniques. It is shown that in the linear viscoelastic limit and in steady state flows the behaviour is similar, except for the FENE-CR dumbbell in shear flow, modelling a Boger fluid. In transient flows larger differences appear.     

The reptation model with segmental stretch

The Doi-Edwards model with segmental stretch and a non-linear finitely extensible spring law is described and examined in simple flow situations where analytic results are derivable; namely oscillatory flow and steady state flow at high deformation rates. The model is shown to be consistent with the Bueche-Ferry hypothesis in fast large strain unidirectional flows but to violate this rule in small strain reversing flows. The discrepancy is identified with a preaveraging approximation used to describe the relative tube-chain velocity. Experimentally verifiable scaling rule for the birefringence as a universal function of a planar flow-type parameter and deformation rate are identified. Sensitivity to the extensional flow character, absent in the original tube model, manifests itself with the introduction of segmental stretch. Although the model generates a non-separable memory function kernel the deformation dependence of the memory function is quantitatively shown to have negligible impact on the predicted theological properties relative to the original Doi-Edwards model. With this simplification, relatively uncomplicated approximations to the segmental stretch model can be deduced.     

Effects of inertia and surface tension on a power-law fluid flowing down a wavy incline

We consider a thin film of a power-law liquid flowing down an inclined wall with sinusoidal topography. Based on the von Kármán–Pohlhausen method an integral boundary-layer model for the film thickness and the flow rate is derived. This allows us to study the influence of the non-Newtonian properties on the steady free surface deformation. For weakly undulated walls we solve the governing equation analytically by a perturbation approach and find a resonant interaction of the free surface with the wavy bottom. Furthermore, the analytical approximation is validated by numerical simulations. Increasing the steepness of the wall reveals that nonlinear effects like the resonance of higher harmonics grow in importance. We find that shear-thickening flows lead to a decrease while shear thinning flows lead to an amplification of the steady free surface. A linear stability analysis of the steady state shows that the bottom undulation has in most cases a stabilizing influence on the free surface. Shear thickening fluids enhance this effect. The open questions which occurred in the linear analysis are then clarified by a nonlinear stability analysis. Finally, we show the important role of capillarity and discuss its influence on the steady solution and on the stability.     

Primary Bjerknes force in a three-dimensional nonlinear resonator: Numerical simulations

The effect of the primary Bjerknes force caused by a three-dimensional nonlinear standing ultrasonic wave on a population of nonlinear oscillating bubbles is studied in this paper by analyzing the results obtained from simulations performed with a numerical model at low and moderate pressure amplitudes. Small air bubbles are evenly distributed in a water filled cavity excited at resonance for which axial symmetry is assumed. Both the bubble oscillation variable and the pressure variable are unknown in the nonlinear set of coupled differential equations that describes the interaction of ultrasound and bubbles. Simulation results show that the three-dimensional primary Bjerknes force field is strongly amplitude dependent. We also analyze whether taking the term in the differential system that defines the nonlinear behavior of the pressure field into account is determinant or not on the computation of the force field. The corresponding results corroborate the one-dimensional conclusions on the fundamental importance of considering this nonlinear acoustic term to obtain an accurate approximation of the force in a cavity.     

Flow of a Weakly Conducting Fluid in a Channel Filled with a Darcy–Brinkman–Forchheimer Porous Medium

We investigate in this article, the fully developed flow in a fluid-saturated channel filled with a Darcy–Brinkman–Forchheimer porous medium, which is conducted with an electrically varying parallel Lorentz force. The Lorentz force varies exponentially in the vertical direction due to low fluid electrical conductivity and the special arrangement of the magnetic and electric fields at the lower plate. With the homotopy analysis method (HAM), a particularly effective technique in solving nonlinear problems, analytical approximation series solutions with high accuracy are derived for fluid velocity and the results are illustrated in form of figures. All these flows are new and are presented for the first time in the literature.     

REDUCED-ORDER MODELS OF UNSTEADY TRANSONIC VISCOUS FLOWS IN TURBOMACHINERY

The proper orthogonal decomposition (POD) technique is applied in the frequency domain to obtain a reduced-order model of the unsteady flow in a transonic turbomachinery cascade of oscillating blades. The flow is described by a inviscid—viscous model, i.e. a full potential equation outer flow model and an integral equation boundary layer model. The nonlinear transonic steady flow is computed first and then the unsteady flow is determined by a small perturbation linearization about the nonlinear steady solution. Solutions are determined for a full range of frequencies and validated. The full model results and the POD method are used to construct a reduced-order model in the frequency domain. A cascade of airfoils forming the Tenth Standard Configuration is investigated to show that the reduced-order model with only 15–75 degrees of freedom accurately predicts the unsteady response of the full system with approximately 15 000 degrees of freedom.     

Numerical modeling of the acoustic wave dynamics on the basis of a kinetic approach

The results of the numerical modeling of a viscous compressible flow on a finite one-dimensional interval with periodic boundary conditions on the basis of a model kinetic equation are presented. The flow is excited by a small time-dependent force, periodic in space and time. It is shown that when the periodicity interval length L ≥ 10³ nonlinear steady oscillations with sharp temporal and spatial variations of the parameters arise. The system dynamics are investigated by means a pseudospectral method using the standard fast Fourier transformation procedure.     

Permeability up-scaling using Haar Wavelets

In the context of flow in porous media, up-scaling is the coarsening of a geological model and it is at the core of water resources research and reservoir simulation. An ideal up-scaling procedure preserves heterogeneities at different length-scales but reduces the computational costs required by dynamic simulations. A number of up-scaling procedures have been proposed. We present a block renormalization algorithm using Haar wavelets which provide a representation of data based on averages and fluctuations. In this work, absolute permeability will be discussed for single-phase incompressible creeping flow in the Darcy regime, leading to a finite difference diffusion type equation for pressure. By transforming the terms in the flow equation, given by Darcy’s law, and assuming that the change in scale does not imply a change in governing physical principles, a new equation is obtained, identical in form to the original. Haar wavelets allow us to relate the pressures to their averages and apply the transformation to the entire equation, exploiting their orthonormal property, thus providing values for the coarse permeabilities. Focusing on the mean-field approximation leads to an up-scaling where the solution to the coarse scale problem well approximates the averaged fine scale pressure profile.     

Robust tracking control method based on composite nonlinear feedback technique for linear systems with time-varying uncertain parameters and disturbances

In this paper, the composite nonlinear feedback control method is considered for robust tracking and model following of uncertain linear systems. The control law guarantees that the tracking error decreases asymptotically to zero in the presence of time varying uncertain parameters and disturbances. For performance improvement of the dynamical system, the proposed robust tracking controller consists of linear and nonlinear feedback parts without any switching element. The linear feedback law is designed to allow the closed loop system have a small damping ratio and a quick response while the nonlinear feedback law increases the damping ratio of the system as the system output approaches the output of the reference model. A new collection of different nonlinear functions used in the control law are offered to improve the reference tracking performance of the system. The proposed robust tracking controller improves the transient performance and steady state accuracy simultaneously. Finally, the simulations are provided to verify the theoretical results.     

Gas-kinetic unified algorithm for plane external force-driven flows covering all flow regimes by modeling of Boltzmann equation

The nonequilibrium steady gas flows under the external forces are essentially associated with some extremely complicated nonlinear dynamics, due to the acceleration or deceleration effects of the external forces on the gas molecules by the velocity distribution function. In this article, the gas-kinetic unified algorithm (GKUA) for rarefied transition to continuum flows under external forces is developed by solving the unified Boltzmann model equation. The computable modeling of the Boltzmann equation with the external force terms is presented at the first time by introducing the gas molecular collision relaxing parameter and the local equilibrium distribution function integrated in the unified expression with the flow state controlling parameter, including the macroscopic flow variables, the gas viscosity transport coefficient, the thermodynamic effect, the molecular power law, and molecular models, covering a full spectrum of flow regimes. The conservative discrete velocity ordinate (DVO) method is utilized to transform the governing equation into the hyperbolic conservation forms at each of the DVO points. The corresponding numerical schemes are constructed, especially the forward-backward MacCormack predictor-corrector method for the convection term in the molecular velocity space, which is unlike the original type. Some typical numerical examples are conducted to test the present new algorithm. The results obtained by the relevant direct simulation Monte Carlo method, Euler/Navier-Stokes solver, unified gas-kinetic scheme, and moment methods are compared with the numerical analysis solutions of the present GKUA, which are in good agreement, demonstrating the high accuracy of the present algorithm. Besides, some anomalous features in these flows are observed and analyzed in detail. The numerical experience indicates that the present GKUA can provide potential applications for the simulations of the nonequilibrium external-force driven flows, such as the gravity, the electric force, and the Lorentz force fields covering all flow regimes.     

Simultaneous measurement of transient stress and flow birefringence in one-sided compression (biaxial extension) of a polymer melt

Summary Investigation of time dependent behaviour of a polystyrene melt is carried out with the aid of a new apparatus for biaxial extension. Use is made of the method of two impinging fluid streams guided by lubricated trumpet shaped metal walls. The flow birefringence is measured in the plane of symmetry and, at the same time, the force is measured which tends to separate the trumpets. The linear stress-optical relation turns out to be valid in this new flow geometry. An accurate value for the stress-optical coefficient can be determined from the relaxation experiments. The stress build-up as calculated from the optical measurements, is compared with the pertinent result of the theory of linear viscoelasticity. For the desired interconversion of dynamic moduli use is made of the approximation by Schwarzl and Struik. The steady state measurements are checked by the results of the non-linear model of Acierno et al.With 16 figures and 2 tables