Boundary value problems numerical solutions

Once the random field involved in the stochastic boundary value problem has been discretized, a solution method has to be adopted in order to solve the boundary value problem numerically. The choice of the solution method depends on the required statistical information of the solution. If only the first two statistical moments of the solution are of interest second moment analysis), the perturbafion method can be applied. However, if a. full probabilistic analysis is necessary, Galerkin schemes can be utilized or one has to resort to Monte Carlo simulations eventually in combination with a von Neumann series expansion. 

The weighted residual method provides a flexible mathematical framework for the construction of a variety of numerical solution schemes for the differential equations arising in engineering problems. In particular, as is shown in the followmg section, its application in conjunction with the finite element discretizations yields powerful solution algorithms for field problems. To outline this technique we consider a steady-state boundary value problem represented by the following mathematical model... 

NUMERICAL SOLUTIONS TO TWO-POINT BOUNDARY VALUE PROBLEMS... 

Finite element methods are one of several approximate numerical techniques available for the solution of engineering boundary value problems. Analysis of materials processing operations lead to equations of this type, and finite element methods have a number of advantages in modeling such processes. This document is intended as an overview of this technique, to include examples relevant to polymer processing technology. 

Difference equations with a symmetric matrix are typical in numerical solution of boundary-value problems associated with self-adjoint differential equations of second order. In what follows we will show that the condition Bi = is necessary and sufficient for the operator [yj] be self-adjoint. As can readily be observed, any difference equation of the form... 

Kubicek, M. Hlavacek, V. In Numerical Solution of Nonlinear Boundary Value Problems with Applications Prentice Hall Englewoods-Cliffs, N.Y., 1983. 

At the same time it is worth to notice that in modern numerical methods of a solution of boundary value problems, based on replacement of differential equations by finite difference, these steps are performed simultaneously. In accordance with the theorem of uniqueness, the field inside the volume V is defined by a distribution of masses inside this volume and boundary conditions, and correspondingly it is natural to derive an equation establishing this link. With this purpose in mind we will again proceed from Gauss s theorem,... 

Other methods can be used in space, such as the finite element method, the orthogonal collocation method, or the method of orthogonal collocation on finite elements. One simply combines the methods for ordinary differential equations (see Ordinary Differential Equations—Boundary Value Problems ) with the methods for initial-value problems (see Numerical Solution of Ordinary Differential Equations as Initial Value Problems ). Fast Fourier transforms can also be used on regular grids (see Fast Fourier Transform ). 

Absorption columns can be modeled in a plate-to-plate fashion (even if it is a packed bed) or as a packed bed. The former model is a set of nonlinear algebraic equations, and the latter model is an ordinary differential equation. Since streams enter at both ends, the differential equation is a two-point boundary value problem, and numerical methods are used (see Numerical Solution of Ordinary Differential Equations as Initial-Value Problems ). 

K. S. Yee, Numerical solution of initial boundary value problems involving Maxwell s equations in isotropic media, IEEE Transactions on Antennas and Propagation 14, 302-307 (1966). 

This is a linear ordinary-differential-equation boundary-value problem that can be solved analytically (see Bird, Stewart, and Lightfoot, Transport Phenomena, Wiley, 1960). Here, however, proceed directly to numerical finite-difference solution, which can be implemented easily in a spreadsheet. Assuming a cone angle of a = 2° and a rotation rate of 2 = 30 rpm, determine f(0) — v /r. 

Numerical Solution Equations 6.40 and 6.41 represent a nonlinear, coupled, boundary-value system. The system is coupled since u and V appear in both equations. The system is nonlinear since there are products of u and V. Numerical solutions can be accomplished with a straightforward finite-difference procedure. Note that Eq. 6.41 is a second-order boundary-value problem with values of V known at each boundary. Equation 6.40 is a first-order initial-value problem, with the initial value u known at z = 0. 

Solution by Shooting Solution of the boundary value problem described by Eq. 6.59 is usually accomplished numerically by a shooting method. To implement a shooting method, the third-order equations is transformed to a system of three first-order equations as... 

Figure E.l represents a highly simplified view of an ideal structure for an application program. The boxes with the rounded borders represent those functions that are problem specific, while the square-comer boxes represent those functions that can be relegated to problem-independent software. This structure is well-suited to problems that are mathematically systems of nonlinear algebraic equations, ordinary differential equation initiator boundary-value problems, or parabolic partial differential equations. In these cases the problem-independent mathematical software is usually written in the form of a subroutine that in turn calls a user-supplied subroutine to define the system of equations. Of course, the analyst must write the subroutine that describes the particular system of equations. Moreover, for most numerical-solution algorithms, the system of equations must be written in a discrete form (e.g., a finite-volume representation). However, the equation-defining sub-...

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]. 

Keller, H. B., 1976, Numerical Solution of Two-Point Boundary Value Problems. Regional Conf. Series in Appl. Math., SIAM. 

Figure 18. The most probable escape path (bottom solid curve) from S5 to the SI, found in the numerical simulations. The stable limit cycle is shown by rombs see Fig. 16 for other symbols. Parameters were h = 0.13, cty = 0.95,o>o 0.597,D = 0.0005. Top optimal force (solid line) corresponding to the optimal path after filtration [169]. The optimal path and optimal force from numerical solution of the boundary-value problem are shown by dots.

This result opens up the possibility of the numerical solution of the corresponding boundary value problem for energy-optimal control formulated above. 

Due to these inner iterations via IVP solvers and due to the need to solve an associated nonlinear systems of equations to match the local solutions globally, boundary value problems are generally much harder to solve and take considerably more time than initial value problems. Typically there are between 30 and 120 I VPs to solve numerous times in each successful run of a numerical BVP solver. 

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. 

Taking into account the aforementioned effects of ice formation in porous materials, a macroscopic quintuple model within the framework of the Theory of Porous Media (TPM) for the numerical simulation of initial and boundary value problems of freezing and thawing processes in saturated porous materials will be investigated. The porous solid is made up of a granular or structured porous matrix (a = S) and ice (a = I), where it will be assumed that both phases have the same motion. Due to the different freezing points of water in the macro and micro pores, the liquid will be distinguished into bulk water ( a = L) in the macro pores and gel water (a = P, pore solution) in the micro pores. With exception of the gas phase (a = G), all constituents will be considered as incompressible. 

The boundary value problem (Eqs. (10), (11)) is usually solved numerically. However, it is also possible to use another approach employing a linearization of this second-order, non-linear problem and a subsequent analytical treatment The analytical solution of the linearized boundary value problem in the film region is obtained in [15] ... 

Equations 1 through 6 can be solved numerically for temperature, pressure and concentration profiles along the length of the reactor, provided that the appropriate initial or boundary conditions are given. Solution methods are now discussed for both the initial and boundary value problems. However, emphasis will be placed on the latter since most industrial applications fall into this category. 

---