SEARCH · Engineering Papers

Results for “perturbation analysis”

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 955 records · Page 53

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Error analysis of numerical methods for thick diffusive neutron transport problems on Shishkin mesh

A thin layer will develop at the boundary if the incoming angular flux is anisotropic in thick diffusive neutron transport problems. Solving such singularly perturbed problems, which have non-smooth solutions with singularity near the boundary, is computationally challenging. Standard finite difference schemes on a uniform mesh cannot yield ε-uniform convergence, where ε is a small parameter, while it can be achieved on a suitable piecewise-uniform Shishkin mesh. We present a formal error analysis of the diamond difference (DD) method and step difference (SD) method for solving the S{sub N} neutron transport equation. The analysis can be extended to other finite difference methods. Numerical results are presented to confirm the error estimates and the advantages of the Shishkin mesh. (author)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Enhancing DESI DR1 full-shape analyses using HOD-informed priors

We present an analysis of DESI Data Release 1 (DR1) that incorporates Halo Occupation Distribution (HOD)-informed priors into Full-Shape (FS) modeling of the power spectrum based on cosmological perturbation theory (PT). By leveraging physical insights from the galaxy-halo connection, these HOD-informed priors on nuisance parameters substantially mitigate projection effects in extended cosmological models that allow for dynamical dark energy. The resulting credible intervals now encompass the posterior maximum from the baseline analysis using gaussian priors, eliminating a significant posterior shift observed in baseline studies. In the ΛCDM framework, a combined DESI DR1 FS information and constraints from the DESI DR1 baryon acoustic oscillations (BAO) — including Big Bang Nucleosynthesis (BBN) constraints and a weak prior on the scalar spectral index — yields Ω m = 0.2994 ± 0.0090 and σ 8 = 0.836$^{+0.024}_{-0.027}$, representing improvements of approximately 4% and 23% over the baseline analysis, respectively. For the w 0 w a CDM model, our results from various data combinations are highly consistent, with all configurations converging to a region with w 0 > -1 and w a < 0. This convergence not only suggests intriguing hints of dynamical dark energy but also underscores the robustness of our HOD-informed prior approach in delivering reliable cosmological constraints.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1994

Galactic cosmic ray modulation and interplanetary medium perturbations due to a long-living active region during October 1989

During October 1989, three very energetic flares were ejected by the same active region at longitudes 9 deg E, 32 deg W, and 57 deg W, respectively. The shape of the galactic cosmic ray variations suggests the presence of large magnetic cloud structures (Nagashima et al., 1990) following the shock-associated perturbations. In spite of long data gaps the interplanetary observations at Interplanetary Monitoring Platform (IMP) 8 (near the Earth) and International Cometary Explorer (ICE)(approximately 1 AU, approximately 65 deg W) confirm this possibility for the event related to the 9 deg E flare; the principal axes analysis shows that the interplanetary magnetic field variations at both spacecraft locations are mainly confined on a meridian plane. This result suggests that the western longitudinal extension of this cloud is indeed very large (greater than or equal to 5 deg). The nonnegligible depression in the cosmic ray intensity observed inside the possible cloud related to the 57 deg W flare indicates that also the eastern extension could be very wide. The analysis of neutron monitor data shows clearly the cosmic ray trapping effect of magnetic clouds; this mechanism seems to be responsible for the enhanced diurnal effect often observed during the recovery phase of Forbush decreases. We give an interpretation for the anisotropic cosmic ray peak occurring in the third event, and, related to that, we suggest that the Forbush decrease modulated region at the Earth's orbit could be somewhat wider than the magnetic cloud, as already anticipated by Nagashima et al. (1990). By this analysis, based mainly on cosmic ray data, we show that it is possible to do reasonable inferences on the large-scale structure of flare-related interplanetary perturbations when interplanetary medium data are not completely present.

Bavassano, B.↗

NASA NTRS1972

Integrated dynamic analysis simulation of space stations with controllable solar array

A methodology is formulated and presented for the integrated structural dynamic analysis of space stations with controllable solar arrays and non-controllable appendages. The structural system flexibility characteristics are considered in the dynamic analysis by a synthesis technique whereby free-free space station modal coordinates and cantilever appendage coordinates are inertially coupled. A digital simulation of this analysis method is described and verified by comparison of interaction load solutions with other methods of solution. Motion equations are simulated for both the zero gravity and artificial gravity (spinning) orbital conditions. Closed loop controlling dynamics for both orientation control of the arrays and attitude control of the space station are provided in the simulation by various generic types of controlling systems. The capability of the simulation as a design tool is demonstrated by utilizing typical space station and solar array structural representations and a specific structural perturbing force. Response and interaction load solutions are presented for this structural configuration and indicate the importance of using an integrated type analysis for the predictions of structural interactions.

Heinrichs, J. A.↗

DOE OSTI2020

Resummation of QED radiative corrections in a strong constant crossed field

By considering radiative corrections of up to 3rd-loop order, Ritus and Narozhny conjectured that the proper expansion parameter for QED in a strong constant crossed field is $g=\alpha\chi^{2/3}$, where the dynamical quantum parameter $\chi=e\sqrt{-(Fp)^2}/m^3$ combines the particle momentum $p$ with the external field strength tensor $F$. Here we present and discuss the first non-perturbative result in this context, the resummed bubble-type polarization corrections to the electron self-energy in a constant crossed field. Our analysis confirms the relevance of the scaling parameter $g$ to the enhancement of bubble-type radiative corrections. This parameter actually represents the characteristic value of the ratio of the 1-loop polarization bubble to the photon virtuality. After an all-order resummation we identify and discuss two contributions to the self-energy with different formation regions and asymptotic behavior for $g\gg1$. Whereas the breakdown of perturbation theory occurs already for $g\gtrsim1$, the leading-order result remains dominant until the asymptotic regime $g\gg 1$ is reached. However, the latter is specific to processes like elastic scattering or photon emission and does not have to remain true for general higher-order QED processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Variational and phase response analysis for limit cycles with hard boundaries, with applications to neuromechanical control problems

Motor systems show an overall robustness, but because they are highly nonlinear, understanding how they achieve robustness is difficult. In many rhythmic systems, robustness against perturbations involves response of both the shape and the timing of the trajectory. This makes the study of robustness even more challenging. To understand how a motor system produces robust behaviors in a variable environment, we consider a neuromechanical model of motor patterns in the feeding apparatus of the marine mollusk Aplysia californica (Shaw et al. in J Comput Neurosci 38(1):25–51, 2015; Lyttle et al. in Biol Cybern 111(1):25–47, 2017). We established in (Wang et al. in SIAM J Appl Dyn Syst 20(2):701–744, 2021. https://doi.org/10.1137/20M1344974) the tools for studying combined shape and timing responses of limit cycle systems under sustained perturbations and here apply them to study robustness of the neuromechanical model against increased mechanical load during swallowing. Interestingly, we discover that nonlinear biomechanical properties confer resilience by immediately increasing resistance to applied loads. In contrast, the effect of changed sensory feedback signal is significantly delayed by the firing rates’ hard boundary properties. Our analysis suggests that sensory feedback contributes to robustness in swallowing primarily by shifting the timing of neural activation involved in the power stroke of the motor cycle (retraction). This effect enables the system to generate stronger retractor muscle forces to compensate for the increased load, and hence achieve strong robustness. The approaches that we are applying to understanding a neuromechanical model in Aplysia, and the results that we have obtained, are likely to provide insights into the function of other motor systems that encounter changing mechanical loads and hard boundaries, both due to mechanical and neuronal firing properties.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1992

Symplectic maps for the n-body problem - Stability analysis

The stability of new symplectic n-body maps is examined from the point of view of nonlinear dynamics. The resonances responsible for the principal artifacts are identified. These are resonances between the stepsize and the difference of mean motions between pairs of planets. For larger stepsizes resonant perturbations are evident in the variation of the energy of the system corresponding to these stepsize resonances. It is shown that the principal instability of the method can be predicted and corresponds to the overlap of the stepsize resonances. It is noted that the analysis suggests that other artifacts will occur. For example, the overlap of a stepsize resonance with a resonance of the actual system may also give a region of chaotic behavior that is an artifact. It is pointed out that the fact that the principal artifacts corresponds to a particular set of stepsize resonances suggests that it may be possible to perturbatively remove the effect when the stepsize resonances are nonoverlapping.

Wisdom, Jack↗

NASA NTRS1990

Instabilities in line-driven stellar winds. IV - Linear perturbations in three dimensions

Nonradial wave propagation in line-driven stellar winds is analyzed including both finite disk effects and the line-drag effect of scattered line radiation. Within the local (WKB) analysis the results apply to wavelengths both longer and shorter than the Sobolev length. The finite disk causes short-wavelength waves with lateral velocity polarization to be unstable in the idealized case of pure absorption; however, the growth rates are smaller than for radially polarized waves, and the instability is damped by a relatively small amount of scattering. Hence, in realistic stellar winds, where the driving is primarily by scattering lines, perturbations with an arbitrary mixture of lateral and radial polarizations at the wind base should quickly become nearly radially polarized farther out in the wind. The implications of these results are discussed, both for interpretation of observational signatures of wind structure and for theoretical calculations aimed at modeling the nonlinear evolution of wind instabilities.

Rybicki, G. B.↗

NASA NTRS2004

Kinematid Parameters of Corrective Postural Responses Differ between Upper and Lower Body Perturbations

Balance control is disrupted following prolonged microgravity exposure, and to better understand this, both upper and lower body perturbations have been used to study postural control in space flight crewmembers. However, differences between several postural response indicators observed using the two techniques suggest that different sensory systems may be involved in organizing responses to these different perturbation approaches. The present study sought to compare differences in parameters of corrective postural responses between upper body perturbations (pushes to the chest) and forward translations of the support surface. Nine subjects participated in this study. Forward translations were performed using a NeuroCom EquiTest(TM) CDP system, which was synchronized with a Northern Digital OptoTrak motion tracking system (3 subjects). Chest pushes were applied using a hand-held force transducer device and were performed using a stabilometric system (6 subjects). Analysis of EMG has shown that: i) the earliest response of the leg muscles was registered significantly later during forward translation of the support surface than during chest pushes, and ii) there was a tendency for the different order of leg muscles activation during the translation tests. Analysis of the kinematic data showed a significant difference in the subject's body segments inclinations during corrective postural responses to upper and lower body perturbations. It appears that upper body perturbations likely engage the vestibular system more rapidly, while lower body perturbations likely engage somatosensory systems more rapidly. These differences must be taken into account when choosing the type of perturbation for testing postural function.

Sayenko, G.↗

DOE OSTI2022

Collisions of localized shocks and quantum circuits

We study collisions between localized shockwaves inside a black hole interior. We give a holographic boundary description of this process in terms of the overlap of two growing perturbations in a shared quantum circuit. The perturbations grow both exponentially as well as ballistically. Due to a competition between different physical effects, the circuit analysis shows dependence on the transverse locations and exhibits four regimes of qualitatively different behaviors. On the gravity side we study properties of the post-collision geometry, using exact calculations in simple setups and estimations in more general circumstances. We show that the circuit analysis offers intuitive and surprisingly accurate predictions about gravity computations involving non-linear features of general relativity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Planck 2018 results. VI. Cosmological parameters

In this paper, we present cosmological parameter results from the final full-mission Planck measurements of the cosmic microwave background (CMB) anisotropies, combining information from the temperature and polarization maps and the lensing reconstruction. Compared to the 2015 results, improved measurements of large-scale polarization allow the reionization optical depth to be measured with higher precision, leading to significant gains in the precision of other correlated parameters. Improved modelling of the small-scale polarization leads to more robust constraints on many parameters, with residual modelling uncertainties estimated to affect them only at the 0.5 σ level. We find good consistency with the standard spatially-flat 6-parameter ΛCDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density Ω c h 2 = 0.120 ± 0.001, baryon density Ω b h 2 = 0.0224 ± 0.0001, scalar spectral index n s = 0.965 ± 0.004, and optical depth τ = 0.054 ± 0.007 (in this abstract we quote 68% confidence regions on measured parameters and 95% on upper limits). The angular acoustic scale is measured to 0.03% precision, with 100 θ * = 1.0411 ± 0.0003. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-ΛCDM cosmology, the inferred (model-dependent) late-Universe parameters are: Hubble constant H 0 = (67.4 ± 0.5) km s –1 Mpc –1 ; matter density parameter Ω m = 0.315 ± 0.007; and matter fluctuation amplitude σ 8 = 0.811 ± 0.006. We find no compelling evidence for extensions to the base-ΛCDM model. Combining with baryon acoustic oscillation (BAO) measurements (and considering single-parameter extensions) we constrain the effective extra relativistic degrees of freedom to be N eff = 2.99 ± 0.17, in agreement with the Standard Model prediction N eff = 3.046, and find that the neutrino mass is tightly constrained to Σ m ν < 0.12 eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base ΛCDM at over 2 σ , which pulls some parameters that affect the lensing amplitude away from the ΛCDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. The joint constraint with BAO measurements on spatial curvature is consistent with a flat universe, Ω K = 0.001 ± 0.002. Also combining with Type Ia supernovae (SNe), the dark-energy equation of state parameter is measured to be w 0 = –1.03 ± 0.03, consistent with a cosmological constant. We find no evidence for deviations from a purely power-law primordial spectrum, and combining with data from BAO, BICEP2, and Keck Array data, we place a limit on the tensor-to-scalar ratio r 0.002 < 0.06. Standard big-bang nucleosynthesis predictions for the helium and deuterium abundances for the base-ΛCDM cosmology are in excellent agreement with observations. The Planck base-ΛCDM results are in good agreement with BAO, SNe, and some galaxy lensing observations, but in slight tension with the Dark Energy Survey’s combined-probe results including galaxy clustering (which prefers lower fluctuation amplitudes or matter density parameters), and in significant, 3.6 σ , tension with local measurements of the Hubble constant (which prefer a higher value). Simple model extensions that can partially resolve these tensions are not favoured by the Planck data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

A transcriptomic atlas of acute stress response to low pH in multiple Issatchenkia orientalis strains

Because of its natural stress tolerance to low pH, Issatchenkia orientalis (a.k.a. Pichia kudriavzevii) is a promising non-model yeast for bio-based production of organic acids. Yet, this organism is relatively unstudied, and specific mechanisms of its tolerance to low pH are poorly understood, limiting commercial use. In this study, we selected 12 I. orientalis strains with varying acid stress tolerance (six tolerant and six susceptible) and profiled their transcriptomes in different pH conditions to study potential mechanisms of pH tolerance in this species. We identified hundreds of genes whose expression response is shared by tolerant strains but not by susceptible strains, or vice versa, as well as genes whose responses are reversed between tolerant and susceptible strains. We mapped regulatory mechanisms of transcriptomic responses via motif analysis as well as differential network reconstruction, identifying several transcription factors, including Stb5, Mac1, and Rtg1/Rtg3, some of which are known for their roles in acid response in Saccharomyces cerevisiae. Functional genomics analysis of short-listed genes and transcription factors suggested significant roles for energy metabolism and translation-related processes, as well as the cell wall integrity pathway and RTG-dependent retrograde signaling pathway. Finally, we conducted additional experiments for two organic acids, 3-hydroxypropionate and citramalate, to eliminate acid-specific effects and found potential roles for glycolysis and trehalose biosynthesis specifically for response to low pH. In summary, our approach of comparative transcriptomics and phenotypic contrasting, along with a multi-pronged bioinformatics analysis, suggests specific mechanisms of tolerance to low pH in I. orientalis that merit further validation through experimental perturbation and engineering.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Data for A Transcriptomic Atlas of Acute Stress Response to Low pH in Multiple Issatchenkia orientalis Strains

Because of its natural stress tolerance to low pH, Issatchenkia orientalis (a.k.a. Pichia kudriavzevii ) is a promising non-model yeast for bio-based production of organic acids. Yet, this organism is relatively unstudied, and specific mechanisms of its tolerance to low pH are poorly understood, limiting commercial use. In this study, we selected 12 I. orientalis strains with varying acid stress tolerance (six tolerant and six susceptible) and profiled their transcriptomes in different pH conditions to study potential mechanisms of pH tolerance in this species. We identified hundreds of genes whose expression response is shared by tolerant strains but not by susceptible strains, or vice versa, as well as genes whose responses are reversed between tolerant and susceptible strains. We mapped regulatory mechanisms of transcriptomic responses via motif analysis as well as differential network reconstruction, identifying several transcription factors, including Stb5, Mac1, and Rtg1/Rtg3, some of which are known for their roles in acid response in Saccharomyces cerevisiae . Functional genomics analysis of short-listed genes and transcription factors suggested significant roles for energy metabolism and translation-related processes, as well as the cell wall integrity pathway and RTG-dependent retrograde signaling pathway. Finally, we conducted additional experiments for two organic acids, 3-hydroxypropionate and citramalate, to eliminate acid-specific effects and found potential roles for glycolysis and trehalose biosynthesis specifically for response to low pH. In summary, our approach of comparative transcriptomics and phenotypic contrasting, along with a multi-pronged bioinformatics analysis, suggests specific mechanisms of tolerance to low pH in I. orientalis that merit further validation through experimental perturbation and engineering.

Conversion↗

DOE OSTI2025

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2007

Earth Global Reference Atmospheric Model (GRAM99): Short Course

Earth-GRAM is a FORTRAN software package that can run on a variety of platforms including PC's. For any time and location in the Earth's atmosphere, Earth-GRAM provides values of atmospheric quantities such as temperature, pressure, density, winds, constituents, etc.. Dispersions (perturbations) of these parameters are also provided and have realistic correlations, means, and variances - useful for Monte Carlo analysis. Earth-GRAM is driven by observations including a tropospheric database available from the National Climatic Data Center. Although Earth-GRAM can be run in a "stand-alone" mode, many users incorporate it into their trajectory codes. The source code is distributed free-of-charge to eligible recipients.

Leslie, Fred W.↗

NASA NTRS2013

Stability of Menisci in Detached Bridgman Growth

Detached growth, also referred to as dewetted growth, is a Bridgman crystal growth process in which the melt is in contact with the crucible wall but the crystal is not. A meniscus bridges the gap between the top of the crystal and the crucible wall. The meniscus shape depends on the contact angle of the melt with the crucible wall, the growth angle of the melt with respect to the solidifying crystal, the gas pressure differential, the Weber number describing the rotation rate of the crucible, and the Bond number. Only some of the meniscus shapes are stable and the stability criterion is the sign of the second variation of the potential energy upon admissible meniscus shape perturbations. The effects of confined gas volumes above and below the melt and crucible rotation are evaluated. The analysis is applicable to the non-stationary case where the crystal radius changes during growth. Static stability maps (crystal radius versus pressure differential) are obtained for a series of Bond numbers, growth angles and Weber numbers. Also, the specific cases of Ge and InSb, in both terrestrial and microgravity conditions, are analyzed. Stability was found to depend significantly on whether the interior surface was considered to be microscopically rough or smooth, corresponding to pinned or unpinned states. It was also found that all meniscus shapes are statically stable in a microgravity environment.

Mazuruk, K.↗

NASA NTRS2013

Stability of Menisci in Detached Bridgman Growth

Mazuruk, Konstantin↗

NASA NTRS2021

Hypersonic Boundary-Layer Instabilities over Ogive-Cylinder Models

Computational investigations of an ogive-cylinder geometry with varying nosetips at zero degrees angle of attack are presented. The model geometry and conditions are selected to match experiments conducted in the Air Force Research Laboratory (AFRL) Mach 6 Ludwieg Tube. Five nosetips of interest were selected for the computational studies herein: two sharp ogives, two blunt ogives, and one hemispherical nosetip. Computations are performed at a freestream Reynolds number of 7.01 × 10 6 m -1 . Each sharp and blunt tip ogives had a 14 and 28 degree version for the tip angles. A cylindrical section follows the nosetip, resulting in a meter long model such as in the experiments. The laminar flow solutions are analyzed. The boundary-layer-edge properties and the velocity and temperature profiles are compared across streamwise locations aft of the ogive-cylinder junction. The blunt nosetips induce an entropy layer that envelopes the boundary-layer profiles. Modal stability analysis identifies most amplified frequencies corresponding to Mack’s second modes that agree with experimental results for the sharp tips. Similar to the experimental measurements based on wall-mounted pressure sensors, no unstable modes are found for the blunt models. Nonmodal analysis revealed a broadband set of disturbances present for the blunter tips, in agreement with experimental observations. Flow perturbation contours of most amplified planar and oblique disturbances are shown to qualitatively match wind tunnel Schlieren images, with a switch from rope-like to elongated structures, i.e., from high frequency Mack’s second modes to low frequency Mack’s first modes, as the nosetip angle is increased for the sharp tip, and from boundary-layer to entropy-layer disturbances as the bluntness is increased.

Boundary layer transition↗