Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data intensive computing”

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

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

At least 1,225 records · Page 68

Woven ceramic matrix composite surrogate model based on physics-informed recurrent neural network

A recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear constitutive behavior of woven ceramic matrix composites (CMCs) driven by matrix damage at multiple length scales. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the multiscale generalized method of cells (MSGMC) approach coupled with a matrix damage model. This coupling permits simulating the nonlinear behavior of woven CMCs based on constituent response at the micro-, meso-, and macroscales. The multiscale repeating unit cell is loaded under non-monotonic conditions including multiple load / unload cycles and tension / compression. The fiber volume fraction as well as the intra- and intertow void volume fractions are also varied in the generation of training data. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input strain sequence and fiber and void volume fractions, the resulting stress versus strain response while satisfying physical constraints such as positive semi-definiteness of the tangent stiffness matrix and linear elastic unloading. Further, the trained surrogate model effectively matches the stress versus strain response and successfully predicts the tangent modulus throughout the loading regime. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex architectures, nonlinear multiaxial material response, and under non-monotonic loading conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation of Proton Beams in the AGS Booster with Space Charge

As higher intensity accelerators are designed, built and operated, the resonances and instabilities caused by space charge can become dominant effects in limiting the intensity of the beam that can be delivered. Particle loss has many causes but space charge can be a major factor in reducing beam intensity and generating emittance growth, particularly at low energies and in the lower energy stages of accelerators such as the AGS Booster. Successful operation relies on a good understanding of basic particle dynamics and it becomes an important challenge to ascertain how the beam evolves when charged particles experience both external forces through the fields of the magnets combined with internal forces through Coulomb interactions. A possible approach to such a study is to include calculations of space charge in computer simulation codes that track the motions of charged particles in 6D phase space. Unfortunately, repeated space charge calculations in 6D phase space are computationally intensive even with the most modern computing facilities. Study with traditional models of space charge is then either limited to a short real-time scale, or sacrifices accuracy and often omits some of the physics. The study described in this report attempts to explore all aspects of the Booster from injection, through RF capture and acceleration to the extraction energy. The simulations are allied to real machine settings. They do not attempt to model effects that are well known and for which accurate experimental data are available, but they try in many cases to extend beyond the normal Booster parameter range so as to understand the limitations on operations and suggest areas for improvement.

43 PARTICLE ACCELERATORS↗

Objective method for analysis and tracking of convective cells as seen by radar

A special method has been developed for the study of cells that are embedded in convective rain systems. This method consists of a package of computer programs that use pattern recognition techniques on three-dimensional digital radar data to identify the rain cells, track them with time, and calculate their properties. The product of the computations is a comprehensive database of physically meaningful properties of rain cells, which can be used to infer the internal structure and the dynamics of convective rain systems. The cell-tracking method has been applied to the summer convective clouds of south Florida for the following purposes: (1) derivation of the relationship between the echo top height and the precipitation characteristics (e.g., area, water yield, rain intensity and duration of the rain cells); (2) study of the microphysical behaviorof cumulus clouds in relation to their cell properties; (3) evaluation of the effect of seeding on cumulus clouds on the cell scale; and (4) examination of cloud-to-ground lightning discharges in relation to convective cell intensity. The cell-tracking method is also currently being used in rain enhancement projects in Texas, Israel, and South Africa. The cell-tracking method, its products and their use in meteorological research are described in this paper.

Rosenfeld, Daniel↗

Simulation of spatially evolving turbulence and the applicability of Taylor's hypothesis in compressible flow

For the numerical simulation of inhomogeneous turbulent flows, a method is developed for generating stochastic inflow boundary conditions with a prescribed power spectrum. Turbulence statistics from spatial simulations using this method with a low fluctuation Mach number are in excellent agreement with the experimental data, which validates the procedure. Turbulence statistics from spatial simulations are also compared to those from temporal simulations using Taylor's hypothesis. Statistics such as turbulence intensity, vorticity, and velocity derivative skewness compare favorably with the temporal simulation. However, the statistics of dilatation show a significant departure from those obtained in the temporal simulation. To directly check the applicability of Taylor's hypothesis, space-time correlations of fluctuations in velocity, vorticity, and dilatation are investigated. Convection velocities based on vorticity and velocity fluctuations are computed as functions of the spatial and temporal separations. The profile of the space-time correlation of dilatation fluctuations is explained via a wave propagation model.

Lee, Sangsan↗

Dynamic Inlet Distortion Prediction with a Combined Computational Fluid Dynamics and Distortion Synthesis Approach

A procedure has been developed for predicting peak dynamic inlet distortion. This procedure combines Computational Fluid Dynamics (CFD) and distortion synthesis analysis to obtain a prediction of peak dynamic distortion intensity and the associated instantaneous total pressure pattern. A prediction of the steady state total pressure pattern at the Aerodynamic Interface Plane is first obtained using an appropriate CFD flow solver. A corresponding inlet turbulence pattern is obtained from the CFD solution via a correlation linking root mean square (RMS) inlet turbulence to a formulation of several CFD parameters representative of flow turbulence intensity. This correlation was derived using flight data obtained from the NASA High Alpha Research Vehicle flight test program and several CFD solutions at conditions matching the flight test data. A distortion synthesis analysis is then performed on the predicted steady state total pressure and RMS turbulence patterns to yield a predicted value of dynamic distortion intensity and the associated instantaneous total pressure pattern.

Norby, W. P.↗

A Composite Diagnosis of Synoptic-Scale Extratropical Cyclone Development over the United States

This paper presents a composite diagnosis of synoptic-scale forcing mechanisms associated with extratropical cyclone evolution. Drawn from 12 cyclone cases that occurred over the continental United States during the cool season months, the diagnosis provides a 'climatology' of development mechanisms for difference categories of cyclone evolution ranging from cyclone weakening through three stages of cyclone intensification. Computational results were obtained using an 'extended' form of the Zwack-Okossi equation applied to routine upper-air and surface data analyzed on a 230 km x 230 km grid. Results show that cyclonic vorticity advection, which maximizes in the upper troposphere, was the primary contributor to cyclone development regardless of the stage of development. A second consistent contributor to development was latent heat release. Horizontal temperature advection, often acknowledged as a development mechanism, was found to contribute to development only during more intense stages. During weakening and weaker development stages, temperature advection opposed development, as the warm-air advection invariably found at upper levels was dominated by cold air advection in the lower half of the troposphere. In the more intense stages, development was moderated by dry-adiabatic cooling associated with the ascending vertical motions.

Rolfson, Donald M.↗

Validation of numerical models for seismic fluid-structure-interaction analysis of nuclear, safety-related equipment

Seismic design, qualification, and risk assessment of nuclear safety-related vessels filled with liquid will have to consider the interaction between the vessel (tank), the contained liquid, and submerged components, if any. Seismic fluid-structure-interaction (FSI) analysis of nuclear vessels will rely on numerical models, which are required to be verified and validated. This study validates previously-verified numerical models using test data generated from earthquake-simulator experiments involving a cylindrical tank. Two solvers in LS-DYNA for FSI simulations are used for the numerical analysis: Arbitrary-Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD). Numerical and test results are compared for responses critical to seismic design of advanced reactor vessels: hydrodynamic pressures on the tank wall, reactions at the support, and wave heights of the contained liquid. Analysis is performed for one-, two-, and three-directional seismic inputs with a range of intensity, and rocking motions. The accuracy of the numerical results and different methods for outputting wave heights are discussed. Recommendations for validation of seismic FSI numerical models of advanced reactors are provided. As a result, the validation exercise presented in this study is broadly applicable to cylindrical tanks, regardless of industry sector.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

GenAI4UQ: A software for forward and inverse uncertainty quantification using conditional generative AI

We introduce GenAI4UQ, a software package for forward and inverse uncertainty quantification in model calibration, parameter estimation, and ensemble forecasting. GenAI4UQ leverages a generative AI-based conditional modeling framework to address limitations of traditional inverse modeling techniques, such as Markov Chain Monte Carlo (MCMC) methods. By replacing computationally intensive iterative processes with a direct, learned mapping, GenAI4UQ enables efficient calibration of input parameters and generation of predictions directly from observations. The software supports rapid ensemble forecasting with robust uncertainty quantification while maintaining computational and storage efficiency. Built-in auto-tuning of hyperparameters simplifies model training, ensuring accessibility for users with varying expertise. Its versatile conditional generative framework is applicable across diverse scientific domains. While GenAI4UQ offers significant advantages in flexibility and efficiency, users should interpret its uncertainty estimates with caution in data-sparse scenarios, as the model may overestimate uncertainty—an effect common to all surrogate-based approaches including MCMC with surrogate models. Despite this, GenAI4UQ transforms inverse modeling by providing a fast, reliable, and user-friendly solution. It empowers researchers and practitioners to quickly estimate parameter distributions and generate model predictions for new observations, facilitating efficient decision-making and advancing the state of uncertainty quantification in computational modeling.

97 MATHEMATICS AND COMPUTING↗

Laboratory infrared spectroscopy of gaseous negatively charged polyaromatic hydrocarbons

Based largely on infrared spectroscopic evidence, polycyclic aromatic hydrocarbon (PAH) molecules are now widely accepted to occur abundantly in the interstellar medium. Laboratory infrared spectra have been obtained for a large variety of neutral and cationic PAHs, but data for anionic PAHs are scarce. Nonetheless, in regions with relatively high electron densities and low UV photon fluxes, PAHs have been suggested to occur predominantly as negatively charged ions (anions), having substantial influence on cloud chemistry. While some matrix spectra have been reported for radical anion PAHs, no data is available for even-electron anions, which are more stable against electron detachment. Here we present the first laboratory infrared spectra of deprotonated PAHs ([PAH-H]{sup –}) in the wavelength ranges between 6 and 16 μm and around 3 μm. Wavelength-dependent infrared multiple-photon electron detachment is employed to obtain spectra for deprotonated naphthalene, anthracene, and pyrene in the gas phase. Spectra are compared with theoretical spectra computed at the density functional theory level. We show that the relative band intensities in different ranges of the IR spectrum deviate significantly from those of neutral and positively charged PAHs, and moreover from those of radical anion PAHs. These relative band intensities are, however, well reproduced by theory. An analysis of the frontier molecular orbitals of the even- and odd-electron anions reveals a high degree of charge localization in the deprotonated systems, qualitatively explaining the observed differences and suggesting unusually high electric dipole moments for this class of PAH molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗

Cloud parameters from GOES visible and infrared radiances during the FIRE Cirrus IFO, October 1986

Visible (VIS, 0.65 micron) and infrared (IR, 10.5 microns) channels on geostationary satellites are the key elements of the International Satellite Cloud Climatology Project (ISCCP). All daytime ISCCP cloud parameters are derived from a combination of VIS and IR data. Validation and improvement of the ISCCP and other cloud retrieval algorithms are important components of the First ISCCP Regional Experiment (FIRE) Intensive Field Observations (IFO). Data from the Cirrus IFO (October 19 to November 2, 1986) over Wisconsin are available for validating cirrus cloud retrievals from satellites. The Geostationary Operational Environmental Satellite (GOES) located over the Equator at approximately 100 deg W provided nearly continuous measurements of VIS and IR radiances over the IFO areas. The preliminary results of cloud parameters derived from the IFO GOES data are presented. Cloud attitudes are first derived using an algorithms without corrections for cloud emissivity. These same parameters will then be computed from the same data relying on an emissivity correction algorithm based on correlative data taken during the Cirrus IFO.

Heck, Patrick W.↗

Automated inspection of turbine blades: Challenges and opportunities

Current inspection methods for complex shapes and contours exemplified by aircraft engine turbine blades are expensive, time-consuming and labor intensive. The logistics support of new manufacturing paradigms such as integrated product-process development (IPPD) for current and future engine technology development necessitates high speed, automated inspection of forged and cast jet engine blades, combined with a capability of retaining and retrieving metrology data for process improvements upstream (designer-level) and downstream (end-user facilities) at commercial and military installations. The paper presents the opportunities emerging from a feasibility study conducted using 3-D holographic laser radar in blade inspection. Requisite developments in computing technologies for systems integration of blade inspection in production are also discussed.

Mehta, Manish↗

Memory-Intensive Benchmarks: IRAM vs. Cache-Based Machines

The increasing gap between processor and memory performance has lead to new architectural models for memory-intensive applications. In this paper, we explore the performance of a set of memory-intensive benchmarks and use them to compare the performance of conventional cache-based microprocessors to a mixed logic and DRAM processor called VIRAM. The benchmarks are based on problem statements, rather than specific implementations, and in each case we explore the fundamental hardware requirements of the problem, as well as alternative algorithms and data structures that can help expose fine-grained parallelism or simplify memory access patterns. The benchmarks are characterized by their memory access patterns, their basic control structures, and the ratio of computation to memory operation.

Biswas, Rupak↗

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗

Colour thresholding and objective quantification in bioimaging

Computer imaging is rapidly becoming an indispensable tool for the quantification of variables in research and medicine. Whilst its use in medicine has largely been limited to qualitative observations, imaging in applied basic sciences, medical research and biotechnology demands objective quantification of the variables in question. In black and white densitometry (0-256 levels of intensity) the separation of subtle differences between closely related hues from stains is sometimes very difficult. True-colour and real-time video microscopy analysis offer choices not previously available with monochrome systems. In this paper we demonstrate the usefulness of colour thresholding, which has so far proven indispensable for proper objective quantification of the products of histochemical reactions and/or subtle differences in tissue and cells. In addition, we provide interested, but untrained readers with basic information that may assist decisions regarding the most suitable set-up for a project under consideration. Data from projects in progress at Tulane are shown to illustrate the advantage of colour thresholding over monochrome densitometry and for objective quantification of subtle colour differences between experimental and control samples.

NASA Discipline Neuroscience↗

Effect of different breathing patterns in the same patient on stereotactic ablative body radiotherapy dosimetry for primary renal cell carcinoma: A case study

Stereotactic ablative body radiotherapy (SABR) for primary renal cell carcinoma (RCC) targets requires motion management strategies to verify dose delivery. This case study highlights the effect of a change in patient breathing amplitude on the dosimetry to organs at risk and target structures. A 73-year-old male patient was planned for receiving 26 Gy of radiation in 1 fraction of SABR for a left primary RCC. The patient was simulated with four-dimensional computed tomography (4DCT) and the tumor internal target volume (ITV) was delineated using the 4DCT maximum intensity projection. However, the initially planned treatment was abandoned at the radiation oncologist's discretion after pretreatment cone-beam CT (CBCT) motion verification identified a greater than 50% reduction in superior to inferior diaphragm motion as compared with the planning 4DCT. This patient was resimulated with respiratory coaching instructions. To assess the effect of the change in breathing on the dosimetry to the target, each plan was recalculated on the data set representing the change in breathing condition. A change from smaller to larger breathing showed a 46% loss in planning target volume (PTV) coverage, whereas a change from larger breathing to smaller breathing resulted in an 8% decrease in PTV coverage. ITV coverage was similarly reduced by 8% in both scenarios. This case study highlights the importance of tools to verify breathing motion prior to treatment delivery. 4D image guided radiation therapy verification strategies should focus on not only verifying ITV margin coverage but also the effect on the surrounding organs at risk.

61 RADIATION PROTECTION AND DOSIMETRY↗

Observation of laser ablation of silicon as a function of pulse length at constant fluence via time-resolved x-ray spectroscopy

We investigate the ablation of silicon as a function of laser pulse length at a constant fluence using time-resolved x-ray spectroscopy data obtained from OMEGA EP experiments at the University of Rochester's Laboratory for Laser Energetics. Our targets consisted of three-layer planar structures composed of Si (50 μm), Cu (25 μm), and SiO 2 (500 μm) layers. The Si layer was irradiated by a 351-nm laser with varying pulse widths of 250 ps, 500 ps, 1 ns, and 10 ns while maintaining a constant fluence of ~27.9 kJ/cm 2 . Electron temperatures and densities of the ablated plasma were determined by analyzing the time-resolved x-ray spectroscopy data through a comparison of experimental measurements with synthetic results obtained from Si atomic calculations in a steady state and non-local thermodynamic equilibrium. These calculations were computed using PrismSPECT. Additionally, radiation-hydrodynamics simulations with FLASH are used to generate simulated plasma-density and plasma-temperature profiles, which are then compared with the experimental measurements. Our analyses reveal that increasing the laser pulse length at a constant fluence results in a decrease in electron temperatures and densities. Furthermore, the longer pulses with lower intensities lead to deeper ablation regions before reaching the peak ablation but lower ionization balances in the silicon layer. Furthermore, these findings emphasize the critical role of laser pulse length in plasma ablation and shock generation for laser-impulse studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation on experimental techniques to detect, locate and quantify gear noise in helicopter transmissions

A robotic system to automate the detection, location, and quantification of gear noise using acoustic intensity measurement techniques has been successfully developed. Major system components fabricated under this grant include an instrumentation robot arm, a robot digital control unit and system software. A commercial, desktop computer, spectrum analyzer and two microphone probe complete the equipment required for the Robotic Acoustic Intensity Measurement System (RAIMS). Large-scale acoustic studies of gear noise in helicopter transmissions cannot be performed accurately and reliably using presently available instrumentation and techniques. Operator safety is a major concern in certain gear noise studies due to the operating environment. The man-hours needed to document a noise field in situ is another shortcoming of present techniques. RAIMS was designed to reduce the labor and hazard in collecting data and to improve the accuracy and repeatability of characterizing the acoustic field by automating the measurement process. Using RAIMS a system operator can remotely control the instrumentation robot to scan surface areas and volumes generating acoustic intensity information using the two microphone technique. Acoustic intensity studies requiring hours of scan time can be performed automatically without operator assistance. During a scan sequence, the acoustic intensity probe is positioned by the robot and acoustic intensity data is collected, processed, and stored.

Flanagan, P. M.↗