Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear phenomena”

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 469 records · Page 26

Unitary Qubit Lattice Algorithms for Plasma Physics

This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An inherent instability in fiber-reinforced composite laminates

Presented is a theory to predict the room-temperature shape of unsymmetric laminates. The theory is an extension of classical lamination theory and includes the effects of geometric nonlinearities. It is shown that the saddle shape predicted by the classical theory can actually be an unstable shape and a laminate will assume a cylindrical shape instead. The phenomena is coupled to the length, width, and thickness of the laminate.

Hyer, M. W.↗

On the relation between 'mixing length' and 'direct interaction approximation' theories of turbulence

The capabilities of direct interaction approximations (DIA) and weak-coupling approximations (WCA) in modeling turbulence in plasma physics studies are examined. DIA and WCA formalisms are reviewed and a DIA expression for homogeneous incompressible hydrodynamics is simplified to derive an energy spectrum in an arbitrary wavenumber space. A solution is defined for the expression when applied to long wavelength phenomena generated by an instability. The equation is further generalized to yield a transport equation for modeling the spectrum generated by homogeneous turbulence in a linearly unstable system undergoing nonlinear interaction. Finally, the technique is applied to two-dimensional turbulence in a low pressure, weakly ionized plasma embedded in a homogeneous magnetic field and experiencing instabilities caused by electron density and electric potential gradients, such as found in the ionospheric E region.

Sudan, R. N.↗

Vertical solidification of dendritic binary alloys

Three numerical techniques are employed to analyze the influence of thermosolutal convection on defect formation in directionally solidified (DS) alloys. The finite-element models are based on the Boussinesq approximation and include the plane-front model and two plane-front models incorporating special dendritic regions. In the second model the dendritic region has a time-independent volume fraction of liquid, and in the last model the dendritic region evolves as local conditions dictate. The finite-element models permit the description of nonlinear thermosolutal convection by treating the dendritic regions as porous media with variable porosities. The models are applied to lead-tin alloys including DS alloys, and severe segregation phenomena such as freckles and channels are found to develop in the DS alloys. The present calculations and the permeability functions selected are shown to predict behavior in the dendritic regions that qualitatively matches that observed experimentally.

Heinrich, J. C.↗

Exploring the Capabilities of a Machine Learning Algorithm to Detect Space Weather-Significant Emerging Active Regions

Active regions are a source of various phenomena responsible for Space Weather disturbances; therefore, developing a technology for early warning about upcoming magnetic activity is crucial to mitigate its impact. However, observational limitations and the high nonlinearity of processes associated with the accumulation of magnetic flux and its interaction with the surrounding plasma during the emergence through the convection zone make early activity detection a challenging problem. To address these challenges, we developed a physics-driven machine learning model that allows us to detect active regions (ARs) before they become visible on the solar surface by analyzing the power spectra of acoustic oscillations observed by the SDO/HMI instrument. This study is based on a time series of Doppler shift maps of 31x31-degree areas tracked with the Carrington rotation rate for four days before and after the emergence. The Doppler shift time series are processed into the oscillation power maps for four frequency ranges and accompanied by line-of-sight magnetograms and the continuum intensity maps from SDO/HMI. The resulting data are converted into a 1D time series representing the mean temporal variations of these quantities. The redacted time series are used as input to predict AR emergence using the Long Short Term Memory (LSTM) method. The training of the LSTM model is based on 40 ARs, which includes an independent analysis for each sub region that exhibits AR emergence or remains quiet. The emergence of magnetic flux (defined as a decrease of the continuum intensity) was detected with the developed LSTM algorithm from 5 to 48 hours before the reported time by NOAA. The developed model is capable of pointing to the time and location of active region formation. In this presentation, we discuss reasons that impact how early in advance the model can identify the upcoming activity and the possibility of improving the current predictive skills and steps to transition to the operational forecast.

Heliophysics↗

Investigating the Role of Accident Tolerant Cladding on Source Term Reduction for High-Burnup PWRs Using MELCOR

The use of accident tolerant fuel (ATF) cladding can increase coping times during and beyond design basis accidents. While such gains may be incremental, they provide a margin that can potentially be recovered to enable high-burnup (HBU) operation. Realizing such a margin requires demonstrating that the combination of HBU and ATF has not led to an overall increase in source term. This study investigates the influence of cladding technology (Zr-based, Cr-coated Zr, and FeCrAl) and fuel cycle length (18 and 24 months) on radiological dose at the boundary of the exclusion zone for a four-loop pressurized water reactor to investigate whether ATF claddings can provide such benefits. We analyze a recovered large break loss-of-coolant accident scenario to investigate the impact of transient timescale on the benefits of such coping time increases. The simulations have been performed using the MELCOR and MELCOR Accident Consequence Code System codes. For the cases analyzed, increased fuel cycle length did not necessarily increase radionuclide release and hydrogen generation, as these were found to be sensitive to the core power distribution. Similarly, off-site dose consequence is dominated by short-lived radionuclides that tend to saturate earlier in the burnup, so higher burnup operation did not necessarily increase the source term for the phenomena and transients analyzed here. Delays in recovery of the lowpressure safety injection system increase hydrogen production and radionuclide release, especially between 780 s and 1620 s, due to the nonlinear oxidation and core degradation behavior. Results show that Cr-coated Zr enhances safety by delaying heatup and gap release. Here, when uncertainty propagation on oxidation properties is considered, FeCrAl exhibits the lowest overall radionuclide release and off-site dose throughout the spectrum. However, while the considered “base model” performance is superior under delayed injection scenarios, upper-bound cases display hydrogen generation risk comparable to the Zr-based cladding.

Accident Tolerant Fuel↗

A simulation study of cold electron heating by loss cone instabilities

Results are presented of a computer simulation study of electron cyclotron harmonic waves; an electrostatic finite size particle simulation code is used. The initial electron velocity distribution function is modeled as a hot ring perpendicular to the ambient magnetic field and a cold Maxwellian component. Nonlinear cyclotron resonance is found to be the saturation mechanism of the unstable waves and the heating mechanism of the cold electrons. The results can be related to certain ionospheric phenomena.

Ashour-Abdalla, M.↗

Fluid Flow in Crystal Growth: Analysis of the Verticle Bridgman and Floating Zone Process

Research is aimed at the fundamental understanding of the interactions of heat and mass transport and fluid mechanics in crystal growth from the melt, especially in how it pertains to experiments supported by the Microgravity Science and Applications program. Emphasis has been on studies of the vertical Bridgman system, on small-scale floating zone systems proposed for space flight, and on quantitative prediction of nonlinear transitions in melt/solid interface morphologies leading to the formation of dendrites in directional solidification. The analysis of vertical Bridgman growth has lead to the complete numerical analysis of transport phenomena in directional solidification and to a detailed comparison of calculations to the GaGe growth experiments of Wang and Witt. Calculations show the central role of the ampoule material in determining the radial temperature gradients in the melt and thus, the intensity of convection. Results indicate that diffusion-controlled growth can not be achieved on Earth with conventional growth systems, but are feasible in microgravity

Brown, R. A.↗

Ion mixing in the plasma sheet boundary layer by drift instabilities

The linear stability properties of collisionless drift instabilities are analyzed in a Harris equilibrium model of the plasma sheet boundary layer (PSBL). The strearmng ions with drift-type instabilities driven in the PSBL are considered. The fluid approximation leads to growth but predicts that the mode width approaches the gyroradius of the energetic ions. Thus an integral equation theory for the modes is developed taking into account that in the PSBL the curvature drift is weak compared with the grad-B drift. The exact wave particle resonance is kept in the nonlocal response functions. Plasma density, temperature, and magnetic gradient drift motions are taken into account. The drift modes produce an anomalous cross-field momentum transport mixing the PSBL ions on the time scale of tens of seconds. A nonlinear simulation is performed which shows the coalescence of the small scale, fast growing modes into large-scale vortices. The relation between these collective modes and plasma sheet transport phenomena is discussed including the comparison with the competing plasma mixing from single-particle stochasticity.

Horton, W.↗

Application of non-linear dynamics to the characterization of cardiac electrical instability

Beat-to-beat alternation in the morphology of the ECG has been previously observed in hearts susceptible to fibrillation. In addition, fibrillation has been characterized by some as a chaotic state. Period doubling phenomena, such as alternation, and the onset of chaos have been connected by non-linear dynamical systems theory. In this paper, we describe the use of a technique from nonlinear dynamics theory, the construction of a first return nap, to assess the susceptibility to fibrillation threshhold in canine experiments.

NASA Discipline Cardiopulmonary↗

On the Recoverability of Reactor Dynamics from Point Kinetics Data using SINDYc

Data driven models for reactor dynamics tend to face a few notable challenges. First, the broad range of timescales present across the various feedback mechanisms. Next, high standards for safety and importance of model performance in all operational domains. Finally, the correlations between state variables observed in most transients leads to difficulties when methods attempt to attribute certain dynamic phenomena to a particular cause. The current paper seeks to support efforts towards incorporating physics knowledge into one particular data-driven method for finding reactor dynamics called ``Sparse Identification of Nonlinear Dynamics with Control (SINDYc). The incorporation of physical knowledge into SINDYc will help address the challenges listed above by giving the model an initial understanding of the system being modeled. The demonstration of an exact representation of a set of point reactor dynamics equations in SINDYc is provided. Then, this allows for further discussion with mathematical justification as to why SINDYc, and other data-driven methods, may be difficult to apply to nuclear reactor dynamics due to correlations in the state variables. A particularly useful result of this work is the set of candidate functions required in SINDYc to exactly represent the reactor dynamics under a point approximation. In the future, SINDYc can be integrated with point models for reactor dynamics before being applied to the physical system to yield higher accuracy.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Phenomenological Model of Nonlinear Dynamics and Deterministic Chaotic Gas Migration in Bentonite: Experimental Evidence and Diagnostic Parameters

Understanding gas migration in compacted clay materials, e.g., bentonite and claystone, is important for the design and performance assessment of an engineered barrier system of a radioactive waste repository system, as well as many practical applications. Existing field and laboratory data on gas migration processes in low-permeability clay materials demonstrate the complexity of flow and transport processes, including various types of instabilities, caused by nonlinear dynamics of coupled processes of liquid–gas exchange, dilation, fracturing, fracture healing, etc., which cannot be described by classical models of fluid dynamics in porous media. We here show that the complexity of gas migration processes can be explained using a phenomenological concept of nonlinear dynamics and deterministic chaos theory. To do so, we analyzed gas pressure and gas influx (i.e., input) and outflux (i.e., output), recorded during the gas injection experiment in the compact Mx80-D bentonite sample, and calculated a set of the diagnostic parameters of nonlinear dynamics and chaos, such a global embedding dimension, a correlation dimension, an information dimension, and a spectrum of Lyapunov exponents, as well as plotted 2D and 3D pseudo-phase-space strange attractors, based on the univariate influx and outflux time series data. These results indicate the presence of phenomena of low-dimensional deterministic chaotic behavior of gas migration in bentonite. In particular, during the onset of gas influx in the bentonite core, before the breakthrough, the development of gas flow pathways is characterized by the process of chaotic gas diffusion. After the breakthrough, with inlet-to-outlet movement of gas, the prevailing process is chaotic advection. During the final phase of the experiment, with no influx to the sample, the relaxation pattern of gas outflux is resumed back to a process of chaotic diffusion. The types of data analysis and a proposed phenomenological model can be used to establish the basic principles of experimental data-gathering, modeling predictions, and a research design.

36 MATERIALS SCIENCE↗

Deterministic symbolic regression with derivative information: General methodology and application to equations of state

Symbolic regression methods simultaneously determine the model functional form and the regression parameter values by generating expression trees. Symbolic regression can capture the complexity of real–world phenomena but the use of deterministic optimization for symbolic regression has been limited due to the complexity of the search space of existing formulations. Herein we present a novel deterministic mixed–integer nonlinear programming formulation for symbolic regression that incorporates derivative constraints through auxiliary expression trees. By applying the chain rule to mathematical operations, binary expression trees are capable of representing the calculation of first and second derivatives. We apply this formulation to illustrative examples using derivative information to show increased model discrimination capability. In addition, we perform a case study of a thermodynamic equation of state to gain insight on valid functional forms with thermodynamics–based constraints on the first and second derivatives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetic effects of self-organization in low-temperature magnetized plasmas (Final Report for DE-SC0019045: 9/1/2018-8/31/2019)

The research goal of this Early Career Research Program project is to advance the science of low temperature magnetized plasmas using first-principles kinetic, e.g., solving the Boltzmann equation, models. A better understanding of the plasma transport properties, such as electrical resistivity and thermal conductivity must be obtained to better control and design magnetized plasma sources. A key physical phenomenon in low-temperature magnetized plasmas is the self-organization of plasma flows, e.g., coherent plasma structures and transition between different discharge modes. Understanding the physics of coherent plasma structures has been identified as one of the four important frontier topics in a recent workshop by the Office of Fusion Energy Science. These self-emerging plasma patterns results from the nonlinear coupling between plasma constituents (ions, electrons, and neutral atoms), electromagnetic fields, and materials. The physics of low-temperature magnetized plasmas is particularly complex because (i) the collisionless instabilities affect collisional phenomena (ionization, transport, etc.) and vice versa and (ii) the presence of materials (sheath, electron injection, etc.) influences the bulk plasma properties in low-temperature plasma (LTP) devices. Computational models of these plasma flows remain challenging because the plasma density can vary several orders of magnitudes, a wide range of temporal and spatial scales must be resolved, the phenomena are inherently three-dimensional, and various physical processes, such as wall interaction, instabilities, inelastic collisions, etc., must be simultaneously and self-consistently taken into account. Understanding the fundamental transport mechanisms in cross-field configurations enables additional controllability of plasma properties and electron energy distribution functions (EEDFs).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Forced response of a cantilever beam with a dry friction damper attached. I - Theory. II - Experiment

A theoretical and experimental study of the forced vibration response of a cantilevered beam with Coulomb damping nonlinearity is described. Viscous damping in the beam is neglected. Beam and dry friction damper configurations of interest for applications to turbine blade vibrations are considered. It is shown that the basic phenomena found by Dowell (1983) for a simply supported beam with an attached dry friction damper of specific geometry also apply to a cantilevered beam and a more general representation of the dry friction damper and its associated mass and stiffness.

Dowell, E. H.↗

Effects of Initial Geometric Imperfections On the Non-Linear Response of the Space Shuttle Superlightweight Liquid-Oxygen Tank

The results of an analytical study of the elastic buckling and nonlinear behavior of the liquid-oxygen tank for the new Space Shuttle superlightweight external fuel tank are presented. Selected results that illustrate three distinctly different types of non-linear response phenomena for thin-walled shells which are subjected to combined mechanical and thermal loads are presented. These response phenomena consist of a bifurcation-type buckling response, a short-wavelength non-linear bending response and a non-linear collapse or "snap-through" response associated with a limit point. The effects of initial geometric imperfections on the response characteristics are emphasized. The results illustrate that the buckling and non-linear response of a geometrically imperfect shell structure subjected to complex loading conditions may not be adequately characterized by an elastic linear bifurcation buckling analysis, and that the traditional industry practice of applying a buckling-load knock-down factor can result in an ultraconservative design. Results are also presented that show that a fluid-filled shell can be highly sensitive to initial geometric imperfections, and that the use a buckling-load knock-down factor is needed for this case.

Nemeth, Michael P.↗

The dynamic aspects of thermo-elasto-viscoplastic snap-through and creep buckling phenomena

Use of a mathematical model and solution methodology, to examine dynamic buckling and dynamic postbuckling behavior of shallow arches and spherical caps made of a realistic material and undergoing non-isothermal, elasto-viscoplastic deformation was examined. Thus, geometric as well as material type nonlinearities of higher order are included in this analysis. The dynamic stability problem is studied under impulsive loading and suddenly applied loading with loads of constant magnitude and infinite duration. A finite element model was derived directly from the incrementally formulated nonlinear shell equations, by using a tensor-oriented procedure. As an example of the results, the time history of the midspan displacement of a damped shallow circular arch is presented.

Riff, R.↗