Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational simulation”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Xenon Gas Separation and Storage using Metal-Organic Frameworks

The global demand for Xe, a noble gas with applications in electronics, lighting, and the medical industry, is expected to rise significantly over the coming decades. However, the low abundance of Xe in the earth’s atmosphere and the costly cryogenic distillation process that is used to obtain Xe commercially via air separation have limited the scale of applications of Xe. A physisorption-based separation using porous materials may be a viable and cost-effective alternative to cryogenic distillation. In particular, metal-organic frameworks (MOFs) have shown promise as highly Xe-selective, porous solids. In this review, we detail the recent advances of MOFs as adsorbents for noble gas adsorption/separation and the role of computer simulation in finding optimal materials for Xe adsorption.

MOFs, Xenon, Krypton↗
DOE OSTI2016

Temperature of maximum density and excess thermodynamics of aqueous mixtures of methanol

In this work, we present a study of representative excess thermodynamic properties of aqueous mixtures of methanol over the complete concentration range, based on extensive computer simulation calculations. In addition to test various existing united atom model potentials, we have developed a new force-field which accurately reproduces the excess thermodynamics of this system. Moreover, we have paid particular attention to the behavior of the temperature of maximum density (TMD) in dilute methanol mixtures. The presence of a temperature of maximum density is one of the essential anomalies exhibited by water. This anomalous behavior is modified in a non-monotonous fashion by the presence of fully miscible solutes that partly disrupt the hydrogen bond network of water, such as methanol (and other short chain alcohols). In order to obtain a better insight into the phenomenology of the changes in the TMD of water induced by small amounts of methanol, we have performed a new series of experimental measurements and computer simulations using various force fields. We observe that none of the force-fields tested capture the non-monotonous concentration dependence of the TMD for highly diluted methanol solutions.

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

On the nature of the molecular ordering of water in aqueous DMSO mixtures

Computer simulation studies of aqueous dimethyl sulfoxyde (DMSO) mixtures show micro-heterogeneous structures, just like aqueous alcohol mixtures. However, there is a marked difference in the aggregate structure of water between the two types of systems. While water molecules form multiconnected globular clusters in alcohols, we report herein that the typical water aggregates in aqueous DMSO mixtures are linear, favouring a 2 hydrogen bond structure per water molecule, and for all DMSO mole fractions ranging from 0.1 to 0.9. This linear-aggregate structure produces a particular signature in the water site-site structure factors, in the form of a pre-peak at k ≈ 0.2–0.8 Å{sup −1}, depending on DMSO concentration. This pre-peak is either absent in other aqueous mixtures, such as aqueous methanol mixtures, or very difficult to see through computer simulations, such as in aqueous-t-butanol mixtures. This difference in the topology of the aggregates explains why the Kirkwood-Buff integrals of aqueous-DMSO mixture look nearly ideal, in contrast with those of aqueous alcohol mixtures, suggesting a connection between the shape of the water aggregates, its fluctuations, and the concentration fluctuations. In order to further study this discrepancy between aqueous DMSO and aqueous alcohol mixture, two models of pseudo-DMSO are introduced, where the size of the sulfur atom is increased by a factor 1.6 and 1.7, respectively, hence increasing the hydrophobicity of the molecule. The study shows that these mixtures become closer to the emulsion type seen in aqueous alcohol mixtures, with more globular clustering of the water molecules, long range domain oscillations in the water-water correlations and increased water-water Kirkwood-Buff integrals. It demonstrates that the local ordering of the water molecules is influenced by the nature of the solute molecules, with very different consequences for structural properties and related thermodynamic quantities. This study illustrates the unique plasticity of water in presence of different types of solutes.

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

Geometrical structures for radiation biology research as implemented in the TOPAS-nBio toolkit

Computational simulations, such as Monte Carlo track structure simulations, offer a powerful tool for quantitatively investigating radiation interactions within cells. Additionally, the modelling of the spatial distribution of energy deposition events as well as diffusion of chemical free radical species, within realistic biological geometries, can help provide a comprehensive understanding of the effects of radiation on cells. Track structure simulations, however, generally require advanced computing skills to implement. The TOPAS-nBio toolkit, an extension to TOPAS (TOol for PArticle Simulation), aims to provide users with a comprehensive framework for radiobiology simulations, without the need for advanced computing skills. This includes providing users with an extensive library of advanced, realistic, biological geometries ranging from the micrometer scale (e.g. cells and organelles) down to the nanometer scale (e.g. DNA molecules and proteins). Finally, were we present the geometries available in TOPAS-nBio.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Best Practices for Body Temperature Measurement with Infrared Thermography: External Factors Affecting Accuracy

Infrared thermographs (IRTs) are commonly used during disease pandemics to screen individuals with elevated body temperature (EBT). To address the limited research on external factors affecting IRT accuracy, we conducted benchtop measurements and computer simulations with two IRTs, with or without an external temperature reference source (ETRS) for temperature compensation. The combination of an IRT and an ETRS forms a screening thermograph (ST). We investigated the effects of viewing angle (θ, 0–75°), ETRS set temperature (T ETRS , 30–40 °C), ambient temperature (T atm , 18–32 °C), relative humidity (RH, 15–80%), and working distance (d, 0.4–2.8 m). We discovered that STs exhibited higher accuracy compared to IRTs alone. Across the tested ranges of T atm and RH, both IRTs exhibited absolute measurement errors of less than 0.97 °C, while both STs maintained absolute measurement errors of less than 0.12 °C. The optimal T ETRS for EBT detection was 36–37 °C. When θ was below 30°, the two STs underestimated calibration source (CS) temperature (T CS ) of less than 0.05 °C. The computer simulations showed absolute temperature differences of up to 0.28 °C and 0.04 °C between estimated and theoretical temperatures for IRTs and STs, respectively, considering d of 0.2–3.0 m, T atm of 15–35 °C, and RH of 5–95%. The results highlight the importance of precise calibration and environmental control for reliable temperature readings and suggest proper ranges for these factors, aiming to enhance current standard documents and best practice guidelines. These insights enhance our understanding of IRT performance and their sensitivity to various factors, thereby facilitating the development of best practices for accurate EBT measurement.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Importance of hydrophobic traps for proton diffusion in lyotropic liquid crystals

The diffusion of protons in self-assembled systems is potentially important for the design of efficient proton exchange membranes. In this work, we study proton dynamics in a low-water content, lamellar phase of a sodium-carboxylate gemini surfactant/water system using computer simulations. The hopping of protons via the Grotthuss mechanism is explicitly allowed through the multi-state empirical valence bond method. We find that the hydronium ion is trapped on the hydrophobic side of the surfactant-water interface, and proton diffusion then proceeds by hopping between surface sites. The importance of hydrophobic traps is surprising because one would expect the hydronium ions to be trapped at the charged headgroups. The physics illustrated in this system should be relevant to the proton dynamics in other amphiphilic membrane systems, whenever there exist exposed hydrophobic surface regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2020

Free-Space Optical Communication for CubeSats in Low Lunar Orbit (LLO)

The NASA ARTEMIS Program will include LunaNet, a highly extensible, open architecture, lunar communications and navigation network. A constellation of CubeSats in Low Lunar Orbit (LLO), 100 km, could form an optical communications and navigation network as part of LunaNet, with terminals on the lunar surface, including mobile ones such as with astronauts and rovers. The proposed CubeSat nodes should provide data relay and navigational aid services. The proposed effort herein is to develop a fine pointing capability for laser beam pointing to augment body pointing by CubeSats. Body pointing was used by Aerospace Corporation for the CubeSats in LEO in NASA’s Optical Communications and Sensors Demonstration (OCSD) program [1]. Previously, this fine pointing capability was computer simulated for the OCSD program [2,3]. With fine pointing, the spot size on the Earth was reduced by a factor of eight with a reduction in laser output power by a factor of sixty-four, thereby mitigating the thermal load challenge on the CubeSats. The same reductions in spot size and laser output power can be achieved for CubeSats in LLO. A new method is described for optical data transmissions from satellites, which uses laser arrays for laser beam pointing. It combines a lens system and an array of vertical-cavity, surface-emitting lasers and photodetectors, an VCSEL/Photodetector Array, (both mature technologies), in a novel way. This system is applied to CubeSats in low lunar orbit, (LLO), which use body pointing. Also, It may be able to replace current architectures which use dynamical systems, (i.e., moving parts) to point the laser, and which may also use vibration isolation platforms. The computer simulations used the optics code, OpticStudio, from Zemax, LLC, which has the capabilities to model the laser source and diffraction effects from wave optics. These capabilities make it possible to model laser beam propagation over long space communication distances.

Goorjian, P. M.↗
DOE OSTI2015

Electron-Beam Dynamics for an Advanced Flash-Radiography Accelerator

Beam dynamics issues were assessed for a new linear induction electron accelerator being designed for multipulse flash radiography of large explosively driven hydrodynamic experiments. Special attention was paid to equilibrium beam transport, possible emittance growth, and beam stability. Especially problematic would be high-frequency beam instabilities that could blur individual radiographic source spots, low-frequency beam motion that could cause pulse-to-pulse spot displacement, and emittance growth that could enlarge the source spots. Furthermore, beam physics issues were examined through theoretical analysis and computer simulations, including particle-in-cell codes. Beam instabilities investigated included beam breakup, image displacement, diocotron, parametric envelope, ion hose, and the resistive wall instability. The beam corkscrew motion and emittance growth from beam mismatch were also studied. It was concluded that a beam with radiographic quality equivalent to the present accelerators at Los Alamos National Laboratory will result if the same engineering standards and construction details are upheld.

43 PARTICLE ACCELERATORS↗
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↗