Applied Mathematics Commons^™

Electrodynamical Modeling For Light Transport Simulation, Michael G. Saunders

Undergraduate Honors Theses

Modernity in the computer graphics community is characterized by a burgeoning interest in physically based rendering techniques. That is to say that mathematical reasoning from first principles is widely preferred to ad hoc, approximate reasoning in blind pursuit of photorealism. Thereby, the purpose of our research is to investigate the efficacy of explicit electrodynamical modeling by means of the generalized Jones vector given by Azzam [1] and the generalized Jones matrix given by Ortega-Quijano & Arce-Diego [2] in the context of stochastic light transport simulation for computer graphics. To augment the status quo path tracing framework with such a modeling …

A Coupled Localized Rbf Meshless/Drbem Formulation For Accurate Modeling Of Incompressible Fluid Flows, Leonardo Bueno, Eduardo Divo, Alain J. Kassab

Publications

Velocity-pressure coupling schemes for the solution of incompressible fluid flow problems in Computational Fluid Dynamics (CFD) rely on the formulation of Poisson-like equations through projection methods. The solution of these Poisson-like equations represent the pressure correction and the velocity correction to ensure proper satisfaction of the conservation of mass equation at each step of a time-marching scheme or at each level of an iteration process. Inaccurate solutions of these Poisson-like equations result in meaningless instantaneous or intermediate approximations that do not represent the proper time-accurate behavior of the flow. The fact that these equations must be solved to convergence at …

An Rbf Interpolation Blending Scheme For Effective Shock-Capturing, M. Harris, Eduardo Divo, Alain J. Kassab

Publications

In recent years significant focus has been given to the study of Radial basis functions (RBF), especially in their use on solving partial differential equations (PDE). RBF have an impressive capability of inter- polating scattered data, even when this data presents localized discontinuities. However, for infinitely smooth RBF such as the Multiquadrics, inverse Multiquadrics, and Gaussian, the shape parameter must be chosen properly to obtain accurate approximations while avoiding ill-conditioning of the interpolating matrices. The optimum shape parameter can vary significantly depending on the field, particularly in locations of steep gradients, shocks, or discontinuities. Typically, the shape parameter is chosen …

Code For "Noise-Enhanced Coding In Phasic Neuron Spike Trains", Cheng Ly, Brent D. Doiron

Statistical Sciences and Operations Research Data

This zip file contains Matlab scripts and ode (XPP) files to calculate the statistics of the models in "Noise-Enhanced Coding in Phasic Neuron Spike Trains". This article is published in PLoS ONE.

Multiplicative Noise Removal With A Sparsity-Aware Optimization Model, Jian Lu, Lixin Shen, Chen Xu, Yuesheng Xu

Mathematics & Statistics Faculty Publications

Restoration of images contaminated by multiplicative noise (also known as speckle noise) is a key issue in coherent image processing. Notice that images under consideration are often highly compressible in certain suitably chosen transform domains. By exploring this intrinsic feature embedded in images, this paper introduces a variational restoration model for multiplicative noise reduction that consists of a term reflecting the observed image and multiplicative noise, a quadratic term measuring the closeness of the underlying image in a transform domain to a sparse vector, and a sparse regularizer for removing multiplicative noise. Being different from popular existing models which focus …

Orbital Stability Results For Soliton Solutions To Nonlinear Schrödinger Equations With External Potentials, Joseph B. Lindgren

Theses and Dissertations--Mathematics

For certain nonlinear Schroedinger equations there exist solutions which are called solitary waves. Addition of a potential $V$ changes the dynamics, but for small enough $||V||_{L^\infty}$ we can still obtain stability (and approximately Newtonian motion of the solitary wave's center of mass) for soliton-like solutions up to a finite time that depends on the size and scale of the potential $V$. Our method is an adaptation of the well-known Lyapunov method.

For the sake of completeness, we also prove long-time stability of traveling solitons in the case $V=0$.

A Physics-Based Approach To Modeling Wildland Fire Spread Through Porous Fuel Beds, Tingting Tang

Theses and Dissertations--Mechanical and Aerospace Engineering

Wildfires are becoming increasingly erratic nowadays at least in part because of climate change. CFD (computational fluid dynamics)-based models with the potential of simulating extreme behaviors are gaining increasing attention as a means to predict such behavior in order to aid firefighting efforts. This dissertation describes a wildfire model based on the current understanding of wildfire physics. The model includes physics of turbulence, inhomogeneous porous fuel beds, heat release, ignition, and firebrands. A discrete dynamical system for flow in porous media is derived and incorporated into the subgrid-scale model for synthetic-velocity large-eddy simulation (LES), and a general porosity-permeability model is …

Fire, Ice, Water, And Dirt: A Simple Climate Model, John Kroll

Mathematics & Statistics Faculty Publications

A simple paleoclimate model was developed as a modeling exercise. The model is a lumped parameter system consisting of an ocean (water), land (dirt), glacier, and sea ice (ice) and driven by the sun (fire). In comparison with other such models, its uniqueness lies in its relative simplicity yet yielding good results. For nominal values of parameters, the system is very sensitive to small changes in the parameters, yielding equilibrium, steady oscillations, and catastrophes such as freezing or boiling oceans. However, stable solutions can be found, especially naturally oscillating solutions. For nominally realistic conditions, natural periods of order 100kyrs are …

Complex Solutions Of The Time Fractional Gross-Pitaevskii (Gp) Equation With External Potential By Using A Reliable Method, Nasir Taghizadeh, Mona N. Foumani

Applications and Applied Mathematics: An International Journal (AAM)

In this article, modified (G'/G )-expansion method is presented to establish the exact complex solutions of the time fractional Gross-Pitaevskii (GP) equation in the sense of the conformable fractional derivative. This method is an effective method in finding exact traveling wave solutions of nonlinear evolution equations (NLEEs) in mathematical physics. The present approach has the potential to be applied to other nonlinear fractional differential equations. Based on two transformations, fractional GP equation can be converted into nonlinear ordinary differential equation of integer orders. In the end, we will discuss the solutions of the fractional GP equation with external potentials.

Impact Of Permeable Lining Of The Wall On The Peristaltic Flow Of Herschel Bulkley Fluid, G. C. Sankad, Asha Patil

Applications and Applied Mathematics: An International Journal (AAM)

The peristaltic motion is modeled for the Herschel Bulkley fluid, considered to flow in a non-uniform inclined channel. The channel wall is supposed to be lined with a non-erodible porous material. The flow is considered to be moving in a wave frame of reference moving with same velocity as of the sinusoidal wave. Low Reynolds number and long wave length assumptions are made to solve the model. Analytical solution is obtained for the pressure difference and also for the frictional force. Graphs are plotted, using Mathematica software, for both the results of pressure difference and frictional force against time average …

On The Propagation Of Atmospheric Gravity Waves In A Non-Uniform Wind Field: Introducing A Modified Acoustic-Gravity Wave Equation, Ahmad Talaei

All Graduate Plan B and other Reports, Spring 1920 to Spring 2023

Atmospheric gravity waves play fundamental roles in a broad-range of dynamical processes extending throughout the Earth’s neutral atmosphere and ionosphere. In this paper, we present a modified form for the acoustic-gravity wave equation and its dispersion relationships for a compressible and non-stationary atmosphere in hydrostatic balance. Importantly, the solutions have been achieved without the use of the well-known Boussinesq approximation which have been used extensively in previous studies.

We utilize the complete set of governing equations for a compressible atmosphere with non-uniform airflows to determine an equation for vertical velocity of possible atmospheric waves. This intricate wave equation is simplified …

Making Waves: A Computational Swimming Lamprey With Sensory Feedback, Christina Hamlet, Eric D. Tytell, Lisa J. Fauci, Kathleen A. Hoffman

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Interdisciplinary Undergraduate Research In Biofluids, Eva M. Strawbridge

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

A Fast Model For The Simulation Of External Gear Pumps, Zechao Lu, Xinran Zhao, Andrea Vacca

The Summer Undergraduate Research Fellowship (SURF) Symposium

External gear pump is an important category of positive displacement fluid machines used to perform the mechanical–hydraulic energy conversions in many fluid power applications. An efficient numerical simulation program is needed to simulate the system in order to provide a direction for design purpose. The model consists of a lumped parameter fluid dynamic model and a model that simulates the radial micro-motions of the gear’s axes of rotation. The system consists of a set of ordinary differential equations related to the conservation on mass of the internal control volumes of the pump, which are given by the tooth space volumes …

Model For Computing Kinetics Of The Graphene Edge Epitaxial Growth On Copper, Mikhail Khenner

Mathematics Faculty Publications

A basic kinetic model that incorporates a coupled dynamics of the carbon atoms and dimers ona copper surface is used to compute growth of a single-layer graphene island. The speed of theisland's edge advancement on Cu[111] and Cu[100] surfaces is computed as a function of the growthtemperature and pressure. Spatially resolved concentration pro les of the atoms and dimers aredetermined, and the contributions provided by these species to the growth speed are discussed.Island growth in the conditions of a thermal cycling is studied.

Approximate Analytical Solution Of Boussinesq Equation In Homogeneous Medium With Leaky Base, Rajeev K. Bansal

Applications and Applied Mathematics: An International Journal (AAM)

Approximate analytical solutions of Boussinesq equation are widely used for approximation of subsurface seepage flow in confined and unconfined aquifers under varying hydrological conditions. In this paper, we use a 2-dimensional linearized Boussinesq equation to simulate the water table fluctuations in an isotropic aquifer overlying a semi pervious bed under multiple localized recharge and withdrawal. The unconfined aquifer is considered to be in contact with two water bodies of constant water head along opposite cost lines, while the remaining two faces have no flow condition. The mathematical model is solved analytically using finite Fourier sine transform and the application of …

On The Global Existence And Boundedness Of Solutions Of Nonlinear Vector Differential Equations Of Third Order, Timur Ayhan, Cemil Tunç

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we give some criteria to ensure the global existence and boundedness of solutions to a kind of third order nonlinear vector differential equations. By using the Lyapunov's direct method, we obtain a new result on the topic and give an example for the illustrations. Our result includes, completes and improves some earlier results in the literature.

Non-Newtonian Prandtl Fluid Over Stretching Permeable Surface, N. R. Jain, M. G. Timol

Applications and Applied Mathematics: An International Journal (AAM)

An analysis is made of the velocity and temperature distribution in the flow of a viscous incompressible fluid caused by the stretching permeable surface which issues in the Prandtl fluid. Parandtl fluid is a pseudoplastic visco-inelastic non-Newtonian fluid. The governing partial differential equations are reduced to ordinary differential equations using deductive group transformation and similarity solution is derived. Numerical solutions to the reduced non-linear similarity equations are then obtained by adopting shooting method using the Nachtsheim-Swigert iteration technique. The results of the numerical solution are then presented graphically in the form of velocity and temperature profiles. The corresponding skin friction …

Nonlinear Harmonic Modes Of Steel Strings On An Electric Guitar, Joel Wenrich

Senior Theses

Steel strings used on electric and acoustic guitars are non-ideal oscillators that can produce imperfect intonation. According to theory, this intonation should be a function of the bending stiffness of the string, which is related to the dimensions of length and thickness of the string. To test this theory, solid steel strings of three different linear densities were analyzed using an oscilloscope and a Fast Fourier Transform function. We found that strings exhibited more drastic nonlinear harmonic behavior as their effective length was shortened and as linear density increased.

Spontaneous Parametric Down Conversion Of Photons Through Β-Barium Borate, Luke Horowitz

Physics

An apparatus for detecting pairs of entangled 405nm photons that have undergone Spontaneous Parametric Down Conversion through β-Barium Borate is described. By using avalanche photo-diodes to detect the low-intensity converted beam and a coincidence module to register coincident photons, it is possible to create an apparatus than can be used to perform quantum information experiments under a budget appropriate for an undergraduate physics lab.

Theoretical Modeling Of The Effect Of Polymer Chain Immobilization Rates On Holographic Recording In Photopolymers, Dana Mackey, Paul O'Reilly, Izabela Naydenova

Articles

This paper introduces an improved mathematical model for holographic grating formation in an acrylamide-based photopolymer, which consists of partial differential equations derived from physical laws. The model is based on the two-way diffusion theory of \cite{izabela}, which assumes short polymer chains are free to diffuse, and generalizes a similar model presented in \cite{josab} by introducing an immobilization rate governed by chain growth and cross-linking. Numerical simulations were carried out in order to investigate the behaviour of the photopolymer system for short and long exposures, with particular emphasis on the effect of recording parameters (such as illumination frequency and intensity), as …

Early Afternoon Concurrent Panel Sessions: Commercial Space Industry Snapshot: Presentation: Small Catapult-Assisted Horizontal-Launch Reusable Rbcc Ssto Spaceplane For Economical Short-Duration Leo Access, Nihad E. Daidzic,

Aviation / Aeronautics / Aerospace International Research Conference

This article discusses the conceptual design, flight trajectory calculations, and utilization of the possible future horizontally-launched reusable Single-Stage-to-Orbit (SSTO) spaceplane for small payload short-duration manned/unmanned access to Low-Earth-Orbit (LEO). The 10,000 lb spaceplane would use 5,000 ft catapult-assist horizontal-launch facility and conduct powered approach and landing on conventional horizontal paved runways following the gliding atmospheric re-entry. To increase the economy of operation, the launch facility located at high elevations (4,000+ ft) equatorial region is required, such as, the plateaus in Kenya and Tanzania in Africa and/or Ecuador in South America. A 500-lb payload, including pilot-commander, is envisioned. The propulsion cycle …

Effect Of Spalled Particles Thermal Degradation On A Hypersonic Flow Field Environment, Raghava S. C. Davuluri, Huaibao Zhang, Alexandre Martin

Mechanical Engineering Faculty Publications

Two-way coupling is performed between a spallation code and a hypersonic aerothermodynamics CFD solver to evaluate the effect of spalled particles on the flow field. Time accurate solutions are computed in argon and air flow fields. A single particle simulations and multiple particles simulations are performed and studied. The results show that the carbon vapor released by spalled particles tend to change the composition of the flow field, particularly the upstream region of the shock.

The Global Stability Of The Solution To The Morse Potential In A Catastrophic Regime, Weerapat Pittayakanchit

HMC Senior Theses

Swarms of animals exhibit aggregations whose behavior is a challenge for mathematicians to understand. We analyze this behavior numerically and analytically by using the pairwise interaction model known as the Morse potential. Our goal is to prove the global stability of the candidate local minimizer in 1D found in A Primer of Swarm Equilibria. Using the calculus of variations and eigenvalues analysis, we conclude that the candidate local minimizer is a global minimum with respect to all solution smaller than its support. In addition, we manage to extend the global stability condition to any solutions whose support has a single …

The Dynamics Of Flat Surface Internal Geophysical Waves With Currents, Alan Compelli, Rossen Ivanov

Articles

A two-dimensional water wave system is examined consisting of two discrete incompressible fluid domains separated by a free common interface. In a geophysical context this is a model of an internal wave, formed at a pycnocline or thermocline in the ocean. The system is considered as being bounded at the bottom and top by a flatbed and wave-free surface respectively. A current profile with depth-dependent currents in each domain is considered. The Hamiltonian of the system is determined and expressed in terms of canonical wave-related variables. Limiting behavior is examined and compared to that of other known models. The linearised …

Stereographic Visualization Of Bose-Einstein Condensate Clouds To Measure The Gravitational Constant, Ed Wesley Wells

College of Graduate Studies: Theses & Dissertations

This thesis describes a set of tools that can be used for the rapid design of atom interferometer schemes suitable for measuring Newton's Universal Gravitation constant also known as "Big G". This tool set is especially applicable to Bose--Einstein--condensed systems present in NASA's Cold Atom Laboratory experiment to be deployed to the International Space Station in 2017. These tools include a method of approximating the solutions of the nonlinear Schrödinger or Gross--Pitaevskii equation (GPE) using the Lagrangian Variational Method. They also include a set of software tools for translating the approximate solutions of the GPE into images of the optical …

Simulation Of Nuclear Fusion Using A One Dimensional Particle In Cell Method, Steven T. Margell

Cal Poly Humboldt theses and projects

In this thesis several novel techniques are developed to simulate fusion events in an isotropic, electrostatic three-dimensional Deuterium-Tritium plasma. These techniques allow us to accurately predict three-dimensional collision events with a one-dimensional model while simultaneously reducing compute time via a nearest neighbor algorithm. Furthermore, a fusion model based on first principles is developed that yields an average fusion reactivity which correlates well with empirical results.

Evaluation Of Ray-Path Integrals In Geometrical Optics, John A. Adam, Michael Pohrivchak

Mathematics & Statistics Faculty Publications

A brief summary of the physical context to this paper is provided, and the deviation angle undergone by an incident ray after k internal reflections inside a transparent unit sphere is formulated. For radially inhomogeneous spheres (in particular) this angle is related to a ray-path integral; an improper integral for which there are relatively few known exact analytical forms, even for simple refractive index profiles n(r). Thus for a linear profile the integral is a combination of incomplete elliptic integrals of the first and third kinds (though not all are as complicated as this). The ray-path integral is evaluated …

Approximate Solutions For The Flow And Heat Transfer Due To A Stretching Sheet Embedded In A Porous Medium With Variable Thickness, Variable Thermal Conductivity And Thermal Radiation Using Laguerre Collocation Method, M. M. Khader, Ahmed M. Megahed

Applications and Applied Mathematics: An International Journal (AAM)

In this article, a numerical approach is given for studying the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with a power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by a non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing PDEs are transformed into a system of coupled non-linear ODEs which are using appropriate boundary conditions for various physical parameters. The proposed method is based on replacement of the unknown function by truncated series …

Exact Implicit Solution Of Nonlinear Heat Transfer In Rectangular Straight Fin Using Symmetry Reduction Methods, M. S. Abdel Latif, A. H. Abdel Kader, H. M. Nour

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, the exact implicit solution of the second order nonlinear ordinary differential equation which governing heat transfer in rectangular fin is obtained using symmetry reduction methods. General relationship among the temperature at the fin tip, the temperature gradient at the fin base, the mode of heat transfer, 𝑛 and the fin parameters 𝑁 and ℰ is obtained. Some numerical examples are discussed and it is shown that the temperature of fin increases when approaching from the heat source. The relationship between the fin efficiency and the temperature of fin tip is obtained for any value of the mode …