Numerical Analysis and Computation Commons^™

Validation Of Weak Form Thermal Analysis Algorithms Supporting Thermal Signature Generation, Elton Lewis Freeman

Masters Theses

Extremization of a weak form for the continuum energy conservation principle differential equation naturally implements fluid convection and radiation as flux Robin boundary conditions associated with unsteady heat transfer. Combining a spatial semi-discretization via finite element trial space basis functions with time-accurate integration generates a totally node-based algebraic statement for computing. Closure for gray body radiation is a newly derived node-based radiosity formulation generating piecewise discontinuous solutions, while that for natural-forced-mixed convection heat transfer is extracted from the literature. Algorithm performance, mathematically predicted by asymptotic convergence theory, is subsequently validated with data obtained in 24 hour diurnal field experiments for ...

Stable, Tunable, Quasimonoenergetic Electron Beams Produced In A Laser Wakefield Near The Threshold For Self-Injection, Sudeep Banerjee, Serguei Y. Kalmykov, Nathan D. Powers, Gregory Golovin, Vidiya Ramanathan, Nathan J. Cunningham, Kevin J. Brown, Shouyuan Chen, Isaac Ghebregziabher, Bradley A. Shadwick, Donald P. Umstadter, Benjamin A. Cowan, David L. Bruhwiler, Arnaud Beck, Erik Lefebvre

Serguei Y. Kalmykov

Stable operation of a laser-plasma accelerator near the threshold for electron self-injection in the blowout regime has been demonstrated with 25–60 TW, 30 fs laser pulses focused into a 3–4 millimeter length gas jet. Nearly Gaussian shape and high nanosecond contrast of the focused pulse appear to be critically important for controllable, tunable generation of 250–430 MeV electron bunches with a low energy spread, ~ 10 pC charge, a few-mrad divergence and pointing stability, and a vanishingly small low-energy background. The physical nature of the near-threshold behavior is examined using three-dimensional particle-in-cell simulations. Simulations indicate that properly locating ...

Two-Dimensional Hydrodynamic Modeling Of Two-Phase Flow For Understanding Geyser Phenomena In Urban Stormwater System, Zhiyu S. Shao

Theses and Dissertations--Civil Engineering

During intense rain events a stormwater system can fill rapidly and undergo a transition from open channel flow to pressurized flow. This transition can create large discrete pockets of trapped air in the system. These pockets are pressurized in the horizontal reaches of the system and then are released through vertical vents. In extreme cases, the transition and release of air pockets can create a geyser feature.

The current models are inadequate for simulating mixed flows with complicated air-water interactions, such as geysers. Additionally, the simulation of air escaping in the vertical dropshaft is greatly simplified, or completely ignored, in ...

Propeller, Joel Kahn

The STEAM Journal

This image is based on several different algorithms interconnected within a single program in the language BASIC-256. The fundamental structure involves a tightly wound spiral working outwards from the center of the image. As the spiral is drawn, different values of red, green and blue are modified through separate but related processes, producing the changing appearance. Algebra, trigonometry, geometry, and analytic geometry are all utilized in overlapping ways within the program. As with many works of algorithmic art, small changes in the program can produce dramatic alterations of the visual output, which makes lots of variations possible.

Sloane’S Gap: Do Mathematical And Social Factors Explain The Distribution Of Numbers In The Oeis?, Nicolas J.-P. Gauvrit, Jean-Paul Delahaye, Hector Zenil

Journal of Humanistic Mathematics

The Online Encyclopedia of Integer Sequences (OEIS) is a catalog of integer sequences. We are particularly interested in the number of occurrences of N(n) of an integer n in the database. This number N(n) marks the importance of n and it varies noticeably from one number to another, and from one number to the next in a series. “Importance” can be mathematically objective (2^10 is an example of an “important” number in this sense) or as the result of a shared mathematical culture (10^9 is more important than 9^10 because we use a decimal notation ...

Petawatt-Laser-Driven Wakefield Acceleration Of Electrons To 2 Gev In 10^{17} Cm^{-3} Plasma, Xiaoming Wang, Rafal B. Zgadzaj, Neil Fazel, Sunghwan A. Yi, X. Zhang, Watson Henderson, Yen-Yu Zhang, Rick Korzekwa, Hai-En Tsai, C.-H. Pai, Zhengyan Li, Hernan Quevedo, Gilliss Dyer, Erhard W. Gaul, Mikael Martinez, Aaron Bernstein, Ted Borger, M. Spinks, M. Donovan, Serguei Y. Kalmykov, Vladimir N. Khudik, Gennady Shvets, Todd Ditmire, Michael C. Downer

Serguei Y. Kalmykov

Electron self-injection into a laser-plasma accelerator (LPA) driven by the Texas Petawatt (TPW) laser is reported at plasma densities 1.7 - 6.2 x 10^{17} cm^{-3}. Energy and charge of the electron beam, ranging from 0.5 GeV to 2 GeV and tens to hundreds of pC, respectively, depended strongly on laser beam quality and plasma density. Angular beam divergence was consistently around 0.5 mrad (FWHM), while shot-to-shot pointing fluctuations were limited to ±1.4 mrad rms. Betatron x-rays with tens of keV photon energy are also clearly observed.

Sub-Millimeter-Scale, 100-Mev-Class Quasi-Monoenergetic Laser Plasma Accelerator Based On All-Optical Control Of Dark Current In The Blowout Regime, Serguei Y. Kalmykov, Xavier Davoine, Bradley A. Shadwick

Serguei Y. Kalmykov

It is demonstrated that by negatively chirping the frequency of a 20-fs, 15-TW driving laser pulse with an ultrabroad bandwidth (corresponding to a sub-2-cycle transform-limited duration it is possible to prevent early compression of the pulse into an optical shock, thus reducing expansion of the accelerating plasma bucket (electron density "bubble") and delaying dephasing of self-injected and accelerated electrons. These features help suppress unwanted continuous self-injection (dark current) in the blowout regime, making possible to use the entire dephasing length to generate low-background, quasi-monoenergetic 200-MeV-scale electron beams from sub-mm-length, dense plasmas (n_{e0} = 1.3 x 10^{19} cm^{−3}).

Dark-Current-Free Laser-Plasma Acceleration In Blowout Regime Using Nonlinear Plasma Lens, Serguei Y. Kalmykov

Serguei Y. Kalmykov

It is demonstrated that a thin dense plasma slab (lens), placed before a multi-centimeter-length, low-density plasma (accelerator), overfocuses an incident petawatt laser pulse at a controlled location inside the accelerator, creating an expanding electron density bubble that traps plasma electrons over a brief time interval. As soon as the pulse stabilizes and self-guiding begins, the bubble stabilizes and transforms into the first (non-broken) bucket of a conventional three-dimensional nonlinear plasma wave, eliminating any chance for further injection. A well collimated, quasi-monoenergetic electron bunch with a zero low-energy background further accelerates to a multi-GeV energy.

Computationally Efficient Methods For Modelling Laser Wakefield Acceleration In The Blowout Regime, Benjamin M. Cowan, Serguei Y. Kalmykov, Arnaud Beck, Xavier Davoine, Kyle Bunkers, Agustin F. Lifschitz, Erik Lefebvre, David L. Bruhwiler, Bradley A. Shadwick, Donald P. Umstadter

Serguei Y. Kalmykov

Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100-terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, 3D particle-in-cell modelling are examined. First, the Cartesian code VORPAL (Nieter, C. and Cary, J. R. 2004 VORPAL: a versatile plasma simulation code. J. Comput. Phys. 196, 538) using a perfect-dispersion electromagnetic solver precisely describes the laser pulse and bubble dynamics, taking advantage of coarser resolution in the propagation direction, with a proportionally larger time step. Using third-order splines for macroparticles helps suppress the sampling noise while ...

The Reasonable Effectiveness Of Mathematics In The Natural Sciences, Nicolas Fillion

Electronic Thesis and Dissertation Repository

One of the most unsettling problems in the history of philosophy examines how mathematics can be used to adequately represent the world. An influential thesis, stated by Eugene Wigner in his paper entitled "The Unreasonable Effectiveness of Mathematics in the Natural Sciences," claims that "the miracle of the appropriateness of the language of mathematics for the formulation of the laws of physics is a wonderful gift which we neither understand nor deserve." Contrary to this view, this thesis delineates and implements a strategy to show that the applicability of mathematics is very reasonable indeed.

I distinguish three forms of the ...

Computational Investigation Of Steady Navier-Stokes Flows Past A Circular Obstacle In Two--Dimensional Unbounded Domain, Carl Fredrik Jonathan Gustafsson

Open Access Dissertations and Theses

This thesis is a numerical investigation of two-dimensional steady flows past a circular obstacle. In the fluid dynamics research there are few computational results concerning the structure of the steady wake flows at Reynolds numbers larger than 100, and the state-of-the-art results go back to the work of Fornberg (1980) Fornberg (1985). The radial velocity component approaches its asymptotic value relatively slowly if the solution is ``physically reasonable''. This presents a difficulty when using the standard approach such as domain truncation. To get around this problem, in the present research we will develop a spectral technique for the solution of ...

Hard And Soft Error Resilience For One-Sided Dense Linear Algebra Algorithms, Peng Du

Doctoral Dissertations

Dense matrix factorizations, such as LU, Cholesky and QR, are widely used by scientific applications that require solving systems of linear equations, eigenvalues and linear least squares problems. Such computations are normally carried out on supercomputers, whose ever-growing scale induces a fast decline of the Mean Time To Failure (MTTF). This dissertation develops fault tolerance algorithms for one-sided dense matrix factorizations, which handles Both hard and soft errors.

For hard errors, we propose methods based on diskless checkpointing and Algorithm Based Fault Tolerance (ABFT) to provide full matrix protection, including the left and right factor that are normally seen in ...

Approximate Methods For Dynamic Portfolio Allocation Under Transaction Costs, Nabeel Butt

Electronic Thesis and Dissertation Repository

The thesis provides robust and efficient lattice based algorithms for solving dynamic portfolio allocation problems under transaction costs. The early part of the thesis concentrates upon developing a toolbox based on multinomial trees. The multinomial trees are shown to provide a reasonable approximation for most popular transaction cost models in the academic literature. The tool, once forged, is implemented in the powerful Mathematica based parallel computing environment. In the second part of the thesis we provide applications of our framework to real world problems. We show re-balancing portfolios is more valuable in an investment environment where the growth and volatility ...

The Octonions And The Exceptional Lie Algebra G2, Ian M. Anderson

Research Applications

The octonions O are an 8-dimensional non-commutative, non-associative normed real algebra. The set of all derivations of O form a real Lie algebra. It is remarkable fact, first proved by E. Cartan in 1908, that the the derivation algebra of O is the compact form of the exceptional Lie algebra G2. In this worksheet we shall verify this result of Cartan and also show that the derivation algebra of the split octonions is the split real form of G2.

PDF and Maple worksheets can be downloaded from the links below.

Modeling And Mathematical Analysis Of Plant Models In Ecology, Eric A. Eager

Dissertations, Theses, and Student Research Papers in Mathematics

Population dynamics tries to explain in a simple mechanistic way the variations of the size and structure of biological populations. In this dissertation we use mathematical modeling and analysis to study the various aspects of the dynamics of plant populations and their seed banks.

In Chapter 2 we investigate the impact of structural model uncertainty by considering different nonlinear recruitment functions in an integral projection model for Cirsium canescens. We show that, while having identical equilibrium populations, these two models can elicit drastically different transient dynamics. We then derive a formula for the sensitivity of the equilibrium population to changes ...