Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear Analysis”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pyomo.DOE: An open-source package for model-based design of experiments in Python

Predictive mathematical models are a cornerstone of science and engineering. Yet selecting, calibrating, and validating said science-based models often remains an art in practice. Model-based design of experiments (MBDoE) provides a systematic framework to maximize information gain from experiments while minimizing time and resource costs. But MBDoE remains limited to niche application areas, in part because practitioners must integrate expertise in statistics, computational optimization, and modeling. To help reduce this barrier, we introduce Pyomo.DOE, an open-source package for MBDoE. Pyomo.DOE uses a nonlinear sensitivity analysis code k_aug to quickly approximate the Fisher information matrix and leverages a new stochastic programming abstraction. We demonstrate Pyomo.DOE with the first application of MBDoE to fixed-bed breakthrough experiments, which highlights the power of Pyomo.DOE to quantify the value of experimental modifications a priori for large-scale partial differential-algebraic equation (PDAE) models. Here we also provide a mathematical primer on MBDoE targeted at general chemical engineers.

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

Convergence and asymptotical stability of numerical solutions for neutral stochastic delay differential equations driven by G-Brownian motion

There are few numerical analysis results for nonlinear neutral stochastic delay differential equations driven by G-Brownian motion (G-NSDDEs). This paper is concerned with the numerical solutions of the G-NSDDEs to fill this gap. In this paper, we first devote to show that the stochastic theta numerical solution converges to the exact solution for the G-NSDDE. We then prove that the backward Euler–Maruyama numerical solution for the G-NSDDE is asymptotically mean-square stable under suitable conditions. Moreover, a numerical example is demonstrated to illustrate the effectiveness of our obtained results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Shape optimization under uncertainty for rotor blades of horizontal axis wind turbines

We present a computational framework for the shape optimization of a Horizontal-Axis Wind Turbine (HAWT) rotor blade under uncertainty. Our framework integrates aerodynamic simulations based on the blade element method which utilizes reduced order models of the blade structure and wind load with design sensitivity analysis and nonlinear programming. The wind velocity is modeled as a stochastic process to account for variations in time and space. An additional stochastic process accounts for uncertainties in the structural material properties. The uncertainty propagation is based on a non-intrusive polynomial chaos expansion that allows accurate estimation of stochastic performance metrics such as generated power and structural compliance. Sensitivities of cost and constraint functions with respect to shape parameters namely twist angles are computed with an efficient scheme to enable gradient based optimization. To demonstrate the effect of uncertainty, designs obtained from optimization under uncertainty are compared to those obtained from deterministic optimization.

42 ENGINEERING↗
DOE OSTI2023

Dynamical responses of constrained pipe conveying fluids and its dependence on the modeling of the contact force

Accurately modeling the impact force used in the analysis of loosely constrained cantilevered pipes conveying fluid is imperative. If little information is known of the motion-limiting constraints used in experiments, the analysis of the system may yield inaccurate predictions. Here in this work, multiple forcing representations of the impact force are defined and analyzed for a cantilevered pipe that conveys fluid. Depending on the representation of the impact force, the dynamics of the pipe can vary greatly when only the stiffness of the constraints is known from experiments. Three gap sizes of the constraints are analyzed, and the representation of the impact force used to analyze the system is found to significantly affect the response of the pipe at each gap size. An investigation on the effects of the vibro-impact force representation is performed through using basin of attraction analysis and nonlinear characterization of the system’s response.

42 ENGINEERING↗
DOE OSTI2025

A Scale‐Adaptive Urban Hydrologic Framework: Incorporating Network‐Level Storm Drainage Pipes Representation

Abstract Below‐ground urban stormwater networks (BUSNs) significantly influence urban flood dynamics, yet their representation at the watershed or larger scales remains challenging. We introduce a scalable urban hydrologic framework that centers on a novel network‐level BUSN representation, balancing the needs for physical basis, parameter parsimony, and computational efficiency. Our framework conceptualizes an urban watershed into four interacting zones: hillslopes (natural), storm‐sewersheds (urban), a sub‐network channel (tributaries), and a main channel. We develop an innovative Graph Theory‐based algorithm to derive network‐level BUSN parameters from publicly available datasets, enabling efficient, scalable parameterization. We demonstrate this framework's applicability at nine representative watersheds in the Houston metropolitan region, USA, with urban imperviousness ranging from 0% to 64% and drainage areas ranging from 24 to 302 . Our model achieves satisfying computational efficiency, completing hourly time step simulations for 18 years in less than 5 sec per watershed on a standard PC. Validation against observed daily streamflow confirms that the model can capture small‐to‐large flood peaks and seasonal and annual water balance over these watersheds. Comparisons with the National Water Model show better performance in predicting flood peaks and overall water balance, underscoring the promises of our new framework for urban hydrologic modeling at large scales. Furthermore, analysis reveals nonlinear relationships between BUSNs' designed capacities and flood reduction effects. Our approach bridges the gap between detailed hydraulic and large‐scale hydrologic models, providing a valuable tool for urban flood prediction and management across broader spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Magnetic reconnection mediated by hyper-resistive plasmoid instability

Magnetic reconnection mediated by the hyper-resistive plasmoid instability is studied with both linear analysis and nonlinear simulations. The linear growth rate is found to scale as S{sub H}{sup 1/6} with respect to the hyper-resistive Lundquist number S{sub H}≡L{sup 3}V{sub A}/η{sub H}, where L is the system size, V{sub A} is the Alfvén velocity, and η{sub H} is the hyper-resistivity. In the nonlinear regime, reconnection rate becomes nearly independent of S{sub H}, the number of plasmoids scales as S{sub H}{sup 1/2}, and the secondary current sheet length and width both scale as S{sub H}{sup −1/2}. These scalings are consistent with a heuristic argument assuming secondary current sheets are close to marginal stability. The distribution of plasmoids as a function of the enclosed flux ψ is found to obey a ψ{sup −1} power law over an extended range, followed by a rapid fall off for large plasmoids. These results are compared with those from resistive magnetohydrodynamic studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Attenuation due to hysteretic damage in the free vibration of a beam

We present an asymptotic analysis of nonlinear free vibration of a beam with a damage plane represented by nonlinear hysteretic bending and shear springs. The perturbation parameter is the product of the ratio of the nonlinear to linear parts of the stiffness times the amplitude of the free vibration. The loss of energy and ensuing attenuation due to hysteresis is accounted for by reducing the amplitude of vibration after each cycle by an amount such that the loss in total system energy equals the work done to traverse the hysteresis loop. A new Fourier representation for each cycle of the hysteresis and the deflection solution is used for this purpose and leads to higher harmonics, an evolving complex stiffness and corrected natural frequency that are linked to the attenuation. The frequency increases to its linear value from an initially reduced value. The damage parameter, frequency shift and fundamental amplitudes are presented as functions of the initial damage parameter and time (cycles of vibration). The amplitudes of several of the higher harmonics are also presented as functions of time. Many of the results exhibit sufficient sensitivity with respect to the damage parameter that they should be able to be used to characterize the damage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Augmenting the one-shot framework by additional constraints

The (multistep) one-shot method for design optimization problems has been successfully implemented for various applications. To this end, a slowly convergent primal fixed-point iteration of the state equation is augmented by an adjoint iteration and a corresponding preconditioned design update. In this paper we present a modification of the method that allows for additional equality constraints besides the usual state equation. Finally, a retardation analysis and the local convergence of the method in terms of necessary and sufficient conditions are given, which depend on key characteristics of the underlying problem and the quality of the utilized preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Impurity transport, turbulence transitions and intrinsic rotation in Alcator C-Mod plasmas

Linear and nonlinear gyrokinetic simulations are used to probe turbulent impurity transport in intrinsically rotating tokamak plasmas. For this simulation-based study, experimental input parameters are taken from a pair of ICRF heated Alcator C-Mod discharges exhibiting a change in the sign of the normalized toroidal rotation gradient at mid-radius (i.e. a change from hollow to peaked intrinsic rotation profiles). The simulations show that there is no change in the peaking of the calcium impurity between the plasmas with peaked and hollow rotation profiles, suggesting that the impurity transport and the shape of the rotation do not always change together. Furthermore, near mid-radius, r/a = 0.5 (normalized midplane minor radius), the linear and nonlinear gyrokinetic simulations exhibit no evidence of a transition from ion temperature gradient (ITG) to trapped electron mode dominance when the intrinsic rotation profile changes from peaked to hollow. Finally, extensive nonlinear sensitivity analysis is performed, and there is no change in the ITG critical gradient or in the stiffness of ion heat transport with the change in the intrinsic toroidal rotation profile shape, which suggests that the shape of the rotation profile is not dominated by the ITG onset in these cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

DEMNUni: the clustering of large-scale structures in the presence of massive neutrinos

We analyse the clustering features of Large Scale Structures (LSS) in the presence of massive neutrinos, employing a set of large-volume, high-resolution cosmological N-body simulations, where neutrinos are treated as separate collisionless particles. The volume of 8 h{sup -3} Gpc{sup 3}, combined with a resolution of about 8×10{sup 10}h{sup -1}M{sub ⊚} for the cold dark matter (CDM) component, represents a significant improvement over previous N-body simulations in massive neutrino cosmologies. In this work we focus, in the first place, on the analysis of nonlinear effects in CDM and neutrinos perturbations contributing to the total matter power spectrum. We show that most of the nonlinear evolution is generated exclusively by the CDM component. We therefore compare mildly nonlinear predictions from Eulerian Perturbation Theory (PT), and fully nonlinear prescriptions (HALOFIT) with the measurements obtained from the simulations. We find that accounting only for the nonlinear evolution of the CDM power spectrum allows to recover the total matter power spectrum with the same accuracy as the massless case. Indeed, we show that, the most recent version of the (HALOFIT) formula calibrated on ΛCDM simulations can be applied directly to the linear CDM power spectrum without requiring additional fitting parameters in the massive case. As a second step, we study the abundance and clustering properties of CDM halos, confirming that, in massive neutrino cosmologies, the proper definition of the halo bias should be made with respect to the cold rather than the total matter distribution, as recently shown in the literature. Here we extend these results to the redshift space, finding that, when accounting for massive neutrinos, an improper definition of the linear bias can lead to a systematic error of about 1-2 % in the determination of the linear growth rate from anisotropic clustering. This result is quite important if we consider that future spectroscopic galaxy surveys, as e.g. Euclid, are expected to measure the linear growth-rate with statistical errors less than about 3 % at z∼<1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Analysis of the linear and nonlinear stability of Alfven eigenmodes and fish-bones in JET DT discharges: mode identification and shear flows generation

The plasma in future nuclear fusion reactors will be heated by neutral beam injectors (NBIs) and high frequency electromagnetic waves as well as fusion born alpha particles. Energetic particles (EPs), with energies up to two orders of magnitude larger than the thermal plasma, can trigger EP driven modes and induce harmful EP losses, reducing the plasma heating efficiency and the economical viability of the reactor. The present study is dedicated to analyze the Alfven Eigenmode (AE) activity in JET D–T discharges, the closest experiment to reactor-like operation performed until now. There, EP driven modes are induced by the combined effect of tangential NBIs and ion cyclotron resonance heating (ICRH) driven EP. Linear and nonlinear simulations are performed with the gyro-fluid FAR3d code to analyze the AE activity observed in the discharge 99896. The linear simulations reproduce the unstable n = 3 to 5 toroidal AEs (TAE) at the inner plasma region observed in the experiment, triggered by highly energetic passing deuterium populations injected by the tangential NBIs, further accelerated by the effect of the ICRH up to 1 MeV. In addition, fish-bones triggered by energetic trapped hydrogen induced by the ICRH are also reproduced. On the other hand, the alpha particles density is too small to destabilize AEs in the experiment. Nonetheless, increasing artificially the alpha density by one order of magnitude, an n = 1 beta induced AE can be destabilized in the inner plasma region. Nonlinear simulations indicate the generation of zonal structures during the AE/fish-bone saturation phase. TAE and fish-bones causes a rather weak increase of the passing D and trapped H EP (around 2%), respectively. Shear flows and zonal currents are generated during the saturation of TAE and fish-bones. Nonlinear simulations performed for D–T and pure deuterium thermal plasma indicate AE/fish-bone activity is weaker and shear flows are less intense in the pure deuterium case, trends consistent with the experimental observations that also indicates a deterioration of the thermal plasma confinement. Therefore, both numerical studies and experimental evidence indicate the generation of shear flows by AE/fish-bones could be connected with an improvement of the thermal plasma confinement.

AE↗
DOE OSTI2025

Efficient Streaming Dynamic Mode Decomposition

We propose a reformulation of the streaming dynamic mode decomposition method that requires maintaining a single orthonormal basis, thereby reducing computational redundancy. The proposed efficient streaming dynamic mode decomposition method results in a constant-factor reduction in computational complexity and memory storage requirements. Numerical experiments on representative canonical dynamical systems show that the enhanced computational efficiency does not compromise the accuracy of the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

High Field Transport in (Ultra) Wide Bandgap Semiconductors: Diamond Versus Cubic GaN

We provide an analysis of nonlinear transport in diamond and cubic GaN (c-GaN) with emphasis on the different types of optical phonon scattering [i.e., optical deformation potential (ODP) scattering versus polar optic phonon (POP) scattering] limiting the carrier velocity. Both types of carrier mobilities and carrier saturation (peak) velocities in diamond and c-GaN as functions of different doping types and concentrations are obtained by directly solving the Boltzmann equation. Our model indicates that the nonrandomizing nature of POP scattering causes carrier temperature cooling, resulting in higher carrier drift velocity than with ODP scattering. This effect, in addition to the small carrier effective masses and large optical phonon energies, is responsible for the higher peak velocities in c-GaN, compared to carrier drift velocity in diamond.

42 ENGINEERING↗
DOE OSTI2015

Shear effects on energy dissipation from an elastic beam on a rigid foundation

This paper describes the energy dissipation arising from microslip for an elastic shell incorporating shear and longitudinal deformation resting on a rough-rigid foundation. This phenomenon is investigated using finite element (FE) analysis and nonlinear geometrically exact shell theory. Both approaches illustrate the effect of shear within the shell and observe a reduction in the energy dissipated from microslip as compared to a similar system neglecting shear deformation. In particular, it is found that the shear deformation allows for load to be transmitted beyond the region of slip so that the entire interface contributes to the load carrying capability of the shell. The energy dissipation resulting from the shell model is shown to agree well with that arising from the FE model, and this representation can be used as a basis for reduced order models that capture the microslip phenomenon.

42 ENGINEERING↗