Exploring the Interplay between CAD and FreeFem++ as an Energy Decision-Making Tool for Architectural Design

The energy modelling software tools commonly used for architectural purposes do not allow a straightforward real-time implementation within the architectural design programs. In addition, the surrounding exterior spaces of the building, including the inner courtyards, hardly present a specific treatment distinguishing these spaces from the general external temperature in the thermal simulations. This is a clear disadvantage when it comes to streamlining the design process in relation to the whole-building energy optimization. In this context, the present study aims to demonstrate the advantages of the FreeFem++ open source program for performing simulations in architectural environments. These simulations include microclimate tests that describe the interactions between a building architecture and its local exterior. The great potential of this mathematical tool can be realized through its complete system integration within CAD (Computer-Aided Design) software such as SketchUp or AutoCAD. In order to establish the suitability of FreeFem++ for the performance of simulations, the most widely employed energy simulation tools able to consider a proposed architectural geometry in a specific environment are compared. On the basis of this analysis, it can be concluded that FreeFem++ is the only program displaying the best features for the thermal performance simulation of these specific outdoor spaces, excluding the currently unavailable easy interaction with architectural drawing programs. The main contribution of this research is, in fact, the enhancement of FreeFem++ usability by proposing a simple intuitive method for the creation of building geometries and their respective meshing (pre-processing). FreeFem++ is also considered a tool for data analysis (post-processing) able to help engineers and architects with building energy-efficiency-related tasks

Automated classification of three-dimensional reconstructions of coral reefs using convolutional neural networks

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Hopkinson, B. M., King, A. C., Owen, D. P., Johnson-Roberson, M., Long, M. H., & Bhandarkar, S. M. Automated classification of three-dimensional reconstructions of coral reefs using convolutional neural networks. PLoS One, 15(3), (2020): e0230671, doi: 10.1371/journal.pone.0230671.Coral reefs are biologically diverse and structurally complex ecosystems, which have been severally affected by human actions. Consequently, there is a need for rapid ecological assessment of coral reefs, but current approaches require time consuming manual analysis, either during a dive survey or on images collected during a survey. Reef structural complexity is essential for ecological function but is challenging to measure and often relegated to simple metrics such as rugosity. Recent advances in computer vision and machine learning offer the potential to alleviate some of these limitations. We developed an approach to automatically classify 3D reconstructions of reef sections and assessed the accuracy of this approach. 3D reconstructions of reef sections were generated using commercial Structure-from-Motion software with images extracted from video surveys. To generate a 3D classified map, locations on the 3D reconstruction were mapped back into the original images to extract multiple views of the location. Several approaches were tested to merge information from multiple views of a point into a single classification, all of which used convolutional neural networks to classify or extract features from the images, but differ in the strategy employed for merging information. Approaches to merging information entailed voting, probability averaging, and a learned neural-network layer. All approaches performed similarly achieving overall classification accuracies of ~96% and >90% accuracy on most classes. With this high classification accuracy, these approaches are suitable for many ecological applications.This study was funded by grants from the Alfred P. Sloan Foundation (BMH, BR2014-049; https://sloan.org), and the National Science Foundation (MHL, OCE-1657727; https://www.nsf.gov). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript

A 8-neighbor model lattice Boltzmann method applied to mathematical-physical equations

© 2016. This version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/A lattice Boltzmann method (LBM) 9-bit model is presented to solve mathematical-physical equations, such as, Laplace equation, Poisson equation, Wave equation and Burgers equation. The 9-bit model has been verified by several test cases. Numerical simulations, including 1D and 2D cases, of each problem are shown respectively. Comparisons are made between numerical predictions and analytic solutions or available numerical results from previous researchers. It turned out that the 9-bit model is computationally effective and accurate for all different mathematical-physical equations studied. The main benefits of the new model proposed is that it is faster than the previous existing models and has a better accuracy.Peer ReviewedPostprint (author's final draft

A mesh adaptivity scheme on the Landau-de Gennes functional minimization case in 3D, and its driving efficiency

This paper presents a 3D mesh adaptivity strategy on unstructured tetrahedral meshes by a posteriori error estimates based on metrics, studied on the case of a nonlinear finite element minimization scheme for the Landau-de Gennes free energy functional of nematic liquid crystals. Newton's iteration for tensor fields is employed with steepest descent method possibly stepping in. Aspects relating the driving of mesh adaptivity within the nonlinear scheme are considered. The algorithmic performance is found to depend on at least two factors: when to trigger each single mesh adaptation, and the precision of the correlated remeshing. Each factor is represented by a parameter, with its values possibly varying for every new mesh adaptation. We empirically show that the time of the overall algorithm convergence can vary considerably when different sequences of parameters are used, thus posing a question about optimality. The extensive testings and debugging done within this work on the simulation of systems of nematic colloids substantially contributed to the upgrade of an open source finite element-oriented programming language to its 3D meshing possibilities, as also to an outer 3D remeshing module

Structure model index does not measure rods and plates in trabecular bone

Structure model index (SMI) is widely used to measure rods and plates in trabecular bone. It exploits the change in surface curvature that occurs as a structure varies from spherical (SMI = 4), to cylindrical (SMI = 3) to planar (SMI = 0). The most important assumption underlying SMI is that the entire bone surface is convex and that the curvature differential is positive at all points on the surface. The intricate connections within the trabecular continuum suggest that a high proportion of the surface could be concave, violating the assumption of convexity and producing regions of negative differential. We implemented SMI in the BoneJ plugin and included the ability to measure the amounts of surface that increased or decreased in area after surface mesh dilation, and the ability to visualize concave and convex regions. We measured SMI and its positive (SMI+) and negative (SMI-) components, bone volume fraction (BV/TV), the fraction of the surface that is concave (CF), and mean ellipsoid factor (EF) in trabecular bone using 38 X-ray microtomography (XMT) images from a rat ovariectomy model of sex steroid rescue of bone loss, and 169 XMT images from a broad selection of 87 species' femora (mammals, birds, and a crocodile). We simulated bone resorption by eroding an image of elephant trabeculae and recording SMI and BV/TV at each erosion step. Up to 70%, and rarely less than 20%, of the trabecular surface is concave (CF 0.155 – 0.700). SMI is unavoidably influenced by aberrations from SMI-, which is strongly correlated with BV/TV and CF. The plate-to-rod transition in bone loss is an erroneous observation resulting from SMI's close and artefactual relationship with BV/TV. SMI cannot discern between the distinctive trabecular geometries typical of mammalian and avian bone, whereas EF clearly detects birds' more plate-like trabeculae. EF is free from confounding relationships with BV/TV and CF. SMI results reported in the literature should be treated with suspicion. We propose that EF should be used instead of SMI for measurements of rods and plates in trabecular bone

Wetting of anisotropic sinusoidal surfaces - experimental and numerical study of directional spreading

Directional wettability, i.e. variation of wetting properties depending on the surface orientation, can be achieved by anisotropic surface texturing. A new high precision process can produce homogeneous sinusoidal surfaces (in particular parallel grooves) at the micro-scale, with a nano-scale residual roughness five orders of magnitude smaller than the texture features. Static wetting experiments have shown that this pattern, even with a very small aspect ratio, can induce a strong variation of contact angle depending on the direction of observation. A comparison with numerical simulations (using Surface Evolver software) shows good agreement and could be used to predict the fluid-solid interaction and droplet behaviour on textured surfaces. Two primary mechanisms of directional spreading of water droplets on textured stainless steel surface have been identified. The first one is the mechanical barrier created by the textured surface peaks, this limits spreading in perpendicular direction to the surface anisotropy. The second one is the capillary action inside the sinusoidal grooves accelerating spreading along the grooves. Spreading has been shown to depend strongly on the history of wetting and internal drop dynamics