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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A fully coupled two-level Schwarz preconditioner based on smoothed aggregation for the transient multigroup neutron diffusion equations

The neutron multigroup diffusion equations is an important approximation of the neutron transport equation, and has been widely used for studying the motion of neutrons and their interactions with the stationary background materials. Solving the multigroup diffusion equations faces challenging because many variables are coupled together and a lot of computing resources are required. In this paper, we focus on the development of a scalable parallel algorithm framework for the system of equations arising from the discretization of the multigroup diffusion equations using a finite element method in space and using an implicit scheme in time. The parallel agorithm framweork consists of an inexact Jacobian-free Newton for the nonlinear equations, a Krylov subspace method for the Jacobian system and a parallel Schwarz preconditioner for the acceleration of the algorithm convergence. A coarse space based on smoothed aggregation is introduced to construct a two-level method to improve the parallel performance of the Schwarz preconditioner. We numerically shown that compared with the one-level method, the proposed two-level method works significantly better in terms of the number of linear iterations and the compute time for a system of equations with millions of unknowns and a supercomputer with thousands of processors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics

Quantum computers can be used to simulate nonlinear non-Hamiltonian classical dynamics on phase space by using the generalized Koopman–von Neumann formulation of classical mechanics. The Koopman–von Neumann formulation implies that the conservation of the probability distribution function on phase space, as expressed by the Liouville equation, can be recast as an equivalent Schrödinger equation on Hilbert space with a Hermitian Hamiltonian operator and a unitary propagator. This Schrödinger equation is linear in the momenta because it derives from a constrained Hamiltonian system with twice the classical phase-space dimension. A quantum computer with finite resources can be used to simulate a finite-dimensional approximation of this unitary evolution operator. Quantum simulation of classical dynamics is exponentially more efficient than a deterministic Eulerian discretization of the Liouville equation if the Koopman–von Neumann Hamiltonian is sparse. Utilizing quantum walk techniques for state preparation and amplitude estimation for the calculation of observables leads to a quadratic improvement over classical probabilistic Monte Carlo algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A mathematical and numerical framework for the analysis of compressible thermal convection in gases at very high temperatures

The relevance of non-equilibrium phenomena, nonlinear behavior, gravitational effects and fluid compressibility in a wide range of problems related to high-temperature gas-dynamics, especially in thermal, mechanical and nuclear engineering, calls for a concerted approach using the tools of the kinetic theory of gases, statistical physics, quantum mechanics, thermodynamics and mathematical modeling in synergy with advanced numerical strategies for the solution of the Navier–Stokes equations. The reason behind such a need is that in many instances of relevance in this field one witnesses a departure from canonical models and the resulting inadequacy of standard CFD approaches, especially those traditionally used to deal with thermal (buoyancy) convection problems. Starting from microscopic considerations and typical concepts of molecular dynamics, passing through the Boltzmann equation and its known solutions, we show how it is possible to remove past assumptions and elaborate an algorithm capable of targeting the broadest range of applications. Moving beyond the Boussinesq approximation, the Sutherland law and the principle of energy equipartition, the resulting method allows most of the fluid properties (density, viscosity, thermal conductivity, heat capacity and diffusivity, etc.) to be derived in a rational and natural way while keeping empirical contamination to the minimum. Special attention is deserved as well to the well-known pressure issue. With the application of the socalled multiple pressure variables concept and a projection-like numerical approach, difficulties with such a term in the momentum equation are circumvented by allowing the hydrodynamic pressure to decouple from its thermodynamic counterpart. The final result is a flexible and modular framework that on the one hand is able to account for all the molecule (translational, rotational and vibrational) degrees of freedom and their effective excitation, and on the other hand can guarantee adequate interplay between molecular and macroscopic-level entities and processes. Performances are demonstrated by computing some incompressible and compressible benchmark test cases for thermal (gravitational) convection, which are then extended to the high-temperature regime taking advantage of the newly developed features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Central Compact Finite‐Difference Scheme With High Spectral Resolution for KdV Equation

This work presents a combination of cell‐node and cell‐centered compact finite difference scheme for the approximation of third derivatives involved in Korteweg–de Vries (KdV) equations. This approach employs a half‐shifted derivative construction at cell centers, avoiding the need for compact interpolation, thereby removing transfer errors; hence, it improves spectral resolution and maintains high‐order accuracy. Fourier analysis is performed to show the spectral properties of the proposed formulation, which provides higher spectral resolutions as compared to node‐based compact schemes. A filtering strategy is incorporated to suppress high‐frequency oscillations without compromising the accuracy of the numerical scheme, and the total variation diminishing Runge Kutta (TVDRK3) method is applied for time integration. Numerical experiments on linear, nonlinear, and coupled KdV systems are conducted, and a comparative analysis with cell‐node compact schemes confirms that the proposed scheme consistently reduces errors by up to an order of magnitude and achieves high spectral resolution properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Projection-based multifidelity linear regression for data-scarce applications

Surrogate modeling for systems with high-dimensional quantities of interest remains challenging, particularly when training data are costly to acquire. This work develops multifidelity methods for multiple-input multiple-output linear regression targeting data-limited applications with high-dimensional outputs. Multifidelity methods integrate many inexpensive low-fidelity model evaluations with limited, costly high-fidelity evaluations. We introduce two projection-based multifidelity linear regression approaches with linear and nonlinear features that leverage principal component basis vectors for dimensionality reduction and combine multifidelity data through: (i) a direct data augmentation using low-fidelity data, and (ii) a data augmentation incorporating explicit linear corrections between low-fidelity and high-fidelity data. The data augmentation approaches combine high-fidelity and low-fidelity data into a unified training set and train the linear regression model through weighted least squares with fidelity-specific weights. We introduce a proximity-based weighting scheme with automatic weight selection strategy through cross-validation. Here, the proposed multifidelity linear regression methods are demonstrated on approximating the surface pressure field of a hypersonic vehicle in flight and the temperature field on an aircraft disc braking system. In an ultra low-data regime of no more than twelve high-fidelity samples, multifidelity linear regression achieves approximately 2% – 12% improvement in median accuracy and a higher R 2 score relative to single-fidelity methods at comparable computational cost.

data augmentation↗
DOE OSTI2018

Wilson Prize article: Reflections on our experiences with developing the theory of intrabeam scattering

We discuss our experiences with developing detailed theoretical descriptions of intrabeam scattering in particle accelerators. We focus on the historical importance of understanding intrabeam scattering for the successful operation of a variety of accelerators around the world. In doing so, we highlight the fact that the theoretical understanding of intrabeam scattering played a crucial role in the discovery of the top quark at Fermilab, intermediate vector bosons W ± , Z and the Higgs particle at CERN, and the perfect liquid quark-gluon plasma at Brookhaven’s Relativistic Heavy Ion Collider. We describe several useful high energy approximations to intrabeam scattering, including those that utilize a Modified Piwinski high energy approximation by Karl Bane that has gained wide usage in applications to electron damping rings and advanced light sources. Finally, we comment on the fact that a detailed understanding of intrabeam scattering at synchrotron-based advanced light sources is empowering many transformational discoveries in a myriad of disciplines.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

In vivo dosimetry with optically stimulated dosimeters and RTQA2 radiochromic film for intraoperative radiotherapy of the breast

Purpose: Measurements were taken with optically stimulated luminescence dosimeters (OSLDs) and with RTQA2 radiochromic film to evaluate the use of each for in vivo dosimetry with intraoperative radiotherapy of the breast.Methods: Nonlinear calibration curves were established for OSLDs and RTQA2 radiochromic film using the Intrabeam 50 kV{sub p} source. Measurements were taken in a water phantom and compared to absolute dose measurements taken with an ionization chamber to investigate the characteristics of both types of dosimeters, including energy response and radiative absorption. In vivo readings were taken on the skin and in the tumor cavity using OSLDs and RTQA2 radiochromic film for 10 patients and 20 patients respectively. A prescription of 20 Gy to the surface of the applicator was used for all in vivo measurements in this study.Results: OSLDs were found to have an approximate uncertainty of ±7% for readings near the surface of the applicator and ±17% for readings at distances typical to the skin. The radiative absorption by OSLD was negative, indicating that this type of dosimeter absorbs less radiation than water in the targeted intraoperative radiotherapy energy range. RTQA2 film exhibited no energy dependence and all film readings were within ±8% of the delivered dose. The maximum radiative absorption in film was 8.5%. Radiochromic film measurements were found to be on average 18.2 ± 3.3 Gy for the tumor cavity and 2.1 ± 0.8 Gy for positions on the skin superior and inferior to the Intrabeam applicator. Average cavity measurements taken with OSLDs were 15.9 ± 3.9 Gy and average skin doses were 1.4 ± 0.8 Gy.Conclusions: OSLDs produce results with an uncertainty comparable to other dosimeters near the surface of the applicator but the uncertainty increases to an unacceptably high level with distance from the applicator. RTQA2 radiochromic film is shown to be accurate both at the surface of the applicator and at distances of 1–2 cm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Strongly nonlinear waves in locally resonant granular chains

In this paper, we explore a recently proposed locally resonant granular system bearing harmonic internal resonators in a chain of beads interacting via Hertzian elastic contacts. In this system, we propose the existence of two types of configurations: (a) small-amplitude periodic traveling waves and (b) dark-breather solutions, i.e. exponentially localized, time-periodic states mounted on top of a non-vanishing background. A remarkable feature distinguishing our results from other settings where dark breathers are observed is the complete absence of precompression in the system, i.e. the absence of a linear spectral band. We also identify conditions under which the system admits long-lived bright breather solutions. Our results are obtained by means of an asymptotic reduction to a suitably modified version of the so-called discrete p-Schrödinger (DpS) equation, which is established as controllably approximating the solutions of the original system for large but finite times (under suitable assumptions on the solution amplitude and the resonator mass). The findings are also corroborated by detailed numerical computations. Long-lived bright breathers are proved to exist over long but finite times, after which numerical simulations indicate that the breathers disintegrate. Finally, in line with these results, we prove that the only exact time-periodic bright breathers consist of trivial linear oscillations, without contact interactions between discrete elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fragme∩t: An Open‐Source Framework for Multiscale Quantum Chemistry Based on Fragmentation

Fragment-based quantum chemistry offers a means to circumvent the nonlinear computational scaling of conventional electronic structure calculations, by partitioning a large calculation into smaller subsystems then considering the many-body interactions between them. Variants of this approach have been used to parameterize classical force fields and machine learning potentials, applications that benefit from interoperability between quantum chemistry codes. However, there is a dearth of software that provides interoperability yet is purpose-built to handle the combinatorial complexity of fragment-based calculations. To fill this void we introduce “Fragme∩t”, an open-source software application that provides a tool for community validation of fragment-based methods, a platform for developing new approximations, and a framework for analyzing many-body interactions. Fragme∩t includes algorithms for automatic fragment generation and structure modification, and for distance- and energy-based screening of the requisite subsystems. Checkpointing, database management, and parallelization are handled internally and results are archived in a portable database. Interfaces to various quantum chemistry engines are easy to write and exist already for Q-Chem, PySCF, xTB, Orca, CP2K, MRCC, Psi4, NWChem, GAMESS, and MOPAC. Applications reported here demonstrate parallel efficiencies around 96% on more than 1000 processors but also showcase that the code can handle large-scale protein fragmentation using only workstation hardware, all with a codebase that is designed to be usable by non-experts. Fragme∩t conforms to modern software engineering best practices and is built upon well established technologies including Python, SQLite, and Ray. The source code is available under the Apache 2.0 license.

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

Spontaneous Polarization Suppression of Exciton–Exciton Annihilation in Rhombohedral-Stacked Bilayer Molybdenum Disulfide

Rapid exciton–exciton annihilation (EEA) in two-dimensional semiconductors limits access to high-density excitonic regimes essential for efficient optoelectronic operation under strong excitation. Here, we show that EEA is strongly suppressed in rhombohedral (3R)-stacked MoS 2 bilayers relative to nonpolar 2H bilayers, and that this suppression is quantitatively consistent with repulsive dipole–dipole interactions between layer-polarized excitons induced by spontaneous polarization. Using ultrafast pump–probe spectroscopy, we measure an EEA rate of γ EEA = (5.03 ± 0.99) × 10 –3 cm 2 s –1 in 3R bilayers, which is approximately 18.2-fold smaller than that in monolayers and 2.9-fold smaller than that in nonpolar 2H bilayers. Despite the higher exciton diffusivity recently reported for 3R relative to 2H bilayers, the reduced EEA rate in 3R indicates a rate-limited regime governed by the close-encounter annihilation probability rather than diffusion. A rate-limited annihilation model incorporating a dipole–dipole repulsive potential captures the observed ratio γ EEA,3R /γ EEA,2H ≈ 0.35 for an exciton–exciton encounter distance of ∼1.3 nm, consistent with the bilayer exciton Bohr radius. These results indicate that spontaneous polarization in 3R-stacked bilayers can suppress nonlinear excitonic losses and provide a route toward high-density excitonics.

dipolar excitons↗
DOE OSTI2022

Stability of exact solutions of the (2 + 1)-dimensional nonlinear Schrödinger equation with arbitrary nonlinearity parameter κ

In this work, we consider the nonlinear Schrödinger equation (NLSE) in 2+1 dimensions with arbitrary nonlinearity exponent κ in the presence of an external confining potential. Exact solutions to the system are constructed, and their stability as we increase the 'mass' (i.e., the L 2 norm) and the nonlinearity parameter κ is explored. Here we observe both theoretically and numerically that the presence of the confining potential leads to wider domains of stability over the parameter space compared to the unconfined case. Our analysis suggests the existence of a stable regime of solutions for all κ as long as their mass is less than a critical value M*(κ). Furthermore, we find that there are two different critical masses, one corresponding to width perturbations and the other one to translational perturbations. The results of Derrick's theorem are also obtained by studying the small amplitude regime of a four-parameter collective coordinate (4CC) approximation. A numerical stability analysis of the NLSE shows that the instability curve M*(κ) versus κ lies below the two curves found by Derrick's theorem and the 4CC approximation. In the absence of the external potential, κ = 1 demarcates the separation between the blowup regime and the stable regime. In this 4CC approximation, for κ < 1, when the mass is above the critical mass for the translational instability, quite complicated motions of the collective coordinates are possible. Energy conservation prevents the blowup of the solution as well as confines the center of the solution to a finite spatial domain. We call this regime the 'frustrated' blowup regime and give some illustrations. In an appendix, we show how to extend these results to arbitrary initial ground state solution data and arbitrary spatial dimension d.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Optical, physical, and chemical properties of springtime aerosol over Barrow Alaska in 2008

Airborne observations from four flights during the 2008 Indirect and Semi-Direct Aerosol Campaign (ISDAC) are used to examine some cloud-free optical, physical, and chemical properties of aerosol particles in the springtime Arctic troposphere. The number concentrations of particles larger than 0.12 μm (Na>120), important for light extinction and cloud droplet formation, ranged from 15 to 2260 cm -3 , with the higher Na>120 cases dominated by measurements from two flights of long-range transported biomass burning (BB) aerosols. The two other flights examined here document a relatively clean aerosol and an Arctic Haze aerosol impacted by larger particles largely composed of dust. For observations from the cleaner case and the BB cases, the particle light scattering coefficients at low relative humidity (RH<20%) increased nonlinearly with increasing Na>120, driven mostly by an increase in mean sizes of particles with increasing Na>120 (BB cases). For those three cases, particle light absorption coefficients also increased nonlinearly with increasing Na>120 and linearly with increasing submicron particle volume concentration. In addition to black carbon, brown carbon was estimated to have increased light absorption coefficients by 27% (450 nm wavelength) and 14% (550 nm) in the BB cases. For the case with strong dust influence, the absorption relative to submicron particle volume was small compared with the other cases. There was a slight gradient of Passive Cavity Aerosol Spectrometer Probe (PCASP) mean volume diameter (MVD) towards smaller sizes with increasing height, which suggests more scavenging of the more elevated particles, consistent with a typically longer lifetime of particles higher in the atmosphere. However, in approximately 10% of the cases, the MVD increased (>0.4 μm) with increasing altitude, suggesting transport of larger fine particle mass (possibly coarse particle mass) at high levels over the Arctic. This may be because of transport of larger particles at higher elevations and relatively slow deposition to the surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling information flow in a computer processor with a multi-stage queuing model

In this paper, we introduce a nonlinear stochastic model to describe the propagation of information inside a computer processor. In this model, a computational task is divided into stages, and information can flow from one stage to another. The model is formulated as a spatially-extended, continuous-time Markov chain where space represents different stages. This model is equivalent to a spatially-extended version of the M/M/s queue. The main modeling feature is the throttling function which describes the processor slowdown when the amount of information falls below a certain threshold. We derive the stationary distribution for this stochastic model and develop a closure for a deterministic ODE system that approximates the evolution of the mean and variance of the stochastic model. In conclusion, we demonstrate the validity of the closure with numerical simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Magnetoacoustic Waves in a Stratified Atmosphere with a Magnetic Null Point

We perform nonlinear MHD simulations to study the propagation of magnetoacoustic waves from the photosphere to the low corona. We focus on a 2D system with a gravitationally stratified atmosphere and three photospheric concentrations of magnetic flux that produce a magnetic null point with a magnetic dome topology. We find that a single wavepacket introduced at the lower boundary splits into multiple secondary wavepackets. A portion of the packet refracts toward the null owing to the varying Alfvén speed. Waves incident on the equipartition contour surrounding the null, where the sound and Alfvén speeds coincide, partially transmit, reflect, and mode-convert between branches of the local dispersion relation. Approximately 15.5% of the wavepacket’s initial energy ( E {sub input}) converges on the null, mostly as a fast magnetoacoustic wave. Conversion is very efficient: 70% of the energy incident on the null is converted to slow modes propagating away from the null, 7% leaves as a fast wave, and the remaining 23% (0.036 E {sub input}) is locally dissipated. The acoustic energy leaving the null is strongly concentrated along field lines near each of the null’s four separatrices. The portion of the wavepacket that refracts toward the null, and the amount of current accumulation, depends on the vertical and horizontal wavenumbers and the centroid position of the wavepacket as it crosses the photosphere. Regions that refract toward or away from the null do not simply coincide with regions of open versus closed magnetic field or regions of particular field orientation. We also model wavepacket propagation using a WKB method and find that it agrees qualitatively, though not quantitatively, with the results of the numerical simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A robust return-map algorithm for general multisurface plasticity

Three new contributions to the field of multisurface plasticity are presented for general situations with an arbitrary number of nonlinear yield surfaces with hardening or softening. A method for handling linearly dependent flow directions is described. A residual that can be used in a line search is defined. An algorithm that has been implemented and comprehensively tested is discussed in detail. Examples are presented to illustrate the computational cost of various components of the algorithm. The overall result is that a single Newton-Raphson iteration of the algorithm costs between 1.5 and 2 times that of an elastic calculation. Examples also illustrate the successful convergence of the algorithm in complicated situations. For example, without using the new contributions presented here, the algorithm fails to converge for approximately 50% of the trial stresses for a common geomechanical model of sedementary rocks, while the current algorithm results in complete success. Since it involves no approximations, the algorithm is used to quantify the accuracy of an efficient, pragmatic, but approximate, algorithm used for sedimentary-rock plasticity in a commercial software package. Furthermore, the main weakness of the algorithm is identified as the difficulty of correctly choosing the set of initially active constraints in the general setting.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Longitudinal Beam Dynamics with a Higher-Harmonic Cavity for Bunch Lengthening

In Sect.1 we review the conditions of operation of an HHC system assuming its stability. We first discuss the longitudinal dynamics with an active HHC and the conditions that lead to a quartic potential in the limit of small oscillations, together with formulae giving the bunch lengthening factor and the synchrotron frequency dependence on the amplitude of the synchrotron oscillations. We discuss next passive HHC operations with uniform fillings, and distinguish the case of normal-conducting HHC, where optimal conditions of operation can be met, from the case of superconducting HHC, where optimal conditions of operation can be satisfied only approximately. In Sect.2 we address the stability of the passive superconducting 3HC system of the NSLS-II storage ring via self-consistent simulations with code SPACE. The case of nominal fractional filling of 80% is compared with the case of fractional filling of 90% and with the uniform filling case. It is shown that for values of the detuning frequency in the neighborhood of the good working point discussed in Sect.2, unstable regimes of oscillation are found, with detuning frequency threshold dependent on the gap in the filling. It is shown that the instability regime is characterized by dipole and quadrupole modes of oscillation exhibiting asymptotically nonlinear saturation. A discussion of the average bunch lengthening, together with the degree of uniformity across the bunch train, achievable under stable conditions as a function of the fractional filling concludes the paper.

36 MATERIALS SCIENCE↗