Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinearly Consistent Implementation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Surrogate-Based Asynchronous Decomposition Technique for Realistic Security-Constrained Optimal Power Flow Problems

Here we present a decomposition approach for obtaining good feasible solutions for the security-constrained, alternating-current, optimal power flow (SC-AC-OPF) problem at an industrial scale and under real-world time and computational limits. The approach was designed while preparing and participating in ARPA-E’s Grid Optimization Competition (GOC) Challenge 1. The challenge focused on a near-real-time version of the SC-AC-OPF problem, where a base operating point is optimized, taking into account possible single-element contingencies, after which the system adapts its operating point following the response of automatic frequency droop controllers and voltage regulators. Our solution approach for this problem relies on state-of-the-art nonlinear programming algorithms, and it employs nonconvex relaxations for complementarity constraints, a specialized two-stage decomposition technique with sparse approximations of recourse terms and contingency ranking and prescreening. The paper describes and justifies our approach and outlines the features of its implementation, including functions and derivatives evaluation, warm-starting strategies, and asynchronous parallelism. We discuss the results of the independent benchmark of our approach by ARPA-E’s GOC team in Challenge 1, where it was found to consistently produce high-quality solutions across a wide range of network sizes and difficulty, and conclude by outlining future extensions of the approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Transitional Markov Chain Monte Carlo Sampler in UQTk

Transitional Markov Chain Monte Carlo (TMCMC) is a variant of a class of Markov Chain Monte Carlo algorithms known as tempering-based methods. In this report, the implementation of TMCMC in the Uncertainty Quantification Toolkit is investigated through the sampling of high-dimensional distributions, multi-modal distributions, and nonlinear manifolds. Furthermore, the Bayesian model evidence estimates obtained from TMCMC are tested on problems with known analytical solutions and shown to provide consistent results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Spectral analysis of shielded gamma ray sources using precalculated library data

In this paper, an approach has been developed for determining the intensity of a shielded source by first determining the thicknesses of three different shielding materials from a passively collected gamma-ray spectrum by making comparisons with predetermined shielded spectra. These evaluations are dependent on the accuracy and validity of the predetermined library spectra which were created by changing the thicknesses of the three chosen materials lead, aluminum and wood that are used to simulate any actual shielding. Each of the spectra produced was generated using MCNP5 with a sufficiently large number of histories to ensure a low relative error at each channel. The materials were held in the same respective order from source to detector, where each material consisted of three individual thicknesses and a null condition. This then produced two separate data sets of 27 total shielding material situations and subsequent predetermined libraries that were created for each radionuclide source used. The technique used to calculate the thicknesses of the materials implements a Levenberg–Marquardt nonlinear search that employs a tri-linear interpolation with the respective predetermined libraries within each channel for the supplied input unknown spectrum. Given that the nonlinear parameters require an initial guess for the calculations, the approach demonstrates first that when the correct values are input, the correct thicknesses are found. It then demonstrates that when multiple trials of random values are input for each of the nonlinear parameters, the average of the calculated solutions that successfully converges also produced the correct thicknesses. Under situations with sufficient information known about the detection situation at hand, the method was shown to behave in a manner that produces reasonable results and can serve as a good preliminary solution. Finally, this technique has the capability to be used in a variety of full spectrum inverse analysis problems including homeland security issues.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1999

Development of An Intelligent Flight Propulsion Control System

The initial design and demonstration of an Intelligent Flight Propulsion and Control System (IFPCS) is documented. The design is based on the implementation of a nonlinear adaptive flight control architecture. This initial design of the IFPCS enhances flight safety by using propulsion sources to provide redundancy in flight control. The IFPCS enhances the conventional gain scheduled approach in significant ways: (1) The IFPCS provides a back up flight control system that results in consistent responses over a wide range of unanticipated failures. (2) The IFPCS is applicable to a variety of aircraft models without redesign and,(3) significantly reduces the laborious research and design necessary in a gain scheduled approach. The control augmentation is detailed within an approximate Input-Output Linearization setting. The availability of propulsion only provides two control inputs, symmetric and differential thrust. Earlier Propulsion Control Augmentation (PCA) work performed by NASA provided for a trajectory controller with pilot command input of glidepath and heading. This work is aimed at demonstrating the flexibility of the IFPCS in providing consistency in flying qualities under a variety of failure scenarios. This report documents the initial design phase where propulsion only is used. Results confirm that the engine dynamics and associated hard nonlineaaities result in poor handling qualities at best. However, as demonstrated in simulation, the IFPCS is capable of results similar to the gain scheduled designs of the NASA PCA work. The IFPCS design uses crude estimates of aircraft behaviour. The adaptive control architecture demonstrates robust stability and provides robust performance. In this work, robust stability means that all states, errors, and adaptive parameters remain bounded under a wide class of uncertainties and input and output disturbances. Robust performance is measured in the quality of the tracking. The results demonstrate the flexibility of the IFPCS architecture and the ability to provide robust performance under a broad range of uncertainty. Robust stability is proved using Lyapunov like analysis. Future development of the IFPCS will include integration of conventional control surfaces with the use of propulsion augmentation, and utilization of available lift and drag devices, to demonstrate adaptive control capability under a greater variety of failure scenarios. Further work will specifically address the effects of actuator saturation.

Calise, A. J.↗
NASA NTRS1998

A Nonlinear Physics-Based Optimal Control Method for Magnetostrictive Actuators

This paper addresses the development of a nonlinear optimal control methodology for magnetostrictive actuators. At moderate to high drive levels, the output from these actuators is highly nonlinear and contains significant magnetic and magnetomechanical hysteresis. These dynamics must be accommodated by models and control laws to utilize the full capabilities of the actuators. A characterization based upon ferromagnetic mean field theory provides a model which accurately quantifies both transient and steady state actuator dynamics under a variety of operating conditions. The control method consists of a linear perturbation feedback law used in combination with an optimal open loop nonlinear control. The nonlinear control incorporates the hysteresis and nonlinearities inherent to the transducer and can be computed offline. The feedback control is constructed through linearization of the perturbed system about the optimal system and is efficient for online implementation. As demonstrated through numerical examples, the combined hybrid control is robust and can be readily implemented in linear PDE-based structural models.

Smith, Ralph C.↗
DOE OSTI2017

Anderson acceleration and application to the three-temperature energy equations

The Anderson acceleration method is an algorithm for accelerating the convergence of fixed-point iterations, including the Picard method. Anderson acceleration was first proposed in 1965 and, for some years, has been used successfully to accelerate the convergence of self-consistent field iterations in electronic-structure computations. Recently, the method has attracted growing attention in other application areas and among numerical analysts. Compared with a Newton-like method, an advantage of Anderson acceleration is that there is no need to form the Jacobian matrix. Thus the method is easy to implement. In this paper, an Anderson-accelerated Picard method is employed to solve the three-temperature energy equations, which are a type of strong nonlinear radiation-diffusion equations. Two strategies are used to improve the robustness of the Anderson acceleration method. One strategy is to adjust the iterates when necessary to satisfy the physical constraint. Another strategy is to monitor and, if necessary, reduce the matrix condition number of the least-squares problem in the Anderson-acceleration implementation so that numerical stability can be guaranteed. Numerical results show that the Anderson-accelerated Picard method can solve the three-temperature energy equations efficiently. Compared with the Picard method without acceleration, Anderson acceleration can reduce the number of iterations by at least half. A comparison between a Jacobian-free Newton–Krylov method, the Picard method, and the Anderson-accelerated Picard method is conducted in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Lattice cluster theory for dense, thin polymer films

While the application of the lattice cluster theory (LCT) to study the miscibility of polymer blends has greatly expanded our understanding of the monomer scale molecular details influencing miscibility, the corresponding theory for inhomogeneous systems has not yet emerged because of considerable technical difficulties and much greater complexity. Here, we present a general formulation enabling the extension of the LCT to describe the thermodynamic properties of dense, thin polymer films using a high dimension, high temperature expansion. Whereas the leading order of the LCT for bulk polymer systems is essentially simple Flory-Huggins theory, the highly non-trivial leading order inhomogeneous LCT (ILCT) for a film with L layers already involves the numerical solution of 3(L − 1) coupled, highly nonlinear equations for the various density profiles in the film. The new theory incorporates the essential “transport” constraints of Helfand and focuses on the strict imposition of excluded volume constraints, appropriate to dense polymer systems, rather than the maintenance of chain connectivity as appropriate for lower densities and as implemented in self-consistent theories of polymer adsorption at interfaces. The ILCT is illustrated by presenting examples of the computed profiles of the density, the parallel and perpendicular bonds, and the chain ends for free standing and supported films as a function of average film density, chain length, temperature, interaction with support, and chain stiffness. The results generally agree with expected general trends.

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

Experimental evidence for dynamic friction on rock fractures from frequency-dependent nonlinear hysteresis and harmonic generation: Low-Frequency Torsional Dynamic Friction

Frictional properties affect the propagation of high-amplitude seismic waves across rock fractures and faults. Laboratory evidence suggests that these properties can be measured in active seismic surveys, potentially offering a route to characterizing friction in situ. Here, we present experimental results from a subresonance torsional modulus and attenuation apparatus that utilizes micron-scale sinusoidal oscillations to probe the nonlinear stress-strain relation at a range of strain amplitudes and rates. Nonlinear effects are further quantified using harmonic distortion; however, time series data best illuminate underlying physical processes. The low-frequency stress-strain hysteretic loops show stiffening at the sinusoid's static ends, but stiffening is reduced above a threshold frequency. This shape is determined by harmonic generation in the strain; the stress signal has no harmonics, confirming that the fractured sample is the source of the nonlinearity. These qualitative observations suggest the presence of rate-dependent friction and are consistent between fractures in three different rock types. We propose that static friction at the low strain rate part of the cycle, when given sufficient “healing” time at low oscillation frequencies, causes this stiffening cusp shape in the hysteresis loop. While rate-and-state friction is commonly used to represent dynamic friction, it cannot capture static friction or negative slip velocities. So we implement another dynamic friction model, based on the work of Dahl, which describes this process and produces similar results. Since the two models have a similar form, parameterizations of field data could constraint fault model inputs, such as specific location velocity strengthening or weakening properties.

58 GEOSCIENCES↗
DOE OSTI2025

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Schwinger-Dyson equations in large-N quantum field theories and nonlinear random processes

We propose a stochastic method for solving Schwinger-Dyson equations in large-N quantum field theories. Expectation values of single-trace operators are sampled by stationary probability distributions of the so-called nonlinear random processes. The set of all the histories of such processes corresponds to the set of all planar diagrams in the perturbative expansions of the expectation values of singlet operators. We illustrate the method on examples of the matrix-valued scalar field theory and the Weingarten model of random planar surfaces on the lattice. For theories with compact field variables, such as sigma models or non-Abelian lattice gauge theories, the method does not converge in the physically most interesting weak-coupling limit. In this case one can absorb the divergences into a self-consistent redefinition of expansion parameters. A stochastic solution of the self-consistency conditions can be implemented as a 'memory' of the random process, so that some parameters of the process are estimated from its previous history. We illustrate this idea on the two-dimensional O(N) sigma model. The extension to non-Abelian lattice gauge theories is discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Coupling CTH to Linear Acoustic Propagation across an Air-Earth Interface

The interface between the Earth and the atmosphere forms a strong contrast in material properties. As such, numerical issues can arise when simulating an elastic wavefield across such a boundary when using a numerical simulation scheme. This is exacerbated when two different simulation codes are coupled straddling that interface. In this report we document how we implement the coupling of CTH, a nonlinear shock physics code, to a linearized elastic/acoustic wave propagation algorithm, axiElasti, across the air-earth interface. We first qualitatively verify that this stable coupling between the two algorithms produces expected results with no visible effects of the coupling interface. We then verify the coupling interface quantitatively by checking consistency with results from previous work and with coupled acoustic-elastic seismo-acoustic source inversions in three earth materials.

58 GEOSCIENCES↗
NASA NTRS1980

A variational technique for smoothing flight-test and accident data

The problem of determining aircraft motions along a trajectory is solved using a variational algorithm that generates unmeasured states and forcing functions, and estimates instrument bias and scale-factor errors. The problem is formulated as a nonlinear fixed-interval smoothing problem, and is solved as a sequence of linear two-point boundary value problems, using a sweep method. The algorithm has been implemented for use in flight-test and accident analysis. Aircraft motions are assumed to be governed by a six-degree-of-freedom kinematic model; forcing functions consist of body accelerations and winds, and the measurement model includes aerodynamic and radar data. Examples of the determination of aircraft motions from typical flight-test and accident data are presented.

Bach, R. E., Jr.↗
NASA NTRS1982

A Variational Technique for Smoothing Flight-Test and Accident Data

The problem of determining aircraft motions along a trajectory is solved using a variational algorithm that generates unmeasured states and forcing functions, and estimates instrument bias and scale-factor errors. The problem is formulated as a nonlinear fixed-interval smoothing problem, and is solved as a sequence or linear two-point boundary value problems, using a sweep method. The algorithm has been implemented for use in flight-test and accident analysis. Aircraft motions are assumed to be governed by a six-degree-of-freedom kinematic model; forcing functions consist of body accelerations and winds, and the measurement model includes aerodynamic and radar data. Examples of the determination of aircraft motions from typical flight-test and accident data are presented.

Bach, Ralph E. , Jr.↗
NASA NTRS1994

Parameter estimation in a structural acoustic system with fully nonlinear coupling conditions

A methodology for estimating physical parameters in a class of structural acoustic systems is presented. The general model under consideration consists of an interior cavity which is separated from an exterior noise source by an enclosing elastic structure. Piezoceramic patches are bonded to or embedded in the structure; these can be used both as actuators and sensors in applications ranging from the control of interior noise levels to the determination of structural flaws through nondestructive evaluation techniques. The presence and excitation of patches, however, changes the geometry and material properties of the structure as well as involves unknown patch parameters, thus necessitating the development of parameter estimation techniques which are applicable in this coupled setting. In developing a framework for approximation, parameter estimation and implementation, strong consideration is given to the fact that the input operator is unbonded due to the discrete nature of the patches. Moreover, the model is weakly nonlinear. As a result of the coupling mechanism between the structural vibrations and the interior acoustic dynamics. Within this context, an illustrating model is given, well-posedness and approximations results are discussed and an applicable parameter estimation methodology is presented. The scheme is then illustrated through several numerical examples with simulations modeling a variety of commonly used structural acoustic techniques for systems excitations and data collection.

Banks, H. T.↗
DOE OSTI2021

Tusas: A fully implicit parallel approach for coupled phase-field equations

In this study, we develop a fully-coupled, fully-implicit approach for phase-field modeling of solidification in metals and alloys. Predictive simulation of solidification in pure metals and metal alloys remains a significant challenge in the field of materials science, as microstructure formation during the solidification process plays a critical role in the properties and performance of the solid material. Our simulation approach consists of a finite element spatial discretization of the fully-coupled nonlinear system of partial differential equations at the microscale, which is treated implicitly in time with a preconditioned Jacobian-free Newton-Krylov method. The approach is algorithmically scalable as well as efficient due to an effective preconditioning strategy based on algebraic multigrid and block factorization. We implement this approach in the open-source Tusas framework, which is a general, flexible tool developed in C++ for solving coupled systems of nonlinear partial differential equations. The performance of our approach is analyzed in terms of algorithmic scalability and efficiency, while the computational performance of Tusas is presented in terms of parallel scalability and efficiency on emerging heterogeneous architectures. We demonstrate that modern algorithms, discretizations, and computational science, and heterogeneous hardware provide a robust route for predictive phase-field simulation of microstructure evolution during additive manufacturing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Plasma equilibria with multiple ion species: Equations and algorithm

Axisymmetric equilibrium of a magnetically confined plasma with multiple ion species is considered. To describe hot plasmas with isothermal surfaces, we adopt a formulation consistent with zero poloidal ion flow. This formulation includes all magnetic field components and also correctly includes all effects of toroidal ion rotation. There are two free surface functions for each species and a third which is determined by a differential equation relating surface functions per species. We have developed and implemented an algorithm for the solution of the resulting nonlinear equations and found solutions with large charge and mass contrast among the ion species for both compact (r = 0 included) and annular (r = 0 excluded) domains. Our solution method allows for arbitrary domain shapes, includes far-field conditions, and treats any combination of electrically conducting or insulating walls. Appropriate surface functions are used to describe the transition from closed to open field plasma in a reasonable manner. Solutions for advanced fuel cycle fusion systems (both D-{sup 3}He and p-{sup 11}B) are presented to illustrate the power of the method. Finally, we briefly discuss the special issues associated with obtaining very elongated solutions and describe the algorithm for implementing these features.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

polymec

Polymec is a set of tools that can be used by a computational scientist/researcher to construct application software that solves nonlinear differential equations on unstructured grids (including point meshes) on parallel computers. It consists of a set of libraries written in C99 that can be linked into an application to give it access to parallel solvers and visualization capabilities. Interfaces to several DOE-funded open source libraries are provided as a convenience to re-searchers to aid the implementation of application software.

Johnson, Jeffery↗
DOE OSTI2013

Nonlinear diffusion acceleration for the multigroup transport equation discretized with S{sub N} and continuous FEM with rattlesnake

Nonlinear diffusion acceleration (NDA) can improve the performance of a neutron transport solver significantly especially for the multigroup eigenvalue problems. The high-order transport equation and the transport-corrected low-order diffusion equation form a nonlinear system in NDA, which can be solved via a Picard iteration. The consistency of the correction of the low-order equation is important to ensure the stabilization and effectiveness of the iteration. It also makes the low-order equation preserve the scalar flux of the high-order equation. In this paper, the consistent correction for a particular discretization scheme, self-adjoint angular flux (SAAF) formulation with discrete ordinates method (S{sub N}) and continuous finite element method (CFEM) is proposed for the multigroup neutron transport equation. Equations with the anisotropic scatterings and a void treatment are included. The Picard iteration with this scheme has been implemented and tested with RattleS{sub N}ake, a MOOSE-based application at INL. Convergence results are presented. (authors)

97 MATHEMATICS AND COMPUTING↗