Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix theories”

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 109 records · Page 6

Fluctuations in Hill’s equation parameters and application to cosmic reheating

Cosmic inflation provides a compelling framework for explaining several observed features of our Universe, but its viability depends on an efficient reheating phase that converts the inflaton’s energy into Standard Model particles. This conversion often proceeds through nonperturbative mechanisms such as parametric resonance, which is described by Hill’s equation. In this work, we investigate how stochastic fluctuations in the parameters of Hill’s equation can influence particle production during reheating. We show that such fluctuations can arise from couplings to light scalar fields and can significantly alter the stability bands in the resonance structure, thereby enhancing the growth of fluctuations and broadening the region of efficient energy transfer. Using random matrix theory and stochastic differential equations, we decompose the particle growth rate into deterministic and noise-induced components and demonstrate analytically and numerically that even modest noise leads to substantial particle production in otherwise stable regimes. Furthermore, these results suggest that stochastic effects can robustly enhance the efficacy of reheating across a wide swath of parameter space, with implications for early Universe cosmology, UV completions involving multiple scalar fields, and the resolution of the cosmological moduli problem.

Cosmology↗

Simulation of enhanced light extraction from periodic, disordered, and quasi-periodic OLED structures

Here we have developed a rigorous scattering matrix theory of light emission from periodically structured media using a Green’s function approach. We computationally simulate the spectral power inside the structure, incorporating Purcell factor enhancements, and find internal waveguiding modes and plasmonic losses. We simulate the light-outcoupling factor η out and describe how corrugations and structured media can enhance η out . We have extended our framework to describe non-periodic disordered arrays as well as aperiodic quasi-crystalline arrays. Flat organic light-emitting diodes (OLEDs) have low outcoupling with η out ~20% since most of the light is trapped in waveguided modes in higher index layers or lost to plasmonic excitations at the metal cathode. Periodically corrugated OLEDs can achieve highly enhanced η out ~65%-70% for pitch values between 1000 and 2000 nm, representing an enhancement factor >3 over planar structures. Periodic corrugations strongly diffract trapped waveguided and plasmonic modes to the emissive air cone. Disordered templates also can cause significant enhanced outcoupling with η out ~50%-55% for smaller nearest-neighbor separations of 300–400 nm. Quasi-crystalline tilings can also lead to enhancements of η out ~50%-55%. This framework can be utilized to design novel structured media that can generate high light extraction.

42 ENGINEERING↗

Machine learning applied to classifying neutron resonances

The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit resonance structure whose shape is determined in part by the angular momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is therefore of paramount importance. In this project, we apply a machine learning technique, namely decision trees, to automate the quantum number assignments. The tree is trained from simulated data generated to mimic the errors found in real data. We explore the use of several physics-motivated features for training our trees, including the nearest neighbor spacing distribution, cumulative level distribution, and channel width distributions. Initial results using random matrix theory motivated fits which demonstrated that we can determine resonance spin groups somewhat reliably. If we use these fits as features in our trees, we can train them to spot outliers corresponding to misassigned resonances. We found that with the large number of features used in this project that the decision tree tended to over t training data resulting in poor performance with respect to the test data. By reducing the number of features, we can achieve nearly perfect assignment of quantum numbers with our training data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Statistical Hauser-Feshbach model calculation for (α, n) reaction cross section and discrete level population

We calculate the α-particle induced reactions on 17,18 O, 19 F, and 23 Na, in the energy range 0 ≤ E α ≤ 10 MeV, with a particular attention to the branching ratios of populated discrete levels in the (α,n) channels. Since there are too many open channels to employ the R-matrix theory, we apply the statistical Hauser-Feshbach model to calculate these reaction cross sections and branching ratios. The branching ratio is defined as b n = $\frac{σ_n}{Σ_iσ_i + σ_c}$ , where σ n is the n-th level production cross section after the neutron emission, and σ c is the production of the continuum state. In the cases of our target nuclei and the energy range of interest, the residual nuclei of the neutron emission channel are always in their discrete states, so that σ c can be negligible.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

EDA R-Matrix Evaluations Report on n+ 9 Be, n+ 16 O [Slides]

This presentation discusses new developments in light-element R-matrix evaluation and recent light-element work. The presentation also reviewed R-matrix theory, EDAf90 code implementation, including interface EDAf90 with NJOY, ENDFtk, EXFOR/CSISRS, and a new evaluation for n+ 9 Be. The presentation concludes by stating that further development is needed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

R-matrix Resolved Resonance Region Evaluation of 140,142 Ce

Oak Ridge National Laboratory completed the resolved resonance region (RRR) evaluation of the two most abundant cerium isotopes, 140 Ce (88.45%) and 142 Ce (11.11%), as requested by the US Nuclear Criticality Safety Program. These evaluations are based on recent high-resolution transmission and capture measurements performed on nat Ce and highly enriched 142 Ce samples at the JRC-Geel Linear Accelerator facility, as well as measured thermal constants available from the EXFOR database. Starting from the resonance parameters of the ENDF/B-VIII.0 library followed by a preliminary R-matrix analysis, an updated set of resonance parameters and corresponding covariance information were derived by fitting these measured data using the Reich–Moore approximation of the R -matrix theory, as implemented in the SAMMY code system. The 140 Ce RRR upper energy limit was kept at 200 keV, whereas the 142 Ce resonance region was extended from 13 to 26 keV. Updated statistical properties were obtained for the new evaluations and compared to those derived from the ENDF/B-VIII.0 nuclear data library. The new evaluation work improved some of the discrepancies found in previous work, such as the capture resonance integral and stellar Maxwellian-averaged cross sections. These integral quantities were mainly derived from the fit of the latest measured data, especially the neutron capture yield data for 142 Ce isotope.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Potential absence of observed π2 linear-chain structures in 14O via 10C(α,α) resonances scattering

Background: The preference for light nuclear systems to coagulate into -particle clusters has been well-studied. The possibility of a linear chain configuration of -particles would allow for a new way to study this phenomenon. Purpose: A rotational band of states in C has been claimed showing a linear chain structure. The mirror system, O, has been studied here to examine how this linear chain structure is affected by replacing the valence neutrons with protons. Method: A beam of C was incident into a chamber filled with He:CO gas with the tracks recorded inside the TexAT Time Projection Chamber and the recoil -particles detected by a silicon detector array to measure the cross section. Results: The experimental cross section was compared with previous studies and fit using R-Matrix theory with the previously-observed O states being transformed to the C using mirror symmetry. The measured cross section does not replicate the claimed states, with the predicted cross section exceeding that observed at several energies and angles. Conclusion: A series of possibilities are highlighted with the most likely being that the originally-seen C states did not constitute a rotational band with a potentially incorrect spin assignment due to the limitations of the angular correlation method with non-zero spin particles. The work highlights the difficulties in measuring broad resonances corresponding to a linear chain state in a high level density.

Bishop, J.↗

Robustness tests utilizing the structure of modelling error

The present investigation is essentially concerned with the extension of results presented by Lehtomaki et al. (1981) on the robustness of multivariable linear time invariant feedback control systems. The work reported by Lehtomaki et al. is based on a multivariable version of Nyquist's theorem from which several robustness theorems were derived. In connection with the current investigation a slightly more general approach based on Nyquist's theorem is given in a fundamental robustness theorem from which various robustness tests may be obtained. A fundamental characterization of robustness is considered, and important tools from matrix theory are introduced. Attention is given to robustness tests and unstructured model error, and a robustness analysis for linear systems with structured model error.

Lehtomaki, N. A.↗

Interaction potential between a helium atom and metal surfaces

By employing an S-matrix theory for evanescent waves, the repulsive potential between a helium atom and corrugated metal surfaces has been calculated. P-wave interactions and intra-atomic correlation effects were found to be very important. The corrugation part of the interaction potential is much weaker than predicted by the effective-medium theory. Application to Cu, Ni, and Ag (110) surfaces gives good agreement with experiment without any adjustable parameters.

Takada, Y.↗

Stability sensitivity analysis for the aeroelastic optimization of a helicopter rotor

A sensitivity study of blade stability in forward flight for a hingeless rotor with respect to structural design variables is carried out using a direct analytical method. Structural design variables include nonstructural mass distribution (spanwise and chordwise), chordwise offset of center of gravity, and blade bending stiffnesses (flap, lag and torsion). The formulation for blade steady response is based on a finite element method in space and time. The vehicle trim and blade steady response are calculated iteratively as one coupled solution using a modified Newton method. Eigenvalues corresponding to different blade modes are calculated using Floquet transition matrix theory. The formulation for derivatives of the eigenvalues with respect to design variables is implemented using a direct analytical approach (chain rule differentiation), and constitutes an integral part of the regular stability analysis. The stability sensitivity derivatives were obtained at a fraction of computation time compared to the frequently adopted finite difference method. A parametric study showed that nonstructural mass and chordwise cg offset of outboarad elements, and lag bending stiffness of inboard elements, have powerful influence on blade stability.

Lim, Joon W.↗

Bearingless rotor aeromechanical stability measurements and correlations using nonlinear aerodynamics

The aeromechanical stability of a 1/8th Froude scaled bearingless rotor model was investigated experimentally in a wind tunnel. Both shaft-fixed and shaft-free conditions were examined to study the aeroelastic stability of a bearingless rotor without the incorporation of auxiliary dampers. This wind tunnel investigation generated a set of stability data for four different advance ratios, and a wide range of collective pitch settings. Theoretical analysis was performed using the newly developed University of Maryland Advanced Rotorcraft Code (UMARC). For analysis, the blade is modeled as an elastic beam undergoing flap bending, lag bending, elastic twist, and axial deformation. Blade response is calculated using a finite element method in time. Nonlinear aerodynamic effects are included by using a semiempirical stall modeling. The linearized periodic rotor perturbation equations in the nonrotating frame are solved for stability roots using Floquet transition matrix theory, as well as constant coefficient approximation. The predicted results are compared with the experimental data.

Wang, James M.↗

Numerical analysis of right-half plane zeros for a single-link manipulator

The purpose of this research is to further develop an understanding of how nonminimum phase zero location is affected by structural link design. As the demand for light-weight robots that can operate in a large workspace increases, the structural flexibility of the links become more of an issue in controls problems. When the objective is to accurately position the tip while the robot is actuated at the base, the system is nonminimum phase. One important characteristic of nonminimum phase systems is system zeros in the right half of the Laplace plane. The ability to pick the location of these nonminimum phase zeros would give the designer a new freedom similar to pole placement. The research targets a single-link manipulator operating in the horizontal plane and modeled as a Euler-Bernoulli beam with pinned-free end conditions. Using transfer matrix theory, one can consider link designs that have variable cross-sections along the length of the beam. A FORTRAN program was developed to determine the location of poles and zeros given the system model. The program was used to confirm previous research on nonminimum phase systems, and develop a relationship for designing linearly tapered links. The method allows the designer to choose the location of the first pole and zero and then defines the appropriate taper to match the desired locations. With the pole and zero location fixes, the designer can independently change the link's moment of inertia about its axis of rotation by adjusting the height of the beam. These results can be applied to inverse dynamic algorithms currently under development at Georgia Tech.

Girvin, Douglas Lynn↗

Numerical analysis of nonminimum phase zero for nonuniform link design

As the demand for light-weight robots that can operate in a large workspace increases, the structural flexibility of the links becomes more of an issue in control. When the objective is to accurately position the tip while the robot is actuated at the base, the system is nonminimum phase. One important characteristic of nonminimum phase systems is system zeros in the right half of the Laplace plane. The ability to pick the location of these nonminimum phase zeros would give the designer a new freedom similar to pole placement. This research targets a single-link manipulator operating in the horizontal plane and modeled as a Euler-Bernoulli beam with pinned-free end conditions. Using transfer matrix theory, one can consider link designs that have variable cross-sections along the length of the beam. A FORTRAN program was developed to determine the location of poles and zeros given the system model. The program was used to confirm previous research on nonminimum phase systems, and develop a relationship for designing linearly tapered links. The method allows the designer to choose the location of the first pole and zero and then defines the appropriate taper to match the desired locations. With the pole and zero location fixed, the designer can independently change the link's moment of inertia about its axis of rotation by adjusting the height of the beam. These results can be applied to the inverse dynamic algorithms that are currently under development.

Girvin, Douglas L.↗

Dynamic Reconstruction and Multivariable Control for Force-Actuated, Thin Facesheet Adaptive Optics

The Multiple Mirror Telescope (MMT) under development at the University of Arizona takes a new approach in adaptive optics placing a large (0.65 m) force-actuated, thin facesheet deformable mirror at the secondary of an astronomical telescope, thus reducing the effects of emissivity which are important in IR astronomy. However, The large size of the mirror and low stiffness actuators used drive the natural frequencies of the mirror down into the bandwidth of the atmospheric distortion. Conventional adaptive optics takes a quasi-static approach to controlling the, deformable mirror. However, flexibility within the control bandwidth calls for a new approach to adaptive optics. Dynamic influence functions are used to characterize the influence of each actuator on the surface of the deformable mirror. A linearized model of atmospheric distortion is combined with dynamic influence functions to produce a dynamic reconstructor. This dynamic reconstructor is recognized as an optimal control problem. Solving the optimal control problem for a system with hundreds of actuators and sensors is formidable. Exploiting the circularly symmetric geometry of the mirror, and a suitable model of atmospheric distortion, the control problem is divided into a number of smaller decoupled control problems using circulant matrix theory. A hierarchic control scheme which seeks to emulate the quasi-static control approach that is generally used in adaptive optics is compared to the proposed dynamic reconstruction technique. Although dynamic reconstruction requires somewhat more computational power to implement, it achieves better performance with less power usage, and is less sensitive than the hierarchic technique.

Grocott, Simon C. O.↗