Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear modeling”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Separation of current density and electric field domains caused by nonlinear electronic instabilities

In 1963 Ridley postulated that under certain bias conditions circuit elements exhibiting a current- or voltage-controlled negative differential resistance will separate into coexisting domains with different current densities or electric fields, respectively, in a process similar to spinodal decomposition of a homogeneous liquid or disproportionation of a metastable chemical compound. The ensuing debate, however, failed to agree on the existence or causes of such electronic decomposition. Using thermal and chemical spectro-microscopy, we directly imaged signatures of current-density and electric-field domains in several metal oxides. The concept of local activity successfully predicts initiation and occurrence of spontaneous electronic decomposition, accompanied by a reduction in internal energy, despite unchanged power input and heat output. This reveals a thermodynamic constraint required to properly model nonlinear circuit elements. Our results explain the electroforming process that initiates information storage via resistance switching in metal oxides and has significant implications for improving neuromorphic computing based on nonlinear dynamical devices.

42 ENGINEERING↗
DOE OSTI2021

A nonlinear elasto-plastic bond model for the discrete element modeling of woody biomass particles

This work presents an experiment-informed, semi-empirical, elasto-plastic bond model for the discrete element modeling of woody biomass particles. The model renders nonlinear plastic deformation of materials when subjected to compression/tension, bending, and twisting, essential for accurately simulating the behavior of biomass in comminution. The model is implemented in an open-source DEM package LIGGGHTS and is assessed in a number of verification tests. The model is applied to simulate the fracturing test of notched loblolly pine blocks. Reduced bond strength is prescribed to the weak planes along the growth rings. It is found that the square-root formulation for describing the elasto-plastic normal bond contact delivers the best agreement with the experimental data over its linear and quadratic counterparts. Furthermore, the model predicts a strong anisotropic behavior of wood blocks with regard to the grain orientations, indicating that the cutting angle should be an important attribute for the optimization of biomass comminution.

42 ENGINEERING↗
DOE OSTI2016

How mesoscopic staircases condense to macroscopic barriers in confined plasma turbulence

This Rapid Communication sets forth the mechanism by which mesoscale staircase structures condense to form macroscopic states of enhanced confinement. Density, vorticity, and turbulent potential enstrophy are the variables for this model. Formation of the staircase structures is due to inhomogeneous mixing of (generalized) potential vorticity (PV). Such mixing results in the local sharpening of density and vorticity gradients. When PV gradients steepen, the density staircase structure develops into a lattice of mesoscale “jumps” and “steps,” which are, respectively, regions of local gradient steepening and flattening. The jumps then merge and migrate in radius, leading to the emergence of a new macroscale profile structure, so indicating that profile self-organization is a global process, which may be described by a local, but nonlinear model. This work predicts and demonstrates how mesoscale condensation of staircases leads to global states of enhanced confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

Element-specific modal formulations for large-displacement multibody dynamics

Large dispacement assumed-mode modeling techniques are examined in the context of multibody elastodynamics. The range of both general and element-specific approaches are studied with the aid of examples involving beams, plates, and shells. For systems undergoing primarily structural bending and twisting with little or no membrane distortion, it is found that fully-linear, element-specific, modal formulations provide the most accurate time history solutions at the least expense. When membrane effects become dominant in structural problems due to loading and boundary conditions, one must naturally resort to a formulation involving a nonlinear stress-strain relationship in addition to nonlinear terms associated with large overall system motion. Such nonlinear models were investigated using assumed modes and found to lead to modal convergence difficulties when standard free-free structural modes are employed. A constrained mode formulation aimed at addressing the convergence problem is proposed.

Ryan, R. R.↗
NASA NTRS2013

Sharp Eccentric Rings in Planetless Hydrodynamical Models of Debris Disks

Exoplanets are often associated with disks of dust and debris, analogs of the Kuiper Belt in our solar system. These "debris disks" show a variety of non-trivial structures attributed to planetary perturbations and utilized to constrain the properties of the planets. However, analyses of these systems have largely ignored the fact that, increasingly, debris disks are found to contain small quantities of gas, a component all debris disks should contain at some level. Several debris disks have been measured with a dust-to-gas ratio around unity where the effect of hydrodynamics on the structure of the disk cannot be ignored. Here we report that dust-gas interactions can produce some of the key patterns seen in debris disks that were previously attributed to planets. Through linear and nonlinear modeling of the hydrodynamical problem, we find that a robust clumping instability exists in this configuration, organizing the dust into narrow, eccentric rings, similar to the Fomalhaut debris disk. The hypothesis that these disks might contain planets, though thrilling, is not necessarily required to explain these systems.

debris disks↗
DOE OSTI2020

Rapid data-driven model reduction of nonlinear dynamical systems including chemical reaction networks using ℓ 1 -regularization

We develop a new data-driven paradigm for efficient model reduction of a broad class of nonlinear dynamical systems. Our model reduction method directly enables the interpretation of key components of the dynamical system, unlike traditional projection-based model reduction methods that focus on reducing computational complexity more than interpretability. Our method is not application specific and is simple to implement on nonlinear dynamical systems arising from a variety of different fields. It requires minimal parameterization using a single parameter to trade-off between model complexity and estimation error. We use a data-driven paradigm to formulate model reduction as an efficient convex optimization problem that scales polynomially in the original size of the complex system, enabling systems with as many as thousands of components to be reduced in a matter of minutes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1974

Modal analysis and control of flexible manipulator arms

The possibility of modeling and controlling flexible manipulator arms was examined. A modal approach was used for obtaining the mathematical model and control techniques. The arm model was represented mathematically by a state space description defined in terms of joint angles and mode amplitudes obtained from truncation on the distributed systems, and included the motion of a two link two joint arm. Three basic techniques were used for controlling the system: pole allocation with gains obtained from the rigid system with interjoint feedbacks, Simon-Mitter algorithm for pole allocation, and sensitivity analysis with respect to parameter variations. An improvement in arm bandwidth was obtained. Optimization of some geometric parameters was undertaken to maximize bandwidth for various payload sizes and programmed tasks. The controlled system is examined under constant gains and using the nonlinear model for simulations following a time varying state trajectory.

Neto, O. M.↗
DOE OSTI2017

Diffusion of New Products with Recovering Consumers

We consider the diffusion of new products in the discrete Bass-SIR model, in which consumers who adopt the product can later “recover” and stop influencing their peers to adopt the product. To gain insight into the effect of the social network structure on the diffusion, we focus on two extreme cases. In the “most-connected” configuration where all consumers are inter-connected (complete network), averaging over all consumers leads to an aggregate model, which combines the Bass model for diffusion of new products with the SIR model for epidemics. In the “least-connected” configuration where consumers are arranged on a circle and each consumer can only be influenced by his left neighbor (one-sided 1D network), averaging over all consumers leads to a different aggregate model which is linear, and can be solved explicitly. We conjecture that for any other network, the diffusion is bounded from below and from above by that on a one-sided 1D network and on a complete network, respectively. When consumers are arranged on a circle and each consumer can be influenced by his left and right neighbors (two-sided 1D network), the diffusion is strictly faster than on a one-sided 1D network. This is different from the case of non-recovering adopters, where the diffusion on one-sided and on two-sided 1D networks is identical. We also propose a nonlinear model for recoveries, and show that consumers’ heterogeneity has a negligible effect on the aggregate diffusion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Observations of modal coupling due to bolted joints in an experimental benchmark structure

The goal of this paper is to present a set of measurements from a benchmark structure containing two bolted joints to support future efforts to predict the damping due to the joints and to model nonlinear coupling between the first two elastic modes. Bolted joints introduce nonlinearities in structures, typically causing a softening in the natural frequency and an increase in damping because of frictional slip between the contact interfaces within the joint. These nonlinearities pose significant challenges when characterizing the response of the structure under a large range of load amplitudes, especially when the modal responses become coupled, causing the effective damping and natural frequency to not only depend on the excitation amplitude of the targeted mode, but also the relative amplitudes of other modes. In this work, two nominally identical benchmark structures, known in some prior works as the S4 beam, are tested to characterize their nonlinear properties for the first two elastic modes. Detailed surface measurements are presented and validated through finite element analysis and reveal distinct contact interactions between the two sets of beams. The free-free test structures are excited with an impact hammer and the transient response is analyzed to extract the damping and frequency backbone curves. A range of impact amplitudes and drive points are used to isolate a single mode or to excite both modes simultaneously. Differences in the nonlinear response correlate with the relative strength of the modes that are excited, allowing one to characterize mode coupling. Each of the beams shows different nonlinear properties for each mode, which is attributed to the different contact pressure distributions between the parts, although the mode coupling relationship is found to be consistent between the two. Lastly, the test data key finding are presented in this paper and the supporting data is available on a public repository for interested researchers.

42 ENGINEERING↗
DOE OSTI2015

Microwave Imaging Reflectometry for the study of Edge Harmonic Oscillations on DIII-D [Microwave Imaging Reflectometry (MIR) for the study of Edge Harmonic Oscillations (EHOs) on DIII-D]

Quiescent H-mode (QH) is an ELM free mode of operation in which edge-localized harmonic oscillations (EHOs) are believed to enhance particle transport, thereby stabilizing ELMs and preventing damage to the divertor and plasma facing components. Microwave Imaging Reflectometer (MIR) enabling direct comparison between the measured and simulated 2D images of density fluctuations near the edge can determine the 2D structure of density oscillation which can help to explain the physics behind EHO modes. MIR data sometimes indicates a counter-propagation between higher (n>1) and dominant (n=1) harmonics of coherent EHOs in the steep gradient regions of the pedestal. To preclude diagnostic artifacts, we have performed forward modeling that includes possible optical misalignments to show that offsets between transmitting and receiving antennas do not account for this feature. We have also simulated the non-uniform rotation of the EHO structure, which induces multiple harmonics that are properly characterized in the synthetic diagnostic. Excluding these possible explanations for the data, the counter-propagation observed in MIR data, which is not corroborated by external Mirnov coil array measurements, may be due to subtleties of the eigenmode structure, such as an inversion radius consistent with a magnetic island. Similar effects are observed in analysis of internal ECE-Imaging and BES data. Furthermore, the identification of a non-ideal structure motivates further exploration of nonlinear models of this instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

On the Boundary Conditions at an Oscillating Contact Line: A Physical/Numerical Experimental Program

We will pursue an improved physical understanding and mathematical model for the boundary condition at an oscillating contact line at high Reynolds number. We expect that the body force is locally unimportant for earth-based systems, and that the local behavior may dominate the mechanics of partially-filled reservoirs in the microgravity environment. One important space-based application for this contact-line study is for Faraday-waves. Oscillations in the direction of gravity (or acceleration) can dominate the fluid motion during take-off and reentry with large steady-state accelerations and in orbit, where fluctuations on the order of 10(exp -4)g occur about a zero mean. Our experience with Faraday waves has shown them to be 'cleaner' than those produced by vertical or horizontal oscillation of walls. They are easier to model analytically or computationally, and they do not have strong vortex formation at the bottom of the plate. Hence many, if not most, of the experiments will be performed in this manner. The importance of contact lines in the microgravity environment is well established. We will compare high resolution measurements of the velocity field (lO micro-m resolution) using particle-tracking and particle-image velocimetry as the fluid/fluid interface is approached from the lower fluid. The spatial gradients in the deviation provide additional means to determine an improved boundary condition and a measure of the slip region. Dissipation, the size of the eddy near the contact line, and hysteresis will be measured and compare to linear and nonlinear models of viscous and irrotational but dissipative models.

Perlin, Marc↗
NASA NTRS1994

Cyclogenesis in a saturated environment

The dynamics of baroclinic wave growth in a saturated environment is examined using linear and nonlinear models employing a parameterization of latent heat release that assumes all rising air is saturated, and saturation equivalent potential temperature is conserved on ascent. Piecewise potential vorticity (PV) diagnostics are used to interpret the results. When the stability to vertical displacements in saturated air is allowed to increase with height, as it must in an atmosphere with a constant, positive lapse rate of potential temperature, the growth rates of the most unstable modes of the Eady problem grow only marginally faster than the modes of the dry problem. The vertical variation of moist static stability produces a gradient of moist potential vorticity in the rising air, eliminating the short-wave cutoff present in the dry Eady problem. The destabilization of the short waves is shown to be associated with the interaction between surface potential temperature anomalies and diabatically generated lower-tropospheric potential vorticity anomalies. Nonlinear primitive equation simulations, starting from normal-mode initial conditions, show that while the dry wave grows at nearly the linear growth rate until maximum amplitude is reached, the moist wave grows significantly faster than the linear growth rate at finite amplitude. This enhanced growth is associated with the rapid amplification of a mesoscale PV anomaly generated by latent heat release at the warm front. The rapid amplification of the surface cyclone results from the superposition of the circulation associated with this misoscale PV anomaly upon the circulation associated with the surface and upper boundary potential temperature anomalies. Additional integrations with finite-amplitude initial conditions more typical of atmospheric conditions exhibit similar behavior. It is suggested that many of the rapid cyclogenesis events that occur as upper-troposheric PV anomalies cross the east coasts of continents may arise from the rapid generation of PV anomalies by condensational heating in the moist maritime lower troposphere.

Whitaker, Jeffrey S.↗
DOE OSTI2014

Communication: Modeling charge-sign asymmetric solvation free energies with nonlinear boundary conditions

We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory after replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for: (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley “bracelet” and “rod” test problems [D. L. Mobley, A. E. Barber II, C. J. Fennell, and K. A. Dill, “Charge asymmetries in hydration of polar solutes,” J. Phys. Chem. B 112, 2405–2414 (2008)]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry.

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

Physics constrained learning for data-driven inverse modeling from sparse observations

Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Furthermore, our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Nonlinear Reynolds stress models and the renormalization group

The renormalization group is applied to derive a nonlinear algebraic Reynolds stress model of anisotropic turbulence in which the Reynolds stresses are quadratic functions of the mean velocity gradients. The model results from a perturbation expansion that is truncated systematically at second order with subsequent terms contributing no further information. The resulting turbulence model applied to both low and high Reynolds number flows without requiring wall functions or ad hoc modifications of the equations. All constants are derived from the renormalization group procedure; no adjustable constants arise. The model permits inequality of the Reynolds normal stresses, a necessary condition for calculating turbulence-driven secondary flows in noncircular ducts.

Rubinstein, Robert↗
DOE OSTI2012

Chaotic saddles in nonlinear modulational interactions in a plasma

A nonlinear model of modulational processes in the subsonic regime involving a linearly unstable wave and two linearly damped waves with different damping rates in a plasma is studied numerically. We compute the maximum Lyapunov exponent as a function of the damping rates in a two-parameter space, and identify shrimp-shaped self-similar structures in the parameter space. By varying the damping rate of the low-frequency wave, we construct bifurcation diagrams and focus on a saddle-node bifurcation and an interior crisis associated with a periodic window. We detect chaotic saddles and their stable and unstable manifolds, and demonstrate how the connection between two chaotic saddles via coupling unstable periodic orbits can result in a crisis-induced intermittency. The relevance of this work for the understanding of modulational processes observed in plasmas and fluids is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Analysis of current-driven oscillatory dynamics of single-layer homoepitaxial islands on crystalline conducting substrates

Here we report results of a systematic study on the complex oscillatory current-driven dynamics of single-layer homoepitaxial islands on crystalline substrate surfaces and the dependence of this driven dynamical behavior on important physical parameters, including island size, substrate surface orientation, and direction of externally applied electric field. The analysis is based on a nonlinear model of driven island edge morphological evolution that accounts for curvature-driven edge diffusion, edge electromigration, and edge diffusional anisotropy. Using a linear theory of island edge morphological stability, we calculate a critical island size at which the island’s equilibrium edge shape becomes unstable, which sets a lower bound for the onset of time-periodic oscillatory dynamical response. Using direct dynamical simulations, we study the edge morphological dynamics of current-driven single-layer islands at larger-than-critical size, and determine the actual island size at which the migrating islands undergo a transition from steady to time-periodic asymptotic states through a subcritical Hopf bifurcation. At the highest symmetry of diffusional anisotropy examined, on {111} surfaces of face-centered cubic crystalline substrates, we find that more complex stable oscillatory states can be reached through period-doubling bifurcation at island sizes larger than those at the Hopf points. We characterize in detail the island morphology and dynamical response at the stable time-periodic asymptotic states, determine the range of stability of these oscillatory states terminated by island breakup, and explain the morphological features of the stable oscillating islands on the basis of linear stability theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗