Search NASA⌕ Search

SEARCH · Search NASA

Results for “randomized methods”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

A new method for predicting response in complex linear systems. II

The paper describes convergence and response studies in the low frequency range of complex systems, particularly with low values of damping of different distributions, and reports on the modification of the relaxation procedure required under these conditions. A new method is presented for response estimation in complex lumped parameter linear systems under random or deterministic steady state excitation. The essence of the method is the use of relaxation procedures with a suitable error function to find the estimated response; natural frequencies and normal modes are not computed. For a 45 degree of freedom system, and two relaxation procedures, convergence studies and frequency response estimates were performed. The low frequency studies are considered in the framework of earlier studies (Kayser and Bogdanoff, 1975) involving the mid to high frequency range.

Bogdanoff, J. L.↗
DOE OSTI2017

Second-Order Uncertainty Analysis for Radiation Shielding Problems Using MCNP

Standard uncertainty propagation is based on a first-order Taylor expansion. This approximation is usually justified by assuming that 'uncertainties are very small numbers compared with the values of their corresponding quantities'. Furthermore, the second derivatives needed for a second-order analysis are not usually available, and a second-order analysis requires more knowledge regarding the distribution of the random variables. For radiation shielding problems, the MCNP6 Monte Carlo code computes first and second derivatives of a response with respect to nuclide density or reaction cross section. This paper explores the application of these derivatives to second-order uncertainty analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A fast Monte Carlo cell-by-cell simulation for radiobiological effects in targeted radionuclide therapy using pre-calculated single-particle track standard DNA damage data

Introduction: We developed a new method that drastically speeds up radiobiological Monte Carlo radiation-track-structure (MC-RTS) calculations on a cell-by-cell basis. Methods: The technique is based on random sampling and superposition of single-particle track (SPT) standard DNA damage (SDD) files from a “pre-calculated” data library, constructed using the RTS code TOPAS-nBio, with “time stamps” manually added to incorporate dose-rate effects. This time-stamped SDD file can then be input into MEDRAS, a mechanistic kinetic model that calculates various radiation-induced biological endpoints, such as DNA double-strand breaks (DSBs), misrepairs and chromosomal aberrations, and cell death. As a benchmark validation of the approach, we calculated the predicted energy-dependent DSB yield and the ratio of direct-to-total DNA damage, both of which agreed with published in vitro experimental data. We subsequently applied the method to perform a superfast cell-by-cell simulation of an experimental in vitro system consisting of neuroendocrine tumor cells uniformly incubated with 177 Lu. Results and discussion: The results for residual DSBs, both at 24 and 48 h post-irradiation, are in line with the published literature values. Our work serves as a proof-of-concept demonstration of the feasibility of a cost-effective “in silico clonogenic cell survival assay” for the computational design and development of radiopharmaceuticals and novel radiotherapy treatments more generally.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Fe Vacancy Order and Domain Distribution in A{sub x}Fe{sub 2−y}Se{sub 2}

The Fe vacancy and its ordering pattern is a unique property of the A{sub x}Fe{sub 2−y}Se{sub 2} (A = K, Rb, Cs) superconductors. Below the structural transition temperature T{sub S} ∼ 550 K, the system phase separates into the I4/mmm phase with no Fe vacancy and the I4/m phase with the √5×√5×1 vacancy order. In this work, an Ising-like Hamiltonian defined on the Fe-sublattice is constructed and employed in a Metropolis Monte Carlo simulation to study the Fe vacancy order and the domain distribution under different thermal treatments. A third phase with randomly distributed Fe vacancies and space group I4/mmm is observed in the simulation. The response of the volume fraction of this simulated phase under thermal treatment shows a direct connection to the change of the superconducting shielding fraction under annealing and quenching processes, indicating a close relation between the vacancy disordered I4/mmm phase and the superconducting phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Performance characterization of freeform finished surfaces of potassium dihydrogen phosphate using fluid jet polishing with a nonaqueous slurry

Abstract Potassium dihydrogen phosphate (KDP) and its deuterated analog (DKDP) are unique nonlinear optical materials for high power laser systems. They are used widely for frequency conversion and polarization control by virtue of the ability to grow optical-quality crystals at apertures suitable for fusion-class laser systems. Existing methods for freeform figuring of KDP/DKDP optics do not produce surfaces with sufficient laser-induced–damage thresholds (LIDT’s) for operation in the ultraviolet portion of high-peak-power laser systems. In this work, we investigate fluid jet polishing (FJP) using a nonaqueous slurry as a sub-aperture finishing method for producing freeform KDP surfaces. This method was used to selectively polish surface areas to different depths on the same substrate with removals ranging from 0.16 μm to 5.13 μm. The finished surfaces demonstrated a slight increase in roughness as the removal depth increased along with a small number of fracture pits. Laser damage testing with 351 nm, 1 ns pulses demonstrated excellent surface damage thresholds, with the highest values in areas devoid of fracture pits. This work demonstrates, for the first time, a method that enables fabrication of a waveplate that provides tailored polarization randomization that can be scaled to meter-sized optics. Furthermore, this method is based on FJP technology that incorporates a nonaqueous slurry specially designed for use with KDP. This novel nonaqueous FJP process can be also used for figuring other types of materials that exhibit similar challenging inherent properties such as softness, brittleness, water-solubility, and temperature sensitivity.

36 MATERIALS SCIENCE↗
NASA NTRS2020

Electromagnetic Scattering by Discrete Random Media Illuminated by a Gaussian Beam II: Solution of the Radiative Transfer Equation

In this paper, we present numerical methods for solving the phenomenological scalar radiative transfer equation for a discrete random medium illuminated by a Gaussian beam. These rely on the Fourier transform method for the horizontal variables and the discrete ordinate method with matrix exponential for solving the underlying one-dimensional radiative transfer equation in the wavenumber domain. The problem of a Gaussian beam at oblique and normal incidence, as well as, the searchlight problem are treated. A complete description of the methods and the numerical algorithms is provided.

Gaussian beam↗
NASA NTRS2016

Multiaxis Rainflow Fatigue Methods for Nonstationary Vibration

Mechanical structures and components may be subjected to cyclical loading conditions, including sine and random vibration. Such systems must be designed and tested accordingly. Rainflow cycle counting is the standard method for reducing a stress time history to a table of amplitude-cycle pairings prior to the Palmgren-Miner cumulative damage calculation. The damage calculation is straightforward for sinusoidal stress but very complicated for random stress, particularly for nonstationary vibration. This paper evaluates candidate methods and makes a recommendation for further study of a hybrid technique.

Irvine, T.↗
DOE OSTI2019

Charge transfer excitation energies from ground state density functional theory calculations

Calculating charge transfer (CT) excitation energies with high accuracy and low computational cost is a challenging task. Kohn-Sham density functional theory (KS-DFT), due to its efficiency and accuracy, has achieved great success in describing ground state problems. To extend to excited state problems, our group recently demonstrated an approach with good numerical results to calculate low-lying and Rydberg excitation energies of an N-electron system from a ground state KS or generalized KS calculations of an (N-1)-electron system via its orbital energies. In the present work, we explore further the same methodology to describe CT excitations. Numerical results from this work show that performance of conventional density functional approximations (DFAs) is not as good for CT excitations as for other excitations due to the delocalization error. Applying localized orbital scaling correction (LOSC) to conventional DFAs, a recently developed method in our group to effectively reduce the delocalization error, can improve the results. Overall, the performance of this methodology is better than time dependent DFT (TDDFT) with conventional DFAs. In addition, it shows that results from LOSC-DFAs in this method reach similar accuracy to other methods, such as ΔSCF, G 0 W 0 with Bethe-Salpeter equations, particle-particle random phase approximation, and even high-level wavefunction methods like CC2. Our analysis shows that the correct 1/R trend for CT excitation can be captured from LOSC-DFA calculations, stressing that the application of DFAs with the minimal delocalization error is essential within this methodology. This work provides an efficient way to calculate CT excitation energies from ground state DFT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Computational methods of robust controller design for aerodynamic flutter suppression

The development of Riccati iteration, a tool for the design and analysis of linear control systems is examined. First, Riccati iteration is applied to the problem of pole placement and order reduction in two-time scale control systems. Order reduction, yielding a good approximation to the original system, is demonstrated using a 16th order linear model of a turbofan engine. Next, a numerical method for solving the Riccati equation is presented and demonstrated for a set of eighth order random examples. A literature review of robust controller design methods follows which includes a number of methods for reducing the trajectory and performance index sensitivity in linear regulators. Lastly, robust controller design for large parameter variations is discussed.

Anderson, L. R.↗
NASA NTRS1972

On coherent oscillations of a string.

Vibrations of an elastic string when the separation between the ends varies randomly are studied. The emphasis is on the evolution of the coherent, or ordered, oscillations of the string. Using a perturbation technique borrowed from quantum field theory and the modified Kryloff-Bogoliuboff method, the 'multiple scattering' effect of the random separation between the ends on the linear and nonlinear coherent oscillations are investigated. It is found that due to the random interactions the coherent fundamental oscillation as well as the harmonies are damped. Their frequencies are also modified.

Liu, C. H.↗
DOE OSTI2019

Photon Diffusion in Random Media and Anisotropy of Scattering in the Henyey–Greenstein and Rayleigh–Gans Models

Monte Carlo (MC) numerical simulation is carried out for the intensity of multiply backscattered radiation as a function of the source–receiver distance for the Henyey–Greenstein (HG) and Rayleigh–Gans (RG) anisotropic scattering models. It is found that, in spite of the standard diffusive regime assumption on the universality of the description of scattering in terms of the extinction coefficient and the transport length, the multiple scattering intensity depends on the form of the phase function and on the degree of anisotropy even for systems with identical extinction coefficients and transport lengths. In this case, the results obtained by the RG phase function are more sensitive to the scattering ansisotropy than those obtained by the HG function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1974

Buckling of cylindrical shells with random imperfections

Buckling of long cylindrical shells with random imperfections subjected to axial loads is treated using the method of truncated hierarchy. A system of homogeneous variational equations is set up to examine the existence of bifurcation in the neighborhood of the equilibrium state. These equations are linear and involve stationary random coefficients. The truncated-hierarchy method is applied, and characteristic equations are obtained. Various exponential cosine correlation functions associated with asymmetric imperfections are examined numerically. Qualitatively, the results obtained are as anticipated.

Fersht, R. S.↗
NASA NTRS1992

Recognition and characterization of hierarchical interstellar structure. II - Structure tree statistics

A new method of image analysis is described, in which images partitioned into 'clouds' are represented by simplified skeleton images, called structure trees, that preserve the spatial relations of the component clouds while disregarding information concerning their sizes and shapes. The method can be used to discriminate between images of projected hierarchical (multiply nested) and random three-dimensional simulated collections of clouds constructed on the basis of observed interstellar properties, and even intermediate systems formed by combining random and hierarchical simulations. For a given structure type, the method can distinguish between different subclasses of models with different parameters and reliably estimate their hierarchical parameters: average number of children per parent, scale reduction factor per level of hierarchy, density contrast, and number of resolved levels. An application to a column density image of the Taurus complex constructed from IRAS data is given. Moderately strong evidence for a hierarchical structural component is found, and parameters of the hierarchy, as well as the average volume filling factor and mass efficiency of fragmentation per level of hierarchy, are estimated. The existence of nested structure contradicts models in which large molecular clouds are supposed to fragment, in a single stage, into roughly stellar-mass cores.

Houlahan, Padraig↗
DOE OSTI2025

DROP DURABILITY ASSESSMENT OF ELECTRONIC ASSEMBLIES UNDER OFF-AXIS LOADING WITH SKEWED FIXTURES

This thesis studies drop durability of electronic assemblies when the acceleration vector is oriented at 45° to the out-of-plane direction of the circuit card. The off-axis drop tests are accomplished with a skewed fixture and are conducted as a proxy for multiaxial drop testing. Advanced shock testing and vibration test methods have been developed over the last few decades to better represent real-world field environments during ground-based laboratory testing. However, many of these test methods require expensive and specialized equipment not available in most laboratories. An alternative approach for approximating simultaneous loading along multiple axes on conventional equipment utilizes skewed fixtures which have seen use in off-axis random vibration and drop impact testing. These methods generally rely on the conversion of a uniaxial input load from the test equipment (using a uniaxial drop tower or shaker) into a multiaxial load when resolved in the reference frame of the test article (mounted on a skewed fixture). Skewed fixture design is presented and recommendations for conducting skewed angle drop testing are introduced based on local measurements along the skewed face of the fixture to accurately monitor the impact event. Characterization tests were performed with a skewed fixture, at simultaneous acceleration loads from 500 to 3,000 g in two (in-plane and out-of-plane) directions, while meeting standard time domain tolerances. Upon experimental characterization, drop shock durability tests were conducted on a printed circuit assembly (PCA). Mean drops-to-failure were measured and quantified with Weibull statistics. Dominant solder joint failure modes were identified via failure analysis. Prior work on inclined angle impact testing is limited, and the majority of solder joint interconnect level fatigue studies are conducted considering perpendicular loading normal the circuit card. Low-cycle fatigue curves are generated based on plastic strain and plastic work density within the solder joint. A multiscale nonlinear finite element model is used to relate board-level flexure to solder joint interconnect level plastic strain. A high strain rate solder constitutive model allows for accurate modeling of solder plasticity resulting from high-impact drop shock. Fatigue parameters are computed from the Coffin-Manson relation and Palmgren-Miner damage accumulation. This work serves to apply established low-cycle fatigue methods for conventional drop shock loading (impact normal to circuit card) to non-perpendicular loading with a skewed fixture.

Hower, Jonathan [Kansas City National Security Cam↗
DOE OSTI2025

Increasing the hardness of posiform planting using random QUBOs for programmable quantum annealer benchmarking

Posiform planting is a method for constructing QUBO instances with a unique planted solution that can be tailored to arbitrary connectivity graphs. In this study we investigate making posiform planted QUBOs computationally harder by fusing many smaller random Ising models, whose global minimum is computed classically, with posiform planted QUBOs. The unique ground state of the resulting QUBO is the concatenation of (exactly one of) the ground states of each smaller problem. Our method generates QUBO instances that have a unique solution, are native to the hardware graph, and have tunable computational hardness. We use our QUBOs to benchmark three D-Wave quantum annealing processors (with 563–5627 qubits), and compare them against simulated annealing and Gurobi. Surprisingly, we find that the D-Wave ground state sampling success rate is not dependent on the glued random QUBO size, and that some QUBO classes are solved at high success rates at short annealing times on the Zephyr processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗