Showing papers on "Boundary value problem published in 1998"

Lattice boltzmann method for fluid flows

Abstract: We present an overview of the lattice Boltzmann method (LBM), a parallel and efficient algorithm for simulating single-phase and multiphase fluid flows and for incorporating additional physical complexities. The LBM is especially useful for modeling complicated boundary conditions and multiphase interfaces. Recent extensions of this method are described, including simulations of fluid turbulence, suspension flows, and reaction diffusion systems.

Inverse problems for partial differential equations

Abstract: Inverse Problems- Ill-Posed Problems and Regularization- Uniqueness and Stability in the Cauchy Problem- Elliptic Equations: Single Boundary Measurements- Elliptic Equations: Many Boundary Measurements- Scattering Problems- Integral Geometry and Tomography- Hyperbolic Problems- Inverse parabolic problems- Some Numerical Methods

PREMIX :A F ORTRAN Program for Modeling Steady Laminar One-Dimensional Premixed Flames

Abstract: This report documents a Fortran computer program that computes species and temperature profiles in steady-state burner-stabilized and freely propagating premixed laminar flames. The program accounts for finite rate chemical kinetics and multicomponent molecular transport. After stating the appropriate governing equations and boundary conditions, we discuss the finite difference discretization and the Newton method for solving the boundary value problem. Global convergence of this algorithm is aided by invoking time integration procedures when the Newton method has convergence difficulties. The program runs in conjunction with preprocessors for the chemical reaction mechanism and the transport properties. Transport property formulations include the option of using multicomponent or mixtureaveraged formulas for molecular diffusion. Discussion of two example problems illustrates many of the program's capabilities.

Thermomechanical analysis of functionally graded cylinders and plates

Abstract: The dynamic thermoelastic response of functionally graded cylinders and plates is studied. Thermomechanical coupling is included in the formulation, and a finite element model of the formulation is developed. The heat conduction and the thermoelastic equations are solved for a functionally graded axisymmetric cylinder subjected to thermal loading. In addition, a thermoelastic boundary value problem using the first-order shear deformation plate theory (FSDT) that accounts for the transverse shear strains and the rotations, coupled with a three-dimensional heat conduction equation, is formulated for a functionally graded plate. Both problems are studied by varying the volume fraction of a ceramic and a metal using a power law distribution.

The spectral element method: An efficient tool to simulate the seismic response of 2D and 3D geological structures

Abstract: We present the spectral element method to simulate elastic-wave propagation in realistic geological structures involving complieated free-surface topography and material interfaces for two- and three-dimensional geometries. The spectral element method introduced here is a high-order variational method for the spatial approximation of elastic-wave equations. The mass matrix is diagonal by construction in this method, which drastically reduces the computational cost and allows an efficient parallel implementation. Absorbing boundary conditions are introduced in variational form to simulate unbounded physical domains. The time discretization is based on an energy-momentum conserving scheme that can be put into a classical explicit-implicit predictor/multi-corrector format. Long-term energy conservation and stability properties are illustrated as well as the efficiency of the absorbing conditions. The associated Courant condition behaves as Δ tC < O ( nel−1/nd N −2), with nel the number of elements, nd the spatial dimension, and N the polynomial order. In practice, a spatial sampling of approximately 5 points per wavelength is found to be very accurate when working with a polynomial degree of N = 8. The accuracy of the method is shown by comparing the spectral element solution to analytical solutions of the classical two-dimensional (2D) problems of Lamb and Garvin. The flexibility of the method is then illustrated by studying more realistic 2D models involving realistic geometries and complex free-boundary conditions. Very accurate modeling of Rayleigh-wave propagation, surface diffraction, and Rayleigh-to-body-wave mode conversion associated with the free-surface curvature are obtained at low computational cost. The method is shown to provide an efficient tool to study the diffraction of elastic waves by three-dimensional (3D) surface topographies and the associated local effects on strong ground motion. Complex amplification patterns, both in space and time, are shown to occur even for a gentle hill topography. Extension to a heterogeneous hill structure is considered. The efficient implementation on parallel distributed memory architectures will allow to perform real-time visualization and interactive physical investigations of 3D amplification phenomena for seismic risk assessment.

The method of fundamental solutions for elliptic boundary value problems

Abstract: The aim of this paper is to describe the development of the method of fundamental solutions (MFS) and related methods over the last three decades. Several applications of MFS-type methods are presented. Techniques by which such methods are extended to certain classes of non-trivial problems and adapted for the solution of inhomogeneous problems are also outlined.

Positive Solutions of Differential, Difference and Integral Equations

Abstract: Preface. Ordinary Differential Equations. 1. First Order Initial Value Problems. 2. Second Order Initial Value Problems. 3. Positone Boundary Value Problems. 4. Semi-positone Boundary Value Problems. 5. Semi-Infinite Interval Problems. 6. Mixed Boundary Value Problems. 7. Singular Boundary Value Problems. 8. General Singular and Nonsingular Boundary Value Problems. 9. Quasilinear Boundary Value Problems. 10. Delay Boundary Value Problems. 11. Coupled System of Boundary Value Problems. 12. Higher Order Sturm-Liouville Boundary Value Problems. 13. (n,p) Boundary Value Problems. 14. Focal Boundary Value Problems. 15. General Focal Boundary Value Problems. 16. Conjugate Boundary Value Problems. Difference Equations. 17. Discrete Second Order Boundary Value Problems. 18. Discrete Higher Order Sturm-Liouville Boundary Value Problems. 19. Discrete (n,p) Boundary Value Problems. 20. Discrete Focal Boundary Value Problems. 21. Discrete Conjugate Boundary Value Problems. Integral and Integrodifferential Equations. 22. Volterra Integral Equations. 23. Hammerstein Integral Equations. 24. First Order Integrodifferential Equations. References. Authors Index. Subject Index.

The natural element method in solid mechanics

Abstract: The application of the Natural Element Method (NEM) 1; 2 to boundary value problems in two-dimensional small displacement elastostatics is presented. The discrete model of the domain consists of a set of distinct nodes N, and a polygonal description of the boundary @. In the Natural Element Method, the trial and test functions are constructed using natural neighbour interpolants. These interpolants are based on the Voronoi tessellation of the set of nodes N. The interpolants are smooth (C 1 ) everywhere, except at the nodes where they are C 0 . In one-dimension, NEM is identical to linear nite elements. The NEM interpolant is strictly linear between adjacent nodes on the boundary of the convex hull, which facilitates imposition of essential boundary conditions. A methodology to model material discontinuities and non-convex bodies (cracks) using NEM is also described. A standard displacement-based Galerkin procedure is used to obtain the discrete system of linear equations. Application of NEM to various problems in solid mechanics, which include, the patch test, gradient problems, bimaterial interface, and a static crack problem are presented. Excellent agreement with exact (analytical) solutions is obtained, which exemplies the accuracy and robustness of NEM and suggests its potential application in the context of other classes of problems|crack growth, plates, and large deformations to name a few. ? 1998 John Wiley & Sons, Ltd.

Direct analysis of particulate suspensions with inertia using the discrete Boltzmann equation

Abstract: An efficient and robust computational method, based on the lattice-Boltzmann method, is presented for analysis of impermeable solid particle(s) suspended in fluid with inertia. In contrast to previous lattice-Boltzmann approaches, the present method can be used for any solid-to-fluid density ratio. The details of the numerical technique and implementation of the boundary conditions are presented. The accuracy and robustness of the method is demonstrated by simulating the flow over a circular cylinder in a two-dimensional channel, a circular cylinder in simple shear flow, sedimentation of a circular cylinder in a two-dimensional channel, and sedimentation of a sphere in a three-dimensional channel. With a solid-to-fluid density ratio close to one, new results from two-dimensional and three-dimensional computational analysis of dynamics of an ellipse and an ellipsoid in a simple shear flow, as well as two-dimensional and three-dimensional results for sedimenting ellipses and prolate spheroids, are presented.

Boundary Value Problems of Mathematical Physics

Abstract: The goal of this final chapter is to show how the boundary value problems of mathematical physics can be solved by the methods of the preceding chapters. This will be done by solving a variety of specific problems that illustrate the principal types of problems that were formulated in Chapter 7. Additional applications are developed in the Exercises. The primary solution method is Fourier’s method of separation of variables and the associated Sturm-Liouville theory of Chapter 8.

Quantum codes on a lattice with boundary

Abstract: A new type of local-check additive quantum code is presented. Qubits are associated with edges of a 2-dimensional lattice whereas the stabilizer operators correspond to the faces and the vertices. The boundary of the lattice consists of alternating pieces with two different types of boundary conditions. Logical operators are described in terms of relative homology groups.

A local boundary integral equation (LBIE) method in computational mechanics, and a meshless discretization approach

Abstract: The Galerkin finite element method (GFEM) owes its popularity to the local nature of nodal basis functions, i.e., the nodal basis function, when viewed globally, is non-zero only over a patch of elements connecting the node in question to its immediately neighboring nodes. The boundary element method (BEM), on the other hand, reduces the dimensionality of the problem by one, through involving the trial functions and their derivatives, only in the integrals over the global boundary of the domain; whereas, the GFEM involves the integration of the “energy” corresponding to the trial function over a patch of elements immediately surrounding the node. The GFEM leads to banded, sparse and symmetric matrices; the BEM based on the global boundary integral equation (GBIE) leads to full and unsymmetrical matrices. Because of the seemingly insurmountable difficulties associated with the automatic generation of element-meshes in GFEM, especially for 3-D problems, there has been a considerable interest in element free Galerkin methods (EFGM) in recent literature. However, the EFGMs still involve domain integrals over shadow elements and lead to difficulties in enforcing essential boundary conditions and in treating nonlinear problems. The object of the present paper is to present a new method that combines the advantageous features of all the three methods: GFEM, BEM and EFGM. It is a meshless method. It involves only boundary integration, however, over a local boundary centered at the node in question; it poses no difficulties in satisfying essential boundary conditions; it leads to banded and sparse system matrices; it uses the moving least squares (MLS) approximations. The method is based on a Local Boundary Integral Equation (LBIE) approach, which is quite general and easily applicable to nonlinear problems, and non-homogeneous domains. The concept of a “companion solution” is introduced so that the LBIE for the value of trial solution at the source point, inside the domain Ω of the given problem, involves only the trial function in the integral over the local boundary Ω s of a sub-domain Ω s centered at the node in question. This is in contrast to the traditional GBIE which involves the trial function as well as its gradient over the global boundary Γ of Ω. For source points that lie on Γ, the integrals over Ω s involve, on the other hand, both the trial function and its gradient. It is shown that the satisfaction of the essential as well as natural boundary conditions is quite simple and algorithmically very efficient in the present LBIE approach. In the example problems dealing with Laplace and Poisson's equations, high rates of convergence for the Sobolev norms ||·||0 and ||·||1 have been found. In essence, the present EF-LBIE (Element Free-Local Boundary Integral Equation) approach is found to be a simple, efficient, and attractive alternative to the EFG methods that have been extensively popularized in recent literature.

A numerical investigation on the effect of the inflow conditions on the self-similar region of a round jet

Abstract: In this paper we consider the direct numerical simulation (DNS) of a spatially developing free round jet at low Reynolds numbers. Simulation of a spatially evolving flow such as the jet requires boundary conditions, which allow entrainment into the turbulent flow across the lateral boundaries of the computational domain. The boundary conditions which satisfy this requirement are so-called traction free boundary conditions. After showing that these boundary conditions lead to a correct behavior of the velocity near the lateral boundary of the jet, we will consider the DNS of the jet flow at a Reynolds number of 2.4×103 and compare the results with experimental data obtained by Hussein et al. [J. Fluid Mech. 258, 31 (1994)] and by Panchapakesan and Lumley [J. Fluid Mech. 246, 197 (1993)]. The results of our numerical simulations agree very well with the experimental data. Next we use the DNS to investigate the influence of the shape of the velocity profile at the jet orifice on the self-similarity scaling f...

Noninvasive determination of the optical properties of two-layered turbid media

Abstract: Light propagation in two-layered turbid media having an infinitely thick second layer is investigated in the steady-state, frequency, and time domains. A solution of the diffusion approximation to the transport equation is derived by employing the extrapolated boundary condition. We compare the reflectance calculated from this solution with that computed with Monte Carlo simulations and show good agreement. To investigate if it is possible to determine the optical coefficients of the two layers and the thickness of the first layer, the solution of the diffusion equation is fitted to reflectance data obtained from both the diffusion equation and the Monte Carlo simulations. Although it is found that it is, in principle, possible to derive the optical coefficients of the two layers and the thickness of the first layer, we concentrate on the determination of the optical coefficients, knowing the thickness of the first layer. In the frequency domain, for example, it is shown that it is sufficient to make relative measurements of the phase and the steady-state reflectance at three distances from the illumination point to obtain useful estimates of the optical coefficients. Measurements of the absolute steady-state spatially resolved reflectance performed on two-layered solid phantoms confirm the theoretical results.

Analyses of slow high-concentration flows of granular materials

Abstract: A theory intended for slow, dense flows of cohesionless granular materials is developed for the case of planar deformations. By considering granular flows on very fine scales, one can conveniently split the individual particle velocities into fluctuating and mean transport components, and employ the notion of granular temperature that plays a central role in rapid granular flows. On somewhat larger scales, one can think of analogous fluctuations in strain rates. Both kinds of fluctuations are utilized in the present paper. Following the standard continuum approach, the conservation equations for mass, momentum and particle translational fluctuation energy are presented. The latter two equations involve constitutive coefficients, whose determination is one of the main concerns of the present paper. We begin with an associated flow rule for the case of a compressible, frictional, plastic continuum. The functional dependence of the flow rule is chosen so that the limiting behaviours of the resulting constitutive relations are consistent with the results of the kinetic theories developed for rapid flow regimes. Following Hibler (1977) and assuming that there are fluctuations in the strain rates that have, for example, a Gaussian distribution function, it is possible to obtain a relationship between the mean stress and the mean strain rate. It turns out, perhaps surprisingly, that this relationship has a viscous-like character. For low shear rates, the constitutive behaviour is similar to that of a liquid in the sense that the effective viscosity decreases with increasing granular temperature, whereas for rapid granular flows, the viscosity increases with increasing granular temperature as in a gas. The rate of energy dissipation can be determined in a manner similar to that used to derive the viscosity coefficients. After assuming that the magnitude of the strain-rate fluctuations can be related to the granular temperature, we obtain a closed system of equations that can be used to solve boundary value problems. The theory is used to consider the case of a simple shear flow. The resulting expressions for the stress components are similar to models previously proposed on a more ad hoc basis in which quasi-static stress contributions were directly patched to rate-dependent stresses. The problem of slow granular flow in rough-walled vertical chutes is then considered and the velocity, concentration and granular temperature profiles are determined. Thin boundary layers next to the vertical sidewalls arise with the concentration boundary layer being thicker than the velocity boundary layer. This kind of behaviour is observed in both laboratory experiments and in granular dynamics simulations of vertical chute flows.

Realization of Fluid Boundary Conditions via Discrete Boltzmann Dynamics

Abstract: We describe a novel way based on lattice-Boltzmann representation for realizing hydrodynamic boundary conditions at a solid surface. It is shown that using this approach the resulting physics properties are independent of the position and the orientation of the surface with respect to the lattice mesh. The fluxes of mass, energy as well as both normal and tangential momenta can be accurately controlled to correspond to various fluid dynamics situations.

A systematic development of 'solid-shell' element formulations for linear and non-linear analyses employing only displacement degrees of freedom

Abstract: In the present contribution we propose a so-called solid-shell concept which incorporates only displacement degrees of freedom. Thus, some major disadvantages of the usually used degenerated shell concept are overcome. These disadvantages are related to boundary conditions—the handling of soft and hard support, the need for special co-ordinate systems at boundaries, the connection with continuum elements—and, in geometrically non-linear analyses, to a complicated update of the rotation vector. First, the kinematics of the so-called solid-shell concept in analogy to the degenerated shell concept are introduced. Then several modifications of the solid-shell concept are proposed to obtain locking-free solid-shell elements, leading also to formulations which allow the use of general three-dimensional material laws and which are also able to represent the normal stresses and strains in thickness direction. Numerical analyses of geometrically linear and non-linear problems are finally performed using solely assumed natural shear strain elements with a linear approximation in in-plane direction. Although some considerations are needed to get comparable boundary conditions in the examples analysed, the solid-shell elements prove to work as good as the degenerated shell elements. The numerical examples show that neither thickness nor shear locking are present even for distorted element shapes. © 1998 John Wiley & Sons, Ltd.