CALIBRATION OF RIVER BED ROUGHNESS

Size: px
Start display at page:

Download "CALIBRATION OF RIVER BED ROUGHNESS"

Transcription

1 CALIBRATION OF RIVER BED ROUGHNESS A. M. Wasantha La 1, M. ASCE Singuar vaue decomposition (SVD) is used to caibrate the Manning s roughness coefficients in a 1-D unsteady fow mode of the Upper Niagara River. The method is used to sove for the parameters after formuating the caibration probem as a generaized inear inverse probem. SVD is usefu in soving under-determined, overdetermined or even-determined probems, and can provide information to compute matrices describing parameter resoution, covariance and correation. This information is usefu in identifying the important parameter groups in the mode. Caibration is repeated with different numbers of parameter groups to determine the variation of the output error and uncertainty of the parameters with the parameter dimension. For purposes of comparison, the mode with a seected group of parameters is caibrated using Gauss-Newton method and minimax methods. The study shows the reationship of the parameters to the geometric ayout of the river and the gaging stations. INTRODUCTION Probems invoving mathematica modes are categorized as direct probems and inverse probems. In direct probems, everything about the mode is known, and the objective is to find the system output for a given input. Inverse probems are cassified as identification, detection and state reconstruction probems. Caibration is an inverse probem associated with identification, and is used to determine unknown constants or parameters in a mode. Direct or expicit parameter determination is not possibe in many noninear probems, incuding unsteady river fow networks. Indirect methods based on an output error criterion are commony used to determine bed roughness parameters iterativey. Agorithms used in caibration generay incude methods based on mathematica programming, optima contro or tria and error techniques. A summary of caibration methods used in groundwater probems is given in the text by Wiis and Yeh (1987). Athough there is theoreticay an infinite number of spatiay dependent 1 Senior Civi Engineer, South Forida Water Management District, 3301 Gun Cub Rd., West Pam Beach, FL 33416; formay Research Engineer, Dept. of Civi and Env. Engg., Carkson University, Potsdam, NY

2 bed roughness parameters in a river network, ony a finite number exists in a computer mode due to spatia discretization used with finite difference or finite eement methods. Out of this finite number, ony a imited number can be identified in most inverse probems (Beck, 1987). The current study investigates the use of Singuar Vaue Decomposition (SVD) to study parameters in a river network mode. The resuts are aso used to form parameter groups to simpify the probem and avoid getting a severey underdetermined systems. The study shows that some of the parameter groups can aso be associated with the geometric ayout of the river network and the pacement of observation stations. When identifiabiity is not a probem, popuar methods of caibration used in the past for open channe probems incude methods by Becker and Yeh (1972), Yeh and Becker (1973), and Fread, et a. (1978). Reated appications incude caibration of surface irrigation parameters by Katopodes et a. (1990). Error criteria commony used in the caibration are based on minimization of the sum of square of the errors in cacuations, the maximum absoute error or the bias. Fread (1978) used a method based on minimization of the bias, which required breaking down of the river into a number of singe channe reaches before caibrating each reach separatey. The Kaman fitering method used by Chiu, et a. (1978) is aso a vauabe too in caibration, adopted from optima contro theory. The SVD method has been previousy used by Wiggins (1972) and Uhrhammer (1980) to caibrate seismoogic parameters. Text by Meneke (1984) gives a description of inverse methods used in soving geophysica probems. The objectives of this study incudes determination of the structure, accuracy, and the optimum parameter dimension in the river network mode. Parameter covariance and correation matrices are aso determined in the study. The caibration is repeated using optimization methods based on the Gauss-Newton method and the minimax methods for comparison of resuts. 2

3 THEORETICAL CONSIDERATIONS The theoretica derivation of the equations required for caibration assumes that the numerica mode to be caibrated can be represented as a continuous function of the mode parameters ocay within the usefu range of the variabes. If y k j are the vaues of state variabes or water eves simuated by the mode at times k at ocations j, and if Y k j same variabes, the objective of the caibration is to find parameters x i, i are the physica observations of the 1 2 n that minimize the errors ε k j y k j Y k j, j 1 2 m k 1 2, using a specified criterion. n number of parameters; m number of observation stations; and number of time steps. The infuence coefficient method (eg., Becker and Yeh, 1972) uses parameter perturbations to define infuence or sensitivity coefficients. a k ji y k j x i y k j im x i 0 x i x i y k j x i for each i 1 2 n j 1 2 m k 1 2 (1) in which, a k ji sensitivity of the j th water eve with respect to i th parameter at time k; x i i th parameter; k time step number; number of time steps. Using matrix notation, (1) can aso be expressed as in which, the matrix A k m n time step. y k m 1 A k m n x k n 1 (2) is caed the infuence matrix. The superscript k is used to identify the Numerica computation of infuence coefficients is carried out by perturbing each of the parameters by a sma amount x i. For the current caibration of Manning s coefficients, x is If n parameters are caibrated, the mode is run n 1 times to obtain a the m n eements of the sensitivity matrix. Generaized inear inverse probem The method used by Wiggins (1972), Ward et a. (1973) and Uhrhammer (1980) to sove inear inverse probems is used in the current study. Since there is a time series of errors at every gage, bias is used as the error indicator, and there are m error indicators for the probem. ε s j 1 y k j Yj k for j 1 2 m (3) k 1 3

4 in which, ε s j average error or bias at gage j, which is ε s in vector form. The corrections in parameters x needed to make ε s 0 can be determined by equating average y of (2) to -ε s and soving for et a. 1973). x. The resuting system of inear equations can be expressed as (Wiggins 1972, Ward 1 k 1 or, in matrix form, in which, a k jp x p 1 y k j Yj k for j 1 2 m p 1 2 n (4) k 1 A s 1 k 1 A s x ε s (5) A k ε s 1 The corrected set of parameters to be used in the next iteration are k 1 ε k (6) x r 1 x r ρ r x r (7) subjected to x r b x r 1 x r ub (8) in which r is the iteration number; ρ r parameter used to contro the step size; ub, b are specified upper and ower bound vaues. ρ r continued unti ε s 2 converges, and x 2 becomes sma 0 8 is used in the current study. The iterative procedure is x x δ in which δ machine precision of the computer or a arger vaue depending on the required parameter toerance. The m equations shown in (5) can be soved easiy using methods such as Gaussian eimination if m n and deta 0. When m n, the probem is overdetermined, and the east squares method eads to the soution x A T A 1 A T ε s as ong as deta T A 0. When m n, the probem is underdetermined, and the generaized or Penrose inverse soution is x A T AA T 1 ε s as ong as detaa T 0 (Nobe and Danie, 1975). None of these soution techniques can be used with rank deficient or i-conditioned matrices. The sensitivity matrix becomes rank deficient if there is at east one parameter that has no significant infuence on the any of the seected observations, or if there is at east one observation station that is not sufficienty affected by any of the parameters. Singuar Vaue Decomposition (SVD) is capabe of soving under, over, even or mixed determined probems as expained in detai by Meneke (1984). 4

5 SVD is based on the key theorem that a matrix A s can be decomposed into three matrices V, Λ and U such that A s UΛV T (9) in which U and V are m m and n n matrices of orthogona singuar vectors and Λ m n diagona matrix of singuar vaues of A s. Texts by Nobe and Danie (1975), Forsythe et a. (1977) and Press, et at., (1989) give detaied information about the method. With symmetric matrices, SVD gives eigenvaues and eigenvectors. The physica meaning of the terms and expressions of the SVD method are expained using (9) and (10). SVD essentiay diagonaizes the sensitivity matrix. After the decomposition, singuar vectors of ength n formed by the coumns of V give coefficients of the inear combinations of od parameters giving rise to new independent parameter groups. Singuar vectors of ength m formed by coumns of U give coefficients of the inear combinations of observations forming new observations groups. The new independent parameter groups are reated to groups of observations through the diagona eements in Λ. Both newy formed observation and parameter groups are independent of each other because of the orthogonaity of matrices V and U. Singuar vaues forming the diagona sensitivity matrix Λ reate these parameter groups to observations groups. The number of nonzero singuar vaues q in Λ is the rank of A s. It gives the maximum number of independent parameter groups that can be identified in the mode. SVD can be used to sove the system of equations shown in (5). The soution method repaces the previousy mentioned methods for m n, m n, m n or mixed determined probems. It works even if the matrices invoved are singuar (Meneke, 1984). The soution to the system of equations shown in (5) is obtained by decomposing A s into A s UΛV T. (5) then becomes Λ z d (10) in which, z V T x and d U T ε s (11) The soution to (10) can be written as x V 1 λ U T ε s (12) 5

6 in which λ are the diagona eements of Λ. If at east one λ i is zero or cose to zero such that the vaue is dominated by the roundoff error, the matrix is singuar. When this happens, 1 λ vaues corresponding to λ faing beow a sma cutoff vaue are a repaced by zeros (Forsythe, et a. 1975, Press, 1989). It can be shown that SVD determines the x for a probem that minimizes x 2 and ε s 2. The soution under these conditions is unique, and incudes the Gaussian, east square and Penrose cases mentioned before. The cutoff eve or the smaest singuar vaue λ min in (12) is usefu in controing the data errors in the soution. If it is sma, 1 λ becomes too arge, and errors in the soution wi be bown out of proportion. In a non singuar matrix, λ min is reated to the condition number, λ max λ min. A matrix is i-conditioned if it is too arge, or the reciproca, λ min λ max approaches foating point precision δ. For singe precision operations, δ 10 6 for most computers and the minimum cutoff possibe is δ λ max. Since data error can cause unnecessary error magnifications at sma cutoffs, a δ vaue of 10 3 was used in the current probem. Even arger vaues (10 2 ) were found to be sufficient during initia iterations. Methods of smoothing the cutoff have been discussed Wiggins (1972) and Uhrhammer (1980). The U and V matrices created by SVD can give additiona information about the parameter behavior. The foowing expression gives an idea about the parameter resoution in the computer mode (Wiggins, 1972): Parameter space resoution, R VV T (13) Reative size of the eements of R show reative resoution or independence of the corrections to the parameters. If the parameters are finey resoved, this matrix is expected to be cose to unity. If the parameters are not individuay resovabe but are resovabe as groups, the probem is characterized by compact resoution, showing square bocks or groups of dominant eements in the matrix. Covariance of the estimated parameters is defined as (Wiis and Yeh, 1987) ϒ cov ˆx E x ˆx x T ˆx (14) in which ϒis the parameter covariance matrix which can be computed as σ 2 A T s A s! determinant of A T s A s 1 when the 1 is not zero. σ 2 error variance at the output. When the matrix is iconditioned, this method fais. However, the foowing expression derived using the resuts of SVD 6

7 gives the covariance matrix for both singuar and non-singuar sensitivity matrices (Wiggins, 1972, Uhrhammer, 1980): Covariance ϒ 2 V σ2 Λ 2 VT (15) A vaue of 1 is used for σ 2 to obtain the reative vaues in the matrix eements. Uhrhammer (1980) referred to the covariance matrix as the uncertainty matrix, and used a damping parameter to damp out the effects of sma singuar vaues on 1 λ during a sharp cutoff. In the current study, a sharp eary cutoff λ min ( 10 3 ) as suggested by Wiggins (1972) is used for simpicity. The covariance matrix can be used to obtain the correation matrix (Uhrhammer, 1980) Correation among parameter corrections, ρ 2 i j ϒ 2 i j ϒ i iϒ j j (16) Correation among parameter corrections is usefu in determining the dependence of parameters and parameter groups. If there are parameters that are dependent as a group, and as a group if they are independent from other groups refecting a bock form in the matrix, caibration can be simpified by casting these groups as singe parameters. The correation matrix, resoution matrix, and V T are used in the current study to understand and isoate the groups of parameters. These groups can aso be reated to the ayout of the river network and gage positions. Gauss-Newton method The Gauss-Newton method as discussed by Wiis and Yeh (1987) considers the identification probem as an optimization probem aimed at minimizing output errors. This method cannot be used with underdetermined systems m n or systems with rank q n. It is used in the study as an aternative method to caibrate groups of parameters seected to satisfy these conditions. When using the method, the objective is to determine parameters x i i 1 2 n such that the summation of the squares of the errors is a minimum, subjected to parameter constraints given by (8). The sum of squares of errors is defined as m S j 1 k 1 w j k y k j Yj k 2 (17) over a the gages and a the time steps. The w j k are weighing factors, which are assumed to be 1 for the current probem. The parameter corrections that minimize S can be shown to be (Yoon and 7

8 Yeh 1976). C x D (18) in which eements i j of C and D are given by c i j n a k p i k 1 p ak p j ; 1 d i n k 1 p 1 a k p i y k p Yp k (19) The n n system of inear equations is soved to obtain the corrections x i, i 1 2 n. The parameters are updated using (7), in which ρ r can be assumed. Optima vaues of ρ r can be determined by the quadratic interpoation method suggested by Yoon and Yeh (1976). Iterations can be terminated when S converges and the parameter corrections are too sma to make physica sense, or when x x machine precision. Minimax method The objective of the minimax method used in the current study is to determine the parameters x i i 1 2 " n that minimize the sum of the absoute vaues of the maximum errors at each of the gages. The maximum errors e 1 e 2 # e m at gages 1 2 m at any time during the entire period are e i max The objective function seected is ε k i k 1 2 i 1 2 m (20) min e 1 e 2 e m (21) If the parameters x i i 1 2 " n are subjected to sma changes x i i 1 2 " n, then simuated vaues y k i i 1 2 " n k 1 2, change to k y$i y$i y k i i 1 2 " n k 1 2. Expressing y k i using (1), y k i can be expressed as y$i k y k i n a i j x j i 1 2 m k 1 2 " (22) j 1 The constraints used to determine maximum error can be set up by confining the new errors after a parameter correction to be e i for a stations. This condition is given by k y$i Y i e i for i 1 2 m and k 1 2 " k (23) 8

9 Substituting for y$i k using (22), (23) can be written as n a k i j x j e i Y k i y k i for i 1 2 " m and k 1 2 (24) j 1 n a k i j j j e i 1 Yi k y k i for i 1 2 m and k 1 2 " (25) Other conditions can aso be incuded to keep the soution within imits. Upper and ower bounds can be set to the parameters using % x i x i x i ub i 1 2 n (26) x i x i x i b i 1 2 n (27) in which, x i ub and x i b are the upper and ower bounds of parameters in (8). The objective function with 2 m 2 n constraints form a inear programming probem, which is soved for the corrections x i, i 1 2 n using the simpex agorithm. Parameters are updated using (7) unti the corrections are too sma to make physica sense or x x machine precision. APPLICATION TO THE UPPER NIAGARA RIVER The mode to be caibrated soves the Saint Venant equations using the four point impicit method. The numerica method is simiar to that used by Potok and Quinn (1979). The Saint Venant equations consist of the foowing continuity and momentum equations: Q t A t Q x Q 2 x A & ga h x 0 (28) p b τ b ρ 0 (29) in which, Q discharge; A fow area; h water surface eevation; x distance aong the river; t time; ρ density of water; p b wetted perimeter; τ b shear stress at channe bottom ρgq 2 n 2 b A 2 R' 1( 3) ; R hydrauic radius A p b ; p b wetted perimeter; n b Manning s coefficient. n b vaues for the discretized river sections are the unknown parameters in the mode to be caibrated. The map of the Upper Niagara River is shown in Fig. 1. In the map, the term NYPA stands for New York Power Authority hydropower intake, and SAB stands for the Sir Adams Beck power 9

10 pant intake on the Canadian side. The river network is discretized into 27 sections. The probem invoves caibration of the (n 27) unknown Mannings coefficients to minimize the errors at a the gaging stations. Houry water eves between Feb 13, 1989, and Feb 16, 1989, are used in the foowing sampe caibration to obtain Manning s coefficients for the range of water eves considered. A compete caibration requires data for a onger period of time, covering a wider range of water eves. The ist of 3 stations with known discharges and 9 stations with known water eves avaiabe to run the mode is shown in Fig. 1. Water eve gages are marked as G1... G9 and known discharges marked as Q10... Q12. The upstream end water eve at Fort Erie gage G1 is specified as the upstream boundary condition, and the downstream end discharges Q10, Q11 and Q12 at the contro dam and intakes are specified as the downstream boundary conditions. Water eve gage G2 is not used because its vaue is very cose (* 0.05 ft) to the vaue of G3, and both are represented by ony one node in the 1-D mode. The remaining 7 gages marked G3-G9 are used in the caibration. In the study, the observations and simuations at these gages are simpy referred to as Y i and y i i 1 m in which m 7. Caibration of individua mode parameters The generaized inear inverse formuation and singuar vaue decomposition are used to caibrate 27 bed roughness vaues n 1 n 27 of the 27 river sections in the 1-D hydrodynamic mode. SVD is needed in this case because there are ony 7 observation stations or 7 equations to sove (5) to determine 27 unknown parameters. A the Manning s coefficients are assumed as in the first iteration. Figure 2 shows the sensitivity matrix A s indicating arge sensitivities of n 2 and n 3 to amost a the water eve gaging stations. In the figure, the circe diameters correspond to the sizes of the eements. n 2 and n 3 correspond to the rapids sections of the upper Niagara river where the drop in water eve is about 2 m and fow veocities reach approximatey 2.6 m/s. Figure 3 shows the convergence of a few seected stabe parameters n 2, n 6, n 11, n 13, and n 16, and the rapid reduction of the error variance in water eves with increasing number of iterations. The rapid convergence indicates that the method is capabe of minimizing the bias error as intended. Figure 3 shows that the method aso reduces the mean square error during caibration. Since the 10

11 probem is underdetermined, the soution shown in the second row of Tabe 1 is nonunique and has a arge random error. Parameter error variance is arge for reativey insensitive parameters as shown ater in Fig. 7. Even if the insensitive parameters have a arge variance, the mode gives a good error reduction as shown in Fig. 4 because the mode output is not affected by them. The set of 27 parameters with arge errors has a imited use. However, the by-products of the method such as the singuar vaues and singuar vectors can be used to understand the parameters of the mode in great detai as expained before. Three of the new ineary independent parameter corrections obtained using z V T x of (11) are shown beow. n n n n n 11 λ (30) n n n n n 10 λ (31) n n n n n 24 λ (32) in which, n+1 n+ 2 are the components of the vector z representing new independent parameter variations. These variations are made up of inear combinations of od parameter variations n 1 n 2. The coefficients in front of n 1 n 2, which form n+i, are obtained from the coumns of V. Variations of the new parameters, n+i are reated to variations in the mode output, h+i, through the reationship h+ i λ i n+i, which is simiar to (10) when the mode output variation is repaced by a negative output error. The singuar vaues λ i of the mode 81210, 2269, 1572, 827, 782, 476 and 22, (nearest integer) show the domination of the caibration by a sma number of parameter groups. The singuar vaues show that the rank of the (7 27) sensitivity matrix is 7 and that ony 7 independent parameter groups can be identified in the probem. The singuar vectors of V written as (30)...(32) show the three most important parameter groups. Singuar vaues λ i provide estimates of the sensitivities. The equations show that adjustment of the group of parameters, consisting of mainy n 3 and n 2 as indicated by the strong coefficients and 0.687, can take care of about 93 % of the overa caibration. Parameters n 2 and n 3 represent the narrow rapids section of the river that has a controing effect on the overa water eves. Parameters n 11 n 16 form another independent group. The inear combinations of variations of the mode output are reated to n+1 n+ 2 through the equation h+ U T h. h h h h h 6 λ (33) 11

12 h h h h h 8 λ (34) h h h h h 9 λ (35) in which, h+1 h+2 are the eements of vector h+ made of inear combinations of variations of water eves h 1 h 2 h 7 which correspond to gages G3-G9. (33)-(35) show that the new parameter n+1 uses simuation errors in water eves h 9 h 7 h 8 h 6 for caibration, and parameter n+2 uses errors in h 4 h 3 for caibration. Approximate equaity of coefficients in (33) shows that a the gages have the same information required to caibrate n+1. Parameter n+ 2 however requires information mosty from h 4 and h 3 (G6 and G5) as shown by coefficients and Figure 5 shows the parameter resoution matrix obtained using (13). In the figure, parameter corrections in groups n 6 n 14 and n 16 n 24 are the east resoved. n 27 is the most resoved, because there are two gages at both ends of the reach. Parameters within a ong river reach without gaging stations in the midde are not sufficienty resoved, and behave as a group. Any further resoution of parameters within these zones is not possibe because observation points cannot distinguish between roughnesses of any two sections within a group if the group ies between two gages. It is aso unnecessary to resove the parameters any further uness the mode is used to estimate the water eves at the interior points. The accuracy of the parameter corrections is of the order 1 λ min as shown by (12). If a sma cutoff vaue λ min is seected, the parameters wi have arger errors. If a arge cutoff is seected, ony a few parameter groups wi be seected as expained before. With a very arge cutoff, the number of parameters may not be adequate to represent the physica system, and sufficienty minimize the errors. The next sections incudes determination of an optima parameter dimension based on the parameter uncertainty and the east square output error. The correation matrix in Fig. 6 shows the group behavior of the parameters. The figure shows the strong correation among parameter groups n 6 n 8, n 9 n 12, n 13 n 15, n 16 n 18, n 19 n 24, n 25 n 26. The map in Fig. 1 shows that these groups correspond to river reaches with water eve gages at both ends. Next sections describe the caibration using parameter groups. Caibration of parameter groups Caibration of the individua parameters in the previous section shows that they tend to group together as indicated by the correation, resoution and V matrices. In the SVD method expained 12

13 earier, singuar vectors in V coud impicity make use of these groups. However, 27+1 computer runs were needed for one iteration of the caibration. In this section, in order to reduce the number of runs and obtain more stabe parameters with ow error variances, 7 new groups of parameters are created from the origina parameters using the information avaiabe from the previous section. A singe parameter vaue is assigned to a the parameters in a group. The parameter groups created in this manner are shown in Tabe 2 and are marked as m 1 m 7. The number of groups is imited to 7 because it is the rank of A s. To study the effect of the number of parameter groups, or parameter dimension on the output error and parameter uncertainty, the 7 parameter groups are combined again with each other to form smaer numbers of parameter groups. Coumn 4 of Tabe 2 for exampe shows how parameters n 1 n 27 are combined to form 4 groups. Figure 7 shows the variation of the output error and parameter uncertainty with different numbers of parameter groups. The mean square bias error is used to measure the output error in which bias is the average difference between the observed and simuated water eves. An approximate estimate of the parameter error or uncertainty is computed using the expression σ 2 λ 2 min, which is approximatey equa to the infinite norm of the covariance matrix. Figure 7 shows that the parameter uncertainty increases and the output error decreases with increasing number of parameter groups. The figure aso shows that the output error variance does not reduce much when there are more than 3 parameters. With more than 3 parameters, the parameter uncertainty increases. The figure shows the existence of an optima parameter dimension of about 3 for the probem where both the output error and parameter uncertainty are ow. SVD with 7 parameter groups shown in Tabe 2 gave a mean square bias error of after 7 iterations showing that the caibration is successfu even if these handpicked parameter groups are not exacty the singuar vectors in (30)-(32). The condition number of the 7 7 group sensitivity matrix is /23.8, which shows that the matrix is nonsinguar, and it is possibe to use the Gauss-Newton method for comparison. The parameter resoution and correation matrices of the 7 parameter group caibration (not shown) indicate that the handpicked parameter groups are sufficienty resoved and independent. The third and fourth rows of Tabe 1 show the vaues of parameters obtained by caibrating the handpicked 7 and 3 groups of parameters using SVD. A these are even determined probems that can aso be soved using Gaussian eimination. The vaue of a singe parameter that fits the 13 4 m 2

14 entire river network is found to be , with an output mean square bias error of m 2. Caibration using Gauss-Newton and Minimax methods The Gauss-Newton and minimax methods are aso used to caibrate the same 7 handpicked parameter groups for comparison. The fina Mannings coefficients after 5 iterations of caibration using both methods are shown in Tabe 3 aong with the sum of error variances and mean square bias errors. The argest correction required on the Mannings coefficients after 5 iterations of the Gauss- Newton method is approximatey , and the argest correction required after 8 iterations of the minimax method is approximatey The iterations were terminated at this point due to the smaness of the adjustment. During the caibration using the Gauss-Newton method, ρ was assumed as 0.75 during the first two iterations, and 1.0 afterwards. Larger vaues of ρ were avoided during first few iterations to avoid crashing the mode. The Gauss-Newton method required the east amount of computer time and the fastest rate of convergence. The minimax method is the most time consuming because of the inear programming (LP) agorithm. The computer time requirement increased with the number of time steps in the simuation runs because it increased the size of the inear programming probem. A Sparc 1 computer using the IMSL package needed minutes of computer time for a singe iteration of the minimax, when compared with spit seconds needed by other methods. Tabe 3 shows that parameter vaues are approximatey the same for the most important parameter group m 1 consisting of n 1 n 6. Pots of the observed and simuated water eves using a three methods do not show any visuay detectabe difference. CONCLUSIONS SVD is usefu in identifying and understanding an arbitrary number of parameters when the equations invoved are overdetermined, even-determined, underdetermined or singuar. Caibration of the 27 parameters in the Niagara river mode using different numbers of parameter groups show that the mode has a maximum of 7 independent parameter groups, and the optimum parameter dimension is about 3. 14

15 The study shows that the underdetermined identification probem with 27 parameters can be reduced to an even-determined or an over-determined probem by creating a number of parameter groups equa to the rank of the sensitivity matrix (7) or ess. The groups can be created using the information from parameter space resoution matrices, correation matrices, singuar vectors or the geometrica ayout of the river system. Parameters corresponding to river reaches between gages appear as groups in the correation and resoution matrices, and therefore can be considered as groups. Once parameter groups are created such that the probem is neither underdetermined nor singuar, optimization methods such as the Gauss Newton and minimax methods can be used to carry out the same caibration. 15

16 ACKNOWLEDGEMENTS The study was partiay supported by the New York Power Authority through contract No entited Upper Niagara River Ice Dynamics Simuation. Mr. Randy Crissman is the Contracting Officer s representative. The author wishes to thank Prof. H. T. Shen for having him in his research team. 16

17 APPENDIX I. REFERENCES Beck, M. B. (1987). Water Quaity Modeing: A Review of the Anaysis of Uncertainty, Water Resources Research, 23(8), Becker, L., Yeh, W. W-G. (1972). Identification of parameters in unsteady open channe fows, Water Resources Research, 8(4), Becker, L., Yeh, W. W-G. (1973). Identification of mutipe reach channe parameters, Water Resources Research, 9(2), Chiu, Chao-Lin, (1978). Kaman Fiter in Open Channe Fow Estimation, J. of Hydr. Engg., ASCE, 104(8). Forsythe, E. George, Macom, A. Michae, and Moer, B. Ceve, (1977), Computer Methods for Mathematica Computations, Prentice Ha Inc. Forsythe, E. G., Macom A. M., and Moer, B. C. (1977). Computer methods for mathematica computations, Prentice Ha, Inc. Fread, D. L., and Smith, G. F. (1977). Caibration techniques for 1-D unsteady fow modes, J. of the Hydr. Div., ASCE, 104, Juy, Katopodes, N. D., Jyh-Haw Tang, and Cemments, Abert J. (1990). Estimation of Surface Irrigation parameters, J. Hydr. Engg., ASCE, 116(5), Meneke, Wiiam, (1984). Geophysica data anaysis: discrete inverse theory, Academic Press, Inc., New York. Nobe, B., and Danie, J. (1975). Appied Linear Agebra, Prentice Ha Pubishers, Engewood Ciff, NJ. Potok, A. J., and Quinn, F. H. (1979), A Hydrodynamic Mode for the Upper St. Lawrence River 17

18 for Water Resources Studies, Water Res. Bu., 15(6). Press, W. H., Fannery, B. P., Teukosky, S. A., and Verrering, W. T. (1989). Numerica Recipes, the Art of Scientific Computing, Cambridge University Press, Cambridge. Uhrhammer, R. A. (1980). Anaysis of Seismographic Station Networks, Bu. of the Seis. Soc. of America, 70(4), Ward, S. H., Peepes, W. J., and Ryu, J. (1973). Anaysis of Geoeectromagnetic Data, Methods in Computationa Physics, 13, Academic Press, New York, 13, Wiggins, R. A. (1972). The Genera Linear Inverse Probem: impications of surface waves and free osciations for surface structures, Rev. Geophysics and Space Physics, 10, Wiis, Robert and Yeh, W-G. Wiiam (1987). Groundwater Systems Panning and Management, Prentice Ha, Inc. Yeh, W. W-G., and Becker, L. (1973). Linear programming and channe fow identification, Journa of the Hydrauics Division, ASCE, 99(11), Yoon, Y. S., and Yeh, W. W-G. (1976). Parameter Identification in an Inhomogeneous Medium with the Finite-Eement method, Soc. Pet. Eng. J.,

19 APPENDIX II. NOTATION The foowing symbos are used in the paper A k A sensitivity matrix a k j i cross sectiona area of a channe A s k 1 Ak a k ji sensitivity of j th water eve with respect to i th parameter at time k. d intermediate variabe defined as U T ε s which give inear combinations of differences in water eves h i k m m j n n b n j n+j Q q r R R water eve at gage i time step number number of time steps number of observation stations group vaues of handpicked parameters number of parameters Manning s coefficient Manning s coefficient of river section j new parameter combination obtained using the SVD method discharge in a channe rank of the matrix, or number of singuar vaues used number of iterations hydrauic radius parameter resoution matrix U V m m and n n matrices of orthogona singuar vectors x i Yj k y k j x ε s ε s j a generic symbo for parameter i observed water eve at node j at time k (m) simuated water eve at node j at time k (m) vector of parameter corrections T n 1 n 2 vector consisting of summations of errors summation of errors at gage j over time 19

20 Λ diagona matrix of singuar vaues λ i i th singuar vaue ρ i j correation among parameter corrections τ b shear stress at channe bottom ϒ 2 covariance of parameter corrections Subscripts s summation over time r iteration number Superscripts k vaue at time step k 20

21 Figures Fig 1: Map of the upper Niagara River. Fig 2: Sensitivity matrix (Max diameter=31519); gages G3-G9 are marked as 1-7. Fig 3: Convergence of parameters and reduction of errors with iterations. Fig 4: Observed and simuated water eves after caibration. Fig 5: Parameter resoution matrix (Max. Diameter = 0.96). Fig 6: Correation matrix (Max diameter = 1.0). Fig 7: Variation of output error and parameter uncertainty with parameter dimension. 21

22 Tabe 1: Manning s Roughness Coefficients obtained using SVD Param Indiv groups groups Param Indiv groups groups Param Indiv groups groups KEY WORDS Bed roughness caibration, cana network, Open chane fow Hydrauics, Inverse probem, parameter identification, Singuar vaue decomposition, east square, Gauss-Newton, Linear Programming, 22

23 Tabe 2: Parameter groups formed by combining individua reach parameters; a groups are combinations of the 7 basic groups in coumn 2 Groups of Number of par. groups Origina and formation Parameters n 1 n 5 m 1 m 1 m 1 m 1 m 1 m 1 n 6 n 8 m 2 m 2 m 2 m 2 m 2 m 1 n 9 n 12 m 3 m 3 m 3 m 3 m 2 m 1 n 13 n 15 m 4 m 4 m 3 m 3 m 2 m 1 n 16 n 18 m 5 m 5 m 4 m 2 m 2 m 1 n 19 n 26 m 6 m 5 m 1 m 3 m 2 m 1 n 27 m 7 m 5 m 3 m 3 m 2 m 1 Tabe 3: Resuts of the caibration of 7 handpicked parameter groups using different methods Group Min. Bias Least Sqr. Minimax Err. Variance (m 2 ) Mean Sq. Bias (m

24 Parameter Individua 7 groups 3 groups

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

Haar Decomposition and Reconstruction Algorithms

Haar Decomposition and Reconstruction Algorithms Jim Lambers MAT 773 Fa Semester 018-19 Lecture 15 and 16 Notes These notes correspond to Sections 4.3 and 4.4 in the text. Haar Decomposition and Reconstruction Agorithms Decomposition Suppose we approximate

More information

Stochastic Variational Inference with Gradient Linearization

Stochastic Variational Inference with Gradient Linearization Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,

More information

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain CORRECTIONS TO CLASSICAL PROCEDURES FOR ESTIMATING THURSTONE S CASE V MODEL FOR RANKING DATA Aberto Maydeu Oivares Instituto de Empresa Marketing Dept. C/Maria de Moina -5 28006 Madrid Spain Aberto.Maydeu@ie.edu

More information

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA) 1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using

More information

Copyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU

Copyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water

More information

ASummaryofGaussianProcesses Coryn A.L. Bailer-Jones

ASummaryofGaussianProcesses Coryn A.L. Bailer-Jones ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Theory and implementation behind: Universal surface creation - smallest unitcell

Theory and implementation behind: Universal surface creation - smallest unitcell Teory and impementation beind: Universa surface creation - smaest unitce Bjare Brin Buus, Jaob Howat & Tomas Bigaard September 15, 218 1 Construction of surface sabs Te aim for tis part of te project is

More information

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation

More information

IE 361 Exam 1. b) Give *&% confidence limits for the bias of this viscometer. (No need to simplify.)

IE 361 Exam 1. b) Give *&% confidence limits for the bias of this viscometer. (No need to simplify.) October 9, 00 IE 6 Exam Prof. Vardeman. The viscosity of paint is measured with a "viscometer" in units of "Krebs." First, a standard iquid of "known" viscosity *# Krebs is tested with a company viscometer

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

BALANCING REGULAR MATRIX PENCILS

BALANCING REGULAR MATRIX PENCILS BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,

More information

Tracking Control of Multiple Mobile Robots

Tracking Control of Multiple Mobile Robots Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart

More information

International Journal of Mass Spectrometry

International Journal of Mass Spectrometry Internationa Journa of Mass Spectrometry 280 (2009) 179 183 Contents ists avaiabe at ScienceDirect Internationa Journa of Mass Spectrometry journa homepage: www.esevier.com/ocate/ijms Stark mixing by ion-rydberg

More information

Primal and dual active-set methods for convex quadratic programming

Primal and dual active-set methods for convex quadratic programming Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:

More information

Sequential Decoding of Polar Codes with Arbitrary Binary Kernel

Sequential Decoding of Polar Codes with Arbitrary Binary Kernel Sequentia Decoding of Poar Codes with Arbitrary Binary Kerne Vera Miosavskaya, Peter Trifonov Saint-Petersburg State Poytechnic University Emai: veram,petert}@dcn.icc.spbstu.ru Abstract The probem of efficient

More information

arxiv: v4 [math.na] 25 Aug 2014

arxiv: v4 [math.na] 25 Aug 2014 BIT manuscript No. (wi be inserted by the editor) A muti-eve spectra deferred correction method Robert Speck Danie Ruprecht Matthew Emmett Michae Minion Matthias Boten Rof Krause arxiv:1307.1312v4 [math.na]

More information

Chemical Kinetics Part 2

Chemical Kinetics Part 2 Integrated Rate Laws Chemica Kinetics Part 2 The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates the rate

More information

Statistics for Applications. Chapter 7: Regression 1/43

Statistics for Applications. Chapter 7: Regression 1/43 Statistics for Appications Chapter 7: Regression 1/43 Heuristics of the inear regression (1) Consider a coud of i.i.d. random points (X i,y i ),i =1,...,n : 2/43 Heuristics of the inear regression (2)

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

Introduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR 1 Maximizing Sum Rate and Minimizing MSE on Mutiuser Downink: Optimaity, Fast Agorithms and Equivaence via Max-min SIR Chee Wei Tan 1,2, Mung Chiang 2 and R. Srikant 3 1 Caifornia Institute of Technoogy,

More information

Research of Data Fusion Method of Multi-Sensor Based on Correlation Coefficient of Confidence Distance

Research of Data Fusion Method of Multi-Sensor Based on Correlation Coefficient of Confidence Distance Send Orders for Reprints to reprints@benthamscience.ae 340 The Open Cybernetics & Systemics Journa, 015, 9, 340-344 Open Access Research of Data Fusion Method of Muti-Sensor Based on Correation Coefficient

More information

A. Distribution of the test statistic

A. Distribution of the test statistic A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch

More information

Related Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage

Related Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is

More information

The EM Algorithm applied to determining new limit points of Mahler measures

The EM Algorithm applied to determining new limit points of Mahler measures Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,

More information

arxiv: v1 [math.ca] 6 Mar 2017

arxiv: v1 [math.ca] 6 Mar 2017 Indefinite Integras of Spherica Besse Functions MIT-CTP/487 arxiv:703.0648v [math.ca] 6 Mar 07 Joyon K. Boomfied,, Stephen H. P. Face,, and Zander Moss, Center for Theoretica Physics, Laboratory for Nucear

More information

An explicit Jordan Decomposition of Companion matrices

An explicit Jordan Decomposition of Companion matrices An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057

More information

Source and Relay Matrices Optimization for Multiuser Multi-Hop MIMO Relay Systems

Source and Relay Matrices Optimization for Multiuser Multi-Hop MIMO Relay Systems Source and Reay Matrices Optimization for Mutiuser Muti-Hop MIMO Reay Systems Yue Rong Department of Eectrica and Computer Engineering, Curtin University, Bentey, WA 6102, Austraia Abstract In this paper,

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

Nonlinear Analysis of Spatial Trusses

Nonlinear Analysis of Spatial Trusses Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes

More information

(This is a sample cover image for this issue. The actual cover is not yet available at this time.)

(This is a sample cover image for this issue. The actual cover is not yet available at this time.) (This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna

More information

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems Componentwise Determination of the Interva Hu Soution for Linear Interva Parameter Systems L. V. Koev Dept. of Theoretica Eectrotechnics, Facuty of Automatics, Technica University of Sofia, 1000 Sofia,

More information

Fitting Algorithms for MMPP ATM Traffic Models

Fitting Algorithms for MMPP ATM Traffic Models Fitting Agorithms for PP AT Traffic odes A. Nogueira, P. Savador, R. Vaadas University of Aveiro / Institute of Teecommunications, 38-93 Aveiro, Portuga; e-mai: (nogueira, savador, rv)@av.it.pt ABSTRACT

More information

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p

More information

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport

More information

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing

More information

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,

More information

Traffic data collection

Traffic data collection Chapter 32 Traffic data coection 32.1 Overview Unike many other discipines of the engineering, the situations that are interesting to a traffic engineer cannot be reproduced in a aboratory. Even if road

More information

Radar/ESM Tracking of Constant Velocity Target : Comparison of Batch (MLE) and EKF Performance

Radar/ESM Tracking of Constant Velocity Target : Comparison of Batch (MLE) and EKF Performance adar/ racing of Constant Veocity arget : Comparison of Batch (LE) and EKF Performance I. Leibowicz homson-csf Deteis/IISA La cef de Saint-Pierre 1 Bd Jean ouin 7885 Eancourt Cede France Isabee.Leibowicz

More information

Computational studies of discrete breathers. Sergej Flach MPIPKS Dresden January 2003

Computational studies of discrete breathers. Sergej Flach MPIPKS Dresden January 2003 Computationa studies of discrete breathers Sergej Fach MPIPKS Dresden January 2003 CONTENT: 0. A bit on numerics of soving ODEs 1. How to observe breathers in simpe numerica runs 2. Obtaining breathers

More information

BP neural network-based sports performance prediction model applied research

BP neural network-based sports performance prediction model applied research Avaiabe onine www.jocpr.com Journa of Chemica and Pharmaceutica Research, 204, 6(7:93-936 Research Artice ISSN : 0975-7384 CODEN(USA : JCPRC5 BP neura networ-based sports performance prediction mode appied

More information

High Spectral Resolution Infrared Radiance Modeling Using Optimal Spectral Sampling (OSS) Method

High Spectral Resolution Infrared Radiance Modeling Using Optimal Spectral Sampling (OSS) Method High Spectra Resoution Infrared Radiance Modeing Using Optima Spectra Samping (OSS) Method J.-L. Moncet and G. Uymin Background Optima Spectra Samping (OSS) method is a fast and accurate monochromatic

More information

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017 In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative

More information

Statistical Learning Theory: A Primer

Statistical Learning Theory: A Primer Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO

More information

SydU STAT3014 (2015) Second semester Dr. J. Chan 18

SydU STAT3014 (2015) Second semester Dr. J. Chan 18 STAT3014/3914 Appied Stat.-Samping C-Stratified rand. sampe Stratified Random Samping.1 Introduction Description The popuation of size N is divided into mutuay excusive and exhaustive subpopuations caed

More information

Explicit overall risk minimization transductive bound

Explicit overall risk minimization transductive bound 1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,

More information

Chemical Kinetics Part 2. Chapter 16

Chemical Kinetics Part 2. Chapter 16 Chemica Kinetics Part 2 Chapter 16 Integrated Rate Laws The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates

More information

A Novel Learning Method for Elman Neural Network Using Local Search

A Novel Learning Method for Elman Neural Network Using Local Search Neura Information Processing Letters and Reviews Vo. 11, No. 8, August 2007 LETTER A Nove Learning Method for Eman Neura Networ Using Loca Search Facuty of Engineering, Toyama University, Gofuu 3190 Toyama

More information

17 Lecture 17: Recombination and Dark Matter Production

17 Lecture 17: Recombination and Dark Matter Production PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was

More information

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University Turbo Codes Coding and Communication Laboratory Dept. of Eectrica Engineering, Nationa Chung Hsing University Turbo codes 1 Chapter 12: Turbo Codes 1. Introduction 2. Turbo code encoder 3. Design of intereaver

More information

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SIAM J. NUMER. ANAL. Vo. 0, No. 0, pp. 000 000 c 200X Society for Industria and Appied Mathematics VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW

More information

https://doi.org/ /epjconf/

https://doi.org/ /epjconf/ HOW TO APPLY THE OPTIMAL ESTIMATION METHOD TO YOUR LIDAR MEASUREMENTS FOR IMPROVED RETRIEVALS OF TEMPERATURE AND COMPOSITION R. J. Sica 1,2,*, A. Haefee 2,1, A. Jaai 1, S. Gamage 1 and G. Farhani 1 1 Department

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

Cryptanalysis of PKP: A New Approach

Cryptanalysis of PKP: A New Approach Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in

More information

Methods for Ordinary Differential Equations. Jacob White

Methods for Ordinary Differential Equations. Jacob White Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa

More information

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled. imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,

More information

APPENDIX C FLEXING OF LENGTH BARS

APPENDIX C FLEXING OF LENGTH BARS Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or

More information

Published in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics

Published in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of

More information

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW Abstract. One of the most efficient methods for determining the equiibria of a continuous parameterized

More information

СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS

СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА Милко Стоянов Милошев 1, Константин Савков Казаков 2 Висше Строително Училище Л. Каравелов - София COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS

More information

Mode in Output Participation Factors for Linear Systems

Mode in Output Participation Factors for Linear Systems 2010 American ontro onference Marriott Waterfront, Batimore, MD, USA June 30-Juy 02, 2010 WeB05.5 Mode in Output Participation Factors for Linear Systems Li Sheng, yad H. Abed, Munther A. Hassouneh, Huizhong

More information

NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION

NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION Hsiao-Chang Chen Dept. of Systems Engineering University of Pennsyvania Phiadephia, PA 904-635, U.S.A. Chun-Hung Chen

More information

Absolute Value Preconditioning for Symmetric Indefinite Linear Systems

Absolute Value Preconditioning for Symmetric Indefinite Linear Systems MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Absoute Vaue Preconditioning for Symmetric Indefinite Linear Systems Vecharynski, E.; Knyazev, A.V. TR2013-016 March 2013 Abstract We introduce

More information

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, as we as various

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

Solving Maxwell s Equations Using the Ultra Weak Variational Formulation

Solving Maxwell s Equations Using the Ultra Weak Variational Formulation Soving Maxwe s Equations Using the Utra Weak Variationa Formuation T. Huttunen, M. Mainen and P. Monk Department of Appied Physics, University of Kuopio, P.O.Box 1627, 7211 Kuopio, Finand Department of

More information

Problem Set 6: Solutions

Problem Set 6: Solutions University of Aabama Department of Physics and Astronomy PH 102 / LeCair Summer II 2010 Probem Set 6: Soutions 1. A conducting rectanguar oop of mass M, resistance R, and dimensions w by fas from rest

More information

Laboratory Exercise 1: Pendulum Acceleration Measurement and Prediction Laboratory Handout AME 20213: Fundamentals of Measurements and Data Analysis

Laboratory Exercise 1: Pendulum Acceleration Measurement and Prediction Laboratory Handout AME 20213: Fundamentals of Measurements and Data Analysis Laboratory Exercise 1: Penduum Acceeration Measurement and Prediction Laboratory Handout AME 20213: Fundamentas of Measurements and Data Anaysis Prepared by: Danie Van Ness Date exercises to be performed:

More information

Paper presented at the Workshop on Space Charge Physics in High Intensity Hadron Rings, sponsored by Brookhaven National Laboratory, May 4-7,1998

Paper presented at the Workshop on Space Charge Physics in High Intensity Hadron Rings, sponsored by Brookhaven National Laboratory, May 4-7,1998 Paper presented at the Workshop on Space Charge Physics in High ntensity Hadron Rings, sponsored by Brookhaven Nationa Laboratory, May 4-7,998 Noninear Sef Consistent High Resoution Beam Hao Agorithm in

More information

SASIMI: Sparsity-Aware Simulation of Interconnect-Dominated Circuits with Non-Linear Devices

SASIMI: Sparsity-Aware Simulation of Interconnect-Dominated Circuits with Non-Linear Devices SASIMI: Sparsity-Aware Simuation of Interconnect-Dominated Circuits with Non-Linear Devices Jitesh Jain, Stephen Cauey, Cheng-Kok Koh, and Venkataramanan Baakrishnan Schoo of Eectrica and Computer Engineering

More information

Inductive Bias: How to generalize on novel data. CS Inductive Bias 1

Inductive Bias: How to generalize on novel data. CS Inductive Bias 1 Inductive Bias: How to generaize on nove data CS 478 - Inductive Bias 1 Overfitting Noise vs. Exceptions CS 478 - Inductive Bias 2 Non-Linear Tasks Linear Regression wi not generaize we to the task beow

More information

Kernel pea and De-Noising in Feature Spaces

Kernel pea and De-Noising in Feature Spaces Kerne pea and De-Noising in Feature Spaces Sebastian Mika, Bernhard Schokopf, Aex Smoa Kaus-Robert Muer, Matthias Schoz, Gunnar Riitsch GMD FIRST, Rudower Chaussee 5, 12489 Berin, Germany {mika, bs, smoa,

More information

Partial permutation decoding for MacDonald codes

Partial permutation decoding for MacDonald codes Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics

More information

Scalable Spectrum Allocation for Large Networks Based on Sparse Optimization

Scalable Spectrum Allocation for Large Networks Based on Sparse Optimization Scaabe Spectrum ocation for Large Networks ased on Sparse Optimization innan Zhuang Modem R&D Lab Samsung Semiconductor, Inc. San Diego, C Dongning Guo, Ermin Wei, and Michae L. Honig Department of Eectrica

More information

On a geometrical approach in contact mechanics

On a geometrical approach in contact mechanics Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128

More information

Crystallisation of a supercooled spherical nodule in a flow

Crystallisation of a supercooled spherical nodule in a flow EUROTHERM 69 Heat and Mass Transfer in Soid-Liquid Phase hange Processes June 25-27, 2003, Bistra caste, Ljubjana, Sovenia Eds.: B. Sarer, D. Gobin rystaisation of a supercooed spherica nodue in a fow

More information

Multilevel Monte Carlo Methods and Applications to Elliptic PDEs with Random Coefficients

Multilevel Monte Carlo Methods and Applications to Elliptic PDEs with Random Coefficients CVS manuscript No. wi be inserted by the editor) Mutieve Monte Caro Methods and Appications to Eiptic PDEs with Random Coefficients K.A. Ciffe M.B. Gies R. Scheich A.L. Teckentrup Received: date / Accepted:

More information

Finite element method for structural dynamic and stability analyses

Finite element method for structural dynamic and stability analyses Finite eement method for structura dynamic and stabiity anayses Modue-9 Structura stabiity anaysis Lecture-33 Dynamic anaysis of stabiity and anaysis of time varying systems Prof C S Manohar Department

More information

Numerical simulation of javelin best throwing angle based on biomechanical model

Numerical simulation of javelin best throwing angle based on biomechanical model ISSN : 0974-7435 Voume 8 Issue 8 Numerica simuation of javein best throwing ange based on biomechanica mode Xia Zeng*, Xiongwei Zuo Department of Physica Education, Changsha Medica University, Changsha

More information

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness 1 Scheduabiity Anaysis of Deferrabe Scheduing Agorithms for Maintaining Rea-Time Data Freshness Song Han, Deji Chen, Ming Xiong, Kam-yiu Lam, Aoysius K. Mok, Krithi Ramamritham UT Austin, Emerson Process

More information

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................

More information

Two-Stage Least Squares as Minimum Distance

Two-Stage Least Squares as Minimum Distance Two-Stage Least Squares as Minimum Distance Frank Windmeijer Discussion Paper 17 / 683 7 June 2017 Department of Economics University of Bristo Priory Road Compex Bristo BS8 1TU United Kingdom Two-Stage

More information

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity

More information

Florida State University Libraries

Florida State University Libraries Forida State University Libraries Eectronic Theses, Treatises and Dissertations The Graduate Schoo Construction of Efficient Fractiona Factoria Mixed-Leve Designs Yong Guo Foow this and additiona works

More information

Available online at ScienceDirect. Procedia Computer Science 96 (2016 )

Available online at  ScienceDirect. Procedia Computer Science 96 (2016 ) Avaiabe onine at www.sciencedirect.com ScienceDirect Procedia Computer Science 96 (206 92 99 20th Internationa Conference on Knowedge Based and Inteigent Information and Engineering Systems Connected categorica

More information