Search NASA⌕ Search

SEARCH · Search NASA

Results for “RESPONSE MATRIX METHOD”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Error Propagation for Response Matrix Methods in the Many Node Limit

The response matrix method has been a useful solution method for radiation transport problems for decades. The method is desirable for many reasons. It can be extended to both fixed source and eigenvalue problems. In addition, a key aspect of the method is a shift of computational expense to a precomputation phase. Methods and codes employing response matrix or response matrix-like methods have enjoyed formidable computational efficiency while solving large problems. However, current applications of the method rely on large node sizes and few nodes to maximize computational efficiency. As part of a desire to make these methods more amenable to applications such as parallel computing in the future, the number of nodes in problems may increase, and node size will decrease. The effect of this change has not been well studied except for a few cases. In the past, it was determined that a poor approximation of angular dependence of flux on node interfaces could lead to large error propagation for fixed source problems, but a newer study indicates that changing node sizes as well as number of nodes poses little challenge for response methods in the eigenvalue case. This paper presents the unintuitive result that for response matrix methods increasing the spatial resolution of the problem may dramatically decrease accuracy - this effect is demonstrated on a simple analytic 1-D benchmark problem. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT

We present a special symmetric Lanczos algorithm and a kernel polynomial method (KPM) for approximating the absorption spectrum of molecules within the linear response time-dependent density functional theory (TDDFT) framework in the product form. In contrast to existing algorithms, the new algorithms are based on reformulating the original non-Hermitian eigenvalue problem as a product eigenvalue problem and the observation that the product eigenvalue problem is self-adjoint with respect to an appropriately chosen inner product. This allows a simple symmetric Lanczos algorithm to be used to compute the desired absorption spectrum. The use of a symmetric Lanczos algorithm only requires half of the memory compared with the nonsymmetric variant of the Lanczos algorithm. The symmetric Lanczos algorithm is also numerically more stable than the nonsymmetric version. The KPM algorithm is also presented as a low-memory alternative to the Lanczos approach, but the algorithm may require more matrix-vector multiplications in practice. We discuss the pros and cons of these methods in terms of their accuracy as well as their computational and storage cost. Applications to a set of small and medium-sized molecules are also presented.

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

Beam manipulation techniques, nonlinear beam dynamics, and space charge effect in high energy high power accelerators

We had carried out a design of an ultimate storage ring with beam emittance less than 10 picometer for the feasibility of coherent light source at X-ray wavelength. The accelerator has an inherent small dynamic aperture. We study method to improve the dynamic aperture and collective instability for an ultimate storage ring. Beam measurement and accelerator modeling are an integral part of accelerator physics. We develop the independent component analysis (ICA) and the orbit response matrix method for improving accelerator reliability and performance. In collaboration with scientists in National Laboratories, we also carry out experimental and theoretical studies on beam dynamics. Our proposed research topics are relevant to nuclear and particle physics using high brightness particle and photon beams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Compensation of Third-order Resonances in the High Intensity Regime

As the Fermilab Accelerator Complex enters the high-intensity era, the Recycler Ring (RR) needs to mitigate the detrimental effect of third-order resonance crossing. Third-order resonance lines can be compensated to first order by cancelling out the global Resonance Driving Terms (RDTs) using the response matrix method. This compensation scheme has been proven to work at low intensities, i.e., in the single-particle regime. In order to evaluate the effectiveness of this compensation scheme at higher intensities, this study looks at dynamic and static tune scans, with and without resonance compensation, and different space charge tune shifts. Special care was taken in order to disentangle effects from space charge tune shift, structure resonances and space charge driven resonances.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

A Bayesian method to estimate the neutron response matrix of a single crystal CVD diamond detector

Detectors made from artificial chemical vapor deposition (CVD) single crystal diamond are very promising candidates for applications where high resolution neutron spectrometry in very high neutron fluxes is required, for example in fusion research. We propose a Bayesian method to estimate the neutron response function of the detector for a continuous range of neutron energies (in our case, 10 MeV ≤ E{sub n} ≤ 16 MeV) based on a few measurements with quasi-monoenergetic neutrons. This method is needed because a complete set of measurements is not available and the alternative approach of using responses based on Monte Carlo calculations is not feasible. Our approach uses Bayesian signal-background separation techniques and radial basis function interpolation methods. We present the analysis of data measured at the PTB accelerator facility PIAF. The method is quite general and it can be applied to other particle detectors with similar characteristics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1952

Matrix Method of Determining the Longitudinal-Stability Coefficients and Frequency Response of an Aircraft from Transient Flight Data

A matrix method is presented for determining the longitudinal-stability coefficients and frequency response of an aircraft from arbitrary maneuvers. The method is devised so that it can be applied to time-history measurements of combinations of such simple quantities as angle of attack, pitching velocity, load factor, elevator angle, and hinge moment to obtain the over-all coefficients. Although the method has been devised primarily for the evaluation of stability coefficients which are of primary interest in most aircraft loads and stability studies, it can be used also, with a simple additional computation, to determine the frequency-response characteristics. The entire procedure can be applied or extended to other problems which can be expressed by linear differential equations.

Donegan, James J↗
NASA NTRS2024

Analytical Simulation of Effects of Local Mechanisms and Microstructure on Creep Response of Unidirectional Ceramic Matrix Composites

A micromechanics-based method has been developed to analyze the creep response of uncoated ceramic matrix minicomposites. Although the global stress level at which the creep response is analyzed is lower than the composite proportional limit, local stresses are assumed to be high enough that localized damage is present in the composite in the form of matrix microcracks. To model the composite creep response, including the effects of matrix cracking, a fiber shear lag-based methodology is employed. In this approach, stresses are assumed to vary in the fiber as a function of time and distance from the crack plane. The varying stresses are then used to compute the overall creep strain for the composite. Various assumptions regarding the level of matrix microcracking in the composite and the level of creep in the fiber and matrix in various portions of the composite unit cell are also examined. The creep response of the fiber is modeled using a linear Burgers model. The model is applied to a SiCf/SiC unidirectional minicomposite system. The computed creep results are compared to experimentally obtained values. The effects of the local fiber volume fraction on the overall creep response of the composite are also studied. This work will allow increased understanding of the key material damage mechanisms and load sharing that take place during creep conditions and can be expanded to provide improved analysis methods for full macrocomposites.

Micromechanics↗
NASA NTRS2024

Analytical Simulation of the Effects of Local Mechanisms and Microstructure on the Creep Response of Unidirectional Ceramic Matrix Composites

A micromechanics-based method has been developed to analyze the creep response of uncoated ceramic matrix minicomposites. Although the global stress level at which the creep response is analyzed is lower than the composite proportional limit, local stresses are assumed to be high enough that localized damage is present in the composite in the form of matrix microcracks. To model the composite creep response including the effects of matrix cracking, a fiber shear lag-based methodology is employed. In this approach, stresses are assumed to vary in the fiber as a function of time and distance from the crack plane. The varying stresses are then used to compute the overall creep strain for the composite. Various assumptions regarding the level of matrix microcracking in the composite, and the level of creep in the fiber and matrix in various portions of the composite unit cell, are also examined. The creep response of the fiber is modeled using a linear Burgers model. The model is applied to a SiC f /SiC unidirectional minicomposite system. The computed creep results are compared to experimentally obtained values. The effects of the local fiber volume fraction on the overall creep response of the composite are also studied. This work will allow increased understanding of the key material damage mechanisms and load sharing that take place during creep conditions and can be expanded to provide improved analysis methods for full macrocomposites.

Micromechanics↗
DOE OSTI2020

Quantitative and three-dimensional assessment of holdup material

Nuclear material deposited in equipment, transfer lines, and ventilation systems of a processing facility is usually referred to as holdup. Operators need to know the location and amount of holdup in their facility for radiation safety and material accountability. In this work, we propose to use an array of detectors co-axial to the inspected pipe to measure the holdup material. This method is implementable into an automated system capable of crawling on surfaces and pipes of various curvatures, which would enable faster, easier, and more accurate holdup safeguards measurements. We first demonstrated that the current holdup assay procedure could lead to a non-negligible bias in the estimate of special nuclear material mass, due to the simplified assumption of deposited geometry introduced by the Generalized Geometry Holdup (GGH) model. The new approach consists of imaging the inner holdup material by characterizing the detector array’s response and unfolding it from the measured light output. Our experimental proof of principle consists of three NaI(Tl) detectors surrounding an aluminum pipe containing two cesium-137 ( 137 Cs) sources. We derived the source distribution inside the pipe by first calculating the detector response matrix using a method adaptive to the surface geometry of the object containing the measured holdup material. Creating a matrix of the detector array’s measured counts, we then proceed to solve an inverse problem, resulting in an accurately located source position and activity distribution within the response matrix’s spatial resolution. Furthermore, we then developed a simulated model of the envisioned experimental setup, which accurately described both the activity and position of the source in 2D. Finally, we extended our model onto a discretized three-dimensional model of the system, encompassing 36 detectors. The model was simulated to measure the source distribution and activity of multiple sources in pipes of varied geometries, which could then be conveniently translated to a real-time holdup measurement protocol. For the 3D simulation of four different source geometries, the model accurately localized the source position in 3D, while the activity retained a maximum relative error of ±5.32%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Interpreting plasmonic response of epitaxial Ag/Si(100) island ensembles

Associating features in the experimentally measured optical response of epitaxial Ag islands grown on Si(100) with the localized surface plasmon resonances (LSPRs) hosted by the Ag islands is challenging due to the variation of the Si dielectric function over the energy range under consideration. However, it is possible to conclusively identify features in the experimental spectra with LSPR modes oscillating both parallel and perpendicular to the epitaxial interface by simulating the optical response. The Abeles matrix method is used to describe the composite layered system and the Ag islands are modeled using the thin island film model developed by Bedeaux and Vlieger. By incorporating island morphology parameters determined by quantitative analysis of electron micrographs, the simulation faithfully reproduces the main features of the experimental spectra. Individually zeroing the dipoles associated with the LSPR modes enables conclusive identification of their contribution to the optical response of the composite system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

Coupled bending-torsion steady-state response of pretwisted, nonuniform rotating beams using a transfer-matrix method

Using the Newtonian method, the equations of motion are developed for the coupled bending-torsion steady-state response of beams rotating at constant angular velocity in a fixed plane. The resulting equations are valid to first order strain-displacement relationships for a long beam with all other nonlinear terms retained. In addition, the equations are valid for beams with the mass centroidal axis offset (eccentric) from the elastic axis, nonuniform mass and section properties, and variable twist. The solution of these coupled, nonlinear, nonhomogeneous, differential equations is obtained by modifying a Hunter linear second-order transfer-matrix solution procedure to solve the nonlinear differential equations and programming the solution for a desk-top personal computer. The modified transfer-matrix method was verified by comparing the solution for a rotating beam with a geometric, nonlinear, finite-element computer code solution; and for a simple rotating beam problem, the modified method demonstrated a significant advantage over the finite-element solution in accuracy, ease of solution, and actual computer processing time required to effect a solution.

Gray, Carl E., Jr.↗
DOE OSTI2014

Linear optical response of finite systems using multishift linear system solvers

We discuss the application of multishift linear system solvers to linear-response time-dependent density functional theory. Using this technique the complete frequency-dependent electronic density response of finite systems to an external perturbation can be calculated at the cost of a single solution of a linear system via conjugate gradients. We show that multishift time-dependent density functional theory yields excitation energies and oscillator strengths in perfect agreement with the standard diagonalization of the response matrix (Casida's method), while being computationally advantageous. We present test calculations for benzene, porphin, and chlorophyll molecules. We argue that multishift solvers may find broad applicability in the context of excited-state calculations within density-functional theory and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Utilization of the Generalized Method of Cells to Analyze the Deformation Response of Laminated Ceramic Matrix Composites

In order to practically utilize ceramic matrix composites in aircraft engine components, robust analysis tools are required that can simulate the material response in a computationally efficient manner. The MAC/GMC software developed at NASA Glenn Research Center, based on the Generalized Method of Cells micromechanics method, has the potential to meet this need. Utilizing MAC/GMC, the effective stiffness properties, proportional limit stress and ultimate strength can be predicted based on the properties and response of the individual constituents. In this paper, the effective stiffness and strength properties for a representative laminated ceramic matrix composite with a large diameter fiber are predicted for a variety of fiber orientation angles and laminate orientations. As part of the analytical study, methods to determine the in-situ stiffness and strength properties of the constituents required to appropriately simulate the effective composite response are developed. The stiffness properties of the representative composite have been adequately predicted for all of the fiber orientations and laminate configurations examined in this study. The proportional limit stresses and strains and ultimate stresses and strains were predicted with varying levels of accuracy, depending on the laminate orientation. However, for the cases where the predictions did not have the desired level of accuracy, the specific issues related to the micromechanics theory were identified which could lead to difficulties that were encountered that could be addressed in future work.

Goldberg, Robert K.↗
NASA NTRS2021

Analytical Simulation of Effects of Local Mechanisms on Tensile Response of Ceramic Matrix Minicomposites

A micromechanics-based method based on a fiber shear lag analysis was developed to analyze the fast fracture response of uncoated ceramic matrix minicomposites. The model was applied to two SiCf/SiC minicomposite systems that were fabricated by Rolls-Royce and the University of Connecticut. The analysis approach accounts in an average sense for the fact that the cracks in an actual minicomposite often have an irregular geometry. The effects of local variations in the fiber volume ratio on the composite response was also investigated. Parametric studies were performed to investigate the effect of the interfacial shear stress, global fiber volume fraction and percentage of the composite that remains uncracked on the proportional limit stress and the composite secondary modulus. This work will facilitate an increased understanding of the key material mechanisms that take place during fast fracture loading.

Micromechanics Ceramic Matrix Composites Damage↗
NASA NTRS1990

Microlevel thermal effects in metal matrix composites

A method for studying the influence of thermal effects on the inelastic response of metal matrix composites is reviewed. A micromechanics approach based upon the method of cells is shown to be quite versatile for studying a variety of materials response phenomena. Yielding and inelastic response of the composite are predicted as functions of thermal stresses, yielding of the matrix, and imperfect fiber/matrix bonding. Results are presented in the form of yield surfaces and nonlinear stress-strain curves for unidirectional and laminated boron/aluminum and silicon-carbide/titanium.

Herakovich, Carl T.↗
NASA NTRS2000

A Study of Influencing Factors on the Tensile Response of a Titanium Matrix Composite With Weak Interfacial Bonding

The generalized method of cells micromechanics model is utilized to analyze the tensile stress-strain response of a representative titanium matrix composite with weak interfacial bonding. The fiber/matrix interface is modeled through application of a displacement discontinuity between the fiber and matrix once a critical debonding stress has been exceeded. Unidirectional composites with loading parallel and perpendicular to the fibers are examined, as well as a cross-ply laminate. For each of the laminates studied, analytically obtained results are compared to experimental data. The application of residual stresses through a cool-down process was found to have a significant effect on the tensile response. For the unidirectional laminate with loading applied perpendicular to the fibers, fiber packing and fiber shape were shown to have a significant effect on the predicted tensile response. Furthermore, the interface was characterized through the use of semi-emperical parameters including an interfacial compliance and a "debond stress;" defined as the stress level across the interface which activates fiber/matrix debonding. The results in this paper demonstrate that if architectural factors are correctly accounted for and the interface is appropriately characterized, the macro-level composite behavior can be correctly predicted without modifying any of the fiber or matrix constituent properties.

Goldberg, Robert K.↗
NASA NTRS2012

Investigation of Effects of Material Architecture on the Elastic Response of a Woven Ceramic Matrix Composite

To develop methods for quantifying the effects of the microstructural variations of woven ceramic matrix composites on the effective properties and response of the material, a research program has been undertaken which is described in this paper. In order to characterize and quantify the variations in the microstructure of a five harness satin weave, CVI SiC/SiC, composite material, specimens were serially sectioned and polished to capture images that detailed the fiber tows, matrix, and porosity. Open source quantitative image analysis tools were then used to isolate the constituents and collect relevant statistics such as within ply tow spacing. This information was then used to build two dimensional finite element models that approximated the observed section geometry. With the aid of geometrical models generated by the microstructural characterization process, finite element models were generated and analyses were performed to quantify the effects of the microstructure and its variation on the effective stiffness and areas of stress concentration of the material. The results indicated that the geometry and distribution of the porosity appear to have significant effects on the through-thickness modulus. Similarly, stress concentrations on the outer surface of the composite appear to correlate to regions where the transverse tows are separated by a critical amount.

Goldberg, Robert K.↗