Search NASA⌕ Search

SEARCH · Search NASA

Results for “COUPLING NETWORK”

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 1,045 records · Page 58

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

15 GEOTHERMAL ENERGY↗

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗

Integrated Commercial Carbon Capture and Storage (CCS) Prefeasibility Study at Rock Springs Uplift, Wyoming (Final Report)

To establish the potential of safe, long-term commercial-scale carbon capture & storage (CCS) in southwestern Wyoming, the University of Wyoming formed a CCS coordination team (CCT) to complete a Phase I CarbonSAFE pre-feasibility study adjacent to the region’s largest coal-fired power plant, the Jim Bridger Plant (JBP). The study’s primary objective was to evaluate the potential of storing 50+ million metric tons of anthropogenic carbon dioxide (CO 2 ) in stacked reservoirs over a period of 25 years. In addition, this integrated study evaluated related aspects of carbon capture utilization & storage (CCUS) at the study site, including: (1) the CO 2 character of post-combustion flue gas from JBP and the technical capacity for retrofitting proven commercial-scale CO 2 capture technology relative to CO 2 source character and technological constraints; (2) utilization of the existing CO 2 pipeline network in the immediate vicinity of JBP for enhanced oil recovery (CO 2 -EOR) and utilizing that network as part of a regional hub for CO 2 storage from other anthropogenic sources; (3) assessment of the challenges and benefits of meeting all of CarbonSAFE programmatic Phase I goals relative to Wyoming’s carbon management regulatory framework; and (4) assessment of the community and environmental factors that may impact CCS at JBP. With respect to storage, we analyzed the potential of stacked Mesozoic reservoirs at the study site, focusing on the Entrada and Nugget sandstones. The reservoir response to injection simulations varied greatly between formations, though they both were capable of CO 2 storage. The Entrada Sandstone is relative thin (~40’) and heterogeneous; storage using one injection and production well was limited to 1.2 Mt over 25 years at the study site. The Nugget Sandstone is much thicker (> 400’) and generally homogeneous; storage using one injection and production well approached 15.0 Mt over 25 years; this scenario suggests that the storage potential in the Nugget Sandstone is close to the estimated P90 calculation of 9.6 Mt/mi 2 as models indicate a storage potential of 8.3 Mt/mi 2 at saturation. Both reservoirs’ storage capacity benefited greatly from coupling an injection well with an offset production well; total storage capacity increased and reservoir pressure stayed below critical thresholds. This indicates that implementing a pressure management strategy at the study site would optimize storage and decrease risk. Seals associated with the targeted reservoirs were shown to be capable of retaining the potential storage volumes, and the site benefits from multiple seal redundancies (i.e. > 7,000’ of seal formations). Furthermore, assessment of different reservoir variables (e.g., fluid composition, pressure, etc.) indicates all reservoirs at the study site are isolated. Storage assessments from Mesozoic reservoirs were coupled with simulations from the site’s Paleozoic reservoirs (the Weber Sandstone and Madison Limestone), which were the focus of previous study. The fully integrated stacked reservoir pre-feasibility studies suggest that the Nugget and Madison formations alone could safely store 50 million metric tons of CO 2 within the study site. Furthermore, both formations are true saline reservoirs near JBP, making them ideal target reservoirs for CO 2 sequestration. Analysis of Wyoming’s laws and regulations confirm that the State has a favorable policy environment for CCS/CCUS. The project’s economic models estimated that the following revenues collectively are sufficient to finance the project over its lifespan if CO 2 utilization is coupled with CCS via the following strategies: (1) sales of CO 2 for EOR (approximately $69 million/year); (2) use of CO 2 tax credits such as amended §45Q and §48A (approximately $484 million total); and (3) sales of low-carbon electricity and marketable carbon offset/credits into carbon-constrained markets (approximately $11-$17 million/year) to West Coast consumers. These revenue estimates are broadly consistent with other studies that concluded arbitrage in JBP retail electricity sales between California and Wyoming could help to support the cost of deploying CCS at JBP. The models also indicate that maintaining a successful business case for long-term CCS -- wherein associated capture, transportation and operational costs are fully addressed at JBP -- is reliant on numerous external factors that complicate predictions. This study suggests that this site meets the CarbonSAFE program’s requirements of being able to feasibly store 50+ million metric tons of CO 2 over 25 years within the site’s stacked reservoirs, especially if coupled with pressure management strategies. This study suggests that the Nugget Sandstone and the Madison Limestone are two of the best reservoirs in the State with respect to overall storage capacity. In addition, the site’s proximity to existing CO 2 transportation networks, CO 2 -EOR opportunities, ability to sell to markets that value low-carbon electricity, Wyoming’s existing carbon regulatory framework, and a public that is well-educated with respect to energy markets and issues collectively increase the potential for implementing commercial-scale CCS adjacent to the JBP. With respect to the available technical and non-technical site data, and the reservoir quality at the study site, we suggest that this site is capable of moving directly into Phase III of the CarbonSAFE program.

01 COAL, LIGNITE, AND PEAT↗

Model Data Archive Associated with Manuscript "Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon"

This data package supports the publication “Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon” by Li et al. (2026). The package contains processed model inputs, configuration files, restart files, simulation outputs, scripts, and visualization products used to evaluate post-fire dissolved organic carbon (DOC) dynamics in the Naches River Watershed, Washington, USA, following the 2021 Schneider Springs Fire. The modeling workflow couples ELM-BGC, the biogeochemistry-enabled Energy Exascale Earth System Model Land Model; ATS, the Advanced Terrestrial Simulator for integrated surface-subsurface hydrology; and PFLOTRAN, a reactive transport model for multicomponent aqueous geochemistry. Together, these models simulate how wildfire-induced changes in vegetation, litter, coarse woody debris, and soil organic matter influence DOC production, transport, and reaction from burned hillslopes to stream networks. The archive includes preprocessed meteorological, geospatial, hydrologic, and biogeochemical forcing data; ELM-BGC-derived DOC source terms; ATS mesh files; PFLOTRAN reactive-transport inputs; model configuration files; spin-up and transient restart files; watershed-scale diagnostic outputs; stream concentration time series; and figures or visualization files used to inspect and reproduce key results. File types include Hierarchical Data Format 5 (HDF5) files for gridded forcing and model-coupling data, model input and configuration files for ELM-BGC, ATS, and PFLOTRAN, restart and simulation-output files generated by the modeling workflow, tabular or time-series diagnostic outputs, scripts for post-processing and figure generation, and image or visualization products associated with the manuscript. Use of the package depends on the intended task. Re-running the simulations requires the relevant modeling software, including ELM-BGC, ATS, and PFLOTRAN as ATS's geochemical engine. Inspecting outputs and reproducing figures requires Python with scientific plotting libraries such as Matplotlib, and three-dimensional model outputs may be viewed with ParaView. Geographic information system files or maps may be inspected with ArcGIS Pro or comparable GIS software. The data package is intended to enable traceability, reuse, and partial reproduction of the coupled land-to-watershed hydro-biogeochemical modeling workflow used to test how wildfire disturbance affects terrestrial carbon pools and downstream DOC dynamics.

ATS↗

Run Time Improvement Efforts for the Roman Space Telescope Thermal Analysis

"Observatory thermal models for large, complex missions, such as the Wide Field InfraRed Survey Telescope (WFIRST) mission, produce an immense amount of data to be processed. Configuration management of the model throughout the project life cycle has mainly focused on which versions of the subsystem models form the current observatory level configuration. However, the results produced by the model are not nearly as well The Roman Space Telescope (RST), formerly known as the Wide Field InfraRed Survey Telescope, is the next great astrophysics observatory mission to follow the James Webb Space Telescope with a planned launch in 2026. As a large scale, flagship mission for NASA with challenging wave front error stability requirements, a single model approach for both thermal discipline analysis and thermo-optical distortion analysis has been used since the early days of the project. In alleviating the need to maintain two separate models for different analysis types, it imposes run time penalties on the thermal analysis with a large model with significant radiation heat exchange. Throughout the lifecycle of the project, the component models have steadily grown in size, resulting in a continuous growth of the overall observatory model with each update and consequently a considerable increase in the model run time. While ongoing efforts to reduce run time are continuously investigated, previous efforts had primarily focused on timestep size and total simulation time to reach quasi-equilibrium. More recently, studies were performed on the total number of radiation couplings (radks) included in the model, which has a nearly linear impact on run time, but increases exponentially with node count. As standard practice for spacecraft analysis, small radks were excluded from the temperature solution based on the assumption that their interchange/view factors have a negligible impact on heat flow. Four approaches were investigated to reduce the model run time while minimizing the impact on accuracy: (1) the Equivalent Radiation Network node, (2) Progressive Radk Inclusion as solution proceeds, (3) Targeted Radk Filtering for critical/non critical areas, and lastly (4) Representation of culled radks with Backloads. Furthermore, the investigation of model run time also revealed that cold cases took noticeably longer to run than hot cases; the root computational inefficiencies were explored along with the computation penalty of linearization of the external radks and recalculation of temperature dependent linear couplings at each timestep. This paper outlines the details of each of the above approaches and their impact on run time and model accuracy.

Thermal Analysis↗

Run Time Improvement Efforts for the Roman Space Telescope Thermal Analysis

Observatory thermal models for large, complex missions, such as the Wide Field InfraRed Survey Telescope (WFIRST) mission, produce an immense amount of data to be processed. Configuration management of the model throughout the project life cycle has mainly focused on which versions of the subsystem models form the current observatory level configuration. However, the results produced by the model are not nearly as well The Roman Space Telescope (RST), formerly known as the Wide Field InfraRed Survey Telescope, is the next great astrophysics observatory mission to follow the James Webb Space Telescope with a planned launch in 2026. As a large scale, flagship mission for NASA with challenging wave front error stability requirements, a single model approach for both thermal discipline analysis and thermo-optical distortion analysis has been used since the early days of the project. In alleviating the need to maintain two separate models for different analysis types, it imposes run time penalties on the thermal analysis with a large model with significant radiation heat exchange. Throughout the lifecycle of the project, the component models have steadily grown in size, resulting in a continuous growth of the overall observatory model with each update and consequently a considerable increase in the model run time. While ongoing efforts to reduce run time are continuously investigated, previous efforts had primarily focused on timestep size and total simulation time to reach quasi-equilibrium. More recently, studies were performed on the total number of radiation couplings (radks) included in the model, which has a nearly linear impact on run time, but increases exponentially with node count. As standard practice for spacecraft analysis, small radks were excluded from the temperature solution based on the assumption that their interchange/view factors have a negligible impact on heat flow. Four approaches were investigated to reduce the model run time while minimizing the impact on accuracy: (1) the Equivalent Radiation Network node, (2) Progressive Radk Inclusion as solution proceeds, (3) Targeted Radk Filtering for critical/non critical areas, and lastly (4) Representation of culled radks with Backloads. Furthermore, the investigation of model run time also revealed that cold cases took noticeably longer to run than hot cases; the root computational inefficiencies were explored along with the computation penalty of linearization of the external radks and recalculation of temperature dependent linear couplings at each timestep. This paper outlines the details of each of the above approaches and their impact on run time and model accuracy.

Thermal Analysis↗

Neural Network Reflectance Prediction Model for Both Open Ocean and Coastal Waters

Remote sensing of global ocean color is a valuable tool for understanding the ecology and biogeochemistry of the worlds oceans, and provides critical input to our knowledge of the global carbon cycle and the impacts of climate change. Ocean polarized reflectance contains information about the constituents of the upper ocean euphotic zone, such as colored dissolved organic matter (CDOM), sediments, phytoplankton, and pollutants. In order to retrieve the information on these constituents, remote sensing algorithms typically rely on radiative transfer models to interpret water color or remote-sensing reflectance; however, this can be resource-prohibitive for operational use due to the extensive CPU time involved in radiative transfer solutions. In this work, we report a fast model based on machine learning techniques, called Neural Network Reflectance Prediction Model (NNRPM), which can be used to predict ocean bidirectional polarized reflectance given inherent optical properties of ocean waters. This supervised model is trained using a large volume of data derived from radiative transfer simulations for coupled atmosphere and ocean systems using the successive order of scattering technique (SOS-CAOS). The performance of the model is validated against another large independent test dataset generated from SOS-CAOS. The model is able to predict both polarized and unpolarized reflectances with an absolute error (AE) less than 0.004 for 99% of test cases. We have also shown that the degree of linear polarization (DoLP) for unpolarized incident light can be predicted with an AE less than 0.002 for 99% of test cases. In general, the simulation time of SOS-CAOS depends on optical depth, and required accuracy. When comparing the average speeds of the NNRPM against the SOS-CAOS model for the same parameters, we see that the NNRPM is able to predict the Ocean BRDF 6000 times faster than SOS-CAOS. Both ultraviolet and visible wavelengths are included in the model to help differentiate between dissolved organic material and chlorophyll in the study of the open ocean and the coastal zone. The incorporation of this model into the retrieval algorithm will make the retrieval process more efficient, and thus applicable for operational use with global satellite observations.

radiative transfer↗

A multifidelity deep operator network approach to closure for multiscale systems

Projection-based reduced order models (PROMs) have shown promise in representing the behavior of multiscale systems using a small set of generalized (or latent) variables. Despite their success, PROMs can be susceptible to inaccuracies, even instabilities, due to the improper accounting of the interaction between the resolved and unresolved scales of the multiscale system (known as the closure problem). In the current work, we interpret closure as a multifidelity problem and use a multifidelity deep operator network (DeepONet) framework to address it. In addition, to enhance the stability and/or accuracy of the multifidelity-based closure, we employ the recently developed ``in-the-loop" training approach from the literature on coupling physics and machine learning models. The resulting approach is tested on shock advection for the one-dimensional viscous Burgers equation and vortex merging for the two-dimensional Navier-Stokes equations. Furthermore, the numerical experiments show significant improvement of the predictive ability of the closure-corrected PROM over the un-corrected one both in the interpolative and the extrapolative regimes.

42 ENGINEERING↗

Coupled cluster and dislocation dynamics modeling of microstructure evolution in irradiated materials

We develop here a coupled cluster and dislocation dynamics framework to study the microstructure evolution of irradiated materials. The framework not only accounts for the three dimensional diffusion of radiation-generated clusters, but also their interaction with dislocation networks and the resultant climb motion of discrete dislocations within finite crystals. The framework is solved with a superposition solution scheme, and is applied to investigate the evolution of the irradiation-induced dislocation loops in zirconium (Zr), considering the effects of various bias factors including the diffusion anisotropy difference (DAD) of interstitials and interstitial clusters, the dislocation bias of defects to discrete dislocation segments, and the production bias of defects from the radiation cascade. We find that the DAD is the most critical factor influencing the kinetics of the loop evolution in Zr, while the recombination/interaction of mobile defects can induce a strong spatial dependence of the loop evolution together with the DAD. Here, the method is also adopted to study the evolution of interstitial $\langle$a$\rangle$ and vacancy $\langle$c$\rangle$ dislocation loop ensembles consistent with the microstructure observed during irradiation-induced growth of Zr. Our findings not only reveal the spatial dependence of the size and ellipticity of the dislocation loops, but also suggest a limit on the anisotropy factor of interstitials to reproduce the co-growth of $\langle$a$\rangle$ and $\langle$c$\rangle$ loops in zirconium, in good agreement with experimental observations and other simulation results.

Bias factors↗

Capturing intrinsic site-dependent spectral signatures and lifetimes of isolated OH oscillators in extended water networks

The extremely broad infrared spectrum of water in the OH stretching region is a manifestation of how profoundly a water molecule is distorted when embedded in its extended hydrogen-bonding network. Many effects contribute to this breadth in solution at room temperature, which raises the question as to what the spectrum of a single OH oscillator would be in the absence of thermal fluctuations and coupling to nearby OH groups. We report the intrinsic spectral responses of isolated OH oscillators embedded in two cold (~20 K), hydrogen-bonded water cages adopted by the Cs + ·(HDO)(D 2 O) 19 and D 3 O + ·(HDO)(D 2 O) 19 clusters. Most OH oscillators yield single, isolated features that occur with linewidths that increase approximately linearly with their redshifts. Oscillators near 3,400 cm -1 , however, occur with a second feature, which indicates that OH stretch excitation of these molecules drives low-frequency, phonon-type motions of the cage. Finally, the excited state lifetimes inferred from the broadening are considered in the context of fluctuations in the local electric fields that are available even at low temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical controls on the propagation rate of fracture in calcite

Calcite (CaCO 3 ) is one of the most abundant minerals in the Earth’s crust, and it is susceptible to subcritical chemically-driven fracturing. Understanding chemical processes at individual fracture tips, and how they control the development of fractures and fracture networks in the subsurface, is critical for carbon and nuclear waste storage, resource extraction, and predicting earthquakes. Chemical processes controlling subcritical fracture in calcite are poorly understood. We demonstrate a novel approach to quantify the coupled chemical-mechanical effects on subcritical fracture. The calcite surface was indented using a Vickers-geometry indenter tip, which resulted in repeatable micron-scale fractures propagating from the indent. Individual indented samples were submerged in an array of aqueous fluids and an optical microscope was used to track the fracture growth in situ. The fracture propagation rate varied from 1.6 × 10 -8 m s -1 to 2.4 × 10 -10 m s -1 . The rate depended on the type of aqueous ligand present, and did not correlate with the measured dissolution rate of calcite or trends in zeta-potential. We postulate that chemical complexation at the fracture tip in calcite controls the growth of subcritical fracture. Previous studies indirectly pointed to the zeta-potential being the most critical factor, while our work indicates that variation in the zeta-potential has a secondary effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conflation of Geospatial POI Data and Ground-level Imagery via Link Prediction on Joint Semantic Graph

With the proliferation of smartphone cameras and social networks, we have rich, multi-modal data about points of interest (POIs) - like cultural landmarks, institutions, businesses, etc. - within a given areas of interest (AOI) (e.g., a county, city or a neighborhood) available to us. Data conflation across multiple modalities of data sources is one of the key challenges in maintaining a geographical information system (GIS) which accumulate data about POIs. Given POI data from nine different sources, and ground-level geo-tagged and scene-captioned images from two different image hosting platforms, in this work we explore the application of graph neural networks (GNNs) to perform data conflation, while leveraging a natural graph structure evident in geospatial data. The preliminary results demonstrate the capacity of a GNN operation to learn distributions of entity (POIs and images) features, coupled with topological structure of entity's local neighborhood in a semantic nearest neighbor graph, in order to predict links between a pair of entities.

Gurav, Rutuja↗

Quantifying Variability and Controls of Riverine Dissolved Organic Carbon Exported to Arctic Coastal Margins of North America (Final Report)

This project involved implementation and application of a coupled permafrost hydrology and dissolved organic carbon process models to investigate how spatial and seasonal variations in terrestrial hydrology and soil freeze/thaw dynamics influence the mobilization, loading, and export of organic carbon to the stream network for selected arctic basins across northern Alaska and northwest Canada. The model simulations were constrained by detailed observations of in- stream chemistry, soil active layer profile moisture and temperature dynamics, streamflow, soil carbon inventories and satellite microwave remote sensing based assessments of surface soil freeze-thaw dynamics. We developed and applied the numerical modeling and data analysis, incorporating observed data for calibration and validation, in order to investigate the terrestrial hydrology, permafrost dynamics, and associated DOC production and loading to rivers across a region encompassing watersheds draining to the coast. The project produced six publications and three datasets archived in public repositories.

54 ENVIRONMENTAL SCIENCES↗

A Graphics Processing Unit (GPU) Approach to Large Eddy Simulation (LES) for Transport and Contaminant Dispersion

Recent advances in the development of large eddy simulation (LES) atmospheric models with corresponding atmospheric transport and dispersion (AT&D) modeling capabilities have made it possible to simulate short, time-averaged, single realizations of pollutant dispersion at the spatial and temporal resolution necessary for common atmospheric dispersion needs, such as designing air sampling networks, assessing pollutant sensor system performance, and characterizing the impact of airborne materials on human health. The high computational burden required to form an ensemble of single-realization dispersion solutions using an LES and coupled AT&D model has, until recently, limited its use to a few proof-of-concept studies. An example of an LES model that can meet the temporal and spatial resolution and computational requirements of these applications is the joint outdoor-indoor urban large eddy simulation (JOULES). A key enabling element within JOULES is the computationally efficient graphics processing unit (GPU)-based LES, which is on the order of 150 times faster than if the LES contaminant dispersion simulations were executed on a central processing unit (CPU) computing platform. JOULES is capable of resolving the turbulence components at a suitable scale for both open terrain and urban landscapes, e.g., owing to varying environmental conditions and a diverse building topology. In this paper, we describe the JOULES modeling system, prior efforts to validate the accuracy of its meteorological simulations, and current results from an evaluation that uses ensembles of dispersion solutions for unstable, neutral, and stable static stability conditions in an open terrain environment.

54 ENVIRONMENTAL SCIENCES↗

NAS - Supercomputing master tool for aeronautics

Features, performance objectives and applications planned for the NASA National Aerodynamics Simulator (NAS) are outlined. NAS was conceived in the 1970s as a means to performing numerical aerodynamic simulations beyond the scope of wind tunnel testing for high-speed flight. Present supercomputers cannot deal with problems exhibiting strongly coupled viscous effects, which are being increasingly more accurately represented by approximations to the full Navier-Stokes equations. Located at the NASA-Ames Center, the NAS will by 1990 comprise a distributed computer network capable of a 4 Gflop computing rate and have a memory capacity of 1 billion 64-bit words. Remote access to the system through UNIX-based microcomputers will be available through land lines and satellite links. New supercomputers will be tested on the system without disturbing ongoing work. The core machine will be a Cray-2 with a 2 Gflop rate.

Bailey, F. R.↗

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems, task 1: Ducted propfan analysis

The time-dependent three-dimensional Euler equations of gas dynamics were solved numerically to study the steady compressible transonic flow about ducted propfan propulsion systems. Aerodynamic calculations were based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. An implicit residual smoothing operator was used to aid convergence. Two calculation grids were employed in this study. The first grid utilized an H-type mesh network with a branch cut opening to represent the axisymmetric cowl. The second grid utilized a multiple-block mesh system with a C-type grid about the cowl. The individual blocks were numerically coupled in the Euler solver. Grid systems were generated by a combined algebraic/elliptic algortihm developed specifically for ducted propfans. Numerical calculations were initially performed for unducted propfans to verify the accuracy of the three-dimensional Euler formulation. The Euler analyses were then applied for the calculation of ducted propfan flows, and predicted results were compared with experimental data for two cases. The three-dimensional Euler analyses displayed exceptional accuracy, although certain parameters were observed to be very sensitive to geometric deflections. Both solution schemes were found to be very robust and demonstrated nearly equal efficiency and accuracy, although it was observed that the multi-block C-grid formulation provided somewhat better resolution of the cowl leading edge region.

Hall, Edward J.↗