Sheath vaporization of a monodisperse fuel-spray jet

The group vaporization of a monodisperse fuel-spray jet discharging into a hot coflowing gaseous stream is investigated for steady flow by numerical and asymptotic methods with a two-continua formulation used for the description of the gas and liquid phases. The jet is assumed to be slender and laminar, as occurs when the Reynolds number is moderately large, so that the boundary-layer form of the conservation equations can be employed in the analysis. Two dimensionless parameters are found to control the flow structure, namely the spray dilution parameter 1, defined as the mass of liquid fuel per unit mass of gas in the spray stream, and the group vaporization parameter e, defined as the ratio of the characteristic time of spray evolution due to droplet vaporization to the characteristic diffusion time across the jet. It is observed that, for the small values of e often encountered in applications, vaporization occurs only in a thin layer separating the spray from the outer droplet-free stream. This regime of sheath vaporization, which is controlled by heat conduction, is amenable to a simplified asymptotic description, independent of ε,in which the location of the vaporization layer is determined numerically as a free boundary in a parabolic problem involving matching of the separate solutions in the external streams, with appropriate jump conditions obtained from analysis of the quasi-steady vaporization front. Separate consideration of dilute and dense sprays, corresponding, respectively, to the asymptotic limits λ<<1 and λ>>1, enables simplified descriptions to be obtained for the different flow variables, including explicit analytic expressions for the spray penetration distance.

Asymptotic stability of a mathematical model of cell population

We consider a simplified system of a growing colony of cells described as a free boundary problem. The system consists of two hyperbolic equations of first order coupled to an ODE to describe the behavior of the boundary. The system for cell populations includes non-local terms of integral type in the coefficients. By introducing a comparison with solutions of an ODE's system, we show that there exists a unique homogeneous steady state which is globally asymptotically stable for a range of parameters under the assumption of radially symmetric initial data.

A gradient-enhanced large-deformation continuum damage model for fibre-reinforced materials

A non-local gradient-based damage formulation within a geometrically non-linear setting is presented. The hyperelastic constitutive response at local material point level is governed by a strain energy which is additively composed of an isotropic matrix and of an anisotropic fibre-reinforced material, respectively. The inelastic constitutive response is governed by a scalar [1–d]-type damage formulation, where only the anisotropic elastic part is assumed to be affected by the damage. Following the concept in Dimitrijević and Hackl [28], the local free energy function is enhanced by a gradient-term. This term essentially contains the gradient of the non-local damage variable which, itself, is introduced as an additional independent variable. In order to guarantee the equivalence between the local and non-local damage variable, a penalisation term is incorporated within the free energy function. Based on the principle of minimum total potential energy, a coupled system of Euler–Lagrange equations, i.e., the balance of linear momentum and the balance of the non-local damage field, is obtained and solved in weak form. The resulting coupled, highly non-linear system of equations is symmetric and can conveniently be solved by a standard incremental-iterative Newton–Raphson-type solution scheme. Several three-dimensional displacement- and force-driven boundary value problems—partially motivated by biomechanical application—highlight the mesh-objective characteristics and constitutive properties of the model and illustratively underline the capabilities of the formulation proposed

Some aspects of lubrication in heavy regimes: thermal effects, stability and turbulence

In this work, a combination of numerical methods applied to thermohydrodynamic lubrication problems with cavitation is presented. It should be emphasized the difficulty of the nonlinear mathematical coupled model involving a free boundary problem, but also the simplicity of the algorithms employed to solve it. So, finite element discretizations for the hydrodynamic and thermal equations combined with upwind techniques for the convection terms and duality methods for nonlinear features are proposed. Additionally, a model describing the movement of the shaft is provided. Considering the shaft as a rigid body this model will consist of an ODE system relating acceleration of the center of gravity and external and pressure loads. The numerical experiments of mechanical stability try to clarify the position of the neutral stability curve. Finally, a rotating machine for ship propulsion involving both axial and radial bearings operating with nonconventional lubricants (seawater to avoid environmental pollution) is analyzed by using laminar and turbulent inertial flows.

Application of the Boundary Method to the determination of the properties of the beam cross-sections

Using the 3-D equations of linear elasticity and the asylllptotic expansion methods in terms of powers of the beam cross-section area as small parameter different beam theories can be obtained, according to the last term kept in the expansion. If it is used only the first two terms of the asymptotic expansion the classical beam theories can be recovered without resort to any "a priori" additional hypotheses. Moreover, some small corrections and extensions of the classical beam theories can be found and also there exists the possibility to use the asymptotic general beam theory as a basis procedure for a straightforward derivation of the stiffness matrix and the equivalent nodal forces of the beam. In order to obtain the above results a set of functions and constants only dependent on the cross-section of the beam it has to be computed them as solutions of different 2-D laplacian boundary value problems over the beam cross section domain. In this paper two main numerical procedures to solve these boundary value pf'oblems have been discussed, namely the Boundary Element Method (BEM) and the Finite Element Method (FEM). Results for some regular and geometrically simple cross-sections are presented and compared with ones computed analytically. Extensions to other arbitrary cross-sections are illustrated.

Application of the homogenization techniques to the determination of the effective flexure constants of plate structural systems

A Mindlin plate with periodically distributed ribs patterns is analyzed by using homogenization techniques based on asymptotic expansion methods. The stiffness matrix of the homogenized plate is found to be dependent on the geometrical characteristics of the periodical cell, i.e. its skewness, plan shape, thickness variation etc. and on the plate material elastic constants. The computation of this plate stiffness matrix is carried out by averaging over the cell domain some solutions of different periodical boundary value problems. These boundary value problems are defined in variational form by linear first order differential operators on the cell domain and the boundary conditions of the variational equation correspond to a periodic structural problem. The elements of the stiffness matrix of homogenized plate are obtained by linear combinations of the averaged solution functions of the above mentioned boundary value problems. Finally, an illustrative example of application of this homogenization technique to hollowed plates and plate structures with ribs patterns regularly arranged over its area is shown. The possibility of using in the profesional practice the present procedure to the actual analysis of floors of typical buildings is also emphasized.

A gradient estimate to a degenerate parabolic equation with a singular absorption term: The global quenching phenomena

We prove global existence of nonnegative solutions to the one dimensional degenerate parabolic problems containing a singular term. We also show the global quenching phenomena for L1 initial datums. Moreover, the free boundary problem is considered in this paper.

Positive solutions for nonlinear m-point Eigenvalue problems

In this paper, we are concerned with determining values of lambda, for which there exist positive solutions of the nonlinear eigenvalue problem [GRAPHICS] where a, b, c, d is an element of [0, infinity), xi(i) is an element of (0, 1), alpha(i), beta(i) is an element of [0 infinity) (for i is an element of {1, ..., m - 2}) are given constants, p, q is an element of C ([0, 1], (0, infinity)), h is an element of C ([0, 1], [0, infinity)), and f is an element of C ([0, infinity), [0, infinity)) satisfying some suitable conditions. Our proofs are based on Guo-Krasnoselskii fixed point theorem. (C) 2004 Elsevier Inc. All rights reserved.

Deflection of jets discharged into a reservoir with a free surface

Deflections of jets discharged into a reservoir with a free surface are investigated numerically. The jets are known to deflect towards either side of the free surface or the bottom, whose direction is not determined uniquely in some experimental conditions, i.e. there are multiple stable states realizable in the same condition. The origin of the multiple stable states is explored by utilizing homotopy transformations in which the top boundary of the reservoir is transformed from a rigid to a free boundary and also the location of the outlet throat is continuously moved from mid-height to the top. We depicted bifurcation diagrams of the flow compiling the data of numerical simulations, from which we identified the origin as an imperfect pitchfork bifurcation, and obtained an insight into the mechanism for the direction to be determined. The parameter region where such multiple stable states are possible is also delimited. © 2011 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.

Determining the temperature from Cauchy data in corner domains

We consider the problem of reconstruction of the temperature from knowledge of the temperature and heat flux on a part of the boundary of a bounded planar domain containing corner points. An iterative method is proposed involving the solution of mixed boundary value problems for the heat equation (with time-dependent conductivity). These mixed problems are shown to be well-posed in a weighted Sobolev space.

A procedure for the reconstruction of a stochastic stationary temperature field

An iterative procedure is proposed for the reconstruction of a temperature field from a linear stationary heat equation with stochastic coefficients, and stochastic Cauchy data given on a part of the boundary of a bounded domain. In each step, a series of mixed well-posed boundary-value problems are solved for the stochastic heat operator and its adjoint. Well-posedness of these problems is shown to hold and convergence in the mean of the procedure is proved. A discretized version of this procedure, based on a Monte Carlo Galerkin finite-element method, suitable for numerical implementation is discussed. It is demonstrated that the solution to the discretized problem converges to the continuous as the mesh size tends to zero.

A variational conjugate gradient method for determining the fluid velocity of a slow viscous flow

The problem considered is that of determining the fluid velocity for linear hydrostatics Stokes flow of slow viscous fluids from measured velocity and fluid stress force on a part of the boundary of a bounded domain. A variational conjugate gradient iterative procedure is proposed based on solving a series of mixed well-posed boundary value problems for the Stokes operator and its adjoint. In order to stabilize the Cauchy problem, the iterations are ceased according to an optimal order discrepancy principle stopping criterion. Numerical results obtained using the boundary element method confirm that the procedure produces a convergent and stable numerical solution.

On the Stabilization of the Wave Equation by the Boundary

* Partially supported by CNPq (Brazil)

Singular Solutions with Exponential Growth for the (3+1)-D Wave Equation

Недю И. Попиванов, Тодор П. Попов, Рудолф Шерер - Разглеждат се четиримерни гранични задачи за нехомогенното вълново уравнение. Те са предложени от М. Протер като многомерни аналози на задачата на Дарбу в равнината. Известно е, че единственото обобщено решение може да има силна степенна особеност само в една гранична точка. Тази сингулярност е изолирана във върха на характеристичния конус и не се разпространява по конуса. Друг аспект на проблема е, че задачата не е фредхолмова, тъй като има безкрайномерно коядро. Предишни резултати сочат, че решението може да има най-много експоненциален ръст, но оставят открит въпроса дали наистина съществуват такива решения. Показваме, че отговора на този въпрос е положителен и строим обобщено решение на задачата на Протер с експоноциална особеност.

Interior Boundaries for Degenerate Elliptic Equations of Second Order Some Theory and Numerical Observations

2010 Mathematics Subject Classification: Primary 35J70; Secondary 35J15, 35D05.