Modeling of multicomponent-fuel drop-laden mixing layers having a multitude of species

A formulation representing multicomponent-fuel (MC-fuel) composition as a probability distribution function (PDF) depending on the molar mass is used to construct a model of a large number of MC-fuel drops evaporating in a gas flow, so as to assess the extent of fuel specificity on the vapor composition. The PDF is a combination of two Gamma PDFs, which was previously shown to duplicate the behavior of a fuel composed of many species during single drop evaporation. The conservation equations are Eulerian for the flow and Lagrangian for the physical drops, all of which are individually followed. The gas conservation equations for mass, momentum, species, and energy are complemented by differential conservation equations for the first four moments of the gas-composition PDF; all coupled to the perfect gas equation of state. Source terms in all conservation equations couple the gas phase to the drops. The drop conservation equations for mass, position, momentum, and energy are complemented by differential equations for four moments of the liquid-composition PDF. The simulations are performed for a three-dimensional mixing layer whose lower stream is initially laden with drops. Initial perturbations excite the layer to promote the double pairing of its four initial spanwise vortices to an ultimate vortex. The drop temperature is initially lower than that of the surrounding gas, initiating drop heating and evaporation. The results focus on both evolution and the state of the drops and gas when layers reach a momentum-thickness maximum past the double vortex pairing; particular emphasis is on the gas composition. Comparisons between simulations with n-decane, diesel, and three kerosenes show that at same initial Reynolds number and Stokes number distribution, a single-component fuel cannot represent MC fuels. Substantial differences among the MC-fuel vapor composition indicate that fuel specificity must be captured for the prediction of combustion.     

MC-Smooth: A mass-conserving,smooth vorticity–velocity formulation for multi-dimensional flows

A novel vorticity–velocity formulation of the Navier–Stokes equations – the Mass-Conserving, Smooth (MC-Smooth) vorticity–velocity formulation – is developed in this work. The governing equations of the MC-Smooth formulation include a new second-order Poisson-like elliptic velocity equation, along with the vorticity transport equation, the energy conservation equation, and N_spec species mass balance equations. In this study, the MC-Smooth formulation is compared to two pre-existing vorticity–velocity formulations by applying each formulation to confined and unconfined axisymmetric laminar diffusion flame problems. For both applications, very good to excellent agreement for the simulation results of the three formulations has been obtained. The MC-Smooth formulation requires the least CPU time and can overcome the limitations of the other two pre-existing vorticity–velocity formulations by ensuring mass conservation and solution smoothness over a broader range of flow conditions. In addition to these benefits, other important features of the MC-Smooth formulation include: (1) it does not require the use of a staggered grid, and (2) it does not require excessive grid refinement to ensure mass conservation. The MC-Smooth formulation is a computationally attractive approach that can effectively extend the applicability of the vorticity–velocity formulation.     

Optical filaments and optical bullets in dispersive nonlinear media

We have investigated the evolution of picosecond and femtosecond optical pulses governed by the amplitude vector equation in the optical and UV domains. We have written this equation in different coordinate frames, namely, in the laboratory frame, the Galilean frame, and the moving-in-time frame and have normalized it for the cases of different and equal transverse and longitudinal sizes of optical pulses or modulated optical waves. For optical pulses with a small transverse size and a large longitudinal size (optical filaments), we obtain the well-known paraxial approximation in all the coordinate frames, while for optical pulses with relatively equal transverse and longitudinal sizes (so-called light bullets), we obtain new non-paraxial nonlinear amplitude equations. In the case of optical fields with low intensity, we have reduced the nonlinear amplitude vector equations governing the light-bullet evolution to the linear amplitude equations. We have solved the linear equations using the method of Fourier transform. An unexpected new result is the relative stability of light bullets and the significant decrease in the diffraction enlargement of light bullets with respect to the case of long pulses in the linear propagation regime.     

Divergence formulation of source term

In this paper, we propose to write a source term in the divergence form. A conservation law with a source term can then be written as a single divergence form. We demonstrate that it enables to discretize both the conservation law and the source term in the same framework, and thus greatly simplifies the construction of numerical schemes. To illustrate the advantage of the divergence formulation, we apply the new formulation to construct a uniformly third-order accurate edge-based finite-volume scheme for conservation laws with a source term. Third-order accuracy is demonstrated for regular and irregular triangular grids for the linear advection and Burgers’ equations with a source term.     

积分梯度法的两种格式

在二维柱坐标系下Lagrange流体力学的计算中,积分梯度法是动量方程的一种有效离散方法.积分梯度法中,IGT(Integral Gradient Total)格式不能保持柱几何下一维球对称性;IGA(Integral Gradient Average)格式可以保持一维球对称性,但当相邻网格质量相差比较大时,会得到远远脱离真实物理现象的加速度.深入研究IGA和IGT格式发现,当相邻网格边界压力取为质量加权时,即使相邻网格质量相差较大,对于一维平面和一维柱问题,IGT与IGA等价;在二维情形下,可以缩小IGT和IGA之间的差异.理论证明,IGA格式不能保持系统的动量守恒,IGT格式能保持系统的动量守恒性.数值模拟结果进一步显示了这两个格式的优缺点.     

A stable high-order Spectral Difference method for hyperbolic conservation laws on triangular elements

Numerical schemes using piecewise polynomial approximation are very popular for high order discretization of conservation laws. While the most widely used numerical scheme under this paradigm appears to be the Discontinuous Galerkin method, the Spectral Difference scheme has often been found attractive as well, because of its simplicity of formulation and implementation. However, recently it has been shown that the scheme is not linearly stable on triangles. In this paper we present an alternate formulation of the scheme, featuring a new flux interpolation technique using Raviart–Thomas spaces, which proves stable under a similar linear analysis in which the standard scheme failed. We demonstrate viability of the concept by showing linear stability both in the semi-discrete sense and for time stepping schemes of the SSP Runge–Kutta type. Furthermore, we present convergence studies, as well as case studies in compressible flow simulation using the Euler equations.     

On the correct representation of bending and axial deformation in the absolute nodal coordinate formulation with an elastic line approach

In finite element methods that are based on position and slope coordinates, a representation of axial and bending deformation by means of an elastic line approach has become popular. Such beam and plate formulations based on the so-called absolute nodal coordinate formulation have not yet been verified sufficiently enough with respect to analytical results or classical nonlinear rod theories. Examining the existing planar absolute nodal coordinate element, which uses a curvature proportional bending strain expression, it turns out that the deformation does not fully agree with the solution of the geometrically exact theory and, even more serious, the normal force is incorrect. A correction based on the classical ideas of the extensible elastica and geometrically exact theories is applied and a consistent strain energy and bending moment relations are derived. The strain energy of the solid finite element formulation of the absolute nodal coordinate beam is based on the St. Venant-Kirchhoff material: therefore, the strain energy is derived for the latter case and compared to classical nonlinear rod theories. The error in the original absolute nodal coordinate formulation is documented by numerical examples. The numerical example of a large deformation cantilever beam shows that the normal force is incorrect when using the previous approach, while a perfect agreement between the absolute nodal coordinate formulation and the extensible elastica can be gained when applying the proposed modifications. The numerical examples show a very good agreement of reference analytical and numerical solutions with the solutions of the proposed beam formulation for the case of large deformation pre-curved static and dynamic problems, including buckling and eigenvalue analysis. The resulting beam formulation does not employ rotational degrees of freedom and therefore has advantages compared to classical beam elements regarding energy-momentum conservation.     

Wave based analysis of the Green's function for a layered cylindrical shell

Cylindrical shells composed of concentric layers may be designed to affect the way that elastic waves are generated and propagated, particularly when some layers are anisotropic. To aid the design process, the present work develops a wave based analysis of the Green's function for a layered cylindrical shell in which the response is given as a sum of waves propagating in the axial coordinate. The analysis assumes linear Hookean materials for each layer. It uses finite element discretizations in the radial coordinate and Fourier series expansions in the circumferential coordinate, leading to linear equations in the axial wavenumber domain that relate shell displacements and forces. Inversion to the axial domain is accomplished via a state-space formulation that is evaluated using residue integration. The resulting expression for the Green's function for each circumferential harmonic is a summation over the natural waves of the shell. The finite element discretization in the radial direction allows the approach to be used for arbitrarily thick shells. The approach is benchmarked to results from an isotropic shell and numerical examples are given for a shell composed of a fiber-reinforced material. The numerical examples illustrate the effect of fiber orientation on the Green's function.     

Surface instability on thin fluid layers of a binary mixture

A long-wave model for thin layers consisting of two miscible fluids is presented. The model is a development of a simplified 2D model variant which considers the temperature and the concentration fields as linear functions on the vertical coordinate and neglect the convective terms from the corresponding equations. Now, the 2D thin film equation is coupled to complete 3D energy and mass conservations equation. We discuss the extended system in the linear approximation and in the initial nonlinear stage.     

Mode selection in weakly unstable two-dimensional detonations

A formulation of the reactive Euler equations in the shock-attached frame is used to study the two-dimensional instability of weakly unstable detonation through direct numerical simulation. The results are shown to agree with the predictions of linear stability analysis. Comparisons are made with linear perturbation growth rates and oscillation frequencies as a function of transverse disturbance wavelength. The perturbation eigenfunctions predicted by linear stability analysis are directly validated through numerical simulation. Three regimes of unstable behavior – linear, weakly nonlinear, and fully nonlinear – are explored and characterized in terms of the power spectrum of the normal shock velocity for a Chapman–Jouguet detonation with weak heat release.     

Multi-dimensional conservative semi-Lagrangian method of characteristics CIP for the shallow water equations

A new characteristic approach that guarantees conservative property is proposed and is applied to the shallow water equations. CIP–CSL (Constrained Interpolation Profile/Conservative Semi-Lagrangian) interpolation is applied to the CIP method of characteristics in order to enhance the mass conservation of the numerical result. Although the characteristic formulation is originally derived from non-conservative form, present scheme achieves complete mass conservation by solving mass conservation simultaneously and reflecting conserving mass in interpolation profile. Present method has less height error compared to the CIP method of characteristics by several orders of magnitude. By the enhanced conservation property, present scheme is applicable to nonlinear problem such as shock. Furthermore, application to two dimensions including the Coriolis term is straightforward with directional splitting technique.     

Modes of a rotating dielectric waveguide

Proceeding from the Maxwell differential equations and the Minkowski constitutive equations of the first order with respect to the ratio of the velocity of motion of a medium to the velocity of light in a vacuum, we derived equations for determining the mode structure of a field in a linear isotropic medium with an axially symmetric distribution of the permittivity and the permeability and with the angular velocity of rotation dependent on the distance from the symmetry axis. A waveguide effect was shown to appear even for a homogeneous medium owing to the spatial nonuniformity of the angular velocity. In the case of the quadratic dependence of the angular velocity on the transverse coordinate, the exact solutions for the modes were found in the form of Gaussian beams.     

Conservation of mass and momentum of the least-squares spectral collocation scheme for the Stokes problem

From the literature it is known that spectral least-squares schemes perform poorly with respect to mass conservation and compensate this lack by a superior conservation of momentum. This should be revised, since the here presented new least-squares spectral collocation scheme leads to an outstanding performance with respect to conservation of momentum and mass. The reasons can be found in using only a few elements, each with high polynomial degree, avoiding normal equations for solving the overdetermined linear systems of equations and by introducing the Clenshaw–Curtis quadrature rule for imposing the average pressure to be zero. Furthermore, we combined the transformation of Gordon and Hall (transfinite mapping) with our least-squares spectral collocation scheme to discretize the internal flow problems.     

A tapered beam finite element for rotor dynamics analysis

The stiffness, mass and gyroscopic matrices of a rotating beam element are developed, a cubic function being used for the transverse displacement. Shear deflection is included by use of end nodal variables of shear strain, along with transverse displacement and cross-section rotation; rotatory inertia effects are included in the energy functional to provide a Timoshenko beam formulation. The gyroscopic effects for small perturbations are linearized as a skew symmetric damping matrix. The formulation is implemented by numerical integration for a linearly tapered circular beam. A technique of reduction of the shear nodal variable prior to global assembly is shown to provide little loss in accuracy with reduced system bandwidth. Numerical comparisons for three previously published beam models are included, with results presented for the case of forward and reverse precession to verify the gyroscopic effects. The utility of the element in a general program for rotor dynamics analysis is identified.     

Spin,torsion, forces

A theory of gravitation with torsion that is derived from a potential is developed. An explicit material action is presented that gives rise to the correct conservation laws and equations of motion. It is shown that the Noether identities yield the same conservation laws as the Bianchi identities, and the Papapetrou method is used to develop the propagation equations and the force law. The equations are put in the nonrelativistic limit and the 3-vector formulation is displayed.     

A nonlinear PSE method for two-fluid shear flows with complex interfacial topology

A nonlinear stability method is developed for laminar two-fluid shear flows which undergo changes in the interface topology. The method is based on the nonlinear parabolized stability equations (PSE) and incorporates a scalar-based interface capturing (IC) scheme in order to track complex deformations of the fluid interface. In doing so, the formulation retains the flexibility and physical insight of instability-wave based methods, while providing hydrodynamic modeling capabilities similar to direct numerical calculations: the new formulation, referred to as the IC-PSE, can capture the nonlinear physical mechanisms responsible for generating large-scale, two-fluid structures, without incurring heavy computational costs. This approach is valid for spatially developing, laminar two-fluid shear flows which are convectively unstable, and can naturally account for the growth of finite amplitude interfacial waves, along with changes to the interfacial topology. We demonstrate the accuracy of the IC-PSE against direct Navier–Stokes calculations for two-fluid mixing layers with density and viscosity stratification. The comparisons show that the IC-PSE can predict the dynamics of the instability waves and capture the formation of Kelvin–Helmholtz vortex rolls and large scale liquid structures, at an order of magnitude less computational cost than direct calculations. The role of surface tension in the IC-PSE formulation is shown to be valid for flows in which Re/We ? 1, and the method accurately predicts the formation and non-linear evolution of flow structures in this limit. This is demonstrated for spatially developing mixing layers which lead to vortex roll-up and ligaments, prior to droplet formation. The pinch-off process itself is a high surface tension phenomenon and in not considered herein. The method also accurately captures the effect of interfacial waves on the mean flow, and the topology changes during the non-linear evolution of the two-fluid structures.     

The effect of an inhomogeneous cladding of a gradient optical waveguide on the mode and envelope characteristics of a soliton pulse

The nonlinear dynamics of a short pulse in a gradient waveguide layer is studied analytically taking into account the occurrence of a cladding layer and the longitudinal inhomogeneity of the two layers. The class of functions is presented in terms of which the transverse profiles of the refractive indices of both the gradient and the cladding layers can be represented, and, at the same time, the mode structure of the pulse can be explicitly described in terms of a hypergeometric equation. Expressions are presented that describe the variations of the propagation constant and transverse distribution of the wave field under the action of the longitudinal inhomogeneities of the gradient and cladding layers. It is shown that the envelope of the pulse satisfies the generalized nonlinear Schrödinger equation the coefficients of which are functions of the longitudinal coordinate and are expressed via the refractive indices of the waveguide layer and cladding.     

Proper-Time Formulation of Relativistic Dynamics

It will be argued that Minkowski's implementation of distances is inconsistent. An alternative implementation will be proposed. In the new model the proper time of an object is taken as its fourth coordinate. Distances will be measured according to a four dimensional Euclidean metric. In the present approach mass is a constant of motion. A mass can therefore be ascribed to photons and neutrinos. Mechanics and dynamics will be reformulated in close correspondence with classical physics. Of particular interest is the equation of motion for the proper time momentum. In the classical limit it reduces to the classical law of conservation of (kinetic+potential) energy. In the relativistic limit it is similar to the conservation of energy of the theory of relativity. The conservation of proper time momentum allows for an alternative explanation for Compton scattering and pair annihilation. On the basis of the proper time formulation of electrodynamics also an alternative explanation will be offered for the spectra of hydrogenic atoms. The proper time formulation of gravitational dynamics leads to the correct predictions of gravitational time dilation, the deflection of light and the precession of the perihelia of planets. For this no curvature will be needed. That is, spacetime is flat everywhere, even in the presence of sources of gravitation. Some cosmological consequences will be discussed. The present approach gives a new notion to energy, antiparticles and the structure of spacetime. The contents of the present paper will have important implications for the foundations of physics in general.     

The photon momentum

The momentum expression of a photon is derived by taking into account the absorption process by an electron. The defect of the actual momentum representation is shown by using the equations of energy and momentum conservation. The phase and group velocities of the photon that represents the particle aspect of light are also obtained with the aid of the new formulation.     

Propagation of surface elastic waves in the Cosserat medium

The problem of surface elastic wave propagation in the Cosserat medium (half-space) is considered. The strained state is characterized by two independent vectors: displacement and rotation. Solutions to the equations of motion are sought in the form of wave packets specified by an arbitrary Fourier spectrum. It is shown that, if the solution is sought in the form of a three-component displacement vector and a three-component rotation vector dependent on time, depth, and longitudinal coordinate, the initial system splits into two systems, one of which describes the Rayleigh wave and the other corresponds to a transverse wave decaying with depth. For both waves, analytical solutions in terms of displacements are obtained. It should be particularly noted that, unlike the Rayleigh wave, the solution for the transverse surface wave has no analogues in the classical elasticity theory. The transverse wave solution is numerically compared with the Rayleigh wave solution.