Matching Items (13)
Filtering by
- All Subjects: Computational Fluid Dynamics
- Creators: Huang, Huei-Ping
Description
Image processing in canals, rivers and other bodies of water has been a very important concern. This research using Image Processing was performed to obtain a photographic evidence of the data of the site which helps in monitoring the conditions of the water body and the surroundings. Images are captured using a digital camera and the images are stored onto a datalogger, these images are retrieved using a cellular/ satellite modem. A MATLAB program was designed to obtain the level of water by just entering the file name into to the program, a curve fit model was created to determine the contrast parameters. The contrast parameters were obtained using the data obtained from the gray scale image mainly the mean and variance of the intensity values. The enhanced images are used to determine the level of water by taking pixel intensity plots along the region of interest. The level of water obtained is accurate to less than 2% of the actual level of water observed from the image. High speed imaging in micro channels have various application in industrial field, medical field etc. In medical field it is tested by using blood samples. The experimental procedure proposed determines the flow duration and the defects observed in these channel using a fluid introduced into the micro channel the fluid being water based dye and whole milk. The viscosity of the fluid shows different types of flow patterns and defects in the micro channel. The defects observed vary from a small effect to the flow pattern to an extreme defect in the channel such as obstruction of flow or deformation in the channel. The sample needs to be further analyzed by SEM to get a better insight on the defects.
ContributorsShasedhara, Abhijeet Bangalore (Author) / Lee, Taewoo (Thesis advisor) / Huang, Huei-Ping (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2011
Description
A new method of adaptive mesh generation for the computation of fluid flows is investigated. The method utilizes gradients of the flow solution to adapt the size and stretching of elements or volumes in the computational mesh as is commonly done in the conventional Hessian approach. However, in the new method, higher-order gradients are used in place of the Hessian. The method is applied to the finite element solution of the incompressible Navier-Stokes equations on model problems. Results indicate that a significant efficiency benefit is realized.
ContributorsShortridge, Randall (Author) / Chen, Kang Ping (Thesis advisor) / Herrmann, Marcus (Thesis advisor) / Wells, Valana (Committee member) / Huang, Huei-Ping (Committee member) / Mittelmann, Hans (Committee member) / Arizona State University (Publisher)
Created2011
Description
The evolution of single hairpin vortices and multiple interacting hairpin vortices are studied in direct numerical simulations of channel flow at Re-tau=395. The purpose of this study is to observe the effects of increased Reynolds number and varying initial conditions on the growth of hairpins and the conditions under which single hairpins autogenerate hairpin packets. The hairpin vortices are believed to provide a unified picture of wall turbulence and play an important role in the production of Reynolds shear stress which is directly related to turbulent drag. The structures of the initial three-dimensional vortices are extracted from the two-point spatial correlation of the fully turbulent direct numerical simulation of the velocity field by linear stochastic estimation and embedded in a mean flow having the profile of the fully turbulent flow. The Reynolds number of the present simulation is more than twice that of the Re-tau=180 flow from earlier literature and the conditional events used to define the stochastically estimated single vortex initial conditions include a number of new types of events such as quasi-streamwise vorticity and Q4 events. The effects of parameters like strength, asymmetry and position are evaluated and compared with existing results in the literature. This study then attempts to answer questions concerning how vortex mergers produce larger scale structures, a process that may contribute to the growth of length scale with increasing distance from the wall in turbulent wall flows. Multiple vortex interactions are studied in detail.
ContributorsParthasarathy, Praveen Kumar (Author) / Adrian, Ronald (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2011
Description
For CFD validation, hypersonic flow fields are simulated and compared with experimental data specifically designed to recreate conditions found by hypersonic vehicles. Simulated flow fields on a cone-ogive with flare at Mach 7.2 are compared with experimental data from NASA Ames Research Center 3.5" hypersonic wind tunnel. A parametric study of turbulence models is presented and concludes that the k-kl-omega transition and SST transition turbulence model have the best correlation. Downstream of the flare's shockwave, good correlation is found for all boundary layer profiles, with some slight discrepancies of the static temperature near the surface. Simulated flow fields on a blunt cone with flare above Mach 10 are compared with experimental data from CUBRC LENS hypervelocity shock tunnel. Lack of vibrational non-equilibrium calculations causes discrepancies in heat flux near the leading edge. Temperature profiles, where non-equilibrium effects are dominant, are compared with the dissociation of molecules to show the effects of dissociation on static temperature. Following the validation studies is a parametric analysis of a hypersonic inlet from Mach 6 to 20. Compressor performance is investigated for numerous cowl leading edge locations up to speeds of Mach 10. The variable cowl study showed positive trends in compressor performance parameters for a range of Mach numbers that arise from maximizing the intake of compressed flow. An interesting phenomenon due to the change in shock wave formation for different Mach numbers developed inside the cowl that had a negative influence on the total pressure recovery. Investigation of the hypersonic inlet at different altitudes is performed to study the effects of Reynolds number, and consequently, turbulent viscous effects on compressor performance. Turbulent boundary layer separation was noted as the cause for a change in compressor performance parameters due to a change in Reynolds number. This effect would not be noticeable if laminar flow was assumed. Mach numbers up to 20 are investigated to study the effects of vibrational and chemical non-equilibrium on compressor performance. A direct impact on the trends on the kinetic energy efficiency and compressor efficiency was found due to dissociation.
ContributorsOliden, Daniel (Author) / Lee, Tae-Woo (Thesis advisor) / Peet, Yulia (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2013
Description
An understanding of aerodynamics is crucial for automobile performance and efficiency. There are many types of “add-on” aerodynamic devices for cars including wings, splitters, and vortex generators. While these have been studied extensively, rear spoilers have not, and their effects are not as widely known. A Computational Fluid Dynamics (CFD) and wind tunnel study was performed to study the effects of spoilers on vehicle aerodynamics and performance. Vehicle aerodynamics is geometry dependent, meaning what applies to one car may or may not apply on another. So, the Scion FRS was chosen as the test vehicle because it is has the “classic” sports car configuration with a long hood, short rear, and 2+2 passenger cabin while also being widely sold with a plethora of aftermarket aerodynamic modifications available. Due to computing and licensing restrictions, only a 2D CFD simulation was performed in ANSYS Fluent 19.1. A surface model of the centerline of the car was created in SolidWorks and imported into ANSYS, where the domain was created. A mesh convergence study was run to determine the optimum mesh size, and Realizable k-epsilon was the chosen physics model. The wind tunnel lacked equipment to record quantifiable data, so the wind tunnel was utilized for flow visualization on a 1/24 scale car model to compare with the CFD.
0° spoilers reduced the wake area behind the car, decreasing pressure drag but also decreasing underbody flow, causing a reduction in drag and downforce. Angled spoilers increased the wake area behind the car, increasing pressure drag but also increasing underbody flow, causing an increase in drag and downforce. Longer spoilers increased these effects compared to shorter spoilers, and short spoilers at different angles did not create significantly different effects. 0° spoilers would be best suited for cases that prioritize fuel economy or straight-line acceleration and speed due to the drag reduction, while angled spoilers would be best suited for cars requiring downforce. The angle and length of spoiler would depend on the downforce needed, which is dependent on the track.
0° spoilers reduced the wake area behind the car, decreasing pressure drag but also decreasing underbody flow, causing a reduction in drag and downforce. Angled spoilers increased the wake area behind the car, increasing pressure drag but also increasing underbody flow, causing an increase in drag and downforce. Longer spoilers increased these effects compared to shorter spoilers, and short spoilers at different angles did not create significantly different effects. 0° spoilers would be best suited for cases that prioritize fuel economy or straight-line acceleration and speed due to the drag reduction, while angled spoilers would be best suited for cars requiring downforce. The angle and length of spoiler would depend on the downforce needed, which is dependent on the track.
ContributorsNie, Alexander (Author) / Wells, Valana (Thesis director) / Huang, Huei-Ping (Committee member) / Mechanical and Aerospace Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2019-12
Description
Advancements to a dual scale Large Eddy Simulation (LES) modeling approach for immiscible turbulent phase interfaces are presented. In the dual scale LES approach, a high resolution auxiliary grid, used to capture a fully resolved interface geometry realization, is linked to an LES grid that solves the filtered Navier-Stokes equations. Exact closure of the sub-filter interface terms is provided by explicitly filtering the fully resolved quantities from the auxiliary grid. Reconstructing a fully resolved velocity field to advance the phase interface requires modeling several sub-filter effects, including shear and accelerational instabilities and phase change. Two sub-filter models were developed to generate these sub-filter hydrodynamic instabilities: an Orr-Sommerfeld model and a Volume-of-Fluid (VoF) vortex sheet method. The Orr-Sommerfeld sub-filter model was found to be incompatible with the dual scale approach, since it is unable to generate interface rollup and a process to separate filtered and sub-filter scales could not be established. A novel VoF vortex sheet method was therefore proposed, since prior vortex methods have demonstrated interface rollup and following the LES methodology, the vortex sheet strength could be decomposed into its filtered and sub-filter components. In the development of the VoF vortex sheet method, it was tested with a variety of classical hydrodynamic instability problems, compared against prior work and linear theory, and verified using Direct Numerical Simulations (DNS). An LES consistent approach to coupling the VoF vortex sheet with the LES filtered equations is presented and compared against DNS. Finally, a sub-filter phase change model is proposed and assessed in the dual scale LES framework with an evaporating interface subjected to decaying homogeneous isotropic turbulence. Results are compared against DNS and the interplay between surface tension forces and evaporation are discussed.
ContributorsGoodrich, Austin Chase (Author) / Herrmann, Marcus (Thesis advisor) / Dahm, Werner (Committee member) / Kim, Jeonglae (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Arizona State University (Publisher)
Created2023
Description
This study identifies the influence that leading-edge shape has on the aerodynamic characteristics of a wing using surface far-field and near-field analysis. It examines if a wake survey is the appropriate means for measuring profile drag and induced drag. The paper unveils the differences between sharp leading-edge and blunt leading-edge wings with the tools of pressure loop, chordwise pressure distribution, span load plots and with wake integral computations. The analysis was performed using Computational Fluid Dynamics (CFD), vortex lattice potential flow code (VORLAX), and a few wind-tunnels runs to acquire data for comparison. This study found that sharp leading-edge wings have less leading-edge suction and higher drag than blunt leading-edge wings.
The blunt leading-edge wings have less drag because the normal vector of the surface in the front section of the airfoil develops forces at opposed skin friction. The shape of the leading edge, in conjunction with the effect of viscosity, slightly alter the span load; both the magnitude of the lift and the transverse distribution. Another goal in this study is to verify the veracity of wake survey theory; the two different leading-edge shapes reveals the shortcoming of Mclean’s equation which is only applicable to blunt leading-edge wings.
The blunt leading-edge wings have less drag because the normal vector of the surface in the front section of the airfoil develops forces at opposed skin friction. The shape of the leading edge, in conjunction with the effect of viscosity, slightly alter the span load; both the magnitude of the lift and the transverse distribution. Another goal in this study is to verify the veracity of wake survey theory; the two different leading-edge shapes reveals the shortcoming of Mclean’s equation which is only applicable to blunt leading-edge wings.
ContributorsOu, Che Wei (Author) / Takahashi, Timothy (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2019
Description
This dissertation introduces FARCOM (Fortran Adaptive Refiner for Cartesian Orthogonal Meshes), a new general library for adaptive mesh refinement (AMR) based on an unstructured hexahedral mesh framework. As a result of the underlying unstructured formulation, the refinement and coarsening operators of the library operate on a single-cell basis and perform in-situ replacement of old mesh elements. This approach allows for h-refinement without the memory and computational expense of calculating masked coarse grid cells, as is done in traditional patch-based AMR approaches, and enables unstructured flow solvers to have access to the automated domain generation capabilities usually only found in tree AMR formulations.
The library is written to let the user determine where to refine and coarsen through custom refinement selector functions for static mesh generation and dynamic mesh refinement, and can handle smooth fields (such as level sets) or localized markers (e.g. density gradients). The library was parallelized with the use of the Zoltan graph-partitioning library, which provides interfaces to both a graph partitioner (PT-Scotch) and a partitioner based on Hilbert space-filling curves. The partitioned adjacency graph, mesh data, and solution variable data is then packed and distributed across all MPI ranks in the simulation, which then regenerate the mesh, generate domain decomposition ghost cells, and create communication caches.
Scalability runs were performed using a Leveque wave propagation scheme for solving the Euler equations. The results of simulations on up to 1536 cores indicate that the parallel performance is highly dependent on the graph partitioner being used, and differences between the partitioners were analyzed. FARCOM is found to have better performance if each MPI rank has more than 60,000 cells.
The library is written to let the user determine where to refine and coarsen through custom refinement selector functions for static mesh generation and dynamic mesh refinement, and can handle smooth fields (such as level sets) or localized markers (e.g. density gradients). The library was parallelized with the use of the Zoltan graph-partitioning library, which provides interfaces to both a graph partitioner (PT-Scotch) and a partitioner based on Hilbert space-filling curves. The partitioned adjacency graph, mesh data, and solution variable data is then packed and distributed across all MPI ranks in the simulation, which then regenerate the mesh, generate domain decomposition ghost cells, and create communication caches.
Scalability runs were performed using a Leveque wave propagation scheme for solving the Euler equations. The results of simulations on up to 1536 cores indicate that the parallel performance is highly dependent on the graph partitioner being used, and differences between the partitioners were analyzed. FARCOM is found to have better performance if each MPI rank has more than 60,000 cells.
ContributorsBallesteros, Carlos Alberto (Author) / Herrmann, Marcus (Thesis advisor) / Adrian, Ronald (Committee member) / Chen, Kangping (Committee member) / Huang, Huei-Ping (Committee member) / Lopez, Juan (Committee member) / Arizona State University (Publisher)
Created2019
Description
The goal of this paper was to do an analysis of two-dimensional unsplit mass and momentum conserving Finite Volume Methods for Advection for Volume of Fluid Fields with interfaces and validating their rates of convergence. Specifically three unsplit transport methods and one split transport method were amalgamated individually with four Piece-wise Linear Reconstruction Schemes (PLIC) i.e. Unsplit Eulerian Advection (UEA) by Owkes and Desjardins (2014), Unsplit Lagrangian Advection (ULA) by Yang et al. (2010), Split Lagrangian Advection (SLA) by Scardovelli and Zaleski (2003) and Unsplit Averaged Eulerian-Lagrangian Advection (UAELA) with two Finite Difference Methods by Parker and Youngs (1992) and two Error Minimization Methods by Pilliod Jr and Puckett (2004). The observed order of accuracy was first order in all cases except when unsplit methods and error minimization methods were used consecutively in each iteration, which resulted in second-order accuracy on the shape error convergence. The Averaged Unsplit Eulerian-Lagrangian Advection (AUELA) did produce first-order accuracy but that was due to a temporal error in the numerical setup. The main unsplit methods, Unsplit Eulerian Advection (UEA) and Unsplit Lagrangian Advection (ULA), preserve mass and momentum and require geometric clipping to solve two-phase fluid flows. The Unsplit Lagrangian Advection (ULA) can allow for small divergence in the velocity field perhaps saving time on the iterative solver of the variable coefficient Poisson System.
ContributorsAnsari, Adil (M.S.) (Author) / Herrmann, Marcus (Thesis advisor) / Peet, Yulia (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2019
Description
First, in a large-scale structure, a 3-D CFD model was built to simulate flow and temperature distributions. The flow patterns and temperature distributions are characterized and validated through spot measurements. The detailed understanding of them then allows for optimization of the HVAC configuration because identification of the problematic flow patterns and temperature mis-distributions leads to some corrective measures. Second, an appropriate form of the viscous dissipation term in the integral form of the conservation equation was considered, and the effects of momentum terms on the computed drop size in pressure-atomized sprays were examined. The Sauter mean diameter (SMD) calculated in this manner agrees well with experimental data of the drop velocities and sizes. Using the suggested equation with the revised treatment of liquid momentum setup, injection parameters can be directly input to the system of equations. Thus, this approach is capable of incorporating the effects of injection parameters for further considerations of the drop and velocity distributions under a wide range of spray geometry and injection conditions. Lastly, groundwater level estimation was investigated using compressed sensing (CS). To satisfy a general property of CS, a random measurement matrix was used, the groundwater network was constructed, and finally the l-1 optimization was run. Through several validation tests, correct estimation of groundwater level by CS was shown. Using this setup, decreasing trends in groundwater level in the southwestern US was shown. The suggested method is effective in that the total measurements of registered wells can be reduced down by approximately 42 %, sparse data can be visualized and a possible approach for groundwater management during extreme weather changes, e.g. in California, was demonstrated.
ContributorsLee, Joon Young (Author) / Lee, Taewoo (Thesis advisor) / Huang, Huei-Ping (Committee member) / Lopez, Juan (Committee member) / Phelan, Patrick (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015