Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 199 records · Page 11

High-Resolution LiDAR Observations for Coupled Effects of Dry-Air Entrainment and Haze-Cloud Interactions on Cloud Vertical Structure

This study demonstrates the high-resolution profiling of cloud microphysics in a laboratory chamber using Time-Correlated Single Photon Counting (TCSPC) LiDAR. We present a novel retrieval method to derive vertical extinction (𝜎) profiles, constrained by in situ measurements, to diagnose responses to dry-air entrainment. In clean clouds, the LiDAR signals and retrieved 𝜎 remain relatively uniform, with entrainment effects confined to the upper layer. In contrast, polluted clouds exhibit strong vertical variability and a transition to a water-vapor-limited state. Entrainment significantly enhances the haze number concentration (𝑁ℎ), particularly near the bottom, creating highly height-dependent extinction profiles for polluted clouds. Our results highlight the capability of high-resolution LiDAR in capturing fine-scale vertical inhomogeneities. This approach provides a robust framework for quantifying how aerosol loading modulates entrainment sensitivity, offering new insights into the transition between buffered and water-vapor-limited regimes.

54 ENVIRONMENTAL SCIENCES↗

Development of a point kinetics subroutine for Molten Salt reactors in RELAP5-3D

RELAP5-3D has historically focused on accurate modeling of design-basis and beyond-design-basis accidents in solid fuel light-water reactors. The assumption used in solid fuel reactors that delayed neutron precursors remain at the location in which they were produced is inherently incorrect for a molten salt reactor with flowing fuel. Therefore, new methods for modeling MSR kinetics were required. Subroutine used in RELAP5-3D were updated to analyze point kinetics in a molten salt reactor during both steady state and transient analyses. The phenomena related to MSR kinetics that were added to RELAP5-3D include arrays to store previous values for independent and dependent variables, arrays to add terms to existing kinetics equations, and a new reactivity bias term. A custom MATLAB assessment code was used to establish a suite of verification tests. This capability is undergoing final documentation and testing for public release with a future version of RELAP5-3D.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

On the use of internal state variables in thermoviscoplastic constitutive equations

The general theory of internal state variables are reviewed to apply it to inelastic metals in use in high temperature environments. In this process, certain constraints and clarifications will be made regarding internal state variables. It is shown that the Helmholtz free energy can be utilized to construct constitutive equations which are appropriate for metallic superalloys. Internal state variables are shown to represent locally averaged measures of dislocation arrangement, dislocation density, and intergranular fracture. The internal state variable model is demonstrated to be a suitable framework for comparison of several currently proposed models for metals and can therefore be used to exhibit history dependence, nonlinearity, and rate as well as temperature sensitivity.

Allen, D. H.↗

Bifurcation properties of a stratospheric vacillation model

Nonlinear properties of the stratospheric vacillation model of Holton and Mass (1976) are studied numerically using bifurcation theory. Severe truncation and vertical differencing are used to obtain 81 nonlinear ordinary differential equations with time-dependent variables. Three branches of the steady solutions are determined using Powell's hybrid method and the pseudoarclength continuation method, and a multiplicity of stable steady-state solutions with different vertical structures are found to exist in some range of the bifurcation parameter. Periodic solutions are found which branch off from a steady solution by a Hopf bifurcation. It is suggested that the interannual variability of the stratospheric circulation in the middle and high latitudes during winter may be explained by the multiplicity of steady and periodic stable solutions.

Yoden, Shigeo↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 1: Theory and numerical solution procedures

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 1 of a series of three reference publications that describe LENS, provide a detailed guide to its usage, and present many example problems. Part 1 derives the governing equations and describes the numerical solution procedures for the types of problems that can be solved. The accuracy and efficiency of LSENS are examined by means of various test problems, and comparisons with other methods and codes are presented. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 3: Illustrative test problems

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 3 of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part 3 explains the kinetics and kinetics-plus-sensitivity analysis problems supplied with LSENS and presents sample results. These problems illustrate the various capabilities of, and reaction models that can be solved by, the code and may provide a convenient starting point for the user to construct the problem data file required to execute LSENS. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Bittker, David A.↗

LSENS, a general chemical kinetics and sensitivity analysis code for homogeneous gas-phase reactions. 2: Code description and usage

LSENS, the Lewis General Chemical Kinetics Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 2 of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part 2 describes the code, how to modify it, and its usage, including preparation of the problem data file required to execute LSENS. Code usage is illustrated by several example problems, which further explain preparation of the problem data file and show how to obtain desired accuracy in the computed results. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions. Part 1 (NASA RP-1328) derives the governing equations describes the numerical solution procedures for the types of problems that can be solved by lSENS. Part 3 (NASA RP-1330) explains the kinetics and kinetics-plus-sensitivity-analysis problems supplied with LSENS and presents sample results.

Radhakrishnan, Krishnan↗

LSENS, A General Chemical Kinetics and Sensitivity Analysis Code for Homogeneous Gas-Phase Reactions: Code Description and Usage - Part 2

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part II of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part II describes the code, how to modify it, and its usage, including preparation of the problem data file required to execute LSENS. Code usage is illustrated by several example problems, which further explain preparation of the problem data file and show how to obtain desired accuracy in the computed results. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions. Part I (NASA RP-1328) derives the governing equations and describes the numerical solution procedures for the types of problems that can be solved by LSENS. Part III (NASA RP-1330) explains the kinetics and kinetics-plus-sensitivity-analysis problems supplied with LSENS and presents sample results.

Radhakrishnan, Krishnan↗

The value of fusion energy to a decarbonized United States electric grid

Fusion could be a part of future decarbonized electricity systems, but it will need to compete with other technologies. In particular, pulsed tokamak plants have a unique operational mode, and evaluating which characteristics make them economically competitive can help select between design pathways. Using a capacity expansion and operations model, here we determined cost thresholds for pulsed tokamaks to reach a range of penetration levels in a future decarbonized US Eastern Interconnection. The required capital cost to reach a fusion capacity of 100 GW varied from $2,700 to $7,500 kW -1 , and the equilibrium penetration increases rapidly with decreasing cost. The value per unit power capacity depends on the variable operational cost and on the cost of its competition, particularly fission, much more than on the pulse cycle parameters. These findings can therefore provide initial cost targets for fusion more generally in the United States.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurement report: Role of organic coating and chemical composition on ice nucleation potential of atmospheric particles in European Arctic

Understanding the ice nucleation (IN) potential of Arctic aerosols is critical for predicting their influence on cloud formation and water cycles in this vulnerable region. This study investigates the role of particle composition, organic coatings, and aerosol sources in modulating ice nucleating particle (INPs) abundance across five aerosol samples collected at the Gruvebadet Observatory Station in Ny-Ålesund, Svalbard. The IN potential of Arctic aerosol particles was studied by investigating chemical, morphological, and INP abundance measurements. Single-particle analyses revealed distinct differences in mixing state, organic volume fraction (OVF), and organic coating morphology across samples. OVF distributions were linked to particle origin, with marine-influenced Na-rich particles often exhibiting thin organic coatings, while long-range transported particles showed thicker organic coatings. Biogenic contributions, though variable, were linked to heat-sensitive INPs, suggesting a role for labile biological macromolecules under certain meteorological conditions. Spearman rank correlation analysis between particle composition and immersion-mode INP concentrations at two freezing temperatures indicated that organic-rich and Na-rich particles were positively associated with enhanced INP abundance. However, discrepancies in INP abundance were observed for particles with thicker organic coatings, where the morphological configuration of the organic material may play a role. The results highlight that Arctic INP variability is governed not only by chemical composition but also by the morphological configuration of organic material, which can either enhance or inhibit ice nucleation depending on its abundance, distribution, thickness, and mixing state. These findings underscore the combined influence of source regions, atmospheric processing, and organic–inorganic interactions in shaping Arctic aerosol freezing behavior.

Lata, Nurun Nahar [Pacific Northwest National Labo↗

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasizes the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M qi) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in turns of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, Derck↗

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasises the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters, but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M (dot) q(sub i)) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in terms of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, D.↗

Space-Time Finite Element Tensor Network Approach for the Time-Dependent Convection–Diffusion–Reaction Equation with Variable Coefficients

In this paper, we present a new space-time Galerkin-like method, where we treat the discretization of spatial and temporal domains simultaneously. This method utilizes a mixed formulation of the tensor-train (TT) and quantized tensor-train (QTT) (please see Section Tensor-Train Decomposition), designed for the finite element discretization (Q1-FEM) of the time-dependent convection–diffusion–reaction (CDR) equation. We reformulate the assembly process of the finite element discretized CDR to enhance its compatibility with tensor operations and introduce a low-rank tensor structure for the finite element operators. Recognizing the banded structure inherent in the finite element framework’s discrete operators, we further exploit the QTT format of the CDR to achieve greater speed and compression. Additionally, we present a comprehensive approach for integrating variable coefficients of CDR into the global discrete operators within the TT/QTT framework. The effectiveness of the proposed method, in terms of memory efficiency and computational complexity, is demonstrated through a series of numerical experiments, including a semi-linear example.

convection–diffusion–reaction equation↗

Continuous conditional generative adversarial networks for data-driven solutions of poroelasticity with heterogeneous material properties

Machine learning-based data-driven modeling can allow computationally efficient time-dependent solutions of PDEs, such as those that describe subsurface multiphysical problems. In this work, our previous approach (Kadeethum et al., 2021d) of conditional generative adversarial networks (cGAN) developed for the solution of steady-state problems involving highly heterogeneous material properties is extended to time-dependent problems by adopting the concept of continuous cGAN (CcGAN). The CcGAN that can condition continuous variables is developed to incorporate the time domain through either element-wise addition or conditional batch normalization. Moreover, this framework can handle training data that contain different timestamps and then predict timestamps that do not exist in the training data. As a numerical example, the transient response of the coupled poroelastic process is studied in two different permeability fields: Zinn & Harvey transformation and a bimodal transformation. The proposed CcGAN uses heterogeneous permeability fields as input parameters while pressure and displacement fields over time are model output. Our results show that the model provides sufficient accuracy with computational speed-up. This robust framework will enable us to perform real-time reservoir management and robust uncertainty quantification in poroelastic problems.

97 MATHEMATICS AND COMPUTING↗

Theoretical study of the reaction mechanism and kinetics of the phenyl + propargyl association

Potential energy surface for the phenyl + propargyl radical recombination reaction has been studied at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G** level of theory for the closed-shell singlet species and at the triplet–singlet gap CASPT2/cc-pVTZ-CCSD(T)-F12/cc-pVTZ-f12//CASSCF/cc-pVTZ level of theory for the diradical species. High-pressure limit rate constants for the barrierless channels were evaluated with variable reaction coordinate transition state theory (VRC-TST). Rice–Ramsperger–Kassel–Marcus Master Equation (RRKM-ME) calculations have been performed to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the entrance channels of the radical association reaction produce 3-phenyl-1-propyne and phenylallene which can further dissociate/isomerize into a variety of unimolecular and bimolecular products. Theoretical evidence is presented that, at combustion relevant conditions, the phenyl + propargyl recombination provides a feasible mechanism for the addition of a second five-member ring to the first six-member aromatic ring producing the prototype two-ring species indene and indenyl. Rate expressions for all important reaction channels in a broad range of temperatures and pressures have been generated for kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Leading-edge receptivity to a vortical freestream disturbance: A numerical analysis

The receptivity to freestream vorticity of the boundary layer over a flat plate with an elliptic leading edge is investigated numerically. The flow is simulated by solving the incompressible Navier-Stokes system in general curvilinear coordinates with the vorticity and stream function as dependent variables. A finite-difference scheme which is second-order accurate in both space and time is used. As a first step, the steady basic-state solution is computed. Then a small amplitude vortical disturbance is introduced at the upstream boundary and the governing equations are solved time-accurately to evaluate the spatial and temporal growth of the perturbations leading to instability waves (Tollmien-Schlichting waves) inside the boundary layer. Preliminary results for a symmetric, 2-D disturbance reveal the presence of Tollmien-Schlichting waves aft of the flat-plate/ellipse juncture.

Buter, Thomas A.↗

Experimental Investigation of Organic Synthesis in Hydrothermal Environments

Seafloor hydrothermal systems may be the most likely locations on the early Earth for the emergence of life. Because of the disequilibrium inherent in such dynamic, mixing environments, abundant chemical energy would have been available for formation of the building blocks of life. In addition, theoretical studies suggest that organic compounds in these conditions would reach metastable states, due to kinetic barriers to the formation of stable equilibrium products (CO2 and methane). The speciation of organic carbon in metastable states is highly dependent on the oxidation state, pH, temperature, pressure and bulk composition of the system. The goal of our research is to investigate the effects of a number external variables on the formation, transformation, and stability of organic compounds at hydrothermal conditions. We have begun experimental work to attempt to control the oxidation state of simulated hydrothermal systems by using buffers composed of mineral powders and gas mixtures. We are also beginning to test the stability of organic compounds under these conditions. The experiments are being performed using the hydrothermal bomb apparatus at the U.S. Geological Survey in Menlo Park, CA and the supercritical water oxidizer (SCWO) at NASA Ames Research Center in Moffet Field, CA. The amino acids decomposed rapidly. Even after the approximately 15 minutes between addition of the amino acids and the first sampling, no amino acids were detected in the PPM system by GC- MS, while in the FeFeO system the amino acids were present at a level of less than 50% of original. Carboxylic acids, ammonia, and CO2 were the main products, along with some unidentified compounds. The ratios of carboxylic acids and concentrations of other products seem to have remained stable during the experiments, consistent with observations of other metastable systems and theoretical predictions.

Shock, Everett L.↗

Experimental Investigation of Organic Synthesis in Hydrothermal Environments

The results of the investigation were presented at a Astrobiology Institute General Meeting. Seafloor hydrothermal systems may be the most likely locations on the early Earth for the emergence of life. Because of the disequilibrium inherent in such dynamic, mixing environments, abundant chemical energy would have been available for formation of the building blocks of life. In addition, theoretical studies suggest that organic compounds in these conditions would reach metastable states, due to kinetic barriers to the formation of stable equilibrium products (CO2 and methane). The speciation of organic carbon in metastable states is highly dependent on the oxidation state, pH, temperature, pressure and bulk composition of the system. The goal of our research is to investigate the effects of a number external variables on the formation, transformation, and stability of organic compounds at hydrothermal conditions. We have begun experimental work to attempt to control the oxidation state of simulated hydrothermal systems by using buffers composed of mineral powders and gas mixtures. We are also beginning to test the stability of organic compounds under these conditions.

Shock, Everett L.↗