Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear physics”

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 127 records · Page 7

Overview of Sensitivity Analysis and Shape Optimization for Complex Aerodynamic Configurations

This paper presents a brief overview of some of the more recent advances in steady aerodynamic shape-design sensitivity analysis and optimization, based on advanced computational fluid dynamics (CFD). The focus here is on those methods particularly well-suited to the study of geometrically complex configurations and their potentially complex associated flow physics. When nonlinear state equations are considered in the optimization process, difficulties are found in the application of sensitivity analysis. Some techniques for circumventing such difficulties are currently being explored and are included here. Attention is directed to methods that utilize automatic differentiation to obtain aerodynamic sensitivity derivatives for both complex configurations and complex flow physics. Various examples of shape-design sensitivity analysis for unstructured-grid CFD algorithms are demonstrated for different formulations of the sensitivity equations. Finally, the use of advanced, unstructured-grid CFDs in multidisciplinary analyses and multidisciplinary sensitivity analyses within future optimization processes is recommended and encouraged.

Newman, James C., III↗

Computer implemented empirical mode decomposition method apparatus, and article of manufacture utilizing curvature extrema

A computer implemented physical signal analysis method is includes two essential steps and the associated presentation techniques of the results. All the steps exist only in a computer: there are no analytic expressions resulting from the method. The first step is a computer implemented Empirical Mode Decomposition to extract a collection of Intrinsic Mode Functions (IMF) from nonlinear, nonstationary physical signals based on local extrema and curvature extrema. The decomposition is based on the direct extraction of the energy associated with various intrinsic time scales in the physical signal. Expressed in the IMF's, they have well-behaved Hilbert Transforms from which instantaneous frequencies can be calculated. The second step is the Hilbert Transform. The final result is the Hilbert Spectrum. Thus, the invention can localize any event on the time as well as the frequency axis. The decomposition can also be viewed as an expansion of the data in terms of the IMF's. Then, these IMF's, based on and derived from the data, can serve as the basis of that expansion. The local energy and the instantaneous frequency derived from the IMF's through the Hilbert transform give a full energy-frequency-time distribution of the data which is designated as the Hilbert Spectrum.

Huang, Norden Eh↗

Computer implemented empirical mode decomposition method, apparatus and article of manufacture

A computer implemented physical signal analysis method is invented. This method includes two essential steps and the associated presentation techniques of the results. All the steps exist only in a computer: there are no analytic expressions resulting from the method. The first step is a computer implemented Empirical Mode Decomposition to extract a collection of Intrinsic Mode Functions (IMF) from nonlinear, nonstationary physical signals. The decomposition is based on the direct extraction of the energy associated with various intrinsic time scales in the physical signal. Expressed in the IMF's, they have well-behaved Hilbert Transforms from which instantaneous frequencies can be calculated. The second step is the Hilbert Transform. The final result is the Hilbert Spectrum. Thus, the invention can localize any event on the time as well as the frequency axis. The decomposition can also be viewed as an expansion of the data in terms of the IMF's. Then, these IMF's, based on and derived from the data, can serve as the basis of that expansion. The local energy and the instantaneous frequency derived from the IMF's through the Hilbert transform give a full energy-frequency-time distribution of the data which is designated as the Hilbert Spectrum.

Huang, Norden E.↗

Gyrokinetic turbulence cascade via predator-prey interactions between different scales

Gyrokinetic simulations in a closed fieldline geometry are presented to explore the physics of nonlinear transfer in plasma turbulence. As spontaneously formed zonal flows and small-scale turbulence demonstrate “predator-prey” dynamics, a particular cascade spectrum emerges. The electrostatic potential and the density spectra appear to be in good agreement with the simple theoretical prediction based on Charney-Hasegawa-Mima equation | ϕ{sup ~}{sub k} |{sup 2}∼| n{sup ~}{sub k} |{sup 2}∝k{sup −3}/(1+k{sup 2}){sup 2}, with the spectra becoming anisotropic at small scales. The results indicate that the disparate scale interactions, in particular, the refraction and shearing of larger scale eddies by the self-consistent zonal flows, dominate over local interactions, and contrary to the common wisdom, the comprehensive scaling relation is created even within the energy injection region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Higher-order Dirac solitons in binary waveguide arrays

We study optical analogues of higher-order Dirac solitons (HODSs) in binary waveguide arrays. Like higher-order solitons obtained from the well-known nonlinear Schrödinger equation governing the pulse propagation in an optical fiber, these HODSs have amplitude profiles which are numerically shown to be periodic over large propagation distances. At the same time, HODSs possess some unique features. Firstly, the period of a HODS depends on its order parameter. Secondly, the discrete nature in binary waveguide arrays imposes the upper limit on the order parameter of HODSs. Thirdly, the order parameter of HODSs can vary continuously in a certain range. - Highlights: • Higher-order Dirac solitons in nonlinear binary waveguide arrays are numerically demonstrated. • Amplitude profiles of higher-order Dirac solitons are periodic during propagation. • The period of higher-order Dirac solitons decreases when the soliton order increases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Non-scalar uncertainty: Uncertainty in dynamic systems

The following point is stated throughout the paper: dynamic systems are usually subject to uncertainty, be it the unavoidable quantic uncertainty when working with sufficiently small scales or when working in large scales uncertainty can be allowed by the researcher in order to simplify the problem, or it can be introduced by nonlinear interactions. Even though non-quantic uncertainty can generally be dealt with by using the ordinary probability formalisms, it can also be studied with the proposed non-scalar formalism. Thus, non-scalar uncertainty is a more general theoretical framework giving insight into the nature of uncertainty and providing a practical tool in those cases in which scalar uncertainty is not enough, such as when studying highly nonlinear dynamic systems. This paper's specific contribution is the general concept of non-scalar uncertainty and a first proposal for a methodology. Applications should be based upon this methodology. The advantage of this approach is to provide simpler mathematical models for prediction of the system states. Present conventional tools for dealing with uncertainty prove insufficient for an effective description of some dynamic systems. The main limitations are overcome abandoning ordinary scalar algebra in the real interval (0, 1) in favor of a tensor field with a much richer structure and generality. This approach gives insight into the interpretation of Quantum Mechanics and will have its most profound consequences in the fields of elementary particle physics and nonlinear dynamic systems. Concepts like 'interfering alternatives' and 'discrete states' have an elegant explanation in this framework in terms of properties of dynamic systems such as strange attractors and chaos. The tensor formalism proves especially useful to describe the mechanics of representing dynamic systems with models that are closer to reality and have relatively much simpler solutions. It was found to be wise to get an approximate solution to an accurate model than to get a precise solution to a model constrained by simplifying assumptions. Precision has a very heavy cost in present physical models, but this formalism allows the trade between uncertainty and simplicity. It was found that modeling reality sometimes requires that state transition probabilities should be manipulated as nonscalar quantities, finding at the end that there is always a transformation to get back to scalar probability.

Martinez, Salvador Gutierrez↗

Phonon-enhanced nonlinearities in hexagonal boron nitride

Abstract Polar crystals can be driven into collective oscillations by optical fields tuned to precise resonance frequencies. As the amplitude of the excited phonon modes increases, novel processes scaling non-linearly with the applied fields begin to contribute to the dynamics of the atomic system. Here we show two such optical nonlinearities that are induced and enhanced by the strong phonon resonance in the van der Waals crystal hexagonal boron nitride (hBN). We predict and observe large sub-picosecond duration signals due to four-wave mixing (FWM) during resonant excitation. The resulting FWM signal allows for time-resolved observation of the crystal motion. In addition, we observe enhancements of third-harmonic generation with resonant pumping at the hBN transverse optical phonon. Phonon-induced nonlinear enhancements are also predicted to yield large increases in high-harmonic efficiencies beyond the third.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Developing nonlinear laser-plasma instability models for high-fidelity, multi-physics simulation capability for ICF/HED (IC Report) [Slides]

Modeling nonlinear effects of laser plasma instabilities (LPI) is critically important for ICF/HED experiments. Our project has the following goals: • Develop physics-based nonlinear LPI models using PIC simulations • Couple the nonlinear LPI effects to macroscopic modeling of ICF/HED exp. through laser ray-tracing (LRT).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Raman-induced mode coupling in temporal waveguides formed by short solitons

We study the propagation of optical pulses trapped inside a temporal waveguide formed by two solitons in a dispersive nonlinear medium such as an optical fiber. The solitons are short enough that they decelerate as their spectra shift continuously toward the red side because of intrapulse Raman scattering. We show that the shape of a probe pulse trapped between the two solitons evolves in a periodic fashion, while its spectrum shifts toward the blue side. We develop a coupled-mode theory showing that such changes occur because of mode coupling induced by the deceleration of the short solitons, resulting in a curved waveguide. Here, a simplified two-mode model is used to introduce a single-parameter governing modal coupling and to find the condition under which coupling becomes weak enough that the probe pulse blue-shifts its spectrum without changes in its pulse shape.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulating non-native cubic interactions on noisy quantum machines

As a milestone for general-purpose computing machines, we demonstrate that quantum processors can be programed to efficiently simulate dynamics that are not native to the hardware. Moreover, on noisy devices without error correction, we show that simulation results are significantly improved when the quantum program is compiled using modular gates instead of a restricted set of standard gates. We demonstrate the general methodology by solving a cubic interaction problem, which appears in nonlinear optics, gauge theories, as well as plasma and fluid dynamics. To encode the non-native Hamiltonian evolution, we decompose the Hilbert space into a direct sum of invariant subspaces in which the nonlinear problem is mapped to a finite-dimensional Hamiltonian simulation problem. Furthermore, in a three-states example, the resultant unitary evolution is realized by a product of approximately 20 standard gates, using which approximately ten simulation steps can be carried out on state-of-the-art quantum hardware before results are corrupted by decoherence. In comparison, the simulation depth is improved by more than an order of magnitude when the unitary evolution is realized as a single cubic gate, which is compiled directly using optimal control. Alternatively, parametric gates may also be compiled by interpolating control pulses. Modular gates thus obtained provide high-fidelity building blocks for quantum Hamiltonian simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

State preparation and measurement in a quantum simulation of the $O(\mathrm{3})$ sigma model

Recently, Singh and Chandrasekharan [Phys. Rev. D 100, 054505 (2019)] showed that fixed points of the nonlinear O ( 3 ) sigma model can be reproduced near a quantum phase transition of a spin model with just two qubits per lattice site. In a paper by the NuQS Collaboration [Phys. Rev. Lett. 123, 090501 (2019)], the proposal is made to simulate such field theories on a quantum computer using the universal properties of a similar model. In this paper, following that direction, we demonstrate how to prepare the ground state of the model from Singh and Chandrasekharan and measure a dynamical quantity of interest, the O ( 3 ) Noether charge, on a quantum computer. In particular, we apply Trotter methods to obtain results for the complexity of adiabatic ground state preparation in both the weak-coupling and quantum-critical regimes and use shadow tomography to measure the dynamics of local observables. Finally, we then present and analyze a quantum algorithm based on nonunitary randomized simulation methods that may yield an approach suitable for intermediate-term noisy quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ultrafast nonlinear absorption of TMDC quantum dots

Here, we study theoretically absorption properties of MoS 2 quantum dots placed in an ultrashort and strong optical pulse with the duration of a few femtoseconds and the amplitude up to 1.5 V/Å. While the monolayer of MoS 2 has a bandgap of ≈ 1.8 eV, the quantum dots of MoS 2 have in-gap edge states, which open additional absorption channels. The absorbance of the system strongly depends on the amplitude of the pulse. At small amplitude, less than 0.5 V/Å, the absorbance varies by more than 20% with the frequency of the pulse. The system shows saturable absorption if the frequency of the pulse is larger than the bulk bandgap and reverse saturable absorption if the frequency of the pulse is in the bandgap. At large amplitudes, > 0.5 V/Å, of the pulse, the ultrafast absorbance of a quantum dot has a weak dependence on the pulse frequency. If the frequency of the pulse is less than the corresponding bulk bandgap, the absorbance has a nonmonotonic dependence on the pulse amplitude with the maximum at ≈ 0.5 V/Å. We also present the results for QDs of different TMDC materials, which show that at large field amplitude, > 0.5 V/Å, the absorbance is universal and has weak dependence on the type of TMDC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Plasma-based beam combiner for very high fluence and energy

Extreme optical fluences, much beyond the damage threshold of conventional optics, are of interest for a range of high-energy-density physics applications. Nonlinear interactions of multiple beams in plasmas have the potential to produce optics that operate at much higher intensity and fluence than is possible in solids. In inertial confinement fusion experiments indirectly driven with lasers, many beams overlap in the plasma inside a hohlraum, and cross-beam energy transfer by Brillouin scattering has been employed to redistribute energy between laser beams within the target. Here in this paper, we show that in a hot, under-dense plasma the energy of many input beams can be combined into a single well-collimated beam. The emerging beam has an energy of 4 kJ (over 1 ns) that is more than triple that of any incident beam, and a fluence that is more than double. Because the optic produced is plasma, and is diffractive, it is inherently capable of generating higher fluences in a single beam than solid-state refractive or reflective optics.

42 ENGINEERING↗

Nonlinear Diffusion Acceleration of the Least-Squares Transport Equation in Geometries with Voids

In this paper we show the extension of the Nonlinear-Diffusion Acceleration (NDA) to geometries containing small voids using a weighted least-squares (WLS) high order equation. Even though the WLS equation is well defined in voids, the low-order drift diffusion equation was not defined in materials with a zero cross section. This paper derives the necessary modifications to the NDA algorithm. We show that a small change to the NDA closure term and a non-local definition of the diffusion coefficient solve the problems for voids regions. These changes do not affect the algorithm for optical thick material regions, while making the algorithm well defined in optically thin ones. We use a Fourier analysis to perform an iterative analysis to confirm that the modifications result in a stable and efficient algorithm. Numerical results of our method will be presented in the second part of the paper. We test this formulation with a small, one-dimensional test problem. Additionally we present results for a modified version of the C5G7 benchmark containing voids as a more complex, reactor like problem. We compared our results to PDT, Texas A&M’s transport code, utilizing a first order discontinuous formulation as reference and the self-adjoint angular flux equation with void treatment (SAAFτ ), a different second order form. The results indicate, that the NDA WLS performed comparably or slightly worse then the asymmetric SAAFτ , while maintaining a symmetric discretization matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

When physics-informed data analytics outperforms black-box machine learning: A case study in thickness control for additive manufacturing

Aerosol jet printing (AJP) has emerged as a promising noncontact additive manufacturing method for high-resolution printing for a wide range of material systems. A key challenge limiting the broader adoption of AJP in the material science community is the lack of methods to precisely control thickness. Herein, we develop a model-based design of experiment (MBDoE) framework that integrates physics-informed models, nonlinear regression, and information criteria to postulate, select and calibrate the best model to describe and optimize the AJP manufacturing process. Starting with already available data from system commissioning (e.g., prior single variable sensitivity analysis), four candidate physics-informed models are postulated and trained. MBDoE identifies a single additional optimal experiment to validate these predictive models with quantified uncertainties, which are then used to determine the best experimental conditions to control printed film thickness. As a comparative benchmark, the analysis is repeated using the same dataset with nonparametric Gaussian process regression (GPR) model that does not incorporate physical information. Using MBDoE principles, we find that only five experiments are necessary to calibrate the nonlinear physics-informed parametric model, and with said limited data, this model outperforms the black-box machine learning GPR model. This key result underscores an emerging trend in the data science community: incorporating physical information into predictive models often drastically reduces the data requirements. Leveraging MBDoE further increased the data efficiency. By design, the proposed data science framework is general in nature and can be easily extended to other experimental and additive manufacturing systems beyond AJP.

Aerosol jet printing↗

Squeeze Film Damper Experimental and Numerical Correlation: Test Setup Description and Parameter Identification of Dry System

The process of characterizing a damper involves both experimental techniques and model simulation. Experiments characterize dynamic behavior of a system through empirical modal analysis using impact and shaker testing. The layout of the test setup is crucial for a low-mass, nonlinear system. Physical arrangement, orientation of the system, and instrumentation affect the results that are discussed in this paper. These results can be used to validate a numerical simulation of the same system. With a validated model, important physical input parameters and modeling techniques can be determined and implemented in future designs. In other words, this modeling can be used with confidence to design alternative dampers that may be difficult, costly, and time-consuming to characterize experimentally.

Cook, Jason↗

The Most Effective Model for Describing the Universal Behavior of a Noisy Kuramoto–Sivashinsky Equation as a Paradigmatic Model

We study a noisy Kuramoto–Sivashinsky (KS) equation which describes unstable surface growth and chemical turbulence. It has been conjectured that the universal long-wavelength behavior of the equation, which is characterized by scale-dependent parameters, is described by a Kardar–Parisi–Zhang (KPZ) equation. We consider this conjecture by analyzing a renormalization-group equation for a class of generalized KPZ equations. We then uniquely determine the parameter values of the KPZ equation that most effectively describes the universal long-wavelength behavior of the noisy KS equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗