Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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

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

At least 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Tuning Stoichiometry to Promote Formation of Binary Colloidal Superlattices

The self-assembly of binary nanoparticle superlattices from colloidal mixtures is a promising method for the fabrication of complex colloidal cocrystal structures. However, binary mixtures often form amorphous or metastable phases instead of the thermodynamically stable phase. Here we show that in binary mixtures of differently sized spherical particles, an excess of the smaller component can promote—and, in some cases, may be necessary for—the self-assembly of a binary cocrystal. Using computer simulations, we identify two mechanisms responsible for this phenomenon. First, excess small particles act like plasticizers and enable systems to reach a greater supersaturation before kinetic arrest occurs. Second, they can disfavor competing structures that may interfere with the growth of the target structure. We find the phase behavior of simulated mixtures of nearly hard spheres closely matches published experimental results. We demonstrate the generality of our findings for mixtures of particles of arbitrary shape by presenting a binary mixture of hard shapes that only self-assembles with an excess of the smaller component.

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

T -matrix analysis of static Wilson line correlators from lattice QCD at finite temperature

Here, we utilize a previously constructed thermodynamic T -matrix approach to the quark-gluon plasma (QGP) to derive constaints on its input by using results on Wilson line correlators (WLCs) of a static quark-antiquark pair from 2 + 1-flavor lattice-QCD (lQCD) computations with realistic pion mass. The self-consistent T -matrix results, which include previous constraints from the lQCD equation of state in the light-parton sector, can describe the lQCD data for WLCs fairly well once refinements of its driving kernel are applied. In particular, the input potential requires less screening than what has been inferred from previous T-matrix analyses. Pertinent predictions for the spectral and transport properties of the QGP are discussed, including the spatial diffusion coefficient for heavy quarks; the latter turns out to have a rather weak temperature dependence, in approximate agreement with recent unquenched lQCD results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1976

Thermodynamics of a solar mixture of molecular hydrogen and helium at high pressure

The thermodynamic properties of a model molecular hydrogen and helium mixture are calculated in the strongly interacting region of 0.005 to 0.3 per cu cm for a range of temperatures that are of interest for the envelopes of the Jovian planets. Computed adiabats fit the gravity data and boundary conditions from model atmospheres of Jupiter.

Slattery, W. L.↗
NASA NTRS1994

Interactions between solidification and compositional convection in mushy layers

Mushy layers are ubiquitous during the solidification of alloys. They are regions of mixed phase wherein solid crystals are bathed in the melt from which they grew. The matrix of crystals forms a porous medium through which the melt can flow, driven either by external forces or by its own buoyancy in a gravitational field. Buoyancy-driven convection of the melt depends both on temperature gradients, which are necessary for solidification, and on compositional gradients, which are generated as certain components of the alloy are preferentially incorporated in the solid phase and the remaining components are expelled into the melt. In fully liquid regions, the combined action of temperature and concentration on the density of the liquid can cause various forms of double-diffusive convection. However, in the interior of mushy regions the temperature and concentration are thermodynamically coupled so only single-diffusive convection can occur. Typically, the effect of composition on the buoyancy of the melt is much greater than the effect of temperature, and thus convection in mushy layers in driven primarily by the computational gradients within them. The rising interstitial liquid is relatively dilute, having come from colder regions of the mushy layer, where the liquidus concentration is lower, and can dissolve the crystal matrix through which it flows. This is the fundamental process by which chimneys are formed. It is a nonlinear process that requires the convective velocities to be sufficiently large, so fully fledged chimneys (narrow channels) might be avoided by means that weaken the flow. Better still would be to prevent convection altogether, since even weak convection will cause lateral, compositional inhomogeneities in castings. This report outlines three studies that examine the onset of convection within mushy layers.

Worster, M. Grae↗
NASA NTRS2021

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA’s exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA’s Entry Systems Modeling (ESM) Project, which serves as the Agency’s only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specificmodels for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures. Results for both configurations have demonstrated the ability to land within the 50-meter precision requirement demanded by Mars human exploration missions, while also reducing propellant requirements by enabling more efficient control through entry and descent. Ongoing research in GN&C is developing mechanical specifications for the systems and establishing engineering feasibility.

aerothermodynamics↗
DOE OSTI2025

Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems

The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2019

Recent Advancements in Modeling and Simulation of Entry Systems at NASA

This paper describes recent development of modeling and simulation technologies for entry systems in support of NASA's exploration missions. Mission-tailored research and development in modeling of entry systems occurs across the Agency (e.g., within the Orion and Mars 2020 Programs), however the aim of this paper is to discuss the broad, cross-mission research conducted by NASA's Entry Systems Modeling (ESM) Project, which serves as the Agency's only concerted effort toward advancing entry systems across a range of technical disciplines. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Predictive material modeling, (2) Shock layer kinetics and radiation, (3) Computational and experimental aerosciences, and (4) Guidance, navigation, and control. Investments in thermal protection material modeling are geared toward high-fidelity, predictive models capable of handling complex structures, with an eye toward optimizing design performance and quantifying thermal protection system reliability. New computational tools have been developed to characterize material properties and behavior at the microstructural level, and experimental techniques (molecular beam scattering, micro-computed tomography, among others) have been developed to measure material kinetics, morphology, and other parameters needed to inform and validate detailed simulations. Advancements have also been made in macrostructural simulation capability to enable 3-D system-scale calculations of material response with complex topological features, including differential recession of tile gaps. Research and development in the area of shock layer kinetics has focused on air and CO2-based atmospheres. Capacity and capability of the NASA Ames Electric Arc Shock Tube (EAST) have been expanded in recent years and analysis of resulting data has led to several improvements in kinetic models, while simultaneously reducing uncertainties associated with radiative heat transfer predictions. First-principles calculations of fundamental kinetic, thermodynamic, and transport data, along with state-specific models for non-equilibrium flow regimes, have also yielded new insights and have the potential to vastly improve model fidelity. Aerosciences is a very broad area of interest in entry systems, yet a number of important challenges are being addressed: Coupled fluid-structure simulations of parachute inflation and dynamics; Experimental and computational studies of vehicle dynamics; Multi-phase flow with dust particles to simulate entry environments at Mars during dust storms; Studies of roughness-induced heating augmentation relevant to tiled and woven thermal protection systems; and Advanced numerical methods to optimize computational analyses for desired accuracy versus cost. Guidance and control in the context of entry systems has focused on development of methods for multi-axis control (i.e. pitch and yaw, rather than bank angle alone) of spacecraft during entry and descent. With precision landing requirements driven by Mars human exploration goals, recent efforts have yielded 6-DOF models of multi-axis control with propulsive descent of both inflatable and rigid ellipsled-like architectures.

Barnhardt, Michael D.↗
DOE OSTI2022

SAXS-guided unbiased coarse-grained Monte Carlo simulation for identification of self-assembly nanostructures and dimensions

Recent studies have shown that solvated amphiphiles can form nanostructured self-assemblies called dynamic binary complexes (DBCs) in the presence of ions. Since the nanostructures of DBCs are directly related to their viscoelastic properties, it is important to understand how the nanostructures change under different solution conditions. However, it is challenging to obtain a three-dimensional molecular description of these nanostructures by utilizing conventional experimental characterization techniques or thermodynamic models. To this end, we combined the structural data from small angle X-ray scattering (SAXS) experiments and thermodynamic knowledge from coarse-grained Monte Carlo (CGMC) simulations to identify the detailed three-dimensional nanostructure of DBCs. Specifically, unbiased CGMC simulations are performed with SAXS-guided initial conditions, which aids us to sample accurate nanostructures in a computationally efficient fashion. As a result, an elliptical bilayer nanostructure is obtained as the most probable nanostructure of DBCs whose dimensions are validated by scanning electron microscope (SEM) images. Then, utilizing the obtained molecular model of DBCs, here we could also explain the pH tunability of the system. Overall, our results from SAXS-guided unbiased CGMC simulations highlight that using potential energy combined with SAXS data, we can distinguish otherwise degenerate nanostructures resulting from the inherent ambiguity of SAXS patterns.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Transferrin receptor targeting by de novo sheet extension

The de novo design of polar protein–protein interactions is challenging because of the thermodynamic cost of stripping water away from the polar groups. Here, we describe a general approach for designing proteins which complement exposed polar backbone groups at the edge of beta sheets with geometrically matched beta strands. We used this approach to computationally design small proteins that bind to an exposed beta sheet on the human transferrin receptor (hTfR), which shuttles interacting proteins across the blood–brain barrier (BBB), opening up avenues for drug delivery into the brain. We describe a design which binds hTfR with a 20 nM K d , is hyperstable, and crosses an in vitro microfluidic organ-on-a-chip model of the human BBB. Our design approach provides a general strategy for creating binders to protein targets with exposed surface beta edge strands.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Bottomonium properties in the quark-gluon plasma from a lattice-QCD informed 𝑇-matrix approach

Recent computations of bottomonium correlation functions with extended sources in lattice-discretized quantum chromodynamics (lQCD) provide new insights into heavy-quark dynamics at distance scales which are of the order of the inverse temperature. We analyze these results employing the thermodynamic 𝑇-matrix framework, in a continued effort to interpret lQCD data for quarkonium correlation functions in a nonperturbative and self-consistently solved quantum-many-body approach to a strongly coupled quark-gluon plasma (QGP). Its key inputs are the in-medium driving kernel (potential) of the scattering equation and an interference function which implements 3-body effects in the quarkonium coupling to the thermal medium. A simultaneous description of lQCD results for the bottomonium correlators with extended operators and the previously analyzed Wilson line correlators only requires minor refinements of the potential but calls for stronger interference effects at larger separation of the bottom quark and antiquark. We then analyze the poles of the self-consistent 𝑇 matrices on the real axis to assess the survival of the various bound states. Here, we estimate the pertinent temperatures where the poles disappear for the various bottomonium states and discuss the relation to the corresponding peaks in the bottomonium spectral functions. We also recalculate the spatial diffusion coefficient of the QGP and find it to be similar to that in our previous study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Hydrodynamic Coupling to a Homogenized Radiation Transport Method based on Young Measures

Resolving radiation transport fields subject to opacity profiles with strong, oscillatory line structure while potentially falling under intermediate optical depth conditions presents a numerical challenge in radiation transport modeling. The Young measure-based homogenization technique formulated by Haut et al. (2017) was investigated as a candidate method for resolving radiation fields under such conditions more accurately. The method was compared against frequently-utilized mean opacity methods as the Rosseland and Planck formulations. In this work, all methods were tested through radiation slab calculations separately comprised of aluminum, copper, and krypton, each for different thermodynamic conditions. Following these offline radiation slab calculations, demonstrations shifted towards the SCEPTRE radiation transport code and, subsequently, the multiphysics ALEGRA code for approximately-coupled radiation-material simulations. Throughout all the simulations shown in this study, for a fixed computational cost, the homogenized method was observed to be more accurate than any of the solutions determined through traditional mean opacity approaches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1974

Computer program for calculation of real gas turbulent boundary layers with variable edge entropy

A user's manual for a computer program which calculates real gas turbulent boundary layers with variable edge entropy on a blunt cone or flat plate at zero angle of attack is presented. An integral method is used. The method includes the effect of real gas in thermodynamic equilibrium and variable edge entropy. A modified Crocco enthalpy velocity relationship is used for the enthalpy profiles and an empirical correlation of the N-power law profile is used for the velocity profile. The skin-friction-coefficient expressions of Spalding and Chi and Van Driest are used in the solution of the momentum equation and in the heat-transfer predictions that use several modified forms of Reynolds analogy.

Boney, L. R.↗
NASA NTRS1979

Modeling of radiative heating of base region of Jovian entry probe

The mechanism of radiative heating of the afterbody region of Jovian entry probe is analyzed. A theoretical model is derived to determine the average thermodynamic properties in the expanding region, recirculating region, recompression region, and neck region through application of one-dimensional conservation equations. Flow parameters are obtained from the shadowgraphs of a free-flight test. Radiative transfer is calculated using spectrally detailed computer codes accounting for nonequilibrium. The results show that the most severe heating occurs immediately behind the frustum, and that the recompression and neck regions are the major sources of radiation that heats the base stagnation point. The radiation flux to the base point is slightly stronger with ablation than without, its value being 0.11(43Pb/Ps)squared times that to the front stagnation point, where Pb/Ps is the ratio of base-to-front stagnation point pressures and its value is in the range 0.023-0.066. The time-integrated heat load to the base point is 18(43Pb/Ps)squared kJ/sq cm. Existing experimental data are shown to agree with the theoretical prediction.

Park, C.↗
NASA NTRS1994

Experimental and numerical study of swept ramp injection into a supersonic flowfield

Time-averaged measurements of pressure, temperature, velocity, and injectant mole fraction are presented using the planar laser-induced iodine fluorescence technique in the complex three-dimensional compressible flowfield around a swept ramp fuel injector. Within the range of thermodynamic conditions present in the test case studied, the technique's accuracy is estimated to be 4% for pressure, temperature, and velocity and 3% for injectant mole fraction. Comparisons with numerical simulations using the SPARK three-dimensional Navier-Stokes computer code with an algebraic turbulence model are made at the centerplane of the flowfield as well as on three crossflow planes downstream of the injector. Calculations and measurements are in good agreement throughout the flowfield, with deviations on the order of 5%; however, in specific regions, such as in the base of the ramp, deviations are larger. A weak asymmetry in the incoming flowfield appears to be amplified by boundary-layer separation occurring when the ramp-generated shock reflects off the tunnel walls. Ramp-generated vortices are weaker in the calculated results due to the effects of numerical viscosity in the vortex cores. This leads to less turning and mixing of the jet plume than observed in the experiments. The rate of decay of the maximum injectant mole fraction with streamwise distance is greater for the present ramp injection scheme than for previously measured transverse injection schemes. In this recirculation region at the base of the injector, laminar calculations show better agreement with the measurements than turbulent calculations.

Donohue, James M.↗
NASA NTRS2003

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Cloud microphysics are inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e.,pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size-distribution functions.A spectral-bin microphysical model is very expensive from a from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region using identical thermodynamic conditions but with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. Besides the initial differences in aerosol concentration, preliminary results indicate that the low CCN concentration case produces rainfall at the surface sooner than the high CCN case but has less cloud water mass aloft. Because the spectral-bin model explicitly calculates and allows for the examination of both the mass and number concentration of species in each size categor, a detailed analysis of the instantaneous size spectrum can be obtained for the two cases. It is shown that since the low CCN case produces fewer droplets, larger sized develop due to the greater condensational and collectional growth, leading to a broader size spectrum in comparison to the high CCN case.

Tao, Wei-Kuo↗
NASA NTRSDate not supplied

Axisymmetric Two-Dimensional Modeling of No Vent Filling of a Cryogenic Tank using Generalized Fluid System Simulation Program

This paper presents preliminary analysis using a finite volume based computational method for network flow analysis for the simulation of chill and fill of a cryogenic tank. An axisymmetric two-dimensional model of a spherical tank was developed using a nodal approach. No Vent Filling (NVF) was accomplished using a Thermodynamic Vent System assisted injector where the injector was cooled by the cold gas generated by Joule-Thompson expansion of the liquid propellant. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around spray droplets, and c) condensation of vapor at the interface of the cooled injector and ullage. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data and one-dimensional (1D) model predictions. It was demonstrated that the axisymmetric 2D model provided a more realistic prediction than the 1D model.

Cryogenics↗
DOE OSTI2020

Comparison of fixed charge and polarizable models for predicting the structural, thermodynamic, and transport properties of molten alkali chlorides

Results from extensive molecular dynamics simulations of molten LiCl, NaCl, KCl, and RbCl over a wide range of temperatures are reported. Comparison is made between the “Polarizable Ion Model” (PIM) and the non-polarizable “Rigid Ion Model” (RIM). Densities, self-diffusivities, shear viscosities, ionic conductivities, and thermal conductivities are computed and compared with experimental data. In addition, radial distribution functions are computed from ab initio molecular dynamics simulations and compared with the two sets of classical simulations as well as experimental data. Here, the two classical models perform reasonably well at capturing structural and dynamic properties of the four molten alkali chlorides, both qualitatively and often quantitatively. With the singular exception of liquid density, for which the PIM is more accurate than the RIM, there are few clear trends to suggest that one model is more accurate than the other for the four alkali halide systems studied here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗