Results for “computation time” — Page 94

NASA NTRSDate not supplied

Simulation of Dynamic Crushing Behavior of Simulated Carbon/Epoxy Stanchions Using MAT58 and MAT213 in LS-DYNA

Carbon fiber-reinforced polymer (CFRP) composites have shown promise as a material for structures designed to absorb energy in crush-style impact loading. In complementing the lightweight and tailorable characteristics of CFRP structures, the benefits of impact energy absorption are particularly interesting in aerospace vehicles that require excellent crashworthiness performance. However, simulating the behavior of CFRP structures in dynamic, crushing-style impact loading is challenging because of the many damage and failure modes that are essential to capture in the model. For the present investigation, CFRP stanchions were tested using a crash sled experimental test rig. The stanchions were designed by the Composite Materials Handbook-17 (CMH-17) Crashworthiness Working Group for the purpose of comparing experimental crash sled tests to corresponding simulations of the tests. The stanchions are C-channel shaped and represent a geometry common in the interior of aerospace vehicle structures. Explicit simulations in LSDYNA were performed using a well-established composite material model (MAT58) and a next-generation material model (MAT213). Fully integrated shell elements were utilized instead of 3D solid elements to limit computation time. In all models, each of the sixteen plies were represented with individual layers of shell elements with tiebreak contact between each adjacent ply pair to simulate interlaminar fracture. Simulating the crushing of the stanchions occurred in two phases. First, the material models were calibrated using flat specimens that were manufactured with the same layup as the stanchions. While the original goal was calibration of material-related properties, meshdependent behavior was observed in simulations with either material model, and an unstructured mesh was selected to remediate undesirable mesh-dependent failure modes. Additionally, for both the MAT58 and MAT213 models, it was found that either the crush force or the failure mode could be modeled accurately, but no set of parameters could be identified to attain both results in the same model. Once satisfactory calibration was achieved, the same material parameters were applied to the stanchion crush simulations. The stanchion simulations showed that MAT213 more accurately predicted the experimentally determined crush force, and both material models predicted key aspects of the experimentally observed failure modes.

Composite Materials↗

NASA NTRSDate not supplied

Energy of Fluid (EOF) Simulation of Large-Scale Densification and Solidification of Hydrogen

The current research develops a finite volume-based Computation Fluid Dynamic (CFD) model utilizing an Energy of Fluid (EOF) approach to generate a simulation of the densification of liquid hydrogen (LH2) in the Integrated Refrigeration and Storage (IRAS) experimental tank for the Ground Operations Demonstration Unit for Liquid Hydrogen Project (GODU-LH2) at the Kennedy Space Center (KSC). The computational code will integrate a commercial software pressure-based model with User Defined Functions (UDF) generated and added to the model. The modified model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. This method will decrease computational time when compared to the built-in commercially available Enthalpy solvers. Currently the research focused on utilizing User Defined Scalars (UDS) compared to the enthalpy model and a lumped node analysis to show the written UDS is as accurate as the built-in solvers, and the models will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

C. Mahony↗

NASA NTRS2023

Real-Time Hardware-in-the-Loop Evaluation of A Partially Turboelectric Propulsion Control Design

In support of aviation fuel burn and emission reduction goals, NASA is pursuing high-payoff research investments that promise to transform aviation. This includes investments in Electrified Aircraft Propulsion (EAP). Multiple technology challenges must be addressed to unlock the full potential of EAP. This includes addressing challenges related to propulsion controls, which will be vital for ensuring efficient coordinated operation of EAP subsystems. This paper presents results from real-time hardware-in-theloop (HIL) testing of a control design for a single aisle partially-turboelectric aircraft propulsion concept conducted at the NASA Electric Aircraft Testbed (NEAT) facility. The control system under test is designed for a propulsion concept consisting of two wing-mounted turbofan engines that produce thrust and generate electrical power to drive a boundary layer ingesting tailfan propulsor via an electrical motor. An integrated control strategy is applied to ensure coordinated operation of the turbofan and tailfan subsystems during steady-state and transient operation throughout the flight envelope. The NEAT test of this integrated control design consists of a partially HIL, partially simulated configuration. A subscale representation of the electrical system design is implemented in hardware and mechanically coupled to electric machines that emulate turbomachinery and propulsor shaft dynamics. The hardware configuration is then operated under the control of a real-time computer application that runs a simulation of the propulsion system and the developed control logic. The NEAT facility test campaign includes a series of experiments that subject the control design to throttle transients conducted throughout the flight envelope and full-flight mission profiles. Testing under simulated performance degradation is also conducted to evaluate control design robustness. This includes constant and abrupt changes in degradation levels. Results from the HIL test are presented and shown to be in good agreement with pretest simulation predictions demonstrating the efficacy of the integrated control design approach.

Electrified Aircraft Propulsion↗

NASA NTRSDate not supplied

Sequential Filtering in the Presence of Uniform Measurement Errors

This paper presents a sequential filtering strategy using observations corrupted with uniform measurement noise. While the Kalman filter remains the best linear estimator of the state, other filtering techniques provide minimum variance optimal estimates, a trait only enjoyed by the Kalman filter when the underlying noises are, in fact, Gaussian. This work develops a new approximate optimal estimator for uniform measurement noises. The resulting recursion requires just slightly more computational time to complete a measurement update than the Kalman filter, which generally cannot be claimed by other optimal strategies such as the particle or Gaussian mixture filters.

James S McCabe↗

NASA NTRSDate not supplied

Forward and Inverse Models for Satellite Remote Sensors using Principal Component Analysis

Satellite remote sensors such as AIRS on Aqua, CrIS on S-NPP, NOAA20 and JPSS-2, IASI on Metop A, B, and C make millions of observations each day with thousands of spectral channels for each observation; this poses challenges for efficiently inversion of the inherently large dataset as needed to retrieve atmospheric and surface properties. This presentation will illustrate the use of Principal Component Analysis (PCA) to speed up radiative transfer forward model calculations and to stabilize the inversion algorithms. A Principal Component-based radiative transfer model (PCRTM) developed at NASA Langley Research Center can simulate top of atmosphere (TOA) radiance or reflectance spectra from 50 cm-1 to 50000 cm-1 (200 m to 0.20 m quickly and accurately. PCRTM demonstrated very high accuracy relative to reference line-by-line radiative transfer models and it saves orders of magnitude computational time. Examples of the PCRTM model developed for hyperspectral sensors such as AIRS, CrIS, IASI, NAST-I, SHIS, CPF, TEMPO, SBG, OMI, and SCIAMACHY will be presented. In addition to using the PCRTM as forward model, the NASA Langley developed inversion algorithm also uses PCA to compress the state vector into a compressed dimension to speed up and stabilize the inversion process. Examples of retrieved atmospheric temperature, water vapor, CO2, CO, CH4, N2O, and O3 profiles, cloud properties (optical depth, size, phase, and height), and surface properties (surface emissivity spectra and skin temperatures) will be presented. This algorithm is being transitioned to the NASA Sounder SIPS and NASA's Goddard Earth Sciences Data and Information Services Center (GES DISC).

forward model↗

NASA NTRS2023

Space-Shielding Radiation Dosage Code Evaluation - Phase 1: SHIELDOSE-2 Radiation-Assessment Code

Radiation-transport codes or radiation-analysis tools are used in the spacecraft-design community to define and assess the local radiation levels. Among various radiation-transport/analysis tools available, SHIELDOSE-2 is commonly used in the early spacecraft-design phases because it has been established that it provides reasonably reliable results with relatively rapid computation time. The NASA Engineering and Safety Center initiated a two-phased study: (1) to better understand and document the SHIELDOSE-2 code capabilities and limitations, and (2) to recommend alternative tools to use for cases where SHIELDOSE-2 is not applicable. This report provides a summary of the Phase 1 study.

SHIELDOSE-2↗

NASA NTRSDate not supplied

Initial Development of A Digital Twin Model for an Electrified Aircraft Propulsion Emulation Rig

In support of aviation fuel burn and emission reduction goals, NASA is pursing high-payoff research investments that promise to transform aviation. This includes investments in Electrified Aircraft Propulsion (EAP), which relies on the generation, storage, transmission, and use of electrical power for producing thrust and optimizing propulsion system efficiency. Multiple technology challenges must be addressed to unlock the full potential of EAP. This includes advances in propulsion controls, which will be vital for ensuring coordinated efficient operation of the complex integrated subsystems that comprise EAP architectures. To support EAP controls research, the NASA Glenn Research Center has developed the Hybrid Propulsion Emulation Rig (HyPER). The HyPER laboratory hardware includes shaft-mounted electric machines, power converters, power supplies, power distribution cables, and an energy storage device that can be reconfigured to represent a variety of EAP architectures. It also includes an integrated real-time computer system that hosts developed EAP control software and turbomachinery simulations. This enables the electrical system and rotating shafts of EAP designs to be implemented in actual hardware and integrated with turbomachinery simulations and system-level EAP control logic implemented in software. In this form, the HyPER laboratory provides a partially simulated, partially hardware-in-the-loop test environment enabling the initial development and evaluation of EAP control technology. A prerequisite for the development of EAP control designs is the availability of a system model that accurately reflects the operation of the electrical system hardware. To support this need, a digital twin model of the HyPER electrical system hardware is under development. This model is being coded in the MATLAB Simulink environment and uses the NASA-developed Electrical Modeling and Thermal Analysis Toolbox (EMTAT) to construct a digital twin framework. EMTAT contains generic electrical component building blocks that are simulated at turbomachinery timescales. Associated inputs and outputs allow the blocks to be combined to model complete electrical systems. The EMTAT blocks also contain adjustable internal maps and parameters that can be set to reflect the operation of a specific electrical component. For the HyPER digital twin, the settings of these EMTAT block internal maps and parameters is determined through machine learning approaches applied to characterization run data collected from the laboratory. During characterization runs the laboratory electrical system hardware is subjected to a full range of torque, speed, and power settings. Acquired data is then used to estimate EMTAT block parameters using a variety of machine learning techniques. The resulting digital twin model is found to match the operation of actual HyPER hardware with an accuracy suitable for control development purposes. It also holds promise for other applications including modeling the performance of HyPER laboratory reconfigurations and model-based anomaly detection. Planned follow-on work to automate post-processing of acquired laboratory data to update the HyPER digital twin model will also be presented and discussed.

Electrified Aircraft Propulsion↗

NASA NTRS2023

Simulations of Zero Boiloff, Densification, and Solidification in a Large-Scale Liquid Hydrogen Tank

The current research utilizes an Energy of Fluid (EOF) approach to develop a finite volume-based Computation Fluid Dynamic (CFD) model to create a benchmark simulation of the densification of liquid hydrogen (LH2) in an experimental Integrated Refrigeration and Storage (IRAS) tank at the Kennedy Space Center (KSC). The computational code will incorporate the FLUENT pressure-based model with User Defined Functions (UDF) to implement the EOF model. The enhanced model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. Specifically, a convection term will be added to an already existing code for a more accurate simulation of LH2 densification. However, currently the research focuses on utilizing User Defined Scalars (UDS) with FLUENT to show the built-in solvers are as accurate as the written UDS. This will reduce the computational time down significantly by utilizing an already built in solver for the momentum equation. A custom UDF has been written to solve transient and steady state conduction, and results obtained using the UDF agree with results obtained using the built-in FLUENT solver. Currently, a convection UDF is being added and will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

Colin Philip Mahony↗

NASA NTRS2023

Numerical Investigation of Occupant Injury Risks During A Realistic Transport Aircraft Crash Conditions

Researchers at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) have conducted a full-scale crash test of a Fokker F28 MK1000 aircraft to investigate the performance of transport aircraft under realistic crash conditions. This crash test was computationally recreated using finite element (FE) human body models (HBMs) to further explore potential injury risks to occupants and analyze the utilization of HBMs in the aerospace crash environment. The Global Human Body Model Consortium (GHBMC) male 50th percentile occupant detailed model (v6.0) and the Toyota Human Model for Safety (THUMS) male 50th percentile occupant model (v6.1) were selected to be used in the crash simulations. The HBMs were simulated in conditions matching those of anthropomorphic test device (ATD) experiments included within the aircraft cabin during the crash test. Seven occupant locations within the cabin were simulated utilizing each of the models. The models were positioned in a neutral upright posture with hands resting on the legs and the feet contacting the floor. Head, brain, neck, and lumbar vertebra injury metrics were calculated for all trials. Both HBM models required minor modifications to stabilize these simulations. The GHBMC model required added erosion for six parts while the THUMS model only required one in order to complete the full simulation. The THUMS model, however, required a much smaller timestep for stability and therefore took significantly more computational time. In addition the GHBMC model includes integrated instrumentation while the THUMS model requires development and implementation of instrumentation. Both models predicted 100% injury risk for lumbar vertebra fracture in all test conditions. This prediction was in family with high lumbar load values measured by the ATDs during the crash test. The THUMS model consistently predicted lower injury risks than the GHBMC model in all three other metrics varying depending on the crash pulse. Overall, the THUMS model required less modifications to allow for this study. However, the GHBMC models significantly faster run time and integrated instrumentation make it a more intuitive model for this research.

Crashworthiness↗

NASA NTRSDate not supplied

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work uses an automated grid adaptation tool recently implemented in VULCAN-CFD, and explores this tool’s ability to solve highspeed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user-generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flowfield parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Furthermore, motivated by the grid convergence index (GCI) used for structured grid simulations, a grid convergence estimate (GCE) was developed and demonstrated for the grid adaptation. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex high-speed mixing problems with only a fraction of the grid cells required to obtain comparable results on a user-generated structured grid.

hypersonics↗

NASA NTRSDate not supplied

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work uses an automated grid adaptation tool recently implemented in VULCAN-CFD, and explores this tool’s ability to solve highspeed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user-generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flowfield parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Furthermore, motivated by the grid convergence index (GCI) used for structured grid simulations, a grid convergence estimate (GCE) was developed and demonstrated for the grid adaptation. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex high-speed mixing problems with only a fraction of the grid cells required to obtain comparable results on a user-generated structured grid.

hypersonics↗

NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗

NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗

NASA NTRS2024

Space-Shielding Radiation Dosage Code Evaluation Phase 2: SHIELDOSE-2 Radiation-Assessment Code

Radiation-transport codes or radiation-analysis tools are used in the spacecraft-design community to define and assess the local radiation levels. Among various radiation-transport/analysis tools available, SHIELDOSE-2 is commonly used in the early spacecraft-design phases because it has been established that it provides reasonably reliable results with relatively rapid computation time. The NASA Engineering and Safety Center initiated a two-phased study: (1) to better understand and document the SHIELDOSE-2 code capabilities and limitations, and (2) to identify alternative tools to use for cases where SHIELDOSE-2 is not applicable. This report provides a summary of the Phase 2 study.

SHIELDOSE-2↗

NASA NTRSDate not supplied

Developing Fast and Accurate Radiative Transfer Models to Meet the Needs of Modern Satellite Remote Sensing Applications

Modern hyperspectral satellite remote sensors provide highly accurate measurements the Earth’s Top-of-Atmosphere (TOA) radiance, reflectance, or polarized spectra with hundreds to thousands of spectral channels and with millions of observations per day. The large data volume and high spectral dimensionality of the data pose challenges for retrieval algorithms. To process the satellite Level-1 data (e.g. calibrated TOA spectra) into Level-2 products (e.g. atmospheric and surface properties) using physical-based retrieval algorithms, accurate and fast Radiative Transfer Models (RTMs) are needed. RTMs are usually the limiting factor in determining the speed of a level-2 algorithm. For example, more than one million Line-by-Line (LBL) radiative transfer (RT) calculations are needed in order to properly capture the spectral contributions of important atmospheric molecules for an IR hyperspectral sensor with a spectral coverage from 3.5 m to 15 m or a solar hyperspectral sensor with spectral coverage from 0.25 m to 2.5 m. In this presentation, we will discuss advantages and disadvantages of different ways (e.g. correlated k and effective transmittance) to accelerate the speed of a fast RTM. We finally describe a Principal Component-based Radiative Transfer Model (PCRTM), which can calculate TOA radiance or reflectance spectra from 50 cm-1 to 40,000 cm-1 (200 m to 0.25 m). It has demonstrated very good accuracy relative to reference LBL RTMs and saves orders of magnitude in computational time. The PCRTM has been used in many satellite remote sensing applications. Examples include forward modeling in Level-2 and Level-3 retrieval algorithms, high fidelity satellite instrument simulators and instrument performance trade studies, spectral and radiometric accuracy characterizations of satellite Level-1 data, tools for inter-satellite calibrations, tools for satellite RTM lookup table generations, and tools for generating physically based training datasets for Artificial Intelligence (AI) algorithms.

climate data record↗

DOE OSTI2020

Dynamic Equivalence of Large-Scale Power Systems Based on Boundary Measurements

Parallel computing helped speed up many tasks that can be done independently. The advance in gaming industry did not produce a technology that helps the power industry to reduce the computation time for dynamic simulation in individual cases. To perform faster than real-time simulation for the purpose of predicting power system dynamic trajectory, power industry continues to struggle to reduce the system size and simulation time of large-scale power systems while keeping its dynamic behavior under various disturbances. There is also an acute need for real-time dynamic model reduction as more renewables enter the generation mix with dramatic changes in the generation outputs. All existing model reduction solutions are based on having access to a detailed dynamic model of the system. The system changes by the minutes, dynamic models are only updated annually. To solve this problem, wide-area measurements obtained by the phasor measurement unit (PMU) at the boundaries between the reduced system and the study system is used to represent the external area. An artificial neuro-fuzzy inference system (ANFIS) is established to perform the mapping of the measurement to the external equivalent model. Here the external area is regarded as a black-box. Model reduction studies are conducted on the Northeast Power Coordinating Council (NPCC) under various types of contingencies, by using a co-simulation approach between PSS/E and MATLAB. The results look very promising and will be discussed in this paper.

Tong, Ning↗

DOE OSTI2020

Theoretical Development of Controller Transfer applied to Dynamical Systems

Data-driven control strategies suffer from the necessity of training and its associated computational time and cost. This motivates the continuous increase of popularity of transfer learning based methods. In this paper, we have investigated transfer learning in the context of controller for dynamical systems. First, we have derived conditions under which controller can be transferred between two linear dynamical systems. Furthermore, we have introduced a novel algorithm to design the transferable and non-transferable control components for linear systems. Secondly, we have derived conditions under which a stabilizing controller can be transferred between two different nonlinear dynamical systems. Furthermore, we have used the notion of zero dynamics to transform the nonlinear dynamics into global normal form, to derive the transferability conditions on the controller. Finally, we have numerically evaluated the performance of our transfer learning based methods for two sets of example linear and nonlinear dynamical systems.

Chakraborty, Indrasis↗

DOE OSTI2021

Improving scalability of electronic structure code for molecular simulations in the presence of environment

A scalable density functional electronic code with Gaussian basis set, called UTEP-NRLMOL, is developed to perform simulations of molecular systems in the presence of the environment with particular attention to the memory requirements. In the electronic structure calculations, the memory and computation time are proportional to the number of atoms. Memory requirements for density functional calculations scale as N*N, where N is the number of atoms. While the recent advances in HPC offer platforms with large numbers of cores, the limited amount of memory available on a given node and poor scalability of the electronic structure codes hinder their efficient usage of these platforms. We have introduced new scaling and parallelization paradigms using MPI-3 shared-memory functionality combined with usage of sparse algebra and storage of matrices in sparse format. This extends the range of applicability of the UTEP-NRLMOL code to large systems over 10,000 atoms, or using up to 67,000 basis functions, and making use of HPC architectures using over 6,000 processors utilizing all available cores. We have also interfaced code with effective fragment potential and polarizable continuum model libraries. The code was used in simulations of several applications which are published in reputed scientific journals.

74 ATOMIC AND MOLECULAR PHYSICS↗