Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation

An implicit sub-grid scale model for large eddy simulation is presented by utilising the concept of a relaxation system for one dimensional Burgers' equation in a novel way. The Burgers' equation is solved for three different unsteady flow situations by varying the ratio of relaxation parameter (epsilon) to time step. The coarse mesh results obtained with a relaxation scheme are compared with the filtered DNS solution of the same problem on a fine mesh using a fourth-order CWENO discretisation in space and third-order TVD Runge-Kutta discretisation in time. The numerical solutions obtained through the relaxation system have the same order of accuracy in space and time and they closely match with the filtered DNS solutions.

Exact analysis of Burgers equation on semiline with flux condition at the origin

The initial boundary value problem for the Burgers equation in the domain x greater-or-equal, slanted 0, t > 0 with flux boundary condition at x = 0 has been solved exactly. The behaviour of the solution as t tends to infinity is studied and the “asymptotic profile at infinity” is obtained. In addition, the uniqueness of the solution of the initial boundary value problem is proved and its inviscid limit as var epsilon → 0 is obtained.

Spectral approach for the soliton and periodic solutions of the nonlinear wave equation

A spectral method that obtains the soliton and periodic solutions to the nonlinear wave equation is presented. The results show that the nonlinear group velocity is a function of the frequency shift as well as of the soliton power. When the frequency shift is a function of time, a solution in terms of the Jacobian elliptic function is obtained. This solution is periodic in nature, and, to generate such an optical pulse train, one must simultaneously amplitude- and frequency-modulate the optical carrier. Finally, we extend the method to include the effect of self-steepening.

A formula for the solution of general Abel integral equation

A new formula for the solution of the general Abel Integral equation is derived, and an important special case is checked with the known result.

Exact N-Wave Solutions for the Non-Planar Burgers Equation

An exact representation of N-wave solutions for the non-planar Burgers equation u(t) + uu(x) + 1/2ju/t = 1/2deltau(xx), j = m/n, m < 2n, where m and n are positive integers with no common factors, is given. This solution is asymptotic to the inviscid solution for Absolute value of x < square-root (2Q0 t), where Q0 is a function of the initial lobe area, as lobe Reynolds number tends to infinity, and is also asymptotic to the old age linear solution, as t tends to infinity; the formulae for the lobe Reynolds numbers are shown to have the correct behaviour in these limits. The general results apply to all j = m/n, m < 2n, and are rather involved; explicit results are written out for j = 0, 1, 1/2, 1/3 and 1/4. The case of spherical symmetry j = 2 is found to be 'singular' and the general approach set forth here does not work; an alternative approach for this case gives the large time behaviour in two different time regimes. The results of this study are compared with those of Crighton & Scott (1979).

Generalised one-dimensional burgers'equation in relativistic fluiddynamics

This work deals with the effects of weak nonlinearity and weak dissipation on a linear wave in relativistic gasdynamics. Using perturbation and asymptotic expansions, a relativistic analogue of generalised one-dimensional Burgers' equation of classical gasdynamics is derived to describe far-field description of the wave. Steady state solution is presented for strict one-dimensional case.

On the solution of the equation ut+unux+H(x, t, u)=0

We consider the equation u(t) + u(n)u(x) + H(x, t, u) = 0 and derive a transformation relating it to u(t) + u(n)u(x) = 0. Special cases of the equation appearing in applications are discussed. Initial value problems and asymptotic behaviour of the solution are studied.

Resonant absorption of alfvén waves near a neutral point

It is shown that, although the mathematical analysis of the Alfven-wave equation does not show any variation at non-zero or zero singular points, the role of surface waves in the physical mechanism of resonant absorption of Alfven waves is very different at these points. This difference becomes even greater when resistivity is taken into account. At the neutral point the zero-frequency surface waves that are symmetric surface modes of the structured neutral layer couple to the tearing mode instability of the layer. The importance of this study for the energy balance in tearing modes and the association of surface waves with driven magnetic reconnection is also pointed out.

Applicability of Butler's equation in interpreting the thermodynamic behavior of surfaces and adsorption in Fe-S-O melts

Expressions for various second-order derivatives of surface tension with respect to composition at infinite dilution in terms of the interaction parameters of the surface and those of the bulk phases of dilute ternary melts have been presented. A method of deducing the parameters, which consists of repeated differentiation of Butler's equations with subsequent application of the appropriate boundary conditions, has been developed. The present investigation calculates the surface tension and adsorption functions of the Fe-S-O melts at 1873 and 1923 K using the modified form of Butler's equations and the derived values for the surface interaction parameters of the system. The calculated values are found to be in good agreement with those of the experimental data of the system. The present analysis indicates that the energetics of the surface phase are considerably different from those of the bulk phase. The present research investigates a critical compositional range beyond which the surface tension increases with temperature. The observed increase in adsorption of sulfur with consequent desorption of oxygen as a function of temperature above the critical compositional range has been ascribed to the increase of activity ratios of oxygen to sulfur in the surface relative to those in the bulk phase of the system.

Conjectures about phase turbulence in the complex Ginzburg-Landau equation

In the complex Ginzburg-Landau equation, we consider possible ''phase turbulent'' regimes, where asymptotic correlations are controlled by phase fluctuations rather than by topological defects. Conjecturing that the decay of such correlations is governed by the Kardar-Parisi-Zhang (KPZ) model of growing interfaces, we derive the following results: (1) A scaling ansatz implies that equal-time spatial correlations in 1d, 2d, and 3d decay like e(-Ax2 zeta), where A is a nonuniversal constant, and zeta=1/2 in 1d. (2) Temporal correlations decay as exp(-t(2 beta)h(t/L(z))), with the scaling law <(beta)over bar> = <(zeta)over bar>/z, where z = 3/2, 1.58..., and 1.66..., for d = 1,2, and 3 respectively. The scaling function h(y) approaches a constant as y --> 0, and behaves like y(2(beta-<(beta)over bar>)), for large y. If in 3d the associated KPZ model turns out to be in its weak-coupling (''smooth'') phase, then, instead of the above behavior, the CGLE exhibits rotating long-range order whose connected correlations decay like 1/x in space or 1/t(1/2) in time. (3) For system sizes, L, and times t respectively less than a crossover length, L(c), and time, t(c), correlations are governed by the free-field or Edwards-Wilkinson (EW) equation, rather than the KPZ model. In 1d, we find that L(c) is large: L(c) similar to 35,000; for L < L(c) we show numerical evidence for stretched exponential decay of temporal correlations with an exponent consistent with the EW value beta(EW)= 1/4.

Determination of the Gibbs energy of diamond using a solid state cell

he chemical potential of carbon in diamond, relative to its value in graphite, has been directly determined using a solid state electrochemical cell incorporating single crystal CaF2 as the solid electrolyte. The cell can be represented as Pt, C(graphite) + CaC2 + CaF2double vertical barCaF2double vertical barCaF2 + CaC2 + C(diamond), Pt The reversible emf of this cell is directly related by the Nernst equation to the Gibbs free energy change for the conversion of diamond to graphite. The difference in the chemical potential of carbon in the two crystal structures varies linearly with temperature in the range 940 to 1260 K ?C(diamond) ? ?C(graphite) = 1100 + 4.64T (±50) J mol?1 On the average, the values given by the equation are 320 J mol?1 less positive than the currently accepted ones based on calorimetric studies. The difference is primarily in the enthalpy term.

Effect of laser processing parameters on the structure of ductile iron

Laser processing of structure sensitive hypereutectic ductile iron, a cast alloy employed for dynamically loaded automative components, was experimentally investigated over a wide range of process parameters: from power (0.5-2.5 kW) and scan rate (7.5-25 mm s(-1)) leading to solid state transformation, all the way through to melting followed by rapid quenching. Superfine dendritic (at 10(5) degrees C s(-1)) or feathery (at 10(4) degrees C s(-1)) ledeburite of 0.2-0.25 mu m lamellar space, gamma-austenite and carbide in the laser melted and martensite in the transformed zone or heat-affected zone were observed, depending on the process parameters. Depth of geometric profiles of laser transformed or melt zone structures, parameters such as dendrile arm spacing, volume fraction of carbide and surface hardness bear a direct relationship with the energy intensity P/UDb2, (10-100 J mm(-3)). There is a minimum energy intensity threshold for solid state transformation hardening (0.2 J mm(-3)) and similarly for the initiation of superficial melting (9 J mm(-3)) and full melting (15 J mm(-3)) in the case of ductile iron. Simulation, modeling and thermal analysis of laser processing as a three-dimensional quasi-steady moving heat source problem by a finite difference method, considering temperature dependent energy absorptivity of the material to laser radiation, thermal and physical properties (kappa, rho, c(p)) and freezing under non-equilibrium conditions employing Scheil's equation to compute the proportion of the solid enabled determination of the thermal history of the laser treated zone. This includes assessment of the peak temperature attained at the surface, temperature gradients, the freezing time and rates as well as the geometric profile of the melted, transformed or heat-affected zone. Computed geometric profiles or depth are in close agreement with the experimental data, validating the numerical scheme.

Solution of a hypersingular integral equation of the second kind

A straightforward analysis involving the complex function-theoretic method is employed to determine the closed-form solution of a special hypersingular integral equation of the second kind, and its known solution is recovered.

An inverse problem for the Helmholtz equation involving two semi-infinite fluids

An analytical method is developed for solving an inverse problem for Helmholtz's equation associated with two semi-infinite incompressible fluids of different variable refractive indices, separated by a plane interface. The unknowns of the inverse problem are: (i) the refractive indices of the two fluids, (ii) the ratio of the densities of the two fluids, and (iii) the strength of an acoustic source assumed to be situated at the interface of the two fluids. These are determined from the pressure on the interface produced by the acoustic source. The effect of the surface tension force at the interface is taken into account in this paper. The application of the proposed analytical method to solve the inverse problem is also illustrated with several examples. In particular, exact solutions of two direct problems are first derived using standard classical methods which are then used in our proposed inverse method to recover the unknowns of the corresponding inverse problems. The results are found to be in excellent agreement.