Loeb solutions of the boltzmann equation

Existence problems for the Boltzmann equation constitute a main area of research within the kinetic theory of gases and transport theory. The present paper considers the spatially periodic case with L¹ initial data. The main result is that the Loeb subsolutions obtained in a preceding paper are shown to be true solutions. The proof relies on the observation that monotone entropy and finite energy imply Loeb integrability of non-standard approximate solutions, and uses estimates from the proof of the H-theorem. Two aspects of the continuity of the solutions are also considered.     

Artificial boundary conditions for Euler–Bernoulli beam equation

In a semi-discretized Euler-Bernoulli beam equa- tion, the non-nearest neighboring interaction and large span of temporal scales for wave propagations pose challenges to the effectiveness and stability for artificial boundary treat- ments. With the discrete equation regarded as an atomic lattice with a three-atom potential, two accurate artificial boundary conditions are first derived here. Reflection co- efficient and numerical tests illustrate the capability of the proposed methods. In particular, the time history treatment gives an exact boundary condition, yet with sensitivity to nu- merical implementations. The ALEX （almost EXact） bound- ary condition is numerically more effective.     

General boundary conditions for the wave equation around non-homogenous scatterers

To solve the wave equation inside a region that contains an inhomogenous dielectric material of arbitrary shape under the influence of an incoming wave, we establish a generalized boundary condition. The solutions inside a finite region resulting form a given incoming wave from the outside, are determined by a linear relation between the normal gradient and the function values on the boundary. This boundary condition is non-local and we show how it can be used in conjunction with the variational principle applied to an open system.     

Solution to the boltzmann equation for monoenergetic neutrons in a slab

Summary This paper presents a physical solution to the linear integro-differential Boltzmann equation, which governs the stationary distribution of the flux of monoenergetic neutrons in the interior and at the boundary of an infinite slab of finite thickness. Results for both the angular and the total fluxes are obtained by properly accounting for the boundary conditions of the problem and by a Fourier transform of the original Boltzmann equation. An integral linear Fredholm equation for the Fourier transform of the total flux is thus derived. The solution of such an equation is given two different representations according to the two decompositions devised for the given Fredholm kernel. The solution obtained in each case is finally inverted termwise to yield the sought total flux in the original space. The result for the angular flux follows, then, directly from the result for the total flux.

---

Sommario In questa nota viene costruita una soluzione fisica dell'equazione integro-differenziale lineare di Boltzmann, che, in regime stazionario, regola la distribuzione del flusso di neutroni monoenergetici all'interno di una lastra piana di spessore finito. Risultati per la distribuzione sia del flusso angolare che totale vengono ottenuti tramite una trasformazione di Fourier dell'equazione integro-differenziale di partenza, in cui le condizioni al contorno del problema vengono opportunamente incorporate. Si perviene in tal modo ad una equazione integrale lineare di Fredholm per la trasformata di Fourier del flusso totale. Della soluzione di questa equazione vengono date due diverse rappresentazioni, che corrispondono a due possibili decomposizioni del nucleo dell'equazione stessa. In ciascuno dei due casi la soluzione ottenuta viene quindi antitrasformata termine a termine per riprodurre il flusso totale nello spazio originale. Dal risultato per il flusso totale segue poi direttamente quello del flusso angolare.

---

---

Work performed under contract C.N.R.     

Process splitting of the boundary conditions for the advection–dispersion equation

Rational strategies are considered for the specification of the intermediate boundary condition at an inflow boundary where process splitting (fractional steps) is adopted in solving the advection–dispersion equation. Three lowest-order methods are initially considered and evaluation is based on comparisons with an analytical solution. For flow and dispersion parameter ranges typical of rivers and estuaries, the given boundary condition for the complete advection–dispersion equation at the end of the complete time step provides a satisfactory estimate of the intermediate boundary value. This was further confirmed by the development and evaluation of two higher-order methods. These required non-centred discrete approximations for spatial derivatives, which offset any special advantages from the higher truncation error order.     

On accurate boundary conditions for a shape sensitivity equation method

This paper studies the application of the continuous sensitivity equation method (CSEM) for the Navier–Stokes equations in the particular case of shape parameters. Boundary conditions for shape parameters involve flow derivatives at the boundary. Thus, accurate flow gradients are critical to the success of the CSEM. A new approach is presented to extract accurate flow derivatives at the boundary. High order Taylor series expansions are used on layered patches in conjunction with a constrained least‐squares procedure to evaluate accurate first and second derivatives of the flow variables at the boundary, required for Dirichlet and Neumann sensitivity boundary conditions. The flow and sensitivity fields are solved using an adaptive finite‐element method. The proposed methodology is first verified on a problem with a closed form solution obtained by the Method of Manufactured Solutions. The ability of the proposed method to provide accurate sensitivity fields for realistic problems is then demonstrated. The flow and sensitivity fields for a NACA 0012 airfoil are used for fast evaluation of the nearby flow over an airfoil of different thickness (NACA 0015). Copyright © 2005 John Wiley & Sons, Ltd.     

Global attraction for the one-dimensional heat equation with nonlinear time-dependent boundary conditions

Archive for Rational Mechanics and Analysis -     

Dirichlet and Neumann boundary conditions for the pressure poisson equation of incompressible flow

In a recent paper Gresho and Sani showed that Dirichlet and Neumann boundary conditions for the pressure Poisson equation give the same solution. The purpose of this paper is to confirm this (for one case at least) by numerically solving the pressure equation with Dirichlet and Neumann boundary conditions for the inviscid stagnation point flow problem. The Dirichlet boundary condition is obtained by integrating the tangential component of the momentum equation along the boundary. The Neumann boundary condition is obtained by applying the normal component of the momentum equation at the boundary. In this work solutions for the Neumann problem exist only if a compatibility condition is satisfied. A consistent finite difference procedure which satisfies this condition on non-staggered grids is used for the solution of the pressure equation with Neumann conditions. Two test cases are computed. In the first case the velocity field is given from the analytical solution and the pressure is recovered from the solution of the associated Poisson equation. The computed results are identical for both Dirichlet and Neumann boundary conditions. However, the Dirichlet problem converges faster than the Neumann case. In the second test case the velocity field is computed from the momentum equations, which are solved iteratively with the pressure Poisson equation. In this case the Neumann problem converges faster than the Dirichlet problem.     

Finite element implementation of boundary conditions for the pressure Poisson equation of incompressible flow

In this paper we address the problem of the implementation of boundary conditions for the derived pressure Poisson equation of incompressible flow. It is shown that the direct Galerkin finite element formulation of the pressure Poisson equation automatically satisfies the inhomogeneous Neumann boundary conditions, thus avoiding the difficulty in specifying boundary conditions for pressure. This ensures that only physically meaningful pressure boundary conditions consistent with the Navier-Stokes equations are imposed. Since second derivatives appear in this formulation, the conforming finite element method requires C¹ continuity. However, for many problems of practical interest (i.e. high Reynolds numbers) the second derivatives need not be included, thus allowing the use of more conventional C⁰ elements. Numerical results using this approach for a wall-driven contained flow within a square cavity verify the validity of the approach. Although the results were obtained for a two-dimensional problem using the p-version of the finite element method, the approach presented here is general and remains valid for the conventional h-version as well as three-dimensional problems.     

Effect of boundary conditions on turbulence characteristics in two-phase jet flows with phase transitions

It is shown that there are no self-similarity and similarity of transverse fields of correlation moments of fluctuating parameters of phases in two-phase jets, in contrast to one-phase jets. The influence of the initial values of a number of parameters of a two-phase jet (gas temperature, volume concentration of droplets in the initial cross section, and radius of the initial cross section of the jet) on turbulence characteristics is analyzed on the basis of numerical simulations. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 29–40, May–June, 2005.     

Supercritical equilibrium solutions of axially moving beams with hybrid boundary conditions

In this paper supercritical equilibria and critical speeds of axially moving beams constrained by sleeves with torsion springs are deduced. Transverse vibration of the beams is governed by a nonlinear integro-partial-differential equation. In the supercritical regime, the corresponding static equilibrium equation for the hybrid boundary conditions is analytically solved for the equilibria and the critical speeds. In the view of the non-trivial equilibrium, comparisons are made among the integro-partial-differential equation, a nonlinear partial-differential equation for transverse vibration, and coupled equations for planar motion under the hybrid boundary conditions.     

Study of stable and unstable jet flows on the basis of the boltzmann equation

Free supersonic underexpanded jets are studied using a direct method conservative splitting scheme for solving the Boltzmann equation. Numerical solutions for a jet flowing into a vacuum and into a fluid-filled space are presented for the following ranges of the parameters: Knudsen number 10⁻⁶<Kn<∞ and pressure ratio 10<n<∞. The solutions are compared with experimental data. Instabilities associated with free turbulence effects in the mixing layer are detected for low Kn numbers. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 153–157, March–April, 1998. The work was carried out with support from the Russian Foundation for Fundamental Research (project No. 96-01-00829).     

On initial,boundary conditions and viscosity coefficient control for Burgers' equation

In order to use the optimal control techniques in models of geophysical flow circulation, an application to a 1D advection–diffusion equation, the so-called Burgers' equation, is described. The aim of optimal control is to find the best parameters of the model which ensure the closest simulation to the observed values. In a more general case, the continuous problem and the corresponding discrete form are formulated. Three kinds of simulation are realized to validate the method. Optimal control processes by initial and boundary conditions require an implicit discretization scheme on the first time step and a decentered one for the non-linear advection term on boundaries. The robustness of the method is tested with a noised dataset and random values of the initial controls. The optimization process of the viscosity coefficient as a time- and space-dependent variable is more difficult. A numerical study of the model sensitivity is carried out. Finally, the numerical application of the simultaneous control by the initial conditions, the boundary conditions and the viscosity coefficient allows a possible influence between controls to be taken into account. These numerical experiments give methodological rules for applications to more complex situations. © 1998 John Wiley & Sons, Ltd.     

High-order boundary conditions for the problems of Laplace equation in infinite region and their application

The high-order boundary conditions for the problems of Laplace equation in infiniteregion have been developed.The improvement in accuracy for numerical solution isachieved by imposing the high-order boundary conditions on the exterior boundary of areduced finite region in which the numerical method is used.So both the computing effortsand the required storage in computer are reduced.The numerical examples show that thelst-order boundary condition approaches to the exact boundary condition and it is clearlysuperior to the traditional boundary condition and the2nd-order boundary condition.