SEARCH · Engineering Papers

Results for “Moving mesh”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS2005

Computing Gravitational Fields of Finite-Sized Bodies

A computer program utilizes the classical theory of gravitation, implemented by means of the finite-element method, to calculate the near gravitational fields of bodies of arbitrary size, shape, and mass distribution. The program was developed for application to a spacecraft and to floating proof masses and associated equipment carried by the spacecraft for detecting gravitational waves. The program can calculate steady or time-dependent gravitational forces, moments, and gradients thereof. Bodies external to a proof mass can be moving around the proof mass and/or deformed under thermoelastic loads. An arbitrarily shaped proof mass is represented by a collection of parallelepiped elements. The gravitational force and moment acting on each parallelepiped element of a proof mass, including those attributable to the self-gravitational field of the proof mass, are computed exactly from the closed-form equation for the gravitational potential of a parallelepiped. The gravitational field of an arbitrary distribution of mass external to a proof mass can be calculated either by summing the fields of suitably many point masses or by higher-order Gauss-Legendre integration over all elements surrounding the proof mass that are part of a finite-element mesh. This computer program is compatible with more general finite-element codes, such as NASTRAN, because it is configured to read a generic input data file, containing the detailed description of the finiteelement mesh.

Quadrelli, Marco↗

NASA NTRS2002

CAPRI: Using a Geometric Foundation for Computational Analysis and Design

CAPRI (Computational Analysis Programming Interface) is a software development tool intended to make computerized design, simulation and analysis faster and more efficient. The computational steps traditionally taken for most engineering analysis (Computational Fluid Dynamics (CFD), structural analysis, etc.) are: Surface Generation, usually by employing a Computer Aided Design (CAD) system; Grid Generation, preparing the volume for the simulation; Flow Solver, producing the results at the specified operational point; Post-processing Visualization, interactively attempting to understand the results. It should be noted that the structures problem is more tractable than CFD; there are fewer mesh topologies used and the grids are not as fine (this problem space does not have the length scaling issues of fluids). For CFD, these steps have worked well in the past for simple steady-state simulations at the expense of much user interaction. The data was transmitted between phases via files. In most cases, the output from a CAD system could go IGES files. The output from Grid Generators and Solvers do not really have standards though there are a couple of file formats that can be used for a subset of the gridding (i.e. PLOT3D) data formats and the upcoming CGNS). The user would have to patch up the data or translate from one format to another to move to the next step. Sometimes this could take days. Instead of the serial approach to analysis, CAPRI takes a geometry centric approach. CAPRI is a software building tool-kit that refers to two ideas: (1) A simplified, object-oriented, hierarchical view of a solid part integrating both geometry and topology definitions, and (2) programming access to this part or assembly and any attached data. The connection to the geometry is made through an Application Programming Interface (API) and not a file system.

Haimes, Robert↗

NASA NTRS2002

Aeroelastic Calculations of Quiet High- Speed Fan Performed

An advanced high-speed fan was recently designed under a cooperative effort between the NASA Glenn Research Center and Honeywell Engines & Systems. The principal design goals were to improve performance and to reduce fan noise at takeoff. Scale models of the Quiet High-Speed Fan were tested for operability, performance, and acoustics. During testing, the fan showed significantly improved noise characteristics, but a self-excited aeroelastic vibration known as flutter was encountered in the operating range. Flutter calculations were carried out for the Quiet High-Speed Fan using a three-dimensional, unsteady aerodynamic, Reynolds-averaged Navier-Stokes turbomachinery code named "TURBO." The TURBO code can accurately model the viscous flow effects that can play an important role in various aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), and flutter in the presence of shock and boundary-layer interaction. Initially, calculations were performed with no blade vibrations. These calculations were at a constant rotational speed and a varying mass flow rate. The mass flow rate was varied by changing the backpressure at the exit boundary of the computational domain. These initial steady calculations were followed by aeroelastic calculations in which the blades were prescribed to vibrate harmonically in a natural mode, at a natural frequency, and with a fixed interblade phase angle between adjacent blades. The AE-prep preprocessor was used to interpolate the in-vacuum mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh and to smoothly propagate the grid deformations from the blade surface to the interior points of the grid. The aeroelastic calculations provided the unsteady aerodynamic forces on the blade surface due to blade vibrations. These forces were vector multiplied with the structural dynamic mode shape to calculate the work done on the blade during one vibration period, then this result was converted to an aerodynamic damping. Flutter occurs when the aerodynamic damping becomes negative, if structural damping is ignored. The results of these aeroelastic calculations are summarized in this plot of aerodynamic damping versus mass flow rate at a constant rotational speed. As the backpressure is increased, the mass flow rate through the fan decreases and the fan operating point moves towards the stall line. The aeroelastic calculations showed that the aerodynamic damping decreases as the stall line is approached, as observed during testing. In addition, the aeroelastic calculations with the TURBO code correctly predicted the aeroelastic parameters: the most unstable vibration mode and interblade phase angle, as observed during testing. The Quiet High-Speed Fan demonstrated significant noise reductions during testing, but flutter imposed limits on its operating range. The accurate calculation of the aeroelastic characteristics using the TURBO code is a significant step toward eliminating flutter from the operating range and toward realizing the benefits of reduced fan noise. The aeroelastic calculations described here were performed under a grant by University of Toledo researchers in collaboration with Glenn's researchers.

Bakhle, Milind A.↗

NASA NTRS1996

Computational modeling of magnetically actuated propellant orientation

Unlike terrestrial applications where gravity positions liquid at the 'bottom' of the tank, the location of liquid propellant in spacecraft tanks is uncertain unless specific actions are taken or special features are built into the tank. Some mission events require knowledge of liquid position prior to a particular action: liquid must be positioned over the tank outlet prior to starting the main engines and must be moved away from the tank vent before vapor can be released overboard to reduce pressure. It may also be desirable to positively position liquid to improve propulsion system performance: moving liquid away from the tank walls will dramatically decrease the rate of heat transfer to the propellant, suppressing the boil-off rate, thereby reducing overall mission propellant requirements. The process of moving propellant to a desired position is referred to as propellant orientation or reorientation. Several techniques have been developed to positively position propellant in spacecraft tanks and each technique imposes additional requirements on vehicle design. Propulsive reorientation relies on small auxiliary thrusters to accelerate the tank. The inertia of the liquid causes it to collect in the aft-end of the tank if the acceleration is forward. This technique requires that additional thrusters be added to the vehicle, that additional propellant be carried in the vehicle, and that an additional operational maneuver be executed. Another technique uses Liquid Acquisition Devices (LAD's) to positively position propellants. These devices rely on surface tension to hold the liquid within special geometries (i.e. vanes, wire-mesh channels, start-baskets). While avoiding some of the penalties of propulsive orientation, this technique requires the addition of complicated hardware inside the propellant tank and performance for long duration missions is uncertain. The subject of the present research is an alternate technique for positively positioning liquid within spacecraft propellant tanks: magnetic fields.

Hochstein, John I.↗

NASA NTRS1998

A Geometry Based Infra-Structure for Computational Analysis and Design

The computational steps traditionally taken for most engineering analysis suites (computational fluid dynamics (CFD), structural analysis, heat transfer and etc.) are: (1) Surface Generation -- usually by employing a Computer Assisted Design (CAD) system; (2) Grid Generation -- preparing the volume for the simulation; (3) Flow Solver -- producing the results at the specified operational point; (4) Post-processing Visualization -- interactively attempting to understand the results. For structural analysis, integrated systems can be obtained from a number of commercial vendors. These vendors couple directly to a number of CAD systems and are executed from within the CAD Graphical User Interface (GUI). It should be noted that the structural analysis problem is more tractable than CFD; there are fewer mesh topologies used and the grids are not as fine (this problem space does not have the length scaling issues of fluids). For CFD, these steps have worked well in the past for simple steady-state simulations at the expense of much user interaction. The data was transmitted between phases via files. In most cases, the output from a CAD system could go to Initial Graphics Exchange Specification (IGES) or Standard Exchange Program (STEP) files. The output from Grid Generators and Solvers do not really have standards though there are a couple of file formats that can be used for a subset of the gridding (i.e. PLOT3D data formats). The user would have to patch up the data or translate from one format to another to move to the next step. Sometimes this could take days. Specifically the problems with this procedure are:(1) File based -- Information flows from one step to the next via data files with formats specified for that procedure. File standards, when they exist, are wholly inadequate. For example, geometry from CAD systems (transmitted via IGES files) is defined as disjoint surfaces and curves (as well as masses of other information of no interest for the Grid Generator). This is particularly onerous for modern CAD systems based on solid modeling. The part was a proper solid and in the translation to IGES has lost this important characteristic. STEP is another standard for CAD data that exists and supports the concept of a solid. The problem with STEP is that a solid modeling geometry kernel is required to query and manipulate the data within this type of file. (2) 'Good' Geometry. A bottleneck in getting results from a solver is the construction of proper geometry to be fed to the grid generator. With 'good' geometry a grid can be constructed in tens of minutes (even with a complex configuration) using unstructured techniques. Adroit multi-block methods are not far behind. This means that a million node steady-state solution can be computed on the order of hours (using current high performance computers) starting from this 'good' geometry. Unfortunately, the geometry usually transmitted from the CAD system is not 'good' in the grid generator sense. The grid generator needs smooth closed solid geometry. It can take a week (or more) of interaction with the CAD output (sometimes by hand) before the process can begin. One way Communication. (3) One-way Communication -- All information travels on from one phase to the next. This makes procedures like node adaptation difficult when attempting to add or move nodes that sit on bounding surfaces (when the actual surface data has been lost after the grid generation phase). Until this process can be automated, more complex problems such as multi-disciplinary analysis or using the above procedure for design becomes prohibitive. There is also no way to easily deal with this system in a modular manner. One can only replace the grid generator, for example, if the software reads and writes the same files. Instead of the serial approach to analysis as described above, CAPRI takes a geometry centric approach. This makes the actual geometry (not a discretized version) accessible to all phases of the analysis. The connection to the geometry is made through an Application Programming Interface (API) and NOT a file system. This API isolates the top-level applications (grid generators, solvers and visualization components) from the geometry engine. Also this allows the replacement of one geometry kernel with another, without effecting these top-level applications. For example, if UniGraphics is used as the CAD package then Parasolid (UG's own geometry engine) can be used for all geometric queries so that no solid geometry information is lost in a translation. This is much better than STEP because when the data is queried, the same software is executed as used in the CAD system. Therefore, one analyzes the exact part that is in the CAD system. CAPRI uses the same idea as the commercial structural analysis codes but does not specify control. Software components of the CAD system are used, but the analysis suite, not the CAD operator, specifies the control of the software session. This also means that the license issues (may be) minimized and individuals need not have to know how to operate a CAD system in order to run the suite.

Haimes, Robert↗

DOE OSTI2020

Experimental Testing of Particle Erosion and Attrition in CSP Systems

In CSP systems using solid particles for heat transfer, erosion from the falling particles can damage the system components upon impact, and from sliding against the component walls. Furthermore, the resulting attrition of particles can change the thermal performance due to alterations in particle morphology. The research on these issues has been limited to industrial applications, and don’t necessarily coincide with the operating conditions in CSP systems. In this work, we investigated the durability of the particles selected as the heat transfer media, as well as special substrate materials that are of interest as lining materials to the particle storage tanks and hoppers. Three different erosion phenomena were studied: impact, abrasion, and attrition erosion. The low-temperature experiments revealed negligible losses from attrition and abrasion erosion. Although the substrates indicated visible surface damage from impact erosion, the actual material mass loss was small even after a throughput of 500kg of particles 63 days of 1MW CSP plant operation assuming 14 hours per day of operation. Erosion is typically classified as impact, abrasion, or attrition erosion. For example, the CSP systems utilizing solid particles as heat transfer media will likely experience impact erosion on the receiver hopper and heat exchanger walls as the solid particles fall from the reservoir. Additionally, there will be sliding erosion along the receiver walls as particles move through the hopper, and lastly, there will also be attrition erosion as the particles rub against each other. The extent of erosion depends on many factors including but not limited to particle shape and size, material hardness/brittleness and ductility, and particle impact speed and angle. Mechanistic erosion equations that are available in the literature are mostly derived from data pertaining to specific experiments and therefore can’t be applied to all applications without suitable adjustments. Survey of the literature relating to particle-based receivers revealed that the most prevalent approach has been to choose low velocities to limit erosion, but to our knowledge, little effort has been put in to fully characterize erosion at the expected operating conditions. A relatively recent work investigated the erosion of wire mesh proposed for use in a particle receiver, where the results actually showed an increase in mass as oxidation occurred on the mesh surface but the work did not investigate the particle attrition rate, and the erosion rate of a solid substrate. The end goal of this work is to develop a comprehensive particle and substrate durability model that will allow a comprehensive understanding of particle-based CSP operating conditions on component durability.

14 SOLAR ENERGY↗

DOE OSTI2024

h5bench: A unified benchmark suite for evaluating HDF5 I/O performance on pre‐exascale platforms

Summary Parallel I/O is a critical technique for moving data between compute and storage subsystems of supercomputers. With massive amounts of data produced or consumed by compute nodes, high‐performant parallel I/O is essential. I/O benchmarks play an important role in this process; however, there is a scarcity of I/O benchmarks representative of current workloads on HPC systems. Toward creating representative I/O kernels from real‐world applications, we have created h5bench , a set of I/O kernels that exercise hierarchical data format version 5 (HDF5) I/O on parallel file systems in numerous dimensions. Our focus on HDF5 is due to the parallel I/O library's heavy usage in various scientific applications running on supercomputing systems. The various tests benchmarked in the h5bench suite include I/O operations (read and write), data locality (arrays of basic data types and arrays of structures), array dimensionality (one‐dimensional arrays, two‐dimensional meshes, three‐dimensional cubes), I/O modes (synchronous and asynchronous). In this paper, we present the observed performance of h5bench executed along several of these dimensions on existing supercomputers (Cori and Summit) and pre‐exascale platforms (Perlmutter, Theta, and Polaris). h5bench measurements can be used to identify performance bottlenecks and their root causes and evaluate I/O optimizations. As the I/O patterns of h5bench are diverse and capture the I/O behaviors of various HPC applications, this study will be helpful to the broader supercomputing and I/O community.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Ion Exchange Processing of AP-105 Hanford Tank Waste through Crystalline Silicotitanate in a Staged 2- then 3-Column System (Rev.1)

The Tank Side Cesium Removal (TSCR) system, under development by Washington River Protection Solutions LLC (WRPS), will send initial low-activity Hanford waste tank supernate feeds to the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility. In addition to entrained solids removal from the supernate, the primary goal of TSCR is to remove cesium-137 ( 137 Cs) by ion exchange, allowing contact handling of the liquid effluent product at the WTP. Crystalline silicotitanate (CST) ion exchange media, manufactured by Honeywell UOP, LLC (product IONSIVTM R9140-B), was selected as the ion exchange media at TSCR. CST is a non-elutable inorganic material that has demonstrated robust chemical, physical, and radiation tolerance while maintaining functionality. However, exchange kinetics of Cs onto CST is slow, resulting in low utilization of the CST Cs load capacity before unacceptable Cs breakthrough. Two process flow designs have been tested, as follows. 1. Lead-lag column processing: The lead column was removed after the lag column effluent reached the waste acceptance criteria (WAC) limit, the lag column was moved into the lead position, and a new lag column was installed. This format used ~52% Cs load capacity on the lead column. 2. Lead-lag-polish column processing: The processing was stopped when the polish column effluent reached the WAC limit. This format resulted in 81% Cs load capacity on the lead column. Testing with diluted feed from Hanford tank AP-105 (AP-105DF) incorporated a nuanced change to the lead-lag-polish column system where the polish column was inserted when the lag column effluent reached WAC limit. A 10.9-L volume of AP-105DF (diluted to 5.6 M Na) was processed through the Direct Feed Test Platform system, established at Pacific Northwest National Laboratory to support small-scale waste qualification efforts. The columns consisted of 10-mL CST beds (CST Lot 2002009604, sieved to screen out >30-mesh particles) placed in 1.5-cm-inner-diameter columns. Feed was processed at 1.83 bed volumes (BV) per hour; the flowrate, in terms of contact time with the CST bed, matched the expected flowrate at TSCR. The <30-mesh CST sieve cut was expected to provide appropriate performance scaling to a full-height column. The installation of the polish column later in processing (after processing 523 BVs) did not appear to fundamentally change the utilization of the lead column for Cs exchange nor did it extend the total feed processing volume when compared to the previous test with AP-107 feed. Table ES.1 and Figure ES.1 summarize the measured AP-105DF Cs load performance.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2015

Los Alamos RAGE Simulations of the HAIV Mission Concept

The mitigation of potentially hazardous objects (PHOs) can be accomplished by a variety of methods including kinetic impactors, gravity tractors and several nuclear explosion options. Depending on the available lead time prior to Earth impact, non- nuclear options can be very effective at altering a PHOs orbit. However if the warning time is short nuclear options are generally deemed most effective at mitigating the hazard. The NIAC mission concept for a nuclear mission has been presented at several meetings, including the last PDC (2013).We use the adaptive mesh hydrocode RAGE to perform detailed simulations of this Hypervelocity Asteroid Intercept Vehicle (HAIV) mission concept. We use the RAGE code to simulate the crater formation by the kinetic impactor as well as the explosion and energy coupling from the follower nuclear explosive device (NED) timed to detonate below the original surface to enhance the energy coupling. The RAGE code has been well validated for a wide variety of applications. A parametric study will be shown of the energy and momentum transfer to the target 100 m diameter object: 1) the HAIV mission as planned; 2) a surface explosion and 3) a subsurface (contained) explosion; both 2) and 3) use the same source energy as 1).Preliminary RAGE simulations show that the kinetic impactor will carve out a surface crater on the object and the subsequent NED explosion at the bottom of the crater transfers energy and momentum to the target effectively moving it off its Earth crossing orbit. Figure 1 shows the initial (simplified) RAGE 2D setup geometry for this study. Figure 2 shows the crater created by the kinetic impactor and Figure 3 shows the time sequence of the energy transfer to the target by the NED.

Asteroid Mitigation↗

DOE OSTI2023

Niowave Neutron Source Converter: Lead-Bismuth-Eutectic (LBE) Windowless Target Design and Evaluation

Los Alamos National Laboratory (LANL) is working with Niowave on the design and evaluation of their lead-bismuth-eutectic (LBE) windowless target (i.e., neutron source converter). Niowave plans to use 200 kW electron beam at 40 MeV beam energy to produce neutrons by photonuclear reaction with LBE. Then, the neutrons undergo fission at the surrounding uranium target assembly (UTA) to produce Molybdenum 99 (Mo-99) as a fission product, which eventually decays to Technetium-99 (Tc-99m). Tc-99m is one of important radioisotopes that is used for medical diagnostics. LANL conducted 3D multiphysics analysis for the Niowave neutron converter design and provided design assessment in thermohydraulic aspects. LANL conducted radiation transport calculations using Monte-Carlo N-Particle (MCNP) code with unstructured meshing scheme. The 3D volumetric heating profiles in the LBE and Stainless-Steel (SS) housing were imported into multiphase computational fluid dynamics (CFD) to obtain 3D temperature profiles of LBE and SS through conjugate heat transfer (CHT) analysis. The key findings are: LBE film thickness at the center of the beam is approximately 1.6 cm with a maximum LBE velocity of approximately 1.8 m/s, which is below a 2 m/s limit to avoid erosion issues on supporting structures; Heat deposition in the LBE peaks at ~1 cm depth from the LBE free-surface because of the forward interactions of electron, photon, and neutron with LBE; LBE maximum temperature is ~360 °C, which is below LBE evaporation initiative temperature, ~450 °C; LBE-SS interface temperature is ~350 °C, which is below the safety thermal limit to prevent severe corrosion on SS. The results indicate that Niowave’s neutron converter design satisfies both hydraulic and thermal criteria for safe operation. By virtue of such computational analysis, Niowave can move toward establishing an experimental setup to experimentally test their LBE neutron converter. The following sections describe the detailed work done by LANL.

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗

NASA NTRS2019

Analytical Techniques for Assessing Gateway and Other Spacecraft Antenna Line-of-Sight

NASA’s Artemis program is committed to landing the next man and the first woman on the Moon by 2024. A critical piece of infrastructure for the mission, the Gateway will inaugurate a sustainable crewed presence beyond low-Earth orbit in cis-lunar space and serve as a staging point for a lunar landing system. The Near Rectilinear Halo Orbit of the Gateway offers numerous operational benefits in terms of its Earth access via the Orion Multipurpose Crew Vehicle, the degree of lunar visibility and surface access it provides, and favorable conditions for station-keeping. To minimize burns during station-keeping and conserve propellant, the Gateway will continually align itself to a solar pressure equilibrium attitude with the Sun; this is convenient for solar array orientation and power generation, but causes difficulties in terms of communication line-of-sight with the Earth, Moon, and visiting vehicles, which appear to be continually rotating about the Gateway’s local frame. A fixed antenna on the Gateway is continually subject to a moving target and can find itself obstructed by the Gateway’s own geometry frequently if not placed carefully. To understand the performance of either Earth-pointing or Moon-pointing antenna placements for a given Gateway geometry, a transient solar system simulation of the Gateway is used to step through orbital ephemeris data and determine the periods in which line-of-sight to a chosen target is achieved and lost. Antenna placements can be tested in the simulation and optimized to increase the average line-of-sight exposure to the target, while also minimizing the maximum duration cut-out in line-of-sight communication experienced; long durations without a successful communication link could pose threats to operations, safety, and mission success. It may be necessary to deviate from the Gateway’s solar pressure equilibrium attitude for short periods of time in order to continue a successful line-of-sight link to a chosen target. Communication coverage spheres can offer insight into what these necessary attitude adjustments may be and are generated by using a 3D mesh of the Gateway and an implementation of the Möller-Trumbore intersection algorithm. This paper outlines the analytical techniques used to perform such antenna placements on the Gateway, offering examples of how line-of-sight strength can be enhanced with an understanding of the Gateway geometry and attitude constraints. These techniques have had direct impact on the writing of communication requirements for NASA and international partner-owned Gateway elements.

Kennedy, William J.↗

NASA NTRS2019

Analytical Techniques for Assessing Gateway and Other Spacecraft Antenna Line-of-Sight

Kennedy, William J.↗

NASA NTRSDate not supplied

2020 International Astronautical Congress Conference Presentation - ANALYTICAL TECHNIQUES FOR ASSESSING GATEWAY AND OTHER SPACECRAFT ANTENNA LINE-OF-SIGHT FOR THE ARTEMIS PROGRAM

NASA’s Artemis program is committed to landing the next man and the first woman on the Moon by 2024. The Gateway, a critical piece of infrastructure for the long-term Artemis mission profiles, will inaugurate a sustainable crewed presence beyond low-Earth orbit in cis-lunar space and serve as a staging point for a lunar landing system. The Near Rectilinear Halo Orbit of the Gateway offers numerous operational benefits in terms of its Earth access via the Orion Multipurpose Crew Vehicle, its degree of Earth and lunar visibility, the versatile surface access it provides, and favorable conditions for station-keeping. To minimize burns during station-keeping and conserve propellant, the Gateway will continually align itself to a solar pressure equilibrium attitude with the Sun; this is convenient for solar array orientation and power generation, but causes difficulties in terms of communication line-of-sight with the Earth, Moon, and visiting vehicles, which appear to be continually rotating about the Gateway’s local frame. A fixed antenna on the Gateway is continually subject to a moving target and can find itself obstructed by the Gateway’s own geometry frequently if not placed carefully. To understand the performance of either Earth-pointing or Moon-pointing antenna placements for a given Gateway geometry, a transient solar system simulation of the Gateway is used to step through orbital ephemeris data and determine the periods in which line-of-sight to a chosen target is achieved and lost. Antenna placements can be tested in the simulation and optimized to increase the average line-of-sight exposure to the target, while also minimizing the maximum duration drop-out in line-of-sight communication experienced; long durations without a successful communication link could pose threats to operations, safety, and mission success. It may be necessary to deviate from the Gateway’s solar pressure equilibrium attitude for short periods of time in order to continue a successful line-of-sight link to a chosen target. Communication coverage spheres can offer insight into what these necessary attitude adjustments may be and are generated by using a 3D mesh of the Gateway and an implementation of the Möller-Trumbore intersection algorithm. This paper outlines the analytical techniques used to perform such antenna placements on the Gateway, offering examples of how line-of-sight strength can be enhanced with an understanding of the Gateway geometry and attitude constraints. These techniques have had direct impact on the writing of communication requirements for NASA and international partner-owned Gateway elements.

Gateway↗

NASA NTRS2002

Molten-Metal Droplet Deposition on a Moving Substrate in Microgravity: Aiding the Development of Novel Technologies for Microelectronic Assembly

Driven by advancements in microelectronics manufacturing, this research investigates the oblique (non-axisymmetric) impact of liquid-metal droplets on flat substrates. The problem of interest is relevant to the development of the novel technology of on-demand dispension (printing) of microscopic solder deposits for the surface mounting of microelectronic devices. The technology, known as solder jetting, features on-demand deposition of miniature solder droplets (30 to 120 microns in diameter) in very fine, very accurate patterns using techniques analogous to those developed for the ink-jet printing industry. Despite its promise, severe limitations exist currently with regards to the throughput rates of the technology; some of these limitations are largely due to the lack of the capability for reliable prediction of solder bump positioning and shapes, especially under ballistic deposition conditions where the droplet impact phenomena are inherently three-dimensional. The study consists of a theoretical and an experimental component. The theoretical work uses a finite element formulation to simulate numerically the non-axisymmetric (3-D) fluid mechanics and heat transfer phenomena of a liquid solder droplet impacting at an angle alpha on a flat substrate. The work focuses on the pre-solidification regime. The modeling of the most challenging fluid mechanics part of the process has been completed successfully. It is based upon the full laminar Navier-Stokes equations employing a Lagrangian frame of reference. Due to the large droplet deformation, the surface (skin) as well as the volumetric mesh have to be regenerated during the calculations in order to maintain the high accuracy of the numerical scheme. The pressure and velocity fields are then interpolated on the newly created mesh. The numerical predictions are being tested against experiments, for cases where wetting phenomena are not important. For the impact parameters used in the example shown (We = 2.38, Fr = 16300, Re = 157), the droplet rolls along the substrate, but its shape remains practically axisymmetric for all impact angles within the range from 0 to 60 deg. Interestingly, the substrate/droplet contact area during the recoiling phase of the impact is not a monotonically decreasing function of time. The experimental component of the research tests the numerical predictions and provides necessary input data (contact angles) for the theoretical model. The experiments are performed in microgravity (2.2s drop tower of the NASA GRC) in order to allow for the use of mm-size solder droplets, which make feasible the performance of accurate measurements, while maintaining similitude of the relevant fluid dynamic groups (Re, Fr, We, Ste). Preliminary oblique impact experiments have been performed using water droplets in normal gravity.

Megaridis, C. M.↗

NASA NTRS2017

Progress of High Efficiency Centrifugal Compressor Simulations Using TURBO

Three-dimensional, time-accurate, and phase-lagged computational fluid dynamics (CFD) simulations of the High Efficiency Centrifugal Compressor (HECC) stage were generated using the TURBO solver. Changes to the TURBO Parallel Version 4 source code were made in order to properly model the no-slip boundary condition along the spinning hub region for centrifugal impellers. A startup procedure was developed to generate a converged flow field in TURBO. This procedure initialized computations on a coarsened mesh generated by the Turbomachinery Gridding System (TGS) and relied on a method of systematically increasing wheel speed and backpressure. Baseline design-speed TURBO results generally overpredicted total pressure ratio, adiabatic efficiency, and the choking flow rate of the HECC stage as compared with the design-intent CFD results of Code Leo. Including diffuser fillet geometry in the TURBO computation resulted in a 0.6 percent reduction in the choking flow rate and led to a better match with design-intent CFD. Diffuser fillets reduced annulus cross-sectional area but also reduced corner separation, and thus blockage, in the diffuser passage. It was found that the TURBO computations are somewhat insensitive to inlet total pressure changing from the TURBO default inlet pressure of 14.7 pounds per square inch (101.35 kilopascals) down to 11.0 pounds per square inch (75.83 kilopascals), the inlet pressure of the component test. Off-design tip clearance was modeled in TURBO in two computations: one in which the blade tip geometry was trimmed by 12 mils (0.3048 millimeters), and another in which the hub flow path was moved to reflect a 12-mil axial shift in the impeller hub, creating a step at the hub. The one-dimensional results of these two computations indicate non-negligible differences between the two modeling approaches.

turbomachinery↗

NASA NTRS2000

Molecular Simulations in Astrobiology

One of the main goals of astrobiology is to understand the origin of cellular life. The most direct approach to this problem is to construct laboratory models of protocells. Such efforts, currently underway in the NASA Astrobiology Program, are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing energy and responding to changes in the environment, are carried out by proteins bound to membranes. We use computer simulations to address the following questions about these proteins: (1) How do small proteins self-organize into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides form membrane-spanning structures (e.g. channels)? (3) By what mechanisms do such structures perform their functions? The simulations are performed using the molecular dynamics method. In this method, Newton's equations of motion for each atom in the system are solved iteratively. At each time step, the forces exerted on each atom by the remaining atoms are evaluated by dividing them into two parts. Short-range forces are calculated in real space while long-range forces are evaluated in reciprocal space, using a particle-mesh algorithm which is of order O(NInN). With a time step of 2 femtoseconds, problems occurring on multi-nanosecond time scales (10(exp 6)-10(exp 8) time steps) are accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude, which requires algorithmic improvements and codes scalable to a large number of processors. Work in this direction is in progress. Two series of simulations are discussed. In one series, it is shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and fold into helical structures (see Figure). Once in the membrane, the peptides exhibit orientational flexibility with changing conditions, which may have provided a mechanism of transmitting signals between the protocell and its environment. In another series of simulations, the mechanism by which a simple protein channel efficiently mediates proton transport across membranes was investigated. This process is a key step in cellular bioenergetics. In the channel under study, proton transport is gated by four histidines that occlude the channel pore. The simulations identify the mechanisms by which protons move through the gate.

Pohorille, Andrew↗