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

Evaluating route to impact convergence of the harmonic balance method for piecewise-smooth systems

Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.

42 ENGINEERING↗
DOE OSTI2023

Physics informed deep neural network embedded in a chemical transport model for the Amazon rainforest

Secondary organic aerosols (SOA) are fine particles in the atmosphere, which interact with clouds, radiation and affect the Earth’s energy budget. SOA formation involves chemistry in gas phase, aqueous aerosols, and clouds. Simulating these chemical processes involve solving a stiff set of differential equations, which are computationally expensive steps for three-dimensional chemical transport models. Deep neural networks (DNNs) are universal function approximators that could be used to represent the complex nonlinear changes in aerosol physical and chemical processes; however, key challenges such as generalizability to extended time periods, preservation of mass balance, simulating sparse model outputs, and maintaining physical constraints have limited their use in atmospheric chemistry. Here, we develop an approach of using a physics-informed DNN that overcomes previous such challenges and demonstrates its applicability for the chemical formation processes of isoprene epoxydiol SOA (IEPOX-SOA) over the Amazon rainforest. The DNN is trained with data generated by simulating IEPOX-SOA over the entire atmospheric column, using the Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem). The trained DNN is then embedded within WRF-Chem to replace the computationally expensive default solver of IEPOX-SOA formation. The trained DNN predictions generalizes well with the default model simulation of the IEPOX-SOA mass concentrations and its size distribution (20 size bins) over several days of simulations in both dry and wet seasons. The embedded DNN reduces the computational expense of WRF-Chem by a factor of 2. Our approach shows promise in terms of application to other computationally expensive chemistry solvers in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Solving Eco-Driving Problems Using Indirect Collocation Method and Smooth Representation

Here, we discuss the eco-driving problem, considering both electric and conventional powertrains, and presents a pathway to solving it numerically using an indirect collocation method. Despite the low-order system dynamics, the piecewise fuel/efficiency map, gear shifting, and real-world traffic/road situations bring system discontinuities/switchings and pure state constraints into the problem formulation, which make the problem highly nonlinear and nontrivial to solve. This letter introduces smooth approximations to convert the original problem to an unconstrained (and penalized) smooth boundary-value problem. This approach eliminates the discussion of the switching structure and leads to a lightweight Newton-method-based solution procedure.

42 ENGINEERING↗
DOE OSTI2026

Robust Optimal Experimental Design of Infinite-Dimensional Bayesian Nonlinear Inverse Problems

Abstract. We consider robust optimal experimental design (ROED) for nonlinear Bayesian inverse problems governed by partial differential equations (PDEs). An optimal design is one that maximizes some utility quantifying the quality of the solution of an inverse problem. However, the optimal design is dependent on elements of the inverse problem such as the simulation model, the prior, or the measurement error model. ROED aims to produce an optimal design that is aware of the additional uncertainties encoded in the inverse problem and remains optimal even after variations in them. We follow a worst-case scenario approach to develop a new framework for robust optimal design of nonlinear Bayesian inverse problems. The proposed framework (a) is scalable and designed for infinite-dimensional Bayesian nonlinear inverse problems constrained by PDEs; (b) develops efficient approximations of the utility, namely the expected information gain; (c) employs eigenvalue sensitivity techniques to develop analytical forms and efficient evaluation methods of the gradient of the utility with respect to the uncertainties against which we wish to be robust; and (d) employs a probabilistic optimization paradigm that properly defines and efficiently solves the resulting combinatorial max-min optimization problem. The effectiveness of the proposed approach is illustrated for optimal sensor placement problem in an inverse problem governed by an elliptic PDE.

Chowdhary, Abhijit↗
DOE OSTI2020

A Computationally Efficient Linear Semi-Lagrangian Scheme for the Advection of Microphysical Variables in Cloud-Resolving Models

An efficient semi-Lagrangian (SL) numerical scheme for the advection of positive-definite microphysical variables in cloud-resolving models is proposed. The scheme is a linear combination of SL schemes of the first and second order. The proposed scheme is compared with two high-order, monotonic, nonlinear Eulerian schemes, in two 2D tests. In the first test, advection of a positive scalar field with large gradients is performed within complicated idealized steering flows. In the second test, a hail storm is simulated using the 2D Hebrew University Cloud Model, describing the interactions between 15 size distributions of different microphysical quantities (each distribution is described by 43 mass bins). The proposed scheme is computationally efficient, has a low numerical diffusion and a high level of mass conservation accuracy, and preserves the sum of multiple advected variables as well as the linear relationships between the advected scalar variables. Although the proposed SL scheme is not exactly positive definite, the negative values and their number, which may appear only in the case of very strong gradients, are negligible. The scheme produces results similar to those of nonlinear Eulerian advection schemes while reducing computation time by approximately a factor of 10.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1966

Interaction of plasmons.

Nonlinear interaction between plasma oscillation normal modes in hydrodynamic approximation

WAVE INTERACTION↗
NASA NTRS1974

System design optimization for a Mars-roving vehicle and perturbed-optimal solutions in nonlinear programming

Work in two somewhat distinct areas is presented. First, the optimal system design problem for a Mars-roving vehicle is attacked by creating static system models and a system evaluation function and optimizing via nonlinear programming techniques. The second area concerns the problem of perturbed-optimal solutions. Given an initial perturbation in an element of the solution to a nonlinear programming problem, a linear method is determined to approximate the optimal readjustments of the other elements of the solution. Then, the sensitivity of the Mars rover designs is described by application of this method.

Pavarini, C.↗
NASA NTRS1974

Analytical studies of ultrashort optical pulse propagation in a resonant absorbing medium

Resonant absorption in a laser pulse can be reduced to zero if the trailing part induces stimulated emission and returns the absorbed energy back to the pulse. This kind of effect is called Self-Induced Transparency (SIT). Analytic descriptions of short pulse propagation in a resonant medium, based on Slowly Varying Envelope Approximation (SVEA) and in linear media, have been extended in studies of SIT along two principal directions: First, to improve the SVEA by introducing Short Pulse Approximation (SPA); and second, to consider the effect of nonlinear refractive index.

Lee, C. T.↗
NASA NTRS1975

Thermodynamics and phase separation of dense fully-ionized hydrogen-helium fluid mixtures

The free energy of a hydrogen-helium fluid mixture is evaluated for the temperatures and densities appropriate to the deep interior of a giant planet such as Jupiter. The electrons are assumed to be fully pressure-ionized and degenerate. In this regime, an appropriate first approximation to the ionic distribution functions can be found by assuming hard sphere interactions. Corrections to this approximation are incorporated by means of the perturbation theory of Anderson and Chandler. Approximations for the three-body interactions and the nonlinear response of the electron gas to the ions are included. It is predicted that a hydrogen-helium mixture, containing 10% by number of helium ions, separates into hydrogen-rich and helium-rich phases below about 8000 K, at the pressures relevant to Jupiter (4-40 Megabars). It is also predicted that the alloy occupies less volume per ion than the separated phases. The equations of state and other thermodynamic derivatives are tabulated. Implications of these results are discussed.

Stevenson, D. J.↗
NASA NTRS1975

Rotorcraft derivative identification from analytical models and flight test data

A general procedure is presented for systematic development of rotorcraft models for use in systems identification, which includes fuselage and rotor degrees of freedom (DOF). Formulations for rigid blade flap and lag as well as the normal mode representation of an elastic blade are developed for hingeless and articulated rotor systems. The method of multiblade coordinates is used to obtain linear constant coefficient state variable models of various levels of approximation. Two of the approximate models, a 6 DOF, are identified from a nonlinear articulated helicopter computer simulation. The results demonstrate the accuracy attainable for each model. Advanced results outline the status of rotorcraft modeling and systems identification and indicate areas that require further investigation.

Molusis, J. A.↗
NASA NTRS1976

Gaussian ideal impulsive loading of rigid viscoplastic plates

The response of a thin, rigid, viscoplastic plate subjected to a spatially axisymmetric Gaussian ideal impulse loading was studied analytically. The Gaussian ideal impulse distribution instantaneously imparts a Gaussian initial velocity distribution to the plate, except at the fixed boundary. The plate deforms with monotonically increasing deflections until the initial dynamic energy is completely dissipated in plastic work. The simply supported plate of uniform thickness obeys the von Mises yield criterion and a generalized constitutive equation for rigid, viscoplastic materials. For the small deflection bending response of the plate, neglecting the transverse shear stress in the yield condition and rotary inertia in the equations of dynamic equilibrium, the governing system of equations is essentially nonlinear. A proportional loading technique, known to give excellent approximations of the exact solution for the uniform load case, was used to linearize the problem and obtain analytical solution in the form of eigenvalue expansions. The linearized governing equations required the knowledge of the collapse load of the corresponding static problem.

Hayduk, R. J.↗
NASA NTRS1981

The three dimensional spherical model for the AGCE

The development of accurate numerical model of the atmospheric general circulation experiment (AGCE) is discussed. The model will serve both as a design and diagnostic tool for the AGCE, as well as for conducting numerical experiments which otherwise cannot be performed by AGCE. The code in its final form will solve the complete three dimensional nonlinear Navier-Stokes energy equations with the Boussinesq approximation. The code will allow for any thermal boundary conditions and any external forcing in the form of rotation and body forces and will allow for variable thermodynamic coefficients.

Antar, B. N.↗
NASA NTRS1981

Multigrid solvers on parallel computers

Massively parallel computers, as considered in this investigation, are not yet available. However, a large-scale parallel computer cannot usefully be designed before the hypothetical algorithms which will employ it are studied. Most of the studies of parallel partial differential equations (PDE) solvers are based on solution techniques much slower (on sequential machines) than multigrid methods. Multigrid methods are highly parallelizable. Each of their processes can simultaneously be performed at all grid points. The present investigation is concerned with a preliminary exploration of the potential of multigrid, or, more generally, Multi-Level Adaptive Techniques (MLAT) on computers with many processors. Basic processes are considered, taking into account coarse-grid approximation, relaxation, coarse-grid corrections, full multigrid algorithms, nonlinear problems and eigenvalue problems, fine-to-coarse correction, and chains of problems. Details of parallel multigrid processing are also examined.

Brandt, A.↗
NASA NTRS1982

Coupled flap-lag-torsional dynamics of hingeless rotor blades in forward flight

The aeroelastic stability and response problem of the coupled flap-lag-torsional dynamics of an isolated hingeless rotor blade in forward flight is considered. Linear, quasi-steady aerodynamics below stall is included. The spatial dependence of the partial differential nonliner equations of motion is discretized using a multimodal Galerkin method. The nonlinear time dependent equilibrium position (steady state response) about which the equations are linearized is obtained by solving a sequence of linear periodic response problems (quasilinearization). Results illustrating blade behavior in forward flight for both soft-in-plane and stiff-in-plane designs are presented. Quasilinearization provides a clear indication of the cases when nonlinear terms due to moderate deflections are important. The results indicate that nonlinearities affect system stability much more than system response. This implies that for obtaining blade vibrations and loads, approximate analytical models based on linearized formulations could be employed, whereas for determining blade stability a more accurate nonlinear analysis would be required.

Friedmann, P. P.↗
NASA NTRS1983

Computational methods for ideal compressible flow

Conservative dissipative difference schemes for computing one dimensional flow are introduced, and the recognition and representation of flow discontinuities are discussed. Multidimensional methods are outlined. Second order finite volume schemes are introduced. Conversion of difference schemes for a single linear convection equation into schemes for the hyperbolic system of the nonlinear conservation laws of ideal compressible flow is explained. Approximate Riemann solvers are presented. Monotone initial value interpolation; and limiters, switches, and artificial dissipation are considered.

Vanleer, B.↗
NASA NTRS1986

Statistical approach to cubic Langmuir turbulence

Previous work on the cubic direct interaction approximation applied to the truncated (in Fourier space) cubically nonlinear Schroedinger equation model of Langmuir turbulence is extended to more modes. In the undriven, undamped case, excellent agreement between the statistical theory and a numerical ensemble of solutions of the dynamic equations is obtained. In the driven, damped case, satisfactory agreement is obtained provided the dynamic ensemble is limited to initial conditions in the basin of attraction.

Sun, G.-Z.↗