Numerical Analysis and Computation Commons^™

Numerical Modeling Of Thermo-Poroelasticity Using Finite Element Method, Maya Mckean

Discovery Day - Daytona Beach

We propose a numerical method for solving and modeling thermo-poroelasticity problems using a finite element formulation. Thermo-poroelasticity models describe the coupled interaction between mechanical deformation, fluid flow, and heat transfer in specific materials or environments over time. These models demonstrate the evolution of displacement, pressure, and temperature; we compute these fields in this work using backward Euler time discretization and Enriched Galerkin finite element spatial discretization. For these computations we used FreeFEM, a partial differential equation solver that uses the finite element method, which produced our numerical results. We then compared these values with the expected analytical solution. This was …

Numerical Method For Strongly Variable-Density Flows At Low Mach Number: Flame-Sheet Regularisation And A Mass-Flux Immersed Boundary Method, Matheus P. Severino, Fernando F. Fachini, Elmer M. Gennaro, Daniel Rodríguez, Leandro F. Souza

Mathematical Modelling and Numerical Simulation with Applications

A low-Mach-number flow, in the laminar regime, has intrinsically two characteristic spatial scales for a given time scale, or two characteristic temporal scales for a given spatial scale, and these dual scales are very different due to the disparity between the flow and acoustic speed. Therefore, low-Mach-number flows impose mathematical and computational challenges in their description. Standard numerical methods for compressible flows, which are typically designed for problems with a single dominant spatial and temporal scale, require alternative approaches, such as preconditioning techniques or solvers tailored for low-Mach-number equations. The present work introduces a simplified fluid dynamics model for flows …

Intermediate-Scale Outflow Dynamics Of Eta Carinae, Edmund J. Garcia, Matthew C. Fleenor

Departmental Honors & Graduate Capstone Projects

η Carinae (η Car) is a binary system, with the larger star being an extremely massive, luminous blue variable (LBV) beyond the Eddington Limit. Surrounding the η Car system, numerous multi-wavelength imaging campaigns reveal axisymmetric structures with the expanding bipolar Homunculus Nebula (¡1 pc). In combination with the episodic eruptive history of the η Car system, our initial intermediate-scale imaging revealed further axisym- metric structures (1-5 pc). To gain a more expansive view of how the small scale structure connects to panoramic imaging of the η Car region, we constructed a deep, optical, narrowband mosaic of 189 images utilizing the …

Microgravity-Induced Alterations In Left Atrial Hemodynamics And Thrombogenic Risk: Insights From Healthy And Atrial Fibrillation Models, Grace M. Hoeppner

Dissertations, Master's Theses and Master's Reports

Background: Microgravity exposure alters cardiovascular loading, yet its impact on left atrial flow dynamics and thrombotic risk remains poorly understood. This study investigates how spaceflight-relevant microgravity-induced changes in cardiac outflow affect left atrial hemodynamics in healthy individuals and patients with atrial fibrillation.

Methods: Patient-specific left atrial models were generated for three healthy individuals and three AF patients. Computational fluid dynamics (CFD) simulations were performed using each patient’s baseline mitral outflow waveform and two modified waveforms representing short- and long-duration post-flight cardiac loading changes derived from echocardiographic observations. Hemodynamic metrics included left atrial velocity, time averaged wall shear stress, oscillatory shear …

Coupled Machine Learning Models: Combining Observations And Numerical Analysis In A Physics-Regularized Approach, Austin B. Schmidt

LSU New Orleans Theses and Dissertations

This dissertation investigates surrogate modeling for fixed-location environmental forecasting using novel data-combination techniques. The work surveys the landscape of observational measurements and numerically generated data, identifying similar research and gaps in current methodologies. The ratio-coupled training framework is introduced to combine two data sources per predicted feature through a tunable parameter that weights training signal strength. An optimization scheme is developed to simultaneously tune surrogate weights and the coupled signal ratio, allowing relative influence between signals to act as an explicit regularizer. Three case studies demonstrate the methodology and approach in a variety of contexts. The first study is based …

Acoustic-Gravity Wave Propagation Based On Solutions To The Generalized Multi-Component Transport Equations, Benedict PiñEyro

Doctoral Dissertations and Master's Theses

Nonlinear atmospheric models have provided important insight into acoustic waves generated by natural and man-made hazards, which may steepen into shocks or N-waves while also dissipating when propagating in the thermosphere. Although models have yielded results that agree with observations of ionospheric perturbations, dynamical models for the diffusive and stratified lower thermosphere often use single gas approximations with height-dependent physical properties that omit the dynamics of the major and minor constituents. Thus, the inter-species diffusion associated with these flows (e.g. variations of mean molecular weight, and specific heat) are not accounted for. This approximation is simpler and less computationally expensive …

Eulerian Smoke Simulation With Multiple Fields, Diyang Zhang

Dartmouth College Master’s Theses

Fluid simulation is a cornerstone of computer graphics, enabling the realistic depiction of dynamic phenomena such as smoke, fire, and other gaseous behaviours. This thesis focuses on advancing Eulerian smoke simulation techniques, with a particular emphasis on grid-based simulations that capture intricate vortical structures and fine visual details.

We propose several detail-preserving frameworks that incorporate various scalar and vector fields within the simulation pipeline, including velocity, impulse, and Lamb vectors, along with their decompositions and transformed representations. By mathematically analyzing the properties of impulse, we derive its scalar fields decomposition (ImpSFD), which introduces an alternative numerical interpretation, and Vortex-Particles in …

Cfd Analysis Of Hydrodynamic Cavitation Through An Orifice: Influence Of Different Inlet Pressures And Number Of Orifice Holes, Lemthong Chanphavong, Vongsavanh Chanthaboune, Keophousone Phonhalath

ASEAN Journal on Science and Technology for Development

Hydrodynamic cavitation (HC) is considered an energy-efficient process with high potential for utilization in many chemical processes. This study presents a computational fluid dynamics (CFD) analysis of cavitating flow through an orifice with a constant flow area. The Reynolds-Averaged Navier-Stokes (RANS) equations, coupled with turbulence and cavitation models, are employed to capture the complex flow behaviors. The effects of inlet pressures and number of orifice-holes on cavitation behavior are investigated. Result of the numerical simulation is validated with the existing experimental data from the literature. The CFD study revealed that cavitation initiates just behind the inlet edge of the orifice …

Les-C Turbulence Models And Fluid Flow Modeling: Analysis And Application To Incompressible Turbulence And Fluid-Fluid Interaction, Kyle J. Schwiebert

Dissertations, Master's Theses and Master's Reports

In the first chapter of this dissertation, we give some background on the Navier-Stokes equations and turbulence modeling. The next two chapters in this dissertation focus on two important numerical difficulties arising in fluid flow modeling: poor mass-conservation and nonphysical oscillations. We investigate two different formulations of the Crank-Nicolson method for the Navier-Stokes equations. The most attractive implementation, second order accurate for both velocity and pressure, is shown to introduce non-physical oscillations. We then propose two options which are shown to avoid the poor behavior. Next, we show that grad-div stabilization, previously assumed to have no effect on the target …

Controlled Manipulation And Transport By Microswimmers In Stokes Flows, Jake Buzhardt

All Dissertations

Remotely actuated microscale swimming robots have the potential to revolutionize many aspects of biomedicine. However, for the longterm goals of this field of research to be achievable, it is necessary to develop modelling, simulation, and control strategies which effectively and efficiently account for not only the motion of individual swimmers, but also the complex interactions of such swimmers with their environment including other nearby swimmers, boundaries, other cargo and passive particles, and the fluid medium itself. The aim of this thesis is to study these problems in simulation from the perspective of controls and dynamical systems, with a particular focus …

Modeling Single And Multiple Pacemaker Interaction In Jellyfish Locomotion, Alexander Hoover

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Boundary Integral Equation Methods For Superhydrophobic Flow And Integrated Photonics, Kosuke Sugita

Dissertations

This dissertation presents fast integral equation methods (FIEMs) for solving two important problems encountered in practical engineering applications.

The first problem involves the mixed boundary value problem in two-dimensional Stokes flow, which appears commonly in computational fluid mechanics. This problem is particularly relevant to the design of microfluidic devices, especially those involving superhydrophobic (SH) flows over surfaces made of composite solid materials with alternating solid portions, grooves, or air pockets, leading to enhanced slip.

The second problem addresses waveguide devices in two dimensions, governed by the Helmholtz equation with Dirichlet conditions imposed on the boundary. This problem serves as a …

Hydrodynamic And Physicochemical Interactions Between An Active Janus Particle And An Inactive Particle, Jessica S. Rosenberg

Dissertations, Theses, and Capstone Projects

Active matter is an area of soft matter science in which units consume energy and turn it into autonomous motion. Groups of these units – whether flocks of birds, bacterial colonies, or even collections of synthetically-made active particles – may exhibit complex behavior on large scales. While the large-scale picture is of great importance, so is the microscopic scale. Studying the individual particles that make up active matter will allow us to understand how they move, and whether and under what circumstances their activity can be controlled.

Here we delve into the world of active matter by studying colloidal-sized (100 …

Mathematical Modeling Of Egain Droplets Sliding Down An Inclined Plane, Shawn Koohy, Kathryn Massey, Luis Schneegans, Megan Vezzetti, Keith Hillaire, Souradip Chattopadhyay, Jessie Chen, Carmen Lee, Karen Daniels

Undergraduate Research Symposium

Eutectic Gallium-Indium (EGaIn) is a room-temperature liquid metal alloy that dramatically changes its surface tension and dynamics under an applied electric field. EGaIn has been used heavily in soft electronics engineering due to its high conductivity, malleability, and safety. However, the absence of mathematical modeling in the current literature makes its behavior difficult to understand and predict. In this study, we present a one-dimensional lubrication model for the dynamics of an EGaIn droplet moving along an inclined plane. Our model incorporates essential physical effects and parameters including oxidation, capillary action, diffusion, gravity, and Marangoni effects. In particular, we incorporate effects …

Practical Implementation Of The Immersed Interface Method With Triangular Meshes For 3d Rigid Solids In A Fluid Flow, Norah Hakami

Mathematics Theses and Dissertations

When employing the immersed interface method (IIM) to simulate a fluid flow around a moving rigid object, the immersed object can be replaced by a virtual fluid enclosed by singular forces on the interface between the real and virtual fluids. These forces represent the impact of the rigid motion on the fluid flow and cause jump discontinuities across the interface in the whole flow field. Then, the IIM resolves the fluid flow on a fixed computational domain by directly incorporating the jump conditions across the interface into numerical schemes. Previous development of the method is limited to simple smooth boundaries. …

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

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

An Examination Of Fontan Circulation Using Differential Equation Models And Numerical Methods, Vanessa Maybruck

Honors Student Research

Certain congenital heart defects can lead to the development of only a single pumping chamber, or ventricle, in the heart instead of the usual two ventricles. Individuals with this defect undergo a corrective, three-part surgery, the third step of which is the Fontan procedure, but as the patients age, their cardiovascular health will likely deteriorate. Using computational fluid dynamics and differential equations, Fontan circulation can be modeled to investigate why the procedure fails and how Fontan failure can be maximally prevented. Borrowing from well-established literature on RC circuits, the differential equation models simulate systemic blood flow in a piecewise, switch-like …

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 …

Five-Wave Resonances In Deep Water Gravity Waves: Integrability, Numerical Simulations And Experiments, Dan Lucas, Marc Perlin, Dian-Yong Liu, Shane Walsh, Rossen Ivanov, Miguel D. Bustamante

Articles

In this work we consider the problem of finding the simplest arrangement of resonant deep water gravity waves in one-dimensional propagation, from three perspectives: Theoretical, numerical and experimental. Theoretically this requires using a normal-form Hamiltonian that focuses on 5-wave resonances. The simplest arrangement is based on a triad of wave vectors K₁ + K₂ = K₃ (satisfying specific ratios) along with their negatives, corresponding to a scenario of encountering wave packets, amenable to experiments and numerical simulations. The normal-form equations for these encountering waves in resonance are shown to be non-integrable, but they admit an integrable reduction …

Improving The Temporal Accuracy Of Turbulence Models And Resolving The Implementation Issues Of Fluid Flow Modeling, Kyle J. Schwiebert

Dissertations, Master's Theses and Master's Reports

A sizeable proportion of the work in this thesis focuses on a new turbulence model, dubbed ADC (the approximate deconvolution model with defect correction). The ADC is improved upon using spectral deferred correction, a means of constructing a higher order ODE solver. Since both the ADC and SDC are based on a predictor-corrector approach, SDC is incorporated with essentially no additional computational cost. We will show theoretically and using numerical tests that the new scheme is indeed higher order in time than the original, and that the benefits of defect correction, on which the ADC is based, are preserved.

The …

Accurate Spectral Algorithms For Solving Variable-Order Fractional Percolation Equations, M. A. Abdelkawy

Applications and Applied Mathematics: An International Journal (AAM)

A high accurate spectral algorithm for one-dimensional variable-order fractional percolation equations (VO-FPEs) is considered.We propose a shifted Legendre Gauss-Lobatto collocation (SL-GLC) method in conjunction with shifted Chebyshev Gauss-Radau collocation (SC-GR-C) method to solve the proposed problem. Firstly, the solution and its space fractional derivatives are expanded as shifted Legendre polynomials series. Then, we determine the expansion coefficients by reducing the VO-FPEs and its conditions to a system of ordinary differential equations (SODEs) in time. The numerical approximation of SODEs is achieved by means of the SC-GR-C method. The under-study’s problem subjected to the Dirichlet or non-local boundary conditions is presented …

Numerical Simulations Of Nonlinear Waves And Their Stability: Stokes Waves And Nonlinear Schroedinger Equation, Anastassiya Semenova

Mathematics & Statistics ETDs

The present work offers an investigation of dynamics and stability of nonlinear waves in Hamiltonian systems. The first part of the manuscript discusses the classical problem of water waves on the surface of an ideal fluid in 2D. We demonstrate how to construct the Stokes waves, and how to apply a continuation method to find waves in close vicinity to the limiting Stokes wave. We provide new insight into the stability of the Stokes waves by identifying previously inaccessible branches of instability in the equations of motion for the fluid. We provide numerical evidence that pairs of unstable eigenvalues of …

A Study Of The Design Of Adaptive Camber Winglets, Justin J. Rosescu

Master's Theses

A numerical study was conducted to determine the effect of changing the camber of a winglet on the efficiency of a wing in two distinct flight conditions. Camber was altered via a simple plain flap deflection in the winglet, which produced a constant camber change over the winglet span. Hinge points were located at 20%, 50% and 80% of the chord and the trailing edge was deflected between -5° and +5°. Analysis was performed using a combination of three-dimensional vortex lattice method and two-dimensional panel method to obtain aerodynamic forces for the entire wing, based on different winglet camber configurations. …

Local Non-Similar Solution Of Powell-Eyring Fluid Flow Over A Vertical Flat Plate, Hemangini Shukla, Hema C. Surati, M. G. Timol

Applications and Applied Mathematics: An International Journal (AAM)

Our objective is to obtain the non-similarity solution of non-Newtonian fluid for Powell-Eyring model by a local non-similarity method. Here, free stream velocity is considered in power-law form (𝑈=𝑥^m). The governing equations are transformed using non-similar transformations and derived equations are treated as ordinary differential equations. Non-similar solutions are obtained for different values of power-law index 𝑚 and stream-wise location 𝜉. Influence of various parameters on velocity and temperature field are presented graphically using MATLAB bvp4c solver.

Meshless Modeling Of Flow Dispersion And Progressive Piping In Poroelastic Levees, Anthony Khoury, Eduardo Divo, Alain J. Kassab, Sai Kakuturu, Lakshmi Reddi

Publications

Performance data on earth dams and levees continue to indicate that piping is one of the major causes of failure. Current criteria for prevention of piping in earth dams and levees have remained largely empirical. This paper aims at developing a mechanistic understanding of the conditions necessary to prevent piping and to enhance the likelihood of self-healing of cracks in levees subjected to hydrodynamic loading from astronomical and meteorological (including hurricane storm surge-induced) forces. Systematic experimental investigations are performed to evaluate erosion in finite-length cracks as a result of transient hydrodynamic loading. Here, a novel application of the localized collocation …

Solitary And Periodic Exact Solutions Of The Viscosity-Capillarity Van Der Waals Gas Equations, Emad A. Az-Zo’Bi

Applications and Applied Mathematics: An International Journal (AAM)

Periodic and soliton solutions are derived for the (1+1)-dimensional van der Waals gas system in the viscosity-capillarity regularization form. The system is handled via the e-^φ(ξ) -expansion method. The obtained solutions have been articulated by the hyperbolic, trigonometric, exponential and rational functions with arbitrary constants. Mathematical analysis and numerical graphs are provided for some solitons, periodic and kink solitary wave solutions to visualize the dynamics of equations. Obtained results reveal that the method is very influential and effective tool for solving nonlinear partial differential equations in applied mathematics.

Numerical Solution Of 3rd Order Ode Using Fdm: On A Moving Surface In Mhd Flow Of Sisko Fluid, Manisha Patel, Jayshri Patel, M. G. Timol

Applications and Applied Mathematics: An International Journal (AAM)

A Similarity group theoretical technique is used to transform the governing nonlinear partial differential equations of two dimensional MHD boundary layer flow of Sisko fluid into nonlinear ordinary differential equations. Then the resulting third order nonlinear ordinary differential equation with corresponding boundary conditions is linearised by Quasi linearization method. Numerical solution of the linearised third order ODE is obtained using Finite Difference method (FDM). Graphical presentation of the solution is given.

Surface Waves Over Currents And Uneven Bottom, Alan Compelli, Rossen Ivanov, Calin I. Martin, Michail D. Todorov

Articles

The propagation of surface water waves interacting with a current and an uneven bottom is studied. Such a situation is typical for ocean waves where the winds generate currents in the top layer of the ocean. The role of the bottom topography is taken into account since it also influences the local wave and current patterns. Specific scaling of the variables is selected which leads to approximations of Boussinesq and KdV types. The arising KdV equation with variable coefficients, dependent on the bottom topography, is studied numerically when the initial condition is in the form of the one soliton solution …

Quasi-Linearization Method With Rational Legendre Collocation Method For Solving Mhd Flow Over A Stretching Sheet With Variable Thickness And Slip Velocity Which Embedded In A Porous Medium, M. M. Khader

Applications and Applied Mathematics: An International Journal (AAM)

The quasi-linearization method (QLM) and the rational Legendre functions are introduced here to present the numerical solution for the Newtonian fluid flow past an impermeable stretching sheet which embedded in a porous medium with a power-law surface velocity, variable thickness and slip velocity. Firstly, due to the high nonlinearity which yielded from the ordinary differential equation which describes the proposed physical problem, we construct a sequence of linear ODEs by using the QLM, hence the resulted equations become a system of linear algebraic equations. The comparison with the available results in the literature review proves that the obtained results via …

Numerical Studies For Mhd Flow And Gradient Heat Transport Past A Stretching Sheet With Radiation And Heat Production Via Dtm, Khadijah M. Abualnaja

Applications and Applied Mathematics: An International Journal (AAM)

This paper presents a numerically study for the effect of the internal heat generation, magnetic field and thermal radiation effects on the flow and gradient heat transfer of a Newtonian fluid over a stretching sheet. By using a similarity transformation, the governing PDEs can be transformed into a coupled non-linear system of ODEs with variable coefficients. Numerical solutions for these equations subject to appropriate boundary conditions are obtained by using the differential transformation method (DTM). The effects of various physical parameters such as viscosity parameter, the suction parameter, the radiation parameter, internal heat generation or absorption parameter and the Prandtl …