• Bridson, R., 2008: Fluid Simulation for Computer Graphics.A K Peters, 228 pp.

  • Choi, J.-I., , Oberoi R. C. , , Edwards J. R. , , and Rosati J. A. , 2007: An immersed boundary method for complex incompressible flows. J. Comput. Phys., 224, 757784, doi:10.1016/j.jcp.2006.10.032.

    • Search Google Scholar
    • Export Citation
  • DeLeon, R., , Felzien K. , , and Senocak I. , 2012: Toward a GPU-accelerated immersed boundary method for wind forecasting over complex terrain. Proc. ASME Fluids Engineering Division Summer Meeting, Rio Grande, Puerto Rico, ASME, FEDSM2012-72145.

  • Eberly, D., 2008: Distance between point and line, ray, or line segment. Geometric Tools, 2 pp. [Available online at http://www.geometrictools.com/Documentation/DistancePointLine.pdf.]

  • Edwards, J. R., , Choi J.-L. , , Ghosh S. , , Gieseking D. A. , , and Eischen J. D. , 2010: An immersed boundary method for general flow applications. Proc. Third Joint US-European Fluids Engineering Summer Meeting, Montreal, QC, Canada, FEDSM-ICNMM2010-31097.

  • Fadlun, E. A., , Verzicco R. , , Orlandi P. , , and Mohd-Yusof J. , 2000: Combined immersed-boundary finite-difference methods for three-dimensional complex flow simulations. J. Comput. Phys., 161, 3560, doi:10.1006/jcph.2000.6484.

    • Search Google Scholar
    • Export Citation
  • Gilmanov, A., , and Sotiropoulos F. , 2005: A hybrid Cartesian/immersed boundary method for simulating flows with 3D, geometrically complex, moving bodies. J. Comput. Phys., 207, 457492, doi:10.1016/j.jcp.2005.01.020.

    • Search Google Scholar
    • Export Citation
  • Guth, P. L., 2013: MICRODEM. [Available online at http://www.usna.edu/Users/oceano/pguth/website/microdem/microdem.htm.]

  • Holcombe, S., 2012: Inpolyhedron—Are points inside a triangulated volume? MATLAB Central File Exchange. [Available online at http://www.mathworks.com/matlabcentral/fileexchange/37856-inpolyhedron-are-points-inside-a-triangulated-volume-.]

  • Kang, S., , Lightbody A. , , Hill C. , , and Sotiropoulos F. , 2011: High-resolution numerical simulation of turbulence in natural waterways. Adv. Water Resour., 34, 98113, doi:10.1016/j.advwatres.2010.09.018.

    • Search Google Scholar
    • Export Citation
  • Lundquist, K. A., , Chow F. K. , , and Lundquist J. K. , 2010: An immersed boundary method for the Weather Research and Forecasting Model. Mon. Wea. Rev., 138, 796817, doi:10.1175/2009MWR2990.1.

    • Search Google Scholar
    • Export Citation
  • Lundquist, K. A., , Chow F. K. , , and Lundquist J. K. , 2012: An immersed boundary method enabling large-eddy simulations of flow over complex terrain in the WRF Model. Mon. Wea. Rev., 140, 39363955, doi:10.1175/MWR-D-11-00311.1.

    • Search Google Scholar
    • Export Citation
  • McDonald, B., 2004: SURF2STL. MATLAB Central File Exchange. [Available online at http://www.mathworks.com/matlabcentral/fileexchange/4512-surf2stl.]

  • Mittal, R., , and Iaccarino G. , 2005: Immersed boundary methods. Annu. Rev. Fluid Mech., 37, 239261, doi:10.1146/annurev.fluid.37.061903.175743.

    • Search Google Scholar
    • Export Citation
  • Mittal, R., , Dong H. , , Bozkurttas M. , , Najjar F. , , Vargas A. , , and von Loebbeck A. , 2008: A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries. J. Comput. Phys., 227, 48254852, doi:10.1016/j.jcp.2008.01.028.

    • Search Google Scholar
    • Export Citation
  • Mohd-Yusof, J., 1997: Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries. Annual research briefs 1997, Center for Turbulence Research Annual Research Briefs, 317–327.

  • O’Rourke, J., 1994: Computational Geometry in C.Cambridge University Press, 346 pp.

  • Osher, S., , and Fedkiw R. , 2003: Level Set Methods and Dynamic Implicit Surfaces.Springer, 269 pp.

  • Peskin, C., 1972: Flow patterns around heart valves: A numerical method. J. Comput. Phys., 10, 252271, doi:10.1016/0021-9991(72)90065-4.

    • Search Google Scholar
    • Export Citation
  • Pu, J., , Zheng B. , , Leader J. K. , , Wang X.-H. , , and Gur D. , 2008: An automated CT based lung nodule detection scheme using geometric analysis of signed distance field. Med. Phys., 35, 34533461, doi:10.1118/1.2948349.

    • Search Google Scholar
    • Export Citation
  • Senocak, I., , Ackerman A. S. , , Stevens D. E. , , and Mansour N. N. , 2004: Topography modeling in atmospheric flows using the immersed boundary method. Annual research briefs 2004, Center for Turbulence Research, 331–341.

  • Sethian, J., 1999: Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science.2nd ed. Cambridge University Press, 378 pp.

  • Sunday, D., 2001a: Intersection of a ray/segment with a triangle. [Available online at http://geomalgorithms.com/a06-_intersect-2.html#Segment-Triangle.]

  • Sunday, D., 2001b: Lines and distance of a point to a line. [Available online at http://geomalgorithms.com/a02-_lines.html.]

  • Zhao, H., 2004: A fast sweeping method for Eikonal equations. Math. Comput., 74, 603627, doi:10.1090/S0025-5718-04-01678-3.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 38 38 3
PDF Downloads 29 29 3

An Immersed Boundary Geometric Preprocessor for Arbitrarily Complex Terrain and Geometry

View More View Less
  • 1 Department of Mechanical and Biomedical Engineering, Boise State University, Boise, Idaho
  • | 2 Department of Computer Science, Boise State University, Boise, Idaho
© Get Permissions
Restricted access

Abstract

There is a growing interest to apply the immersed boundary method to compute wind fields over arbitrarily complex terrain. The computer implementation of an immersed boundary module into an existing flow solver can be accomplished with minor modifications to the rest of the computer program. However, a versatile preprocessor is needed at the first place to extract the essential geometric information pertinent to the immersion of an arbitrarily complex terrain inside a 3D Cartesian mesh. Errors in the geometric information can negatively impact the correct implementation of the immersed boundary method as part of the solution algorithm. Additionally, the distance field from the terrain is needed to implement various subgrid-scale turbulence models and to initialize wind fields over complex terrain. Despite the popularity of the immersed boundary method, procedures used in the geometric preprocessing stage have received less attention. The present study found that concave and convex regions of complex terrain are particularly challenging to process with existing procedures discussed in the literature. To address this issue, a geometric preprocessor with a distance field solver was presented, and the solver demonstrated its versatility for arbitrarily complex geometry, terrain, and urban environments. The distance field solver uses the initial distance field at the immersed boundaries and propagates it to the rest of the domain by solving the Eikonal equation with the fast sweeping method.

Denotes Open Access content.

Corresponding author address: Inanc Senocak, Department of Mechanical and Biomedical Engineering, Boise State University, 1910 University Drive, Boise, ID 83725-2085. E-mail: senocak@boisestate.edu

Abstract

There is a growing interest to apply the immersed boundary method to compute wind fields over arbitrarily complex terrain. The computer implementation of an immersed boundary module into an existing flow solver can be accomplished with minor modifications to the rest of the computer program. However, a versatile preprocessor is needed at the first place to extract the essential geometric information pertinent to the immersion of an arbitrarily complex terrain inside a 3D Cartesian mesh. Errors in the geometric information can negatively impact the correct implementation of the immersed boundary method as part of the solution algorithm. Additionally, the distance field from the terrain is needed to implement various subgrid-scale turbulence models and to initialize wind fields over complex terrain. Despite the popularity of the immersed boundary method, procedures used in the geometric preprocessing stage have received less attention. The present study found that concave and convex regions of complex terrain are particularly challenging to process with existing procedures discussed in the literature. To address this issue, a geometric preprocessor with a distance field solver was presented, and the solver demonstrated its versatility for arbitrarily complex geometry, terrain, and urban environments. The distance field solver uses the initial distance field at the immersed boundaries and propagates it to the rest of the domain by solving the Eikonal equation with the fast sweeping method.

Denotes Open Access content.

Corresponding author address: Inanc Senocak, Department of Mechanical and Biomedical Engineering, Boise State University, 1910 University Drive, Boise, ID 83725-2085. E-mail: senocak@boisestate.edu
Save