Solving Algebraic Equations

The term in parentheses in Eq. (8-17) is zero at steady state and thus it can be dropped. Next the Laplace transform is taken, and the resulting algebraic equation solved. Denoting X s) as the Laplace transform of and X,(.s) as the transform of 4, the final transfer Function can be written as ... 

Equation (7-54) allows calculation of the residence time required to achieve a given conversion or effluent composition. In the case of a network of reactions, knowing the reaction rates as a function of volumetric concentrations allows solution of the set of often nonlinear algebraic material balance equations using an implicit solver such as the multi variable Newton-Raphson method to determine the CSTR effluent concentration as a function of the residence time. As for batch reactors, for a single reaction all compositions can be expressed in terms of a component conversion or volumetric concentration, and Eq. (7-54) then becomes a single nonlinear algebraic equation solved by the Newton-Raphson method (for more details on this method see the relevant section this handbook). 

The Solve and NSolve commands are for algebraic equation solving. The Solve provides a symbolic result and NSolve numerically evaluates for the variable that is sought. These are used either for single or sets of equations. They are best illustrated by example. We can begin with Solve. 

APPENDIX A Review of Methods for Nonlinear Algebraic Equations Solving for a from Eq. A. 19, we obtain... 

The preceding examples show that the differentiation is equivalent to the multiplication by the parameter s and the integration is equivalent to the division by s in the Laplace domain. This allows for an easy transformation of differential or integral equations into algebraic equations, solving them in the Laplace domain and then carrying an inverse transformation into the time domain. This is schematically shown below ... 

If the given and wanted quantities are related by one or more Pea expressions, solve the problem by If the given and wanted quantities are related by an algebraic equation, solve the problem by... 

Equation (6.57) is a nonlinear algebraic equation. Solving it at each time step A t would require much CPU time, and so it is converted to an approximate linear algebraic equation. 

By contrast, a numerical computer program for solving such integration problems would depend on approximating the mathematical expression by a series of algebraic equations over expHcit integration limits. 

To solve the general problem using the backward Euler method, replace the nonlinear differential equation with the nonhuear algebraic equation for one step. 

This represents a set of nonlinear algebraic equations that can he solved with the Newton-Raphson method. However, in this case, a viable iterative strategy is to evaluate the transport coefficients at the last value and then solve... 

Elliptic Equations Elhptic equations can be solved with both finite difference and finite element methods. One-dimensional elhptic problems are two-point boundary value problems. Two- and three-dimensional elliptic problems are often solved with iterative methods when the finite difference method is used and direct methods when the finite element method is used. So there are two aspects to consider howthe equations are discretized to form sets of algebraic equations and howthe algebraic equations are then solved. 

Discretization of the governing equations. In this step, the exact partial differential equations to be solved are replaced by approximate algebraic equations written in terms of the nodal values of the dependent variables. Among the numerous discretization methods, finite difference, finite volume, and finite element methods are the most common. Tlxe finite difference method estimates spatial derivatives in terms of the nodal values and spacing between nodes. The governing equations are then written in terms of... 

Solution of the algebraic equations. For creeping flows, the algebraic equations are hnear and a linear matrix equation is to be solved. Both direct and iterative solvers have been used. For most flows, the nonlinear inertial terms in the momentum equation are important and the algebraic discretized equations are therefore nonlinear. Solution yields the nodal values of the unknowns. 

The UCKR.ON test problem assumes the simplest uniform surface implicitly, because adsorbed hydrogen coverage is directly proportional to the partial pressure of gaseous hydrogen and adversely affected by the partial pressure of the final products. Such a simple mechanism still amounts to a complex and unaccustomed rate expression of the type solved by second order algebraic equations. 

In the finite-difference appntach, the partial differential equation for the conduction of heat in solids is replaced by a set of algebraic equations of temperature differences between discrete points in the slab. Actually, the wall is divided into a number of individual layers, and for each, the energy conserva-tk>n equation is applied. This leads to a set of linear equations, which are explicitly or implicitly solved. This approach allows the calculation of the time evolution of temperatures in the wall, surface temperatures, and heat fluxes. The temporal and spatial resolution can be selected individually, although the computation time increa.ses linearly for high resolutions. The method easily can be expanded to the two- and three-dimensional cases by dividing the wall into individual elements rather than layers. 

More advanced models, for example the algebraic stress model (ASM) and the Reynolds stress model (RSM), are not based on the eddy-viscosity concept and can thus account for anisotropic turbulence thereby giving still better predictions of flows. In addition to the transport equations, however, the algebraic equations for the Reynolds stress tensor also have to be solved. These models are therefore computationally far more complex than simple closure models (Kuipers and van Swaaij, 1997). 

We continue by substituting the eigenvalues, in turn, into the algebraic equations (3-110). Because these equations are not independent, it is not possible to solve uniquely for the individual 7a, 7z values only ratios can be obtained, as follows ... 

The Laplace transformation is based upon the Laplace integral which transforms a differential equation expressed in terms of time to an equation expressed in terms of a complex variable a + jco. The new equation may be manipulated algebraically to solve for the desired quantity as an explicit function of the complex variable. 

Applications of Newton s Second Law. Problems involving no unbalanced couples can often be solved with the second law and the principles of kinematics. As in statics, it is appropriate to start with a free-body diagram showing all forces, decompose the forces into their components along a convenient set of orthogonal coordinate axes, and then solve a set of algebraic equations in each coordinate direction. If the accelerations are known, the solution will be for an unknown force or forces, and if the forces are known the solution will be for an unknown acceleration or accelerations. 

Substitute the equilibrium terms into the expression for K. This gives an algebraic equation that must be solved for x. 

---