SEARCH · Engineering Papers

Results for “Non linear dynamics”

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

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS2021

Truncated ARQ Statistical Link Analysis for Dynamic Links

The future deep space links are migrating towards higher frequency bands such as Ka band and optical. These links are susceptible to non Gaussian and non linear effects such as atmospheric turbulence, scintillation, antenna mis-pointing, jitter, etc. These dynamic links thus will experience various degrees of fading loss, and some of these link disruptions cannot be effectively mitigated by forward error correction coding and/or interleaving. One effective way to ensure reliable communication is by using Automatic Repeat Request (ARQ) protocol, where the receiver acknowledges to the transmitter whether or not a data unit is successfully received. If a data unit is not successfully received (such as after a pre-set time-out), the transmitter would then re-transmit the lost data unit to the receiver. In a previous paper, we derived a statistical link analysis method of finding the optimal operating Signal-to-Noise Ratio (SNR) and estimating the latency of an ARQ scheme. In a more recent paper, we demonstrated the above method using the SNR distribution constructed from the Ka-band (32 GHz) flight data. To simplify the discussion, we considered the academic approach that the ARQ scheme allows for an infinite number of retransmissions. In this paper, we consider the more practical case of a truncated ARQ scheme, where there is a limit on the number of retransmissions. We derive the error probability, the optimal SNR setting, and the latency statistics of the correctly received frames of the truncated ARQ schemes. We first discuss the truncated ARQ link analysis principles using the Gaussian assumption for SNR distribution with a large variance. Next, we demonstrate the statistical truncated ARQ link analysis using the SNR distribution constructed from the Ka-band flight data. The results in this paper can be applied in the design of reliable communication systems such as the Consultative Committee for Space Data System (CCSDS) File Transfer Protocol (CFTP) and the Delay Tolerant Network (DTN).

Morabito, David↗

DOE OSTI2021

Plasma Lens in Parametric Resonance Ionization Cooling

We present a concept of a plasma lens in a parametric resonance ionization cooling (PIC) for extra cooling of muon colliders. The PIC concept has been developed to overcome most aberrations. However, parasitic non-linear aberrations which is induced in an ionization cooling material still reduce a dynamic aperture of the PIC. The plasma lens generates an azimuthally symmetric focusing field in gas-filled RF cavities that is several orders of magnitude stronger than the conventional superconducting magnets. It allows for reduction of the beam size at the beam expansion points. This reduces the size of the aberrations and therefore greatly simplifies their compensation.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Assessing Dynamic Behaviors in Converter- Dominated Power Systems via RMS and EMT Simulations: A Study of Hawaii’s NELHA Microgrid

The transition from conventional power systems to converter-based microgrids has significantly advanced sustainability, clean energy integration, and operational reliability. However, this paradigm shift introduces operational challenges due to the intermittent nature of renewable energy sources and the non-linear characteristics of power electronic loads, inducing voltage fluctuations and harmonic distortions that complicate voltage and frequency regulation. Accurate dynamic modeling is hypothesized to be critical for capturing such effects, enabling reliable simulation and control strategy development. This study introduces an innovative dynamic modeling framework for a real-world converter-based microgrid, utilizing both root mean square (RMS) and electromagnetic transient (EMT) simulation methods. The microgrid was modeled in DIgSILENT PowerFactory, with simulations calibrated against high-resolution field measurements from SEL-735 power quality meters. Results show that RMS simulations effectively characterize steady-state dynamics, while EMT simulations are essential for capturing high-frequency transients and non-linear effects from photovoltaic inverters and variable frequency drives (VFDs). This complementary approach provides a comprehensive understanding of microgrid behavior, providing critical insights for improving simulation accuracy, advancing protection schemes, and improving resilience in future low-inertia power networks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Polariton spectra under the collective coupling regime. II. 2D non-linear spectra

In our previous work [Mondal et al., J. Chem. Phys. 162, 014114 (2025)], we developed several efficient computational approaches to simulate exciton–polariton dynamics described by the Holstein–Tavis–Cummings (HTC) Hamiltonian under the collective coupling regime. Here, we incorporated these strategies into the previously developed Lindblad-partially linearized density matrix (⁠$\mathscr{L}$-PLDM) approach for simulating 2D electronic spectroscopy (2DES) of exciton–polariton under the collective coupling regime. In particular, we apply the efficient quantum dynamics propagation scheme developed in Paper I to both the forward and the backward propagations in the PLDM and develop an efficient importance sampling scheme and graphics processing unit vectorization scheme that allow us to reduce the computational costs from $\mathscr{O}$($\mathscr{K}$ 2 )$\mathscr{O}$(T 3 ) to $\mathscr{O}$($\mathscr{K}$)$\mathscr{O}$(T 0 ) for the 2DES simulation, where $\mathscr{K}$ is the number of states and T is the number of time steps of propagation. As a result, we further simulated the 2DES for an HTC Hamiltonian under the collective coupling regime and analyzed the signal from both rephasing and non-rephasing contributions of the ground state bleaching, excited state emission, and stimulated emission pathways.

2D non-linear spectra↗

DOE OSTI2023

Direct Estimation of Parameters in ODE Models Using WENDy: Weak-Form Estimation of Nonlinear Dynamics

Abstract We introduce the Weak-form Estimation of Nonlinear Dynamics (WENDy) method for estimating model parameters for non-linear systems of ODEs. Without relying on any numerical differential equation solvers, WENDy computes accurate estimates and is robust to large (biologically relevant) levels of measurement noise. For low dimensional systems with modest amounts of data, WENDy is competitive with conventional forward solver-based nonlinear least squares methods in terms of speed and accuracy. For both higher dimensional systems and stiff systems, WENDy is typically both faster (often by orders of magnitude) and more accurate than forward solver-based approaches. The core mathematical idea involves an efficient conversion of the strong form representation of a model to its weak form, and then solving a regression problem to perform parameter inference. The core statistical idea rests on the Errors-In-Variables framework, which necessitates the use of the iteratively reweighted least squares algorithm. Further improvements are obtained by using orthonormal test functions, created from a set of $$C^{\infty }$$ C ∞ bump functions of varying support sizes.We demonstrate the high robustness and computational efficiency by applying WENDy to estimate parameters in some common models from population biology, neuroscience, and biochemistry, including logistic growth, Lotka-Volterra, FitzHugh-Nagumo, Hindmarsh-Rose, and a Protein Transduction Benchmark model. Software and code for reproducing the examples is available at https://github.com/MathBioCU/WENDy .

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Diffusion in intact secondary cell wall models of plants at different equilibrium moisture content

Secondary plant cell walls are composed of carbohydrate and lignin polymers, and collectively represent a significant renewable resource. Leveraging these resources depends in part on a mechanistic understanding for diffusive processes within plant cell walls. Common wood protection treatments and biomass conversion processes to create biorefinery feedstocks feature ion or solvent diffusion within the cell wall. X-ray fluorescence microscopy experiments have determined that ionic diffusion rates are dependent on cell wall hydration as well as the ionic species through non-linear relationships. In this work, we use classical molecular dynamics simulations to map the diffusion behavior of different plant cell wall components (cellulose, hemicellulose, lignin), ions (Na + , K + , Cu2 + , Cl - ) and water within a model for an intact plant cell wall at various hydration states (3-30 wt% water). From these simulations, we analyze the contacts between different plant cell wall components with each other and their interaction with the ions. Generally, diffusion increases with increasing hydration, with lignin and hemicellulose components increasing diffusion by an order of magnitude over the tested hydration range. Ion diffusion depends on charge. Positively charged cations preferentially interact with hemicellulose components, which include negatively charged carboxylates. As a result, positive ions diffuse more slowly than negatively charged ions. Measured diffusion coefficients are largely observed to best fit piecewise linear trends, with an inflection point between 10 and 15% hydration. These observations shed light onto the molecular mechanisms for diffusive processes within secondary plant cell walls at atomic resolution.

09 BIOMASS FUELS↗

DOE OSTI2021

Spectroscopy of GR1 centers in synthetic diamonds

We report on the spectroscopic characterization, absorption saturation, and excited-state dynamics of GR1 centers in synthetic diamonds. The non-linear optical measurements reveal an efficient bleaching of the GR1 center’s ground level under ns-pulsed 633 nm excitation. The maxima of absorption and emission cross sections were estimated to be 4.5 × 10 −17 cm 2 and 9 × 10 −17 cm 2 at 630 nm and 780 nm, respectively. The radiative lifetime of the excited level was estimated to be 8.5 ns. The 658 nm probe kinetics uncovered relaxation of 1 T 2 excited level going predominantly to a metastable state with a lifetime of 220 μs. An induced absorption detected with the use of a highly concentrated diamond sample could be due to up-conversion and photoionization processes in additional impurity-vacancy center with absorption at the short-wavelength tail of GR1. The results presented here indicate that synthetic diamonds with GR1 centers could serve as near infrared gain media or passive Q-switchers for laser cavities over the 633–750 nm spectral range. Optimization of the center concentration and the crystal parameters are required to minimize the induced absorption at the lasing wavelengths.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Mitigation of the Brighter-Fatter Effect in the LSST Camera

Thick, fully depleted charge-coupled devices (CCDs) are known to exhibit non-linear behavior at high signal levels due to the dynamic behavior of charges collecting in the potential wells of pixels, called the brighter-fatter effect (BFE). This particularly impacts bright calibration stars, which appear larger than their intrinsic shape, creating a flux-dependent point-spread function (PSF) that if left unmitigated, could make up a large fraction of the error budget in Stage IV weak-lensing (WL) surveys such as the Legacy Survey of Space and Time (LSST). In this paper, we analyze image measurements of flat fields and artificial stars taken at different illumination levels with the LSST Camera (LSSTCam) at SLAC National Accelerator Laboratory in order to quantify this effect in the LSST Camera before and after a previously introduced correction technique. We observe that the BFE evolves anisotropically as a function of flux due to higher-order BFEs, which violates the fundamental assumption of this correction method. We then introduce a new sampling method based on a physically motivated model to account these higher-order terms in the correction, and then we test the modified correction on both datasets. We find that the new method corrects the effect in flat fields better than it corrects the effect in artificial stars which we conclude is the result of a unmodeled curl component of the deflection field by the correction. We use these results to define a new metric for the full-well capacity of our sensors and advise image processing strategies to further limit the impact of the effect on LSST WL science pathways.

47 OTHER INSTRUMENTATION↗

NASA NTRS1975

Non-linear resonances in the forced responses of plates. I - Symmetric responses of circular plates

The dynamic analogue of the von Karman equations is used to study the symmetric response of a circular plate to a harmonic excitation when the frequency of the excitation is near one of the natural frequencies. It is shown that, in general, when there is no internal resonance (i.e., the natural frequencies are not commensurable), only the mode having a frequency near that of the excitation is strongly excited (i.e., is needed to represent the response in the first approximation). A clamped, circular plate is used as a numerical example to show that, when there is an internal resonance, more than one of the modes involved in this resonance can be strongly excited; moreover, when more than one mode is strongly excited, the lower modes can dominate the response, even when the frequency of the excitation is near that of the highest mode. This possibility was not revealed by any of the earlier studies which were based on the same governing equations.

Sridhar, S.↗

NASA NTRS1999

The Role of Gravity Waves in Modulating Atmospheric Tides

We discuss results for the diurnal and semidiurnal tides obtained from our 3-D, time dependent numerical spectral model (NMS), extending from the ground up into the thermosphere, which incorporates Hines' Doppler spread parameterization of small scale gravity waves (GW). In the DSP, GW momentum (and energy) are conserved as the waves modulate the background flow and are filtered by the flow.As a consequence, the GW interaction tightly couples the dynamic components of the middle atmosphere with strong non-linear interactions between mean zonal circulation, tides and planetary waves to produce complicated patterns of variability much like those observed. The major conclusions are: (1) Since GW momentum is deposited in the altitude regime of increasing winds, the amplitude of the diurnal tide is amplified and its vertical wavelength is reduced at altitudes between 80 and 120 km. Wave filtering by the mean zonal circulation (with peak velocities during solstice) causes the GW flux to peak during equinox, and this produces a large semi-annual variation in the tide that has been observed on UARS. (2) Without the diurnal tide, the semidiurnal tide would also be modulated in this way. But the diurnal tide filters out the GW preferentially during equinox, so that the semidiurnal tide, at higher altitudes, tends to peak during solstice. (3) Under the influence of GW, the tides are modulated also significantly by planetary waves, with periods between 2 and 30 days, which are generated preferentially during solstice in part due to baroclinic instability.

Mayr, H. G.↗

NASA NTRS2017

Wave-Particle Energy Exchange Directly Observed in a Kinetic Alfven-Branch Wave

Alfven waves are fundamental plasma wave modes that permeate the universe. At small kinetic scales they provide a critical mechanism for the transfer of energy between electromagnetic fields and charged particles. These waves are important not only in planetary magnetospheres, heliospheres, and astrophysical systems, but also in laboratory plasma experiments and fusion reactors. Through measurement of charged particles and electromagnetic fields with NASAs Magnetospheric Multiscale (MMS) mission, we utilize Earths magnetosphere as a plasma physics laboratory. Here we confirm the conservative energy exchange between the electromagnetic field fluctuations and the charged particles that comprise an undamped kinetic Alfven wave. Electrons confined between adjacent wave peaks may have contributed to saturation of damping effects via non-linear particle trapping. The investigation of these detailed wave dynamics has been unexplored territory in experimental plasma physics and is only recently enabled by high-resolution MMS observations.

magnetohydrodynamic↗

NASA NTRS1993

Dynamics and control of coherent structures in the turbulent wall layer: An overview

We expand the velocity field in the vicinity of the wall in empirical eigenfunctions obtained from experiment. Truncating our system and using Galerkin projection, we obtain a closed set of non-linear ordinary differential equations with ten degrees of freedom. We find a rich dynamical behavior, including in particular a heteroclinic attracting orbit giving rise to intermittency. The intermittent jump from one attracting point to the other resembles in many respects the bursts observed in experiments. Specifically, the time between jumps and the duration of the jumps, is approximately that observed in a burst; the jump begins with the formation of a narrowed and intensified updraft, like the ejection phase of a burst, and is followed by a gentle, diffuse downdraft like the sweep phase of a burst. The magnitude of the Reynolds stress spike produced during a burst is limited by our truncation. The behavior is quite robust, much of it being due to the symmetries present (Aubry's group has examined dimensions up to 128 with persistence of the global behavior). We have examined eigenvalues and coefficients obtained from experiment, and from exact simulation, which differ in magnitude. Similar behavior is obtained in both cases; in the latter case, the heteroclinic orbits connect limit cycles instead of fixed points, corresponding to cross-stream waving of the streamwise rolls. The bifurcation diagram remains structurally similar, but somewhat distorted. The role of the pressure term is made clear - it triggers the intermittent jumps, which otherwise would occur at longer and longer intervals, as the system trajectory is attracted closer and closer to the heteroclinic cycle. The pressure term results in the jumps occurring at essentially random times, and the magnitude of the signal determines the average timing. Stretching of the wall region shows that the model is consistent with observations of polymer drag reduction. Change of the third order coefficients, corresponding to acceleration or deceleration of the mean flow, changes the heteroclinic cycles from attracting to repelling, increasing or decreasing the stability, in agreement with observations. The existence of fixed points is an artifact introduced by the projection; however, a decoupled model still displays the rich dynamics. Numerous assumptions made in Aubry et al. (1988) can now be proved exactly. Feeding back eigenfuncitons with the proper phase can delay the bursting, (the heteroclinic jump to the other fixed point), decreasing the drag. It is also possible to speed up the bursting, increasing mixing to control separation. Our approach is optimal for short time tracking in control.

Berkooz, Gal↗

DOE OSTI2022

Probing ultrafast excited-state dynamics using EUV-IR six-wave-mixing emission spectroscopy

Non-linear processes such as four-wave-mixing have become instrumental in attosecond EUV spectroscopy. Using EUV high harmonics in conjunction with collinear near-infrared and mid-infrared fields, we extended this approach to high-order-mixing between three colors. Specifically, we find that atomic resonances in neon exhibit a significant cross section for six-wave mixing. The MIR probe frequency tunability in our multicolor scheme is employed to access several optically dark resonances and probe the quantum beat of a coherent electronic wavepacket using background-free EUV emission as a diagnostic. This technique can be easily extended to other atomic and molecular systems, and opens the door to multi-dimensional non-linear spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Unifying Energy-Management Problems for Inverter-Based Power Networks

This paper outlines a suite of energy-management problems for inverter-based power networks from the vantage point of optimal control and (non)linear optimization. The problems are categorized based on timescales dictated by the network dynamics, and organized methodologically based on the problem complexity. A growing body of literature has addressed problems in this domain, albeit, with poorly motivated assumptions and behavioral models that obscure precise device behavior. With a combination of circuit- and control-theoretic lenses, we establish appropriate dynamic models for the networked resources, illustrate how common engineering assumptions arise, uncover how problems are linked, and postulate open challenges.

(non)linear optimization↗

DOE OSTI2023

Quantum dynamics simulations of the 2D spectroscopy for exciton polaritons

We develop an accurate and numerically efficient non-adiabatic path-integral approach to simulate the non-linear spectroscopy of exciton–polariton systems. This approach is based on the partial linearized density matrix approach to model the exciton dynamics with explicit propagation of the phonon bath environment, combined with a stochastic Lindblad dynamics approach to model the cavity loss dynamics. Through simulating both linear and polariton two-dimensional electronic spectra, we systematically investigate how light–matter coupling strength and cavity loss rate influence the optical response signal. Our results confirm the polaron decoupling effect, which is the reduced exciton–phonon coupling among polariton states due to the strong light–matter interactions. We further demonstrate that the polariton coherence time can be significantly prolonged compared to the electronic coherence outside the cavity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Modelica-based system modeling for studying control-related faults in chiller plants and boiler plants serving large office buildings

System modeling is critical when studying operation faults in chiller plants and boiler plants. However, current fault models have difficulties faithfully representing the operation of chiller plants and boiler plants under the effects of those faults, especially the control-related ones. In this study, we present a systematic method to develop high-fidelity models for approximating the behaviors of chiller plants and boiler plants under faulty conditions. Compared to existing ones, the resulting fault models have two advantages: first, they better characterize the dynamic patterns in the system operation. In those models, control architecture and control logic are faithfully implemented. Thus, they can be used to study control-related faults, such as incorrect staging control due to sensor bias and mistuned feedback control. Second, they are readily extensible and can support large-scale investigations to explore different faulty conditions/scenarios. Those models are established in a hierarchical structure while modules in each layer can be redeclared and parameterized at upper layers. In such a case, modifications to the models can be realized through model modifiers and the process can be easily streamlined with scripts. We applied the proposed models in a comprehensive fault impact evaluation of the 13 control-related faults of chiller and boiler plants. In this evaluation, the proposed plant model is coupled with the EnergyPlus thermal load model to study the impact of various faulty scenarios. Based on the evaluation results, we identified the faults that have the most significant impacts on the operation of the chiller and boiler plants, respectively. We also found that the relationship between the impacts of the studied faults and the severity level of the faults are highly non-linear

Huang, Sen↗

NASA NTRS2004

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗