Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel coordinates”

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 217 records · Page 12

Space Station Human Factors Research Review. Volume 1: EVA Research and Development

An overview is presented of extravehicular activity (EVA) research and development activities at Ames. The majority of the program was devoted to presentations by the three contractors working in parallel on the EVA System Phase A Study, focusing on Implications for Man-Systems Design. Overhead visuals are included for a mission results summary, space station EVA requirements and interface accommodations summary, human productivity study cross-task coordination, and advanced EVAS Phase A study implications for man-systems design. Articles are also included on subsea approach to work systems development and advanced EVA system design requirements.

Cohen, Marc M.↗

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict 3-D ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗

Gigatraj: An Atmospheric Trajectory Model

Atmospheric trajectory models have a long history of success in tracking air motions in the lower stratosphere and upper troposphere over periods of up to a few days. Parcels have been traced backwards from observations to identify whatever phenomena (strong convection, volcanic eruptions, rocket launches, etc.) put their signature on them. Parcels have also been initialized at a known event and traced forward to examine their subsequent physical and chemical evolution. We describe a new trajectory model, "gigatraj," that aims to increase exibility by (a) making it straightforward to use new meteorological data sources, including those not based on regular latitude-longitude grids; (b) enabling a run-time choice of vertical coordinate system for kinematic and/or quasi-isentropic calculations; (c) allowing for the output of arbitrary meteorological products, selectable by the user and interpolated to the parcels' locations and times. The model can be run in a serial or parallel processing environment, so that large numbers of parcels can be traced in a reasonable time. Information is presented on model accuracy and performance. The former is demonstrated by runs using both test-pattern winds (comparing expected paths with actual output) and real-world winds (comparing forward and backward runs to characterize how well parcels retrace their paths). Sample cases are also shown, including a reverse domain lling (RDF) calculation illustrating a tropopause fold event. Model output can be displayed using the new Visualization And Lagrangian dynamics Immersive eXtended Reality (VALIXR) system, and an example will be shown. In addition, we describe work to incorporate a version of gigatraj into the Goddard Earth Observing System (GEOS) of NASA's Global Modeling and Assimilation O ce (GMAO) at Goddard Space Flight Center. This enables trajectory calculations to be performed within the running GEOS model at the latter's native time resolution, instead of using the winds from every few hours. It also provides access to all of GEOS's internal variables as they are calculated. This module may be useful, for example, for tracking rapid chemical changes in a Lagrangian framework.

dynamics↗

Uranium (VI) Adsorbate Structures on Portlandite [Ca(OH)2] Type Surfaces Determined by Computational Modelling and X-ray Absorption Spectroscopy

Portlandite [Ca(OH) 2 ] is a potentially dominant solid phase in the high pH fluids expected within the cementitious engineered barriers of Geological Disposal Facilities (GDF). This study combined X-ray Absorption Spectroscopy with computational modelling in order to provide atomic-scale data which improves our understanding of how a critically important radionuclide (U) will be adsorbed onto this phase under conditions relevant to a GDF environment. Such data are fundamental for predicting radionuclide mass transfer. Surface coordination chemistry and speciation of uranium with portlandite [Ca(OH) 2 ] under alkaline groundwater conditions (ca. pH 12) were determined by both in situ and ex situ grazing incidence extended X-ray absorption fine structure analysis (EXAFS) and by computational modelling at the atomic level. Free energies of sorption of aqueous uranyl hydroxides, [UO 2 (OH) n ] 2–n (n = 0–5) with the (001), (100) and (203) or (101) surfaces of portlandite are predicted from the potential of mean force using classical molecular umbrella sampling simulation methods and the structural interactions are further explored using fully periodic density functional theory computations. Although uranyl is predicted to only weakly adsorb to the (001) and (100) clean surfaces, there should be significantly stronger interactions with the (203/101) surface or at hydroxyl vacancies, both prevalent under groundwater conditions. The uranyl surface complex is typically found to include four equatorially coordinated hydroxyl ligands, forming an inner-sphere sorbate by direct interaction of a uranyl oxygen with surface calcium ions in both the (001) and (203/101) cases. In contrast, on the (100) surface, uranyl is sorbed with its axis more parallel to the surface plane. The EXAFS data are largely consistent with a surface structural layer or film similar to calcium uranate, but also show distinct uranyl characteristics, with the uranyl ion exhibiting the classic dioxygenyl oxygens at 1.8 Å and between four and five equatorial oxygen atoms at distances between 2.28 and 2.35 Å from the central U absorber. These experimental data are wholly consistent with the adsorbate configuration predicted by the computational models. These findings suggest that, under the strongly alkaline conditions of a cementitious backfill engineered barrier, there would be significant uptake of uranyl by portlandite to inhibit the mobility of U(VI) from the near field of a geological disposal facility.

36 MATERIALS SCIENCE↗

Noise generation by a finite span swept airfoil

An analysis is developed for the noise generated by a convected gust interacting with a single swept airfoil. The airfoil spans a channel with infinite parallel walls which contains a uniform subsonic mean flow. High frequency gusts, for which the noise generation is concentrated at the airfoil leading edge, are considered. The solution utilizes the Wiener-Hopf technique applied to the equations in nonorthogonal coordinates. Closed form expressions for the farfield acoustic pressure are obtained. The convected gust is assumed to be a plane wave, i.e., the lowest order spanwise mode. Because of the swept leading edge, this convected gust generates acoustic modes of all orders. At high frequencies and low Mach numbers, moderate sweep angles produce substantial noise reductions. The residual acoustic energy is carried mainly by the higher modes. The airfoil sweep also modulates the basic high frequency farfield directivity, producing lobed patterns. The angular modulation of the directivity patterns is most pronounced for the lower order modes.

Kerschen, E. J.↗

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict three-dimensional ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however, this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗

A Procedure for 3-D Contact Stress Analysis of Spiral Bevel Gears

Contact stress distribution of spiral bevel gears using nonlinear finite element static analysis is presented. Procedures have been developed to solve the nonlinear equations that identify the gear and pinion surface coordinates based on the kinematics of the cutting process and orientate the pinion and the gear in space to mesh with each other. Contact is simulated by connecting GAP elements along the intersection of a line from each pinion point (parallel to the normal at the contact point) with the gear surface. A three dimensional model with four gear teeth and three pinion teeth is used to determine the contact stresses at two different contact positions in a spiral bevel gearset. A summary of the elliptical contact stress distribution is given. This information will be helpful to helicopter and aircraft transmission designers who need to minimize weight of the transmission and maximize reliability.

Kumar, A.↗

Computational aeroacoustics on massively parallel computers

A computational fluid dynamics (CFD) technique is employed to solve aeroacoustics problems on massively parallel computers. The algorithm is based on a 4th order accurate central finite difference and 4th order accurate 4-stage Runge-Kutta time integration method that solves the 3-D full Navier-Stokes/Euler equations on a curvilinear coordinate system. The code has been developed to predict noise radiation from ducted fans. A specific application of this code is made here to the oscillating circular piston problem (Workshop Category 4, Problem 2) by solving the full Euler equations with nonreflecting boundary conditions in the axisymmetric mode. The features of the code that apply to this particular problem are described in this paper.

Oezyoeruek, Yusuf↗

LIS-Hydro: Authoritative Source for OCONUS Hydro-Intelligence

U.S. military forces are often tasked to participate in a variety of transboundary water-related decision-making activities, including humanitarian assistance operations through Department of State tasking, support of in-country infrastructure development activities that help develop or improve diplomatic relationships, and support of transboundary water treaty negotiations or disputes to reduce risk of conflict caused by water security issues. The U.S. intelligence communities have identified the coordination over shared water resources as an area of significant concern to U.S. national security (U.S. National Defense Strategy, 2018). Such transboundary water issues are projected to intensify in the future under increasingly complex population dynamics, political tensions due to parallel issues, and a changing climate. A 2017 joint NASA, USACE/RDC, and U.S. Air Force co-sponsored workshop revealed a lack of sufficient decision support tools and access to timely technical and contextual information needed to assess and respond to potential water-related threats around the world. The need for an integrated operational service, with the capacity to combine and synthesize hydrological modeling, assimilation, forecasting, and visualization capabilities across the U.S. Government, was highlighted as a key recommendation. In direct response, a subset of the U.S. Department of Defense, National Intelligence Community, and Oak Ridge National Laboratory are collaborating on the development of a fully integrated hydro-modeling and streamflow prediction system (i.e., LIS-Hydro). Completion of the project and sustainment of the operational capability by Air Force Weather will establish a national asset to assist federal agencies implement government-wide strategies around water resources (U.S. Global Water Strategy, 2017). The hydrological products and services will, for the very first time, establish a routinely available authoritative source of global water intelligence information supporting war-fighters, planners, and decision makers at all echelons and services of the U.S. military, Federal government, and intelligence community. A summary of the interagency scientific collaboration in addressing some of the key gaps and needs identified during the 2017 workshop will be presented.

Land Information System↗

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Scalable Comparative Visualization of Ensembles of Call Graphs

Optimizing the performance of large-scale parallel codes is critical for efficient utilization of computing resources. Code developers often explore various execution parameters, such as hardware configurations, system software choices, and application parameters, and are interested in detecting and understanding bottlenecks in different executions. They often collect hierarchical performance profiles represented as call graphs, which combine performance metrics with their execution contexts. The crucial task of exploring multiple call graphs together is tedious and challenging because of the many structural differences in the execution contexts and significant variability in the collected performance metrics (e.g., execution runtime). In this paper, we present Ensemble CallFlow to support the exploration of ensembles of call graphs using new types of visualizations, analysis, graph operations, and features. We introduce ensemble-Sankey , a new visual design that combines the strengths of resource-flow (Sankey) and box-plot visualization techniques. Whereas the resource-flow visualization can easily and intuitively describe the graphical nature of the call graph, the box plots overlaid on the nodes of Sankey convey the performance variability within the ensemble. Our interactive visual interface provides linked views to help explore ensembles of call graphs, e.g., by facilitating the analysis of structural differences, and identifying similar or distinct call graphs. Finally, we demonstrate the effectiveness and usefulness of our design through case studies on large-scale parallel codes.

97 MATHEMATICS AND COMPUTING↗

Manufactured Porous Ambient Surface Simulants

The planetary science decadal survey for 2013-2022 (Vision and Voyages, NRC 2011) has promoted mission concepts for sample acquisition from small solar system bodies. Numerous comet-sampling tools are in development to meet this standard. Manufactured Porous Ambient Surface Simulants (MPASS) materials provide an opportunity to simulate variable features at ambient temperatures and pressures to appropriately test potential sample acquisition systems for comets, asteroids, and planetary surfaces. The original "flavor" of MPASS materials is known as Manufactured Porous Ambient Comet Simulants (MPACS), which was developed in parallel with the development of the Biblade Comet Sampling System (Backes et al., in review). The current suite of MPACS materials was developed through research of the physical and mechanical properties of comets from past comet missions results and modeling efforts, coordination with the science community at the Jet Propulsion Laboratory and testing of a wide range of materials and formulations. These simulants were required to represent the physical and mechanical properties of cometary nuclei, based on the current understanding of the science community. Working with cryogenic simulants can be tedious and costly; thus MPACS is a suite of ambient simulants that yields a brittle failure mode similar to that of cryogenic icy materials. Here we describe our suite of comet simulants known as MPACS that will be used to test and validate the Biblade Comet Sampling System (Backes et al., in review).

Carey, Elizabeth M.↗

Overview of NASA ISRU Plans, Priorites, and Activities

Introduction:The National Aeronautics andSpace Administration (NASA) of the United States ofAmerica (US) has initiated the Artemis Moon to Marsprogram to send astronauts (the first woman andperson of color) back to the lunar surface, create asustainable human lunar exploration program, andlead the first human exploration mission to the Marssurface in the 2030’s [1]. A major objective of thisprogram is to characterize the resources that exist onthe Moon and Mars, and learn how to utilize them forsustained and affordable exploration. Commonlyknown as In Situ Resource Utilization (ISRU), thesearch for, acquisition, and processing of resources inspace has the potential to greatly reduce thedependency on transporting mission consumables andinfrastructure from Earth, thereby reducing missioncosts, risks, and dependency on Earth.ISRU is Enabling: Through the extraction andprocessing of resources into mission commoditiessuch as rocket propellants, life support consumables,and fuel cell reactants, ISRU enhances and evolvesthe cis-lunar, lander, and surface transportationsystems required for human exploration; expandingand enhancing HOW humans can explore and returnfrom the Moon. Through the extraction andprocessing of resources into metals, silicon, and othermanufacturing and construction feedstock, ISRUenhances and allows for the expansion of criticalinfrastructure using in situ manufacturing andconstruction capabilities that influence WHAT humanscan do on the Moon and in cis-lunar space. Becauseof this, ISRU supports and enables commercialinvolvement beyond NASA and governmentalagencies by both lowering the cost of sustainedtransportation to/from/on the Moon as well assupporting the market required for needing thesetransportation systems. Strategic Framework:To achieve this vision,NASA’s Space Technology Mission Directorate(STMD) ensures the coordinated development ofISRU and other critical space and surfaceinfrastructure elements such as propulsion, power,manufacturing, construction, and robotics through theStrategic Technology Architecture Roundtable(STAR) process. Through STAR, an integratedframework and process has been created allowing forcapabilities and technologies to be linked andassessed, gaps to be identified, specifications andmetrics to be established, and provide a means toprioritize and implement technology development andmissions. A critical part of the STAR effort has beenthe establishment of the Strategic Framework thatorganizes all work under four major Thrusts (Go,Land, Live, and Explore) and identifies the drivingOutcomes for each of these Thrusts. From the Thrustsand Outcomes, all work can be categorized and linkedbetween Capability Areas, and Technology Gaps canbe identified and addressed (Figure 1.)Figure 1. Strategic Framework and STAR FrameworkISRU Envisioned Future: To drive thedevelopment of technologies and capabilities, theSTAR process starts with establishing a ‘grand vision’of where each Outcome and Capability is aiming tobe considered complete. For ISRU, the EnvisionedFuture is “Scalable ISRU production/utilizationcapabilities including sustainable commodities on thelunar and Mars Surface”. This involves starting with10’s of metric tons of products, but evolves into 100’sto 1000’s of metric tons of water, oxygen, propellants,construction and manufacturing feedstock, andcommodities for habitat and food production andoperations. For ISRU, the ‘Prospect to Product’philosophy starts with Destination Reconnaissance &Resource Assessment, followed by ResourceAcquisition, Isolation, and Preparation, leading intoResource Processing (which is further subdivided intomission consumables and feedstocks for constructionand manufacturing). The ISRU Envisioned Futurealso considers what resources are available andattempts to address what and when these resourceswill be evaluated and harnessed, as well asconsidering which products/commodities can beobtained for early use and which ones require moretime and/or users of refined products.It Takes an Architecture: ISRU does not existon its own. By definition, it requires customers/users SHORT TITLE HERE: A. B. Author and C. D. Authorto use the products/commodities produced by ISRUsystems. Also, for an ISRU capability to exist, itmust obtain products and services from other systemsand infrastructure. An important aspect of the STARprocess and the ISRU Envisioned Futures Prioritiesstrategy is to identify and link all of these systems andcapabilities to achieve the desired end state (Figure2).Figure 2. ISRU as Part of a Larger ArchitectureISRU Capability Drivers: The guidingprinciples for NASA’s Space TechnologyDevelopment for Artemis are to develop criticaltechnologies and capabilities that enable (i) asustainable Lunar surface presence, (ii) the future goalof sending humans to Mars, and (iii) promotingcritical technologies to enable future science andcommercial missions. It is a major goal of theArtemis campaign to establish some sort of base campat the lunar South Pole by approximately the end ofthe decade. The ISRU Envisioned Futures Prioritiesstrategy is aligned with the Artemis campaign todevelop and demonstrate ISRU capabilities in thistimeframe that could lead to sustained surfaceoperations, infrastructure growth, and commercialoperations in the next decade (Figure 3).Figure 3. ISRU Dual Path to Full Implementation and CommercializationState of the Art and Gaps: To achieve theenvisioned future, an extensive effort was performedto understand the State of the Art (SOA) for ISRUgoing back decades, and to assess the SOA against thenear and long-term goals and objectives of the ISRUStrategic Outcome objectives. While the releasedISRU Envisioned Futures Priorities only includes atop-level definition of both the SOA and Gaps, furtherinformation on these for ISRU can be found in theISRU Gap Assessment Study performed for theInternational Space Exploration Coordination Group(ISECG) [2]. To provide further guidance to industryand academia, a top level assessment was performedand provide that divides critical areas of ISRUcapabilities and technologies into 3 categories:Significant Funding, Partially Covered/MoreRequired, and Limited/No Funded Activities.Envisioned Future Priorities- Next Steps forISRU: While a significant amount of work over abroad range of technology areas has been performedover the last several years for lunar ISRU, to reach theenvisioned future for ISRU, a lot more work isrequired at the technology level leading to bothsystems and technology demonstrations in the nearfuture. To guide investments within NASA, industry,and academia, 5 specific areas of high priority wereidentified. These are:1.Complete development of the Water and Oxygen Mining Paths and close technology gaps, with emphasis on oxygen extraction from Highland regolith and parallel paths for polar water mining.2.Expand development of metal extraction and feedstock for manufacturing and construction, with emphasis on aluminum and initial/easy to obtain/make construction feedstocks leading to more refined metals and other regolith resources. Also, evaluate biologically inspired/derived technologies in bio-mining, bio-plastic, and other feedstock commodities.3.Ensure the resource assessment needed for future ISRU commercial operations is coordinated with both near/long-term science objectives as well as Artemis mission locations of interest.4.Initiate NASA and industry-led system-level analyses, integration, and testing activities for ISRU capabilities. While significant work has been performed at the technology and subsystemlevel, it is now important to understand how these technology investments can be leveraged and utilized in actual systems and applications5.Initiate lunar ISRU technology flight demonstrations leading to initial ‘Pilot Plant’ end-to-end production capability demonstrations, led by industry

ISRU↗

Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation

The coordinate-space representation of the Hartree-Fock-Bogoliubov theory is the method of choice to study weakly bound nuclei whose properties are affected by the quasiparticle continuum space. To describe such systems, we developed a three-dimensional Skyrme-Hartree-Fock-Bogoliubov solver HFBFFT based on the existing, highly optimized and parallelized Skyrme-Hartree-Fock code Sky3D. The code does not impose any self-consistent spatial symmetries such as mirror inversions or parity. The underlying equations are solved in HFBFFT directly in the canonical basis using the fast Fourier transform. To remedy the problems with pairing collapse, we implemented the soft energy cutoff and pairing annealing. The convergence of HFB solutions was improved by a sub-iteration method. The Hermiticity violation of differential operators brought by Fourier-transform-based differentiation has also been solved. Furthermore, the accuracy and performance of HFBFFT were tested by benchmarking it against other HFB codes, both spherical and deformed, for a set of nuclei, both well-bound and weakly-bound.

3D coordinate-space representation↗

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗

Application of the multigrid method to grid generation

The multigrid method (MGM), used to numerically solve the pair of nonlinear elliptic equations commonly used to generate two dimensional boundary-fitted coordinate systems is discussed. Two different geometries are considered: one involving a coordinate system fitted about a circle and the other selected for an impinging jet flow problem. Two different relaxation schemes are tried: one is successive point overrelaxation and the other is a four-color scheme vectorizeable to take advantage of a parallel processor computer for greater computational speed. Results using MGM are compared with those using SOR (doing successive overrelaxations with the corresponding relaxation scheme on the fine grid only). It is found that MGM becomes significantly more effective than SOR as more accuracy is demanded and as more corrective grids, or more grid points, are used. For the accuracy required, it is found that MGM is two to three times faster than SOR in computing time. With the four-color relaxation scheme as applied to the impinging jet problem, the advantage of MGM over SOR is not as great. This may be due to the effect of a poor initial guess on MGM for this problem.

Ohring, S.↗

Distributed Stochastic Optimization of a Neural Representation Network for Time-Space Tomography Reconstruction

4D time-space reconstruction of dynamic events or deforming objects using X-ray computed tomography (CT) is an important inverse problem in non-destructive evaluation. Conventional back-projection based reconstruction methods assume that the object remains static for the duration of several tens or hundreds of X-ray projection measurement images (reconstruction of consecutive limited-angle CT scans). However, this is an unrealistic assumption for many in-situ experiments that causes spurious artifacts and inaccurate morphological reconstructions of the object. To solve this problem, we propose to perform a 4D time-space reconstruction using a distributed implicit neural representation (DINR) network that is trained using a novel distributed stochastic training algorithm. Our DINR network learns to reconstruct the object at its output by iterative optimization of its network parameters such that the measured projection images best match the output of the CT forward measurement model. Here, we use a forward measurement model that is a function of the DINR outputs at a sparsely sampled set of continuous valued 4D object coordinates. Unlike previous neural representation architectures that forward and back propagate through dense voxel grids that sample the object's entire time-space coordinates, we only propagate through the DINR at a small subset of object coordinates in each iteration resulting in an order-of-magnitude reduction in memory and compute for training. DINR leverages distributed computation across several compute nodes and GPUs to produce high-fidelity 4D time-space reconstructions. We use both simulated parallel-beam and experimental cone-beam X-ray CT datasets to demonstrate the superior performance of our approach.

36 MATERIALS SCIENCE↗

Bounce-averaged Hamiltonian for charged particles in an axisymmetric but nondipolar model magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M,J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar vector-B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. Our formulation provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗