Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORTHOGONAL FUNCTION”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Effect of Multi-Axial Loading on Residual Strain Tensor for 12L14 Steel Alloy

Evaluating the state of residual strain or stress is critically important for structural materials and for reliable design of complex shape components that need to function in extreme environment subjected to large thermo-mechanical loading. When residual stress state is superposed to external loads, it can lead to reduction or increase in failure strength. Past diffraction studies for evaluating the residual strain state involved measuring lattice spacings in three orthogonal directions and do not often correspond to principal directions. To completely resolve the state of strain at a given location, a full strain tensor must be determined. This is especially important when characterizing materials or metallic components exposed to biaxial or complex loading. Neutron diffraction at the second Generation Neutron Residual Stress Facility (NRSF2) at Oak Ridge National Laboratory is used in this study to measure strain tensors associated with different modes of stress path. Hollow cylinder steel samples with 2 mm wall thickness are subjected to either pure axial extension or pure torsion to simulate multi-axial loading conditions. A virgin sample that is not subjected to any deformation, but subjected to identical manufacturing conditions and machining steps involved to obtain hollow cylinder geometry is used for obtaining reference d-spacing for given hkl planes at target spatial location(s). The two samples which are subjected to either pure tension or torsion are loaded to a deformation state that corresponded to equal amount of octahedral shear strain which is an invariant. This procedure is used so that a basis for comparison between the two samples can be made to isolate the stress path effects. A 2-circle Huber orienter is used to obtain strain measurements on identical gauge volume at a series of chi and psi values. The residual state of stress tensor corresponding to ex situ (upon unloading) conditions is presented for three lattice planes (211, 110, 200) for a bcc ferritic system exposed to tension and pure torsion.

HB-2B↗
DOE OSTI2017

Hybrid grid/basis set discretizations of the Schrödinger equation

We present a new kind of basis function for discretizing the Schrödinger equation in electronic structure calculations, called a gausslet, which has wavelet-like features but is composed of a sum of Gaussians. Gausslets are placed on a grid and combine advantages of both grid and basis set approaches. They are orthogonal, infinitely smooth, symmetric, polynomially complete, and with a high degree of locality. Because they are formed from Gaussians, they are easily combined with traditional atom-centered Gaussian bases. As a result, we also introduce diagonal approximations that dramatically reduce the computational scaling of two-electron Coulomb terms in the Hamiltonian.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗
NASA NTRS1988

Modal identification using single-mode projection filters and comparison with ERA and MLE results

The Single-Mode Projection Filter (SPF) is a newly developed algorithm for eigensystem parameter identification from both analytical results and test data. The SPF is formulated with a single mode only and practical for parallel processing implementation. Explicit formulations of SPF are derived for the multi-input multi-output (MIMO) system by using the orthogonal matrices of the controllability and observability matrices in the general sense. The modal parameters of SPF are initially obtained from an analytical model in modal space. The experimental data are then processed through SPF to update its modal parameters and to minimize a cost function defined by the norm of an error matrix. The updated modal parameters represent the characteristics of the test data. A two-dimensional global minimum optimization algorithm is developed and applied for the filter update by using the interval analysis method. The SPF is developed based on a single-mode subsystem and identifies only one modal frequency and one modal damping within a specified region. For an n-modes structure, n SPF can be implemented for parallel processing to reduce the computational burden. The SPF is applied to analyze the simulated data for the MAST beam structure. The estimated modal parameters are comparable to those from the Eigensystem Realization Algorithm (ERA) and repeated modal frequencies are identified. The modal analysis of the Spacecraft Control Laboratory Experiment (SCOLE) data is also performed by using the ERA and the Maximum Likelihood Estimate (MLE). The result shows that the first five modal frequencies are very close from ERA and MLE. However, there are slight disparities in the damping rates and the computational burdens are quite different among these two algorithms.

Huang, Jen-Kuang↗
NASA NTRS1996

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗
NASA NTRS2003

Aeroelastic System Development Using Proper Orthogonal Decomposition and Volterra Theory

This research combines Volterra theory and proper orthogonal decomposition (POD) into a hybrid methodology for reduced-order modeling of aeroelastic systems. The out-come of the method is a set of linear ordinary differential equations (ODEs) describing the modal amplitudes associated with both the structural modes and the POD basis functions for the uid. For this research, the structural modes are sine waves of varying frequency, and the Volterra-POD approach is applied to the fluid dynamics equations. The structural modes are treated as forcing terms which are impulsed as part of the uid model realization. Using this approach, structural and uid operators are coupled into a single aeroelastic operator. This coupling converts a free boundary uid problem into an initial value problem, while preserving the parameter (or parameters) of interest for sensitivity analysis. The approach is applied to an elastic panel in supersonic cross ow. The hybrid Volterra-POD approach provides a low-order uid model in state-space form. The linear uid model is tightly coupled with a nonlinear panel model using an implicit integration scheme. The resulting aeroelastic model provides correct limit-cycle oscillation prediction over a wide range of panel dynamic pressure values. Time integration of the reduced-order aeroelastic model is four orders of magnitude faster than the high-order solution procedure developed for this research using traditional uid and structural solvers.

Lucia, David J.↗
NASA NTRS2010

Calibration Designs for Non-Monolithic Wind Tunnel Force Balances

This research paper investigates current experimental designs and regression models for calibrating internal wind tunnel force balances of non-monolithic design. Such calibration methods are necessary for this class of balance because it has an electrical response that is dependent upon the sign of the applied forces and moments. This dependency gives rise to discontinuities in the response surfaces that are not easily modeled using traditional response surface methodologies. An analysis of current recommended calibration models is shown to lead to correlated response model terms. Alternative modeling methods are explored which feature orthogonal or near-orthogonal terms.

Wind tunnel models↗
DOE OSTI2013

MAGNETIC NANOPARTICLES IN THE INTERSTELLAR MEDIUM: EMISSION SPECTRUM AND POLARIZATION

The presence of ferromagnetic or ferrimagnetic nanoparticles in the interstellar medium would give rise to magnetic dipole radiation at microwave and submillimeter frequencies. Such grains may account for the strong millimeter-wavelength emission observed from a number of low-metallicity galaxies, including the Small Magellanic Cloud. We calculate the absorption and scattering cross sections for such grains, with particular attention to metallic Fe, magnetite Fe{sub 3}O{sub 4}, and maghemite {gamma}-Fe{sub 2}O{sub 3}, all potentially present in the interstellar medium. The rate of Davis-Greenstein alignment by magnetic dissipation is also estimated. We determine the temperature of free-flying magnetic grains heated by starlight and calculate the polarization of the magnetic dipole emission from both free-fliers and inclusions. For inclusions, the magnetic dipole emission is expected to be polarized orthogonally relative to the normal electric dipole radiation. Magnetic dipole radiation will contribute significantly to the 20-40 GHz anomalous microwave emission only if a large fraction of the Fe is in metallic Fe iron nanoparticles with extreme elongations. Finally, we present self-consistent dielectric functions for metallic Fe, magnetite Fe{sub 3}O{sub 4}, and maghemite {gamma}-Fe{sub 2}O{sub 3}, enabling calculation of absorption and scattering cross sections from microwave to X-ray wavelengths.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A prototype, glassless densitometer traceable to primary optical standards for quantitative radiochromic film dosimetry

Purpose: To evaluate a prototype densitometer traceable to primary optical standards and compare its performance to an EPSON Expression{sup ®} 10000XL flatbed scanner (the Epson) for quantitative radiochromic film (RCF) dosimetry. Methods: A prototype traceable laser densitometry system (LDS) was developed to mitigate common film scanning artifacts, such as positional scan dependence and high noise in low-dose regions, by performing point-based measurements of RCF suspended in free-space using coherent light. The LDS and the Epson optical absorbance scales were calibrated up to 3 AU, using reference materials calibrated at a primary standards laboratory and a scanner calibration factor (SCF). Calibrated optical density (OD) was determined for 96 Gafchromic{sup ®} EBT3 film segments before and after irradiation to one of 16 dose levels between 0 and 10 Gy, exposed to {sup 60}Co in a polymethyl-methacrylate (PMMA) phantom. The sensitivity was determined at each dose level and at two rotationally orthogonal readout orientations to obtain the sensitometric response of each RCF dosimetry system. LDS rotational scanning dependence was measured at nine angles between 0°and 180°, due to the expected interference between coherent light and polarizing EBT3 material. The response curves were fit to the analytic functions predicted by two physical response models: the two-parameter single-hit model and the four-parameter percolation model. Results: The LDS and the Epson absorbance measurements were linear to primary optical standards to within 0.2% and 0.3% up to 2 and 1 AU, respectively. At higher densities, the LDS had an over-response (2.5% at 3 AU) and the Epson an under-response (3.1% and 9.8% at 2 and 3 AU, respectively). The LDS and the Epson SCF over the applicable range were 0.968% ± 0.2% and 1.561% ± 0.3%, respectively. The positional scan dependence was evaluated on each digitizer and shown to be mitigated on the LDS, as compared to the Epson. Maximum EBT3 rotational dependence was found to have a strong dependence on dose (0.1% and 34% at 30 mGy and 5 Gy, respectively). The preferred EBT3 polymerization axis angle was constant within experimental uncertainties. In its most sensitive orientation, the LDS-measured EBT3 sensitivity was 7.13 × 10{sup −4} ± 9.2 × 10{sup −6} AU/mGy, which represented a 4.5 fold increase over the Epson of 1.58 × 10{sup −4} ± 9.8 × 10{sup −6} AU/mGy. To first order approximations, EBT3 response was linear up to 500 mGy to within 0.80% and to within 7.5% for the most sensitive LDS and the Epson orientations, respectively. The corresponding single-hit and percolation model relative residual norms were 0.082 and 0.074 for LDS as compared to 0.29 and 0.18 for the Epson, which represented a significant increase in LDS-measured agreement with the simple physical model. Less sensitive LDS and the Epson orientations showed a marked decrease in the physical model agreement, which suggested that suboptimal readout device characteristics may be the origin of the complex sensitometric functional forms currently required for accurate RCF dosimetry. Conclusions: The prototype densitometer was shown to be superior to a conventional scanner for quantitative RCF dosimetry based on physical models of film response. The Epson was shown to be a reliable tool for routine RCF dosimetry in a clinical setting, yet calibration to primary optical standards did not mitigate the necessity for complex, empirical functional form fitting.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1976

On differential transformations between Cartesian and curvilinear (geodetic) coordinates

Differential transformations are developed between Cartesian and curvilinear orthogonal coordinates. Only matrix algebra is used for the presentation of the basic concepts. After defining the reference systems used the rotation (R), metric (H), and Jacobian (J) matrices of the transformations between cartesian and curvilinear coordinate systems are introduced. A value of R as a function of H and J is presented. Likewise an analytical expression for J(-1) as a function of H(-2) and R is obtained. Emphasis is placed on showing that differential equations are equivalent to conventional similarity transformations. Scaling methods are discussed along with ellipsoidal coordinates. Differential transformations between elipsoidal and geodetic coordinates are established.

Soler, T.↗
NASA NTRS2002

Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation

The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.

Chi, X.↗
DOE OSTI2025

Imine Reductase-Catalyzed, Radical-Mediated Asymmetric Cyano Group Migration

Functional group migration (FGM) reactions represent a fundamental class of transformations in organic chemistry, enabling the repositioning of functional moieties in nonobvious ways. However, catalytic asymmetric radical-mediated FGMs remain rare due to the inherent challenges of achieving catalyst-controlled enantioselectivity over free radical intermediates. Herein, we repurpose imine reductases (IREDs), a class of biotechnologically important enzymes known for their substrate promiscuity, to enable the first examples of catalytic asymmetric cyano group migration via a radical mechanism. An orthogonal set of radical enzymes, including PbaIREDCym and SmiIREDCym, was engineered, allowing both 1,4- and 1,5-cyano group migration reactions to occur in an enantiodivergent fashion. The use of the nonionic surfactant TPGS-1000 was found to improve both the yield and enantioselectivity of these cyano migration reactions. Furthermore, this biocatalytic process exhibited a broad substrate scope and is readily scalable, affording a rare example of chiral nonamine product assembly with imine reductases. More broadly, stereoselective radical biocatalysis with engineered IREDs and other versatile enzymes provides a potentially general solution to challenging asymmetric FGM reactions.

Biocatalysis↗
DOE OSTI2023

An efficient virus–induced gene silencing (VIGS) system for gene functional studies in Miscanthus

Virus-induced gene silencing (VIGS) is a powerful tool for transient gene functional analysis in plants, especially for monocot species (e.g., grasses) that are recalcitrant to transformation. Despite various VIGS systems that have been developed in different plant species, none was previously available for the bioenergy crop Miscanthus. Here, we report the establishment of an efficient and robust VIGS system mediated by Tobacco Rattle Virus (TRV) in Miscanthus. We first investigated the impact of various factors that may affect gene silencing efficiency using the Miscanthus sinensis Phytoene Desaturase (MsPDS) gene as a visual indicator of photobleaching. Then, we optimized the TRV-elicited VIGS procedure using an orthogonal experimental design with four factors (sprout size, Agrobacterium concentration, vacuum infiltration time, and co-incubation time) each at three levels. The following led to the highest silencing efficiency (~76%): inoculation of germinating seedlings (1.0–2.0 mm), Agrobacterium tumefaciens culture grown to optical density at 600 nm (OD 600 ) of 0.4, vacuum infiltration for 90 min, and co-incubation for 5 h. The VIGS system established was applicable for both M. sinensis and M. lutarioriparius, with comparable gene silencing efficiency. We verified the efficacy of the VIGS system via the functional characterization of the role of a MYB transcription factor, MsMYB112, in salt stress tolerance. Expression of MsMYB112 was successfully knocked down using the VIGS system, and this led to compromised salt tolerance in the silenced Miscanthus plants. The TRV-based VIGS system established may, therefore, substantially facilitate functional genomic studies in Miscanthus.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

On the interaction of a wind turbine wake with a conventionally neutral atmospheric boundary layer

In this work, we investigate the dynamics of wind turbine tip-vortex breakdown in a conventionally neutral atmospheric boundary layer (ABL). To this end, high-resolution data are collected from large-eddy simulations of a wind turbine operating within a neutral ABL and studied by means of proper orthogonal decomposition (POD) and Fourier analysis. The high resolution of the generated data in both space and time allows us to gain insight into the tip-vortex breakdown mechanisms by (i) capturing the energy modes of the coherent structures, (ii) studying their contribution to the tip-vortex breakdown through their power spectra functions and mean kinetic energy (MKE) flux, and (iii) analysing the growth rate of each contributing perturbation frequency along tip vortices. Our analysis shows that under a fully turbulent scenario, the growth rate of perturbations along the tip vortices is largest for low wave numbers, i.e. long-wave perturbations. Additionally, the MKE flux reaches its highest value at two diameters downstream of the rotor plane, a behaviour that can be attributed to the coexistence of multiple interacting POD modes, with the streamwise vortex roller mode being the primary contributor to the total MKE flux budget, contributing approximately . Finally, comparisons with a laminar, uniform flow scenario subject to a single-frequency perturbation highlight the differences between the two ambient flow conditions. In the non-turbulent, uniform flow scenario, the growth rate attains its maximum value at a wave number corresponding to the out-of-phase mutual-inductance mechanism, whereas the MKE flux exhibits local minima and maxima along the wake and at different downstream locations depending on the perturbation frequency. Our analyses suggest that the breakdown of the wind turbine tip vortices under a fully turbulent neutral ABL inflow is due to complex interactions across a range of excitation frequencies, in which the mutual-inductance instability may not be the dominant one.

17 WIND ENERGY↗
DOE OSTI2020

The completed SDSS-IV extended baryon oscillation spectroscopic survey: geometry and growth from the anisotropic void–galaxy correlation function in the luminous red galaxy sample

ABSTRACT We present an analysis of the anisotropic redshift-space void–galaxy correlation in configuration space using the Sloan Digital Sky Survey extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 luminous red galaxy (LRG) sample. This sample consists of LRGs between redshifts 0.6 and 1.0, combined with the high redshift z > 0.6 tail of the Baryon Oscillation Spectroscopic Survey Data Release 12 CMASS sample. We use a reconstruction method to undo redshift-space distortion (RSD) effects from the galaxy field before applying a watershed void-finding algorithm to remove bias from the void selection. We then perform a joint fit to the multipole moments of the correlation function for the growth rate fσ8 and the geometrical distance ratio DM/DH, finding $f\sigma _8(z_\rm {eff})=0.356\pm 0.079$ and $D_M/D_H(z_\rm {eff})=0.868\pm 0.017$ at the effective redshift $z_\rm {eff}=0.69$ of the sample. The posterior parameter degeneracies are orthogonal to those from galaxy clustering analyses applied to the same data, and the constraint achieved on DM/DH is significantly tighter. In combination with the consensus galaxy BAO and full-shape analyses of the same sample, we obtain fσ8 = 0.447 ± 0.039, DM/rd = 17.48 ± 0.23, and DH/rd = 20.10 ± 0.34. These values are in good agreement with the ΛCDM model predictions and represent reductions in the uncertainties of $13{{\ \rm per\ cent}}$, $23{{\ \rm per\ cent}}$, and $28{{\ \rm per\ cent}}$, respectively, compared to the combined results from galaxy clustering, or an overall reduction of 55 per cent in the allowed volume of parameter space.

Nadathur, Seshadri↗
DOE OSTI2020

Temperature-dependent angle-resolved photoemission spectroscopy study of the Ce 4f states in a possible topological Kondo insulator CeRhAs

The electronic structure of a possible topological Kondo insulator of CeRhAs has been investigated by employing temperature (T)-dependent angle-resolved photoemission spectroscopy (ARPES). Fermi surfaces (FSs) and band structures are successfully measured for three orthogonal crystallographic planes. The measured Fermi-edge states are found to have a three-dimensional (3D) character, contradictory to the proposed topological surface states of the 2 D character. The measured FSs due to the Ce 4f electrons agree well with those from the density functional theory band structures unfolded into the reduced Ce-only unit cell. The theoretically predicted hourglass-type bulk bands along $\bar{X}\bar{S}\bar{X}$ are not clearly resolved, in spite of the evident existence of those band features. Finally, T -dependent ARPES measurements reveal that the coherent Ce 4f states, having two pseudogap structures of Δ 1 ~ 80 meV and Δ 2 ~ 30 meV, persist to remain above 200 K, in agreement with the high Kondo temperature of CeRhAs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

MO-FG-CAMPUS-TeP1-05: Rapid and Efficient 3D Dosimetry for End-To-End Patient-Specific QA of Rotational SBRT Deliveries Using a High-Resolution EPID

Purpose: EPID-based patient-specific quality assurance provides verification of the planning setup and delivery process that phantomless QA and log-file based virtual dosimetry methods cannot achieve. We present a method for EPID-based QA utilizing spatially-variant EPID response kernels that allows for direct calculation of the entrance fluence and 3D phantom dose. Methods: An EPID dosimetry system was utilized for 3D dose reconstruction in a cylindrical phantom for the purposes of end-to-end QA. Monte Carlo (MC) methods were used to generate pixel-specific point-spread functions (PSFs) characterizing the spatially non-uniform EPID portal response in the presence of phantom scatter. The spatially-variant PSFs were decomposed into spatially-invariant basis PSFs with the symmetric central-axis kernel as the primary basis kernel and off-axis representing orthogonal perturbations in pixel-space. This compact and accurate characterization enables the use of a modified Richardson-Lucy deconvolution algorithm to directly reconstruct entrance fluence from EPID images without iterative scatter subtraction. High-resolution phantom dose kernels were cogenerated in MC with the PSFs enabling direct recalculation of the resulting phantom dose by rapid forward convolution once the entrance fluence was calculated. A Delta4 QA phantom was used to validate the dose reconstructed in this approach. Results: The spatially-invariant representation of the EPID response accurately reproduced the entrance fluence with >99.5% fidelity with a simultaneous reduction of >60% in computational overhead. 3D dose for 10{sub 6} voxels was reconstructed for the entire phantom geometry. A 3D global gamma analysis demonstrated a >95% pass rate at 3%/3mm. Conclusion: Our approach demonstrates the capabilities of an EPID-based end-to-end QA methodology that is more efficient than traditional EPID dosimetry methods. Displacing the point of measurement external to the QA phantom reduces the necessary complexity of the phantom itself while offering a method that is highly scalable and inherently generalizable to rotational and trajectory based deliveries. This research was partially supported by Varian.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

RUScal : Software for the analysis of resonant ultrasound spectroscopy measurements

Resonant ultrasound spectroscopy is used to nondestructively measure the elastic resonances of small solids to elucidate the material's elastic properties or other qualities like size, shape, or composition. In this work, we introduce the software RUScal for the purpose of determining elastic properties by analyzing the eigenfrequencies of solid specimens with common shapes, such as rectangular parallelepipeds, cylinders (solid and hollow tube), ellipsoids, and octahedrons, as well as irregularly shaped ellipsoids that can be described analytically. All symmetry classes are supported, from isotropic to triclinic, along with the option to add or remove up to three orthogonal mirror planes as well as the ability to reorient the crystal axes with respect the sample edges via Euler angles. Additional features include tools to help find initial sets of elastic constants, including grid exploration and Monte Carlo methods, a tool to analyze frequencies as a function of sample length or crystal orientation, an error analysis tool to assess fit quality, and formatting of the input and output files for batch fitting, e.g., as a function of temperature. This software was validated with published resonant ultrasound spectroscopy data for various materials, shapes, and symmetries with noted improvements in calculation time compared to finite element methods.

47 OTHER INSTRUMENTATION↗