Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational methods”

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
NASA NTRS1989

LaRC computational dynamics overview

Present research centers on the development of advanced computational methods for transient simulation analyses. Aircraft, launch vehicles and space structure components are potential applications, but primary focus is presently on large space structures. There are both in-house and out-of-house activities. The in-house activity centers around the development of a multibody simulation tool for truss-like structures called LATDYN for Large Angle Transient DYNamics. Multibody analysis involves articulation of structural components as well as robotic maneuvers. These items are necessary for construction (erection or deployment) of large space structures in orbit and the carrying out of certain operations on board the space station. Thus, part of the in-house activity involves the development of methods which treat the changing mass, stiffness and constraints associated with articulating systems. The out-of-house activity involves subcycling, development of large deformation/motion beam formulation, constraint stabilization and direct time integration transient algorithms in parallel computing.

Husner, J. M.↗
NASA NTRS1998

Dynamic Load Balancing for Adaptive Unstructured Grids

Dynamic mesh adaptation on unstructured grids is a powerful tool for computing unsteady three-dimensional problems that require grid modifications to efficiently resolve solution features. By locally refining and coarsening the mesh to capture phenomena of interest, such procedures make standard computational methods more cost effective. Highly refined meshes are required to accurately capture shock waves, contact discontinuities, vortices, and shear layers in fluid flow problems. Adaptive meshes have also proved to be useful in several other areas of computational science and engineering like computer vision and graphics, semiconductor device modeling, and structural mechanics. Local mesh adaptation provides the opportunity to obtain solutions that are comparable to those obtained on globally-refined grids but at a much lower cost. Additional information is contained in the original extended abstract.

Biswas, Rupak↗
DOE OSTI2018

Overcoming the Memory Bottleneck in Auxiliary Field Quantum Monte Carlo Simulations with Interpolative Separable Density Fitting

We investigate the use of interpolative separable density fitting (ISDF) as a means to reduce the memory bottleneck in auxiliary field quantum Monte Carlo (AFQMC) simulations of real materials in Gaussian basis sets. Here, we find that ISDF can reduce the memory scaling of AFQMC simulations from O(M 4 ) to O(M 2 ) . We test these developments by computing the structural properties of carbon in the diamond phase, comparing to results from existing computational methods and experiment. (LLNL-JRNL-759024)

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

Method and computer-readable model for additively manufacturing ducting arrangement for a gas turbine engine

Method and computer-readable model for additively manufacturing a ducting arrangement (10) for a gas turbine engine are provided. Ducting arrangement (10) may include a duct (18) to be fluidly coupled to receive a cross-flow of combustion gases from a main combustion stage. Duct (18) includes a duct segment (23) with an expanding cross-sectional area (24) where one or more injector assemblies (26) are disposed. Injector assembly (26) includes one or more reactant-guiding structures (27) arranged to deliver a flow of reactants to be mixed with the cross-flow of combustion gases. The ducting arrangement is effective to reduce total pressure loss while providing an effective level of mixing of the injected reactants with the passing cross-flow. Respective duct components or the entire ducting arrangement may be formed as a unitized structure, such as a single piece using a rapid manufacturing technology, such as 3D Printing/Additive Manufacturing (AM) technologies.

Chen, Yuntao↗
DOE OSTI2025

Silver diamine fluoride differentially affects dentin and hypomineralized enamel permeabilities

OBJECTIVES: To investigate the physicochemical effect of silver diamine fluoride (SDF) by correlating permeability with mineral density and elemental composition of hypomineralized enamel and carious dentin. METHODS: Enamel and dentin from human carious primary teeth with and without SDF treatment in-vivo, and hypomineralized enamel from permanent molars with and without SDF treatment in-vitro were scanned using micro X-ray computed tomography. Spatial maps of biometals (calcium, zinc), phosphorus, and silver were generated using X-ray fluorescence microprobe. Permeabilities were computed using Porous Microstructure Analysis software. RESULTS: The intrinsic permeability of SDF-treated carious dentin was 14.3 % lower than untreated sound dentin (6.39e-15 ± 3.01e-15 m² vs 7.46e-15 ± 1.82e-15 m²; P < 0.0001), while untreated carious dentin was 98.4 % higher (1.48e-14 ± 7.11e-15 m²; P < 0.0001). SDF-treated and untreated transparent dentin showed similar reduced permeabilities (75.6 % and 78.4 % lower than untreated sound dentin, respectively; P = 0.93). Severely hypomineralized enamel showed permeability reaching 108.1 % of adjacent sound dentin (5.71e-15 ± 2.04e-15 m² vs 5.28e-15 ± 1.30e-15 m²; P = 0.1409) and was significantly higher than mildly hypomineralized enamel (1.39e-15 ± 1.04e-15 m²; P < 0.0001). SDF treatment did not significantly impact the permeability of severely hypomineralized enamel (12.4 % reduction; P = 0.07). Principal component regression identified Zn level as a significant effector of tissue permeabilities in carious primary teeth (P < 0.0001). SIGNIFICANCE: This study introduces a computational method to measure dental tissue permeability, and demonstrates that SDF significantly reduces permeability in carious dentin but not intact hypomineralized enamel. The study reveals biometal Zn localization can alter dentin and enamel permeabilities, providing new insights into pathobiological mechanisms underlying caries and hypomineralization.

Chou, Conrad↗
DOE OSTI2019

Hydrodynamic Analysis of a Suspended Cylinder Under Regular Wave Loading Based on Computational Fluid Dynamics

An investigation into the computation of hydrodynamic loads on a suspended cylinder in regular waves is presented. The primary goal was to perform a three-way validation of the loads between experimental measurements and simulations from two computational methods. Experimental measurements of the longitudinal in-line force on a cylinder suspended at a fixed position were available from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project, Phase Ia. These measurements were compared to computational fluid dynamics (CFD) simulations based on the solution of Reynolds-averaged Navier-Stokes (RANS) equations, as implemented in STAR-CCM+. The study encompassed a sensitivity analysis of the loads computed in STAR-CCM+ based on wave modeling, boundary conditions, turbulence modeling, and spatial and temporal discretization. The analysis was supplemented by results generated with the offshore wind turbine engineering software OpenFAST, based on a hybrid combination of second-order potential flow and viscous drag from Morison's equation. The focus of the investigation was on the assessment of the accuracy of the computation of first- and higher-order hydrodynamic loads. Substantial differences were observed in the numerical prediction of the second and third harmonic force contribution. Local flow field analysis with CFD was applied to study the physics of wave run-up and diffraction dynamics to identify the causes.

16 TIDAL AND WAVE POWER↗
NASA NTRS1984

Calculation of boundary layers of oscillating airfoils

A two-point finite difference unsteady laminar and turbulent boundary layer computational method was used to investigate the properties of the flow around an airfoil (NACA 0012) oscillating through angles of attack up to 18 degrees, for reduced frequencies of 0.01 and 0.20. The unsteady potential flow was determined using the unsteady potential flow method of Geissler. The influence of transition location on stal behavior was investigated, using both experimentally determined transition information, and transition located at the pressure peak; the results show the need for viscous-inviscid interaction in future computation of such flows.

Cebeci, T.↗
NASA NTRS1988

Generation of dynamo magnetic fields in protoplanetary and other astrophysical accretion disks

A computational method for treating the generation of dynamo magnetic fields in astrophysical disks is presented. The numerical difficulty of handling the boundary condition at infinity in the cylindrical disk geometry is overcome by embedding the disk in a spherical computational space and matching the solutions to analytically tractable spherical functions in the surrounding space. The lowest lying dynamo normal modes for a 'thick' astrophysical disk are calculated. The generated modes found are all oscillatory and spatially localized. Tha potential implications of the results for the properties of dynamo magnetic fields in real astrophysical disks are discussed.

Stepinski, T. F.↗
NASA NTRS1993

Computational/experimental studies of isolated, single component droplet combustion

Isolated droplet combustion processes have been the subject of extensive experimental and theoretical investigations for nearly 40 years. The gross features of droplet burning are qualitatively embodied by simple theories and are relatively well understood. However, there remain significant aspects of droplet burning, particularly its dynamics, for which additional basic knowledge is needed for thorough interpretations and quantitative explanations of transient phenomena. Spherically-symmetric droplet combustion, which can only be approximated under conditions of both low Reynolds and Grashof numbers, represents the simplest geometrical configuration in which to study the coupled chemical/transport processes inherent within non-premixed flames. The research summarized here, concerns recent results on isolated, single component, droplet combustion under microgravity conditions, a program pursued jointly with F.A. Williams of the University of California, San Diego. The overall program involves developing and applying experimental methods to study the burning of isolated, single component droplets, in various atmospheres, primarily at atmospheric pressure and below, in both drop towers and aboard space-based platforms such as the Space Shuttle or Space Station. Both computational methods and asymptotic methods, the latter pursued mainly at UCSD, are used in developing the experimental test matrix, in analyzing results, and for extending theoretical understanding. Methanol, and the normal alkanes, n-heptane, and n-decane, have been selected as test fuels to study time-dependent droplet burning phenomena. The following sections summarizes the Princeton efforts on this program, describe work in progress, and briefly delineate future research directions.

Dryer, Frederick L.↗
NASA NTRS2024

On the Role of Friction and Particle Size Distribution in Granular Packings

Packing of particles in a disordered arrangement has tremendous significance in both condensed matter physics and engineering applications. The last three decades have seen remarkable progress in our understanding of the physics of granular packings that has been largely facilitated by a rapid growth in the power of modern computers. Although granular packings are ubiquitous in diverse natural settings, from clogging of powders in hoppers to the crowding of living cells, a significant motivation for modeling granular packings has emerged from a proposal that the transition from a fluid-like state to a solid-like state of a granular material upon increasing volume fraction, called jamming, is intimately related to the origins of glass transition in thermal systems. As such, a majority of modeling efforts have focused on the jamming behavior of an idealized granular material: frictionless, monodisperse sphere. While such studies have illuminated the rich physics of jamming, granular materials in nature and engineering practice are rarely frictionless or monodisperse. The analogous research on the packing of these ‘real-world’ granular materials is still not fully developed. Besides requiring the exploration of a huge parameter space, three key computational considerations have inhibited their modeling: (i) traditional computational methods are not adept at simulating mechanically-stable packings of frictional particles near the jamming transition; (ii) standard algorithms of contact detection in discrete element methods are impractical to simulate granular packings with a wide distribution of particle sizes; (iii) a lack of well-established contact mechanics models of friction that can accurately reproduce experimental data. This chapter will review the latest computational advances to simulate the jamming of size-dispersed frictional particles, and describe the rich microstructural diversity that emerges in their packings.

granular↗
DOE OSTI2019

Hydrodynamic Analysis of a Suspended Cylinder Under Regular Wave Loading Based on Computational Fluid Dynamics: Preprint

An investigation into the computation of hydrodynamic loads on a suspended cylinder in regular waves is presented. The primary goal was to perform a three-way-validation of the loads between experimental measurements and simulations from two computational methods. Experimental measurements of the longitudinal inline force on a cylinder suspended at a fixed po- sition were available from the Offshore Code Comparison Col- laboration, Continued, with Correlation (OC5) project, Phase Ia. These measurements were compared to Computational Fluid Dynamics (CFD) simulations based on the solution of Reynolds- averaged Navier-Stokes (RANS) equations, as implemented in STAR-CCM plus. The study encompassed a sensitivity analysis of the loads computed in STAR-CCM plus based on wave modeling, boundary conditions, turbulence modeling, and spatial and tem- poral discretization. The analysis was supplemented by results generated with the offshore wind turbine engineering software OpenFAST, based on a hybrid combination of second-order po- tential flow and viscous drag fromMorison's equation. The focus of the investigation was on the assessment of the accuracy of the computation of first- and higher-order hydrodynamic loads. Sub- stantial differences were observed in the numerical prediction of the second and third harmonic force contribution. Local flow field analysis with CFD was applied to study the physics of wave run-up and diffraction dynamics to identify the causes.

16 TIDAL AND WAVE POWER↗
DOE OSTI2019

Long non-coding subgenomic flavivirus RNAs have extended 3D structures and are flexible in solution

Most mosquito-borne flaviviruses, including Zika virus (ZIKV), Dengue virus (DENV), and West Nile virus (WNV), produce long non-coding subgenomic RNAs (sfRNAs) in infected cells that link to pathogenicity and immune evasion. Until now, the structural characterization of these lncRNAs remains limited. Here, we studied the 3D structures of individual and combined subdomains of sfRNAs, and visualized the accessible 3D conformational spaces of complete sfRNAs from DENV2, ZIKV, and WNV by small angle X-ray scattering (SAXS) and computational modeling. The individual xrRNA1s and xrRNA2s adopt similar structures in solution as the crystal structure of ZIKV xrRNA1, and all xrRNA1-2s form compact structures with reduced flexibility. While the DB12 of DENV2 is extended, the DB12s of ZIKV and WNV are compact due to the formation of intertwined double pseudoknots. All 3' stem-loops (3'SLs) share similar rod-like structures. Complete sfRNAs are extended and sample a large conformational space in solution. Our work not only provides structural insight into the function of flavivirus sfRNAs, but also highlights strategies of visualizing other lncRNAs in solution by SAXS and computational methods.

3D structure↗
NASA NTRS1985

Boundary layers on oscillating airfoils

A two-point finite-difference unsteady laminar and turbulent boundary-layer computational method has been used to investigate the properties of the flow around an airfoil (NACA 0012) oscillating through angles of attack up to 18 degrees, and for reduced frequencies of 0.01 and 0.20. The unsteady potential flow was determined using the method of Geissler. The influence of transition location on stall behavior has been investigated, using both experimentally determined transition information, and transition located at the pressure peak; the results show the need for viscous-inviscid interaction in future computations of such flows.

Carr, L. W.↗
DOE OSTI2019

On Xenon Fission Product Poisoning

We all appreciate the special role Xenon-135 plays in reactor transients. Unless there is sufficient reactivity override, on shutdown a conventional reactor cannot immediately restart because of the buildup of strongly absorbing Xenon from Iodine decay and the absence of neutron flux to continue its elimination. We are also aware of how the flux imbalance in a large reactor core gives rise to Xenon induced power precession around the core in unpredictable ways. In addition, with regard to computational methods, we lack readily accessible numerical algorithms enabling reliable characterization of Xenon oscillations with simple models to explain this fundamentally important phenomenon. Nowhere is this more apparent than in the classroom. In this presentation, we first consider a new 12-element comprehensive decay to model Xenon and its corresponding numerical implementation and verification to establish a predictive model of Xenon buildup in a constant flux. We then follow with a reduced self-consistent flux model in an infinite reactor- a model not found in textbooks. This presentation should appeal to those interested in developing high order numerical methods for reactor transients.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Magnetic structure of UO 2 and NpO 2 by first-principle methods

Here, the magnetic structure of the actinide dioxides (AnO 2 ) remains a field of intense research. A low-temperature experimental investigation of the magnetic ground-state is complicated by thermal energy released from the radioactive decay of the actinide nuclei. To establish the magnetic ground-state, we have employed high-accuracy computational methods to systematically probe different magnetic structures. A transverse 1k antiferromagnetic ground-state with Fmmm (No. 69) crystal symmetry has been established for UO 2 , whereas a ferromagnetic (111) ground-state with R¯3 with combining macron]m (No. 166) has been established for NpO 2 . Band structure calculations have been performed to analyse these results.

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

Validation of Floating Node Method Using Three-Point Bend Doubler Under Quasi-Static Loading

The NASA Advanced Composite Project (ACP), an industry/government/university partnership, has embarked upon the task of developing technology that can aid in reducing the time line for structural certification of aircraft composite parts using a combination of technologies, one of which is high fidelity damage progression computational methods. Phase II of this project included a task for validating an approach based on the Floating Node Method combined with Directional Cohesive Elements (FNM-DCZE). This paper discusses predicted damage onset and growth in a three-point bend doubler specimen compared to experimental results. Sensitivity of the simulations to mesh refinement as well as key material properties and thermal effects are studied and reported. Overall, qualitative results suggest the main aspects of the damage progression have been captured, with the simulated damage morphology and sequence of events resembling closely what was observed experimentally. Quantitatively, the first load-peak is predicted. However, the re-loading observed in the experiments, after the first load peak, is not captured numerically, suggesting further investigation may be worth pursuing.

Selvarathinam, Alex S.↗
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Nonlinear Aerodynamics and the Design of Wing Tips

The analysis and design of wing tips for fixed wing and rotary wing aircraft still remains part art, part science. Although the design of airfoil sections and basic planform geometry is well developed, the tip regions require more detailed consideration. This is important because of the strong impact of wing tip flow on wing drag; although the tip region constitutes a small portion of the wing, its effect on the drag can be significant. The induced drag of a wing is, for a given lift and speed, inversely proportional to the square of the wing span. Concepts are proposed as a means of reducing drag. Modern computational methods provide a tool for studying these issues in greater detail. The purpose of the current research program is to improve the understanding of the fundamental issues involved in the design of wing tips and to develop the range of computational and experimental tools needed for further study of these ideas.

Kroo, Ilan↗