Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer simulations”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Accelerated fast iterative shrinkage thresholding algorithms for sparsity-regularized cone-beam CT image reconstruction

Purpose: The development of iterative image reconstruction algorithms for cone-beam computed tomography (CBCT) remains an active and important research area. Even with hardware acceleration, the overwhelming majority of the available 3D iterative algorithms that implement nonsmooth regularizers remain computationally burdensome and have not been translated for routine use in time-sensitive applications such as image-guided radiation therapy (IGRT). In this work, two variants of the fast iterative shrinkage thresholding algorithm (FISTA) are proposed and investigated for accelerated iterative image reconstruction in CBCT. Methods: Algorithm acceleration was achieved by replacing the original gradient-descent step in the FISTAs by a subproblem that is solved by use of the ordered subset simultaneous algebraic reconstruction technique (OS-SART). Due to the preconditioning matrix adopted in the OS-SART method, two new weighted proximal problems were introduced and corresponding fast gradient projection-type algorithms were developed for solving them. We also provided efficient numerical implementations of the proposed algorithms that exploit the massive data parallelism of multiple graphics processing units. Results: The improved rates of convergence of the proposed algorithms were quantified in computer-simulation studies and by use of clinical projection data corresponding to an IGRT study. The accelerated FISTAs were shown to possess dramatically improved convergence properties as compared to the standard FISTAs. For example, the number of iterations to achieve a specified reconstruction error could be reduced by an order of magnitude. Volumetric images reconstructed from clinical data were produced in under 4 min. Conclusions: The FISTA achieves a quadratic convergence rate and can therefore potentially reduce the number of iterations required to produce an image of a specified image quality as compared to first-order methods. We have proposed and investigated accelerated FISTAs for use with two nonsmooth penalty functions that will lead to further reductions in image reconstruction times while preserving image quality. Moreover, with the help of a mixed sparsity-regularization, better preservation of soft-tissue structures can be potentially obtained. The algorithms were systematically evaluated by use of computer-simulated and clinical data sets.

60 APPLIED LIFE SCIENCES↗
NASA NTRSDate not supplied

Single-Scattering Properties of Melting Precipitation for GPM Passive Microwave and Radar Remote Sensing Applications

Over the past two decades, detailed computational simulations of the intricate three dimensional structures of ice-phase crystals and aggregates of those crystals have been developed. The microwave single-scattering properties of these simulated ice particles have been computed and used to improve quantitative estimates of snow rates and to better deRne the vertical structure of snow water contents in deep convective systems, as derived from satellite-borne passive microwave and/or radar remote sensing measurements. The same icephase particles have more recently been used as the starting point for simulations of melting precipitation using computational melting methods. In the current study, a heuristic melting method, as well as a physically-based melting procedure based on smoothed-particle hydrodynamics, are applied to ice particle models to describe the full evolution of the particles from dry snow to liquid drops. The discrete dipole approximation is utilized to calculate the single-scattering properties of the mixed-phase particles throughout the melting process. Then, the properties of the particles are “mapped” into simpliRed microphysical simulations of particle size spectra in the melting layers of stratiform, precipitating clouds. The bulk single-scattering properties of the melting layers and the sensitivity of their properties to modeling assumptions are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗
NASA NTRS2004

Composite Fan Blade Design for Advanced Engine Concepts

The aerodynamic and structural viability of composite fan blades of the revolutionary Exo-Skeletal engine are assessed for an advanced subsonic mission using the NASA EST/BEST computational simulation system. The Exo-Skeletal Engine (ESE) calls for the elimination of the shafts and disks completely from the engine center and the attachment of the rotor blades in spanwise compression to a rotating casing. The fan rotor overall adiabatic efficiency obtained from aerodynamic analysis is estimated at 91.6 percent. The flow is supersonic near the blade leading edge but quickly transitions into a subsonic flow without any turbulent boundary layer separation on the blade. The structural evaluation of the composite fan blade indicates that the blade would buckle at a rotor speed that is 3.5 times the design speed of 2000 rpm. The progressive damage analysis of the composite fan blade shows that ply damage is initiated at a speed of 4870 rpm while blade fracture takes place at 7640 rpm. This paper describes and discusses the results for the composite blade that are obtained from aerodynamic, displacement, stress, buckling, modal, and progressive damage analyses. It will be demonstrated that a computational simulation capability is readily available to evaluate new and revolutionary technology such as the ESE.

Abumeri, Galib H.↗
NASA NTRS1993

Experimental and numerical analysis of the wing rock characteristics of a 'wing-body-tail' configuration

Free-to-roll wind tunnel tests were conducted and a computer simulation exercise was performed in an effort to investigate in detail the mechanism of wing rock on a configuration that consisted of a highly-slender forebody and a 78 deg swept delta wing. In the wind tunnel test, the roll angle and wing surface pressures were measured during the wing rock motion. A limit cycle oscillation was observed for angles of attack between 22 deg and 30 deg. In general, the wind tunnel test confirmed that the main flow phenomena responsible for the wing-body-tail wing rock are the interactions between the forebody and the wing vortices. The variation of roll acceleration (determined from the second derivative of the roll angle time history) with roll angle clearly showed the energy balance necessary to sustain the limit cycle oscillation. Pressure measurements on the wing revealed the hysteresis of the wing rock process. First, second and nth order models for the aerodynamic damping were developed and examined with a one degree of freedom computer simulation. Very good agreement with the observed behavior from the wind tunnel was obtained.

Suarez, Carlos J.↗
DOE OSTI2021

Effects of Al:Si and (Al+Na):Si Ratios on the Properties of the International Simple Glass, Part II: Structure

High-alumina containing high-level waste (HLW) will be vitrified at the Waste Treatment Plant at the Hanford Site. The resulting glasses, high in alumina, will have distinct composition-structure-property (C-S-P) relationships compared to previously studied HLW glasses. These C-S-P relationships determine the processability and product durability of glasses and therefore must be understood. The main purpose of this study is to understand the detailed structural changes caused by Al:Si and (Al+Na):Si substitutions in a simplified nuclear waste model glass (ISG, International Simple Glass) by combining experimental structural characterizations and molecular dynamics (MD) simulations. The structures of these two series of glasses were characterized by neutron total scattering and 27Al, 23Na, 29Si, and 11B solid-state nuclear magnetic resonance (NMR) spectroscopy. Additionally, MD simulations were used to generate atomistic structural models of the borosilicate glasses and simulation results were validated by the experimental structural data. Short-range (e.g., bond distance, coordination number, etc.) and medium-range (e.g., oxygen speciation, network connectivity, polyhedral linkages) structural features of the borosilicate glasses were systematically investigated as a function of the degree of substitution. The results show that bond distance and coordination number of the cation-oxygen pairs are relatively insensitive to Al:Si and (Al+Na):Si substitutions with the exception of the B-O pair. Additionally, the Al:Si substitution results in an increase of tri-bridging oxygen species, while (Al+Na):Si substitution creates non-bridging oxygen species. Charge compensator preferences were found for Si-[NBO] (Na+), [3]B-[NBO] (Na+), [4]B (mostly Ca2+), [4]Al (nearly equally split Na+ and Ca2+), and [6]Zr (mostly Ca2+). The network former-BO-network former linkages preferences were also tabulated; Si-O-Al and Al-O-Al were preferred at the expense of lower Si-O-[3]B and [3]B-O-[3]B linkages. These results provide insights on the structural origins of property changes such as glass transition temperature caused by the substitutions, providing a basis for future improvements of theoretical and computer simulation models.

Lu, Xiaonan↗
DOE OSTI2022

Internal Model-Based Active Damping Strategy for a Back-to-Back Modular Multilevel Converter System for Advanced Grid Support: Preprint

The proposed work focuses on the possibility of achieving the reduction in the size of the interfacing filter for the grid connection for a back-to-back modular multilevel converter system with the use of third order LCL filters. To achieve the same attenuation, it is possible to reduce the size of the interfacing filters by the usage of third order filters. However, this kind of filtering comes with the limitation of having sustained oscillations due to lack of damping especially during transient changes. Therefore, in this work, a methodology has been proposed based on the principle of internal model to cater for this unwanted oscillations. The third order system is modeled inside the microcontroller with the damping enabled. The error between the output from the model and the actual are compared and the error is utilized to accomplish the active damping strategy. The proposed architecture is verified via computer simulations based on MATLAB/Simulink domain and various case study results along with their discussion is presented in this paper.

dq control↗
DOE OSTI2016

TH-CD-202-07: A Methodology for Generating Numerical Phantoms for Radiation Therapy Using Geometric Attribute Distribution Models

Purpose: A persistent challenge for the quality assessment of radiation therapy treatments (e.g. contouring accuracy) is the absence of the known, ground truth for patient data. Moreover, assessment results are often patient-dependent. Computer simulation studies utilizing numerical phantoms can be performed for quality assessment with a known ground truth. However, previously reported numerical phantoms do not include the statistical properties of inter-patient variations, as their models are based on only one patient. In addition, these models do not incorporate tumor data. In this study, a methodology was developed for generating numerical phantoms which encapsulate the statistical variations of patients within radiation therapy, including tumors. Methods: Based on previous work in contouring assessment, geometric attribute distribution (GAD) models were employed to model both the deterministic and stochastic properties of individual organs via principle component analysis. Using pre-existing radiation therapy contour data, the GAD models are trained to model the shape and centroid distributions of each organ. Then, organs with different shapes and positions can be generated by assigning statistically sound weights to the GAD model parameters. Organ contour data from 20 retrospective prostate patient cases were manually extracted and utilized to train the GAD models. As a demonstration, computer-simulated CT images of generated numerical phantoms were calculated and assessed subjectively and objectively for realism. Results: A cohort of numerical phantoms of the male human pelvis was generated. CT images were deemed realistic both subjectively and objectively in terms of image noise power spectrum. Conclusion: A methodology has been developed to generate realistic numerical anthropomorphic phantoms using pre-existing radiation therapy data. The GAD models guarantee that generated organs span the statistical distribution of observed radiation therapy patients, according to the training dataset. The methodology enables radiation therapy treatment assessment with multi-modality imaging and a known ground truth, and without patient-dependent bias.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Positioning colloids at the surfaces of cholesteric liquid crystal droplets

We report on the internal configurations of aqueous dispersions of droplets of cholesteric liquid crystals (LCs; 5–50 μm-in-diameter; comprised of 4-cyano-4'-pentylbiphenyl and 4-(1-methylheptyloxycarbonyl)phenyl-4-hexyloxybenzoate) and their influence on the positioning of surface-adsorbed colloids (0.2 or 1 μm-in-diameter polystyrene (PS)). When N = 2D/P was less than 4, where D is the droplet diameter and P is the cholesteric pitch, the droplets adopted a twisted bipolar structure (TBS) and colloids were observed to assume positions at either the poles or equator of the droplets. A statistical analysis of the distribution of locations of the colloids revealed a potential well of depth 2.7 k B T near the equator, a conclusion that was supported by computer simulations performed via the minimization of the Landau–de Gennes free energy (well depth of 7 k B T from simulation). In contrast, for N > 4, a majority of the droplets exhibited a radial spherical structure (RSS) characterized by a pair of closely spaced surface defects (angle of separation with respect to the center of the droplet θ < 5°) connected by a disclination winding to/from the droplet center, which led to the positioning of pairs of colloids with well-defined spacing at these surface defects. The separation of the pairs of surface-adsorbed colloids was colloid size-dependent, ranging from 1.11 ± 0.04 μm for 1 μm-in-diameter colloids to 1.7 ± 0.2 μm for 200 nm-in-diameter colloids. We also observed long-lived metastable configurations in which the two surface point defects were separated by much larger distances (corresponding to populations with angles of θ = 20 ± 10° and 85 ± 10° with respect to the center), and observed these pairs of defects to also position pairs of colloids. A third configuration, the diametrical spherical structure (DSS) was also observed. Consistent with the predictions of computer simulations, we found experimentally that the DSS is indeed composed of disconnected defect rings positioned along the diameter of the droplet. Finally, these results reveal that the rich palette of defects exhibited by confined cholesteric LC systems (equilibrium and metastable) provide the basis of a versatile class of templates that enable the surface positioning of colloids in ways that are not possible with achiral LC droplets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

SuperNeuro: A Fast and Scalable Simulator for Neuromorphic Computing

In many neuromorphic workflows, simulators play a vital role for important tasks such as training spiking neural networks, running neuroscience simulations, and designing, implementing, and testing neuromorphic algorithms. Currently available simulators cater to either neuroscience workflows (e.g., NEST and Brian2) or deep learning workflows (e.g., BindsNET). Problematically, the neuroscience-based simulators are slow and not very scalable, and the deep learning-based simulators do not support certain functionalities that are typical of neuromorphic workloads (e.g., synaptic delay). In this paper, we address this gap in the literature and present SuperNeuro, which is a fast and scalable simulator for neuromorphic computing capable of both homogeneous and heterogeneous simulations as well as GPU acceleration. We also present preliminary results that compare SuperNeuro to widely used neuromorphic simulators such as NEST, Brian2, and BindsNET in terms of computation times. We demonstrate that SuperNeuro can be approximately 10×--300× faster than some of the other simulators for small sparse networks. On large sparse and large dense networks, SuperNeuro can be approximately 2.2×--3.4× faster than the other simulators, respectively.

Date, Prasanna↗
DOE OSTI2026

Data for "Photoinduced Chemomimetic Biocatalysis for Enantioselective Intermolecular Radical Conjugate Addition"

Exploiting nature’s catalysts for non-natural transformations that are inaccessible to chemocatalysis is highly desirable but challenging. On the one hand, the widespread nicotinamide-dependent oxidoreductases have not been utilized for single-electron-transfer-induced bimolecular cross-couplings; on the other, the addition of catalytic asymmetric radical conjugate to terminal alkenes remains a challenge owing to strong racemic background reaction and unselective termination of prochiral radical species. Here we report a chemomimetic biocatalysitic approach for construction of alpha-carbonyl stereocentres via an unnatural intermolecular conjugate addition of N-(acyloxy)phthalimides-derived radicals with acceptor-substituted terminal alkenes, by combination of visible-light excitation and nicotinamide-dependent ketoreductases (KREDs). Based on protein crystal structure, we engineered KREDs via a semi-rational mutagenesis strategy to improve reaction outcomes with a small and high-quality variants library. Mechanistic investigations combining wet experiments, crystallographic studies and computational simulations demonstrate that the repurposed biocatalyst can suppress racemic background reaction and unselected side reactions, yielding enantioselectivity that is challenging to achieve by chemocatalysis.

Catalysis↗
NASA NTRS1990

Influence of linear profile modification and loading conditions on the dynamic tooth load and stress of high contact ratio gears

A computer simulation for the dynamic response of high-contact-ratio spur gear transmissions is presented. High contact ratio gears have the potential to produce lower dynamic tooth loads and minimum root stress but they can be sensitive to tooth profile errors. The analysis presented examines various profile modifications under realistic loading conditions. The effect of these modifications on the dynamic load (force) between mating gear teeth and the dynamic root stress is presented. Since the contact stress is dependent on the dynamic load, minimizing dynamic loads will also minimize contact stresses. It is shown that the combination of profile modification and the applied load (torque) carried by a gear system has a significant influence on gear dynamics. The ideal modification at one value of applied load will not be the best solution for a different load. High-contact-ratio gears were found to require less modification than standard low-contact-ratio gears. High-contact-ratio gears are more adversely affected by excess modification than by under modification. In addition, the optimal profile modification required to minimize the dynamic load (hence the contact stress) on a gear tooth differs from the optimal modification required to minimize the dynamic root (bending) stress. Computer simulation can help find the design tradeoffs to determine the best profile modification to satisfy the conflicting constraints of minimizing both the load and root stress in gears which must operate over a range of applied loads.

Lee, Chinwai↗
NASA NTRSDate not supplied

Validation of Kestrel IDDES Simulations for SLS Transition Analysis

Complex computational simulations are needed to support the Space Launch System (SLS) program, and the fidelity of the computational results must be defended. More specifically, a number of databases are produced for the transition phase of flight, which occurs after the rocket clears the tower but before reaching transonic speeds. In an effort to reduce computational uncertainty, many of the computational parameters were altered to determine the sensitivity of the results to the value of the parameter and then updating the best practice procedures. The baseline routines were developed over many years through the maturation of the SLS program, and this paper delivers a detailed discussion of these baseline results. This section is followed by demonstration of perturbing some of the more significant parameters, including time step and turbulence model, and concluded with a summary of the herein determined best practices for such analysis.

CFD↗
NASA NTRS1975

Long period nodal motion of sun synchronous orbits

An approximative model is formulated for assessing these perturbations that significantly affect long term modal motion of sun synchronous orbits. Computer simulations with several independent computer programs consider zonal and tesseral gravitational harmonics, third body gravitational disturbances induced by the sun and the moon, and atmospheric drag. A pendulum model consisting of evenzonal harmonics through order 4 and solar gravity dominated nodal motion approximation. This pendulum motion results from solar gravity inducing an inclination oscillation which couples into the nodal precession induced by the earth's oblateness. The pendulum model correlated well with simulations observed flight data.

Duck, K. I.↗
NASA NTRS2010

A Review of Discrete Element Method (DEM) Particle Shapes and Size Distributions for Lunar Soil

As part of ongoing efforts to develop models of lunar soil mechanics, this report reviews two topics that are important to discrete element method (DEM) modeling the behavior of soils (such as lunar soils): (1) methods of modeling particle shapes and (2) analytical representations of particle size distribution. The choice of particle shape complexity is driven primarily by opposing tradeoffs with total number of particles, computer memory, and total simulation computer processing time. The choice is also dependent on available DEM software capabilities. For example, PFC2D/PFC3D and EDEM support clustering of spheres; MIMES incorporates superquadric particle shapes; and BLOKS3D provides polyhedra shapes. Most commercial and custom DEM software supports some type of complex particle shape beyond the standard sphere. Convex polyhedra, clusters of spheres and single parametric particle shapes such as the ellipsoid, polyellipsoid, and superquadric, are all motivated by the desire to introduce asymmetry into the particle shape, as well as edges and corners, in order to better simulate actual granular particle shapes and behavior. An empirical particle size distribution (PSD) formula is shown to fit desert sand data from Bagnold. Particle size data of JSC-1a obtained from a fine particle analyzer at the NASA Kennedy Space Center is also fitted to a similar empirical PSD function.

Lane, John E.↗
DOE OSTI2019

Preface to the special volume on the third Sandia Fracture Challenge

The mounting reliance on computational simulations to predict all aspects of the lifecycle of a mechanical system, from fabrication to failure, has prompted the mechanics community to selfassess its abilities to perform those predictions. Benchmark problems in mechanics that compare simulations that use different computational approaches with experiments have sprung up lately, including the NIST AM-Bench looking at additively manufactured (AM) materials (https://www.nist.gov/ambench),the Contact-Mechanics Challenge (Miiser, 2017) considering adhesion between two nominally flat surfaces, Numisheet providing semiannual benchmarking activities in sheet metal forming (http://numisheet2018.org),and the Sandia Fracture Challenge (SFC) (Boyce, 2014 and Boyce, 2016) investigating ductile failure. The previous SFCs have shown that progress has been made in computations of ductile failure, but improvements still can be made, hence the third Sandia Fracture Challenge (SFC3), the subject of this Special Volume. The most recent installment of SFC is building on previous successes and tackling the difficult problem of fracture in an AM 316L stainless steel structure.

42 ENGINEERING↗
DOE OSTI2016

Program optimizations: The interplay between power, performance, and energy

Practical considerations for future supercomputer designs will impose limits on both instantaneous power consumption and total energy consumption. Working within these constraints while providing the maximum possible performance, application developers will need to optimize their code for speed alongside power and energy concerns. This paper analyzes the effectiveness of several code optimizations including loop fusion, data structure transformations, and global allocations. A per component measurement and analysis of different architectures is performed, enabling the examination of code optimizations on different compute subsystems. Using an explicit hydrodynamics proxy application from the U.S. Department of Energy, LULESH, we show how code optimizations impact different computational phases of the simulation. This provides insight for simulation developers into the best optimizations to use during particular simulation compute phases when optimizing code for future supercomputing platforms. Here, we examine and contrast both x86 and Blue Gene architectures with respect to these optimizations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Sandia National Laboratories Advanced Simulation and Computing (ASC) Appraisal Method for the Implementation of the ASC Software Quality Engineering Practices (V.2.0)

This document provides a guide to the process of conducting software appraisals under the Sandia National Laboratories (SNL) ASC Program. The goal of this document is to describe a common methodology for planning, conducting, and reporting results of software appraisals thereby enabling: 1) Identification of improvements in implementation of the software quality engineering (SQE) practices identified in the ASC Software Quality Plan across the ASC Program against objective baselines. 2) Feedback from project teams on SQE opportunities for improvement. 3) Identification of strengths and opportunities for improvement for individual project teams. 4) Guidance to the ASC Program on the focus of future SQE activities. Document contents include process descriptions, templates to promote consistent conduct of appraisals, and an explanation of the relationship of this procedure to the SNL ASC software program. The activities described by this document are under the oversight of SNL ASC Program Director and SNL ASC Program Manager.

42 ENGINEERING↗