On Mathematical Modeling and Simulation of the Pressing Section of a Paper Machine Including Dynamic Capillary Effects: One-Dimensional Model

This study presents the dynamic capillary pressure model (Hassanizadeh and Gray, Adv Water Resour 13:169–186, 1990; Water Resour Res 29:3389–3405, 1993) adapted for the needs of paper manufacturing process simulations. The dynamic capillary pressure–saturation relation is included in a one-dimensional simulation model for the pressing section of a paper machine. The one-dimensional model is derived from a two-dimensional model by averaging with respect to the vertical direction. Then, the model is discretized by the finite volume method and solved by Newton’s method. The numerical experiments are carried out for parameters typical for the paper layer. The dynamic capillary pressure–saturation relation shows significant influence on the distribution of water pressure. The behavior of the solution agrees with laboratory experiments (Beck, Fluid pressure in a press nip: measurements and conclusions, 1983).     

Calculation of stagnation streamline quantities in hypersonic blunt body flows

An approximate method for the efficient calculation of stagnation-streamline quantities in hypersonic flows about spheres or cylinders is suggested. Based on the local similarity of the flow field the two-dimensional Navier-Stokes equations are simplified to a one-dimensional approximation for the stagnation streamline. These equations are solved with an implicit finite-volume scheme. Comparisons with fully two–dimensional Euler and Navier–Stokes calculations for flows about spheres are presented, that include perfect gas flows and flows in chemical non-equilibrium. Comparisons with a number of experiments conclude this report. Received 8 May 1996 / Accepted 31 October 1996     

Global Solvability in W 2 2,1 -Weighted Spaces of the Two-Dimensional Navier–Stokes Problem in Domains with Strip-Like Outlets to Infinity

The time-dependent Navier–Stokes system is studied in a two-dimensional domain with strip-like outlets to infinity in weighted Sobolev function spaces. It is proved that under natural compatibility conditions there exists a unique solution with prescribed fluxes over cross-sections of outlets to infinity which tends in each outlet to the corresponding time-dependent Poiseuille flow. The obtained results are proved for arbitrary large norms of the data (in particular, for arbitrary fluxes) and globally in time. The authors are supported by EC FP6 MC–ToK programme SPADE2, MTKD–CT–2004–014508.     

Second-Order Fredholm Equations for the First Boundary-Value Problem in the Two-Dimensional Anisotropic Theory of Elasticity

A complete potential theory is constructed for the first boundary-value problem in the two-dimensional anisotropic theory of elasticity (the force vector is specified on the boundary) in a bounded domain on a plane with a Lyapunov boundary. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 2, pp. 85–94, March–April, 2006.     

Stress-strain analysis of closed nonthin orthotropic conical shells of varying thickness

The approach developed to solve two-dimensional static problems for nonthin conical shells of varying thickness is used to examine the effect of the geometrical parameters on the stress-strain state of shells. The approach is based on spline-approximation and a stable numerical method of solving one-dimensional problems __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 6, pp. 46–58, June 2008.     

Stabilization and evolution of parameters of a symmetric laminar flow in a plane channel with sudden expansion

A problem is formulated for computing the fields of parameters of a stationary laminar symmetric flow. A two-dimensional flow in a channel with a sudden change in the cross-sectional area is computed. The evolution of a three-dimensional perturbation inserted into the channel at the initial stage of computations is analyzed. It is demonstrated that the parameters of a two-dimensional flow in the channel at a Reynolds number Re = 50 become stabilized at a dimensionless time t > 20, whereas the steady state is reached under the same conditions at t ≈ 100. At a distance of approximately 10h (h is the channel width at the entrance), the flow becomes one-dimensional, but the streamwise component of the velocity vector remains a function of the streamwise coordinate owing to flow compressibility. __________ Translated from PrikladnayaMekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 35–42, January–February, 2007.     

On a Discrete Fourier Series Solution to Static Problem for Conical Shells of Circumferentially Varying Thickness

An approach is developed to solve the two-dimensional boundary-value problems of the stress-strain state of conical shells with circumferentially varying thickness. The approach employs discrete Fourier series to separate variables and make the problem one-dimensional. The one-dimensional boundary-value problem is solved by the stable discrete-orthogonalization method. The results obtained are presented as plots and tables __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 9, pp. 26–37, September 2005.     

Numerical solution of solidification in a square prism using an algebraic grid generation technique

The solidification of an infinitely long square prism was analyzed numerically. A front fixing technique along with an algebraic grid generation scheme was used, where the finite difference form of the energy equation is solved for the temperature distribution in the solid phase and the solid–liquid interface energy balance is integrated for the new position of the moving solidification front. Results are given for the moving solidification boundary with a circular phase change interface. An algebraic grid generation scheme was developed for two-dimensional domains, which generates grid points separated by equal distances in the physical domain. The current scheme also allows the implementation of a finer grid structure at desired locations in the domain. The method is based on fitting a constant arc length mesh in the two computational directions in the physical domain. The resulting simultaneous, nonlinear algebraic equations for the grid locations are solved using the Newton-Raphson method for a system of equations. The approach is used in a two-dimensional solidification problem, in which the liquid phase is initially at the melting temperature, solved by using a front-fixing approach. The difference of the current study lies in the fact that front fixing is applied to problems, where the solid–liquid interface is curved such that the position of the interface, when expressed in terms of one of the coordinates is a double valued function. This requires a coordinate transformation in both coordinate directions to transform the complex physical solidification domain to a Cartesian, square computational domain. Due to the motion of the solid–liquid interface in time, the computational grid structure is regenerated at every time step.     

Refined design of longitudinally corrugated cylindrical shells

A refined Timoshenko-type model based on the straight-line hypothesis is used to develop an approach to analyzing the stress state of longitudinally corrugated cylindrical shells with elliptic cross-section. The approach is to reduce the two-dimensional boundary-value problem that describes the stress–strain state of the shell to a one-dimensional one and to solve it with the stable numerical discrete-orthogonalization method. The solutions obtained using the straight-line hypothesis and the equations of three-dimensional elasticity are compared. The dependence of the stress–strain state of the shell on the number and amplitude of corrugations and the aspect ratio of the cross-section is analyzed     

2D Navier–Stokes Equations in a Time Dependent Domain with Neumann Type Boundary Conditions

In this paper the two-dimensional Navier–Stokes system for incompressible fluid coupled with a parabolic equation through the Neumann type boundary condition for the second component of the velocity is considered. Navier–Stokes equations are defined on a given time dependent domain. We prove the existence of a weak solution for this system. In addition, we prove the continuous dependence of solutions on the data for a regularized version of this system. For a special case of this regularized system also a problem with an unknown interface is solved.     

Domain decomposition method with hybrid approximations applied to solve problems of elasticity

To solve two-dimensional boundary-value problems of elasticity, two iteration algorithms of the domain decomposition method are proposed: parallel Neumann–Neumann and sequential Dirichlet–Neumann. They are based on the hybrid boundary–finite-element approximations. The algorithms are proved to converge. The optimal parameters are selected using the minimum-residual and steepest-descent methods. Some plane problems of elasticity are solved as examples, and stationary and nonstationary iteration algorithms in these examples are analyzed for efficiency Translated from Prikladnaya Mekhanika, Vol. 44, No. 11, pp. 18–29, November 2008.     

Solving the torsion problem for an anisotropic prism by the advanced Kantorovich–Vlasov method

The torsion problem for a rectangular prism with general anisotropy loaded on the lateral surface is solved using the advanced Kantorovich–Vlasov method, which reduces the original three-dimensional problem to three coupled one-dimensional problems, each for one of the variables of the domain. The warping of the cross-section and the deformation of the axis of the prism for different types of anisotropy are analyzed     

Spline-approximation solution of two-dimensional problems of statics for orthotropic conical shells in a refined formulation

An approach to solving two-dimensional stress-strain problems for orthotropic conical shells of variable thickness in a refined formulation is developed. The approach is based on the spline-approximation and the stable discrete-orthogonalization method, which is used to solve the one-dimensional problem. The dependence of the deflection and stresses on the thickness and apex angle of a shell of constant volume is given as an example __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 11, pp. 43–54, November 2007.     

Determining modes and Grashof number in 2D turbulence: a numerical case study

We study how the number of numerically determining modes in the Navier–Stokes equations depends on the Grashof number. Consider the two-dimensional incompressible Navier–Stokes equations in a periodic domain with a fixed time-independent forcing function. We increase the Grashof number by rescaling the forcing and observe through numerical computation that the number of numerically determining modes stabilizes at some finite value as the Grashof number increases. This unexpected result implies that our theoretical understanding of continuous data assimilation is incomplete until an analytic proof which makes use of the non-linear term in the Navier–Stokes equations is found.      

A return toward equilibrium in a 2D Rayleigh–Taylor instability for compressible fluids with a multidomain adaptive Chebyshev method

We numerically simulate a single-mode Rayleigh–Taylor instability between compressible miscible fluids with a highly accurate self-adaptive pseudospectral Chebyshev multidomain method in two two-dimensional boxes at small aspect ratios. The simulations are started from rest and pursued until the return toward mechanical equilibrium of the mixing. Four regimes—linear and weakly nonlinear, nonlinear steady bubble rise, return toward equilibrium, and finally a system of acoustic waves—can be identified. We show that this one-dimensional system of stationary acoustic waves is damped by the physical viscosity. This provides a reference solution.      

Spline-approximation solution of problems of the statics of orthotropic shallow shells with variable parameters

An approach to the solution of problems of the statics of shallow orthotropic shells is proposed. It is based on reducing a two-dimensional boundary value problem to a one-dimensional one using the spline-collocation method and solution of the problem by the stable numerical method of discrete orthogonalization. Solutions are presented for problems on the stress state of orthotropic shells of double curvature for several values of the elastic constants of the material. Translated from Prikladnaya Mekhanika, Vol. 36, No. 7, pp. 60–66, July, 2000.     

Another Proof of Stability for Global Weak Solutions of 2D Degenerated Shallow Water Models

We consider non-linear viscous shallow water models with varying topography, extra friction terms and capillary effects, in a two-dimensional framework. Water-depth dependent laminar and turbulent friction coefficients issued from an asymptotic analysis of the three-dimensional free-surface Navier–Stokes equations are considered here. A new proof of stability for global weak solutions is given in periodic domain Ω = T², adapting the method introduced by J. Simon in [15] for the non-homogeneous Navier–Stokes equations. Existence results for such solutions can be obtained from this stability analysis.     

Static analysis of two-dimensional elastic structures using global collocation

Based on previous findings concerning the numerical solution of one-dimensional elastodynamical problems [Provatidis in Arch Appl Mech 78(4):241–250, 2008] this paper extends the methodology to the static analysis of two-dimensional problems in quadrilateral domains. This target is achieved by replacing the Galerkin/Ritz procedure involved in Lagrangian (or Gordon–Coons) type finite elements by a global collocation scheme. In brief, the boundary conditions are fulfilled at all boundary nodes, while the governing equation is fulfilled at internal points. The theory is supported by four test cases concerning rectangular and curvilinear structures under plane-stress or plane-strain conditions, where the convergence rate is successfully compared with that of conventional bilinear finite elements with the same mesh density.     

Ligament formation in sheared liquid–gas layers

We perform numerical simulations of two-phase liquid–gas sheared layers, with the objective of studying atomization. The Navier–Stokes equations for two-dimensional incompressible flow are solved in a periodic domain. A volume-of-fluid method is used to track the interface. The density ratio is kept around 10. The calculations show good agreement with a fully viscous Orr–Sommerfeld linear theory over several orders of magnitude of interface growth. The nonlinear development shows the growth of finger-like structures, or ligaments, and the detachment of droplets. The effect of the Weber and Reynolds numbers, the boundary layer width and the initial perturbation amplitude are discussed through a number of typical cases. Inversion of the liquid boundary layer is shown to yield more readily ligaments bending upwards and is thus more likely to produce droplets.     

Optimal Results for the Nonhomogeneous Initial-Boundary Value Problem for the Two-Dimensional Navier–Stokes Equations

In this work we study the fully nonhomogeneous initial boundary value problem for the two-dimensional time-dependent Navier–Stokes equations in a general open space domain in R² with low regularity assumptions on the initial and the boundary value data. We show that the perturbed Navier–Stokes operator is a diffeomorphism from a suitable function space onto its own dual and as a corollary we get that the Navier–Stokes equations are uniquely solvable in these spaces and that the solution depends smoothly on all involved data. Our source data space and solution space are in complete natural duality and in this sense, without any smallness assumptions on the data, we solve the equations for data with optimally low regularity in both space and time.