Solidification and Flow of a Binary Alloy Over a Moving Substrate

We study the solidification and flow of a binary alloy over a horizontally moving substrate. A situation in which the solid, liquid and mushy regions are separated by the stationary two-dimensional interfaces is considered. The self-similar solutions of the governing boundary layer equations are obtained, and their parametric dependence is analysed asymptotically. The effect of the boundary layer flow on the physical characteristics is determined. It is found that the horizontal pulling and the resulting flow in the liquid enhance the formation of the mushy region.     

On Nonlinear Instability and Stability for Stratified Shear Flow

An example of stratified shear flow is presented in which an explicit construction is given for unstable eigenvalues with smooth eigenfunctions for the Taylor--Goldstein equation. It is proved for any stratified, plane parallel shear flow that the unstable spectrum of the linear operator is purely discrete. A general theorem is then invoked to prove that the specific example is nonlinearly unstable. A sufficient condition for nonlinear stability for stratified shear flow is discussed.     

A Surrogate Modelling Approach Based on Nonlinear Dimension Reduction for Uncertainty Quantification in Groundwater Flow Models

Transport in Porous Media - In this paper, we develop a surrogate modelling approach for capturing the output field (e.g. the pressure head) from groundwater flow models involving a stochastic...     

Evolution Equations for Plane Cubically Nonlinear Elastic Waves

Consideration is given to the nonlinear theory of elastic waves with cubic nonlinearity. This nonlinearity is separated out, and the interaction of four harmonic waves is studied. The method of slowly varying amplitudes is used. The shortened and evolution equations, the first integrals of these equations (Manley–Rowe relations), and energy balance law for a set of four interacting waves (quadruplet) are derived. The interaction of waves is described using the wavefront reversal scheme     

Metal Flow during Friction Stir Welding of 7075-T651 Aluminum Alloy

Bronze foil of 0.1 mm thickness was placed between faying surfaces of two plates to be butt-welded as marker material to reveal the flow behavior of weld metal during friction stir welding of 7075-T651 aluminum alloy. By tracing the bronze foil fragments in the weld after welding, the metal flow behavior during the welding process was revealed. Besides, the tool forces in the welding process were measured by the octagonal loop resistance turning dynamometer to expound the periodic variation of metal flow pattern. Results show that the flow behavior of the weld metal is different along the thickness direction. The flow pattern presents a periodic variation, and a formula has been proposed to calculate the periodicity of the metal flow. In addition, the weld nugget zone presents a “spoon” shape and the fine grains at the spoon handle and those at the spoon bowl are originated from different zones. A plastic metal flow model in FSW was proposed based on the results. Furthermore, the formation of defects was explained by researching the weld metal flow behavior.     

Attractors for Second-Order Evolution Equations with a Nonlinear Damping

We study long-time dynamics of abstract nonlinear second-order evolution equations with a nonlinear damping. Under suitable hypotheses we prove existence of a compact global attractor and finiteness of its fractal dimension. We also show that any solution is stabilized to an equilibrium and estimate the rate of the convergence which, in turn, depends on the behaviour at the origin of the function describing the dissipation. If the damping is bounded below by a linear function, this rate is exponential. Our approach is based on far reaching generalizations of the Ceron–Lopes theorem on asymptotic compactness and Ladyzhenskayas theorem on the dimension of invariant sets. An application of our results to nonlinear damped wave and plate equations allow us to obtain new results pertaining to structure and properties of global attractors for nonlinear waves and plates.     

A Steady-State Upscaling Approach for Immiscible Two-Phase Flow

The paper presents a model for computing rate-dependent effective capillary pressure and relative permeabilities for two-phase flow, in 2 and 3 space-dimensions. The model is based on solving the equations for immiscible two-phase flow at steady-state, accounting for viscous and capillary forces, at a given external pressure drop. The computational performance of the steady-state model and its accuracy is evaluated through comparison with a commercial simulator ECLIPSE. The properties of the rate-dependent effective relative permeabilities are studied by way of computations using the developed steady-state model. Examples presented show the dependence of the effective relative permeabilities and capillary pressures, which incorporate the effects of fine scale wettability heterogeneity, on the external pressure drop, and thereby on the dimensionless macro-scale capillary number. The effective relative permeabilities converge towards the viscous limit functions as the capillary number tends to infinity. Special cases, when the effective relative permeabilities are rate-invariant, are also studied. The applicability of the steady-state upscaling algorithm in dynamic displacement situations is validated by comparing fine-gridded simulations in heterogeneous reservoirs against their homogenized counterparts. It is concluded that the steady-state upscaling method is able to accurately predict the dynamic behavior of a heterogeneous reservoir, including small scale heterogeneities in both the absolute permeability and the wettability.     

On the Evolution of the Flow Field in a Spark Ignition Engine

The development of the turbulent flow field inside a spark ignition engine is examined by large-eddy simulation (LES), from the intake flow to the tumble break-down. Ten consecutive cold flow engine cycles on a coarse and twenty cycles on a fine grid are simulated and compared to experiments of the same engine. The turbulent subgrid scales are modeled by the standard Smagorinsky and by the recently developed Sigma model. A comparison of the intake flow is made against Particle Image Velocimetry (PIV) measurements along horizontal and vertical lines and to an LES simulation performed by the Darmstadt group. Furthermore, we show the first LES comparison to Magnetic Resonance Velocimetry (MRV conducted by Freudenhammer et al.) measurements, which provided the 3D flow field inside a full scale dummy of the entire upper cylinder head including the valve seat region, at a time which mimics inflow conditions of the corresponding engine. Our LES is in good qualitative and quantitative agreement with the simulation and the experiments, with the notable exception of the measured in-cylinder pressure, which is discussed in detail and compared to 0D simulations and simulations from other groups. A criterion is proposed for estimating the number of cycles needed in a simulation, if experimental data is available. We put emphasis on the flow in the valve seat region, where turbulence is generated, and discuss the formation of the large scale tumble motion, including a comparison of the radial velocity fields on rolled-up planes around the valve seat. Here, spots of high velocities were found in the under flow region, which cannot been seen by the ensemble averaged MRV measurement. Within the compression stroke, a 2D vortex center identification algorithm is applied on slices inside the combustion chamber, yielding a 3D visualization of the tumble vortex, which is found to have a “croissant-like” shape. The tumble vortex trajectory is plotted on the symmetry plane and compared to measurements. Finally, we consider a modified definition of the (turbulent) integral length scale that provided further insight to the tumble break-down process.     

A Galerkin Approach for Power Spectrum Determination of Nonlinear Oscillators

A numerical method to estimate spectral properties of nonlinear oscillators with random input is presented. The stationary system response is expanded into a trigonometric Fourier series. A set of nonlinear algebraic equations, solved by Newtons method, leads to the determination of the unknown Fourier series coefficients of single samples of the response process. For cubic polynomial nonlinearities, closed-form expressions are used to find the nonlinear terms at each step of the solution scheme. Further, a simple procedure yields an approximation of an arbitrary nonlinearity by a cubic polynomial. Power spectral density estimates for the response process are constructed by averaging the square modulus of the computed Fourier coefficients over various samples or by means of well-established smoothing techniques of spectral analysis. Two applications are presented illustrating the effectiveness of the method as compared to statistical linearization and digital Monte-Carlo simulation.     

A Multiscale Approach for Geologically and Flow Consistent Modeling

Subsurface geological models are usually constructed on high-resolution grids in a way that various complexities and heterogeneities are depicted properly. Such models, however, cannot be used directly in the current flow simulators, as they are tied with high computational cost. Thus, using upscaling, by which one can produce flow consistent models that can alleviate the computational burden of flow simulators, is inevitable. Although the upscaling methods are able to reproduce the flow responses, they might not retain the initial geological assumptions. The reservoir models are initially constructed on uniform and high-resolution grids and then, if necessary, are upscaled to be used for flow simulations. A subsurface modeling approach that not only preserves the geological heterogeneity but also provides models that can be used, straight or with a small level of upscaling, in the flow simulators is desirable. In this paper, a new multiresolution method based on (1) the importance of conditioning well data and (2) being geologically and flow consistent is presented. This method discretizes the initial model into several regions based on the available data. Then, the initial assumed geological model is converted into, for example, various high- and low-resolution models. Next, the high-resolution model is used for regions with high-quality data (e.g., well locations), while the low-resolution model is used for the remaining areas. Finally, the patterns of these areas are interlocked, which result in a multiresolution geologically and flow consistent subsurface model. The accuracy of this method is demonstrated using two-phase flow simulation on four complex subsurface systems. The results indicate that the same flow responses, in a much less time, are reproduced using the multiscale models. The speed-up factor gained using the proposed method is also several orders of magnitude.     

Negative Saturation Approach for Non-Isothermal Compositional Two-Phase Flow Simulations

This article deals with developing a solution approach, called the non-isothermal negative saturation (NegSat) solution approach. The NegSat solution approach solves efficiently any non-isothermal compositional flow problem that involves phase disappearance, phase appearance, and phase transition. The advantage of the solution approach is that it circumvents using different equations for single-phase and two-phase regions and the ensuing unstable procedure. This paper shows that the NegSat solution approach can also be used for non-isothermal systems. The NegSat solution approach can be implemented efficiently in numerical simulators to tackle modeling issues for mixed CO₂–water injection in geothermal reservoirs, thermal recovery processes, and for multicontact miscible and immiscible gas injection in oil reservoirs. We illustrate the approach by way of example to cold mixed CO₂–water injection in a 1D geothermal reservoir. The solution is compared with an analytical solution obtained with the wave-curve method (method of characteristics) and shows excellent agreement. A complete set of simulations is carried out, which identifies six bifurcations. The two main bifurcations are (1) when the most downstream compositional wave is replaced by a compositional shock and (2) when an extra Buckley–Leverett rarefaction appears. The plot of the useful energy (exergy) versus the CO₂ storage capacity shows a Z-shape. The top horizontal part represents a branch of high exergy recovery/relatively lower storage capacity, whereas the bottom part represents a branch of lower exergy recovery/higher storage capacity.     

同心旋转圆柱间粘弹性流的非线性动力学模型

基于小间隙假设,将速度场和应力场有杉Ｆｏｕｒｉｅｒ展开式截断近似,推导了同心旋转圆柱间Ｏｌｄｒｏｙｄ－Ｂ型流体的六维动力系统,探讨了高分子添加对滑动轴承间油膜稳定性的影响。结果表明,少量的高分子添加剂具有推迟流体层流失稳的作用。     

同心旋转圆柱间强弹性流的非线性稳定性分析

基于Oldroyd-B型粘弹性流体模型，采用同心旋转圆柱间非线性动力系统分析了流体的弹性对轴对称Taylor涡稳定性的影响．分析结果表明，对于弱弹性流体，Taylor涡出现时，系统存在超临界分岔；而对于强弹性流体则出现亚临界分岔．在小间隙大扰动条件下，采用有限差分法分析了非线性效应对系统稳定性的影响．数值计算结果表明，随着流动速度的增加，润滑油膜的失稳结构与流体的弹性有关，对于弱弹性流，流体以同宿轨道分岔失稳；强弹性流则出现倍周期分岔，直至发生混沌，流场最终发展为湍流．     

Modular Nonlinear Filter Stabilization of Methods for Higher Reynolds Numbers Flow

Stabilization using filters is intended to model and extract the energy lost to resolved scales due to nonlinearity breaking down resolved scales to unresolved scales. This process is highly nonlinear and yet current models for it use linear filters to select the eddies that will be damped. In this report we consider for the first time nonlinear filters which select eddies for damping (simulating breakdown) based on knowledge of how nonlinearity acts in real flow problems. The particular form of the nonlinear filter allows for easy incorporation of more knowledge into the filter process and its computational complexity is comparable to calculating a linear filter of similar form. We then analyze nonlinear filter based stabilization for the Navier?CStokes equations. We give a precise analysis of the numerical diffusion and error in this process.     

Study on Exact Analytical Solutions for Two Systems of Nonlinear Evolution Equations

ＩｎｔｒｏｄｕｃｔｉｏｎＤｕｒｉｎｇｔｈｅｃｏｕｒｓｅｏｆｓｔｕｄｙｉｎｇｔｈｅｗａｔｅｒｗａｖｅ,ｍａｎｙｃｏｍｐｌｅｔｅｌｙｉｎｔｅｇｒａｂｌｅｍｏｄｅｌｓｗｅｒｅｏｂｔａｉｎｅｄ ,ｓｕｃｈａｓＫｄＶｅｑｕａｔｉｏｎ ,ｍＫｄＶｅｑｕａｔｉｏｎ ,(2 1 )＿ｄｉｍｅｎｓｉｏｎａｌＫＰｅｑｕａｔｉｏｎ ,ｃｏｕｐｌｅｄＫｄＶｅｑｕａｔｉｏｎｓ,ｖａｒｉａｎｔＢｏｕｓｓｉｎｅｓｑｅｑｕａｔｉｏｎｓ ,ＷＫＢｅｑｕａｔｉｏｎｓｅｔｃ .[1- 13 ].Ｉｎｏｒｄｅｒｔｏｆｉｎｄｅｘｐｌｉｔｉｃｅｘａｃｔｓｏｌｕｔｉｏ…     

Comparison between Lever and Scheil Rules for Modeling of Microporosity Formation during Solidification

Segregation and microporosity formation are two important physicalphenomena that occur during solidification of binary alloys. The aim ofthis study is to investigate the effect of the model of solute diffusionat the local scale (which means at the local scale of the microscopicrepresentative elemental average volume (REV)) on solute transport andthe microporosity formation during this process. The Scheil rule and thelever rule are used to describe the solute diffusion at the local scale.Results indicate that solute diffusion at the local scale is animportant factor in microporosity formation. Also, microporosityformation slightly reduces inverse segregation because it partiallycompensates for shrinkage. The increase of the external pressure at thefree surface or the decrease of the initial hydrogen concentration inthe molten alloy can be effectively utilized to control microporosityformation.     

A General Data-Based Approach for Developing Reduced-Order Models of Nonlinear MDOF Systems

A general procedure is presented for analyzing dynamic response measurements from complex multi-degree-of-freedom nonlinear systems incorporating arbitrary types of nonlinear elements. The analysis procedure develops a reduced-order, nonlinear model whose format is convenient for numerical simulation studies. No information about the systems mass properties is needed, and only the applied excitations and corresponding response are needed to develop the model whose dimension is compatible with the number of available sensors. The utility of the approach is demonstrated by means of a three-degree-of-freedom system incorporating polynomial-type nonlinear features with hardening as well as softening characteristics.     

Global Nonlinear Stability for Steady Ideal Fluid Flow in Bounded Planar Domains

We prove stability of steady flows of an ideal fluid in a bounded, simply connected, planar region, that are strict maximisers or minimisers of kinetic energy on an isovortical surface. The proof uses conservation of energy and transport of vorticity for solutions of the vorticity equation with initial data in L^p for p>4/3. A related stability theorem using conservation of angular momentum in a circular domain is also proved.