Computational Geosciences 2 (1998) 1, 23-36
|
|
- Alban Matthews
- 5 years ago
- Views:
Transcription
1 A STUDY OF THE MODELLING ERROR IN TWO OPERATOR SPLITTING ALGORITHMS FOR POROUS MEDIA FLOW K. BRUSDAL, H. K. DAHLE, K. HVISTENDAHL KARLSEN, T. MANNSETH Comptational Geosciences 2 (998), Abstract. Operator splitting methods are often sed to solve convection-dision problems of convection dominated natre. However, it is well known that sch methods can prodce signicant (splitting) errors in regions containing self sharpening fronts. To amend this shortcoming, corrected operator splitting methods have been developed. These approaches se the wave strctre from the convection step to identify the splitting error. This error is then compensated for in the dision step. The main prpose of the present work is to illstrate the importance of the correction step in the context of an inverse problem. The inverse problem will consist of estimating the fractional ow fnction in a one-dimensional satration eqation.. Introdction We are here interested in the initial-bondary vale problem associated with nonlinear convection-dision eqations of f() ; () = (x t) 2 (a b) where = (x t) denotes the nknown, f = f() is the x fnction, d = d() is the dision fnction and " > is a scaling parameter. Convection-dision eqations arise in a variety of applications. We here consider two-phase, immiscible and incompressible ow of oil and water in poros media. In this context, is the water satration, f() is the fractional ow fnction and d() is the capillary dision fnction. De to the nonlinear natre of the dierential eqation (), there will in general be no analytical expression for (x t). Hence we mst rely on nmerically techniqes to obtain the soltion. When () is convection dominated, that is, when " is small compared with Key words and phrases. Operator splitting, two-phase ow, satration eqation, modelling error, inverse problem. Brsdal and Mannseth have been spported by the Research Concil of Norwaynder grant 765/42. Karlsen has been spported by VISTA, a research cooperation between the Norwegian Academy of Science and Letters and Den norske stats oljeselskap a.s. (Statoil).
2 2 BRUSDAL ET AL. other scales in (), conventional methods sally exhibit some combination of diclties, ranging from non-physical oscillations to severe nmerical dision at the trailing end of moving fronts. To overcome sch diclties we make se of operator splitting. In the present context, operator splitting means to split the convection-dision eqation () into ahyperbolic eqation and a parabolic eqation, each of which is solved separately in an alternate fashion. Frthermore, each step of sch an algorithm is flly discretized by accrate and ecient nmerical schemes developed especially for hyperbolic conservation laws and parabolic heat type eqations, respectively. The operator splitting approach has been taken by many athors, [2, 4, 5, 6,7,8, 9,, 2, 5, 6,7] (see the references cited therein for a more complete list of relevant papers). In the case where d(), the physical front widths shold be of order O("). However, it is easy to see that disive splitting errors will lead to front widths of size O( p "t) when the splitting time step, t, is large ( "). This type of splitting error comes from thefactthattheentropy condition forces the hyperbolic solver to throw away certain \information" that controls the strctre of the self-sharpening fronts. It is however possible to compensate for this loss of information, which manifests itself in the form of residal x terms. For example, assme that the soltion of () is a moving front. If the time step is large enogh, the convection step will generate a shock with left and right limits, say, l and r, respectively. We can then identify a residal x term associated with this discontinity f res () = f() ; f c (), where f c () denotes the correct envelope (dictated by the entropy condition) of f() on the interval bonded by l and r. There are several ways to take the residal x term f res () into accont. We can perform a separate correction step after the dision step. Correction is then calclated by solving the \residal" conservation law v t + f res (v) x = over a time interval ( ], where > is some parameter that has to be chosen (see [6] for details). Another approach is to inclde the residal term in the eqation modelling dision, that is, instead of solving the eqation w t ; "(d(w)w x ) x =,we solve w t +(f res (w) ; "d(w)w x ) x =. This eqation contains the necessary information needed to prodce the correct strctre of the front. We will here rely on the latter approach, which we shall refer to as corrected operator splitting (COS). Splitting methods which do not accont for residal x terms will be denoted by OS. The idea of sing a residal x term in the dision step was introdced in [] (and frther developed in [4, 6, 7, 8]) in the context of two-phase ow with an established front as initial data. In this setting the residal term can be determined a priori, which
3 OPERATOR SPLITTING AND CORRECTED OPERATOR SPLITTING 3 means that the convection step is to solve v t + f c (v) x = (where f c is linear in the front region). This fact was exploited in [4, 6, 7, 8, ] in the sense that the modied method of characteristics was sed for the convection step. For general initial data this approach may pt severe restrictions on the time step. However, in [5, 6] itwas observed that by sing front tracking [4] to solve the nonlinear conservation law v t + f(v) x = it is possible to dynamically constrct residal x terms. Conseqently, the COS approach makes sense in general. The calclation of these residal terms is a direct conseqence of the fact that front tracking is based on solving Riemann problems. Instead of front tracking, a second order Godnov method will here be sed in the convection step to illstrate that other techniqes can be tilised in the constrction of residal x terms, see also [5]. In what follows, dierences in the Godnov based OS and COS methods will be illstrated in the context of an inverse problem. The inverse problem will consist of recovering the x fnction sing COS and OS as forward models in the inversion procedre. Note that in general it is not possible to dene the residal x term a priori in the inversion procedre since the representation of the x fnction is changing throghot the procedre. This motivates the se of COS where the residal x terms are constrcted dynamically. The main prpose of the present work is to stdy the inence that the modelling errors in the OS and COS forwards models have on the recovery of the x fnction. The paper is organised as follows. In section 2 we give a description of the forward problem and the operator splitting techniqes sed to solve it. In section 3we discss the setp of the inverse problem. Frthermore, we present and discss nmerical experiments. Finally, in section 4 we make some conclding remarks. 2. The forward problem 2.. Nmerical algorithms. In this section we describe in some detail the operator splitting methods that are sed to solve the forward eqation (). This eqation mst be eqipped with appropriate initial and bondary conditions conditions t= = (x) x=a = a 2 R x=b = b 2 R where we assme that is a piecewise smooth fnction. Let n (x) (x t n ) denote the piecewise linear nmerical soltion of () at time level t n, where =t <t < t N = T is a time discretization of [ T]. Here the fnction n = n (x) is piecewise linear with
4 4 BRUSDAL ET AL. respect to a niform grid with grid cell size x >. We determine n+ from n via the following OS algorithm. Convection: Let v(x t n ), t n = t n+ ; t n >, denote the (entropy weak) soltion of the hyperbolic conservation f(v) v(x ) = n (x): The second order Godnov (slope limiter) method described in [3] is sed to calclate this soltion. The Godnov method ses a niform (local) grid with grid cell size x c and a time step < t c t n chosen according to the CFL condition, i.e., t c is chosen so that max jf j t c x c =. Dision: Let w(x t n )bethe soltion of the parabolic heat type @x = w(x ) = v(x t n ) with the Dirichlet bondary data imposed at x = a b. Nmerically this soltion is fond sing the Galerkin method on a niform grid with grid cell x d x and the sal (piecewise linear) \hat" basis fnctions. The time discretization is done by the backward Eler method with a single time step of length t n. Finally, we set n+ = w(x t n ): In the constant dision case, d(), convergence of the OS algorithm is proved in [7]. Convergence reslts in the case where d() is nonlinear and possibly strongly degenerate are obtained in []. Althogh this algorithm converges as varios discretization parameters tend to zero, it is not diclt to constrct an explicit example (see [6]) which shows that OS is too disive in regions containing self-sharpening fronts, at least when the time step is large ( "). This potential shortcoming has motivated the corrected operator splitting (COS) algorithm, which diers from OS in the way that the dision step is carried ot. Dision revisited (residal x terms): To obtain correct strctre of fronts, a residal x term has to be constrcted. To thisend, we introdce the envelope fnction 8 < the lower convex envelope of f between v l and v r if v l <v r, (4) f c (v v l v r )= : the pper concave envelope of f between v r and v l if v l >v r. Let v = v(x t n ) be the tre soltion of (2) at time t = t n, which is known to be piecewise smooth. Let fy i g denote the discontinity points of v. For a xed i, let ~v i and
5 OPERATOR SPLITTING AND CORRECTED OPERATOR SPLITTING Figre. The left plot shows an exact soltion (solid line) from the convection step, which contains for discontinities (left and right shock states are shown by circles). The thick solid lines indicate where the corresponding residal x terms are dened when sed in the dision step. The right plots show the residal x terms (solid line) associated with the for discontinities (the pper left plot corresponds to the rst discontinity, the pper right plot to the second and so on). The x fnction is shown as the thick dashed line and the envelope fnction as dotted. ~v i+ denote the left and right limits of v at x = y i, respectively. Choose positions x i and x i+ sch that y i is located somewhere in the interval (x i x i+ ). Then dene the residal x fnction associated with this (ith) discontinity by (see Figre ) 8 < f(v) ; fres(v i f c (v ~v i ~v i+ ) for v 2 [~v i ~v i+ ] t n )= : for v 62 [~v i ~v i+ ]: For convenience, introdce the globally dened residal x term f res (x v t n )=f i res(v t n ), x 2 [x i x i+ ) for some i. Let now w(x t n ) be the weak soltion of the following f res (x w t n ) ; = w(x ) = v(x t
6 6 BRUSDAL ET AL. with the Dirichlet bondary data imposed at x = a b. Nmerically this soltion is fond by a Petrov-Galerkin method with qadratic test fnctions, see e.g. [4, 6, 7,, 5] for details. To nd a residal x term, we need to locate the associated discontinities (shocks) based on the approximate soltion generated by the Godnov method. This is done by introdcing a sitable parameter > and identify a shock whenever the distance between two sbseqent satration vales in the Godnov soltion is greater than, see [5] for details. Remark. The loss of symmetry de to the residal x term does not aect the overrall comptational time for the COS method since the tridiagonal linear systems are solved by direct methods. This claim is conrmed in [5] where the nmerical experiments show that COS' rntime (with one or a few time steps) is less than that of OS (with many time steps). In higher dimensions this loss of symmetry may become an isse since iterative methods have to be sed to solve the reslting linear systems. However, since the eqations are (partially) symmetrized by the se of a Petrov-Galerkin method this isse seems to be of mch less consern than the restriction on the time step imposed for the OS method. Remark. A sitable vale of is problem dependent and the parameter is adjsted by considering the nmber of grid blocks sed in the nmerical comptations. Therefore, changing cold reslt in a less accrate constrction of the residal x fnction. However, since the nmber of grid blocks is the same in each nmerical experiment, the parameter is kept xed Nmerical experiments. For two phase ow in poros media [3] the advective x fnction in () is given by (5) f() = w () w ()+ o () where i = k ri i, is the mobility of phase i, i = w o. Here, k ri, denotes the relative permeability and i the viscosity. The relative permeabilities are represented by third order normalised B-spline expansions [9] in which the coecient vectors are given by ~c w = (: : :52 :) and ~c o = (: :25 :2 :). The knot vector for the B-spline basis is ~y = (: : : :5 : : :). Frthermore, we set d(), = :5 and " =:.
7 OPERATOR SPLITTING AND CORRECTED OPERATOR SPLITTING 7 The satration proles are simlated sing axfnction, f (), which corresponds to the B-spline expansion given above (see Figre 2a)). The initial satration (see Figre 2b)) is 8 >< : ; x : x :3 (x) = :7 ; 7 (x ; :3) :3 x :3 >: : :3 <x:: The Dirichlet bondary data are = at x = and = at x =. In the nmerical experiments the size of the grid blocks is given by x = and x c =. The grids are 2 xed since we here want to stdy the eect of varying the time step. Only one shock frontis present in the soltion, therefore we let the corresponding residal x term be dened on the whole interval [ ]. In the COS simlation we have sed time step to reach T =:2 since with the COS method the correct width of the front is obtained independently of the time step. In the OS simlations the nmber of time steps is and. We refer to these simlations as COS(), OS() and OS(), respectively. In Figre 3 the COS and OS soltions are shown together with the reference soltion, which is generated by a nite dierence scheme with very ne space and time discretization (t = T 48 To measre the relative errors in the (C)OS satration proles we dene (6) r (C)OS = jj~ ref ; ~ (C)OS jj 2 jj~ ref jj 2 : and x = 5 ). We have r COS() = :6, r OS() = :5 and r OS() = :36. Ths the relative errors measred in the L 2 -norm are abot eqally large for the COS() and OS() soltions. From Figre 3 and the relative errors it is seen that COS is the most ecient and accrate algorithm for the problem nder consideration. We refer to [5, 6] for a more detailed comparison of COS and OS for the forward problem. One may ask whether the relative merits of OS and COS will be similar also when the models are applied to solve an associated inverse problem. The importance of neglecting some specic physics may be dierent for the inverse problem than for the forward problem. 3. The inverse problem 3.. Problem formlation. We will attempt to recover the x fnction, f(), from discrete vales of, which are denoted by ~ obs. The forward model is given by eqation () (recall that d() ) with an appropriate soltion algorithm for (x t). The OS()
8 8 BRUSDAL ET AL..9 a).9 b) f().5 (x,) x Figre 2. The x fnction, f (), in plot a) and the initial satration, (x), in plot b) a) Reference (x) COS() OS() b) Reference (x) OS() x x Figre 3. The reference, COS() and OS() soltions of the forward problem in plot a). The reference soltion and OS() soltion in plot b). and COS() soltion algorithms are selected since these gave similar accracy when applied to the forward problem. The discrete observed vales will be sampled from the ne-grid soltion of the forward problem at time T = :2. The x fnction sed when calclating the ne-grid soltion is f () (seesection 2.2). Objective fnctions, which measre the least sqares distance between observed and corresponding calclated vales of ~, are formed J (C)OS (f) =k~ obs ;~ (C)OS (f)k 2 2. Recovery
9 OPERATOR SPLITTING AND CORRECTED OPERATOR SPLITTING 9 of f () is carried ot throgh minimisation of the objective fnctions. Representing the x fnction by a B-spline expansion leads to a nite dimensional parameter estimation problem for the coecients in the expansion. This is solved by a Levenberg-Marqardt algorithm [8] Soltion criteria. Dierent parametrisations of the x fnction lead to dierent parameter estimation problems with soltions corresponding to dierent fnctional spaces for the estimated x fnction. Hence, having solved the parameter estimation problem does not necessarily mean that one has solved the nderlying inverse problem. Frthermore, inverse problems are known to be inherently ill-posed. Soltion non-existence, soltion non-niqeness and instability of the soltion with respect to data may occr. This shows the need for soltion criteria to decide whether a soltion to the parameter estimation problem is also a soltion to the original inverse problem. Often, inverse problems where the modelling errors are negligible compared to the error in the observed data are considered. When the noise in the observed data is a reslt of additive independent random errors (e.g., reslting from independent measrements), it is possible to apply statistically based soltion criteria [, 2] to test if the soltion to a parameter estimation problem can be accepted as a soltion to the nderlying inverse problem. These criteria concern the vale of the objective fnction and the randomness of the elements in the residal vector at soltion. Frthermore, if a soltion to this kind of inverse problem has been fond, it is possible to give approximate error bonds for the soltion. The inverse problem stdied here is qite dierent. We consider a problem where the modelling error dominates the error in the observed data. In fact, the whole prpose is to stdy the eect of two dierent modelling errors on the recovery of the x fnction. Contrary to measrement errors, modelling errors can not be expected to be additive independent random errors. Hence, neither the above soltion criteria nor the error bonds can be applied. Presently, we do not know enogh abot the natre of the modelling error to pt forth any alternative soltion criteria covering this case. Therefore, the estimated x fnctions presented later will be soltions to specic parameter estimation problems, while we can not claim that any inverse problem has actally been solved. However, it is still of interest to stdy the inence of the modelling errors in OS and COS when attempting to recover the x fnction.
10 BRUSDAL ET AL. We will se the same parametrisations of the x fnction when rnning the OS() model as when rnning the COS() model, thereby ensring that corresponding OS and COS estimates belong to the same fnctional space Nmerical experiments. The observed data are generated at T = :2 by a nite dierence scheme sing the x fnction, f (), and the ne grid given in section 2.2. From the reslting satration prole eqally spaced data are sampled. This constittes the observation data sed in the parameter estimation. The size of the grid cells sed in OS() and COS() is xed x =. For a given order, m, of the B-spline basis fnctions and a given nmber of internal knots, k, in the knot vector, the nmber of coecients in each relative permeability crve is m + k. The nmber of elements in the knot vector is 2m + k. Here the order, m, is xed m = 3. The rst and last coecients in the B-spline expansions of each relative permeability crve are also xed, c w = c m+k o =: and c m+k w = c o =:, so that the total nmber of parameters to be estimated is 2(m + k) ; 4=2(k +). In the nmerical experiments, 2, 3 and 4, the exibility of the estimated fractional ow fnction is systematically increased by adding knots in the knot vector (increasing k). The nmber of knots is k =,k =2,k = 3 and k = 7, respectively. Increased exibility in the x fnction allows the x fnction to compensate for modelling errors in the OS and COS models. Hence, we expect the minimm vale of the objective fnction, J (C)OS (f), to decrease with increased exibility in the x fnction while the relative errors in the estimated x fnctions may increase. The COS and OS estimated x fnctions and the corresponding minimm vales of J (C)OS (f) willbe compared at each level of exibility. The knot vector, initial gess on the nknown parameters, and the estimated parameters from the experiments can be fond in Table, 2, 3 and 4, respectively. The minimm vales of the objective fnction, J (C)OS (f), and the relative errors in the estimated x fnctions are given in Table 5. In Figre 4b) the observed satrations and (C)OS simlated satrations corresponding to the minimm vale of J (C)OS are given for experiment. We omit plotting the satration proles for experiment 2, 3 and 4, since they are visally identical to the plot in Figre 4b). In Figre 5 and 6 the tre x fnction, f (), and the estimated x fnctions, f (C)OS (), are plotted.
11 OPERATOR SPLITTING AND CORRECTED OPERATOR SPLITTING.9 a).9 b).8.7 Flx fnction, f * () Initial f().8.7 Observed (x) COS() OS().6.6 f() x Figre 4. The tre x fnction, f (), and the initial gess on the x fnction in plot a). In plot b) the observed satration, obs, and simlated satrations, (C)OS,inexperiment..9 a).9 b) f().5 f() Tre flx fnction, f * () COS() OS().3.2 Tre flx fnction, f * () COS() OS() Figre 5. The tre x fnction, f (), and estimated x fnctions, f (C)OS (), in experiment (plot a)) and experiment 2 (plot b)). The minimm vales of the objective fnction, J (C)OS (f), are decreasing with increased exibility of the fractional ow fnction, whereas in most cases the relative error of the x fnction is increasing. We also have J COS() < J OS() and that the relative error in the x fnction is greater with OS() than COS() at each level of exibility. We interpret these reslts as follows: In the convection step the speed of a point ^ on a satration crve is given by f (^). The OS() soltion contains nmerical dision
12 2 BRUSDAL ET AL..9 a).9 b) f().5 f() Tre flx fnction, f * () COS() OS().2 Tre flx fnction, f * () COS() OS() Figre 6. The tre x fnction, f (), and estimated x fnctions, f (C)OS (), in experiment 3 (plot a)) and experiment 4 (plot b)). Table Experiment no. Knot vector ~y =(: : : :5 : : :) Initial parameters ~c w =(: :3 :82 :) ~c o =(: :45 :2 :) Estimated parameters ~c w =(: :98 :399 :) COS() ~c o =(: :66 :27 :) Estimated parameters ~c w =(: : :725 :) OS() ~c o =(: :326 :9 :) sch that low-satration points on the front travel too fast. To compensate for this the inversion algorithm seeks to decrease the speed of these points by deforming the estimated x fnction in the lower satration region, cf. Figre 5 and 6. The ability to do so, withot aecting the higher satration regions too mch, increases with increased exibility in the fnctional representation of the x fnction. Observe that not even the deformed x fnction in Figre 6b) did bring the minimm vale of J OS() below that of J COS() in experiment. Ths, the modelling errors in OS() are clearly more serios than those of COS(), as far as the crrent inverse problem is concerned. In an attempt to remove nmerical dision in the OS simlations, the nmber of time steps is increased to 2. With 2 time steps the width of the front in the OS soltion is of
13 OPERATOR SPLITTING AND CORRECTED OPERATOR SPLITTING 3 Table 2 Experiment no. 2 Knot vector ~y =(: : : :3 :5 : : :) Initial parameters ~c w =(: :8 :456 :82 :) ~c o =(: :67 :32 :2 :) Estimated parameters ~c w =(: : :35 :79 :) COS() ~c o =(: :562 :258 :92 :) Estimated parameters ~c w =(: : :29 :67 :) OS() ~c o =(: : :22 :2 :) Table 3 Experiment no. 3 Knot vector ~y =(: : : :3 :5 :7 : : :) Initial parameters ~c w =(: :8 :456 :664 :892 :) ~c o =(: :67 :32 :49 :2 :) Estimated parameters ~c w =(: : :25 :35 :52 :) COS() ~c o =(: :338 :52 :63 :4 :) Estimated parameters ~c w =(: : :2 :255 :587 :) OS() ~c o =(: : :66 :64 :4 :) Table 4 Experiment no. 4 Knot vector ~y =(: : : :5 : :2 :3 :5 :6 :7 : : :) Initial parameters ~c w =(: :3 :89 :77 :29 :456 :62 :7 :892 :) ~c o =(: :945 :842 :697 :53 :32 :92 :22 :2 :) Estimated parameters ~c w =(: : : :34 :89 :89 :269 :39 :477 :) COS() ~c o =(: : : :7 :4 :28 :8 :47 :39 :) Estimated parameters ~c w =(: : : : : : :296 :352 :587 :) OS() ~c o =(: : : : :4 :4 :94 :58 :34 :) the order of the physical front, ". We repeat the nmerical experiments, 2, 3 and 4 sing OS(2). The minimm vales of the objective fnction, J OS (f), and the relative errors of
14 4 BRUSDAL ET AL. Table 5 Experimental reslts Objective fnctions, J (C)OS (f) Experiment no COS() :96 :34 :27.23 OS() :25 :245 : OS(2) :2 :9 :82.8 Relative error in f (C)OS () Experiment no COS() :889 :252 : OS() :449 :728 : OS(2).365 :3 : x 3.9 a) b).7 2 f() Objective fnction.5 COS() OS() OS(2).2 Tre flx fnction, f * () OS() OS(2) Relative error Figre 7. The tre x fnction, f (), and OS() and OS(2) estimated x fnctions in plot a). The minimm vales of the objective fnction, J (C)OS (f), verss the relative errors in the estimated x fnction, f (C)OS (), in plot b). the estimated x fnctions have decreased compared with the OS() reslts (see Table 5 and Figre 7b)). Figre 7a) shows the tre x fnction, f (), and the OS() and OS(2) estimated x fnctions in experiment 3. The OS(2) estimated x fnction is considerably better than the one achieved with OS(). This spports or view that the deformation in the OS() estimated x fnctions is de to nmerical dision.
15 OPERATOR SPLITTING AND CORRECTED OPERATOR SPLITTING 5 4. Conclsion Wehave compared COS and OS models with respect to accracy in the context of solving the one-dimensional satration eqation of poros media ow. In line with previos work, the nmerical experiments indicate that the COS method is mch more accrate than the OS method. For the problem considered here the OS approach needed time steps to achieve the same accracy in the L 2 -norm as COS did with timestep. We have also stdied the COS and OS models in the context of an inverse problem where we recover the x fnction in the one-dimensional satration eqation from discrete satration data. The nmerical experiments indicate that even thogh the COS() and OS() models solve the forward problem with similar accracy, and ths in that sense contain modelling errors of eqal size, the modelling error in OS() is the most serios one for this inverse problem. References [] Y. Bard. Nonlinear Parameter Estimation. John Wiley and Sons, New York City, 98. [2] J. T. Beale and A. Majda. Rates of convergence for viscos splitting of the Navier-Stokes eqations. Math. Comp., 37(56):243{259, 98. [3] G. Chavent and J. Jare. Mathematical models and nite elements for reservoir simlation, volme 7 of Stdies in mathematics and it's applications. North Holland, Amsterdam, 986. [4] H. K. Dahle. Adaptive characteristic operator splitting techniqes for convection-dominated dision problems in one and two space dimensions. PhD thesis, Department of Mathematics, University of Bergen, 988. [5] H. K. Dahle. ELLAM-based operator splitting for nonlinear advection dision eqations. Technical Report 98, Department of Mathematics, University of Bergen, 995. [6] H. K. Dahle, M. S. Espedal, R. E. Ewing, and O. Svareid. Characteristic adaptive sbdomain methods for reservoir ow problems. Nmerical Methods for PDEs, 6:279{39, 99. [7] H. K. Dahle, M. S. Espedal, and O. Svareid. Characteristic, local grid renement techniqes for reservoir ow problems. International J. for Nmerical Methods in Engineering, 34:5{69, 992. [8] H. K. Dahle, R. E. Ewing, and T. F. Rssell. Elerian-Lagrangian localized adjoint methods for a nonlinear advection-dision eqation. Compt. Methods. Appl. Mech. Engrg., 22:223{25, 995. [9] C. N. Dawson. Godnov-mixed methods for advective ow problems in one space dimension. SIAM J. Nm. Anal., 28(5):282{39, Oct. 99. [] M. S. Espedal and R. E. Ewing. Characteristic Petrov-Galerkin sbdomain methods for two-phase immiscible ow. Compt. Methods. Appl. Mech. Engrg., 64:3{35, 987.
16 6 BRUSDAL ET AL. [] S. Evje and K. H. Karlsen. Viscos splitting approximation of mixed hyperbolic-parabolic convection-dision eqations. Report No. 2, 997/98. Institt Mittag-Leer. Available at the URL [2] R. E. Ewing. Operator splitting and Elerian-Lagrangian localized adjoint methods for mltiphase ow. In Whiteman, editor, The Mathematics of Finite Elements and Applications VII MAFELAP, pages 25{232. Academic press, San Diego, CA, 99. [3] J. B. Goodman and R. J. LeVeqe. A geometric approach to high resoltion tvd schemes. SIAM J. Nm. Anal., 25:268{284, 988. [4] H. Holden, L. Holden, and R. Hegh-Krohn. A nmerical method for rst order nonlinear scalar conservation laws in one-dimension. Compt. Math. Applic., 5(6{8):595{62, 988. [5] K. H. Karlsen, K. Brsdal, H. K. Dahle, S. Evje, and K.-A. Lie. The corrected operator splitting approach applied to a nonlinear advection-dision problem. To appear in Compt. Methods. Appl. Mech. Engrg. [6] K. H. Karlsen and N. H. Risebro. Corrected operator splitting for nonlinear parabolic eqations. Preprint, University of Bergen, Norway, 997. Available at the URL [7] K. H. Karlsen and N. H. Risebro. An operator splitting method for nonlinear convection-dision eqations. Nmer. Math., 77(3):365{382, 997. [8] D.. W. Marqardt. An algorithm for least sqares estimation of nonlinear parameters. SIAM J. Appl. Math., :43{44, 963. [9] L. L. Schmaker. Spline Fnctions: Basic Theory. John Wiley and Sons, New York, 98. [2] S. Weisberg. Applied Linear Regression. John Wiley and Sons, New York, 985. (Kari Brsdal) Department of Mathematics, University of Bergen, Johs. Brnsgt. 2, N-58 Bergen, Norway address: kari.brsdal@mi.ib.no (Helge K. Dahle) Department of Mathematics, University of Bergen, Johs. Brnsgt. 2, N-58 Bergen, Norway address: helge.dale@mi.ib.no (Kenneth Hvistendahl Karlsen) Department of Mathematics, University of Bergen, Johs. Brnsgt. 2, N-58 Bergen, Norway address: kenneth.karlsen@mi.ib.no (Trond Mannseth) RF Rogaland Research,, Thormhlensgt. 55, N-58 Bergen, Norway address: trond.mannseth@rf.no
A study of the modelling error in two operator splitting algorithms for porous media flow
Computational Geosciences 2 (1998) 23 36 23 A study of the modelling error in two operator splitting algorithms for porous media flow K. Brusdal a,h.k.dahle a, K. Hvistendahl Karlsen a and T. Mannseth
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More information4 Exact laminar boundary layer solutions
4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved
More informationA Survey of the Implementation of Numerical Schemes for Linear Advection Equation
Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear
More informationSubcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany
Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2
More informationHomotopy Perturbation Method for Solving Linear Boundary Value Problems
International Jornal of Crrent Engineering and Technolog E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rights Reserved Available at http://inpressco.com/categor/ijcet Research Article Homotop
More informationarxiv: v1 [physics.flu-dyn] 4 Sep 2013
THE THREE-DIMENSIONAL JUMP CONDITIONS FOR THE STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY, SINGULAR FORCES, AND AN INCOMPRESSIBLE INTERFACE PRERNA GERA AND DAVID SALAC arxiv:1309.1728v1 physics.fl-dyn]
More informationThe spreading residue harmonic balance method for nonlinear vibration of an electrostatically actuated microbeam
J.L. Pan W.Y. Zh Nonlinear Sci. Lett. Vol.8 No. pp.- September The spreading reside harmonic balance method for nonlinear vibration of an electrostatically actated microbeam J. L. Pan W. Y. Zh * College
More informationDetermining of temperature field in a L-shaped domain
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 0, (:-8 Determining of temperatre field in a L-shaped domain Oigo M. Zongo, Sié Kam, Kalifa Palm, and Alione Oedraogo
More informationAffine Invariant Total Variation Models
Affine Invariant Total Variation Models Helen Balinsky, Alexander Balinsky Media Technologies aboratory HP aboratories Bristol HP-7-94 Jne 6, 7* Total Variation, affine restoration, Sobolev ineqality,
More informationNonlinear parametric optimization using cylindrical algebraic decomposition
Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationCurves - Foundation of Free-form Surfaces
Crves - Fondation of Free-form Srfaces Why Not Simply Use a Point Matrix to Represent a Crve? Storage isse and limited resoltion Comptation and transformation Difficlties in calclating the intersections
More informationA Single Species in One Spatial Dimension
Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,
More informationDirect linearization method for nonlinear PDE s and the related kernel RBFs
Direct linearization method for nonlinear PDE s and the related kernel BFs W. Chen Department of Informatics, Uniersity of Oslo, P.O.Box 1080, Blindern, 0316 Oslo, Norway Email: wenc@ifi.io.no Abstract
More informationChapter 2 Difficulties associated with corners
Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces
More informationOptimal Control, Statistics and Path Planning
PERGAMON Mathematical and Compter Modelling 33 (21) 237 253 www.elsevier.nl/locate/mcm Optimal Control, Statistics and Path Planning C. F. Martin and Shan Sn Department of Mathematics and Statistics Texas
More informationA Macroscopic Traffic Data Assimilation Framework Based on Fourier-Galerkin Method and Minimax Estimation
A Macroscopic Traffic Data Assimilation Framework Based on Forier-Galerkin Method and Minima Estimation Tigran T. Tchrakian and Sergiy Zhk Abstract In this paper, we propose a new framework for macroscopic
More informationA fundamental inverse problem in geosciences
A fndamental inverse problem in geosciences Predict the vales of a spatial random field (SRF) sing a set of observed vales of the same and/or other SRFs. y i L i ( ) + v, i,..., n i ( P)? L i () : linear
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationActive Flux Schemes for Advection Diffusion
AIAA Aviation - Jne, Dallas, TX nd AIAA Comptational Flid Dynamics Conference AIAA - Active Fl Schemes for Advection Diffsion Hiroaki Nishikawa National Institte of Aerospace, Hampton, VA 3, USA Downloaded
More informationOptimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance
Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for
More informationExperimental Study of an Impinging Round Jet
Marie Crie ay Final Report : Experimental dy of an Impinging Rond Jet BOURDETTE Vincent Ph.D stdent at the Rovira i Virgili University (URV), Mechanical Engineering Department. Work carried ot dring a
More informationIncompressible Viscoelastic Flow of a Generalised Oldroyed-B Fluid through Porous Medium between Two Infinite Parallel Plates in a Rotating System
International Jornal of Compter Applications (97 8887) Volme 79 No., October Incompressible Viscoelastic Flow of a Generalised Oldroed-B Flid throgh Poros Medim between Two Infinite Parallel Plates in
More informationA New Approach to Direct Sequential Simulation that Accounts for the Proportional Effect: Direct Lognormal Simulation
A ew Approach to Direct eqential imlation that Acconts for the Proportional ffect: Direct ognormal imlation John Manchk, Oy eangthong and Clayton Detsch Department of Civil & nvironmental ngineering University
More informationSuyeon Shin* and Woonjae Hwang**
JOURNAL OF THE CHUNGCHEONG MATHEMATICAL SOCIETY Volme 5, No. 3, Agst THE NUMERICAL SOLUTION OF SHALLOW WATER EQUATION BY MOVING MESH METHODS Syeon Sin* and Woonjae Hwang** Abstract. Tis paper presents
More informationTheoretical and Experimental Implementation of DC Motor Nonlinear Controllers
Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,
More informationStudy of the diffusion operator by the SPH method
IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait
More informationJ.A. BURNS AND B.B. KING redced order controllers sensors/actators. The kernels of these integral representations are called fnctional gains. In [4],
Jornal of Mathematical Systems, Estimation, Control Vol. 8, No. 2, 1998, pp. 1{12 c 1998 Birkhaser-Boston A Note on the Mathematical Modelling of Damped Second Order Systems John A. Brns y Belinda B. King
More informationTIME ACCURATE FAST THREE-STEP WAVELET-GALERKIN METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS
International Jornal of Wavelets, Mltiresoltion and Information Processing Vol. 4, No. (26) 65 79 c World Scientific Pblishing Company TIME ACCURATE FAST THREE-STEP WAVELET-GALERKIN METHOD FOR PARTIAL
More informationTechnical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty
Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte
More informationIn the seqel we also sppose that the range R(A) of the operator A is non-closed, in general, and then the main eqation nder consideration, A = f ; (1.
The Lavrentiev-reglarized Galerkin method for linear accretive ill-posed problems Robert Plato z Fachbereich Mathematik, Technische Universit t Berlin, Stra e des 17. Jni 135, D10623 Berlin, Germany Abstract
More informationPC1 PC4 PC2 PC3 PC5 PC6 PC7 PC8 PC
Accracy verss interpretability in exible modeling: implementing a tradeo sing Gassian process models Tony A. Plate (tap@mcs.vw.ac.nz) School of Mathematical and Compting Sciences Victoria University of
More informationPIPELINE MECHANICAL DAMAGE CHARACTERIZATION BY MULTIPLE MAGNETIZATION LEVEL DECOUPLING
PIPELINE MECHANICAL DAMAGE CHARACTERIZATION BY MULTIPLE MAGNETIZATION LEVEL DECOUPLING INTRODUCTION Richard 1. Davis & 1. Brce Nestleroth Battelle 505 King Ave Colmbs, OH 40201 Mechanical damage, cased
More informationSources of Non Stationarity in the Semivariogram
Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary
More informationLecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2
BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE
More informationEVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE
13 th World Conference on Earthqake Engineering Vancover, B.C., Canada Agst 1-6, 2004 Paper No. 3099 EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE Ellen M. RATHJE 1, Wen-Jong CHANG 2, Kenneth
More informationDiscussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli
1 Introdction Discssion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen Department of Economics, University of Copenhagen and CREATES,
More informationAnalytical Investigation of Hyperbolic Equations via He s Methods
American J. of Engineering and Applied Sciences (4): 399-47, 8 ISSN 94-7 8 Science Pblications Analytical Investigation of Hyperbolic Eqations via He s Methods D.D. Ganji, M. Amini and A. Kolahdooz Department
More informationResearch Article Permanence of a Discrete Predator-Prey Systems with Beddington-DeAngelis Functional Response and Feedback Controls
Hindawi Pblishing Corporation Discrete Dynamics in Natre and Society Volme 2008 Article ID 149267 8 pages doi:101155/2008/149267 Research Article Permanence of a Discrete Predator-Prey Systems with Beddington-DeAngelis
More informationTed Pedersen. Southern Methodist University. large sample assumptions implicit in traditional goodness
Appears in the Proceedings of the Soth-Central SAS Users Grop Conference (SCSUG-96), Astin, TX, Oct 27-29, 1996 Fishing for Exactness Ted Pedersen Department of Compter Science & Engineering Sothern Methodist
More informationA Regulator for Continuous Sedimentation in Ideal Clarifier-Thickener Units
A Reglator for Continos Sedimentation in Ideal Clarifier-Thickener Units STEFAN DIEHL Centre for Mathematical Sciences, Lnd University, P.O. Box, SE- Lnd, Sweden e-mail: diehl@maths.lth.se) Abstract. The
More informationDEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS (H/U < 20) AND CONSEQUENCES ON CRITICALITY SAFETY
DEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS ( < 20) AND CONSEQUENCES ON CRITICALITY SAFETY N. Leclaire, S. Evo, I.R.S.N., France Introdction In criticality stdies, the blk density
More informationInvestigation of Haar Wavelet Collocation Method to Solve Ninth Order Boundary Value Problems
Global Jornal of Pre and Applied Mathematics. ISSN 0973-1768 Volme 13, Nmber 5 (017), pp. 1415-148 Research India Pblications http://www.ripblication.com Investigation of Haar Wavelet Collocation Method
More informationStudy on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom
EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical
More informationA Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations
Applied Mathematics, 05, 6, 04-4 Pblished Online November 05 in SciRes. http://www.scirp.org/jornal/am http://d.doi.org/0.46/am.05.685 A Comptational Stdy with Finite Element Method and Finite Difference
More informationA sixth-order dual preserving algorithm for the Camassa-Holm equation
A sith-order dal preserving algorithm for the Camassa-Holm eqation Pao-Hsing Chi Long Lee Tony W. H. She November 6, 29 Abstract The paper presents a sith-order nmerical algorithm for stdying the completely
More informationFOUNTAIN codes [3], [4] provide an efficient solution
Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn
More informationIntrodction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc
A New Architectre for a Parallel Finite Field Mltiplier with Low Complexity Based on Composite Fields Christof Paar y IEEE Transactions on Compters, Jly 996, vol 45, no 7, pp 856-86 Abstract In this paper
More informationFluid-induced rotordynamic forces produced by the uid in an annular seal or in the leakage
Flid Flow Eqations For Rotordynamic Flows In Seals And Leakage Paths C.E. Brennen R.V.Uy California Institte of Technology, Pasadena, California 95 ABSTRACT Flid-indced rotordynamic forces prodced by the
More informationE ect Of Quadrant Bow On Delta Undulator Phase Errors
LCLS-TN-15-1 E ect Of Qadrant Bow On Delta Undlator Phase Errors Zachary Wolf SLAC Febrary 18, 015 Abstract The Delta ndlator qadrants are tned individally and are then assembled to make the tned ndlator.
More informationThe Linear Quadratic Regulator
10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.
More informationDynamics of a Holling-Tanner Model
American Jornal of Engineering Research (AJER) 07 American Jornal of Engineering Research (AJER) e-issn: 30-0847 p-issn : 30-0936 Volme-6 Isse-4 pp-3-40 wwwajerorg Research Paper Open Access Dynamics of
More informationEffects of Soil Spatial Variability on Bearing Capacity of Shallow Foundations
Geotechnical Safety and Risk V T. Schweckendiek et al. (Eds.) 2015 The athors and IOS Press. This article is pblished online with Open Access by IOS Press and distribted nder the terms of the Creative
More informationON THE PERFORMANCE OF LOW
Monografías Matemáticas García de Galdeano, 77 86 (6) ON THE PERFORMANCE OF LOW STORAGE ADDITIVE RUNGE-KUTTA METHODS Inmaclada Higeras and Teo Roldán Abstract. Gien a differential system that inoles terms
More informationSafe Manual Control of the Furuta Pendulum
Safe Manal Control of the Frta Pendlm Johan Åkesson, Karl Johan Åström Department of Atomatic Control, Lnd Institte of Technology (LTH) Box 8, Lnd, Sweden PSfrag {jakesson,kja}@control.lth.se replacements
More informationDiscussion Papers Department of Economics University of Copenhagen
Discssion Papers Department of Economics University of Copenhagen No. 10-06 Discssion of The Forward Search: Theory and Data Analysis, by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen,
More informationBASIC RELATIONS BETWEEN TSUNAMIS CALCULATION AND THEIR PHYSICS II
BASIC RELATIONS BETWEEN TSUNAMIS CALCULATION AND THEIR PHYSICS II Zygmnt Kowalik Institte of Marine Science, University of Alaska Fairbanks, AK 99775, USA ABSTRACT Basic tsnami physics of propagation and
More informationEfficiency Increase and Input Power Decrease of Converted Prototype Pump Performance
International Jornal of Flid Machinery and Systems DOI: http://dx.doi.org/10.593/ijfms.016.9.3.05 Vol. 9, No. 3, Jly-September 016 ISSN (Online): 188-9554 Original Paper Efficiency Increase and Inpt Power
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationSubstructure Finite Element Model Updating of a Space Frame Structure by Minimization of Modal Dynamic Residual
1 Sbstrctre Finite Element Model Updating of a Space Frame Strctre by Minimization of Modal Dynamic esidal Dapeng Zh, Xinn Dong, Yang Wang, Member, IEEE Abstract This research investigates a sbstrctre
More informationIntrodction In this paper we develop a local discontinos Galerkin method for solving KdV type eqations containing third derivative terms in one and ml
A Local Discontinos Galerkin Method for KdV Type Eqations Je Yan and Chi-Wang Sh 3 Division of Applied Mathematics Brown University Providence, Rhode Island 09 ABSTRACT In this paper we develop a local
More informationApplying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB
IOS Jornal of Mathematics (IOS-JM) e-issn: 78-578, p-issn: 319-765X. Volme 13, Isse 6 Ver. II (Nov. - Dec. 17), PP 5-59 www.iosrjornals.org Applying Laminar and Trblent Flow and measring Velocity Profile
More informationAn optimizing reduced order FDS for the tropical Pacific Ocean reduced gravity model
An optimizing redced order FDS for the tropical Pacific Ocean redced gravity model Zhendong Lo a,, Jing Chen b,, Jiang Zh c,, Riwen Wang c,, and I. M. Navon d,, a School of Science, Beijing Jiaotong University,
More informationStep-Size Bounds Analysis of the Generalized Multidelay Adaptive Filter
WCE 007 Jly - 4 007 London UK Step-Size onds Analysis of the Generalized Mltidelay Adaptive Filter Jnghsi Lee and Hs Chang Hang Abstract In this paper we analyze the bonds of the fixed common step-size
More informationrate of convergence. (Here " is the thickness of the layer and p is the degree of the approximating polynomial.) The importance of the reslts in 13 li
athematical odels and ethods in Applied Sciences fc World Scientic Pblishing Company Vol. 3, o.8 (1998) THE hp FIITE ELEET ETHOD FOR SIGULARLY PERTURBED PROES I SOOTH DOAIS CHRISTOS A. XEOPHOTOS Department
More informationAbstract. A front tracking method is used to construct weak solutions to
A Front Tracking Method for Conservation Laws with Boundary Conditions K. Hvistendahl Karlsen, K.{A. Lie, and N. H. Risebro Abstract. A front tracking method is used to construct weak solutions to scalar
More information), σ is a parameter, is the Euclidean norm in R d.
WEIGHTED GRAPH LAPLACIAN AND IMAGE INPAINTING ZUOQIANG SHI, STANLEY OSHER, AND WEI ZHU Abstract. Inspired by the graph Laplacian and the point integral method, we introdce a novel weighted graph Laplacian
More informationMODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION
MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION Rssell Qadros, Krishnend Sinha Department of Aerospace Engineering Indian Institte of Technology Bombay Mmbai, India 476 Johan
More informationNumerical Simulation of Three Dimensional Flow in Water Tank of Marine Fish Larvae
Copyright c 27 ICCES ICCES, vol.4, no.1, pp.19-24, 27 Nmerical Simlation of Three Dimensional Flo in Water Tank of Marine Fish Larvae Shigeaki Shiotani 1, Atsshi Hagiara 2 and Yoshitaka Sakakra 3 Smmary
More informationUNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL
8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a
More informationOn Multiobjective Duality For Variational Problems
The Open Operational Research Jornal, 202, 6, -8 On Mltiobjective Dality For Variational Problems. Hsain *,, Bilal Ahmad 2 and Z. Jabeen 3 Open Access Department of Mathematics, Jaypee University of Engineering
More informationApproximate Solution of Convection- Diffusion Equation by the Homotopy Perturbation Method
Gen. Math. Notes, Vol. 1, No., December 1, pp. 18-114 ISSN 19-7184; Copyright ICSRS Pblication, 1 www.i-csrs.org Available free online at http://www.geman.in Approximate Soltion of Convection- Diffsion
More informationSTABILISATION OF LOCAL PROJECTION TYPE APPLIED TO CONVECTION-DIFFUSION PROBLEMS WITH MIXED BOUNDARY CONDITIONS
STABILISATION OF LOCAL PROJECTION TYPE APPLIED TO CONVECTION-DIFFUSION PROBLEMS WITH MIXED BOUNDARY CONDITIONS GUNAR MATTHIES, PIOTR SRZYPACZ, AND LUTZ TOBISA Abstract. We present the analysis for the
More informationPHASE STEERING AND FOCUSING BEHAVIOR OF ULTRASOUND IN CEMENTITIOUS MATERIALS
PHAS STRING AND FOCUSING BHAVIOR OF ULTRASOUND IN CMNTITIOUS MATRIALS Shi-Chang Wooh and Lawrence Azar Department of Civil and nvironmental ngineering Massachsetts Institte of Technology Cambridge, MA
More informationAn Investigation into Estimating Type B Degrees of Freedom
An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information
More informationInternational Journal of Physical and Mathematical Sciences journal homepage:
64 International Jornal of Physical and Mathematical Sciences Vol 2, No 1 (2011) ISSN: 2010-1791 International Jornal of Physical and Mathematical Sciences jornal homepage: http://icoci.org/ijpms PRELIMINARY
More informationDownloaded 01/04/14 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 46 No. pp. 996 1011 c 008 Society for Indstrial and Applied Mathematics Downloaded 01/04/14 to 147.10.110.133. Redistribtion sbject to SIAM license or copyright; see http://www.siam.org/jornals/ojsa.php
More informationEXPT. 5 DETERMINATION OF pk a OF AN INDICATOR USING SPECTROPHOTOMETRY
EXPT. 5 DETERMITIO OF pk a OF IDICTOR USIG SPECTROPHOTOMETRY Strctre 5.1 Introdction Objectives 5.2 Principle 5.3 Spectrophotometric Determination of pka Vale of Indicator 5.4 Reqirements 5.5 Soltions
More informationEffects of modifications on the hydraulics of Denil fishways
BOREAL ENVIRONMENT RESEARCH 5: 67 79 ISSN 1239-6095 Helsinki 28 March 2000 2000 Effects of modifications on the hydralics of Denil fishways Riitta Kamla 1) and Jan Bärthel 2) 1) Water Resorces and Environmental
More informationDownloaded 07/06/18 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 4, No., pp. A4 A7 c 8 Society for Indstrial and Applied Mathematics Downloaded 7/6/8 to 8.83.63.. Redistribtion sbject to SIAM license or copyright; see http://www.siam.org/jornals/ojsa.php
More informationDesign and Data Acquisition for Thermal Conductivity Matric Suction Sensors
68 TRANSPORTATION RSARCH RCORD 1432 Design and Data Acqisition for Thermal Condctivity Matric Sction Sensors J. K.-M. GAN, D. G. FRDLUND, A. XING, AND W.-X. LI The principles behind sing the thermal condctivity
More informationChapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS
Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system
More information08.06 Shooting Method for Ordinary Differential Equations
8.6 Shooting Method for Ordinary Differential Eqations After reading this chapter, yo shold be able to 1. learn the shooting method algorithm to solve bondary vale problems, and. apply shooting method
More informationNetwork Coding for Multiple Unicasts: An Approach based on Linear Optimization
Network Coding for Mltiple Unicasts: An Approach based on Linear Optimization Danail Traskov, Niranjan Ratnakar, Desmond S. Ln, Ralf Koetter, and Mriel Médard Abstract In this paper we consider the application
More informationAssignment Fall 2014
Assignment 5.086 Fall 04 De: Wednesday, 0 December at 5 PM. Upload yor soltion to corse website as a zip file YOURNAME_ASSIGNMENT_5 which incldes the script for each qestion as well as all Matlab fnctions
More informationRESGen: Renewable Energy Scenario Generation Platform
1 RESGen: Renewable Energy Scenario Generation Platform Emil B. Iversen, Pierre Pinson, Senior Member, IEEE, and Igor Ardin Abstract Space-time scenarios of renewable power generation are increasingly
More informationMagnitude-Recurrence Relationship and its Effect on Uniform Hazard Spectra: A Current Assessment
Magnitde-Recrrence Relationship and its Effect on Uniform Hazard Spectra: A Crrent Assessment T. K. Sen PhD (London), CEng., FICE Granherne, UK, a ssidiary of Kellogg Brown & Root Ltd SUMMARY: Gtenerg
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationarxiv: v1 [hep-ph] 17 Oct 2011
Extrapolation Algorithms for Infrared Divergent Integrals arxiv:1110.587v1 [hep-ph] 17 Oct 2011 Western Michigan University E-mail: elise.dedoncker@wmich.ed Jnpei Fjimoto High Energy Accelerator Research
More informationMechanisms and topology determination of complex chemical and biological network systems from first-passage theoretical approach
Mechanisms and topology determination of complex chemical and biological network systems from first-passage theoretical approach Xin Li and Anatoly B. Kolomeisky Citation: J. Chem. Phys. 39, 4406 (203);
More informationOn the Optimization of Numerical Dispersion and Dissipation of Finite Difference Scheme for Linear Advection Equation
Applied Mathematical Sciences, Vol. 0, 206, no. 48, 238-2389 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/ams.206.6463 On the Optimization of Nmerical Dispersion and Dissipation of Finite Difference
More informationThe Dual of the Maximum Likelihood Method
Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All
More informationFormal Methods for Deriving Element Equations
Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact
More informationSolving a Class of PDEs by a Local Reproducing Kernel Method with An Adaptive Residual Subsampling Technique
Copyright 25 Tech Science Press CMES, vol.8, no.6, pp.375-396, 25 Solving a Class of PDEs by a Local Reprodcing Kernel Method with An Adaptive Residal Sbsampling Techniqe H. Rafieayan Zadeh, M. Mohammadi,2
More informationPrediction of Effective Asphalt Layer Temperature
TRANSPORTATION RESEARCH RECORD 1473 93 Prediction of Effective Asphalt Layer Temperatre EARL H. INGE, JR., AND Y. RICHARD KIM The most widely sed method for evalating deflection measrements for overlay
More informationDevelopment of Second Order Plus Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model
Development of Second Order Pls Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model Lemma D. Tfa*, M. Ramasamy*, Sachin C. Patwardhan **, M. Shhaimi* *Chemical Engineering Department, Universiti
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 02-01 The Krylov accelerated SIMPLER method for incompressible flow C. Vik and A. Saghir ISSN 1389-6520 Reports of the Department of Applied Mathematical Analysis
More informationMulti-Voltage Floorplan Design with Optimal Voltage Assignment
Mlti-Voltage Floorplan Design with Optimal Voltage Assignment ABSTRACT Qian Zaichen Department of CSE The Chinese University of Hong Kong Shatin,N.T., Hong Kong zcqian@cse.chk.ed.hk In this paper, we stdy
More information