Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear problems”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multi-element flow-driven spectral chaos (ME-FSC) method for uncertainty quantification of dynamical systems

The flow-driven spectral chaos (FSC) is a recently developed method for tracking and quantifying uncertainties in the long-time response of stochastic dynamical systems using the spectral approach. The method uses a novel concept called enriched stochastic flow maps as a means to construct an evolving finite-dimensional random function space that is both accurate and computationally efficient in time. In this paper, we present a multi-element version of the FSC method (the ME-FSC method for short) to tackle (mainly) those dynamical systems that are inherently discontinuous over the probability space. In ME-FSC, the random domain is partitioned into several elements, and then the problem is solved separately on each random element using the FSC method. Subsequently, results are aggregated to compute the probability moments of interest using the law of total probability. To demonstrate the effectiveness of the ME-FSC method in dealing with discontinuities and long-time integration of stochastic dynamical systems, four representative numerical examples are presented in this paper, including the Van-der-Pol oscillator problem and the Kraichnan-Orszag three-mode problem. Results show that the ME-FSC method is capable of solving problems that have strong nonlinear dependencies over the probability space, both reliably and at low computational cost.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Electromagnetic full-$f$ gyrokinetics in the tokamak edge with discontinuous Galerkin methods

Here, we present an energy-conserving discontinuous Galerkin scheme for the full-f electromagnetic gyrokinetic system in the long-wavelength limit. We use the symplectic formulation and solve directly for ∂A II /∂t, the inductive component of the parallel electric field, using a generalized Ohm’s law derived directly from the gyrokinetic equation. Linear benchmarks are performed to verify the implementation and show that the scheme avoids the Amp`ere cancellation problem. We perform a nonlinear electromagnetic simulation in a helical open-field-line system as a rough model of the tokamak scrape-off layer using parameters from the National Spherical Torus Experiment (NSTX). This is the first published nonlinear electromagnetic gyrokinetic simulation on open field lines. Comparisons are made to a corresponding electrostatic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

MAPS: the MFEM Anisotropic Plasma Solver

Simulating magnetically confined fusion plasmas presents a uniquely challenging problem due to the nonlinear anisotropic heat conduction. We introduce the MAPS (MFEM Anisotropic Plasma Solver) tool, which uses a high-order finite element method to compute transport solutions on unstructured meshes. We show results for a set of three 2-D verification tests, two of which demonstrate the expected convergence properties for various mesh resolutions and polynomial degrees. We then discuss the convergence rate for the third test.

Barnett, Rhea [ORNL] (ORCID:0000000317527979)↗
DOE OSTI2025

Genetic algorithm-based geometry calibration for dynamic compression x-ray diffraction experiments

An important component of dynamic compression x-ray diffraction (XRD) experiment analysis is geometry calibration: proper data interpretation requires knowledge of the precise detector position and orientation and, if the experiment involves a single-crystal sample, knowledge of the lattice orientation. The determination of these parameters in the arbitrary three-dimensional (3D) scattering geometries often present in dynamic compression facilities is challenging, as the associated optimization problem can be highly nonlinear, nonsmooth, and discontinuous. We present a genetic algorithm-based approach for performing dynamic compression XRD calibrations that overcomes these obstacles. We provide details regarding the image processing, algorithm implementation, and open-source software deployment and demonstrate the capability of the approach to calibrate the detector and crystal parameters in 3D geometries. Notably, we demonstrate the solver’s capacity to find the crystal orientation without a priori rotation constraints.

Brown, Nathan P. [Sandia National Laboratories (SN↗
DOE OSTI2026

Neoclassical transport and profile prediction in transport barriers

Strong gradient regions in tokamaks, such as the pedestal or internal transport barriers, are regions of reduced turbulence where neoclassical transport can play a dominant role. However, standard neoclassical transport theory assumes that the gradient length scales of density, temperature, and potential are of the order of the system size. In the pedestal, gradient length scales are much shorter and are measured to be of the order of the ion poloidal gyroradius. We present an extension of neoclassical theory that is applicable in transport barriers of large aspect ratio tokamaks. We show that particle and momentum transport are connected in such a way that a source of parallel momentum can drive a significant neoclassical ion particle flux. In strong gradient regions, density, electric potential, mean parallel flow, and ion temperature are shown to no longer be flux functions. Instead, they have a small but important poloidally varying piece that modifies the transport equations to lowest order. This introduces a nonlinearity in the transport problem through the coupling with quasineutrality that yields multiple co-existing solutions when solving for the plasma profiles. The different solutions could be connected to low and high transport states and jumps between solutions could provide a new neoclassical explanation for H-L back-transitions.

Trinczek, Silvia [Princeton Plasma Physics Laborat↗
DOE OSTI2021

Efficient ptychographic phase retrieval via a matrix-free Levenberg-Marquardt algorithm

The phase retrieval problem, where one aims to recover a complex-valued image from far-field intensity measurements, is a classic problem encountered in a range of imaging applications. Modern phase retrieval approaches usually rely on gradient descent methods in a nonlinear minimization framework. Calculating closed-form gradients for use in these methods is tedious work, and formulating second order derivatives is even more laborious. Additionally, second order techniques often require the storage and inversion of large matrices of partial derivatives, with memory requirements that can be prohibitive for data-rich imaging modalities. We use a reverse-mode automatic differentiation (AD) framework to implement an efficient matrix-free version of the Levenberg-Marquardt (LM) algorithm, a longstanding method that finds popular use in nonlinear least-square minimization problems but which has seen little use in phase retrieval. Furthermore, we extend the basic LM algorithm so that it can be applied for more general constrained optimization problems (including phase retrieval problems) beyond just the least-square applications. Since we use AD, we only need to specify the physics-based forward model for a specific imaging application; the first and second-order derivative terms are calculated automatically through matrix-vector products, without explicitly forming the large Jacobian or Gauss-Newton matrices typically required for the LM method. We demonstrate that this algorithm can be used to solve both the unconstrained ptychographic object retrieval problem and the constrained “blind” ptychographic object and probe retrieval problems, under the popular Gaussian noise model as well as the Poisson noise model. We compare this algorithm to state-of-the-art first order ptychographic reconstruction methods to demonstrate empirically that this method outperforms best-in-class first-order methods: it provides excellent convergence guarantees with (in many cases) a superlinear rate of convergence, all with a computational cost comparable to, or lower than, the tested first-order algorithms.

36 MATERIALS SCIENCE↗
NASA NTRS1972

The prediction of nonlinear longitudinal combustion instability in liquid propellant rockets

An analytical technique was developed to solve nonlinear longitudinal combustion instability problems. The analysis yields the transient and limit cycle behavior of unstable motors and the threshold amplitude required to trigger a linearly stable motor into unstable operation. The limit cycle waveforms were found to exhibit shock wave characteristics for most unstable engine operating conditions. A method of correlating the analytical solutions with experimental data was developed. Calculated results indicate that a second-order solution adequately describes the behavior of combustion instability oscillations over a broad range of engine operating conditions, but that higher order effects must be accounted for in order to investigate engine triggering.

Lores, M. E.↗
NASA NTRS1973

Nonlinear longitudinal combustion instability in rocket motors.

A new analytical technique for the solution of nonlinear longitudinal combustion instability problems in rocket combustors is developed. Using relatively little computation time, this technique is capable of predicting the transient and limit cycle behavior of the combustion instability oscillations as well as the disturbance amplitude required to trigger an instability in a linearly stable motor. The limit cycle waveforms are found to exhibit shock wave characteristics for most unstable engine operating conditions. It is shown that the characteristics of the resulting instability are independent of the nature of the initial disturbance and they depend solely upon the engine operating conditions and the characteristics of the unsteady combustion process.

Lores, M. E.↗
NASA NTRS1973

A finite element tensor approach to plate buckling and postbuckling.

A practical finite element method is presented for geometrically nonlinear (von Karman) plate problems. Symmetric strain energy tensors provide an efficient formulation and solution, and allow the effects of initial imperfections to be considered. An unbalanced force iteration is employed, with the advantage that the time consuming Gauss decomposition needs to be performed only once. The convergence is controlled by an automatically computed relaxation factor, and this makes the method applicable in advanced stages of non-linearity. Several different triangular elements are used in explicitly deriving the required tensors.

Vos, R. G.↗
NASA NTRS1974

Thermal regime and temperature stresses in bodies during thermoradiational heating

An approach is developed to the application of the perturbation method for the solution of problems with essential external nonlinearities, based on identification in the boundary condition of a small nonlinear complex which is considered a perturbing function. The solutions obtained in the first approximation with error of 1 to 2% in calculating the unsteady temperature fields are then used to determine the temperature stresses and deformations in solid bodies of classical form.

Chistopyanova, N. V.↗
NASA NTRS1975

Analysis of nonlinear structures via mode synthesis

An effective procedure for NASTRAN was developed that permits any number of substructures of any size to be synthesized for the purpose of developing normal modes of vibration of the complete structural system. The technique is extended to permit modal transient analysis of the subdivided system. This latter procedure permits the use of NASTRAN's ability to include nonlinear forces in the problem. The five-phase process is accomplished using standard NASTRAN rigid formats with problem-independent alter packages and DMAP sequences.

Gieseke, R. K.↗
NASA NTRS1974

Nonlinear theory of electron neutralization waves in ions beams with dissipation

An analytical theory of nonlinear neutralization waves generated by injection of electrons from a grid in the direction of a homogeneous ion beam of uniform velocity and infinite extension is presented. The electrons are assumed to interact with the ions through the self-consistent space charge field and by strong collective interactions, while diffusion in the pressure gradient is disregarded (zero-temperature approximation). The associated nonlinear boundary-value problem is solved in closed form by means of a von Mises transformation. It is shown that the electron gas moves into the ion space in the form of a discontinuous neutralization wave, which exhibits a periodic field structure (incomplete neutralization). This periodic wave structure is damped out by intercomponent momentum transfer - i.e., after a few relaxation lengths a quasi-neutral plasma results.

Wilhelm, H. E.↗