SEARCH · Search NASA

Results for “nonlinear reaction”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS2020

Nuclear Fusion Reactions in Deuterated Metals

Nuclear fusion reactions of D-D are examined in an environment comprised of high density cold fuel embedded in metal lattices in which a small fuel portion is activated by hot neutrons. Such an environment provides for enhanced screening of the Coulomb barrier due to conduction and shell electrons of the metal lattice, or by plasma induced by ionizing radiation (γ quanta). We show that neutrons are far more efficient than energetic charged particles, such as light particles (e−, e+) or heavy particles (p, d, α) in transferring kinetic energy to fuel nuclei (D) to initiate fusion processes. It is well-known that screening increases the probability of tunneling through the Coulomb barrier. Electron screening also significantly increases the probability of large- versus small-angle Coulomb scattering of the reacting nuclei to enable subsequent nuclear reactions via tunneling. This probability is incorporated into the astrophysical factor S(E). Aspects of screening effects to enable calculation of nuclear reaction rates are also evaluated, including Coulomb scattering and localized heating of the cold fuel, primary D-D reactions, and subsequent reactions with both the fuel and the lattice nuclei. The effect of screening for enhancement of the total nuclear reaction rate is a function of multiple parameters including fuel temperature and the relative scattering probability between the fuel and lattice metal nuclei. Screening also significantly increases the probability of interaction between hot fuel and lattice nuclei increasing the likelilhood of Oppenheimer-Phillips processes opening a potential route to reaction multiplication. We demonstrate that the screened Coulomb potential of the target ion is determined by the nonlinear Vlasov potential and not by the Debye potential. In general, the effect of screening becomes important at low kinetic energy of the projectile. We examine the range of applicability of both the analytical and asymptotic expressions for the well-known electron screening lattice potential energy Ue, which is valid only for E >> Ue (E is the energy in the center of mass reference frame). We demonstrate that for E ≤ Ue, a direct calculation of Gamow factor for screened Coulomb potential is required to avoid unreasonably high values of the enhancement factor f (E) by the analytical—and more so by the asymptotic—formulas.

Vladimir Pines↗

NASA NTRS1980

The collective excitation of g-modes in the sun

Oscillations of the solar interior (linear g-modes) may be strongly driven by the collective influence of all the modes upon the nuclear reactions in the core. This heretofore neglected effect could couple the modes, reduce the effective amplitudes near the center, and spatially concentrate most of the oscillation energy into just a portion of the radiative interior. If operating at sufficient strength, this can reverse the conventional conclusion, drawn from single mode calculations, that almost all solar g-modes are damped. Furthermore, it would put the theory into rough harmony with three otherwise troubling observations: (1) the 'low' neutrino flux measured by Davis (1978), (2) the high correspondence found by Wolff (1976) between recurrence periods in solar activity and the rotational beat periods of g-modes, and (3) the fluctuations in the sun's diameter which imply g-mode activity at high angular harmonics (Hill and Caudell 1979). A nonlinear expression is derived for the local rate of work done on an array of oscillation modes by the nuclear reactions. Three additional tests of the model are suggested.

Wolff, C. L.↗

DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2015

Design and Stability of an On-Orbit Attitude Control System Using Reaction Control Thrusters

Principles for the design and stability of a spacecraft on-orbit attitude control system employing on-off Reaction Control System (RCS) thrusters is presented. Both the vehicle dynamics and the control system actuators are inherently nonlinear, hence traditional linear control system design approaches are not directly applicable. This paper has three main aspects: It summarizes key RCS control System design principles from the Space Shuttle and Space Station programs, it demonstrates a new approach to develop a linear model of a phase plane control system using describing functions, and applies each of these to the initial development of the NASA's next generation of upper stage vehicles. Topics addressed include thruster hardware specifications, phase plane design and stability, jet selection approaches, filter design metrics, and automaneuver logic.

Hall, Robert A.↗

DOE OSTI2026

Phenomena-based graph representations and applications to chemical process simulation

Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.

Convergence↗

DOE OSTI2023

Microphysically modified magnetosonic modes in collisionless, high-β plasmas

With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-β collisionless plasmas, with particular attention to their excitation of and reaction to kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpendicular pressure over parallel pressure in regions where the plasma β is increased. For mode amplitudes |δB/B 0 |≳0.3, this excess excites the mirror instability. Particle scattering off these micro-scale mirrors frustrates the nonlinear saturation of transit-time damping, ensuring that large-amplitude NP modes continue their decay to small amplitudes. At asymptotically large wavelengths, we predict that the mirror-induced scattering will be large enough to interrupt transit-time damping entirely, isotropizing the pressure perturbations and morphing the collisionless NP mode into the magnetohydrodynamic (MHD) entropy mode. In fast waves, a fluctuating pressure anisotropy drives both mirror and firehose instabilities when the wave amplitude satisfies |δB/B 0 |≳2β -1 . The induced particle scattering leads to delayed shock formation and MHD-like wave dynamics. Taken alongside prior work on self-interrupting Alfvén waves and self-sustaining ion-acoustic waves, our results establish a foundation for new theories of electromagnetic turbulence in low-collisionality, high-β plasmas such as the intracluster medium, radiatively inefficient accretion flows and the near-Earth solar wind.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Inversion-based correction of Double-Torsion (DT) subcritical crack growth tests for crack profile geometry

Because of its simplicity and the ability to produce a stable, slow-propagating crack, the Double-Torsion (DT) method has been used widely for investigating the critical and subcritical propagation of a slow-propagating tensile (mode-I) crack. However, to determine the complex relationship between the crack velocity $\mathcal{v_c}$ vs. the strain energy release rate $\mathscr{G}$ (or the stress intensity factor K) from laboratory measurements, several corrections must be made to account for the impact of sample and crack geometry. Particularly, DT test typically produces a crack with a curved edge profile instead of a straight line, causing the local $\mathcal{v_c}$ and $\mathscr{G}$ vary along the crack front. The experimentally measured $\mathcal{v_c}$ and $\mathscr{G}$ data merely reflect collective, averaged behavior of the crack. This makes inversion for the intrinsic, “true” crack growth kinetics necessary, based upon the knowledge of the crack geometry. Simple and effective correction methods have been proposed and validated for the slow, chemical-reaction-controlled part (Region I) of the $\mathcal{v_c}$-$\mathscr{G}$ curve. However, reliable methods for the highly nonlinear, transport-dominated part (Region II) and its sudden transition to the dynamic propagation part (Region III) are still lacking. Here we propose a method for determining the intrinsic $\mathcal{v_c}$-$\mathscr{G}$ relationship across all three Regions based upon DT test data, using a simple model function and its numerical inversion. The performance of this approach is examined and demonstrated using both synthetic and laboratory data for subcritical crack growth in soda lime glass.

42 ENGINEERING↗

NASA NTRS1986

Direct numerical simulations of a temporally evolving mixing layer subject to forcing

The vortical evolution of mixing layers subject to various types of forcing is numerically simulated using pseudospectral methods. The effect of harmonic forcing and random noise in the initial conditions is examined with some results compared to experimental data. Spanwise forcing is found to enhance streamwise vorticity in a nonlinear process leading to a slow, secondary growth of the shear layer. The effect of forcing on a chemical reaction is favorably compared with experimental data at low Reynolds numbers. Combining harmonic and subharmonic forcing is shown to both augment and later destroy streamwise vorticity.

Claus, Russell W.↗

NASA NTRS1989

Direct numerical simulations of a temporally evolving mixing layer subject to forcing

Claus, R. W.↗

DOE OSTI2024

Nonlinear model predictive control for mode‐switching operation of reversible solid oxide cell systems

Abstract Solid oxide cells (SOCs) are a promising dual‐mode technology for the production of hydrogen through high‐temperature water electrolysis, and the generation of power through a fuel cell reaction that consumes hydrogen. Switching between these two modes as the price of electricity fluctuates requires reversible SOC operation and accurate tracking of hydrogen and power production set points. Moreover, a well‐functioning control system is important to avoid cell degradation during mode‐switching operation. In this article, we apply nonlinear model predictive control (NMPC) to an SOC module and supporting equipment and compare NMPC performance to classical proportional‐integral (PI) control strategies, while switching between the modes of hydrogen and power production. While both control methods provide similar performance across various metrics during mode switching, NMPC demonstrates a significant advantage in reducing cell thermal gradients and curvatures (mixed spatial‐temporal partial derivatives), thereby helping to mitigate long‐term degradation.

08 HYDROGEN↗

NASA NTRS2024

Nonlinear Dynamic Analyses of STS-1 Forward RCS Oxidizer Tank Structural Failure

At liftoff of the first Space Shuttle mission, a strong ignition overpressure pulse excited the stack. The frequency of the experienced pressure wave exceeded pre-launch predictions, and resulted in high, alternating normal accelerations along the length of the vehicle. As result of this unexpected loading, an oxidizer tank support strut in the orbiter’s forward reaction control system module failed in buckling under axial compression, which could have led to loss of the mission, crew, and vehicle. This incident is investigated in more detail using linear and nonlinear finite element analyses, complementing previous analyses with classical structural mechanics. The results of these analyses are examined to provide additional insight into the accident and to better inform future design decisions for modern space vehicle systems and structures.

Space Shuttle↗

NASA NTRS2024

Nonlinear Dynamic Analyses of STS-1 Forward RCS Oxidizer Tank Structural Failure

Space Shuttle↗

NASA NTRS2005

Updated Chemical Kinetics and Sensitivity Analysis Code

An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗

DOE OSTI2025

Simplified calculation of the area specific impedance for solid-state battery design

Simplified algebraic area specific impedance (ASI) correlations have been developed for solid-state composite battery electrodes made of a single ion conducting electrolyte, conductive additive, and intercalation active material. Two ASI expressions were developed, one for short times ( i.e. , pulsed power operation) and another for the pseudo steady state operation (i.e., sustained discharge for energy estimation). A full electrochemical model based on porous electrode theory was developed to examine the accuracy of the simplified ASI expressions. The simplified expressions agree favorably with full model results over a wide range of parameters (i.e. , electrode thicknesses, electrolyte conductivities, solid-state diffusion coefficients, specific surface areas, etc.) and conditions (i.e. , C-rates, states of charge, and pulse times). Under most conditions, the error between the full model and the correlations is well below 7 %. Higher errors were observed for the pseudo steady state expression at high/ low states of charge where the assumption of uniform reaction distributions loses validity. Here, the short time ASI has higher error at low states of charge due to the nonlinearity of the open circuit voltage equation, which is assumed linear in the formulation of the simplified algebraic expression.

25 ENERGY STORAGE↗

DOE OSTI2021

Similarity properties in capacitive radio frequency plasmas with nonlinear collision processes

Abstract Similarity laws (SL) are essential for correlating the characteristics of plasmas at different dimensional scales, which have been validated for radio frequency (rf) discharges at low pressure but under limited conditions. In this work, we evaluate the effects of nonlinear collisions (e.g. stepwise ionization) on the similarity properties in capacitive rf discharges in argon across a wide range of pressure regimes via fully kinetic particle-in-cell simulations. The SL scalings of fundamental discharge parameters, e.g. the electron power absorption and electron energy probability function, are examined with and without nonlinear collisions, respectively. Without the nonlinear collisions, the similarity scalings are found to be rigorously valid. When the nonlinear collisions are considered, the similarity properties in rf discharges still exist approximately, which indicates that the violations caused by the nonlinear collisions are still minor in the studied cases. The reasons for the effectiveness of SL scalings with nonlinear collision processes are also discussed. The results from this study confirmed the validity of similarity transformations with more complex reaction kinetics, which may promote the applicability of the SL scalings for the design and fabrications of plasma devices.

Physics↗

DOE OSTI2023

Long-term ultrasonic monitoring of concrete affected by alkali-silica reaction

This article presents continuous monitoring results of alkali-silica reaction (ASR) development in concrete specimens for over 400 days using ultrasonic testing and expansion measurements. Eight concrete specimens with nonreactive aggregate (Control), reactive coarse aggregate, and reactive fine aggregates were cast with two reinforced confinement conditions. The specimens were conditioned in an environmental chamber with high temperature and humidity (38°C and 90% relative humidity) to accelerate the ASR development. A multichannel ultrasonic monitoring system was developed to collect ultrasonic signals automatically, and the expansions in three directions were measured periodically. Results showed that the relative velocity change could detect the ASR initiation in all reactive specimens and show a correlation with expansion in the early stage. However, these correlations are inconsistent for different ASR specimens, and the velocity change becomes less sensitive to ASR damage in the late stage (after 300 days). Irrecoverable velocity drop was observed during every chamber shutdown period, especially in specimens with higher levels of ASR damage. This phenomenon suggests that the nonlinear ultrasonic response caused by the ambient temperature variation may indicate the ASR damage.

Engineering↗

DOE OSTI2022

Nonlinear burn control in ITER using adaptive allocation of actuators with uncertain dynamics

Abstract ITER will be the first tokamak to sustain a fusion-producing, or burning, plasma. If the plasma temperature were to inadvertently rise in this burning regime, the positive correlation between temperature and the fusion reaction rate would establish a destabilizing positive feedback loop. Careful regulation of the plasma’s temperature and density, or burn control, is required to prevent these potentially reactor-damaging thermal excursions, neutralize disturbances and improve performance. In this work, a Lyapunov-based burn controller is designed using a full zero-dimensional nonlinear model. An adaptive estimator manages destabilizing uncertainties in the plasma confinement properties and the particle recycling conditions (caused by plasma–wall interactions). The controller regulates the plasma density with requests for deuterium and tritium particle injections. In ITER-like plasmas, the fusion-born alpha particles will primarily heat the plasma electrons, resulting in different electron and ion temperatures in the core. By considering separate response models for the electron and ion energies, the proposed controller can independently regulate the electron and ion temperatures by requesting that different amounts of auxiliary power be delivered to the electrons and ions. These two commands for a specific control effort (electron and ion heating) are sent to an actuator allocation module that optimally maps them to the heating actuators available to ITER: an electron cyclotron heating system (20 MW), an ion cyclotron heating system (20 MW), and two neutral beam injectors (16.5 MW each). Two different actuator allocators are presented in this work. The first actuator allocator finds the optimal mapping by solving a convex quadratic program that includes actuator saturation and rate limits. It is nonadaptive and assumes that the mapping between the commanded control efforts and the allocated actuators (i.e. the effector model) contains no uncertainties. The second actuator allocation module has an adaptive estimator to handle uncertainties in the effector model. This uncertainty includes actuator efficiencies, the fractions of neutral beam heating that are deposited into the plasma electrons and ions, and the tritium concentration of the fueling pellets. Furthermore, the adaptive allocator considers actuator dynamics (actuation lag) that contain uncertainty. This adaptive allocation algorithm is more computationally efficient than the aforementioned nonadaptive allocator because it is computed using dynamic update laws so that finding the solution to a static optimization problem is not required at every time step. A simulation study assesses the performance of the proposed adaptive burn controller augmented with each of the actuator allocation modules.

Physics↗

DOE OSTI2023

Unexpected Thermomechanical Behavior of Off-Stoichiometry Epoxy/Amine Materials

Recent studies on off-stoichiometric thermosets reveal unique viscoelastic behavior derived from increased free volume and physical interactions between chain ends. To understand structural characteristics arising from cure and its effect on properties, we developed a Monte Carlo model based on step-growth polymerization. Our model accurately predicted structure–property trends for a two-component system of EPON 828 (EPON) and ethylenediamine. A second epoxy monomer, D.E.R. 732 (DER), was investigated to modulate T g . Binary mixtures of EPON and DER in off-stoichiometric, amine-rich formulations resulted in nonlinear evolution of thermomechanical properties with respect to initial formulation stoichiometry. Modifying our model with kinetic parameters allowing for differential epoxide/amine reaction kinetics only partially accounted for trends in T g , suggesting that spatiotemporal contributions-not captured by our model-were significant determinants of material properties compared to polymer architecture for three-component systems. Here these findings underpin the importance of spatial awareness in modeling to inform the development of dynamic thermosets.

36 MATERIALS SCIENCE↗