Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing continuum”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Aerodynamic Database Development for Mars Smart Lander Vehicle Configurations

An aerodynamic database has been generated for the Mars Smart Lander Shelf-All configuration using computational fluid dynamics (CFD) simulations. Three different CFD codes, USM3D and FELISA, based on unstructured grid technology and LAURA, an established and validated structured CFD code, were used. As part of this database development, the results for the Mars continuum were validated with experimental data and comparisons made where applicable. The validation of USM3D and LAURA with the Unitary experimental data, the use of intermediate LAURA check analyses, as well as the validation of FELISA with the Mach 6 CF(sub 4) experimental data provided a higher confidence in the ability for CFD to provide aerodynamic data in order to determine the static trim characteristics for longitudinal stability. The analyses of the noncontinuum regime showed the existence of multiple trim angles of attack that can be unstable or stable trim points. This information is needed to design guidance controller throughout the trajectory.

Bobskill, Glenn J.↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A comprehensive study on three typical photoacid generators using photoelectron spectroscopy and ab initio calculations

Conducting a comprehensive molecular-level evaluation of a photoacid generator (PAG) and its subsequent impact on lithography performance can facilitate the rational design of a promising 193 nm photoresist tailored to specific requirements. In this study, we integrated spectroscopy and computational techniques to meticulously investigate the pivotal factors of three prototypical PAG anions, p-toluenesulfonate (pTS - ), 2-(trifluoromethyl)benzene-1-sulfonate (TFMBS - ), and triflate (TF - ), in the lithography process. Our findings reveal a significant redshift in the absorption spectra caused by specific PAG anions, attributed to their involvement in electronic transition processes, thereby enhancing the transparency of the standard PAG cation, triphenylsulfonium (TPS + ), particularly at ~193 nm. Furthermore, the electronic stability of PAG anions can be enhanced by solvent effects with varying degrees of strength. Here we observed the lowest vertical detachment energy of 6.6 eV of pTS - in PGMEA solution based on the polarizable continuum model, which prevents anion loss at 193 nm lithography. In addition, our findings indicate gas-phase proton affinity values of 316.4 kcal/mol for pTS - , 308.1 kcal/mol for TFMBS - , and 303.2 kcal/mol for TF - , which suggest the increasing acidity strength, yet even the weakest acid pTS - is still stronger than strong acid HBr. The photolysis of TPS + -based PAG, TPS + ·pTS - , generated an excited state leading to homolysis bond cleavage with the lowest reaction energy of 83 kcal/mol. Overall, the PAG anion pTS - displayed moderate acidity, possessed the lowest photolysis reaction energy, and demonstrated an appropriate redshift. These properties collectively render it a promising candidate for an effective acid producer.

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

A predictive modeling tool for damage analysis and design of hydrogen storage composite pressure vessels

In this study, a predictive modeling tool is developed for damage analysis and design of hydrogen (H 2 ) storage composite pressure vessels. It integrates micromechanics of matrix cracking into a continuum damage mechanics (CDM) description for damage evolution, and three-dimensional (3D) finite element (FE) modeling of the vessel structural response. At the scale of the composite layer (mesoscale), the temperature-dependent stiffness reduction law in terms of the damage variable for transverse matrix cracking is computed using an Eshelby-Mori-Tanaka approach for the initial composite thermoelastic properties and a self-consistent model for the stiffness reduction as a function of the damage variable. While transverse matrix cracking obeying a damage evolution relation can progressively evolve from an initiation to a saturation state, fiber failure is predicted by a micromechanical fiber rupture criterion that accounts for the fiber strength and matrix stress. The implementation of this integrated multiscale modeling model into a 3D FE formulation enables damage analysis and design of H 2 storage composite pressure vessels. The developed tool is illustrated through 3D damage analyses of a cryogenically compressed H 2 storage vessel model subjected to thermomechanical loadings to investigate effects of the helical layer fiber orientation and loading scenario on damage development, vessel integrity and burst pressure.

08 HYDROGEN↗
DOE OSTI2022

Computational Design of a Simple Flyer Plate Launcher

The response of materials to shock loading is important to understand for a variety of applications. When shock physics emerged during and after WWII, direct explosive loading or explosively driven plate impact was the primary tool for these studies. Subsequent decades have seen the widespread use of large caliber guns for plate impact studies, laser-shock facilities and pulsed power facilities. INL currently lacks a gun suitable for plate impact or explosives casting and machining facilities; however, it does possess explosives use and handling capabilities. An option for performing plate impact experiments was needed, therefore continuum scale models were utilized to explore a few simple donor-acceptor explosive plane wave lens designs, one of which could be hand packed with plastic explosives to launch flyer plates. 2D simulations were performed to study different geometries in an effort to minimize the difference in shock arrival across the central portion of a small copper flyer plate. A shock wave arrival time difference under 50 ns across 50 percent of the center of the flyer was achieved with a few designs. This work summarizes the computational models and results.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of the High-Pressure High-Temperature (HPHT) growth of single crystal diamond

A multi-scale, computational model is developed to describe the growth characteristics of single-crystal diamond in the High-Pressure, High-Temperature (HPHT) process. This model is the first to connect phase-change kinetics governing crystal growth to the continuum transport of carbon through the growth cell. Results show the importance of convective transport driven by buoyant flow in the metallic solvent, which increases the growth rate by nearly an order of magnitude over that obtained under diffusion alone. Parametric studies show how crystal growth may be kinetically-limited or transport-limited, depending on the value of the macroscopic kinetic coefficient. Estimating this kinetic coefficient from growth experiments yields a phase-change Damköhler number of unity, indicating a mixed regime where phase-change kinetics and transport are comparable and strongly coupled in this system. Mechanisms responsible for slowing growth as the crystal size increases are explained. Finally, supersaturation inhomogeneities along the facets of larger crystals are predicted, which may be relevant to solvent inclusion formation during growth.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Levels of Digitization, Digitalization, and Automation for Advanced Reactors

Much has been written about levels of automation (LOA), but comparatively little has been written about levels of digitization (LODi) and levels of digitalization (LODa). Digitization is a digital representation of analog information and is typical of migration from analog to digital control systems, digitalization involves enhancing the functionality of digital information, and automation changes control from humans to machines. Each of these technology implementations has its own scales, and each forms a viable type of functionality that should be considered not as a continuum toward automation but rather as separate categories of solutions that meet the needs of advanced reactors. In this paper we develop separate LODi, LODa, and LOA scales and demonstrate how conflation of these technologies, using the example of computer-based procedures, can lead to confusion in the design process. With the race to develop advanced reactors, the surest metric of success and safety is proper consideration of the right technology requirements for different control systems.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

A simple coin for a 2d entangled walk

We analyze the effect of a simple coin operator, built out of Bell pairs, in a 2d Discrete Quantum Random Walk (DQRW) problem. The specific form of the coin enables us to find analytical and closed form solutions to the recursion relations of the DQRW. The coin induces entanglement between the spin and position degrees of freedom, which oscillates with time and reaches a constant value asymptotically. We probe the entangling properties of the coin operator further, by two different measures. First, by integrating over the space of initial tensor product states, we determine the Entangling Power of the coin operator. Secondly, we compute the Generalized Relative Rényi Entropy between the corresponding density matrices for the entangled state and the initial pure unentangled state. Both the Entangling Power and Generalized Relative Rényi Entropy behaves similar to the entanglement with time. Finally, in the continuum limit, the specific coin operator reduces the 2d DQRW into two 1d massive fermions coupled to synthetic gauge fields, where both the mass term and the gauge fields are built out of the coin parameters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1970

Analysis of surface ablation of noncharring materials

Computer program solves combined problem of heat transfer and material response for the stagnation region of blunt bodies experiencing melting and vaporizing or subliming ablation. Program contains formulas for the transitional regime to bridge between the free-molecule and continuum regimes.

Matting, F. W.↗
NASA NTRS1971

Nonelastic nuclear reactions and accompanying gamma radiation

Several aspects of nonelastic nuclear reactions which proceed through the formation of a compound nucleus are dealt with. The full statistical model and the partial statistical model are described and computer programs based on these models are presented along with operating instructions and input and output for sample problems. A theoretical development of the expression for the reaction cross section for the hybrid case which involves a combination of the continuum aspects of the full statistical model with the discrete level aspects of the partial statistical model is presented. Cross sections for level excitation and gamma production by neutron inelastic scattering from the nuclei Al-27, Fe-56, Si-28, and Pb-208 are calculated and compared with avaliable experimental data.

Snow, R.↗
NASA NTRS1972

Spectra of type I supernovae. I - Generalities. II - Interpretation of the spectra of Sn 1960 in NGC 4496 from 8 to 20 days after the maximum of luminosity.

In the late stages of type I supernovae evolution (t greater than 50 days) the spectrum is without doubt an emission line spectrum, but in the early stages it is difficult to tell if an emission line spectrum is observed, or a continuum spectrum modified by the overlaying absorption formed in the outer layers. The mechanism of formation of the observed spectra, the physical conditions in the shell, different possible sources of heating of the shell, identification of the lines, and theoretical computations of the line intensities are discussed. An attempt is made to interpret the set of spectra assuming that a spectrum of widely broadened emission lines truncated by an overlaying absorption is seen. The absorption spectrum is first studied. An identification for the emission lines is proposed, and their intensity on the observed spectra is tentatively estimated. These identifications have to be justified by the theoretical computations of a model of the shell.

Gordon, C.↗
NASA NTRS1996

Computational Materials Research

Computational Materials aims to model and predict thermodynamic, mechanical, and transport properties of polymer matrix composites. This workshop, the second coordinated by NASA Langley, reports progress in measurements and modeling at a number of length scales: atomic, molecular, nano, and continuum. Assembled here are presentations on quantum calculations for force field development, molecular mechanics of interfaces, molecular weight effects on mechanical properties, molecular dynamics applied to poling of polymers for electrets, Monte Carlo simulation of aromatic thermoplastics, thermal pressure coefficients of liquids, ultrasonic elastic constants, group additivity predictions, bulk constitutive models, and viscoplasticity characterization.

Hinkley, Jeffrey A.↗
NASA NTRSDate not supplied

Application of OpenFOAM to Plume Impingement in Space Environments

After 30 years of continuous human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. Travelling beyond low earth orbit requires NASA to learn how humans can live in Deep Space environments – beyond the protection of Earth’s magnetosphere and at distances from Earth that prevent a quick return in case of trouble. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to be put in orbit around the moon to act as a home base for Lunar exploration for NASA astronauts. The Gateway Lunar outpost will be built incrementally, via modules which will arrive at separate times and dock to the existing structure. The incremental addition of Gateway modules, and the docking of visiting vehicles, is achieved via a sequence of firings from the approaching body’s onboard reaction control system (RCS) thrusters to achieve the required approach trajectory. The typical hypergolic chemical RCS thrusters work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking process. The exhaust gas from the RCS thrusters form plumes that expand into the vacuum of space and can impinge onto the outer surfaces of the Lunar Gateway, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion - all mechanisms that can damage the Lunar Gateway and must be minimized. Both permanent and visiting modules will have this RCS thruster exhaust impingement problem. This research aims to establish existing OpenFOAM solvers as a methodology for improving simulation techniques of rocket exhaust plume impingement in space environments. The flow structure of a plume in a space environment is complex; a plume that originates from a hypergolic chemical RCS thruster and expands into a vacuum will experience several regimes of rarefication. This range includes the continuum flow in the rocket nozzle through the fully rarefied free molecular flow further from the nozzle. The flow physics is different at these two extremes, and as such, the simulation approach for plumes is generally divided into a traditional computational fluid dynamics (CFD) simulation in and near the nozzle which is coupled to a subsequent direct simulation Monte Carlo (DSMC) simulation. At this time, the scope of this research is developing, verifying, and validating a method using existing solvers in the OpenFOAM framework for performing coupled CFD/DSMC calculations to determine the extent of plume impingement loading on generic space structures. This presentation will detail code-to-code comparisons between the hyStrath dsmcFoam+ solver, developed using OpenFOAM and available as open-source, and NASA’s in-house DSMC Analysis Code (DAC). Comparisons to several open-source publication findings using DAC [3,4] are presented, and advantages of using an OpenFOAM based solver are also discussed. The presentation concludes with a discussion of future work, and a plan for coupling the dsmcFoam+ solver with CFD simulations of chemical rocket engines for unified coupled plume simulation.

DSMC↗
NASA NTRS2001

CFD Simulation of Liquid Rocket Engine Injectors

Detailed design issues associated with liquid rocket engine injectors and combustion chamber operation require CFD methodology which simulates highly three-dimensional, turbulent, vaporizing, and combusting flows. The primary utility of such simulations involves predicting multi-dimensional effects caused by specific injector configurations. SECA, Inc. and Engineering Sciences, Inc. have been developing appropriate computational methodology for NASA/MSFC for the past decade. CFD tools and computers have improved dramatically during this time period; however, the physical submodels used in these analyses must still remain relatively simple in order to produce useful results. Simulations of clustered coaxial and impinger injector elements for hydrogen and hydrocarbon fuels, which account for real fluid properties, is the immediate goal of this research. The spray combustion codes are based on the FDNS CFD code' and are structured to represent homogeneous and heterogeneous spray combustion. The homogeneous spray model treats the flow as a continuum of multi-phase, multicomponent fluids which move without thermal or velocity lags between the phases. Two heterogeneous models were developed: (1) a volume-of-fluid (VOF) model which represents the liquid core of coaxial or impinger jets and their atomization and vaporization, and (2) a Blob model which represents the injected streams as a cloud of droplets the size of the injector orifice which subsequently exhibit particle interaction, vaporization, and combustion. All of these spray models are computationally intensive, but this is unavoidable to accurately account for the complex physics and combustion which is to be predicted, Work is currently in progress to parallelize these codes to improve their computational efficiency. These spray combustion codes were used to simulate the three test cases which are the subject of the 2nd International Workshop on-Rocket Combustion Modeling. Such test cases are considered by these investigators to be very valuable for code validation because combustion kinetics, turbulence models and atomization models based on low pressure experiments of hydrogen air combustion do not adequately verify analytical or CFD submodels which are necessary to simulate rocket engine combustion. We wish to emphasize that the simulations which we prepared for this meeting are meant to test the accuracy of the approximations used in our general purpose spray combustion models, rather than represent a definitive analysis of each of the experiments which were conducted. Our goal is to accurately predict local temperatures and mixture ratios in rocket engines; hence predicting individual experiments is used only for code validation. To replace the conventional JANNAF standard axisymmetric finite-rate (TDK) computer code 2 for performance prediction with CFD cases, such codes must posses two features. Firstly, they must be as easy to use and of comparable run times for conventional performance predictions. Secondly, they must provide more detailed predictions of the flowfields near the injector face. Specifically, they must accurately predict the convective mixing of injected liquid propellants in terms of the injector element configurations.

Farmer, Richard↗
DOE OSTI2021

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Resonant ultrasound spectroscopy measurement and modeling of additively manufactured octet truss lattice cubes

In this work, the elastic characteristics of additively manufactured Ti-5553 octet truss lattice cubes were examined using resonant ultrasound spectroscopy (RUS) techniques and correlated to X-ray computed tomography measurements of strut continuity. RUS models and measurements identified an outlier lattice possessing lower effective modulus despite larger effective density. Here, experimental measurements of density and effective Young's modulus were treated as input parameters in a homogenous anisotropic continuum model. The continuum model exhibited good agreement with experimental RUS measurements, thus confirming the potential for a simplified approximation of the octet truss lattice.

42 ENGINEERING↗
NASA NTRS1983

Development of a computational model for predicting solar wind flows past nonmagnetic terrestrial planets

A computational model for the determination of the detailed plasma and magnetic field properties of the global interaction of the solar wind with nonmagnetic terrestrial planetary obstacles is described. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super-Alfvenic solar wind flow past terrestrial ionospheres.

Stahara, S. S.↗
NASA NTRS1989

Flare model sensitivity of the Balmer spectrum

Careful studies of various chromospheric spectral signatures are very important in order to explore their possible sensitivity to the modifications of the thermodynamic quantities produced by the flare occurrence. Pioneer work of Canfield and co-workers have shown how the H alpha behavior is able to indicate different changes in the atmospheric parameters structure associated to the flare event. It was decided to study the behavior of the highest Balmer lines and of the Balmer continuum in different solar flare model atmospheres. These spectral features, originating in the deep photosphere in a quiet area, may have a sensitivity different from H alpha to the modification of a flare atmosphere. The details of the method used to compute the Stark profile of the higher Balmer line (n is greater than or equal to 6) and their merging were extensively given elsewhere (Donati-Falchi et al., 1985; Falchi et al., 1989). The models used were developed by Ricchiazzi in his thesis (1982) evaluating the chromospheric response to both the nonthermal electron flux, for energy greater than 20 kev, (F sub 20) and to the thermal conduction, (F sub c). The effect of the coronal pressure values (P sub O) at the apex of the flare loop is also included.

Falchi, A.↗