Newton's method for solving heat transfer,fluid flow and interface shapes in a floating molten zone

Newton's method is applied to the finite volume approximation for the steady state heat transfer, fluid flow and unknown interfaces in a floating molten zone. The streamfunction/vorticity and temperature formulation of the Navier–Stokes and energy equations and their associated boundary conditions are written in generalized curvilinear co-ordinates and conservative law form with the Boussinesq approximation. During Newton iteration the ILU(0) preconditioned GMRES matrix solver is applied for solving the linear system, where the sparse Jacobian matrix is estimated by finite differences. Nearly quadratic convergence of the method is observed. Sample calculations are reported for sodium nitrate, a high-Prandtl-number material (Pr = 9.12). Both natural convection and thermocapillary flow as well as an overall mass balance constraint in the molten zone are considered. The effects of convection and heat input on the flow patterns, zone position and interface shapes are illustrated. After the lens effect due to the molten zone is considered, the calculated flow patterns and interface shapes are compared with the observed ones and are found to be in good agreement.     

A non-iterative transformation method for Newton's free boundary problem

In book II of Newton's Principia Mathematica of 1687 several applicative problems are introduced and solved. There, we can find the formulation of the first calculus of variations problem that leads to the first free boundary problem of history. The general calculus of variations problem is concerned with the optimal shape design for the motion of projectiles subject to air resistance. Here, for Newton's optimal nose cone free boundary problem, we define a non-iterative initial value method which is referred in the literature as a transformation method. To define this method we apply invariance properties of Newton's free boundary problem under a scaling group of point transformations. Finally, we compare our non-iterative numerical results with those available in the literature and obtained via an iterative shooting method. We emphasize that our non-iterative method is faster than shooting or collocation methods and does not need any preliminary computation to test the target function as the iterative method or even provide any initial iterate. Moreover, applying Buckingham Pi-Theorem we get the functional relation between the unknown free boundary and the nose cone radius and height.     

A fully implicit direct Newton's method for the steady-state Navier–Stokes equations

Newton's method and banded Gaussian elimination can be a CPU efficient method for steady-state solutions to two-dimensional Navier–Stokes equations. In this paper we look at techniques that increase the radius of convergence of Newton's method, reduce the number of times the Jacobian must be factored, and simplify evaluation of the Jacobian. The driven cavity and natural convection problems are used as test problems, and finite volume discretization is employed.     

Interactions between a rectangular cylinder and a free-surface flow

The salient features of the interaction between a free-surface flow and a cylinder of rectangular cross-section are investigated and discussed. Laboratory-scale experiments are performed in a water channel under various flow conditions and elevations of the cylinder above the channel floor. The flow field is characterized on the basis of time-averaged and fluctuating local velocity measurements. Dynamic loadings on the cylinder are measured by two water-insulated dynamometers placed inside the cylinder structure. Starting from frequency and spectral analyses of the force signals, insights on the relationship between force dominant frequencies and the Strouhal number of the vortex shedding phenomenon are provided. Experimental results highlight the strong influence of the asymmetric configuration imposed by the two different boundary conditions (free surface and channel floor) on (i) the mean force coefficients and (ii) the vortex shedding frequencies. We provide an analysis of the nature of the dependence of average force coefficients on relevant dimensionless groups, i.e., the Reynolds number, normalized flow depth and cylinder submersion.     

模拟微可压粘弹性流体的WCCBS_SU方法

针对微可压缩粘弹性流动问题,发展了微可压缩流的WCCBS方法,详细推导了基于Oldroyd-B本构模型的WCCBS_SU方法的求解过程。在流场微可压的条件下,分别对平面Poiseuille流和4：1粘弹性收缩流进行了数值模拟。Poiseuille流在不同We数下数值结果与解析解的比较,验证了本文方法具有较高的精度和较好...     

An asymptotic iteration method for the numerical analysis of near-critical free-surface flows

Satisfying the boundary conditions at the free surface may impose severe difficulties to the computation of turbulent open-channel flows with finite-volume or finite-element methods, in particular, when the flow conditions are nearly critical. It is proposed to apply an iteration procedure that is based on an asymptotic expansion for large Reynolds numbers and Froude numbers close to the critical value 1.The iteration procedure starts by prescribing a first approximation for the free surface as it is obtained from solving an ODE that has been derived previously by means of an asymptotic expansion (Grillhofer and Schneider, 2003). The numerical solution of the full equations of motion then gives a surface pressure distribution that differs from the constant value required by the dynamic boundary condition. To determine a correction to the elevation of the free surface we next solve an ODE that is obtained from the asymptotic analysis of the flow with a prescribed pressure disturbance at the free surface. The full equations of motion are then solved for the corrected surface, and the procedure is repeated until criteria of accuracy for surface elevation and surface pressure, respectively, are satisfied.The method is applied to an undular hydraulic jump as a test case.     

Efficient solving method for unsteady incompressible interfacial flow problems

Unsteady interfacial problems, considered in an Eulerian form, are studied. The phenomena are modeled using the incompressible viscous Navier–Stokes equations to get the velocity field and an advection equation to predict interface evolutions. The momentum equation is solved by means of an implicit hybrid augmented Lagrangian–Projection method, whereas an explicit characteristic method coupled with a TVD SUPERBEE scheme is applied to the advection equation. The velocity components and the pressure are discretized on staggered grids with finite volumes. Emphasis is on the accuracy and robustness of the techniques described before. A precise explanation on the validation phase will be given, which uses such tests as the advection of a step function or Zalesak's problem to improve the calculation of the interface. The global approach is used on a physically hard interfacial test with strong disparities between viscosities and densities. Copyright © 1999 John Wiley & Sons, Ltd.     

Effects of the Jacobian evaluation on Newton's solution of the Euler equations

Newton's method is developed for solving the 2‐D Euler equations. The Euler equations are discretized using a finite‐volume method with upwind flux splitting schemes. Both analytical and numerical methods are used for Jacobian calculations. Although the numerical method has the advantage of keeping the Jacobian consistent with the numerical residual vector and avoiding extremely complex analytical differentiations, it may have accuracy problems and need longer execution time. In order to improve the accuracy of numerical Jacobians, detailed error analyses are performed. Results show that the finite‐difference perturbation magnitude and computer precision are the most important parameters that affect the accuracy of numerical Jacobians. A method is developed for calculating an optimal perturbation magnitude that can minimize the error in numerical Jacobians. The accuracy of the numerical Jacobians is improved significantly by using the optimal perturbation magnitude. The effects of the accuracy of numerical Jacobians on the convergence of the flow solver are also investigated. In order to reduce the execution time for numerical Jacobian evaluation, flux vectors with perturbed flow variables are calculated only for neighbouring cells. A sparse matrix solver that is based on LU factorization is used. Effects of different flux splitting methods and higher‐order discretizations on the performance of the solver are analysed. Copyright © 2005 John Wiley & Sons, Ltd.     

A pressure-smoothing scheme for incompressible flow problems

A pressure-smoothing scheme for Stokes and Navier–Stokes flows of Newtonian fluids and for Stokes flow of Maxwell fluids is described. The stress deviator obtained from the calculated velocity field is substituted into the governing equilibrium equation. The resulting equation is then solved to obtain a new, smoothed pressure by a least square finite element method.     

A total linearization method for solving viscous free boundary flow problems by the finite element method

In this paper a total linearization method is derived for solving steady viscous free boundary flow problems (including capillary effects) by the finite element method. It is shown that the influence of the geometrical unknown in the totally linearized weak formulation can be expressed in terms of boundary integrals. This means that the implementation of the method is simple. Numerical experiments show that the iterative method gives accurate results and converges very fast.     

A general methodology for investigating flow instabilities in complex geometries: application to natural convection in enclosures

This paper presents a general methodology for studying instabilities of natural convection flows enclosed in cavities of complex geometry. Different tools have been developed, consisting of time integration of the unsteady equations, steady state solving, and computation of the most unstable eigenmodes of the Jacobian and its adjoint. The methodology is validated in the classical differentially heated cavity, where the steady solution branch is followed for vary large values of the Rayleigh number and most unstable eigenmodes are computed at selected Rayleigh values. Its effectiveness for complex geometries is illustrated on a configuration consisting of a cavity with internal heated partitions. We finally propose to reduce the Navier–Stokes equations to a differential system by expanding the unsteady solution as the sum of the steady state solution and of a linear combination of the leading eigenmodes. The principle of the method is exposed and preliminary results are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Forces on a spherical particle in shear flow at intermediate Reynolds numbers

When particles are submerged in a shear flow, there are lateral (lift) forces on the particles, and these lateral forces affect the dispersion of the particles very much. Recent literature survey indicates that there are large discrepancies among the results from the previous numerical investigations on this subject. A small computational domain ranging between 20–30 sphere radii was used in all the previous numerical investigations. However, the result from the present study reveals that the value of lift coefficient strongly depends on the size of computational domain. To provide correct numerical data and physical interpretation for the forces on a spherical particle in linear shear flow, accurate numerical computations were performed for 5≤Re≤200 using a computational domain of 101 sphere radii.     

A numerical scheme for solving a periodically forced Reynolds equation

A modified Reynolds equation is used to model the air‐film in a high‐speed squeeze‐film bearing. The axial position of the bearing stator is prescribed as a finite amplitude periodic oscillation. A numerical approach is considered for solving the uncoupled and coupled periodic problems associated with this model. The uncoupled problem requires the computation of the squeeze‐film dynamics when the rotor is held at a fixed axial position and the coupled problem incorporates the additional air–rotor interaction since the rotor position is unknown and modelled as a spring‐mass‐damper system. The details of a Fourier spectral collocation scheme are provided for the reduction of the modified Reynolds equation to a system of non‐linear, first‐order ordinary differential equations in space. Using the Matlab boundary value problem solver bvp4c this system of equations is solved to give the periodic pressure distributions and rotor heights. The high degree of accuracy in the spectral collocation scheme is demonstrated through comparison with an appropriate analytical solution. Further analysis indicates that the direct periodic solver is at least 10 times faster than the equivalent Crank–Nicholson finite‐difference scheme. For changing values of a selected physical parameter the method of arc‐length continuation is employed to track branches of solutions computed using the spectral collocation scheme. A selection of results is presented to demonstrate the range of accessible solutions and the robust nature of the numerical scheme. Copyright © 2010 John Wiley & Sons, Ltd.     

A numerical method for solving equilibrium problems of masonry-like solids

This paper proposes a numerical method for the solution of equilibrium problems of solids which do not support tension. Some boundary-value problems are solved numerically and the solution obtained is compared to the exact one.

---

Sommario In questo lavoro viene proposto un metodo numerico per la soluzione di problemi di equilibrio di solidi non resistenti a trazione. Vengono successivamente risolti numericamente alcuni problemi di equilibrio e la soluzione ottenuta è confrontata con quella esatta.

     

Interpolation perturbation method for solving nonlinear problems

I.IntroductionInthispaper,onthebasisofRefll],usingtheinterPotationPerturbationmethod,theauthorseekstosolveseveralnon-Iinearproblems.Itsmainpointsare:Introducinganinterpolationfunction,wedeterminethisfunctionbyusingtheperturbationmethodandthenproceedtoseek…     

A finite element method for high Reynolds number viscous fluid flow using two step explicit scheme

This paper presents the finite element method for the analysis of unsteady viscous flow of fluid at high Reynolds numbers. The method is based on the explicit numerical integration scheme in time and uses three node triangular finite elements. For the convenience of the formulation, slight compressibility is considered. For the explicit scheme, the selective lumping two step scheme has been successfully employed. Vortex shedding behind a cylinder has been computed and compared with the conventional experimental results. The results agree favourably when both schemes are compared.     

自由面势流问题边界元法的数值色散误差

以连续及离散Fourier分析研究自由面势流问题边界元法的数值色散误差,并从理论上探讨有关计算中数值色散误差的改善问题.研究表明:对于该问题的数值色散误差而言,重要的在于以问题相应的离散算子考察计及各种数值手段后的总体色散误差,而非仅考虑该数值手段自身的数值色散误差大小.高阶面元、自由面域外奇点或适当的耦合方法是降低有关问题算子总体色散误差的较好选择.     

两点边值问题的小波配点法

根据多分辨分析,提出用任意连续的尺度函数构造区间上的插值基函数,形成以尺度函数为基础的求解两点边值问题的小波配点法.该方法中,尺度函数不受紧支撑、插值等性质的限制,计算复杂度小,数值解收敛性由多分辨分析理论保证.同时,给出边值条件的积分处理方法,能够方便地处理任意边界条件,当尺度函数不具有高阶导数时,该方法也能有效使用.数值算例表明,该方法是一个高效、高精度的算法.