Matching Items (93)
Filtering by
- Genre: Academic theses
- Creators: Huang, Huei-Ping
- Member of: Theses and Dissertations
- Status: Published
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
Rapid expansion of dense beds of fine, spherical particles subjected to rapid depressurization is studied in a vertical shock tube. As the particle bed is unloaded, a high-speed video camera captures the dramatic evolution of the particle bed structure. Pressure transducers are used to measure the dynamic pressure changes during the particle bed expansion process. Image processing, signal processing, and Particle Image Velocimetry techniques, are used to examine the relationships between particle size, initial bed height, bed expansion rate, and gas velocities.
The gas-particle interface and the particle bed as a whole expand and evolve in stages. First, the bed swells nearly homogeneously for a very brief period of time (< 2ms). Shortly afterward, the interface begins to develop instabilities as it continues to rise, with particles nearest the wall rising more quickly. Meanwhile, the bed fractures into layers and then breaks down further into cellular-like structures. The rate at which the structural evolution occurs is shown to be dependent on particle size. Additionally, the rate of the overall bed expansion is shown to be dependent on particle size and initial bed height.
Taller particle beds and beds composed of smaller-diameter particles are found to be associated with faster bed-expansion rates, as measured by the velocity of the gas-particle interface. However, the expansion wave travels more slowly through these same beds. It was also found that higher gas velocities above the the gas-particle interface measured \textit{via} Particle Image Velocimetry or PIV, were associated with particle beds composed of larger-diameter particles. The gas dilation between the shocktube diaphragm and the particle bed interface is more dramatic when the distance between the gas-particle interface and the diaphragm is decreased-as is the case for taller beds.
To further elucidate the complexities of this multiphase compressible flow, simple OpenFOAM (Weller, 1998) simulations of the shocktube experiment were performed and compared to bed expansion rates, pressure fluctuations, and gas velocities. In all cases, the trends and relationships between bed height, particle diameter, with expansion rates, pressure fluctuations and gas velocities matched well between experiments and simulations. In most cases, the experimentally-measured bed rise rates and the simulated bed rise rates matched reasonably well in early times. The trends and overall values of the pressure fluctuations and gas velocities matched well between the experiments and simulations; shedding light on the effects each parameter has on the overall flow.
The gas-particle interface and the particle bed as a whole expand and evolve in stages. First, the bed swells nearly homogeneously for a very brief period of time (< 2ms). Shortly afterward, the interface begins to develop instabilities as it continues to rise, with particles nearest the wall rising more quickly. Meanwhile, the bed fractures into layers and then breaks down further into cellular-like structures. The rate at which the structural evolution occurs is shown to be dependent on particle size. Additionally, the rate of the overall bed expansion is shown to be dependent on particle size and initial bed height.
Taller particle beds and beds composed of smaller-diameter particles are found to be associated with faster bed-expansion rates, as measured by the velocity of the gas-particle interface. However, the expansion wave travels more slowly through these same beds. It was also found that higher gas velocities above the the gas-particle interface measured \textit{via} Particle Image Velocimetry or PIV, were associated with particle beds composed of larger-diameter particles. The gas dilation between the shocktube diaphragm and the particle bed interface is more dramatic when the distance between the gas-particle interface and the diaphragm is decreased-as is the case for taller beds.
To further elucidate the complexities of this multiphase compressible flow, simple OpenFOAM (Weller, 1998) simulations of the shocktube experiment were performed and compared to bed expansion rates, pressure fluctuations, and gas velocities. In all cases, the trends and relationships between bed height, particle diameter, with expansion rates, pressure fluctuations and gas velocities matched well between experiments and simulations. In most cases, the experimentally-measured bed rise rates and the simulated bed rise rates matched reasonably well in early times. The trends and overall values of the pressure fluctuations and gas velocities matched well between the experiments and simulations; shedding light on the effects each parameter has on the overall flow.
ContributorsZunino, Heather (Author) / Adrian, Ronald J (Thesis advisor) / Clarke, Amanda (Committee member) / Chen, Kangping (Committee member) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2019
Description
Laser radars or lidar’s have been used extensively to remotely study winds within the atmospheric boundary layer and atmospheric transport. Lidar sensors have become an important tool within the meteorology and the wind energy community. For example, Doppler lidars are used frequently in wind resource assessment, wind turbine control as well as in atmospheric science research. A Time of Flight based (ToF) direct detection lidar sensor is used in vehicles to navigate through complex and dynamic environments autonomously. These optical sensors are used to map the environment around the car accurately for perception and localization tasks that help achieve complete autonomy.
This thesis begins with a detailed discussion on the fundamentals of a Doppler lidar system. The laser signal flow path to and from the target, the optics of the system and the core signal processing algorithms used to extract velocity information, were studied to get closer to the hardware of a Doppler lidar sensor. A Doppler lidar simulator was built to study the existing signal processing algorithms to detect and estimate doppler frequency, and radial velocity information. Understanding the sensor and its processing at the hardware level is necessary to develop new algorithms to detect and track specific flow structures in the atmosphere. For example, the aircraft vortices have been a topic of extensive research and doppler lidars have proved to be a valuable sensor to detect and track these coherent flow structures. Using the lidar simulator a physics based doppler lidar vortex algorithm is tested on simulated data to track a pair of counter rotating aircraft vortices.
At a system level the major components of a time of flight lidar is very similar to a Doppler lidar. The fundamental physics of operation is however different. While doppler lidars are used for radial velocity measurement, ToF sensors as the name suggests provides precise depth measurements by measuring time of flight between the transmitted and the received pulses. The second part of this dissertation begins to explore the details of ToF lidar system. A system level design, to build a ToF direct detection lidar system is presented. Different lidar sensor modalities that are currently used with sensors in the market today for automotive applications were evaluated and a 2D MEMS based scanning lidar system was designed using off-the shelf components.
Finally, a range of experiments and tests were completed to evaluate the performance of each sub-component of the lidar sensor prototype. A major portion of the testing was done to align the optics of the system and to ensure maximum field of view overlap for the bi-static laser sensor. As a laser range finder, the system demonstrated capabilities to detect hard targets as far as 32 meters. Time to digital converter (TDC) and an analog to digital converter (ADC) was used for providing accurate timing solutions for the lidar prototype. A Matlab lidar model was built and used to perform trade-off studies that helped choosing components to suit the sensor design specifications.
The size, weight and cost of these lidar sensors are still very high and thus making it harder for automotive manufacturers to integrate these sensors into their vehicles. Ongoing research in this field is determined to find a solution that guarantees very high performance in real time and lower its cost over the next decade as components get cheaper and can be seamlessly integrated with cars to improve on-road safety.
This thesis begins with a detailed discussion on the fundamentals of a Doppler lidar system. The laser signal flow path to and from the target, the optics of the system and the core signal processing algorithms used to extract velocity information, were studied to get closer to the hardware of a Doppler lidar sensor. A Doppler lidar simulator was built to study the existing signal processing algorithms to detect and estimate doppler frequency, and radial velocity information. Understanding the sensor and its processing at the hardware level is necessary to develop new algorithms to detect and track specific flow structures in the atmosphere. For example, the aircraft vortices have been a topic of extensive research and doppler lidars have proved to be a valuable sensor to detect and track these coherent flow structures. Using the lidar simulator a physics based doppler lidar vortex algorithm is tested on simulated data to track a pair of counter rotating aircraft vortices.
At a system level the major components of a time of flight lidar is very similar to a Doppler lidar. The fundamental physics of operation is however different. While doppler lidars are used for radial velocity measurement, ToF sensors as the name suggests provides precise depth measurements by measuring time of flight between the transmitted and the received pulses. The second part of this dissertation begins to explore the details of ToF lidar system. A system level design, to build a ToF direct detection lidar system is presented. Different lidar sensor modalities that are currently used with sensors in the market today for automotive applications were evaluated and a 2D MEMS based scanning lidar system was designed using off-the shelf components.
Finally, a range of experiments and tests were completed to evaluate the performance of each sub-component of the lidar sensor prototype. A major portion of the testing was done to align the optics of the system and to ensure maximum field of view overlap for the bi-static laser sensor. As a laser range finder, the system demonstrated capabilities to detect hard targets as far as 32 meters. Time to digital converter (TDC) and an analog to digital converter (ADC) was used for providing accurate timing solutions for the lidar prototype. A Matlab lidar model was built and used to perform trade-off studies that helped choosing components to suit the sensor design specifications.
The size, weight and cost of these lidar sensors are still very high and thus making it harder for automotive manufacturers to integrate these sensors into their vehicles. Ongoing research in this field is determined to find a solution that guarantees very high performance in real time and lower its cost over the next decade as components get cheaper and can be seamlessly integrated with cars to improve on-road safety.
ContributorsBhaskaran, Sreevatsan (Author) / Calhoun, Ronald J (Thesis advisor) / Dahm, Werner (Committee member) / Huang, Huei-Ping (Committee member) / Chen, Kang Pin (Committee member) / Choukulkar, Aditya (Committee member) / Arizona State University (Publisher)
Created2018
Description
Two urban flows are analyzed, one concerned with pollutant transport in a Phoenix, Arizona neighborhood and the other with windshear detection at the Hong Kong International Airport (HKIA).
Lagrangian measures, identified with finite-time Lyapunov exponents, are first used to characterize transport patterns of inertial pollutant particles. Motivated by actual events the focus is on flows in realistic urban geometry. Both deterministic and stochastic transport patterns are identified, as inertial Lagrangian coherent structures. For the deterministic case, the organizing structures are well defined and are extracted at different hours of a day to reveal the variability of coherent patterns. For the stochastic case, a random displacement model for fluid particles is formulated, and used to derive the governing equations for inertial particles to examine the change in organizing structures due to ``zeroth-order'' random noise. It is found that, (1) the Langevin equation for inertial particles can be reduced to a random displacement model; (2) using random noise based on inhomogeneous turbulence, whose diffusivity is derived from $k$-$\epsilon$ models, major coherent structures survive to organize local flow patterns and weaker structures are smoothed out due to random motion.
A study of three-dimensional Lagrangian coherent structures (LCS) near HKIA is then presented and related to previous developments of two-dimensional (2D) LCS analyses in detecting windshear experienced by landing aircraft. The LCS are contrasted among three independent models and against 2D coherent Doppler light detection and ranging (LIDAR) data. Addition of the velocity information perpendicular to the lidar scanning cone helps solidify flow structures inferred from previous studies; contrast among models reveals the intramodel variability; and comparison with flight data evaluates the performance among models in terms of Lagrangian analyses. It is found that, while the three models and the LIDAR do recover similar features of the windshear experienced by a landing aircraft (along the landing trajectory), their Lagrangian signatures over the entire domain are quite different - a portion of each numerical model captures certain features resembling those LCS extracted from independent 2D LIDAR analyses based on observations. Overall, it was found that the Weather Research and Forecast (WRF) model provides the best agreement with the LIDAR data.
Finally, the three-dimensional variational (3DVAR) data assimilation scheme in WRF is used to incorporate the LIDAR line of sight velocity observations into the WRF model forecast at HKIA. Using two different days as test cases, it is found that the LIDAR data can be successfully and consistently assimilated into WRF. Using the updated model forecast LCS are extracted along the LIDAR scanning cone and compare to onboard flight data. It is found that the LCS generated from the updated WRF forecasts are generally better correlated with the windshear experienced by landing aircraft as compared to the LIDAR extracted LCS alone, which suggests that such a data assimilation scheme could be used for the prediction of windshear events.
Lagrangian measures, identified with finite-time Lyapunov exponents, are first used to characterize transport patterns of inertial pollutant particles. Motivated by actual events the focus is on flows in realistic urban geometry. Both deterministic and stochastic transport patterns are identified, as inertial Lagrangian coherent structures. For the deterministic case, the organizing structures are well defined and are extracted at different hours of a day to reveal the variability of coherent patterns. For the stochastic case, a random displacement model for fluid particles is formulated, and used to derive the governing equations for inertial particles to examine the change in organizing structures due to ``zeroth-order'' random noise. It is found that, (1) the Langevin equation for inertial particles can be reduced to a random displacement model; (2) using random noise based on inhomogeneous turbulence, whose diffusivity is derived from $k$-$\epsilon$ models, major coherent structures survive to organize local flow patterns and weaker structures are smoothed out due to random motion.
A study of three-dimensional Lagrangian coherent structures (LCS) near HKIA is then presented and related to previous developments of two-dimensional (2D) LCS analyses in detecting windshear experienced by landing aircraft. The LCS are contrasted among three independent models and against 2D coherent Doppler light detection and ranging (LIDAR) data. Addition of the velocity information perpendicular to the lidar scanning cone helps solidify flow structures inferred from previous studies; contrast among models reveals the intramodel variability; and comparison with flight data evaluates the performance among models in terms of Lagrangian analyses. It is found that, while the three models and the LIDAR do recover similar features of the windshear experienced by a landing aircraft (along the landing trajectory), their Lagrangian signatures over the entire domain are quite different - a portion of each numerical model captures certain features resembling those LCS extracted from independent 2D LIDAR analyses based on observations. Overall, it was found that the Weather Research and Forecast (WRF) model provides the best agreement with the LIDAR data.
Finally, the three-dimensional variational (3DVAR) data assimilation scheme in WRF is used to incorporate the LIDAR line of sight velocity observations into the WRF model forecast at HKIA. Using two different days as test cases, it is found that the LIDAR data can be successfully and consistently assimilated into WRF. Using the updated model forecast LCS are extracted along the LIDAR scanning cone and compare to onboard flight data. It is found that the LCS generated from the updated WRF forecasts are generally better correlated with the windshear experienced by landing aircraft as compared to the LIDAR extracted LCS alone, which suggests that such a data assimilation scheme could be used for the prediction of windshear events.
ContributorsKnutson, Brent (Author) / Tang, Wenbo (Thesis advisor) / Calhoun, Ronald (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Mahalov, Alex (Committee member) / Arizona State University (Publisher)
Created2018
Numerical simulation of environmental flow over urban landscape for applications to renewable energy
Description
Development of renewable energy solutions has become a major interest among environmental organizations and governments around the world due to an increase in energy consumption and global warming. One fast growing renewable energy solution is the application of wind energy in cities. To qualitative and quantitative predict wind turbine performance in urban areas, CFD simulation is performed on real-life urban geometry and wind velocity profiles are evaluated. Two geometries in Arizona is selected in this thesis to demonstrate the influence of building heights; one of the simulation models, ASU campus, is relatively low rise and without significant tall buildings; the other model, the downtown phoenix model, are high-rise and with greater building height difference. The content of this thesis focuses on using RANS computational fluid dynamics approach to simulate wind acceleration phenomenon in two complex geometries, ASU campus and Phoenix downtown model. Additionally, acceleration ratio and locations are predicted, the results are then used to calculate the best location for small wind turbine installments.
ContributorsYing, Xiaoyan (Author) / Huang, Huei-Ping (Thesis advisor) / Peet, Yulia (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2015
Description
For the increasing concerns of influence on environment by fossil-electricity generation, application of renewable energy becomes one of the most focused issues in society. Based on the limitation on urban environment, wind turbines, which can be mounted on rooftop or between buildings, are regarded as a feasible way for wind energy generation. This study presents wind flow simulations in a large-scale environment with certain dimension buildings. Different inlet velocity boundary conditions are tested firstly, and the non-uniform inlet boundary condition shows better agreement with realistic situation. Turbulence intensity is set to be 10% for comparison consistency. The k-epsilon turbulence model is regarded as a better simulation for this certain condition. After that, three different structures, which include single building, pristine double building and modified circular gap double building systems, are tested in this environment condition. The result shows 18.8% velocity increasing on the top of single building system. Pristine double building systems are tested with 4 different gap distances, and building with 10 meters gap achieved the best velocity condition, which 32.8% velocity increasing and 11.8% improvement comparing to single building system, respectively. But the location of maximum velocity moves to the gap and the maximum velocity on the rooftop of double building system is approximately 5.1% lower than single building system. Based on previous study, modified circular gap double building system is created with 10 meters gap. Comparing result with single building system, modified circular gap system achieves higher improvement for wind flow, whose improvement of velocity increasing in the gap and on the rooftop of building are 47.1% and 3.0%, respectively. As a result, the modified circular gap double building can be regarded as a high efficiency system of environmental wind flow over buildings for renewable energy system.
ContributorsLi, Guoyi (Author) / Huang, Huei-Ping (Thesis advisor) / Lee, Taewoo (Committee member) / Forzani, Erica (Committee member) / Arizona State University (Publisher)
Created2015
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
Description
A moving overlapping mesh methodology that achieves spectral accuracy in space and up to second-order accuracy in time is developed for solution of unsteady incompressible flow equations in three-dimensional domains. The targeted applications are in aerospace and mechanical engineering domains and involve problems in turbomachinery, rotary aircrafts, wind turbines and others. The methodology is built within the dual-session communication framework initially developed for stationary overlapping meshes. The methodology employs semi-implicit spectral element discretization of equations in each subdomain and explicit treatment of subdomain interfaces with spectrally-accurate spatial interpolation and high-order accurate temporal extrapolation, and requires few, if any, iterations, yet maintains the global accuracy and stability of the underlying flow solver. Mesh movement is enabled through the Arbitrary Lagrangian-Eulerian formulation of the governing equations, which allows for prescription of arbitrary velocity values at discrete mesh points.
The stationary and moving overlapping mesh methodologies are thoroughly validated using two- and three-dimensional benchmark problems in laminar and turbulent flows. The spatial and temporal global convergence, for both methods, is documented and is in agreement with the nominal order of accuracy of the underlying solver.
Stationary overlapping mesh methodology was validated to assess the influence of long integration times and inflow-outflow global boundary conditions on the performance. In a turbulent benchmark of fully-developed turbulent pipe flow, the turbulent statistics are validated against the available data.
Moving overlapping mesh simulations are validated on the problems of two-dimensional oscillating cylinder and a three-dimensional rotating sphere. The aerodynamic forces acting on these moving rigid bodies are determined, and all results are compared with published data. Scaling tests, with both methodologies, show near linear strong scaling, even for moderately large processor counts.
The moving overlapping mesh methodology is utilized to investigate the effect of an upstream turbulent wake on a three-dimensional oscillating NACA0012 extruded airfoil. A direct numerical simulation (DNS) at Reynolds Number 44,000 is performed for steady inflow incident upon the airfoil oscillating between angle of attack 5.6 and 25 degrees with reduced frequency k=0.16. Results are contrasted with subsequent DNS of the same oscillating airfoil in a turbulent wake generated by a stationary upstream cylinder.
The stationary and moving overlapping mesh methodologies are thoroughly validated using two- and three-dimensional benchmark problems in laminar and turbulent flows. The spatial and temporal global convergence, for both methods, is documented and is in agreement with the nominal order of accuracy of the underlying solver.
Stationary overlapping mesh methodology was validated to assess the influence of long integration times and inflow-outflow global boundary conditions on the performance. In a turbulent benchmark of fully-developed turbulent pipe flow, the turbulent statistics are validated against the available data.
Moving overlapping mesh simulations are validated on the problems of two-dimensional oscillating cylinder and a three-dimensional rotating sphere. The aerodynamic forces acting on these moving rigid bodies are determined, and all results are compared with published data. Scaling tests, with both methodologies, show near linear strong scaling, even for moderately large processor counts.
The moving overlapping mesh methodology is utilized to investigate the effect of an upstream turbulent wake on a three-dimensional oscillating NACA0012 extruded airfoil. A direct numerical simulation (DNS) at Reynolds Number 44,000 is performed for steady inflow incident upon the airfoil oscillating between angle of attack 5.6 and 25 degrees with reduced frequency k=0.16. Results are contrasted with subsequent DNS of the same oscillating airfoil in a turbulent wake generated by a stationary upstream cylinder.
ContributorsMerrill, Brandon Earl (Author) / Peet, Yulia (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2016
Optical simulation and optimization of light extraction efficiency for organic light emitting diodes
Description
Current organic light emitting diodes (OLEDs) suffer from the low light extraction efficiency. In this thesis, novel OLED structures including photonic crystal, Fabry-Perot resonance cavity and hyperbolic metamaterials were numerically simulated and theoretically investigated. Finite-difference time-domain (FDTD) method was employed to numerically simulate the light extraction efficiency of various 3D OLED structures. With photonic crystal structures, a maximum of 30% extraction efficiency is achieved. A higher external quantum efficiency of 35% is derived after applying Fabry-Perot resonance cavity into OLEDs. Furthermore, different factors such as material properties, layer thicknesses and dipole polarizations and locations have been studied. Moreover, an upper limit for the light extraction efficiency of 80% is reached theoretically with perfect reflector and single dipole polarization and location. To elucidate the physical mechanism, transfer matrix method is introduced to calculate the spectral-hemispherical reflectance of the multilayer OLED structures. In addition, an attempt of using hyperbolic metamaterial in OLED has been made and resulted in 27% external quantum efficiency, due to the similar mechanism of wave interference as Fabry-Perot structure. The simulation and optimization methods and findings would facilitate the design of next generation, high-efficiency OLED devices.
ContributorsSu, Hang (Author) / Wang, Liping (Thesis advisor) / Li, Jian (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2016
Description
Analytical solution of the pressure field for water uptake through a composite root, coupled with fully saturated soil is derived by using the slender body approximation. It is shown that in general, the resistance of the root and soil are not additive. This result can play a very important role in modelling water uptake through plant roots and determination of hydraulic resistances of plant roots. Optimum plant root structure that minimizes a single root’s hydraulic resistance is also studied in this work with the constraint of prescribed root volume. Hydraulic resistances under the slender body approximation and without such a limitation are considered. It is found that for large stele-to-cortex permeability ratio, there exists an optimum root length-to-base-radius ratio that minimizes the hydraulic resistance. A remarkable feature of the optimum root structure is that the optimum dimensionless stele conductivity depends only on a single geometrical parameter, the stele-to-root base-radius ratio. Once the stele-to-root base-radius ratio and the stele-to-cortex permeability ratio are given, the optimum root length-to-radius ratio can be found. While these findings remain to be verified by experiments for real plant roots, they offer theoretical guidance for the design of bio-inspired structures that minimizes hydraulic resistance for fluid production from porous media.
ContributorsChandrashekar, Sriram (Author) / Chen, Kang-Ping (Thesis advisor) / Huang, Huei-Ping (Committee member) / Rykaczewski, Konrad (Committee member) / Arizona State University (Publisher)
Created2016