Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear approximation”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Development of a steady potential solver for use with linearized, unsteady aerodynamic analyses

A full potential steady flow solver (SFLOW) developed explicitly for use with an inviscid unsteady aerodynamic analysis (LINFLO) is described. The steady solver uses the nonconservative form of the nonlinear potential flow equations together with an implicit, least squares, finite difference approximation to solve for the steady flow field. The difference equations were developed on a composite mesh which consists of a C grid embedded in a rectilinear (H grid) cascade mesh. The composite mesh is capable of resolving blade to blade and far field phenomena on the H grid, while accurately resolving local phenomena on the C grid. The resulting system of algebraic equations is arranged in matrix form using a sparse matrix package and solved by Newton's method. Steady and unsteady results are presented for two cascade configurations: a high speed compressor and a turbine with high exit Mach number.

Hoyniak, Daniel↗
NASA NTRS1993

Techniques for designing rotorcraft control systems

This report summarizes the work that was done on the project from 1 Apr. 1992 to 31 Mar. 1993. The main goal of this research is to develop a practical tool for rotorcraft control system design based on interactive optimization tools (CONSOL-OPTCAD) and classical rotorcraft design considerations (ADOCS). This approach enables the designer to combine engineering intuition and experience with parametric optimization. The combination should make it possible to produce a better design faster than would be possible using either pure optimization or pure intuition and experience. We emphasize that the goal of this project is not to develop an algorithm. It is to develop a tool. We want to keep the human designer in the design process to take advantage of his or her experience and creativity. The role of the computer is to perform the calculation necessary to improve and to display the performance of the nominal design. Briefly, during the first year we have connected CONSOL-OPTCAD, an existing software package for optimizing parameters with respect to multiple performance criteria, to a simplified nonlinear simulation of the UH-60 rotorcraft. We have also created mathematical approximations to the Mil-specs for rotorcraft handling qualities and input them into CONSOL-OPTCAD. Finally, we have developed the additional software necessary to use CONSOL-OPTCAD for the design of rotorcraft controllers.

Levine, William S.↗
NASA NTRS1996

Computational Modeling of Ultrafast Pulse Propagation in Nonlinear Optical Materials

There is an emerging technology of photonic (or optoelectronic) integrated circuits (PICs or OEICs). In PICs, optical and electronic components are grown together on the same chip. rib build such devices and subsystems, one needs to model the entire chip. Accurate computer modeling of electromagnetic wave propagation in semiconductors is necessary for the successful development of PICs. More specifically, these computer codes would enable the modeling of such devices, including their subsystems, such as semiconductor lasers and semiconductor amplifiers in which there is femtosecond pulse propagation. Here, the computer simulations are made by solving the full vector, nonlinear, Maxwell's equations, coupled with the semiconductor Bloch equations, without any approximations. The carrier is retained in the description of the optical pulse, (i.e. the envelope approximation is not made in the Maxwell's equations), and the rotating wave approximation is not made in the Bloch equations. These coupled equations are solved to simulate the propagation of femtosecond optical pulses in semiconductor materials. The simulations describe the dynamics of the optical pulses, as well as the interband and intraband.

Goorjian, Peter M.↗
NASA NTRS2004

Shell Buckling Design Criteria Based on Manufacturing Imperfection Signatures

An analysis-based approach .for developing shell-buckling design criteria for laminated-composite cylindrical shells that accurately accounts for the effects of initial geometric imperfections is presented. With this approach, measured initial geometric imperfection data from six graphite-epoxy shells are used to determine a manufacturing-process-specific imperfection signature for these shells. This imperfection signature is then used as input into nonlinear finite-element analyses. The imperfection signature represents a "first-approximation" mean imperfection shape that is suitable for developing preliminary-design data. Comparisons of test data and analytical results obtained by using several different imperfection shapes are presented for selected shells. Overall, the results indicate that the analysis-based approach presented for developing reliable preliminary-design criteria has the potential to provide improved, less conservative buckling-load estimates, and to reduce the weight and cost of developing buckling-resistant shell structures.

Hilburger, Mark W.↗
NASA NTRS2017

Transition Induced by a Streamwise Array of Roughness Elements on a Supersonic Flat Plate

Roughness is unavoidable on practical high-speed vehicles, so it is critical to determine its impact on boundary layer transition. The flow field downstream of a streamwise array of cylindrical roughness elements is probed with hot-wire anemometry in this experiment. Mean flow distortion is examined in several measurement planes in the wake of the cylindrical roughness using the streak strength profiles and contour plots of the mass flux and total temperature. The roughness element heights and spacings were varied and their instability modes were examined. Cylindrical roughness elements approximately 140 micron tall produce an odd instability mode that grows weakly with downstream distance in the measurement range of this experiment. Cylindrical roughness elements approximately 280 micron tall produce an even instability mode that grows, becomes nonlinear, and then breaks down. Transition onset remains constant relative to the most downstream roughness in the streamwise array when the 280 micron roughness elements are spaced 2 diameters apart. Transition onset occurs at an earlier upstream location relative to the most downstream roughness in the streamwise array when the roughness elements are spaced 4 diameters appear to recover before the next downstream roughness element, so the location of transition shifts with the location of the most downstream roughness element in the array. When the rough- apart. The wake behind roughness elements spaced 2 diameters apart do not ness elements are spaced 4 diameters apart, the flow behind the first roughness element has enough space to recover before feeding into the second roughness element, and thus, moves transition forward.

Chou, Amanda↗
DOE OSTI2023

Recent Advances in PyROS: The Pyomo Solver for Two-Stage Nonconvex Robust Optimization

The slides present recent algorithmic and implementation advances of the two-stage robust optimization (RO) solver PyROS, and a benchmarking study which demonstrates the utility of PyROS for two-stage RO problems. The advances include extensions of the scope of PyROS to models with uncertain variable bounds, improvements to the initializations of the subproblems used by the underlying cutting set algorithm, and extensions of the uncertainty set interfaces. The benchmarking study is performed on a library of over 8,500 instances, with variations in the nonlinearities, degree-of-freedom partitioning, uncertainty sets, and polynomial decision rule approximations. Overall, the results highlight the effectiveness of PyROS for obtaining robust solutions to optimization problems with uncertain equality constraints.

Sherman, Jason↗
DOE OSTI2023

Recent Advances in PyROS: The Pyomo Solver for Two-Stage Nonconvex Robust Optimization

The slides present recent algorithmic and implementation advances of the two-stage robust optimization (RO) solver PyROS, and a benchmarking study which demonstrates the utility of PyROS for two-stage RO problems. The advances include extensions of the scope of PyROS to models with uncertain variable bounds, improvements to the initializations of the subproblems used by the underlying cutting set algorithm, and extensions of the uncertainty set interfaces. The benchmarking study is performed on a library of over 8,500 instances, with variations in the nonlinearities, degree-of-freedom partitioning, uncertainty sets, and polynomial decision rule approximations. Overall, the results highlight the effectiveness of PyROS for obtaining robust solutions to optimization problems with uncertain equality constraints.

Sherman, Jason↗
NASA NTRS1990

Nonlinear damping model for flexible structures

The study of nonlinear damping problem of flexible structures is addressed. Both passive and active damping, both finite dimensional and infinite dimensional models are studied. In the first part, the spectral density and the correlation function of a single DOF nonlinear damping model is investigated. A formula for the spectral density is established with O(Gamma(sub 2)) accuracy based upon Fokker-Planck technique and perturbation. The spectral density depends upon certain first order statistics which could be obtained if the stationary density is known. A method is proposed to find the approximate stationary density explicitly. In the second part, the spectral density of a multi-DOF nonlinear damping model is investigated. In the third part, energy type nonlinear damping model in an infinite dimensional setting is studied.

Zang, Weijian↗
NASA NTRS1988

Three-dimensional inelastic approximate analysis code (MOMM)

The Mechanics of Materials Model (MOMM) is one of a series of new stand-alone three dimensional nonlinear structural analysis codes. Incorporation of a general purpose finite element computer code into the hot section design process was severely limited by the high costs involved. MOMM is a stiffness method finite element code that uses an internally generated network of beams to characterize hot section component behavior. The method was proposed as a fast, easy to use, computationally efficient tool for approximate analyses. MOMM incorporates a wide variety of analysis capabilities, material models, and load type specifiers instrumental for the analysis of hot section components.

Meister, Jeffrey P.↗
NASA NTRS1988

On Hilbert-Schmidt norm convergence of Galerkin approximation for operator Riccati equations

An abstract approximation framework for the solution of operator algebraic Riccati equations is developed. The approach taken is based on a formulation of the Riccati equation as an abstract nonlinear operator equation on the space of Hilbert-Schmidt operators. Hilbert-Schmidt norm convergence of solutions to generic finite dimensional Galerkin approximations to the Riccati equation to the solution of the original infinite dimensional problem is argued. The application of the general theory is illustrated via an operator Riccati equation arising in the linear-quadratic design of an optimal feedback control law for a 1-D heat/diffusion equation. Numerical results demonstrating the convergence of the associated Hilbert-Schmidt kernels are included.

Rosen, I. G.↗
DOE OSTI2023

Bayesian High-Rank Hankel Matrix Completion for Nonlinear Synchrophasor Data Recovery

Phasor measurement units (PMUs) provide high temporal-resolution synchrophasor measurements for power system monitoring and control. The frequent data quality issues, such as missing and bad data, prevent the incorporation of synchrophasor data in real-time operations. Most existing data-driven data recovery methods assume the power system dynamics can be approximated by a linear dynamical system, and the recovery performance degrades significantly when the power system is experiencing nonlinear dynamics during significant events. Here, this paper proposes a data-driven Bayesian nonlinear synchrophasor data recovery method (Ba-NSDR) that can recover a consecutive time period of simultaneous data losses or errors across all channels, even when the underlying system is highly nonlinear. The idea is to lift the Hankel matrix of the spatial-temporal synchrophasor data to a higher dimension such that the lifted Hankel matrix is low-rank in that space and can be processed with the kernel trick. Our proposed Bayesian method then infers the probabilistic distributions of synchrophasor from the partial observations. Some distinctive features of Ba-NSDR include an uncertainty index to measure the accuracy of the recovery result and the robustness to parameter selections. Our method is verified on both synthetic and recorded event datasets.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1990

Time development of a perturbed-spherical nucleus in a pure supercooled liquid. I - Power-law growth of morphological instabilities. II - Nonlinear development

The linear growth stage of the morphological instabilities developing at the liquid-solid interface during crystal growth from an axisymmetric spherical nucleus is analyzed. The corresponding growth rate parameters are calculated numerically, and it is shown that morphological instabilities for free growth evolve according to a power law in agreement with WKB results and contrary to an exponential law found in a quasi-stationary approximation. The time evolution of an arbitrary perturbed liquid-solid interface from a linear into the nonlinear stage is studied. The initial perturbations include the eigenfunctions for the linear problem, localized perturbations of the grain-boundary type, and a stochastic noise. It is shown that the perturbations grow and spread in a wavelike manner. The formation of a predendritic growth stage is characterized by establishment of constant values of tip radius and velocity.

Pines, V.↗
DOE OSTI2019

A direct approach to estimate the anisotropy of protein structures from small-angle X-ray scattering

In the field of small-angle X-ray scattering (SAXS), the task of estimating the size of particles in solution is usually synonymous with the Guinier plot. The approximation behind this plot, developed by Guinier in 1939, provides a simple yet accurate characterization of the scattering behavior of particles at low scattering angle or momentum transfer q , together with a computationally efficient way of inferring their radii of gyration R G . Moreover, this approximation is valid beyond spherical scatterers, making its use ubiquitous in the SAXS world. However, when it is important to estimate further particle characteristics, such as the anisotropy of the scatterer's shape, no similar or extended approximations are available. Existing tools to characterize the shape of scatterers rely either on prior knowledge of the scatterers' geometry or on iterative procedures to infer the particle shape ab initio . In this work, a low-angle approximation of the scattering intensity I ( q ) for ellipsoids of revolution is developed and it is shown how the size and anisotropy information can be extracted from the parameters of that approximation. The goal of the approximation is not to estimate a particle's full structure in detail, and thus this approach will be less accurate than well known iterative and ab initio reconstruction tools available in the literature. However, it can be considered as an extension of the Guinier approximation and used to generate initial estimates for the aforementioned iterative techniques, which usually rely on R G and D max for initialization. This formulation also demonstrates that nonlinearity in the Guinier plot can arise from anisotropy in the scattering particles. Beyond ideal ellipsoids of revolution, it is shown that this approximation can be used to estimate the size and shape of molecules in solution, in both computational and experimental scenarios. The limits of the approach are discussed and the impact of a particle's anisotropy in the Guinier estimate of R G is assessed.

catalysis (homogeneous), catalysis (heterogeneous)↗
DOE OSTI2018

Benchmarking semiclassical approaches to strong-field QED: Nonlinear Compton scattering in intense laser pulses

The recoil associated with photon emission is key to the dynamics of ultrarelativistic electrons in strong electromagnetic fields, as found in high-intensity laser-matter interactions and astrophysical environments such as neutron star magnetospheres. When the energy of the photon becomes comparable to that of the electron, it is necessary to use quantum electrodynamics (QED) to describe the dynamics accurately. However, computing the appropriate scattering matrix element from strong-field QED is not generally possible due to multiparticle effects and the complex structure of the electromagnetic fields. Therefore, these interactions are treated semiclassically, coupling probabilistic emission events to classical electrodynamics using rates calculated in the locally constant field approximation. Here, we provide comprehensive benchmarking of this approach against the exact QED calculation for nonlinear Compton scattering of electrons in an intense laser pulse. We find agreement at the percentage level between the photon spectra, as well as between the models' predictions of absorption from the background field, for normalized amplitudes a0 > 5. We discuss possible routes towards improved numerical methods and the implications of our results for the study of QED cascades.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

To What Extent Can Raindrop Size Be Determined by a Multiple-Frequency Radar?

In this paper, an analytical treatment of the atmospheric remote sensing problem of determining the raindrop size distribution (DSD) with a spaceborne multifrequency microwave nadir-looking radar system is presented. It is typically assumed that with two radar measurements at different frequencies one ought to be able to calculate two state variables of the DSD: a bulk quantity, such as the rain rate, and a distribution shape parameter. To determine if this nonlinear problem can indeed be solved, the DSD is modeled as a Gamma distribution and quadratic approximations to the corresponding radar-rain relations are used to examine the invertibility of the resulting system of equations in the case of two as well as three radar frequencies. From the investigation, it is found that for regions of DSD state space multiple solutions exist for two or even three different frequency radar measurements. This should not be surprising given the nonlinear coupled nature of the problem.

rain characteristics↗
NASA NTRS1989

Kinetic effects of Alfven wave nonlinearity. I - Ponderomotive density fluctuations

The Vlasov theory is used to study kinetic corrections to fluid descriptions of Alfven wave nonlinearity. The method is to obtain an expression for the second-order perturbed distribution function produced by a nonlinear Alfven wave. From this distribution function a kinetically correct expression is obtained for the plasma density perturbation associated with an envelope-modulated Alfven wave. This kinetic theory result differs substantially from the fluid expression when the plasma beta is greater than about 1, and the electron and ion temperatures are approximately equal. This result is of interest because density fluctuations are an observationally accessible indicator of wave nonlinearity in solar system Alfven waves. It also will assist in the determination of properties of Alfven waves in the interstellar medium. Finally, this analysis also yields a kinetically correct expression for u, the magnetic field-aligned component of the plasma fluid velocity.

Spangler, Steven R.↗
DOE OSTI2022

Automated identification of dominant physical processes

The identification of processes that locally and approximately dominate dynamical system behavior has enabled significant advances in understanding and modeling nonlinear differential dynamical systems. Conventional methods of dominant process identification involve piecemeal and ad hoc (non-rigorous, informal) scaling analyses to identify dominant balances of governing equation terms and to delineate the spatiotemporal boundaries (boundaries in space and/or time) of each dominant balance. For the first time, we present an objective global measure of the fit of dominant balances to observations, which is desirable for automation, and was previously undefined. Furthermore, we propose a formal definition of the dominant balance identification problem in the form of an optimization problem. Here, we show that the optimization can be performed by various machine learning algorithms, enabling the automatic identification of dominant balances. Our method is algorithm agnostic and it eliminates reliance upon expert knowledge to identify dominant balances which are not known beforehand.

42 ENGINEERING↗