Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear mechanics”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A transient FETI methodology for large-scale parallel implicit computations in structural mechanics, part 2

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallellize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet and perhaps will never be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient than explicit codes when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

Farhat, Charbel↗
DOE OSTI2021

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗
NASA NTRS1998

Analysis and Design of Variable Stiffness Composite Cylinders

An investigation of the possible performance improvements of thin circular cylindrical shells through the use of the variable stiffness concept is presented. The variable stiffness concept implies that the stiffness parameters change spatially throughout the structure. This situation is achieved mainly through the use of curvilinear fibers within a fiber-reinforced composite laminate, though the possibility of thickness variations and discrete stiffening elements is also allowed. These three mechanisms are incorporated into the constitutive laws for thin shells through the use of Classical Lamination Theory. The existence of stiffness variation within the structure warrants a formulation of the static equilibrium equations from the most basic principles. The governing equations include sufficient detail to correctly model several types of nonlinearity, including the formation of a nonlinear shell boundary layer as well as the Brazier effect due to nonlinear bending of long cylinders. Stress analysis and initial buckling estimates are formulated for a general variable stiffness cylinder. Results and comparisons for several simplifications of these highly complex governing equations are presented so that the ensuing numerical solutions are considered reliable and efficient enough for in-depth optimization studies. Four distinct cases of loading and stiffness variation are chosen to investigate possible areas of improvement that the variable stiffness concept may offer over traditional constant stiffness and/or stiffened structures. The initial investigation deals with the simplest solution for cylindrical shells in which all quantities are constant around the circumference of the cylinder. This axisymmetric case includes a stiffness variation exclusively in the axial direction, and the only pertinent loading scenarios include constant loads of axial compression, pressure, and torsion. The results for these cases indicate that little improvement over traditional laminates exists through the use of curvilinear fibers, mainly due to the presence of a weak link area within the stiffness variation that limits the ultimate load that the structure can withstand. Rigorous optimization studies reveal that even though slight increases in the critical loads can be produced for designs with an arbitrary variation of the fiber orientation angle, the improvements are not significant when compared to traditional design techniques that utilize ring stiffeners and frames. The second problem that is studied involves arbitrary loading of a cylinder with a stiffness variation that changes only in the circumferential direction. The end effects of the cylinder are ignored, so that the problem takes the form of an analysis of a cross-section for a short cylinder segment. Various load cases including axial compression, pressure, torsion, bending, and transverse shear forces are investigated. It is found that the most significant improvements in load-carrying capability exist for cases which involve loads that also vary around the circumference of the shell, namely bending and shear forces. The stiffness variation of the optimal designs contribute to the increased performance in two ways: lowering the stresses in the critical areas through redistribution of the stresses; and providing a relatively stiff region that alters the buckling behavior of the structure. These results lead to an in-depth optimization study involving weight optimization of a fuselage structure subjected to typical design constraints. Comparisons of the curvilinear fiber format to traditional stiffened structures constructed of isotropic and composite materials are included. It is found that standard variable stiffness designs are quite comparable in terms of weight and load-carrying capability yet offer the added advantage of tailorability of distinct regions of the structure that experience drastically different loading conditions. The last two problems presented in this work involve the nonlinear phenomenon of long tubes under bending. Though this scenario is not as applicable to fuselage structures as the previous problems, the mechanisms that produce the nonlinear effect are ideally suited to be controlled by the variable stiffness concept. This is due to the fact that the dominating influence for long cylinders under bending is the ovalization of the cross-section, which is governed mainly by the stiffness parameters of the cylindrical shell. Possible improvement of the critical buckling moments for these structures is investigated using either a circumferential or axial stiffness variation. For the circumferential case involving infinite length cylinders, it is found that slight improvements can be observed by designing structures that resist the cross-sectional deformation yet do not detract from the buckling resistance at the critical location. The results also indicate that buckling behavior is extremely dependent on cylinder length. This effect is most easily seen in the solution of finite length cylinders under bending that contain an axial stiffness variation. For these structures, the only mechanism that exhibits improved response are those that effectively shorten the length of the cylinder, thus reducing the cross-sectional deformation due to the forced restraint at the ends. It was found that the use of curvilinear fibers was not able to achieve this effect in sufficient degree to resist the deformation, but that ring stiffeners produced the desired response admirably. Thus, it is shown that the variable stiffness concept is most effective at improving the bending response of long cylinders through the use of a circumferential stiffness variation.

Tatting, Brian F.↗
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

Geometric distortions of autocollimation type scanners

One criterion of image quality is its geometric similarity to the control test. The ability of television, phototelevision, phototelegraphic, and other systems to transmit geometrically similar images is characterized by the geometric distortion coefficient. Linear distortions result from the inaccuracy of scanning mechanism element assembly, and distortions can be minimized during adjustment operations or may be easily compensated during image reception. Brief descriptions are given of the possible sources of nonlinear distortions which arise in these devices due to the scanning mechanism, objective distortion, or the nodding mirror.

Narayeva, M. K.↗
DOE OSTI2021

Post-compression of long-wave infrared 2 picosecond sub-terawatt pulses in bulk materials

We have experimentally demonstrated the post-compression of a long-wave infrared (9.2 μm) 150 GW peak power pulse from 2 ps to less than 500 fs using a sequence of two bulk materials with negative group velocity dispersion (GVD). The compression resulted in up to 1.6-fold increase of the peak power and up to 2.8-fold increase of the intensity in the center of a quasi-Gaussian beam. The partial decoupling of the self-phase modulation and chirp compensation stages by using two materials with significantly different ratios of nonlinear refractive index to GVD provides accurate optimization of the compression mechanism and promises a viable path to scaling peak powers to supra-terawatt levels. During the preparatory study, we measured, for the first time to our knowledge, the nonlinear refractive indices of NaCl, KCl, and BaF 2 for picosecond pulses in the long-wave infrared region.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Towards understanding the robust strengthening of ENSO and more frequent extreme El Niño events in CMIP6 global warming simulations

The El Niño Southern Oscillation (ENSO) has profound impacts on weather patterns across the globe, yet there is no consensus on its response to global warming. Several modelling studies find a stronger ENSO in global warming scenarios, while other studies suggest ENSO weakening. Using a broad range of models from the Coupled Model Intercomparison Project phase 6 (CMIP6) and four types of warming experiments, here we show that the majority of the models predict a stronger ENSO by century-end in Shared Social Pathway (SSP) experiments, and in idealized 1pctCO 2 and abrupt 4xCO 2 experiments. Several models, however, do predict no change or ENSO weakening, especially in the idealized experiments. Critically, the strongest forcing (abrupt-4xCO 2 ) does not induce the strongest ENSO response, while differences between the models are much greater than those between warming scenarios. For the long-term response (over 1000 years) the models disagree even on the sign of change. Furthermore, changes in ENSO sea surface temperature (SST) variability are only modestly correlated with the tropical Pacific mean state change. The highest correlation for ENSO SST amplitude is found with the mean zonal SST gradient in the SSP5-8.5 experiment (R = –0.58). In contrast, changes in ENSO rainfall variability correlate well with changes in the mean state, as well as with changes in ENSO SST variability. When evaluating the Bjerknes Stability Index for a subset of models, we find that it is not a reliable predictor of ENSO strengthening, as this index tends to predict greater stability with warming. We argue that the enhanced ENSO stability is offset by increases in atmospheric noise or/and potential nonlinear effects. However, a robust inter-model mechanism that can explain a stronger ENSO simulated with global warming is still lacking. Therefore, caution should be exercised when considering ENSO changes based on a single model or warming scenario.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Optical Bound States in the Continuum Enabled by Magnetic Resonances Coupled to a Mirror

Dielectric metasurfaces made of high refractive index and low optical loss materials have emerged as promising platforms to achieve high-quality factor modes enabling strong light–matter interaction. Bound states in the continuum have shown potential to demonstrate narrow spectral resonances but often require asymmetric geometry and typically feature strong polarization dependence, complicating fabrication and limiting practical applications. Here, we introduce a novel approach for designing high-quality bound states in the continuum using magnetic dipole resonances coupled to a mirror. The resulting metasurface has simple geometric parameters requiring no broken symmetry. To demonstrate the unique features of our photonic platform we show a record-breaking third harmonic generation efficiency from the metasurface benefiting from the strongly enhanced electric field at high-quality resonances. Our approach mitigates the shortcomings of previous platforms with simple geometry enabling facile and large-area fabrication of metasurfaces paving the way for applications in optical sensing, detection, quantum photonics, and nonlinear devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Temperature-Dependent Properties of Molten Li 2 BeF 4 Salt Using Ab Initio Molecular Dynamics

Molten lithium tetrafluoroberyllate (Li 2 BeF 4 ) salt, also known as FLiBe, with a 2:1 mixture of LiF and BeF 2 is being proposed as a coolant and solvent in advanced nuclear reactor designs, such as the molten salt reactor or the fluoride salt cooled high-temperature reactor. We present the results on the structure and properties of FLiBe over a wide range of temperatures, 0–2000 K, from high-throughput ab initio molecular dynamics simulation using a supercell model of 504 atoms. The variations in the local structures of solid and liquid FLiBe with temperature are discussed in terms of a pair distribution function, coordination number, and bond angle distribution. The temperature-dependent electronic structure and optical and mechanical properties of FLiBe are calculated. The optical and mechanical property results are reported for the first time. The results above and below the melting temperature (~732 K) are compared with the experimental data and with data for crystalline FLiBe. The electronic structure and interatomic bonding results are discussed in correlation with the mechanical strength. A novel concept of total bond order density (TBOD), an important quantum mechanical parameter, is used to characterize the internal cohesion and strength in the simulated models. The results show a variation in the rate of change in properties in solid and liquid phases with anomalous behavior across the melting region. The observed trend is the decrease in mechanical strength, band gap, and TBOD in a nonlinear fashion as a function of temperature. The refractive index shows a surprising minimum at 850 K, among the tested temperatures, which lies above the melting point. These findings provide a new platform to understand the interplay between the temperature-dependent structures and properties of FLiBe salt.

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

Joint Experimental and Computational Characterization of Sum-Frequency Generation between a Continuous Wave Laser and an Ultrafast Frequency Comb Laser for Tunable Laser Development

Ultrafast optical frequency combs allow for both high spectral and temporal resolution in molecular spectroscopy and have become a powerful tool in many areas of chemistry and physics. Ultrafast lasers and frequency combs generated from ultrafast mode-locked lasers often need to be converted to other wavelengths. Commonly used wavelength conversions are optical parametric oscillators, which require an external optical cavity, and supercontinuum generation combined with optical parametric amplifiers. Whether commercial or home-built, these systems are complex and costly. Here, we investigate an alternative, simple, and easy-to-implement approach to tunable frequency comb ultrafast lasers enabled by new continuous-wave laser technology. Sum-frequency generation between an Nd:YAG continuous-wave laser and a Yb:fiber femtosecond frequency comb in a beta-barium borate (BBO) crystal is explored. The resulting sum-frequency beam is a pulsed frequency comb with the same repetition rate as the Yb:fiber source. SNLO simulation software is used to simulate the results and provide benchmarks for designing future systems to achieve wavelength conversion and tunability in otherwise difficult-to-reach spectral regions.

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

A database of refractive indices and dielectric constants auto-generated using ChemDataExtractor

The ability to auto-generate databases of optical properties holds great potential for advancing optical research, especially with regards to the data-driven discovery of optical materials. An optical property database of refractive indices and dielectric constants is presented, which comprises a total of 49,076 refractive index and 60,804 dielectric constant data records on 11,054 unique chemicals. The database was auto-generated using the state-of-the-art natural language processing software, ChemDataExtractor, using a corpus of 388,461 scientific papers. The data repository offers a representative overview of the information on linear optical properties that resides in scientific papers from the past 30 years. Public availability of these data will enable a quick search for the optical property of certain materials. The large size of this repository will accelerate data-driven research on the design and prediction of optical materials and their properties. To the best of our knowledge, this is the first auto-generated database of optical properties from a large number of scientific papers. We provide a web interface to aid the use of our database.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Magnetic field effects and transverse ratchets in charge lattices coupled to asymmetric substrates

We examine a charge lattice coupled to a one-dimensional asymmetric potential in the presence of an applied magnetic field, which induces gyrotropic effects in the charge motion. This system could be realized for Wigner crystals in nanostructured samples, dusty plasmas, or other classical charge-ordered states where gyrotropic motion and damping can arise. For zero magnetic field, an applied external ac drive can produce a ratchet effect in which the particles move along the easy flow direction of the substrate asymmetry. The zero field ratchet effect can only occur when the ac drive is aligned with the substrate asymmetry direction; however, when a magnetic field is added, the gyrotropic forces generate a Hall effect that leads to a variety of new behaviors, including a transverse ratchet motion that occurs when the ac drive is perpendicular to the substrate asymmetry direction. We show that this system exhibits commensuration effects as well as reversals in the ratchet effect and the Hall angle of the motion. The magnetic field also produces a nonmonotonic ratchet efficiency when the particles become localized at high fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Generating and processing optical waveforms using spectral singularities

Here, we show that a laser at threshold can be utilized to generate the class of coherent and transform-limited waveforms (vt — z) m e i(kz—ωt) at optical frequencies. We derive these properties analytically and demonstrate them in semiclassical time-domain laser simulations. We then utilize these waveforms to expand other waveforms with high modulation frequencies and demonstrate theoretically the feasibility of complex-frequency coherent absorption at optical frequencies, with efficient energy transduction and cavity loading. This approach has potential applications in quantum computing, photonic circuits, and biomedicine.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Semidefinite programming algorithm for the quantum mechanical bootstrap

Here, we present a semidefinite programming algorithm to find eigenvalues of Schrödinger operators within the bootstrap approach to quantum mechanics. The bootstrap approach involves two ingredients: a nonlinear set of constraints on the variables (expectation values of operators in an energy eigenstate), plus positivity constraints (unitarity) that need to be satisfied. By fixing the energy we linearize all the constraints and show that the feasibility problem can be presented as an optimization problem for the variables that are not fixed by the constraints and one additional slack variable that measures the failure of positivity. To illustrate the method we are able to obtain high-precision, sharp bounds on eigenenergies for arbitrary confining polynomial potentials in one dimension.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

A three-dimensional model of dynamical processes in the Venus atmosphere

Three-dimensional calculations of the circulation of the Venus atmosphere have resulted in mean zonal winds in the same direction and of the same magnitude as those observed; i.e., retrograde with speeds of the order of 100 m/s. The solutions exhibit other observed properties of the circulation: small horizontal temperature contrasts with the larger variations being between equator and pole, meridional velocities at mid and low latitudes less than 10 m/s, and the existence of planetary waves which at certain times show vertical flow-field contours in a horizontal Y configuration. The mechanism maintaining the large zonal winds is a nonlinear instability involving both the mean meridional circulation and planetary-scale eddies. The meridional circulation is the principal means by which zonal momentum is transported vertically. Planetary-scale eddies are the principal means by which potential energy is released, and they are also significant in transporting angular momentum horizontally. Planetary rotation plays an important role in initially generating the mean zonal winds starting from rest. Initial conditions affect the characteristics of the solutions, including the magnitude of the mean zonal velocity and whether or not planetary waves are generated.

Young, R. E.↗
NASA NTRS1980

Proposed method for determining the thickness of glass in solar collector panels

An analytical method was developed for determining the minimum thickness for simply supported, rectangular glass plates subjected to uniform normal pressure environmental loads such as wind, earthquake, snow, and deadweight. The method consists of comparing an analytical prediction of the stress in the glass panel to a glass breakage stress determined from fracture mechanics considerations. Based on extensive analysis using the nonlinear finite element structural analysis program ARGUS, design curves for the structural analysis of simply supported rectangular plates were developed. These curves yield the center deflection, center stress and corner stress as a function of a dimensionless parameter describing the load intensity. A method of estimating the glass breakage stress as a function of a specified failure rate, degree of glass temper, design life, load duration time, and panel size is also presented.

Moore, D. M.↗