Numerical Analysis and Computation Commons^™

Impact Of Spallation And Internal Radiation On Fibrous Ablative Materials, Raghava Sai Chaitanya Davuluri

Theses and Dissertations--Mechanical Engineering

Space vehicles are equipped with Thermal Protection Systems (TPS) that encounter high heat rates and protect the payload while entering a planetary atmosphere. For most missions that interest NASA, ablative materials are used as TPS. These materials undergo several mass and energy transfer mechanisms to absorb intense heat. The size and construction of the TPS are based on the composition of the planetary atmosphere and the impact of various ablative mechanisms on the flow field and the material. Therefore, it is essential to quantify the rates of different ablative phenomena to model TPS accurately. In this work, the impact of …

Numerical Investigation On The Effect Of Spectral Radiative Heat Transfer Within An Ablative Material, Raghava S. C. Davuluri, Rui Fu, Kaveh A. Tagavi, Alexandre Martin

Mechanical Engineering Faculty Publications

The spectral radiative heat flux could impact the material response. In order to evaluate it, a coupling scheme between KATS - MR and P1 approximation model of radiation transfer equation (RTE) is constructed and used. A Band model is developed that divides the spectral domain into small bands of unequal widths. Two verification studies are conducted: one by comparing the simulation computed by the Band model with pure conduction results and the other by comparing with similar models of RTE. The comparative results from the verification studies indicate that the Band model is computationally efficient and can be used to …

Fully Coupled Internal Radiative Heat Transfer For The 3d Material Response Of Heat Shield, Raghava S. C. Davuluri, Rui Fu, Kaveh A. Tagavi, Alexandre Martin

Mechanical Engineering Faculty Publications

The radiative transfer equation (RTE) is strongly coupled to the material response code KATS. A P-1 approximation model of RTE is used to account for radiation heat transfer within the material. First, the verification of the RTE model is performed by comparing the numerical and analytical solutions. Next, the coupling scheme is validated by comparing the temperature profiles of pure conduction and conduction coupled with radiative emission. The validation study is conducted on Marschall et al. cases (radiant heating, arc-jet heating, and space shuttle entry), 3D Block, 2D IsoQ sample, and Stardust Return Capsule. The validation results agree well for …

Numerical Reconstruction Of Spalled Particle Trajectories In An Arc-Jet Environment, Raghava S. C. Davuluri, Sean C. C. Bailey, Kaveh A. Tagavi, Alexandre Martin

Mechanical Engineering Faculty Publications

To evaluate the effects of spallation on ablative material, it is necessary to evaluate the mass loss. To do so, a Lagrangian particle trajectory code is used to reconstruct trajectories that match the experimental data for all kinematic parameters. The results from spallation experiments conducted at the NASA HYMETS facility over a wedge sample were used. A data-driven adaptive methodology was used to adapts the ejection parameters until the numerical trajectory matches the experimental data. The preliminary reconstruction results show that the size of the particles seemed to be correlated with the location of the ejection event. The size of …

Orthogonal Recurrent Neural Networks And Batch Normalization In Deep Neural Networks, Kyle Eric Helfrich

Theses and Dissertations--Mathematics

Despite the recent success of various machine learning techniques, there are still numerous obstacles that must be overcome. One obstacle is known as the vanishing/exploding gradient problem. This problem refers to gradients that either become zero or unbounded. This is a well known problem that commonly occurs in Recurrent Neural Networks (RNNs). In this work we describe how this problem can be mitigated, establish three different architectures that are designed to avoid this issue, and derive update schemes for each architecture. Another portion of this work focuses on the often used technique of batch normalization. Although found to be successful …

Unitary And Symmetric Structure In Deep Neural Networks, Kehelwala Dewage Gayan Maduranga

Theses and Dissertations--Mathematics

Recurrent neural networks (RNNs) have been successfully used on a wide range of sequential data problems. A well-known difficulty in using RNNs is the vanishing or exploding gradient problem. Recently, there have been several different RNN architectures that try to mitigate this issue by maintaining an orthogonal or unitary recurrent weight matrix. One such architecture is the scaled Cayley orthogonal recurrent neural network (scoRNN), which parameterizes the orthogonal recurrent weight matrix through a scaled Cayley transform. This parametrization contains a diagonal scaling matrix consisting of positive or negative one entries that can not be optimized by gradient descent. Thus the …

Effects Of Aperiodicity And Frustration On The Magnetic Properties Of Artificial Quasicrystals, Barry Farmer

Theses and Dissertations--Physics and Astronomy

Quasicrystals have been shown to exhibit physical properties that are dramatically different from their periodic counterparts. A limited number of magnetic quasicrystals have been fabricated and measured, and they do not exhibit long-range magnetic order, which is in direct conflict with simulations that indicate such a state should be accessible. This dissertation adopts a metamaterials approach in which artificial quasicrystals are fabricated and studied with the specific goal of identifying how aperiodicity affects magnetic long-range order. Electron beam lithography techniques were used to pattern magnetic thin films into two types of aperiodic tilings, the Penrose P2, and Ammann-Beenker tilings. SQUID …

A Detection And Data Acquisition System For Precision Beta Decay Spectroscopy, Aaron P. Jezghani

Theses and Dissertations--Physics and Astronomy

Free neutron and nuclear beta decay spectroscopy serves as a robust laboratory for investigations of the Standard Model of Particle Physics. Observables such as decay product angular correlations and energy spectra overconstrain the Standard Model and serve as a sensitive probe for Beyond the Standard Model physics. Improved measurement of these quantities is necessary to complement the TeV scale physics being conducted at the Large Hadron Collider. The UCNB, ⁴⁵Ca, and Nab experiments aim to improve upon existing measurements of free neutron decay angular correlations and set new limits in the search for exotic couplings in beta decay. To …

On The Role Of Ill-Conditioning: Biharmonic Eigenvalue Problem And Multigrid Algorithms, Kasey Bray

Theses and Dissertations--Mathematics

Very fine discretizations of differential operators often lead to large, sparse matrices A, where the condition number of A is large. Such ill-conditioning has well known effects on both solving linear systems and eigenvalue computations, and, in general, computing solutions with relative accuracy independent of the condition number is highly desirable. This dissertation is divided into two parts.

In the first part, we discuss a method of preconditioning, developed by Ye, which allows solutions of Ax=b to be computed accurately. This, in turn, allows for accurate eigenvalue computations. We then use this method to develop discretizations that yield accurate computations …

High-Order Integral Equation Methods For Quasi-Magnetostatic And Corrosion-Related Field Analysis With Maritime Applications, Robert Pfeiffer

Theses and Dissertations--Electrical and Computer Engineering

This dissertation presents techniques for high-order simulation of electromagnetic fields, particularly for problems involving ships with ferromagnetic hulls and active corrosion-protection systems.

A set of numerically constrained hexahedral basis functions for volume integral equation discretization is presented in a method-of-moments context. Test simulations demonstrate the accuracy achievable with these functions as well as the improvement brought about in system conditioning when compared to other basis sets.

A general method for converting between a locally-corrected Nyström discretization of an integral equation and a method-of-moments discretization is presented next. Several problems involving conducting and magnetic-conducting materials are solved to verify the accuracy …

Investigation Of Volatile Organic Compounds (Vocs) Detected At Vapor Intrusion Sites, Mohammadyousef Roghani

Theses and Dissertations--Civil Engineering

This dissertation investigates unexplained vapor intrusion field data sets that have been observed at hazardous waste sites, including: 1) non-linear soil gas concentration trends between the VOC source (i.e. contaminated groundwater plume) and the ground surface; and, 2) alternative pathways that serve as entry points for vapors to infiltrate into buildings and serve to increase VOC exposure risks as compared to the classic vapor intrusion model, which primarily considered foundation cracks as the route for vapor entry. The overall hypothesis of this research is that theoretical knowledge of fate and transport processes can be systematically applied to vapor intrusion field …

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

Theses and Dissertations--Mechanical 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 …

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.

A Model For Spheroid Versus Monolayer Response Of Sk-N-Sh Neuroblastoma Cells To Treatment With 15-Deoxy-Pgj2, Dorothy I. Wallace, Ann Dunham, Paula X. Chen, Michelle Chen, Milan Huynh, Evan Rheingold, Olivia F. Prosper

Mathematics Faculty Publications

Researchers have observed that response of tumor cells to treatment varies depending on whether the cells are grown in monolayer, as in vitro spheroids or in vivo. This study uses data from the literature on monolayer treatment of SK-N-SH neuroblastoma cells with 15-deoxy-PGJ₂ and couples it with data on growth rates for untreated SK-N-SH neuroblastoma cells grown as multicellular spheroids. A linear model is constructed for untreated and treated monolayer data sets, which is tuned to growth, death, and cell cycle data for the monolayer case for both control and treatment with 15-deoxy-PGJ₂. The monolayer …

The Krylov Subspace Methods For The Computation Of Matrix Exponentials, Hao Wang

Theses and Dissertations--Mathematics

The problem of computing the matrix exponential e^tA arises in many theoretical and practical problems. Many methods have been developed to accurately and efficiently compute this matrix function or its product with a vector, i.e., e^tAv. In the past few decades, with the increasing need of the computation for large sparse matrices, iterative methods such as the Krylov subspace methods have proved to be a powerful class of methods in dealing with many linear algebra problems. The Krylov subspace methods have been introduced for computing matrix exponentials by Gallopoulos and Saad, and the corresponding error bounds …

Singular Value Computation And Subspace Clustering, Qiao Liang

Theses and Dissertations--Mathematics

In this dissertation we discuss two problems. In the first part, we consider the problem of computing a few extreme eigenvalues of a symmetric definite generalized eigenvalue problem or a few extreme singular values of a large and sparse matrix. The standard method of choice of computing a few extreme eigenvalues of a large symmetric matrix is the Lanczos or the implicitly restarted Lanczos method. These methods usually employ a shift-and-invert transformation to accelerate the speed of convergence, which is not practical for truly large problems. With this in mind, Golub and Ye proposes an inverse-free preconditioned Krylov subspace method, …

A Posteriori Error Estimates For Surface Finite Element Methods, Fernando F. Camacho

Theses and Dissertations--Mathematics

Problems involving the solution of partial differential equations over surfaces appear in many engineering and scientific applications. Some of those applications include crystal growth, fluid mechanics and computer graphics. Many times analytic solutions to such problems are not available. Numerical algorithms, such as Finite Element Methods, are used in practice to find approximate solutions in those cases.

In this work we present L2 and pointwise a posteriori error estimates for Adaptive Surface Finite Elements solving the Laplace-Beltrami equation −△_Γ u = f . The two sources of errors for Surface Finite Elements are a Galerkin error, and a …

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 …