Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

NEAMS Workbench and Bison Fuel Performance Remote Application Configuration

The mission of the US Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program is to develop, apply, deploy, and support state-of-the-art predictive modeling and simulation tools for the design and analysis of current and future nuclear energy systems. This is accomplished using computing architectures that range from laptops to leadership class facilities. The NEAMS Workbench [1,2,3] initiative will facilitate the transition from conventional computer codes to high-fidelity ones by providing a common analysis environment and interface for model creation, review, execution, output review, and visualization for integrated codes. The BISON fuel performance code [4] is a finite element method (FEM)-based software being developed by multiple organizations led by the Idaho National Laboratory for the US Department of Energy Office of Nuclear Energy. BISON can simulate a variety of fuel and cladding compositions and geometries in varying dimensionalities, depending upon the problem type and complexity, including one-dimensional (1D), 1.5D, 2D axially symmetric (2D-rz), 2D radial-azimuthal (2D-rθ), or 3D models. The NEAMS Workbench allows execution of the BISON fuel performance code on local and remote compute resources. In this report, the detailed instructions of configuring the NEAMS Workbench for remote execution of the BISON fuel performance code on a portable batch system (PBS)-scheduled compute resource, as well the transfer of associated Bison files, is illustrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Debris Dispersion Model Using Java 3D

This paper describes web based simulation of Shuttle launch operations and debris dispersion. Java 3D graphics provides geometric and visual content with suitable mathematical model and behaviors of Shuttle launch. Because the model is so heterogeneous and interrelated with various factors, 3D graphics combined with physical models provides mechanisms to understand the complexity of launch and range operations. The main focus in the modeling and simulation covers orbital dynamics and range safety. Range safety areas include destruct limit lines, telemetry and tracking and population risk near range. If there is an explosion of Shuttle during launch, debris dispersion is explained. The shuttle launch and range operations in this paper are discussed based on the operations from Kennedy Space Center, Florida, USA.

Thirumalainambi, Rajkumar↗

IC reports (2021)

We requested HPC support to continue research on the seismic waves generated by impacts. We used the following codes: (1) the Hybrid Optimization Software Suite (HOSS), developed at LANL. HOSS is based on a combined Finite and Discrete Element Method (FDEM). New material models are developed for the sedimentary rocks. HOSS has been recently benchmarked to iSale and FLAG codes. SPECFEM3D is an open-source code developed since the last 90s. It won the Gordon Bell award for best performance in 2003, was finalist again in 2008 for a run at 0.16 petaflops on 149,784 cores on the ‘Jaguar’ Cray system at Oak Ridge National Laboratory. It also won the BULL Joseph Fourier supercomputing award in 2010.; SW4 is a 4th-order finite difference code developed at LLNL which is currently actively developed to handle complex 3D models and to be ported on future exascale platforms. We assessed our need to a total of 3.9M CPU-hrs for year 1 and and 3.1 M for year 2.

79 ASTRONOMY AND ASTROPHYSICS↗

Tech-X contributions to the Center for Tokamak Transients (CTTS) Scientific Discovery through Advanced Computing (SciDAC) project (Final Report)

The Center for Tokamak Transients Simulations (CTTS) consists of computational plasma physicists, computer scientists, and applied mathematicians who are together developing and applying simulation software to improve our understanding of and ability to predict and control transient events in tokamak discharges, particularly those that lead to violent disruptions. We have major simulation activities concerning disruptions that are caused by violating ideal MHD stability criteria, vertical displacement events, resistive wall modes, and tearing modes. We are performing comprehensive modeling of disruption mitigation techniques, particularly by shattered pellet injection (SPI). The Center is built around two mature extended MHD codes: NIMROD and M3D-C1. These codes are being extended in a number of ways as required by the demands of disruption modeling, and are being optimized to work well for larger problem size and modern computer architectures as required to model tokamak plasmas with realistic parameters over the required timescales. The codes will be coupled to the FronTier code for state-of-the-art modeling of SPI. Validation with experimental data on DIII-D, NSTX, and JET is being facilitated by rigorous uncertainty quantification (UQ). The high-level objectives of the center are to: (1) Develop better understanding and improved predictive capability of when crossing a linear stability boundary will likely lead to a disruption, (2) Provide a tool for calculating vessel forces, both axisymmetric and non-axisymmetric, for worst case VDE and other disruptions in tokamaks, (3) Develop increased understanding and improved predictive capability for locked mode disruptions and how best to avoid them, (4) Develop, verify, and validate 3D models for disruption mitigation by shattered pellets that can be used to design and optimize a system for future tokamaks, and (5) Increase the efficiency and scalability of the M3D-C1 and NIMROD codes so that they can model disruptions in ITER scale devices using the most powerful DOE HPC hardware, and (6) Demonstrate code-coupling via FronTier coupling with NIMROD and M3D-C1 and by interfacing these codes within a WDM framework.

97 MATHEMATICS AND COMPUTING↗

3D Radiative MHD Modeling of the Solar Atmospheric Dynamics and Structure

Dramatic dynamical phenomena accompanied by strong thermodynamic and magnetic structuring are the primary drivers of great interest in studying the solar atmosphere with high spatial and temporal resolutions. Using current computational capabilities, it has become possible to model the magnetized solar plasma in different regimes with a high degree of realism. To study the fine structuring of the solar atmosphere and dynamics, we use 3D MHD radiative models covering all layers from the upper convection zone to the corona. Realistic 3D radiative MHD modeling of the solar magnetoconvection and atmosphere allows us to generate synthetic observables that directly link the physical properties of the solar plasma to spectroscopic observables. In the presentation, we discuss qualitative and quantitative changes of the atmospheric structure and dynamics at different layers of the solar atmosphere, properties of acoustic and surface gravity waves, sources of local heating in the chromosphere-corona transition region, formation of shocks, and high-frequency oscillations in the corona, as well as manifestation of these phenomena in the modeled observables.

SMD↗

Neural Radiance Fields for Shadow Removal and Relighting of Satellite Imagery

Multi-view stereo photogrammetric techniques have been conventionally utilized to generate Global Digital Elevation Models (GDEM) of planetary and lunar surfaces. However, these methods, relying on conventional feature detectors, are often subject to inaccuracies caused by changes in lighting conditions, including diffuse reflection and harsh shading. This has limited the ability of these methods to accurately reconstruct shadowed regions in orbital imagery, such as highly shaded urban areas and the permanently shadowed regions (PSRs) located in the lunar surface, which are critical targets for NASA's Artemis program. Neural Radiance Fields (NeRFs) offer a novel solution to these limitations by breaking away from traditional photogrammetric assumptions of ridged, opaque surfaces. NeRFs are capable of reconstructing 3D objects with variably transmissive properties, such as fluids, PSRs, and reflective surfaces. In this work, we demonstrate the robustness of NeRFs in generating high-fidelity 3D models of terrain from highly shaded orbital imagery acquired from small satellites in low earth orbit (LEO). Our results showcase the ability of NeRF to overcome the limitations of traditional photogrammetric methods and provide a promising solution for reconstructing complex scenes in challenging lighting conditions.

Caleb Ashmore Adams↗

Benchmark Simulation of Natural Circulation Cooling System with Salt Working Fluid Using SAM

Liquid salt-cooled reactors, such as the Fluoride Salt-Cooled High-Temperature Reactor (FHR), offer passive decay heat removal through natural circulation using Direct Reactor Auxiliary Cooling System (DRACS) loops. The behavior of such systems should be well-understood through performance analysis. The advanced system thermal-hydraulics tool System Analysis Module (SAM) from Argonne National Laboratory has been selected for this purpose. The work presented here is part of a larger study in which SAM modeling capabilities are being enhanced for the system analyses of FHR or Molten Salt Reactors (MSR). Liquid salt thermophysical properties have been implemented in SAM, as well as properties of Dowtherm A, which is used as a simulant fluid for scaled experiments, for future code validation studies. Additional physics modules to represent phenomena specific to salt-cooled reactors, such as freezing of coolant, are being implemented in SAM. This study presents a useful first benchmark for the applicability of SAM to liquid salt-cooled reactors: it provides steady-state and transient comparisons for a salt reactor system. A RELAP5-3D model of the Mark-1 Pebble-Bed FHR (Mk1 PB-FHR), and in particular its DRACS loop for emergency heat removal, provides steady state and transient results for flow rates and temperatures in the system that are used here for code-to-code comparison with SAM. The transient studied is a loss of forced circulation with SCRAM event. To the knowledge of the authors, this is the first application of SAM to FHR or any other molten salt reactors. While building these models in SAM, any gaps in the code’s capability to simulate such systems are identified and addressed immediately, or listed as future improvements to the code.

Ahmed, K. K.↗

Inhibition of TGF-β pathway reverts extracellular matrix remodeling in T. cruzi-infected cardiac spheroids

Highlights: • A model of cardiac spheroids to study fibrosis in pathological processes is proposed. • The study evaluates a TGF-β inhibitor in restoring extracellular matrix balance disrupted by T. cruzi infection. • Inhibition of TGF-β signaling reduces fibrosis and hypertrophy in infected-cardiac spheroids. • The extracellular matrix remodeling is achieved through regulation of MMP-2 activity and TIMP-1 expression. Chagasic cardiomyopathy (CC) is the main manifestation of Chagas Disease (CD). CC is a progressive dysfunctional illness, in which transforming growth factor beta (TGF-β) plays a central role in fibrogenesis and hypertrophy. In the present study, we tested in a three-dimensional (3D) model of cardiac cells culture (named cardiac spheroids), capable of mimicking the aspects of fibrosis and hypertrophy observed in CC, the role of TGF-β pathway inhibition in restoring extracellular matrix (ECM) balance disrupted by T. cruzi infection. Treatment of T. cruzi-infected cardiac spheroids with SB 431542, a selective inhibitor of TGF-β type I receptor, resulted in a reduction in the size of spheroids, which was accompanied by a decrease in parasite load and in fibronectin expression. The inhibition of TGF-β pathway also promoted an increase in the activity of matrix metalloproteinase (MMP)-2 and a decrease in tissue inhibitor of matrix metalloproteinase (TIMP)-1 expression, which may be one of the mechanisms regulating extracellular matrix remodeling. Therefore, our study provides new insights into the molecular mechanisms by which inhibition of TGF-β signaling reverts fibrosis and hypertrophy generated by T. cruzi during CC and also highlights the use of cardiac spheroids as a valuable tool for the study of fibrogenesis and anti-fibrotic compounds.

60 APPLIED LIFE SCIENCES↗

Flow boiling transient critical heat flux tests with stainless steel and FeCrAl: Transient correlation implementation, model calibration, and sensitivity analysis

In this study, steady-state and transient internal-flow critical heat flux (CHF) experiments were carried out under two flow conditions, at atmospheric pressure, on 316 stainless steel (316-SS) and iron-chromium-aluminum (FeCrAl) tubes. Slow power transients at a low mass flow (300 kg/m 2 - s) with very low subcooling (X e =-0.0054) generated premature CHFs, which were prevented by faster power transients and high mass flow (1000 kg/m 2 - s). The measured transient CHFs increased linearly with increasing power transients, compared with the steady-state CHF. Yet, the wall superheat at the CHFs decreased with faster power transients. Transient CHF correlations highlighting heterogeneous spontaneous nucleation were calibrated to the measured CHFs and compared with other existing correlations. Transient CHF multipliers were acquired from pool and flow boiling empirical CHF correlations that were generated. The multipliers were applied to the RELAP5-3D nuclear system code to analyze the discrepancies between the measured data and the predicted CHF and post-CHF behavior, which improved peak cladding temperature predictions by 24.6%. A variance-based global sensitivity study perturbing the experimental uncertainties and cladding material thermal properties showed the diminishing influence of flow boiling heat transfer with increasing power transients, highlighting the significance of the volumetric heat capacity for the cladding integrity during transients. Transient CHF correlations were applied to the most limiting design basis accident: a hot-zero-power reactivity-initiated accident on a generic pressurized water reactor RELAP5-3D model.

36 MATERIALS SCIENCE↗

A mechanistic model of a PWR-based nuclear power plant in response to external hazard-induced station blackout accidents

Natural hazard-induced nuclear accidents, such as the Fukushima Daiichi Accident that occurred in Japan in 2011, have significantly increased reactor safety studies in understanding nuclear power plant (NPP) responses to external hazard events such as earthquakes and floods. Natural hazards could cause the loss of offsite power in nuclear power plants, potentially leading to a Station Blackout (SBO) accident that significantly contributes to the overall risk of nuclear power plant accidents. Despite the fact that extensive research has been conducted on the station blackout accident for nuclear power plant, further understanding of these events is needed, particularly in the context of the dynamic nature of external hazards such as external flooding. This paper estimates the progression of station blackout events for a generic pressurized water reactor (PWR) in response to external flooding events. The original RELAP5-3D model of the Westinghouse four-loop design pressurized water reactor was adopted and modified to simulate the external flood-induced station blackout accident, including the short-term and long-term station blackout scenarios. A sensitivity analysis of long-term station blackout, examining reactor operation times and analyzing key parameters over time, was also conducted in this work. The results of the analyses, especially the critical timing parameters of key event sequences, provide useful insights about the time during the external flooding event, which is important for plant operators to make timely decisions to prevent potential core damage. This paper represents significant progress toward developing an integrated risk assessment framework for further identifying and assessing the effects of the critical sources of uncertainties of nuclear power plant under external hazard-induced events.

Liu, Tao↗

Vadose and Saturated Zone Flow and Transport Calculations for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

The purpose of the fate and transport modeling described in this environmental calculation file (ECF) is to evaluate the impacts to groundwater associated with waste disposal operations at Mixed Waste Low-Level Burial Ground (LLBG) Trenches 31 and 34 to satisfy requirements in DOE O 435.1, Radioactive Waste Management. The model integrates the flow and transport in the vadose zone beneath the active trenches with the saturated zone downgradient of the trenches to predict the radionuclide concentration at the point of assessment (POA). DOE M 435.1-1, Radioactive Waste Management Manual, defines the POA as the point of highest projected dose or concentration beyond a 100 m (328 ft) buffer zone surrounding the disposed waste. The modeling is conducted in accordance with the DOE G 435.1, Implementation Guide for Use with DOE M 435.1-1, performance assessment (PA) guidelines. The modeling involves evaluation of the groundwater concentrations and radionuclide arrival times during the 1,000-year compliance and 10,000-year sensitivity-uncertainty periods per DOE O 435.1 and DOE M 435.1. This analysis does not consider radionuclide release during facility operations, only the post-closure impacts of the radionuclides to the environment. The evaluation of potential radiological dose to groundwater receptors caused by releases from a closed facility containing radioactive waste typically includes the following: (1) Release of radionuclides from that facility (2) Transport of those radionuclides through the environment, and (3) Exposure to humans to environmental concentration levels of those radionuclides The fate and transport three-dimensional (3D) model analysis involves the post-closure impacts to the environment of the technetium-99, iodone-129, and uranium (all isotopes in the waste). The residual inventory estimates include several radionuclides, but technetium-99 is typically responsible for almost all of the beta-gamma dose equivalent associated with groundwater (water resources) protection per 40 CFR 141, “National Primary Drinking Water Regulations” (e.g., see the results in WCH-520, Performance Assessment of Environmental Restoration Disposal Facility, Hanford Site, Washington; hereinafter referred to as the ERDF PA), and iodone-129 can also be a significant dose contributor for some waste (e.g., RPP-RPT-59958, Performance Assessment for the Integrated Disposal Facility, Hanford Site, Washington; hereinafter referred to as the IDF PA). Uranium does not typically factor significantly into the impacts to groundwater, even during the 10,000-year sensitivity-uncertainty period, but always remains of interest as a contaminant. This ECF does not address vadose and saturated zone modeling for Trench 94 of the 200 East Area LLBG. Current information confirms the validity of the low corrosion rate of the naval reactor plant carbon steel (HY-80), and the even lower corrosion rate of the nickel-iron-chromium alloy reactor vessel (Inconel Alloy 600) presented in DOE/EIS-0259, Final Environmental Impact Statement on the Disposal of Decommissioned, Defueled Cruiser, Ohio Class, and Los Angeles Class Naval Reactor Plants. Based on these low corrosion rates, the time to breach the reactor vessel to allow release of radionuclides from the activated metal of the reactor vessel internal structure is at least 10,000 years. This time to breach precludes the need to evaluate the vadose and saturated zone transport of contaminants released from the reactor compartment disposal packages in the 200 East Area LLBGs PA (CP-63826, Waste Release Model Package Report for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Representativity error scaling of models of the high temperature test facility

Error scaling is a critical step toward the validation of modeling capabilities for high-temperature gas-cooled reactors (HTGRs). Extensive effort is being made to bring HTGRs into the validation basis of numerous thermal-hydraulics codes. Here, this paper demonstrates how systems-level codes can be leveraged to perform error scaling analyses between an experimental facility and a plant-to-be facility. Specifically, we focus on two conduction cooldown experiments from the High Temperature Test Facility (HTTF) and the General Atomics 350 MW th modular high-temperature gas-cooled reactor (MHTGR-350). The error scaling methodology employed in this study is representativity, which in the context of this work was used to quantify how well experiments captured the physics of the plant facility by comparing sensitivity vectors between the two facilities. In addition, two different RELAP5–3D models of the HTTF were included in the comparison to determine if modeling methodologies notably impact the representativity results. The key figures of merit are the maximum block temperature and the coolant outlet temperature. The time-dependent maximum block temperature had a low representativity of below 0.1 for both models across both transients. The coolant outlet temperature had a higher representativity of around 0.6 for both models during the pressurized conduction cooldown, but it was below 0.2 for the depressurized conduction cooldown. Overall, the HTTF experiments were not representative of the transients in the MHTGR-350. However, these results can play a significant role in informing future potential experiments for HTGR systems that iterate on what the HTTF accomplished.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

A generalized massively parallel ultra-high order FFT-based Maxwell solver

Dispersion-free ultra-high order FFT-based Maxwell solvers have recently proven to be paramount to a large range of applications, including the high-fidelity modeling of high-intensity laser–matter interactions with Particle-In-Cell (PIC) codes. To enable a massively parallel scaling of these solvers, a novel parallelization technique was recently proposed, which consists in splitting the simulation domain into several processor sub-domains, with guard regions appended at each sub-domain boundary. Maxwell's equations are advanced independently on each sub-domain using local shared-memory FFTs (instead of a single distributed global FFT). This implies small truncation errors at sub-domain boundaries, the amplitude of which depends on guard regions sizes and order of the Maxwell solver. For moderate guard region sizes, this ’local’ technique proved to be highly scalable on up to a million cores and notably enabled the 3D modeling of so-called plasma mirrors, for which 8 guard cells only were enough to prevent truncation error growth. Yet, for other applications, the required number of guard cells might be much higher, which would severely limit the parallel efficiency of this technique due to the large volume of guard cells to be exchanged between sub-domains. In this context, we propose a novel parallelization technique that ensures very good scaling of FFT-based solvers with an arbitrarily high number of guard cells. Our ’hybrid’ technique consists in performing distributed FFTs on local groups of processors with guard regions now appended to boundaries of each group of processors. It uses a dual domain decomposition method for the Maxwell solver and other parts of the PIC cycle to keep the simulation load-balanced. This ’hybrid’ technique was implemented in the open source exascale library PICSAR. Benchmarks show that for a large number of guard cells (>16), the ’hybrid’ technique offers up to ×3 speed-up and ×8 memory savings compared to the ’local’ one.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Casing-integrity mapping using top-casing electrodes and surface based EM fields

Wellbore integrity is of paramount importance to subsurface resource extraction, energy storage and hazardous waste disposal. We introduce a simple non-invasive technology for casing integrity screening, based on the continuity of electrical current flow. Applying low frequency current to a wellhead, with a distant return electrode, produces a casing current dependent on the properties and depth extent of the well casing as well as the background formation. These currents in-turn generate surface electrical fields that can be captured in a radial profile and be used to analyze properties of the well casing. Numerical modeling results reveal a strong relation of the electric field to the casing properties and depth extent of the well. A small breakage in the casing produces a profile coincident to a cased well with a completion depth above the break. A corroded patch, where the casing conductivity is reduced, also alters the field profiles and its depth may be estimated by comparing to the profile expected from the well completion diagrams. The electric field profiles are also strongly dependent on background resistivity distributions and on whether the well was drilled using water or oil-based drilling fluids. We validate the proof of concept in a field experiment, where we applied currents at the wellheads of two wells with different casing lengths. The two profiles were similar in appearance but offset in amplitude by more than a factor of 5, consistent with the theoretical analysis as well as the 3D modeling results. These results demonstrate that our proposed approach has promise for mapping the general casing condition without well intervention. This approach can be a practical and effective tool for rapidly screening a number of wells before expensive logging-based technologies are employed for casing inspection in detail.

58 GEOSCIENCES↗

DIMENSIONAL DEPENDENCE OF THE HYDRODYNAMICS OF CORE-COLLAPSE SUPERNOVAE

A major goal over the last decade has been understanding which multidimensional effects are crucial in facilitating core-collapse supernova (CCSN) explosions. Unfortunately, much of this work has necessarily assumed axisymmetry. In this work, we present analyses of simplified two-dimensional (2D) and three-dimensional (3D) CCSN models with the goal of comparing the hydrodynamics in setups that differ only in dimension. Not surprisingly, we find many differences between 2D and 3D models. While some differences are subtle and perhaps not crucial, others are dramatic and make interpreting 2D models problematic. In particular, axisymmetric models produce excess power at the largest spatial scales, power that has been deemed critical in previous explosion models. Nevertheless, our 3D models, which have an order of magnitude less power than 2D models on large scales, explode earlier. Since explosions occur earlier in 3D than in 2D, the vigorous large-scale sloshing is either not critical in any dimension or the explosion mechanism operates differently in 2D and 3D. On the other hand, we find that the average parcel of matter in the gain region has been exposed to net heating for up to 30% longer in 3D than in 2D, an effect we attribute to the differing characters of turbulence in 2D and 3D. We suggest that this effect plays a prominent role in producing earlier explosions in 3D. Finally, we discuss a simple model for the runaway growth of buoyant bubbles that is able to quantitatively account for the growth of the shock radius and predicts a critical luminosity relation.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of trophoblast cystic structures from human induced pluripotent stem cells in limited-area cell culture

Highlights: • Limited-area culture of iPSs induced trophectoderm-like cyst structures. • Trophectoderm-like cyst structures exhibited trophoblast function. • These cysts provide a 3D model for the study of differentiation and morphogenesis. We developed a novel engineering technique to induce differentiation of human induced pluripotent stem cells (hiPSCs) into organoids mimicking the trophectoderm (TE). Here, hiPSCs were cultured on a limited area of 2–4 mm in diameter. After 15–20 days, spherical cysts appeared on the surface of the limited area. Secretion of human chorionic gonadotrophin (hCG) began to increase after ∼ 20 days and remained dramatically elevated over the next 20 days. Limited-area-cultured cysts exhibited expression of hCG, which was a result of epithelial differentiation. Low expression levels of pluripotent genes and high expression levels of trophoblast lineage-specific genes were detected in the cells of spherical cysts. Multinucleated syncytia trophoblast was observed in the reseeded cystic cells. We observed hiPSC-derived cysts that morphologically resembled trophectoderm in vivo. The limited-area cell culture induced a three-dimensional (3D) trophectoderm organoid, which has potential for use in the study of human trophoblast differentiation and placental morphogenesis.

60 APPLIED LIFE SCIENCES↗

Tightly-coupled camera/LiDAR integration for point cloud generation from GNSS/INS-assisted UAV mapping systems

Unmanned aerial vehicles (UAVs) equipped with integrated global navigation satellite systems/inertial navigation systems (GNSS/INS) together with cameras and/or LiDAR sensors are being widely used for topographic mapping in a variety of applications such as precision agriculture, coastal monitoring, and archaeological documentation. Integration of image-based and LiDAR point clouds can provide a comprehensive 3D model of the area of interest. For such integration, ensuring a good alignment between data from the different sources is critical. Although many works have been conducted on this topic, there is still a need for a rigorous integration approach that minimizes the discrepancy between camera and LiDAR data caused by inaccurate system calibration parameters and/or trajectory artifacts. This study proposes an automated tightly-coupled camera/LiDAR integration workflow for GNSS/INS-assisted UAV systems. The proposed strategy is conducted in three main steps. First, an image-based point cloud is generated using a LiDAR/GNSS/INS-assisted structure from motion (SfM) strategy. Then, feature correspondences between image-based and LiDAR point clouds are automatically identified. Finally, an integrated-bundle adjustment procedure including image points, LiDAR raw measurements, and GNSS/INS information is conducted to minimize the discrepancy between point clouds from different sensors while estimating system calibration parameters and refining the trajectory information. The proposed SfM strategy and integration framework are evaluated using five datasets. The SfM results show that using LiDAR data can facilitate feature matching and further increase the number of reconstructed 3D points. The experimental results also illustrate that the developed automated camera/LiDAR integration strategy is capable of accurately estimating system calibration parameters to achieve good alignment among camera/LiDAR data from single/multiple systems. Finally, an absolute accuracy in the range of 3–5 cm is achieved for the image/LiDAR point clouds after the integration process.

42 ENGINEERING↗

Advancing Understanding of the Role of Belowground Processes in Terrestrial Carbon Sinks trhrough Ground-Penetrating Radar. Final Report

Coarse roots play a significant role in belowground carbon cycling and will likely play an increasingly crucial role in belowground carbon sequestration as atmospheric CO 2 levels continue to rise, yet they are one of the most difficult ecosystem parameters to quantify. Despite promising results with ground-penetrating radar (GPR) as a nondestructive method of quantifying biomass of coarse roots, this application of GPR is in its infancy and neither the complete potential nor limitations of the technology have been fully evaluated. The primary goals and questions of this study fell into four groups: (1) GPR methods: Can GPR detect change in root biomass over time, differentiate live roots from dead roots, differentiate between coarse roots, fine roots bundled together, and a fine root mat, remain effective with varied soil moisture, and detect shadowed roots (roots hidden below larger roots); (2) CO 2 enrichment study at Kennedy Space Center in Brevard County, Florida: Are there post-fire legacy effects of CO 2 fertilization on plant carbon pools following the end of CO 2 application ? (3) Disney Wilderness Study: What is the overall coarse root biomass and potential for belowground carbon storage in a restored longleaf pine flatwoods system? Can GPR effectively quantify coarse roots in soils that are wetter than the previous sites and that have a high percentage of saw palmetto rhizomes present? (4) Can GPR accurately represent root architecture in a three-dimensional model? When the user is familiar with the equipment and software in a setting that minimizes unsuitable conditions, GPR is a relatively precise, non-destructive, useful tool for estimating coarse root biomass. However, there are a number of cautions and guidelines that should be followed to minimize inaccuracies or situations that are untenable for GPR use. GPR appears to be precise as it routinely predicts highly similar values for a given area across multiple scanning events; however, it appears to lack sufficient accuracy at small scales. Knowledge of soil conditions and their effects on GPR wave propagation and reception are paramount for the collection of useful data. Strong familiarity with the software and equipment is both important and necessary for GPR use in estimating coarse root biomass. GPR must be utilized at low soil moisture levels in order to accurately represent existing coarse root structures. Our results from Disney Wilderness Preserve highlight the need for a strong understanding of the limitations of GPR, specifically knowledge of root structures (saw palmetto rhizomes) or environmental factors (low moisture content) that may hinder its application within a given system. The 3D modeling of course roots with GPR appears quite promising, as it has become more accurate and precise as the software has advanced and become more robust, but there is still a need for more precision before it will likely be able to model anything more than simple root systems comprised mostly of large diameter roots. Our results from Kennedy Space Center suggest that there are legacy effects from CO 2 fertilization in the form of more root mass providing a greater capacity for aboveground plant regrowth following fire, even 7 years after treatment ended.

54 ENVIRONMENTAL SCIENCES↗