Open Access. Powered by Scholars. Published by Universities.®
- Discipline
- Institution
- Publication
- Publication Type
Articles 1 - 9 of 9
Full-Text Articles in Analysis
Optimal Control Of Stochastic Systems: A Numerical Study Of The Stochastic Linear Quadratic Regulator Framework, Gülşen Orucova Büyüköz, Yaprak Güldoğan Dericioğlu, Tuğçem Partal
Optimal Control Of Stochastic Systems: A Numerical Study Of The Stochastic Linear Quadratic Regulator Framework, Gülşen Orucova Büyüköz, Yaprak Güldoğan Dericioğlu, Tuğçem Partal
Mathematical Modelling and Numerical Simulation with Applications
Optimal control of stochastic linear systems is fundamental in control theory, with applications in robotics, finance, and engineering. The Stochastic Linear Quadratic Regulator (SLQR) derives optimal feedback laws via the Riccati equation but requires numerical discretization of the resulting stochastic dynamics. Despite extensive studies on numerical methods for stochastic differential equations, their performance within the SLQR framework remains insufficiently explored. This study compares two predictor–corrector schemes of different orders: the Order 1.0 Predictor-Corrector (PC) method and the Order 2.0 Weak PC method. A one-dimensional linear quadratic problem with a closed-form solution enables precise error evaluation against the analytical trajectory. Convergence …
Stable And Convergent Difference Schemes For Weakly Singular Convolution Integrals, Wesley Davis, Richard D. Noren
Stable And Convergent Difference Schemes For Weakly Singular Convolution Integrals, Wesley Davis, Richard D. Noren
Mathematics & Statistics Faculty Publications
We obtain new numerical schemes for weakly singular integrals of convolution type called Caputo fractional order integrals using Taylor and fractional Taylor series expansions and grouping terms in a novel manner. A fractional Taylor series expansion argument is utilized to provide fractional-order approximations for functions with minimal regularity. The resulting schemes allow for the approximation of functions in Cγ [0, T], where 0 < γ <= 5. A mild invertibility criterion is provided for the implicit schemes. Consistency and stability are proven separately for the whole-number-order approximations and the fractional-order approximations. The rate of convergence in the time variable is shown …
Numerical Simulation For A Rising Bubble Interacting With A Solid Wall: Impact, Bounce, And Thin Film Dynamics, Changjuan Zhang, Jie Li, Li-Shi Luo, Tiezheng Qian
Numerical Simulation For A Rising Bubble Interacting With A Solid Wall: Impact, Bounce, And Thin Film Dynamics, Changjuan Zhang, Jie Li, Li-Shi Luo, Tiezheng Qian
Mathematics & Statistics Faculty Publications
Using an arbitrary Lagrangian-Eulerian method on an adaptive moving unstructured mesh, we carry out numerical simulations for a rising bubble interacting with a solid wall. Driven by the buoyancy force, the axisymmetric bubble rises in a viscous liquid toward a horizontal wall, with impact on and possible bounce from the wall. First, our simulation is quantitatively validated through a detailed comparison between numerical results and experimental data. We then investigate the bubble dynamics which exhibits four different behaviors depending on the competition among the inertial, viscous, gravitational, and capillary forces. A phase diagram for bubble dynamics has been produced using …
Numerical Analysis And Design Optimization Of Organic Bulk Hetero-Junction Solar Cells Using The Mixed Finite Element And Arc-Length Methods, James L. Dean
Mechanical & Aerospace Engineering Theses & Dissertations
Due to the high coupling and non-linear nature of the semiconductor equations, numerical stability becomes a major problem when using a purely Newtonian numerical method in solving these equations. Presented is a robust mixed finite element formulation with application for organic solar cell devices that incorporates the arc-length predictor-corrector method to bypass such stability issues. The simulated device represents a multilayered solar cell using organic bulk hetero-junction materials including donor materials Poly[2,1,3-benzothiadiazole-4,7-diyl[4,4-bis(2-ethylhexyl)-4H-cyclopenta[2,1-b:3,4- b']dithiophene-2,6-diyl]] (PCPDTBT) and poly 3-hexylthiophene (P3HT) with acceptor fullerene materials [6,6]-phenyl C61 butyric acid methyl ester (PC61BM) or [6,6]-phenyl C71 butyric acid methyl ester (PC71BM) as the photoactive …
Polynomial Extension Operators. Part Ii, Leszek Demkowicz, Jay Gopalakrishnan, Joachim Schöberl
Polynomial Extension Operators. Part Ii, Leszek Demkowicz, Jay Gopalakrishnan, Joachim Schöberl
Mathematics and Statistics Faculty Publications and Presentations
Consider the tangential trace of a vector polynomial on the surface of a tetrahedron. We construct an extension operator that extends such a trace function into a polynomial on the tetrahedron. This operator can be continuously extended to the trace space of H(curl ). Furthermore, it satisfies a commutativity property with an extension operator we constructed in Part I of this series. Such extensions are a fundamental ingredient of high order finite element analysis.
Multigrid Convergence For Second Order Elliptic Problems With Smooth Complex Coefficients, Jay Gopalakrishnan, Joseph E. Pasciak
Multigrid Convergence For Second Order Elliptic Problems With Smooth Complex Coefficients, Jay Gopalakrishnan, Joseph E. Pasciak
Mathematics and Statistics Faculty Publications and Presentations
The finite element method when applied to a second order partial differential equation in divergence form can generate operators that are neither Hermitian nor definite when the coefficient function is complex valued. For such problems, under a uniqueness assumption, we prove the continuous dependence of the exact solution and its finite element approximations on data provided that the coefficients are smooth and uniformly bounded away from zero. Then we show that a multigrid algorithm converges once the coarse mesh size is smaller than some fixed number, providing an efficient solver for computing discrete approximations. Numerical experiments, while confirming the theory, …
Locally Conservative Fluxes For The Continuous Galerkin Method, Bernardo Cockburn, Jay Gopalakrishnan, Haiying Wang
Locally Conservative Fluxes For The Continuous Galerkin Method, Bernardo Cockburn, Jay Gopalakrishnan, Haiying Wang
Mathematics and Statistics Faculty Publications and Presentations
The standard continuous Galerkin (CG) finite element method for second order elliptic problems suffers from its inability to provide conservative flux approximations, a much needed quantity in many applications. We show how to overcome this shortcoming by using a two step postprocessing. The first step is the computation of a numerical flux trace defined on element inter- faces and is motivated by the structure of the numerical traces of discontinuous Galerkin methods. This computation is non-local in that it requires the solution of a symmetric positive definite system, but the system is well conditioned independently of mesh size, so it …
Error Analysis Of Variable Degree Mixed Methods For Elliptic Problems Via Hybridization, Bernardo Cockburn, Jay Gopalakrishnan
Error Analysis Of Variable Degree Mixed Methods For Elliptic Problems Via Hybridization, Bernardo Cockburn, Jay Gopalakrishnan
Mathematics and Statistics Faculty Publications and Presentations
A new approach to error analysis of hybridized mixed methods is proposed and applied to study a new hybridized variable degree Raviart-Thomas method for second order elliptic problems. The approach gives error estimates for the Lagrange multipliers without using error estimates for the other variables. Error estimates for the primal and flux variables then follow from those for the Lagrange multipliers. In contrast, traditional error analyses obtain error estimates for the flux and primal variables first and then use it to get error estimates for the Lagrange multipliers. The new approach not only gives new error estimates for the new …
A Characterization Of Hybridized Mixed Methods For Second Order Elliptic Problems, Bernardo Cockburn, Jay Gopalakrishnan
A Characterization Of Hybridized Mixed Methods For Second Order Elliptic Problems, Bernardo Cockburn, Jay Gopalakrishnan
Mathematics and Statistics Faculty Publications and Presentations
In this paper, we give a new characterization of the approximate solution given by hybridized mixed methods for second order self-adjoint elliptic problems. We apply this characterization to obtain an explicit formula for the entries of the matrix equation for the Lagrange multiplier unknowns resulting from hybridization. We also obtain necessary and sufficient conditions under which the multipliers of the Raviart–Thomas and the Brezzi–Douglas–Marini methods of similar order are identical.