Search NASA⌕ Search

SEARCH · Search NASA

Results for “tangent linear”

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 109 records · Page 6

Reinforcement Learning via Gaussian Processes with Neural Network Dual Kernels

While deep neural networks (DNNs) and Gaussian Processes (GPs) are both popularly utilized to solve problems in reinforcement learning, both approaches feature undesirable drawbacks for challenging problems. DNNs learn complex non-linear embeddings, but do not naturally quantify uncertainty and are often data-inefficient to train. GPs infer posterior distributions over functions, but popular kernels exhibit limited expressivity on complex and high-dimensional data. Fortunately, recently discovered conjugate and neural tangent kernel functions encode the behavior of overparameterized neural networks in the kernel domain. We demonstrate that these kernels can be efficiently applied to regression and reinforcement learning problems by analyzing a baseline case study.We apply GPs with neural network dual kernels to solve reinforcement learning tasks for the first time. We demonstrate, using the well understood mountain-car problem, that GPs empowered with dual kernels perform at least as well as those using the conventional radial basis function kernel. Finally, we conjecture that by inheriting the probabilistic rigor of GPs and the powerful embedding properties of DNNs, GPs using NN dual kernels will empower future reinforcement learning models on difficult domains.

97 MATHEMATICS AND COMPUTING↗

Analytical performance prediction for large antennas

A linear mathematical modeling of the static and dynamic characteristics of a deployable offset wrap rib antenna was attempted using small motions in an effort to simulate the surface distortion resulting from thermal loading and transient oscillations from control. Both cyclic symmetry and model synthesis were used for parametric analysis. Although the ribs can be modeled as slender beams varying cross sections and equivalent linear material properties, the mesh poses a problem since its equivalent stress strain relations are highly nonlinear and sensitive to the biaxial prestressed state. Consequently, the macroscopic material properties were measured experimentally in the neighborhood of the design prestress. An equivalent tangent modulus technique was then adopted about this equilibrium state and assumed to be valid within a small range of incremental stresses.

El-Raheb, M.↗

A computational procedure to analyze metal matrix laminates with nonlinear lamination residual strains

An approximate computational procedure is described for the analysis of angleplied laminates with residual nonlinear strains. The procedure consists of a combination of linear composite mechanics and incremental linear laminate theory. The procedure accounts for initial nonlinear strains, unloading, and in-situ matrix orthotropic nonlinear behavior. The results obtained in applying the procedure to boron/aluminum angleplied laminates show that this is a convenient means to accurately predict the initial tangent properties of angleplied laminates in which the matrix has been strained nonlinearly by the lamination residual stresses. The procedure predicted initial tangent properties results which were in good agreement with measured data obtained from boron/aluminum angleplied laminates.

Chamis, C. C.↗

A computational procedure to analyze metal matrix laminates with nonlinear lamination residual strains

An approximate computational procedure is described for the analysis of angleplied laminates with residual nonlinear strains. The procedure consists of a combination of linear composite mechanics and incremental linear laminate theory. The procedure accounts for initial nonlinear strains, unloading, and in-situ matrix orthrotropic nonlinear behavior. The results obtained in applying the procedure to boron/aluminum angleplied laminates show that this is a convenient means to accurately predict the initial tangent properties of angleplied laminates in which the matrix has been strained nonlinearly by the lamination residual stresses. The procedure predicted initial tangent properties results which were in good agreement with measured data obtained from boron/aluminum angleplied laminates.

Chamis, C. C.↗

Space-Split Algorithm for Sensitivity Analysis of Discrete Chaotic Systems With Multidimensional Unstable Manifolds

Accurate approximations of the change of a system's output and its statistics with respect to the input are highly desired in computational dynamics. Ruelle's linear response theory provides breakthrough mathematical machinery for computing the linear response of chaotic dynamical systems. In this paper, we propose an algorithm for sensitivity analysis of discrete chaos with an arbitrary number of positive Lyapunov exponents. We combine the concept of perturbation space-splitting, which regularizes Ruelle's original expression, together with measure-based parameterization of the expanding subspace. We use these tools to rigorously derive trajectory-following recursive relations that converge exponentially fast, and construct a memory-efficient Monte Carlo scheme for derivatives of the output statistics. Thanks to the regularization and lack of simplifying assumptions on the system's behavior, our method is immune to the common problems of other popular methods such as the exploding tangent solutions and unphysical shadowing directions. Here, we provide a ready-to-use algorithm, analyze its complexity, and demonstrate several numerical examples of sensitivity computation using physically-inspired low-dimensional systems.

97 MATHEMATICS AND COMPUTING↗

Nonlinear steady incompressible lifting-surface analysis with wake roll-up

The problem of lifting surfaces in steady incompressible flow is considered in terms of an integral equation relating the potential discontinuity on wing and wake to the normal derivative of the potential (normal wash) on the lifting surface. The integral equation is approximated by a system of linear algebraic equations obtained by dividing the surfaces into small quadrilateral elements and by assuming the potential discontinuity and the normal wash to be constant within each element. The wake geometry is obtained through iteration by satisfying the condition that the velocity be tangent to the surface of the wake and that the potential discontinuity be constant along the streamlines. Numerical results are in good agreement with existing ones

Suciu, E. O.↗

Design of a Modular Monolithic Implicit Solver for Multi-Physics Applications

The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.

Carton De Wiart, Corentin↗

Differential surface models for tactile perception of shape and on-line tracking of features

Tactile perception of shape involves an on-line controller and a shape perceptor. The purpose of the on-line controller is to maintain gliding or rolling contact with the surface, and collect information, or track specific features of the surface such as edges of a certain sharpness. The shape perceptor uses the information to perceive, estimate the parameters of, or recognize the shape. The differential surface model depends on the information collected and on the a priori information known about the robot and its physical parameters. These differential models are certain functionals that are projections of the dynamics of the robot onto the surface gradient or onto the tangent plane. A number of differential properties may be directly measured from present day tactile sensors. Others may have to be indirectly computed from measurements. Others may constitute design objectives for distributed tactile sensors of the future. A parameterization of the surface leads to linear and nonlinear sequential parameter estimation techniques for identification of the surface. Many interesting compromises between measurement and computation are possible.

Hemami, H.↗

Stability and control of hypersonic waveriders

The longitudinal and lateral stability derivatives of caret-wing waveriders flying on-design are calculated using linear piston theory. The experimental work of Kipke is used to validate the lift and pitch stiffness coefficients. Comparison is also made to the analytical work done by Hui. The calculations are extended to off-design conditions with an attached shock wave using tangent wedge theory. The stability derivatives of a Mach 6 waverider are computed and conclusions about its stability are drawn.

Tarpley, Christopher↗

Damping in Pitch and Roll of Triangular Wings at Supersonic Speeds

A method is derived for calculating the damping coefficients in pitch and roll for a series of triangular wings and a restricted series of sweptback wings at supersonic speeds. The elementary "supersonic source" solution of the linearized equation of motion is used to find the potential function of a line of doublets, and the flows are obtained by surface distributions of these doublet lines. The damping derivatives for triangular wings are found to be a function of the ratio of the tangent of the apex angle to the tangent of the Mach angle. As this ratio becomes equal to and greater than 1.0 for triangular wings, the damping derivatives, in pitch and in roll, become constant. The damping derivative in roll becomes equal to one-half the value calculated for an infinite rectangular wing, and the damping derivative in pitch for pitching about the apex becomes equal to 3.375 times that of an infinite rectangular wing.

Brown, Clinton E↗

A nonlinear finite-element analysis of wings in steady incompressible flows with wake roll-up

The problem of lifting surfaces and complex aircraft configurations in steady incompressible flow is considered. For lifting surfaces the problem is formulated in terms of an integral equation relating the potential discontinuity on wing and wake to the normal derivative of the potential on the lifting surface. For complex configurations the problem is formulated in terms of an integral equation relating the potential to its normal derivative on the surface of the aircraft. The integral equation is approximated by a system of linear algebraic equations obtained by dividing the surfaces into small quadrilateral elements and by assuming the potential (or the potential discontinuity) and its normal derivative to be constant within each element. The wake geometry is obtained by iteration by satisfying the condition that the velocity be tangent to the surface of the wake and that the potential discontinuity be constant along the streamlines.

Suciu, E. O.↗

Spatially Quasi-Periodic Water Waves of Infinite Depth

Abstract We formulate the two-dimensional gravity-capillary water wave equations in a spatially quasi-periodic setting and present a numerical study of solutions of the initial value problem. We propose a Fourier pseudo-spectral discretization of the equations of motion in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on a torus. We adopt a conformal mapping formulation and employ a quasi-periodic version of the Hilbert transform to determine the normal velocity of the free surface. Two methods of time-stepping the initial value problem are proposed, an explicit Runge–Kutta (ERK) method and an exponential time-differencing (ETD) scheme. The ETD approach makes use of the small-scale decomposition to eliminate stiffness due to surface tension. We perform a convergence study to compare the accuracy and efficiency of the methods on a traveling wave test problem. We also present an example of a periodic wave profile containing vertical tangent lines that is set in motion with a quasi-periodic velocity potential. As time evolves, each wave peak evolves differently, and only some of them overturn. Beyond water waves, we argue that spatial quasi-periodicity is a natural setting to study the dynamics of linear and nonlinear waves, offering a third option to the usual modeling assumption that solutions either evolve on a periodic domain or decay at infinity.

Wilkening, Jon (ORCID:0000000327827596)↗

A simple nonlinear joint model

Hertzian contact theory is applied to a butt joint with specially mismatched bearing surfaces to devise a simple mathematical model of nonlinear axial force-displacement behavior in jointed members. Normalized tangent stiffness-force plots, for several values of a joint imperfection parameter, are presented for the sample case of solid structural members of circular section. The results illustrate the potential problem of high joint compliance at low axial-force levels, as well as the generally desirable stiffening and linearizing effects of preload. A nonlinear oscillator problem based on the static model is also formulated and solved to illustrate the effect of amplitude on natural frequency. As expected, natural frequency is low when amplitude is small. The results call attention to the important roles that tight tolerances and preload are expected to play in the design and fabrication of deployable and erectable truss-type space structures.

Fichter, W. B.↗

Chaotic /strange/ and periodic behavior in instability saturation by the oscillating two-stream instability

The nonlinear Schroedinger equation with linear growth and damping is truncated to three waves. The resulting system of nonlinear ordinary differential equations describes the excitation of linearly damped waves by the oscillating two-stream instability driven by a linearly unstable pump wave. This system represents a simple model for the nonlinear saturation of a linearly unstable wave. The model is examined analytically and numerically as a function of the dimensionless parameters of the system. It is found that the model can exhibit a wealth of characteristic dynamical behavior including stationary equilibria, Hopf bifurcations to periodic orbits, period doubling bifurcations, chaotic solutions characteristic of a strange attractor, tangent bifurcations from chaotic to periodic solutions, transient chaos, and hysteresis. Many of these features are shown to be explainable on the basis of one-dimensional maps. In the case of chaotic solutions, evidence for the presence of a strange attractor is provided by demonstrating Cantor set-like structure (i.e., scale invariance) in the surface of section.

Russell, D. A.↗

Theoretical lift and drag of thin triangular wings at supersonic speeds

A method is derived for calculating the lift and the drag due to lift of point-forward triangular wings and a restricted series of sweptback wings at supersonic speeds. The elementary or "supersonic sources" solution of the linearized equation of motion is used to find the potential function of a line of doublets. The flow about the triangular flat plate is then obtained by a surface distribution of these doublet lines. The lift-curve slope of triangular wings is found to be a function of the ratio of the tangent of the apex angle to the tangent of the Mach angle. As the apex angle approaches and becomes greater than the Mach angle, the lift coefficient of the triangular wing becomes equal to that of a two-dimensional supersonic airfoil at the same Mach number. The drag coefficient due to lift of triangular wings with leading edges well behind the Mach cone is shown to be close to that of elliptically loaded wings of the same aspect ratio in subsonic flight. The resultant force on wings with leading edges outside the Mach cone, however, is shown to act normal to the surfaces and thus an induced drag equal to the lift times the angle of attack is obtained.

Brown, Clinton E↗

Optimum receiver structures for phase-multiplexed modulations

A MAP (maximum a posteriori) estimation loop is derived for a QPSK signal and, by suitable approximation to the nonlinearity which arises as a consequence of the MAP theory, is reconfigured with practical realizations that are valid for high and low SNRs. In particular, it is shown that by approximating the hyperbolic tangent nonlinearity in the MAP estimation loop by the first two terms in its power series, an interesting practical realization of this loop results which applies at low SNRs. The error signal in this loop is formed by multiplying the error signal and lock detector output signal of a conventional biphase Costas loop. A generalized linear in-phase channel configuration is proposed which allows carrier reconstruction from an unbalanced QPSK signal at all ratios of data rates and powers in the two channels, even in the limit as these ratios simultaneously approach unity, i.e., balanced quadriphase.

Simon, M. K.↗

A matrix-free approach for finite-strain hyperelastic problems using geometric multigrid

This paper investigates matrix-free algorithms for problems in quasi-static finite-strain hyperelasticity. Iterative solvers with matrix-free operator evaluation have emerged as an attractive alternative to sparse matrices in the fluid dynamics and wave propagation communities because they significantly reduce the memory traffic, the limiting factor in classical finite element solvers. Specifically, we study different matrix-free realizations of the finite element tangent operator and determine whether generalized methods of incorporating complex constitutive behavior might be feasible. Furthermore, in order to improve the convergence behavior of iterative solvers, we also propose a method by which to construct level tangent operators and employ them to define a geometric multigrid preconditioner. Additionally, the performance of the matrix-free operator and the geometric multigrid preconditioner is compared to the matrix-based implementation with an algebraic multigrid (AMG) preconditioner on a single node for a representative numerical example of a heterogeneous hyperelastic material in two and three dimensions. We find that matrix-free methods for finite-strain solid mechanics are very promising, outperforming linear matrix-based schemes by two to five times, and that it is possible to develop numerically efficient implementations that are independent of the hyperelastic constitutive law.

42 ENGINEERING↗

WORM (Write One, Read Many)

WORM (Write One, Run Many) is an easy to use, cross platform, embedded and extensible, functional programming language designed to facilitate the creation of input-decks for computer codes that use standard ASCII text files for input. WORM makes it easy to create generic (yet, complex and powerful) reusable models. Additionally its nature allows for complex calculations and routines to be coded once and easily reused, further simplifying the creation of input decks. WORM (Write One, Run Many) is a powerful and versatile tool designed to improve the efficiency of today’s criticality safety analyst by allowing: + input decks for parametric studies to be created quickly and easily, + calculations and variables to be imbedded into any input deck, thus allowing for meaningful parameter specifications, + problems to be specified using any combination of units, and + complex mathematically defined models to be created. A very simple syntax is employed, and therefore the WORM is easy to learn. A WORM model is essentially a standard input deck with some of its numerical values replaced by WORM code. WORM code may include and evaluate the following mathematical operators and functions: addition, subtraction, multiplication, division, exponentiation, modulus, sine, cosine, tangent, arcsine, arccosine, arctangent, the natural logarithm, logarithm base 10, integer truncation, absolute value, and random number. Several common constants, e.g., pi, e, and Avogadros’s Number (both as 6.022e23 and 0.6022), are predefined in WORM. Additionally, many unit conversion factors are also predefined: millimeters, meters, inches, feet, yards, and mils to centimeters; kilograms, pounds, and ounces to grams; liters, milliliters, gallons, and fluid ounces to cubic centimeters; and angular degrees to radians. For parametric studies, WORM supports various shorthand list specifications: the explicit step size, linear interpolation, and logarithmic interpolation. The list notation sequentially assigns multiple values to a name. WORM creates an input deck for each value of the name. If multiple lists are used, WORM steps through each list individually, i.e., WORM creates input decks corresponding to each and every permutation of the list values. Additionally, a library of standard material definitions and Perl subroutines are included. Any one of these files can be incorporated into the subject model with a simple WORM read command. WORM is completely written in Perl, the Practical Extraction and Reporting Language. Perl is one of the most portable programming languages available today. As such, the WORM works on practically any computer platform.

Sartor, Raymond↗