Boundary value problems integrals

In integrated photoelasticity it is impossible to achieve a complete reconstruction of stresses in samples by only illuminating a system of parallel planes and using equilibrium equations of the elasticity theory. Theory of the fictitious temperature field allows one to formulate a boundary-value problem which permits to determine all components of the stress tensor field in some cases. If the stress gradient in the axial direction is smooth enough, then perturbation method can be used for the solution of the inverse problem. As an example, distribution of stresses in a bow tie type fiber preforms is shown in Fig. 2 [2]. 

In practice it would not be reasonable to solve the balances at the limit of Knudsen diffusion control by considering the n simultaneous boundary value problems (11.7). All the partial pressures can be expressed in terms of by integrating equations (11,25), with the result... 

Diffusion problems in one dimension lead to boundaiy value problems. The boundaiy conditions are applied at two different spatial locations at one side the concentration may be fixed and at the other side the flux may be fixed. Because the conditions are specified at two different locations, the problems are not initial value in character. It is not possible to begin at one position and integrate directly because at least one of the conditions is specified somewhere else and there are not enough conditions to begin the calculation. Thus, methods have been developed especially for boundary value problems. 

References Courant, R., and D. Hilbert, Methods of Mathematical Physics, vol. I, Interscience, New York (1953) Linz, P., Analytical and Numerical Methods for Volterra Equations, SIAM Publications, Philadelphia (1985) Porter, D., and D. S. G. Stirling, Integral Equations A Practical Treatment from Spectral Theory to Applications, Cambridge University Press (1990) Statgold, I., Greens Functions and Boundary Value Problems, 2d ed., Interscience, New York (1997). 

The limits of integration are fixed, and these problems are analogous to boundary value problems. 

References Brown, J. W., and R. V. Churchill, Fourier Series and Boundary Value Problems, 6th ed., McGraw-Hill, New York (2000) Churchill, R. V, Operational Mathematics, 3d ed., McGraw-Hill, New York (1972) Davies, B., Integral Transforms and Their Applications, 3d ed., Springer (2002) Duffy, D. G., Transform Methods for Solving Partial Differential Equations, Chapman Hall/CRC, New York (2004) Varma, A., and M. Morbidelli, Mathematical Methods in Chemical Engineering, Oxford, New York (1997). 

A differential equation for a function that depends on only one variable, often time, is called an ordinary differential equation. The general solution to the differential equation includes many possibilities the boundary or initial conditions are needed to specify which of those are desired. If all conditions are at one point, then the problem is an initial value problem and can be integrated from that point on. If some of the conditions are available at one point and others at another point, then the ordinary differential equations become two-point boundary value problems, which are treated in the next section. Initial value problems as ordinary differential equations arise in control of lumped parameter models, transient models of stirred tank reactors, and in all models where there are no spatial gradients in the unknowns. 

Packages to solve boundary value problems are available on the Internet. On the NIST web page http //gams.nist.gov/, choose problem decision tree and then differential and integral equations and then ordinary differential equations and multipoint boundary value problems. On the Netlibweb site http //www.netlib.org/, search on boundary value problem. Any spreadsheet that has an iteration capability can be used with the finite difference method. Some packages for partial differential equations also have a capability for solving one-dimensional boundary value problems [e.g. Comsol Multiphysics (formerly FEMLAB)]. 

An alternative (and much simpler) way of solving 20.4-6 is to covert it from a boundary- value problem to an initial-value problem, which may then be numerically integrated. 

Referring to Fig. 4.15, it is seen that the concentration and the concentration gradient are unknown at Z = 0. The above boundary condition relation indicates that if one is known, the other can be calculated. The condition of zero gradient at the outlet (Z=L) does not help to start the integration at Z=0, because, as Fig. 4.17 shows, two initial conditions are necessary. The procedure to solve this split-boundary value problem is therefore as follows ... 

A detailed treatment of the theoretical approach used in treating LSV and CV boundary value problems can be found in the monograph by MacDonald [23], More specific information on the numerical solution of integral equations common to electrochemical methods is available in the chapter by Nicholson [30]. The most commonly used method for the calculation of the theoretical electrochemical response, at the present time, is digital simulation which has been well reviewed by Feldberg [31, 32], Prater [33], Maloy [34], and Britz [35]. 

The axial dispersion model has led to the two-point boundary value problem (5.37) for uj from uj = ujstart = 0 to u> = uJend = 1- DEs are standardly solved by numerical integration over subintervals of the desired interval [uj start, we d. For more on the process of solving BVPs, see Section 1.2.4 or click on the Help line under the View icon on the MATLAB desktop, followed by a click on the Search tab in the Help window and searching for BVP . 

For initial value problems such as the ones we have encountered in Chapter 4, there is only one possible direction of integration, namely from the initial value ujstart onwards. But for two-point BVPs we have function information at both ends ujstart < atend and we could integrate forwards from uj start to ujend, or backwards from uiend to uj start-Regardless of the direction of integration, each boundary value problem on [uJstart, uJend]... 

The first one is based on a classical variation method. This approach is also known as an indirect method as it focuses on obtaining the solution of the necessary conditions rather than solving the optimization directly. Solution of these conditions often results in a two-point boundary value problem (TPBVP), which is accepted that it is difficult to solve [15], Although several numerical techniques have been developed to address the solution of TPBVP, e.g. control vector iteration (CVI) and single/multiple shooting method, these methods are generally based on an iterative integration of the state and adjoint equations and are usually inefficient [16], Another difficulty relies on the fact that it requires an analytical differentiation to derive the necessary conditions. 

For the boundary value problems in this paper, the Lx and L2 are linear Laplacian operators, the R and R2 disappear, the inhomogeneous terms g, (x) and g2(x) also equal to zero. The function f, f2 arising from integrating should be fiX = a0 + axx and... 

An expanded formulation of the steady-state permeation model has been presented. Two numerical problems - stiffness and an ill-conditioned boundary value problem - are encountered in solving the system equations. These problems can be circumvented by matching forward and reverse integrations at a point near the inlet (n = 0) but outside the combustion zone. The model predicts a... 

The approximating scheme converts the system of partial differential equations to a set of ordinary differential equations in the axial spatial coordinate. The detailed equations are contained in Yu et al. (2). The advantage of reduction in this manner is that the transient location of the combustion zone does not have to be known a priori, but can be found in the course of the integrations. Gear integration, which is designed for stiff systems, is used to solve the two point boundary value problem in the axial direction. 

Boundary Value Problem. If the conditions at the reactor outlet, such as the pressure and feed conversion are specified, it becomes a boundary value problem. This situation often arises when a reactor is to be designed for a specified feed conversion and for an outlet pressure dictated by downstream equipment. The solution of this problem involves a trial-and-error procedure and several passes of integration. 

Hanna, R. Fourier Series and Integrals of Boundary Value Problems, Wiley, New York (1982). 

---