Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Optical nonlinear properties and dynamics of interband transitions in multilayer MoS{sub 2} structures under femtosecond excitation at a wavelength of 514 nm

The optical nonlinear absorption and bleaching of aqueous suspensions of multilayer MoS{sub 2} sheets (structural modification 2H) under excitation by a 400-fs pulse at a wavelength of 514 nm is investigated using longitudinal scanning. The sample exhibits nonlinear absorption at intensities up to 15 GW cm{sup −2}, while a further increase in intensity to 70 GW cm{sup −2} causes nonlinear bleaching with a relative change in transmission to 14 %. The dynamics of interband transitions in the picosecond range is studied by femtosecond laser photolysis. The relaxation time of photoexcited excitons is measured to be 20 ± 2 ps. The transition dynamics is calculated in the three-level approximation, and the absorption cross sections of photoinduced electron transitions from the valence band to the conduction band and from the first to the second conduction band are estimated. It is shown that the optical nonlinear properties of suspensions of multilayer 2H MoS{sub 2} sheets are mainly determined by the dynamics of single-photon interband transitions. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A boundary condition to the Khokhlov-Zabolotskaya equation for modeling strongly focused nonlinear ultrasound fields

An equivalent source model was proposed as a boundary condition to the nonlinear parabolic Khokhlov-Zabolotskaya (KZ) equation to simulate high intensity focused ultrasound (HIFU) fields generated by medical ultrasound transducers with the shape of a spherical shell. The boundary condition was set in the initial plane; the aperture, the focal distance, and the initial pressure of the source were chosen based on the best match of the axial pressure amplitude and phase distributions in the Rayleigh integral analytic solution for a spherical transducer and the linear parabolic approximation solution for the equivalent source. Analytic expressions for the equivalent source parameters were derived. It was shown that the proposed approach allowed us to transfer the boundary condition from the spherical surface to the plane and to achieve a very good match between the linear field solutions of the parabolic and full diffraction models even for highly focused sources with F-number less than unity. The proposed method can be further used to expand the capabilities of the KZ nonlinear parabolic equation for efficient modeling of HIFU fields generated by strongly focused sources.

42 ENGINEERING↗
DOE OSTI2025

Experimental observation of nonlinear relation between pressure and water flux is consistent with the solution-diffusion model

In several recent studies, it has been proposed that the fundamental understanding of penetrant transport in dense polymer membranes occurring via the solution-diffusion model, which has been the generally accepted theoretical framework for describing penetrant transport in such materials for the past several decades, is flawed. An alternate mechanistic framework based on the idea of two-phase flow in a porous medium (i.e., pore-flow) has been broadly advanced instead, with proponents of this approach claiming that the pore-flow theoretical framework provides the necessary mechanistic insight to design novel polymeric membrane materials for emerging applications. In this study, we show experimental results for hydraulic permeation of water that are entirely consistent with the solution-diffusion theory, without modification, for three dense polymeric membranes: crosslinked poly(ethylene glycol diacrylate) (XLPEGDA), Nafion 117 ionomer in the sodium counterion form (Nafion 117-Na), and cellulose acetate (CA). By measuring water flux at transmembrane pressures up to 240 bar, we observe a nonlinear relationship between the transmembrane pressure (TMP) and water flux, J w , for XLPEGDA and Nafion 117-Na, while this relationship is linear for CA. We demonstrate that the behavior of these three materials is described via the solution-diffusion model. According to the solution-diffusion model, flux is, to a good approximation, proportional to the transmembrane concentration difference induced by the pressure difference across the membrane, rather than to TMP itself. Water sorption isotherms are reported for all three materials. They further justify the nonlinear relationship between TMP and J w observed in XLPEGDA and Nafion 117-Na, emphasizing that the nonlinearity in the flux/TMP relationship stems from nonlinearities in the sorption isotherm with pressure. Additionally, the relationship between water flux and TMP can be predicted, a priori, with no adjustable parameters when a predictive model for the diffusion coefficient of water is employed in conjunction with the experimental water sorption isotherms in the solution-diffusion model. Furthermore, our results demonstrate the validity of the solution-diffusion model to describe transport of penetrants in dense polymer membranes, while highlighting the sensitivity of the solution-diffusion model to the many physical and mathematical simplifications commonly applied to the theory in literature.

materials↗
DOE OSTI2023

Control of Fractional Diffusion Problems via Dynamic Programming Equations

In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton–Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. Here, the formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling.

42 ENGINEERING↗
DOE OSTI2024

Optimizing spectral phase transfer in four-wave mixing with gas-filled capillaries

Four-wave mixing (FWM) in gas-filled hollow-core capillaries, a nonlinear optical process that mixes signal and pump photon frequencies to generate idler frequency photons, offers a method for precise spectral phase transfer from signal to idler at ultrashort timescales and extreme powers. However, this regime is challenged by competing linear and nonlinear dynamics, leading to significant trade-offs between spectral phase transfer and conversion efficiency. Our computational investigation focuses on the upconversion of femtosecond pulses from the infrared (IR) to the ultraviolet (UV), a range notoriously difficult to manipulate. We explore an intermediate energy regime that strikes an optimal balance between FWM-mediated phase-transfer fidelity and nonlinear conversion efficiency. By adjusting the energy ratios and spectral phase profiles of the input signal, we achieve conversion efficiencies of approximately 5-15% while maintaining an effective quasi-linear spectral phase transfer. These findings will contribute to establishing first-principles and scaling laws essential for applications such as high-precision imaging, spectroscopy, quantum transduction, and distributed entangled interconnects, facilitating advanced control of ultrafast photonic and electronic wavepackets in quantum materials with unprecedented spatial and temporal precision.

Zhang, Hao↗
DOE OSTI2022

Paraxial, Thin-Lens Analysis of Fixed-Tune, Non-Scaling FFAs with Two Magnets per Cell

Fixed-field accelerators (FFAs) have closed orbits that change as a function of beam momentum. It is sometimes useful to avoid various resonances by keeping the tunes constant as these orbits change. A well-known example is the scaling FFA where the entire beam orbit and optics are geometrically scaled as a function of momentum. However, this is stricter than necessary: cells with three or more lenses can have fixed tunes even when the focussing strengths in the lenses change. Recently, Dejan Trbojevic has found a pair of nonlinear magnets that produce fixed tunes and fixed beta functions while not being a scaling FFA. This improves on a more approximate solution in. Notably a scaling FFA would require one magnet to be entirely reverse-bending, whereas the Trbojevic solution does not do this and instead resembles an intermediate point between the traditional scaling field profiles and the nonscaling profiles centred around a momentum with equal and positive fields. This suggests there are at least three levels of stringency that one can apply to a fixed-tune FFA design: 1. Fixed cell tunes (in both planes) as a function of momentum; 2. Fixed optics (beta functions) as a function of momentum; 3. Similarity of all orbits via a scaling symmetry law ↔ traditional scaling FFA. This note studies the interesting case #2 above (#3 being fully characterised by the orbit at a single energy) in the simplest possible example: a cell of two thin lenses in the small angle (paraxial) approximation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Assessing the Optimality of LinDist3Flow for Optimal Tap Selection of Step Voltage Regulators in Unbalanced Distribution Networks

The adoption of distributed energy resources such as photovoltaics (PVs) has increased dramatically during the previous decade. The increased penetration of PVs into distribution networks (DNs) can cause voltage fluctuations that have to be mitigated. One of the key utility assets employed to this end are step-voltage regulators (SVRs). It is desirable to include tap selection of SVRs in optimal power flow (OPF) routines, a task that turns out to be challenging because the resultant OPF problem is nonconvex with added complexities stemming from accurate SVR modeling. While several convex relaxations based on semi-definite programming (SDP) have been presented in the literature for optimal tap selection, SDP based schemes do not scale well and are challenging to implement in large-scale planning or operational frameworks. This paper deals with the optimal tap selection (OPTS) problem for wye-connected SVRs using linear approximations of power flow equations. Specifically, the LinDist3Flow model is adopted and the effective SVR ratio is assumed to be continuous- enabling the formulation of a problem called LinDist3Flow-OPTS, which amounts to a linear program. The scalability and optimality gap of LinDist3Flow-OPTS are evaluated with respect to existing SDP-based and nonlinear programming techniques for optimal tap selection in three standard feeders, namely, the IEEE 13-bus, 123-bus, and 8500-node DNs. For all DNs considered, LinDist3Flow-OPTS achieves an optimality gap of approximately 1% or less while significantly lowering the computational burden.

linear approximations↗
DOE OSTI2022

Assessing the Optimality of LinDist3Flow for Optimal Tap Selection of Step Voltage Regulators in Unbalanced Distribution Networks: Preprint

The adoption of distributed energy resources such as photovoltaics (PVs) has increased dramatically during the previous decade. The increased penetration of PVs into distribution networks (DNs) can cause voltage fluctuations that have to be mitigated. One of the key utility assets employed to this end are step-voltage regulators (SVRs). It is desirable to include tap selection of SVRs in optimal power flow (OPF) routines, a task that turns out to be challenging because the resultant OPF problem is nonconvex with added complexities stemming from accurate SVR modeling. While several convex relaxations based on semi-definite programming (SDP) have been presented in the literature for optimal tap selection, SDP based schemes do not scale well and are challenging to implement in large-scale planning or operational frameworks. This paper deals with the optimal tap selection (OPTS) problem for wye-connected SVRs using linear approximations of power flow equations. Specifically, the LinDist3Flow model is adopted and the effective SVR ratio is assumed to be continuous–enabling the formulation of a problem called LinDist3Flow-OPTS, which amounts to a linear program. The scalability and optimality gap of LinDist3Flow-OPTS are evaluated with respect to existing SDP-based and nonlinear programming techniques for optimal tap selection in three standard feeders, namely, the IEEE 13-bus, 123-bus, and 8500-node DNs. For all DNs considered, LinDist3Flow-OPTS achieves an optimality gap of approximately 1% or less while significantly lowering the computational burden.

linear approximations↗
DOE OSTI2021

Nonlinear proper orthogonal decomposition for convection-dominated flows

Autoencoder techniques find increasingly common use in reduced order modeling as a means to create a latent space. This reduced order representation offers a modular data-driven modeling approach for nonlinear dynamical systems when integrated with a time series predictive model. In this Letter, we put forth a nonlinear proper orthogonal decomposition (POD) framework, which is an end-to-end Galerkin-free model combining autoencoders with long short-term memory networks for dynamics. By eliminating the projection error due to the truncation of Galerkin models, a key enabler of the proposed nonintrusive approach is the kinematic construction of a nonlinear mapping between the full-rank expansion of the POD coefficients and the latent space where the dynamics evolve. We test our framework for model reduction of a convection-dominated system, which is generally challenging for reduced order models. Our approach not only improves the accuracy, but also significantly reduces the computational cost of training and testing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Hyper-Rayleigh scattering in centrosymmetric systems

Hyper-Rayleigh scattering (HRS) is an incoherent mechanism for optical second harmonic generation. The frequency-doubled light that emerges from this mechanism is not emitted in a laser-like manner, in the forward direction; it is scattered in all directions. The underlying theory for this effect involves terms that are quadratic in the incident field and involves an even-order optical susceptibility (for a molecule, its associated hyperpolarizability). In consequence, HRS is often regarded as formally forbidden in centrosymmetric media. However, for the fundamental three-photon interaction, theory based on the standard electric dipole approximation, representable as E1{sup 3}, does not account for all experimental observations. The relevant results emerge upon extending the theory to include E1{sup 2}M1 and E1{sup 2}E2 contributions, incorporating one magnetic dipolar or electric quadrupolar interaction, respectively, to a consistent level of multipolar expansion. Both additional interactions require the deployment of higher orders in the multipole expansion, with the E1{sup 2}E2 interaction analogous in rank and parity to a four-wave susceptibility. To elicit the correct form of response from fluid or disordered media invites a tensor representation which does not oversimplify the molecular components, yet which can produce results to facilitate the interpretation of experimental observations. The detailed derivation in this work leads to results which are summarized for the following: perpendicular detection of polarization components both parallel and perpendicular to the pump radiation, leading to distinct polarization ratio results, as well as a reversal ratio for forward scattered circular polarizations. The results provide a route to handling data with direct physical interpretation, to enable the more sophisticated design of molecules with sought nonlinear optical properties.

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

The Response of Precipitation Extremes to the Twentieth- and Twenty-First-Century Global Temperature Change in a Comprehensive Suite of CESM1 Large Ensemble Simulation: Revisiting the Role of Forcing Agents Vs. the Role of Forcing Magnitudes

The response of precipitation extremes (PEs) to global warming is found to be nonlinear in Community Earth System Model version 1 (CESM1) and other global climate models (Pendergrass et al., 2019), which led to the concern that it is not accurate to approximate the response of PE to a single forcing as the difference between simulations with all forcing agents and those that exclude one specific forcing. This calls into question previous model-based results that the sensitivity of PE with warming due to aerosol forcing is significantly larger than that due to greenhouse gases (GHGs). We reevaluate the PE sensitivity to GHGs and aerosols using available CESM1 ensemble simulations. We show that although the PE response to warming is nonlinear in CESM1, especially for the high warming projected in the twenty-first-century, PE sensitivity to aerosols is still significantly stronger than that due to GHGs when the comparison is made within similar warming regimes to avoid the bias induced by the nonlinear behavior. But the difference is smaller than previously estimated. We also conclude that the additivity assumption is largely valid to isolate the PE response due to aerosol forcing from the paired simulations including the “all forcing” experiment when the warming regime is small (e.g., 1°C–2°C in the next few decades when aerosol forcing is projected to decline and becomes a major source of uncertainty for model projection).

58 GEOSCIENCES↗
DOE OSTI2015

Postprocessing Mixed Finite Element Methods For Solving Cahn–Hilliard Equation: Methods and Error Analysis

A postprocessing technique for mixed finite element methods for the Cahn-Hilliard equation is developed and analyzed. Once the mixed finite element approximations have been computed at a fixed time on the coarser mesh, the approximations are postprocessed by solving two decoupled Poisson equations in an enriched finite element space (either on a finer grid or a higher-order space) for which many fast Poisson solvers can be applied. The nonlinear iteration is only applied to a much smaller size problem and the computational cost using Newton and direct solvers is negligible compared with the cost of the linear problem. The analysis presented here shows that this technique remains the optimal rate of convergence for both the concentration and the chemical potential approximations. The corresponding error estimate obtained in our paper, especially the negative norm error estimates, are non-trivial and different with the existing results in the literatures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Derivative-free stochastic optimization via adaptive sampling strategies

In this paper, we present a novel derivative-free framework for solving unconstrained stochastic optimization problems. Many problems in fields ranging from simulation optimization to reinforcement learning to quantum computing involve settings where only stochastic function values are obtained via a zeroth-order oracle, which has no available gradient information and necessitates the usage of derivative-free optimization methodologies. Our approach includes estimating gradients using stochastic function evaluations and integrating adaptive sampling techniques to control the accuracy in these stochastic approximations. Our framework encapsulates several gradient estimation techniques, including standard finite-difference, Gaussian smoothing, sphere smoothing, randomized coordinate finite-difference, and randomized subspace finite-difference methods. We provide theoretical convergence guarantees for our framework and analyze the worst-case iteration and sample complexities associated with each gradient estimation method. Finally, we demonstrate the empirical performance of the methods on logistic regression and nonlinear least squares problems.

Adaptive sampling↗
DOE OSTI2017

Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method

The Molecular Beam Epitaxial model is derived from the variation of a free energy, that consists of either a fourth order Ginzburg–Landau double well potential or a nonlinear logarithmic potential in terms of the gradient of a height function. One challenge in solving the MBE model numerically is how to develop proper temporal discretization for the nonlinear terms in order to preserve energy stability at the time-discrete level. In this paper, we resolve this issue by developing a first and second order time-stepping scheme based on the “Invariant Energy Quadratization” (IEQ) method. The novelty is that all nonlinear terms are treated semi-explicitly, and the resulted semi-discrete equations form a linear system at each time step. Moreover, the linear operator is symmetric positive definite and thus can be solved efficiently. We then prove that all proposed schemes are unconditionally energy stable. The semi-discrete schemes are further discretized in space using finite difference methods and implemented on GPUs for high-performance computing. Various 2D and 3D numerical examples are presented to demonstrate stability and accuracy of the proposed schemes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Microstructural characterization of a star-linear polymer blend under shear flow by using rheo-SANS

We present an investigation into the dynamic relaxation mechanisms of a polybutadiene blend composed of a four-arm star (10 wt. %) and a linear polymer matrix in the presence of an applied shear flow. Our focus was the response of the star polymer, which cannot be unambiguously assessed via linear viscoelastic measurements since the signature of the star polymer can barely be detected due to the dominant contribution of the linear matrix. By utilizing small-angle neutron scattering (SANS) coupled with a Couette shear device and a deuterated matrix polymer, we investigated the dynamics of the minority star component of the blend. Our results confirm that the stars deform anisotropically with increasing shear rate. We have compared the SANS data with predictions from the well-established scattering adaptation of the state-of-the-art tube model for entangled linear polymer melts undergoing shear, i.e., Graham, Likhtman, Milner, and McLeish (GLaMM) approach, appropriately modified following earlier studies in order to apply to the star. This modified model, GLaMM-R, includes the physics necessary to understand stress relaxation in both the linear and nonlinear flow regimes, i.e., contour length fluctuations, constraint release, convective constraint release, and chain retraction. The full scattering signal is due to the minority star component and, although the contribution of the linear chains is hidden from the neutron scattering, they still influence the star polymer molecular dynamics, with the applied shear rate ranging from approximately 8 to 24 s –1 , below the inverse relaxation time of the linear component. This study provides another confirmation that the combination of rheology and neutron scattering is an indispensable tool for investigating the nonlinear dynamics of complex polymeric systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On unifying randomized methods for inverse problems

This work unifies the analysis of various randomized methods for solving linear and nonlinear inverse problems with Gaussian priors by framing the problem in a stochastic optimization setting. By doing so, we show that many randomized methods are variants of a sample average approximation (SAA). More importantly, we are able to prove a single theoretical result that guarantees the asymptotic convergence for a variety of randomized methods. Additionally, viewing randomized methods as an SAA enables us to prove, for the first time, a single non-asymptotic error result that holds for randomized methods under consideration. Another important consequence of our unified framework is that it allows us to discover new randomization methods. Here, we present various numerical results for linear, nonlinear, algebraic, and PDE-constrained inverse problems that verify the theoretical convergence results and provide a discussion on the apparently different convergence rates and the behavior for various randomized methods.

42 ENGINEERING↗