Hyperbolic Method for Dispersive PDEs: Same High- Order of Accuracy for Solution, Gradient, and Hessian
|
|
- Jade Burke
- 5 years ago
- Views:
Transcription
1 AIAA Aviation 3-7 Jne 6, Washington, D.C. 46th AIAA Flid Dynamics Conference AIAA Hyperbolic Method for Dispersive PDEs: Same High- Order of Accracy for, Gradient, and Hessian Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ Alireza Mazaheri NASA Langley Research Center, Hampton, VA 368 Mario Ricchito INRIA Bordea Sd-Oest and Institt de Mathématiqes de Bordea, 3345 Talence Cede, France Hiroaki Nishikawa National Institte of Aerospace, Hampton, VA 3666 In this paper, we introdce a new hyperbolic first-order system for general dispersive partial differential eqations (PDEs). We then etend the proposed system to general advection-diffsion-dispersion PDEs. We apply the forth-order RD scheme of Ref. to the proposed hyperbolic system, and solve time-dependent dispersive eqations, inclding the classical two-soliton KdV and a dispersive shock case. We demonstrate that the predicted reslts, inclding the gradient and Hessian (second derivative), are in a very good agreement with the eact soltions. We then show that the RD scheme applied to the proposed system accrately captres dispersive shocks withot nmerical oscillations. We also verify that the soltion, gradient and Hessian are predicted with eqal order of accracy. I. Introdction Dispersion effects play a fndamental role in many applications involving hydrodynamics. At the large scale, the flow is dominated by advection, while the dissipative effects are more important at the microscopic level. At the mesoscopic level (the intermediate level) the dispersive effects become important, as it is the case, for eample, in nonlinear optics, electromagnetism, qantm mechanics, 4 relativity and Bose-Einstein condensates, 5, 6 atmospheric, coastal, and flvial hydrodynamics, 7 magma, highly viscos flids and/or capillary effects. 4 The main physical effects associated with dispersion are the appearance of dispersive (or ndlating) shocks, and the eistence of smooth traveling solitary waves, which may prodce comple interactions with one another. These systems, independent of the physical natre of the involved medim, admit a mesoscale hydrodynamic model, which consists of a set of Partial Differential Eqations (PDEs). These PDEs are mostly reglarizations of hyperbolic models that, in one-dimension, may be written as t + F () = D, () where, in a classical sense (endowed with an entropy pair, with a diagonalizable fl Jacobian F (), etc.), the left hand side defines a hyperbolic model. Depending on the application and the physical hypotheses made, the reglarization on the right hand side may take different forms, and have a dissipative and/or dispersive character (see e.g., Refs. 3, 5 9). Therefore, term D, may be a combination of the following forms: (a) dissipative reglarization D = ν, (b) dispersive reglarization with time derivative D = ɛ t, Research Aerospace Engineer, Aerothermodynamics Branch, M/S 48A, alireza.mazaheri@nasa.gov Grop Leader at INRIA Sd-Oest, mario.ricchito@inria.fr Associate Research Fellow, hiro@nianet.org of 3 This material is declared a work of the U.S. Government and is not sbject to copyright protection in the United States.
2 (c) flly dispersive reglarization D = ɛ. The form (a) is a classical viscos reglarization, while form (b) can be recast as a first-order evoltionary PDE with an embedded steady state second-order time independent elliptic problem; i.e., t w = f, + = w. Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ We previosly showed e.g., in Refs.,, how one can constrct very accrate nmerical approimations of problems of form (a) and the elliptic form of (b), by first reformlating the system of PDEs as a first-order hyperbolic system. The schemes proposed with the hyperbolic formlation of PDE systems are pwind and highly accrate for both the soltion and its gradient, and have a natral potential for etension on arbitrary nstrctred meshes as illstrated in Refs., 3. The presence of a third-order derivative term in form (c), however, introdces discretization difficlties, which are often related to the nderstanding of the type of stencil that is reqired to approimate these high-order derivative terms, the stability of the method sed, and the imposition of bondary conditions. In Refs. 4 and 7, possible soltions to some of these isses are proposed, where the athors showed very good reslts with a non-hyperbolic first-order system reformlation of a PDE and carefl discretization of the fles. The hyperbolic reformlation of dispersive PDEs similar to the one initially proposed for diffsion in Ref. alleviate the above mentioned isses. However, it is shown and proved in Ref. 5 that the hyperbolic formlation of a dispersive PDE in the form given in Ref. is not possible. Ths, we are motivated to introdce an alternative hyperbolic formlation that is careflly designed for general dispersive PDEs that are relevant to, for eample, qantm mechanics, relativistic hydrodynamics, and coastal engineering 5, 8, 4 applications, sch as the Korteweg-de Vries (KdV) eqation. In this work, we first present a first-order system for a pre dispersion eqation, and show how the proposed formlation can be made hyperbolic. As an intermediate generalization, we then show that an advective-dispersive PDE (sch as the classical KdV) can also be made flly-hyperbolic as well. The flly-hyperbolic advection-dispersion system cold be sefl in imposing a characteristics bondary condition. A practical etension of the proposed hyperbolic dispersion system for general advective-diffsivedispersive PDEs follows net. We then present some nmerical eamples by applying the high-order residaldistribtion (RD) scheme of Ref. to the proposed system, and solving general dispersive PDEs, inclding the classical KdV eqation, on randomly distribted nodes. We verify the order of accracy of the scheme for both the soltion, the gradient, and the Hessian (second derivative) with the se of method of manfactred soltion, and show that the RD scheme applied to the proposed hyperbolic system prodces accrate soltion, gradient and Hessian with eqal order of accracy. The ability to obtain accrate gradient and Hessian are very important as they are sed in many hydrodynamic dispersive models to define physically relevant qantities (e.g., potentials, energy, flow properties at an arbitrary depth, etc.). We also present soltions for the zero dispersion limit of conservation laws and demonstrate that the RD scheme applied to the proposed hyperbolic system is robst and can captre physical oscillations associated with the generated dispersive shocks. II. Hyperbolic Dispersion In this section, we start with a time-dependent dispersive PDE, and reformlate it to a first-order system that can be sccessflly transformed to a first-order hyperbolic system. Consider the following linear dispersive PDE that is often referred to as the Airy eqation, t = ɛ, () where ɛ is the dispersion coefficient (positive or negative). Following the process we otlined in Ref. (althogh other choices may also be possible), we consider here the semi-discrete form of Eq. () obtained with some implicit time integration scheme: α t = ɛ + s(),, where α and s() depend on the time discretization and the known vales of. We then replace the of 3
3 semi-discrete PDE by the steady limit of the psedo-time dependent system τ = ɛ q α t + s(), τ p = T ɛ ( p γ), (3) τ q = T ɛ (γ + p q), Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ where τ is the psedo time, t is the physical time, and T ɛ and γ are, respectively, the dispersion relaation time and an arbitrary constant, both to be defined later. It is easy to verify that in the above system satisfies the original semi-discrete dispersion eqation at psedo steady state (see also Refs. and ). Note that, the proposed system with γ = redces to the first-order system formlation that is proven in Ref. 5 not to be hyperbolic. The stdy of the characteristic polynomial of the Jacobian matri of the qasi-linear differential operator on the right hand side of Eq. (3) reveals that the form of γ is critical in constrcting a hyperbolic system. For eample, we have fond that the following form of γ cold reslt in a hyperbolic formlation of the first-order system for dispersion γ = β + ɛ β T ɛ, (4) where for dimensional consistency β = κ/l ɛ, κ is an arbitrary constant, and L ɛ is the dispersion length scale, to be defined later. With this choice of γ, the first-order system (3) admits the characteristic speeds = βt ɛ, = λdisp ( + + 4ɛβ 3 T ɛ ), 3 = λdisp with the following corresponding eigenvectors ɛ T ɛ 3 T ɛ R = ɛβ β /(ɛβ) β 3 /(ɛβ) ( + 4ɛβ 3 T ɛ ), (5), (6) which is non-singlar (nless κ =,, or ). We now define the relaation time T ɛ as the ratio between the dispersion length scale, L ɛ, and either of the λ or the λ 3, and arrive at the following relation: T ɛ = L3 ɛ ɛ, (7) where we have chosen κ = κ g (+ 5)/, which is denoted here as the golden ratio, satisfying κ g κ g = by definition. We remark that the proposed system have three real eigenvales with linearly independent eigenvectors for any positive κ (ecept for κ =, which makes the matri of eigenvectors singlar). Ths, the proposed first order system is hyperbolic. With the choice of the golden ratio, we arrive at γ = /L ɛ. With the definitions of T ɛ, κ(= κ g ), β, and γ, we recast the proposed hyperbolic dispersion system in a vector form as τ U + A U = Q, (8) where U = p q, A = A disp = ɛ /T ɛ γ/t ɛ /T ɛ, Q = α t + s() (p + γ)/t ɛ q/t ɛ. (9) The three real eigenvales of the Jacobian matri A disp can be shown to be (see also Fig. ) = L ɛ κ g T ɛ, = κ g, 3 = κ g, () and the corresponding right and left eigenvectors are 3 of 3
4 t t 3 3 O O (a) ɛ < (b) ɛ > Figre : Schematics of the dispersion waves strctres for the proposed hyperbolic dispersion system. Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ R = κ g L ɛ L ɛ κ g L ɛ L ɛ κ g L ɛ L ɛ κ g, L = κ g + /L ɛ κ g/l ɛ κ g (κ g + )/L ɛ (κ g + )/L ɛ κ g + /L ɛ /(κ gl ɛ ) /κ g. () We now define the dispersion length scale, L ɛ, by first inserting the Forier mode of phase angle θ [, π] (i.e., U = U e Iθ/h ), where U is a constant vector, h is a grid spacing, and I =, into the hyperbolic dispersion system. We then evalate the eigenvales of the Forier-transformed operator with the κ g vale, and epand them for small θ to obtain the leading terms: λ = Iɛ h 3 θ3 + IɛL ɛ h 5 θ 5 + ɛl3 ɛ h 6 θ6 + O(θ 7 ), λ /3 = ɛ L 3 ɛ ± I ɛ L ɛh θ + O(θ ). () The first eigenvale is clearly an approimation to the Forier symbol of the operator ɛ, with a secondorder error in terms of θ. The other two represent damping modes with the real part, ɛ/l 3 ɛ, bt also have propagation in the net leading term, ±I ɛ θ/(l ɛh). To determine the dispersion length scale and possibly enhance the effects of error propagation dring transient psedo-time, we eqate the imaginary parts of the ɛ second and the third eigenvales with the leading term of the first eigenvale (i.e., h θ 3 = ɛ 3 L ɛh θ) for the smoothest mode (i.e., lowest freqency error mode) by setting θ = πh. After some algebra, we arrive at the following dispersion length scale: L ɛ = ± 4 π, (3) where the positive and negative length scales are associated with the positive and negative dispersion coefficient, ɛ, respectively. For all the nmerical eamples presented in Sec. VI, we observed eperimentally that both dobling and halving the above dispersion length scale reslt in a larger nmber of linear relaations and therefore, the above dispersion length scale appears to be optimm. The proposed hyperbolic dispersion system is now completely defined. In the net section, we present an etension of this system to an advective-dispersive PDE (sch as KdV) and show that a flly-hyperbolic system formlation of sch PDEs can also be fond. The general etension of the proposed hyperbolic system for a general advection-diffsion-dispersion eqation is discssed in Sec. IV. III. Hyperbolic Advection-Dispersion In this section, we briefly present an immediate etension of the hyperbolic dispersion system to general advective-dispersive PDEs (sch as KdV). Here, we present one possible flly-hyperbolic advective-dispersive system (in one-dimension), which cold be beneficial in imposing characteristic bondary conditions (BCs). The sbject of the characteristic BC is, however, beyond the scope of this paper and therefore, the effects of the formlation on characteristic BC is left for ftre investigations. Consider general advective-dispersive PDEs: t + f = ɛ, where advection speed is defined as 4 of 3
5 a() = f/. One possible hyperbolic formlation for advective-dispersive PDEs is: τ = a + ɛ q α t + s(), ( τ p = + ) ( p γ), T a T ɛ (4) τ q = T ɛ (γ + p q), Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ where T a is the advection relaation time, γ = β + the following three characteristic speeds λ adv-disp = β ( T a + T ɛ ), λ adv-disp,3 = T ɛ ɛ T β a T ɛ ɛ β T, and T = +. This system has T a T ɛ ( a λ adv-disp ± ) (a λ adv-disp ) + 4ɛβ, (5) T ɛ with the corresponding right eigenvectors epressed as ɛ λ adv-disp T λ adv-disp 3 T R = ɛβ. (6) λ adv-disp a βt/t ɛ λ adv-disp /(ɛβ) βt/t ɛ λ adv-disp 3 /(ɛβ) We follow the same procedre as given in Sec. II, and obtain the relaation time as T = (alt ɛ a L T ɛ + 4LT ɛ (βl + ))/(ɛβ). With this relaation time, we can show that the eigenvales are always real, and matri of right eigenvectors is non-singlar with κ > (ecept for κ = ). Ths, the presented system is flly-hyperbolic for κ > and κ. An optimm length scale L can be obtained by a Forier analysis, as discssed in Sec. II. However, this comptation is not performed here, becase for general nmerical approimation we se the hyperbolic strctre of advection, diffsion, and dispersion operators separately (see Sec. IV). The individal treatment of these operators makes the etension to mlti-dimensions and/or more comple eqations straightforward. IV. Generalization: Hyperbolic Advection-Diffsion-Dispersion Here, we present an etension of the proposed hyperbolic dispersion formlation to a general hyperbolic advection-diffsion-dispersion system. This etension is not trivial and therefore, is reported here for completeness. Consider a nonlinear advection-diffsion-dispersion eqation t + (f) = (ν ) + ɛ, (7) where f is a nonlinear fnction of, and the diffsion coefficient, ν, may be a fnction of the soltion variable, and the advection speed is defined as a() = f/. Using the hyperbolic dispersion system introdced in Sec. II, we propose the following hyperbolic system for a general advection-diffsion-dispersion eqation, Eq. (7), τ = f + (νp + γν) + ɛ q α + s(), (8) ( t τ p = + ) ( p γ), (9) T ν T ɛ τ q = T ɛ (γ + p q), () which is formlated sch that it properly redces to a pre hyperbolic advection-dispersion and a pre hyperbolic advection-diffsion system, respectively, in the dispersion limit (ν ) and the diffsion limit (ɛ ). 5 of 3
6 Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ Writing in vector form, Eq. (8), and following the non-nified approach of Ref. 3 for a separate treatment of advection, diffsion, and (in the this case) dispersion components, the fl Jacobian matri A of the proposed hyperbolic advection-diffsion-dispersion formlation becomes A = A adv + A diff + A disp = a γν + ν /T ν + A disp, () where A disp is given in Eq. (9). Note that for simplicity in the discssion, we presented the fl Jacobian for a linear ν, bt the proposed first-order system is also applicable for a nonlinear ν. Note also that we have inclded the fl Jacobian reslting from the presence of the γν term in the first eqation into the advective fl Jacobian, and not to the diffsion fl Jacobian; that is, we have treated γν as an added advection speed, which acts as an artificial advection in the diffsion limit (as γ does not vanish when ɛ ). With the above formlation, we also recover the identical hyperbolic diffsion system (which is in fact a [ ] system) of Ref. and ths, the same eigenvales, λ diff = ν/t ν, λ diff = ν/t ν, the diffsion relaation time, T ν = L ν/ν, and the diffsion length scale, L ν = /π, given in Ref.?, are applicable here. V. Discretization and Implicit Solver The proposed hyperbolic advection-diffsion-dispersion system can now be discretized with a desired scheme, sch as finite volme (FV), Discontinos-Galerkin (DG), RD, etc. Here we describe the discretization within the RD framework. Consider a one-dimensional domain discretized with N randomly distribted nodes. We store soltion vector U at each node denoted by j, where j =,, 3,..., N, and compte the cell residal, Φ E, by integrating the right-hand-side of Eq. (8) over the cell defined by the node j and j + (Fig. ): Φ E = j+ j ( A U + Q) d, () where we follow the techniqe of Ref. to compte cell residals of both the steady and nsteady sorce terms. j L E Figre : Schematic of the cell definition, and the left and the right nodes on each comptational cell. We now distribte the cell residal Φ E to the right and the left nodes of each cell with the following distribtion fnction, which is obtained sch that we can garantee conservation when advection, diffsion, and dispersion are treated separately: Φ L = B Φ E, Φ R = B + Φ E, B ± = I ± Dadv ± D diff ± D disp, (3) where the positive and the negative signs of the advection, D adv, the diffsion, D diff, and the dispersion, D disp, stabilization terms correspond to the distribtions of the cell residal Φ E to the right and the left nodes, respectively (see Fig. ). We define these stabilization terms independently by reformlating the advection, the diffsion, and the dispersion projection matrices in the form of the above distribtion fnction. To do this, we recall that an pwind distribtion matri can be constrcted as (after some algebra): B ± = R ± λ λ ± λ λ ± λ3 λ 3 R j L = I ± A ( R Λ L ) = I ± 3 A which reslts in the following advection, diffsion, and dispersion stabilization terms: D adv = Aadv /( a γν) + ɛ), D diff = Adiff l= λ diff l Πdiff l, D disp = Adisp l= l= 3 l λ l Π l, (4) Πdisp l, (5) 6 of 3
7 where ɛ is added to avoid division by zero when the total advection speed, a γν, is identically zero. Note that the distribtion matri sms p to the identity matri for each cell and therefore, we have conservation. The diffsion and dispersion projection matrices (i.e., Π diff and Π disp ), which are given here for convenience, are easily defined by projecting the diffsion and dispersion fles Jacobians onto their corresponding rnning waves (i.e., A diff. = l= λdiff l Π diff l, and A disp. = 3 l= λdisp l Π disp l ), to arrive at λ diff ν Π diff = /T ν, Π diff = I Π diff, (6) λ diff λ diff λ diff Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ Π disp = κ g + κ g κ 3 gl ɛ κ gl ɛ κ g L ɛ κ 4 g κ 3 gl ɛ L ɛ κ g L ɛ κ g, Πdisp = L ɛ L ɛ L ɛ L ɛ L ɛ L ɛ, Π disp 3 = I Π disp Π disp, (7) where I is an identity matri. We now follow the procedre described in Ref., and constrct an implicit solver with a second-order compact Jacobian matri, which we compte nmerically with an atomatic-differentiation techniqe based on an operator-overloading algorithm. VI. Reslts In this section, we eamine and verify the accracy of the proposed hyperbolic advection-diffsiondispersion system. The presented eamples are solved by discretizing the hyperbolic system sing the RD scheme discssed in Sec. V. We rela the linearized system with a Gass-Seidel (GS) algorithm and redce the linear residals by three orders of magnitde with the maimm of relaations. This typically takes, on average, abot 4 5 GS relaations with an nder relaation parameter of.8. The implicit solver is then contined ntil ten orders of magnitde redction is achieved for all the eqations; this sally takes abot 4 5 Newton iterations (per time step). We first present an eample for solving a nonlinear advection-diffsion-dispersion eqation throgh the method-of-manfactred soltion. This eample is sed to verify the order of accracy of the scheme. We then present two more eamples; one for solving the classical KdV eqation with small dispersion coefficient, and one to demonstrate the capability of the proposed scheme in resolving physical oscillations in the zero dispersion limit. A. Eample : Order of accracy verification Consider the following nonlinear advection-diffsion-dispersion eqation: t + f = (ν ) + ɛ + s(, ) (8) where f = 3, ν =.5, ɛ = and s is the manfactred sorce term. We seek, throgh the method of manfactred soltion, a time-dependent two-soliton soltion of the (generalized) KdV eqation: 6 e (, t) = (η η ) ( η sech [χ(η )] + η csch [χ(η )] ) ( η tanh [χ(η )] η coth [χ(η )] ), (9) where χ(η) = η ( η t ã), (3) and η >, η >, and ã are arbitrary constants. Here we consider η =.5, η =., and ã =., bt similar reslts are also obtained with other vales. After reformlating Eq. (8) in the form of the proposed hyperbolic advection-diffsion-dispersion system discssed in Sec. IV, we apply the efficient forth-order RD scheme of Ref. for the spatial discretization, along with the A-stable second-order Backward-Differencing-Formla (BDF) for the temporal discretization 7 of 3
8 to the hyperbolic system. We then solve the system of eqations with an imposed Dirichlet bondary condition (i.e., the time-depended soltions at the bondary nodes are fied) and e (, ) as an initial condition. Figre 3 shows the predicted soltion and its first and second gradients (i.e., and ) at different times. These soltions are obtained with t =. sing 6 randomly distribted grid points in (, 3). As shown, the predicted forth-order reslts are in a very good agreement with the eact vales even on.5 (eact) (eact) (eact).5 (eact) (eact) (eact) Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ (a) t = (c) t = 8. (eact) (eact) (eact) (b) t = 4. (eact) (eact) (eact) (d) t =. Figre 3: Two-soliton KdV with method of manfactred soltion: comparison with the predicted (forthorder spatial pls BDF) and eact soltions (inclding and ) at different times obtained on a randomly distribted nodes (N = 6) in [, 3] ( t =.). sch a relatively coarse grid. Figre 4 shows a comparison between the forth-order soltions and the second-order soltions for the same problem on the same randomly distribted grid. The presented reslts correspond to the soltion at t = with t =.. This figre shows that the second-order scheme over/nder predicts the peaks and valleys of the gradients, while the forth-order scheme captres them more accrately. We also verified the formal temporal and spatial orders of accracy for the soltion, gradient, and the Hessian (second-derivative) on a series of grids with randomly distribted nodes. The error convergence is 8 of 3
9 .5 Eact nd-order 4th-order.5 Eact nd-order 4th-order (,) (,) Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ (a) at t = (b) at t = Figre 4: Comparison between the predicted second- and forth-order soltion gradients for the two-soliton KdV problem compted on a randomly distribted nodes (N = 6) in [, 3] (t =, t =.). obtained by compting the L norms of the soltions (i.e., L = N j= U j U e j /N, where U is the soltion vector, the sperscript e denotes the eact vale, and N is the total nmber of grid points). Similar reslts are also obtained with the L norm and therefore not shown. For temporal accracy verification, we sed a grid with 64 randomly distribted points, and varied t. For spatial accracy, we fied t =. and sed a series of randomly distribted grids. These reslts are shown in Fig. 5, indicating that all soltion variables, inclding the first two soltion gradients (i.e., and ) are forth-order accrate in space and second-order accrate in time. This is significant becase a sith-order accrate soltion is typically needed with conventional schemes to achieve forth-order accrate second-derivatives. The reconstrcted soltion gradients, which are compted sing a cbic polynomial fit of the forth-order compted soltion, are at least orders of magnitde less accrate than the compted soltion gradients. The orders of accracy of the reconstrcted soltion gradients are also shown in Fig. 5 for comparison. B. Eample : Colliding soliton waves Consider the following classical KdV eqation: t + f = ɛ. (3) where f = / and ɛ <. The following initial condition reslts in a time-dependent colliding two-soliton soltion: ( ) ( ) (, ) = 3 η sech η /ɛ [( ) η t] + 3 η sech η /ɛ [( ) η t] (3) where η and η are arbitrary positive constants, and and are the initial locations of the two soliton waves. Following the eample given in Ref. 7, we consider η =.3, η =., =.4, =.8 and ɛ = with a periodic bondary condition, and solve the classical KdV eqation in [, ] with randomly distribted nodes. Similar reslts are also obtained with other vales and therefore, not shown. Figre 6 shows the predicted soltion and its gradient and Hessian (i.e., and ) at t =,, 3. The initial soltion is also shown as reference. These soltions are obtained with t =. sing the forth-order (spatial) RD scheme (with BDF for time discretization) of Ref.. We note that these reslts are converged and frther grid refinement does not change the reslts qalitatively. 9 of 3
10 Hyperbolic Advection Diffsion Dispersion (Randomly Distribted Grid) Hyperbolic Advection Diffsion Dispersion (Randomly Distribted Grid) Slope Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ log L Slope log t (a) Temporal accracy (BDF), N = 64 log L Slope 3 Slope 4 Slope log h (reconstrcted) (reconstrcted) (b) Spatial accracy (forth-order), t =. Figre 5: Temporal and spatial accracy verifications of the proposed hyperbolic advection-diffsiondispersion system for the soltion variable and its first two derivatives sing the forth-order RD-GT scheme (with BDF for temporal discretization) on randomly distribted nodes. C. Eample 3: Dispersive shock Consider a dispersive nonlinear Brger eqation (i.e., Eq. (3) when ɛ + ). This eqation has a dispersive behavior and prodces continos wavelets in the vicinity of the discontinity (aka dispersive shock waves). In this eample, we illstrate the capability of or code in resolving high-freqency wavelets for very small ɛ. Figre 7 shows the predicted reslts for ɛ = 3 and 4. The predicted high-order gradients and Hessians (second-derivatives) of the soltion are also provided. The reslts show that physical oscillations are captred and soltions are noise-free particlarly before and after the dispersive shocks (others reported soltions with some nmerical oscillations, see e.g., Ref. 4). For comparison, the reconstrcted gradient (third-order) and Hessian (second-order) are also shown for comparison. We note that these soltions are converged (i.e., no noticeable change in reslts are observed with frther grid refinement) and we eperienced no instability in obtaining these soltions. VII. Conclsions We have introdced, for the first time, a first-order system approach for general advection-diffsiondispersion eqations. We showed that the proposed system has real eigenvales with linearly independent eigenvectors, and ths is hyperbolic. We apply the forth-order RD scheme of Ref. to the proposed hyperbolic system, and solved several time-dependent dispersive eqations, inclding the classical two-soliton KdV and a dispersive shock case. We demonstrated that the predicted reslts, inclding the gradient and Hessian (i.e., and ), are in a very good agreement with the eact soltions. We also showed that the RD scheme applied to the proposed system can captre dispersive shocks with no nmerical oscillations. The design order of accracy of the scheme (forth-order spatial and second-order temporal) is also verified, and we achieved eqal order of accracy for the soltion, the gradient, and the Hessian. Acknowledgments The first and last athors wold like to thank the Center Chief Technology Office of NASA Langley Research Center for their spport throgh the Center Innovation Fnd (CIF). of 3
11 (a) t =, (b) t =, (c) t =, Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ (d) t =, (g) t =, (e) t =, (h) t =, (f) t =, (i) t =, (j) t = 3, (k) t = 3, (l) t = 3, Figre 6: Two-soliton soltion to t + = ɛ : Predicted (forth-order spatial pls BDF) reslts with randomly distribted grid points in a periodic domain [, ] (ɛ = , t =.). of 3
12 Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ (a) N = 3, t =.5, ɛ = (d) N = 5, t =.5, ɛ = 5 5 (reconstrcted) (b) N = 3, t =.5, ɛ = 4 5 (reconstrcted) (e) N = 5, t =.5, ɛ = 5 (reconstrcted) (c) N = 3, t =.5, ɛ = 4 5 (reconstrcted) (f) N = 5, t =.5, ɛ = 5 Figre 7: Dispersive shock case (zero dispersion limit): t + = ɛ, ɛ, (, ) = + / sin(π): predicted soltion, its gradient and Hessian (second-derivative) sing the forth-order (pls BDF) RD scheme applied to the proposed hyperbolic dispersion system on randomly distribted grid points in a periodic domain [, ]. of 3
13 References Downloaded by NASA LANGLEY RESEARCH CENTRE on Agst 4, 6 DOI:.54/ A. Mazaheri and H. Nishikawa. Very efficient high-order hyperbolic schemes for time-dependent advection-diffsion problems: Third-, forth-, and sith-order. Compters and Flids, :3 47, 4. H. Leblond. Interaction of two solitary waves in a ferromagnet. J. Phys. A: Math. Gen., 8:3763, W. Wan, S. Jia, and W. Fleischer. Dispersive sperflid-like shock waves in nonlinear optics. Natre Physics, 3():46 5, 7. 4 N. Ghofraniha, C. Conti, G. Rocco, and S. Trillo. Shocks in nonlocal media. Physical Review Letters, 99:4393, 7. 5 M.A. Hoefer, M.J. Ablowitz, I. Coddington, E.A. Cornell, P. Engels, and V. Schweikhard. Dispersive and classical shock waves in Bose-Einstein condensates and gas dynamics. Phys. Rev. A, 74:363, 6. 6 M.A. Hoefer. Shock waves in dispersive Elerian flids. Jornal of Nonlinear Science, 4(3):55 577, 4. 7 R.K. Smith, N. Crook, and G. Roff. The morning glory: an etraordinary atmospheric ndlar bore. Qart. J. R. Met. Soc., 8: , 98. J.W. Rottman and R. Grimshaw. Atmospheric internal solitary waves. In Roger Grimshaw, editor, Environmental Stratified Flows, volme 3 of Topics in Environmental Flid Mechanics, pages Springer US,. 9 S. Glimsdal, G.K. Pedersen, C.B. Harbitz, and F. Lovholt. Dispersion of tsnamis: does it really matter? Nat. Hazards Earth Syst. Sci., 3:57 56, 3. H. Chanson. Crrent knowledge in hydralic jmps and related phenomena. a srvey of eperimental reslts. Eropean Jornal of Mechanics - B/Flids, 8():9, 9. D.L. Wilkinson and M.L. Banner. Undlar bores. 6th Astralian Hydralics and Flid Mechanics Conference, Adelaide (Astralia), December 977. N.K. Lowman and M.A. Hoefer. Dispersive shock waves in viscosly deformable media. J. Flid Mech., 78:54 557, 3. 3 S. Benzoni-Gavage. Planar traveling waves in capillary flids. Differential Integral Eqations, 6(3-4): , 3. 4 C.-Y. Kao, A. Krganov, and Z. Q. A fast eplicit operator splitting method for modified Bckley-Leverett eqations. J. Sci. Compt., doi:.7/s , 4. 5 R. Carles, R. Danchin, and J.-C. Sat. Madelng, Gross-Pitaevskii and Korteweg. Nonlinearity, 5(): ,. 6 C.-K. Lin and K.C. W. Hydrodynamic limits of the nonlinear Klein-Gordon eqation. J. Math. Pres Appl., 98:38 345,. 7 H.-L. Li and J. Yan. A local discontinos Galerkin method for the Korteweg-de Vries eqation with bondary effect. J. Compt. Phys., 5:97 8, 6. 8 D. Lannes. The Water Waves Problem: Mathematical Analysis and Asymptotics. American Mathematical Society, Providence, Rhode Island, 3. 9 A.G. Filippini, S. Bellec, M. Colin, and M. Ricchito. On the nonlinear behavior of Bossinesq type models: amplitdevelocity vs ampltitde-fl forms. Coastal Engng., doi:.6/j.coastaleng , 5. H. Nishikawa. A first-order system approach for diffsion eqation. I: Second order residal distribtion schemes. J. Compt. Phys., 7:35 35, 7. A. Mazaheri and H. Nishikawa. First-order hyperbolic system method for time-dependent advection-diffsion problems. Technical Report TM 4 875, NASA, March 4. H. Nishikawa. Accracy-preserving bondary fl qadratre for finite-volme discretization on nstrctred grids. J. Compt. Phys., 8:58 555, 5. 3 A. Mazaheri and H. Nishikawa. Improved second-order hyperbolic residal-distribtion scheme and its etension to third-order on arbitrary trianglar grids. J. Compt. Phys., 3:455 49, 5. 4 J. Yan and C.-W. Sh. A local discontinos Galerkin method for KdV type eqations. SIAM J. Nmer. Anal., 4:769 79,. 5 E.F. Toro and G.I. Montecinos. Advection-diffsion-reaction eqation: hyperbolization and high-order ADER discritizations. SIAM J. Sci. Compt., 36(5):A43 A457, 4. 6 K. Braer. Accessed: Jly, 5. 7 A. Debssche and J. Printems. Nmerical simlation of the stochastic Kortewegde Vries eqation. Physica D: Nonlinear Phenomena, 34: 6, of 3
Active 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 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 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 informationFREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS
7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,
More informationBurgers Equation. A. Salih. Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram 18 February 2016
Brgers Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 18 Febrary 216 1 The Brgers Eqation Brgers eqation is obtained as a reslt of
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 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 informationfür Mathematik in den Naturwissenschaften Leipzig
Ma-Planck-Institt für Mathematik in den Natrwissenschaften Leipzig Nmerical simlation of generalized KP type eqations with small dispersion by Christian Klein, and Christof Sparber Preprint no.: 2 26 NUMERICAL
More informationSimplified Identification Scheme for Structures on a Flexible Base
Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles
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 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 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 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 information1. Introduction. In this paper, we are interested in accurate numerical approximations to the nonlinear Camassa Holm (CH) equation:
SIAM J. SCI. COMPUT. Vol. 38, No. 4, pp. A99 A934 c 6 Society for Indstrial and Applied Matematics AN INVARIANT PRESERVING DISCONTINUOUS GALERKIN METHOD FOR THE CAMASSA HOLM EQUATION HAILIANG LIU AND YULONG
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 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 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 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 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 informationModelling by Differential Equations from Properties of Phenomenon to its Investigation
Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University
More informationReduction of over-determined systems of differential equations
Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical
More informationOn the scaling of entropy viscosity in high order methods
On te scaling of entropy viscosity in ig order metods Adeline Kornels and Daniel Appelö Abstract In tis work, we otline te entropy viscosity metod and discss ow te coice of scaling inflences te size of
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 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 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 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 informationSolving the Lienard equation by differential transform method
ISSN 1 746-7233, England, U World Jornal of Modelling and Simlation Vol. 8 (2012) No. 2, pp. 142-146 Solving the Lienard eqation by differential transform method Mashallah Matinfar, Saber Rakhshan Bahar,
More informationThe Scalar Conservation Law
The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan
More informationHybrid modelling and model reduction for control & optimisation
Hybrid modelling and model redction for control & optimisation based on research done by RWTH-Aachen and TU Delft presented by Johan Grievink Models for control and optimiation market and environmental
More informationNew Fourth Order Explicit Group Method in the Solution of the Helmholtz Equation Norhashidah Hj. Mohd Ali, Teng Wai Ping
World Academy of Science, Engineering and Tecnology International Jornal of Matematical and Comptational Sciences Vol:9, No:, 05 New Fort Order Eplicit Grop Metod in te Soltion of te elmoltz Eqation Norasida.
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 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 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 informationAdjoint-Based Sensitivity Analysis for Computational Fluid Dynamics
Adjoint-Based Sensitivity Analysis for Comptational Flid Dynamics Dimitri J. Mavriplis Department of Mecanical Engineering niversity of Wyoming Laramie, WY Motivation Comptational flid dynamics analysis
More informationAN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD FOR BUCKLING AND INITIAL POST-BUCKLING ANALYSIS OF COMPOSITE SHELLS
th Eropean Conference on Comptational Mechanics (ECCM ) 7th Eropean Conference on Comptational Flid Dynamics (ECFD 7) 5 Jne 28, Glasgow, UK AN ISOGEOMETRIC SOLID-SHELL FORMULATION OF THE KOITER METHOD
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 informationComputational Geosciences 2 (1998) 1, 23-36
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), 23-36 Abstract.
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 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 informationBridging the Gap Between Multigrid, Hierarchical, and Receding-Horizon Control
Bridging the Gap Between Mltigrid, Hierarchical, and Receding-Horizon Control Victor M. Zavala Mathematics and Compter Science Division Argonne National Laboratory Argonne, IL 60439 USA (e-mail: vzavala@mcs.anl.gov).
More informationLewis number and curvature effects on sound generation by premixed flame annihilation
Center for Trblence Research Proceedings of the Smmer Program 2 28 Lewis nmber and crvatre effects on sond generation by premixed flame annihilation By M. Talei, M. J. Brear AND E. R. Hawkes A nmerical
More informationConditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane
Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching
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 informationA Model-Free Adaptive Control of Pulsed GTAW
A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,
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 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 informationEfficient quadratic penalization through the partial minimization technique
This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC272754474, IEEE Transactions
More information3 2D Elastostatic Problems in Cartesian Coordinates
D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates
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 informationA Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane
Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani
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 informationLecture Notes: Finite Element Analysis, J.E. Akin, Rice University
9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)
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 informationModeling Effort on Chamber Clearing for IFE Liquid Chambers at UCLA
Modeling Effort on Chamber Clearing for IFE Liqid Chambers at UCLA Presented by: P. Calderoni own Meeting on IFE Liqid Wall Chamber Dynamics Livermore CA May 5-6 3 Otline his presentation will address
More informationRestricted Three-Body Problem in Different Coordinate Systems
Applied Mathematics 3 949-953 http://dx.doi.org/.436/am..394 Pblished Online September (http://www.scirp.org/jornal/am) Restricted Three-Body Problem in Different Coordinate Systems II-In Sidereal Spherical
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 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 informationBäcklund transformation, multiple wave solutions and lump solutions to a (3 + 1)-dimensional nonlinear evolution equation
Nonlinear Dyn 7 89:33 4 DOI.7/s7-7-38-3 ORIGINAL PAPER Bäcklnd transformation, mltiple wave soltions and lmp soltions to a 3 + -dimensional nonlinear evoltion eqation Li-Na Gao Yao-Yao Zi Y-Hang Yin Wen-Xi
More information1 Differential Equations for Solid Mechanics
1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to
More informationA local discontinuous Galerkin method for the Camassa-Holm equation
A local discontinos Galerkin method for the Camassa-Holm eqation Yan X and Chi-Wang Sh Janary, 7 Abstract In this paper, we develop, analyze and test a local discontinos Galerkin (LDG) method for solving
More informationComparison of Monte Carlo and deterministic simulations of a silicon diode Jose A. Carrillo Λ, Irene M. Gamba y, Orazio Muscato z and Chi-Wang Shu x N
Comparison of Monte Carlo and deterministic simlations of a silicon diode Jose A. Carrillo Λ, Irene M. Gamba y, Orazio Mscato z and Chi-ang Sh November 7, Abstract lectron transport models in Si transistors
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 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 informationMixed Type Second-order Duality for a Nondifferentiable Continuous Programming Problem
heoretical Mathematics & Applications, vol.3, no.1, 13, 13-144 SSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Mied ype Second-order Dality for a Nondifferentiable Continos Programming Problem.
More informationSimilarity Solution for MHD Flow of Non-Newtonian Fluids
P P P P IJISET - International Jornal of Innovative Science, Engineering & Technology, Vol. Isse 6, Jne 06 ISSN (Online) 48 7968 Impact Factor (05) - 4. Similarity Soltion for MHD Flow of Non-Newtonian
More informationOnline Solution of State Dependent Riccati Equation for Nonlinear System Stabilization
> REPLACE American HIS Control LINE Conference WIH YOUR PAPER IDENIFICAION NUMBER (DOUBLE-CLICK HERE O EDI) FrC3. Marriott Waterfront, Baltimore, MD, USA Jne 3-Jly, Online Soltion of State Dependent Riccati
More informationOn peaked solitary waves of the Degasperis-Procesi equation
. Article. Special Topic: Flid Mechanics SCIENCE CHINA Physics, Mechanics & Astronomy Febrary 213 Vol. 56 No. 2: 418 422 doi: 1.17/s11433-13-4993-9 On peaked solitary waves of the Degasperis-Procesi eqation
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 informationApproximate Solution for the System of Non-linear Volterra Integral Equations of the Second Kind by using Block-by-block Method
Astralian Jornal of Basic and Applied Sciences, (1): 114-14, 008 ISSN 1991-8178 Approximate Soltion for the System of Non-linear Volterra Integral Eqations of the Second Kind by sing Block-by-block Method
More informationFRÉCHET KERNELS AND THE ADJOINT METHOD
PART II FRÉCHET KERNES AND THE ADJOINT METHOD 1. Setp of the tomographic problem: Why gradients? 2. The adjoint method 3. Practical 4. Special topics (sorce imaging and time reversal) Setp of the tomographic
More informationGravity-capillary waves on the free surface of a liquid dielectric in a tangential electric field
Gravity-capillary waves on the free srface of a liqid dielectric in a tangential electric field Evgeny A. Kochrin Institte of Electrophysics, Ural Branch of Rssian Academy of Sciences 66, 6 Amndsen str.,
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 informationPartial Differential Equations with Applications
Universit of Leeds MATH 33 Partial Differential Eqations with Applications Eamples to spplement Chapter on First Order PDEs Eample (Simple linear eqation, k + = 0, (, 0) = ϕ(), k a constant.) The characteristic
More informationChaotic and Hyperchaotic Complex Jerk Equations
International Jornal of Modern Nonlinear Theory and Application 0 6-3 http://dxdoiorg/0436/ijmnta000 Pblished Online March 0 (http://wwwscirporg/jornal/ijmnta) Chaotic and Hyperchaotic Complex Jerk Eqations
More informationNonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator
Proceedings of the 6 IEEE International Conference on Control Applications Mnich, Germany, October 4-6, 6 WeB16 Nonparametric Identification and Robst H Controller Synthesis for a Rotational/Translational
More informationReflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by
Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introdction The transmission line eqations are given by, I z, t V z t l z t I z, t V z, t c z t (1) (2) Where, c is the per-nit-length
More informationLarge time-stepping methods for higher order time-dependent evolution. equations. Xiaoliang Xie. A dissertation submitted to the graduate faculty
Large time-stepping methods for higher order time-dependent evoltion eqations by Xiaoliang Xie A dissertation sbmitted to the gradate faclty in partial flfillment of the reqirements for the degree of DOCTOR
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 informationTime-adaptive non-linear finite-element analysis of contact problems
Proceedings of the 7th GACM Colloqim on Comptational Mechanics for Yong Scientists from Academia and Indstr October -, 7 in Stttgart, German Time-adaptive non-linear finite-element analsis of contact problems
More informationUttam Ghosh (1), Srijan Sengupta (2a), Susmita Sarkar (2b), Shantanu Das (3)
Analytical soltion with tanh-method and fractional sb-eqation method for non-linear partial differential eqations and corresponding fractional differential eqation composed with Jmarie fractional derivative
More informationStability of Model Predictive Control using Markov Chain Monte Carlo Optimisation
Stability of Model Predictive Control sing Markov Chain Monte Carlo Optimisation Elilini Siva, Pal Golart, Jan Maciejowski and Nikolas Kantas Abstract We apply stochastic Lyapnov theory to perform stability
More information3.1 The Basic Two-Level Model - The Formulas
CHAPTER 3 3 THE BASIC MULTILEVEL MODEL AND EXTENSIONS In the previos Chapter we introdced a nmber of models and we cleared ot the advantages of Mltilevel Models in the analysis of hierarchically nested
More informationA Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model
entropy Article A Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model Hyenkyn Woo School of Liberal Arts, Korea University of Technology and Edcation, Cheonan
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 informationInternational Journal of Scientific & Engineering Research, Volume 5, Issue 3, March ISSN
International Jornal of Scientific & Engineering Research, Volme 5, Isse 3, March-4 83 ISSN 9-558 Doble Dispersion effects on free convection along a vertical Wavy Srface in Poros Media with Variable Properties
More informationAN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS
AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS Fang-Ming Y, Hng-Yan Chng* and Chen-Ning Hang Department of Electrical Engineering National Central University, Chngli,
More informationAdvances in Water Resources
Advances in Water Resorces () Contents lists available at ScienceDirect Advances in Water Resorces ornal homepage: www.elsevier.com/locate/advwatres Positivity-preserving high order well-balanced discontinos
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 informationAN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS
AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS Fang-Ming Y, Hng-Yan Chng* and Chen-Ning Hang Department of Electrical Engineering National Central University, Chngli,
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 informationTransient Approach to Radiative Heat Transfer Free Convection Flow with Ramped Wall Temperature
Jornal of Applied Flid Mechanics, Vol. 5, No., pp. 9-1, 1. Available online at www.jafmonline.net, ISSN 175-57, EISSN 175-645. Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped
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 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 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 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 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 informationReducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation
JOURNAL OF AIRCRAFT Vol. 42, No. 4, Jly Agst 2005 Redcing Conservatism in Fltterometer Predictions Using Volterra Modeling with Modal Parameter Estimation Rick Lind and Joao Pedro Mortaga University of
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 informationGravitational Instability of a Nonrotating Galaxy *
SLAC-PUB-536 October 25 Gravitational Instability of a Nonrotating Galaxy * Alexander W. Chao ;) Stanford Linear Accelerator Center Abstract Gravitational instability of the distribtion of stars in a galaxy
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 information