Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Truncated distribution”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Statistical theory of cubic Langmuir turbulence

The cubic direct interaction approximation is applied to a truncated (in Fourier space) version of the cubically nonlinear Schroedinger equation model of Langmuir physics. The results are compared (in the three-mode case) to those for an ensemble of numerical solutions of the dynamical equations with 10,000 different sets of Gaussianly distributed initial conditions. In the undriven, undamped case, the statistical theory (but not the ensemble) evolves to a state of thermal equilibrium. In the driven, damped case, the statistical theory appears to evolve to a state close to that corresponding to one of the limit cycles of the dynamical equations.

Sun, G.-Z.↗
NASA NTRS1988

Estimating residual fault hitting rates by recapture sampling

For the recapture debugging design introduced by Nayak (1988) the problem of estimating the hitting rates of the faults remaining in the system is considered. In the context of a conditional likelihood, moment estimators are derived and are shown to be asymptotically normal and fully efficient. Fixed sample properties of the moment estimators are compared, through simulation, with those of the conditional maximum likelihood estimators. Properties of the conditional model are investigated such as the asymptotic distribution of linear functions of the fault hitting frequencies and a representation of the full data vector in terms of a sequence of independent random vectors. It is assumed that the residual hitting rates follow a log linear rate model and that the testing process is truncated when the gaps between the detection of new errors exceed a fixed amount of time.

Lee, Larry↗
DOE OSTI2022

Computational Design of Miniproteins as SARS-CoV-2 Therapeutic Inhibitors

A rational therapeutic strategy is urgently needed for combating SARS-CoV-2 infection. Viral infection initiates when the SARS-CoV-2 receptor-binding domain (RBD) binds to the ACE2 receptor, and thus, inhibiting RBD is a promising therapeutic for blocking viral entry. In this study, the structure of lead antiviral candidate binder (LCB1), which has three alpha-helices (H1, H2, and H3), is used as a template to design and simulate several miniprotein RBD inhibitors. LCB1 undergoes two modifications: structural modification by truncation of the H3 to reduce its size, followed by single and double amino acid substitutions to enhance its binding with RBD. We use molecular dynamics (MD) simulations supported by ab initio density functional theory (DFT) calculations. Complete binding profiles of all miniproteins with RBD have been determined. The MD investigations reveal that the H3 truncation results in a small inhibitor with a -1.5 kcal/mol tighter binding to RBD than original LCB1, while the best miniprotein with higher binding affinity involves D17R or E11V + D17R mutation. DFT calculations provide atomic-scale details on the role of hydrogen bonding and partial charge distribution in stabilizing the minibinder:RBD complex. This study provides insights into general principles for designing potential therapeutics for SARS-CoV-2.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1972

The stability of motion of satellites with flexible appendages

The mathematical formulation associated with the problem of stability of motion of a satellite consisting of a main rigid body and three (or less) pairs of flexible rods is presented. The rods are capable of flexure in two orthogonal directions. Whereas the rotational motion of the body is described by generalized coordinates depending on time alone, the elastic displacements of the rods depend both on spatial position and time. Assuming no external torques, there exist motion integrals in the form of momentum integrals. These integrals can be regarded as constraint equations relating the system velocities, and used to reduce the number of variables describing the motion. The stability analysis has been carried out by means of an extension of the Liapunov direct method. Since the elastic vibrations result in energy dissipation, it is shown that the equilibrium position is asymptotically stable if the Hamiltonian is positive definite and unstable if it can take negative values in the neighborhood of the equilibrium. Determining the sign definiteness of the Hamiltonian is complicated by the fact that it contains spatial derivatives of the elastic displacements. Two methods are presented to cope with this problem. The first, the standard modal analysis in conjunction with series truncation, develops criteria in terms of infinite series associated with the natural modes and frequencies of the elastic rods. The second, the method of integral coordinates, yields closed-form stability criteria involving the system parameters, such as the body moments of inertia, the length and mass distribution of the rods, the lowest natural frequencies of the rods, and the satellite spin velocity.

Meirovitch, L.↗
NASA NTRS1994

Compression strength of composite primary structural components

The linear elastic response is determined for an internally pressurized, long circular cylindrical shell stiffened on the inside by a regular arrangement of identical stringers and identical rings. Periodicity of this configuration permits the analysis of a portion of the shell wall centered over a generic stringer-ring joint; i.e., a unit cell model. The stiffeners are modeled as discrete beams, and the stringer is assumed to have a symmetrical cross section and the ring an asymmetrical section. Asymmetery causes out-of-plane bending and torsion of the ring. Displacements are assumed as truncated double Fourier series plus simple terms in the axial coordinate to account for the closed and pressure vessel effect (a non-periodic effect). The interacting line loads between the stiffeners and the inside shell wall are Lagrange multipliers in the formulation, and they are also assumed as truncated Fourier series. Displacement continuity constraints between the stiffeners and shell along the contact lines are satisfied point-wise. Equilibrium is imposed by the principle of virtual work. A composite material crown panel from the fuselage of a large transport aircraft is the numerical example. The distributions of the interacting line loads, and the out-of-plane bending moment and torque in the ring, are strongly dependent on modeling the deformations due to transverse shear and cross-sectional warping of the ring in torsion. This paper contains the results from the semiannual report on research on 'Pressure Pillowing of an Orthogonally Stiffened Cylindrical Shell'. The results of the new work are illustrated in the included appendix.

Johnson, Eric R.↗
NASA NTRS1995

Microlensing by the galactic bar

We compute the optical depth and duration distribution of microlensing events towrd Baade's window in a model composed of a Galactic disk and a bar. The bar model is a self-consistent dynamical model built out of individual orbits that has been populated to be consistent with the COBE maps of the Galaxy and kinematic observations of the Galactic bulge. We find that most of the lenses are in the bulge with a line-of-sight distance 6.25 kpc (adopting R(sub 0) = 8 kpc). The microlensing optical depth of a 2 x 10(exp 10) solar mass bar plus a truncated disk is (2.2 +/- 0.45) x 10(exp -6), consistent with the large optical depth (3.2 +/- 1.2) x 10(exp -6) found by Udalski et al. (1994). This model optical depth is enhanced over the predictions of axisymmetric models by Kiraga & Paczynski (1994) by slightly more than a factor of 2, since the bar is elongated along the line of sight. The large Einstein radius and small transverse velocity dispersion also predict a longer event duration in the self-consistent bar model than in the Kiraga-Paczynski model. The event rate and duration distribution also depend on the lower mass cutoff of the lens mass function. With a 0.1 solar mass cutoff, five to seven events (depending on the contribution of disk lenses) with a logarithmic mean duration of 20 days are expected for the Optical Gravitational Lensing Experiment (OGLE) according to our model, while Udalski et al. (1994) observed nine events with durations from 8 to 62 days. On the other hand, if most of the lenses are brown dwarfs, our model predicts too many short-duration events. A Kolmogorov-Smirnov test finds only 7% probability for the model with 0.01 solar mass cutoff to be consistent with current data.

Zhao, Hongsheng↗
DOE OSTI2023

Quantum efficiency enhancement in simulated nanostructured negative electron affinity GaAs photocathodes

We report nanostructured negative electron affinity GaAs photocathodes for a polarized electron source are studied using finite difference time domain optical simulation. The structures studied are nanosquare columns, truncated nanocones, and truncated nanopyramids. Mie-type resonances in the 700–800 nm waveband, suitable for generation of polarized electrons, are identified. At resonance wavelengths, the nanostructures can absorb up to 99% of the incident light. For nanosquare columns and truncated nanocones, the maximum quantum efficiency (QE) at 780 nm obtained from simulation is 27%, whereas for simulated nanopyramids, the QE is ~21%. The high photocathode quantum efficiency is due to the shift of Mie resonance toward the longer wavelength, leading to increased light absorption. The field profile distribution shows the excitation of dipole and quadrupole modes within the nanostructures at resonant frequencies. This leads to enhanced photoabsorption and photoelectron generation closer to emission surfaces than for a flat photocathode. The enhanced photoabsorption and reduced electron transport distance for the nanostructured photocathode enhance its QE compared to that for the flat surface wafer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Trivalent ion overcharging on electrified graphene

We report the structure of the electrical double layer (EDL) formed near graphene in aqueous environments strongly impacts its performance for a plethora of applications, including capacitive deionization. In particular, adsorption and organization of multivalent counterions near the graphene interface can promote nonclassical behaviors of EDL including overcharging followed by co-ion adsorption. In this paper, we characterize the EDL formed near an electrified graphene interface in dilute aqueous YCl 3 solution using in situ high resolution x-ray reflectivity (also known as crystal truncation rod) and resonant anomalous x-ray reflectivity (RAXR). These interface-specific techniques reveal the electron density profiles with molecular-scale resolution. We find that yttrium ions (Y 3+ ) readily adsorb to the negatively charged graphene surface to form an extended ion profile. This ion distribution resembles a classical diffuse layer but with a significantly high ion coverage, i.e., 1 Y 3+ per 11.4 ± 1.6 Å 2 , compared to the value calculated from the capacitance measured by cyclic voltammetry (1 Y 3+ per ~240 Å 2 ). Such overcharging can be explained by co-adsorption of chloride that effectively screens the excess positive charge. The adsorbed Y 3+ profile also shows a molecular-scale gap (≥ 5 Å) from the top graphene surfaces, which is attributed to the presence of intervening water molecules between the adsorbents and adsorbates as well as the lack of inner-sphere surface complexation on chemically inert graphene. We also demonstrate controlled adsorption by varying the applied potential and reveal consistent Y 3+ ion position with respect to the surface and increasing cation coverage with increasing the magnitude of the negative potential. This is the first experimental description of a model graphene-aqueous system with controlled potential and provides important insights into the application of graphene-based systems for enhanced and selective ion separations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

Constrained and unconstrained modes - Some modeling aspects of flexible spacecraft

Spacecraft that are partially rigid and partially flexible may be dynamically modeled in terms of either 'constrained' modes of vibration, for which the rigid part is held motionless, or the 'unconstrained' modes of the whole vehicle. A model 'completeness index' is defined for each expansion and used as a truncation criterion. Using as an example a large flexible platform with a rigid body at the center, it is shown that as the rigid portion of the vehicle becomes very small, a great many constrained modes are needed to achieve satisfactory model completeness. Otherwise, one finds an erroneous distribution of angular momentum among the modes. It is also shown that the importance of unconstrained modes is not necessarily ordered by frequency.

Hablani, H. B.↗
DOE OSTI2024

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗
NASA NTRS1989

Maximum likelihood estimation for distributed parameter models of flexible spacecraft

A distributed-parameter model of the NASA Solar Array Flight Experiment spacecraft structure is constructed on the basis of measurement data and analyzed to generate a priori estimates of modal frequencies and mode shapes. A Newton-Raphson maximum-likelihood algorithm is applied to determine the unknown parameters, using a truncated model for the estimation and the full model for the computation of the higher modes. Numerical results are presented in a series of graphs and briefly discussed, and the significant improvement in computation speed obtained by parallel implementation of the method on a supercomputer is noted.

Taylor, L. W., Jr.↗
NASA NTRS2004

Over-Luminous Elliptical Galaxies

The first paper from our work has been completed and accepted for publication. Another paper presents a study of the ESO 30601 70 galaxy group, combining Chandra, XMM-Newton, and optical observations. We find that the system is a true fossil galaxy group - a group whose optical light is dominated by a single galaxy. The group X-ray emission is composed of a central, dense, cool core (10 kpc in radius) and an isothermal medium beyond the central 10 kpc. The region between 10 and 50 kpc (the cooling radius) has the same temperature as the gas from 50 to 400 kpc, although the gas cooling time between 10 and 50 kpc (2-6 Gyr) is shorter than the Hubble time. Thus, the ESO 3060170 group does not have a group-sized cooling core. We suggest that the group cooling core may have been heated by a central active galactic nucleus (AGN) outburst in the past and that the small, dense, cool core is the truncated relic of a previous cooling core. The Chandra observations also reveal a variety of X-ray features in the central region, including a finger, an edge-like feature, and a small tail, all aligned along a north-south axis, as are the galaxy light and group galaxy distribution. The proposed AGN outburst may cause gas to slosh around the center and produce these asymmetric features. The observed flat temperature profile to 1/3rvir is not consistent with the predicted temperature profile in recent numerical simulations. We compare the entropy profile of the ESO 3060170 group with those of three other groups and find a flatter relation than that predicted by simulations involving only shock heating, S approximately r approximately 0.85. This is direct evidence of the importance of non-gravitational processes in group centers. We derive the mass profiles within 1/3rvir and find that the ESO 3060170 group is the most massive fossil group known.

Forman, William↗
NASA NTRS1993

Applications of computer algebra to distributed parameter systems

In the analysis of vibrations of continuous elastic systems, one often encounters complicated transcendental equations with roots directly related to the system's natural frequencies. Typically, these equations contain system parameters whose values must be specified before a numerical solution can be obtained. The present paper presents a method whereby the fundamental frequency can be obtained in analytical form to any desired degree of accuracy. The method is based upon truncation of rapidly converging series involving inverse powers of the system natural frequencies. A straightforward method to developing these series and summing them in closed form is presented. It is demonstrated how Computer Algebra can be exploited to perform the intricate analytical procedures which otherwise would render the technique difficult to apply in practice. We illustrate the method by developing two analytical approximations to the fundamental frequency of a vibrating cantilever carrying a rigid tip body. The results are compared to the numerical solution of the exact (transcendental) frequency equation over a range of system parameters.

Storch, Joel A.↗
DOE OSTI2024

Solving differential‐algebraic equations in power system dynamic analysis with quantum computing

Abstract Power system dynamics are generally modeled by high dimensional non‐linear differential‐algebraic equations (DAEs) given a large number of components forming the network. These DAEs' complexity can grow exponentially due to the increasing penetration of distributed energy resources, whereas their computation time becomes sensitive due to the increasing interconnection of the power grid with other energy systems. This paper demonstrates the use of quantum computing algorithms to solve DAEs for power system dynamic analysis. We leverage a symbolic programming framework to equivalently convert the power system's DAEs into ordinary differential equations (ODEs) using index reduction methods and then encode their data into qubits using amplitude encoding. The system non‐linearity is captured by Hamiltonian simulation with truncated Taylor expansion so that state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can solve the power system's DAEs accurately with a computational complexity polynomial in the logarithm of the system dimension. We also illustrate the use of recent advanced tools in scientific machine learning for implementing complex computing concepts, that is, Taylor expansion, DAEs/ODEs transformation, and quantum computing solver with abstract representation for power engineering applications.

computational complexity↗
DOE OSTI2020

Ultrasonic waves from radial mode excitation of a piezoelectric disc on the surface of an elastic solid

In this work, we discuss the ultrasonic wave propagation characteristics in an isotropic elastic solid due to the radial mode excitation of a piezoelectric disc actuator bonded to its surface. Finite element simulations using coupled electromechanical modeling is employed to investigate the wave propagation behavior. We find that the radial mode vibrations on the surface of the elastic solid generates all the three types of ultrasonic waves: longitudinal, shear, and surface waves. The waves in the solid are comprised of a central lobe and multiple side lobes based on the frequency of the radial mode excitation. The central lobe is predominantly composed of longitudinal waves and the side lobes are composed of shear waves. While the longitudinal waves have a strong central lobe and weak side lobes that are not fully developed, shear waves consists of fully developed side lobes. In addition, we observe that the longitudinal waves have fewer side lobes inside the solid compared to the shear waves. A semi-analytical approach is presented to explain the above observations. The interface traction boundary condition between the piezoelectric disc and the elastic solid is approximated as a truncated 'Bessel' excitation over the area of the piezoelectric disc in contact with the solid. A Fourier–Bessel series expansion technique is used to investigate the wave propagation behavior from the above traction boundary condition. The results obtained explain the observations from the finite element simulations with regard to the number and distribution of side lobes pertaining to the longitudinal and shear waves generated from the radial modes of the piezoelectric disc. The proposed semi-analytical approach is general and can be applied to any arbitrary axisymmetric excitations on the surface of an elastic solid.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Large space system control technology status and accomplishments

In the area of antenna modeling and control, preliminary structural models were defined for two representative parabolic reflectors. A control system design evolved for attitude control of the reflectors. The controller design was based on a lumped control concept where the control hardware (sensor and actuators) was mounted at the base of the antenna. In the area of distributed control, static shape control techniques were worked out to establish a prescribed vehicle shape. The corresponding estimation process for determination of vehicle shape from selected sensor measurements was also developed. A model order reduction study was conducted to find the best preflight dynamical models for on board controller design. The problems caused by truncation of the vehicle dynamics required to minimize on board computations were investigated. Estimator designs were developed for on board detection of large structure model errors.

Rodriguez, G.↗
NASA NTRS2007

Propagation of Computational Uncertainty Using the Modern Design of Experiments

This paper describes the use of formally designed experiments to aid in the error analysis of a computational experiment. A method is described by which the underlying code is approximated with relatively low-order polynomial graduating functions represented by truncated Taylor series approximations to the true underlying response function. A resource-minimal approach is outlined by which such graduating functions can be estimated from a minimum number of case runs of the underlying computational code. Certain practical considerations are discussed, including ways and means of coping with high-order response functions. The distributional properties of prediction residuals are presented and discussed. A practical method is presented for quantifying that component of the prediction uncertainty of a computational code that can be attributed to imperfect knowledge of independent variable levels. This method is illustrated with a recent assessment of uncertainty in computational estimates of Space Shuttle thermal and structural reentry loads attributable to ice and foam debris impact on ascent.

DeLoach, Richard↗
DOE OSTI2020

A Data-Driven Nonparametric Approach for Probabilistic Load-Margin Assessment Considering Wind Power Penetration

A modern power system is characterized by an increasing penetration of wind power, which results in large uncertainties in its states. These uncertainties must be quantified properly; otherwise, the system security may be threatened. Facing this challenge, here we propose a cost-effective, data-driven approach to assessing a power system's load margin probabilistically. Using actual wind data, a kernel density estimator is applied to infer the nonparametric wind speed distributions, which are further merged into the framework of a vine copula. The latter enables us to simulate complex multivariate and highly dependent model inputs with a variety of bivariate copulae that precisely represent the tail dependence in the correlated samples. Furthermore, to reduce the prohibitive computational time of traditional Monte-Carlo simulations that process a large amount of samples, we propose to use a nonparametric, Gaussian-process-emulator-based reduced-order model to replace the original complicated continuation power-flow model through a Bayesian-learning framework. To accelerate the convergence rate of this Bayesian algorithm, a truncated polynomial chaos surrogate, which serves as a highly efficient, parametric Bayesian prior, is developed. This emulator allows us to execute the time-consuming continuation power-flow solver at the sampled values with a negligible computational cost. Results of simulations that are performed on several test systems reveal the impressive performance of the proposed method in the probabilistic load-margin assessment.

17 WIND ENERGY↗