SEARCH · Engineering Papers

Results for “classic nonlinearity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations

We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Optomechanical interaction in fibre lasers with micro-optomechanical resonance structures

We report the results of studies of fibre lasers with passive modulators based on light-excited micro-optomechanical resonance structures (micro-oscillators). It is shown that in fibre lasers based on active fibres doped with rare-earth elements (Er, Er/Yb, Yb, Nd), the optomechanical interaction of laser radiation with micro-oscillators of various types (fibre-optic, microvolume) leads to self-oscillations of the characteristics of laser radiation at frequencies of relaxation oscillations and intermode beats synchronised with the frequencies of elastic eigenoscillations of micro-oscillators. It is found that in an ultra-long erbium – ytterbium fibre laser with a nonlinear mirror based on a microcantilever, laser photothermal excitation of the second mode of elastic transverse oscillations of the microcantilever makes it possible to perform passive mode locking exclusively due to Q-switching of the laser cavity. Pulsed lasing with a controlled repetition rate (∼76 kHz), a pulse duration of 2 – 5 μs and an output energy of 0.1 μJ per pulse is implemented. Based on a simplified physical model of the indicated fibre lasers with micro-oscillators, we have developed an approximate mathematical model describing the regimes of passive mode locking of fibre lasers with micro-oscillators that play the role of mirrors with a nonlinear reflection coefficient in the laser cavity. The prospects for the development and application of the considered laser systems are discussed. (control of laser radiation parameters)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1974

A Computing Method for Sound Propagation Through a Nonuniform Jet Stream

Understanding the principles of jet noise propagation is an essential ingredient of systematic noise reduction research. High speed computer methods offer a unique potential for dealing with complex real life physical systems whereas analytical solutions are restricted to sophisticated idealized models. The classical formulation of sound propagation through a jet flow was found to be inadequate for computer solutions and a more suitable approach was needed. Previous investigations selected the phase and amplitude of the acoustic pressure as dependent variables requiring the solution of a system of nonlinear algebraic equations. The nonlinearities complicated both the analysis and the computation. A reformulation of the convective wave equation in terms of a new set of dependent variables is developed with a special emphasis on its suitability for numerical solutions on fast computers. The technique is very attractive because the resulting equations are linear in nonwaving variables. The computer solution to such a linear system of algebraic equations may be obtained by well-defined and direct means which are conservative of computer time and storage space. Typical examples are illustrated and computational results are compared with available numerical and experimental data.

Padula, S. L.↗

NASA NTRS1980

Current understanding of the physics of type III solar radio bursts

One of the most exciting plasma physics investigations of recent years has been connected with the understanding of a new strong turbulent plasma state excited by propagating electron beams. This new state is initiated on the linear level by parametric instabilities (OTS, modulational, etc.) and results in a very dynamic state composed of collective clusters of modes called solitons, cavitons, spikons, etc. Introduction of these concepts into the classic beam-plasma interaction problem has rendered quasi-linear and weak turbulence theories inapplicable over most of the interesting parameter range, and helped explain many paradoxes connected with the propagation of beams in the laboratory and space. Following a brief review of these nonlinear notions, the means by which their application to type III solar radiobursts has revolutionized understanding of their propagation, radioemission and scaling properties and has guided the in situ observations towards a more complete understanding are demonstrated. A particular burst (May 16, 1971) is analyzed in detail and compared with numerical predictions.

Papadopoulos, K.↗

NASA NTRS1985

On Fofonoff's Mode

The classical barotropic ocean model of Veronis (1966) is reexamined to verify the possible existence of Fofonoff's (1954) steady free inertial mode by considering the balance where J is the Jacobian, eta is the streamfunction of the horizontal velocity, and R is a parameter of nonlinearity. A linear functional relationship between the potential vorticity and the streamfunction and inertial boundary layers along the ocean basin was suggested by Fofonoff. However, the numerical integration showed that the system is barotropically unstable. The typical instantaneous field is given. The statistical equilibrium state eta is also given. Although it resembles Fofonoff's solution, it is maintained by the Reynolds stress field induced by the eddies in addition to the wind stress and dissipation. The plot of eta versus R delta sq +y is given. It is shown that the Fofonoff's linear relation between the potential vorticity and the streamfunction does not hold uniformly over the entire basin. The above results show that the Fofonoff's solution can not be realized physically.

Merkine, L. O.↗

NASA NTRS2004

Deformations of Flat Unsymmetric Laminates Subjected to Inplane Loads

The geometrically nonlinear deformation response of initially flat unsymmetric cross-ply laminates subjected to an inplane compressive load and two sets of boundary conditions is studied. Stability of the deformations is considered. At issue is whether or not the plate remains flat with increased compressive loading, and whether it buckles. A semi-infinite unsymmetric cross-ply laminate is used to show the combined effects of geometric nonlinearities and bending-stretch coupling. Finite element results for finite laminates are then presented, and it is shown that to a large degree the boundary conditions control the character of the deformation response. It appears that clamped boundary conditions support buckling behavior, in the classic sense of bifurcation, whereas simply-supported conditions do not.

Hyer, Michael W.↗

DOE OSTI2020

Surface Energies and Structure of Salt–Brine Interfaces

Permeability of salt formations is controlled by the equilibrium between the salt–brine and salt–salt interfaces described by the dihedral angle, which can change with the composition of the intergranular brine. In this work, classical molecular dynamics (MD) simulations were used to investigate the structure and properties of the salt–brine interface to provide insight into the stability of salt systems. Mixed NaCl–KCl brines were investigated to explore differences in ion size on the surface energy and interface structure. Nonlinearity was noted in the salt–brine surface energy with increasing KCl concentration, and the addition of 10% KCl increased surface energies by 2–3 times (5.0 M systems). Size differences in Na+ and K+ ions altered the packing of dissolved ions and water molecules at the interface, impacting the surface energy. Additionally, ions at the interface had lower numbers of coordinating water molecules than those in the bulk and increased hydration for ions in systems with 100% NaCl or 100% KCl brines. Ultimately, small changes in brine composition away from pure NaCl altered the structure of the salt–brine interface, impacting the dihedral angle and the predicted equilibrium permeability of salt formations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Analytical WKB theory for high-harmonic generation and its application to massive Dirac electrons

Here we propose an analytical approach to high-harmonic generation (HHG) for nonperturbative low-frequency and high-intensity fields based on the (Jeffreys-)Wentzel-Kramers-Brillouin (WKB) approximation. By properly taking into account Stokes phenomena of WKB solutions, we obtain wave functions that systematically include the repetitive dynamics of production and acceleration of electron-hole pairs and quantum interference due to phase accumulation between different pair production times (Stückelberg phase). Using the obtained wave functions without relying on any phenomenological assumptions, we explicitly compute electric current (including intra- and interband contributions) as the source of HHG for a massive Dirac system in (1+1) dimensions under an ac electric field. We demonstrate that the WKB approximation agrees well with numerical results obtained by solving the time-dependent Schrödinger equation and point out that the quantum interference is important in HHG. We also predict in the deep nonperturbative regime that (1) harmonic intensities oscillate with respect to electric-field amplitude E 0 and frequency Ω with a period determined by the Stückelberg phase, (2) the cutoff order of HHG is determined by 2eE 0 ℏΩ 2 , with e being the electron charge, and that (3) noninteger harmonics, controlled by the Stückelberg phase, appear as a transient effect. Our WKB theory is particularly suited for a parameter regime, where the Keldysh parameter γ = (Δ/2)Ω/eE 0 ,with Δ being the gap size, is small. This parameter regime corresponds to intense lasers in the terahertz regime for realistic massive Dirac materials. Our analysis implies that the so-called HHG plateau can be observed at the terahertz frequency within the current technology.

36 MATERIALS SCIENCE↗

DOE OSTI2021

de Gennes Narrowing and Relationship between Structure and Dynamics in Self-Organized Ion-Beam Nanopatterning

Investigating the relationship between structure and dynamical processes is a central goal in condensed matter physics. Perhaps the most noted relationship between the two is the phenomenon of de Gennes narrowing, in which relaxation times in liquids are proportional to the scattering structure factor. Here, a similar relationship is discovered during the self-organized ion-beam nanopatterning of silicon using coherent x-ray scattering. However, in contrast to the exponential relaxation of fluctuations in classic de Gennes narrowing, the dynamic surface exhibits a wide range of behaviors as a function of the length scale, with a compressed exponential relaxation at lengths corresponding to the dominant structural motif—self-organized nanoscale ripples. These behaviors are reproduced in simulations of a nonlinear model describing the surface evolution. Our team suggests that the compressed exponential behavior observed here is due to the morphological persistence of the self-organized surface ripple patterns which form and evolve during ion-beam nanopatterning.

36 MATERIALS SCIENCE↗

NASA NTRS1988

Local buckling and crippling of composite stiffener sections

Local buckling, postbuckling, and crippling (failure) of channel, zee, and I- and J-section stiffeners made of AS4/3502 graphite-epoxy unidirectional tape are studied by experiment and analysis. Thirty-six stiffener specimens were tested statically to failure in axial compression as intermediate length columns. Web width is 1.25 inches for all specimens, and the flange width-to-thickness ratio ranges from 7 to 28 for the specimens tested. The radius of the stiffener corners is either 0.125 or 0.250 inches. A sixteen-ply orthotropic layup, an eight-ply quasi-isotropic layup, and a sixteen-ply quasi-isotropic layup are examined. Geometrically nonlinear analyses of five specimens were performed with the STAGS finite element code. Analytical results are compared to experimental data. Inplane stresses from STAGS are used to conduct a plane stress failure analysis of these specimens. Also, the development of interlaminar stress equations from equilibrium for classical laminated plate theory is presented. An algorithm to compute high order displacement derivatives required by these equations based on the Discrete Fourier Transform (DFT) is discussed.

Bonanni, David L.↗

NASA NTRS1992

Compression failure of angle-ply laminates

Test results from the compression loading of (+ or - Theta/ - or + Theta)(sub 6s) angle-ply IM7-8551-7a specimens, 0 less than or = Theta less than or = 90 degs, are presented. The observed failure strengths and modes are discussed, and typical stress-strain relations shown. Using classical lamination theory and the maximum stress criterion, an attempt is made to predict failure stress as a function of Theta. This attempt results in poor correlation with test results and thus a more advanced model is used. The model, which is based on a geometrically nonlinear theory, and which was taken from previous work, includes the influence of observed layer waviness. The waviness is described by the wave length and the wave amplitude. The theory is briefly described and results from the theory are correlated with test results. It is shown that by using levels of waviness observed in the specimens, the correlation between predictions and observations is good.

Peel, L. D.↗

DOE OSTI2022

Double copy for massive scalar field theories

We explore extensions of the double copy to massive theories and find a new cubic theory with a local double copy. We consider the nonlinear sigma model and the special galileon theory, massless versions of which are known to be related through the double copy. We show that by performing a Kaluza-Klein reduction of these theories from five dimensions down to four, a double copy relation exists between the resulting massive four-dimensional scalar field theories. This requires the vanishing contribution of new galileon terms arising in high dimensions. We further explore if other interactions that do not arise from a dimensional reduction of the nonlinear sigma model could be double copied and find a new cubic interaction which satisfies the BCJ relations up to 5-point amplitudes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Linear and nonlinear benchmarks between the CLT code and the M3D-C1 code for the 2/1 resistive tearing mode and the 1/1 resistive kink mode

Here, the linear and nonlinear benchmarks between the CLT code and the M3D-C1 code for the 2/1 resistive tearing mode and the 1/1 resistive kink mode are presented. CLT is an explicit finite difference code, while M3D-C1 is an implicit finite element code. Although the implementations of CLT and M3D-C1 are totally different, we find that the simulation results of the resistive-kink mode and the m/n=2/1 tearing mode from M3D-C1 and CLT are almost the same, including the linear and nonlinear growth rates, the mode structures, the nonlinear saturation levels, the Poincare plots, and the scaling laws. This confirms that the nonlinear results for the 1/1 resistive-kink mode and 2/1 tearing mode are accurate and reliable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1991

Combined thermal and bending fatigue of high-temperature metal-matrix composites: Computational simulation

The nonlinear behavior of a high-temperature metal-matrix composite (HT-MMC) was simulated by using the metal matrix composite analyzer (METCAN) computer code. The simulation started with the fabrication process, proceeded to thermomechanical cyclic loading, and ended with the application of a monotonic load. Classical laminate theory and composite micromechanics and macromechanics are used in METCAN, along with a multifactor interaction model for the constituents behavior. The simulation of the stress-strain behavior from the macromechanical and the micromechanical points of view, as well as the initiation and final failure of the constituents and the plies in the composite, were examined in detail. It was shown that, when the fibers and the matrix were perfectly bonded, the fracture started in the matrix and then propagated with increasing load to the fibers. After the fibers fractured, the composite lost its capacity to carry additional load and fractured.

Gotsis, Pascal K.↗

NASA NTRS1992

Combined bending and thermal fatigue of high-temperature metal-matrix composites - Computational simulation

Gotsis, Pascal K.↗

NASA NTRS1985

Elastohydrodynamic lubrication calculations used as a tool to study scuffing

A new Reynolds equation was developed that takes into account the nonlinear viscous behavior of the fluid. The new Reynolds equation considers the nonlinear viscous fluid model of Eyring, the equilibrium equation, the constant mass flow, and the kinematic boundary condition. The new Reynolds equation and the elasticity equation are solved simultaneously by using a system approach and a Newton-Raphson technique. Comparisons are made with results obtained from the classical Reynolds equation. The effects of sliding speed and introducing a bump or a groove within the conjunction are studied. Results are shown for both moderate and heavy loads.

Houpert, L. G.↗

NASA NTRS1988

Helicopter mathematical models and control law development for handling qualities research

Progress made in joint NASA/Army research concerning rotorcraft flight-dynamics modeling, design methodologies for rotorcraft flight-control laws, and rotorcraft parameter identification is reviewed. Research into these interactive disciplines is needed to develop the analytical tools necessary to conduct flying qualities investigations using both the ground-based and in-flight simulators, and to permit an efficient means of performing flight test evaluation of rotorcraft flying qualities for specification compliance. The need for the research is particularly acute for rotorcraft because of their mathematical complexity, high order dynamic characteristics, and demanding mission requirements. The research in rotorcraft flight-dynamics modeling is pursued along two general directions: generic nonlinear models and nonlinear models for specific rotorcraft. In addition, linear models are generated that extend their utilization from 1-g flight to high-g maneuvers and expand their frequency range of validity for the design analysis of high-gain flight control systems. A variety of methods ranging from classical frequency-domain approaches to modern time-domain control methodology that are used in the design of rotorcraft flight control laws is reviewed. Also reviewed is a study conducted to investigate the design details associated with high-gain, digital flight control systems for combat rotorcraft. Parameter identification techniques developed for rotorcraft applications are reviewed.

Chen, Robert T. N.↗