Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equivalence 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 541 records · Page 30

A Novel Approach for Adaptive Signal Processing

Adaptive linear predictors have been used extensively in practice in a wide variety of forms. In the main, their theoretical development is based upon the assumption of stationarity of the signals involved, particularly with respect to the second order statistics. On this basis, the well-known normal equations can be formulated. If high- order statistical stationarity is assumed, then the equivalent normal equations involve high-order signal moments. In either case, the cross moments (second or higher) are needed. This renders the adaptive prediction procedure non-blind. A novel procedure for blind adaptive prediction has been proposed and considerable implementation has been made in our contributions in the past year. The approach is based upon a suitable interpretation of blind equalization methods that satisfy the constant modulus property and offers significant deviations from the standard prediction methods. These blind adaptive algorithms are derived by formulating Lagrange equivalents from mechanisms of constrained optimization. In this report, other new update algorithms are derived from the fundamental concepts of advanced system identification to carry out the proposed blind adaptive prediction. The results of the work can be extended to a number of control-related problems, such as disturbance identification. The basic principles are outlined in this report and differences from other existing methods are discussed. The applications implemented are speech processing, such as coding and synthesis. Simulations are included to verify the novel modelling method.

Chen, Ya-Chin↗

Fine structure of the isovector giant dipole resonance in Nd 142 – 150 and Sm 152

Inelastic proton scattering at energies of a few hundred MeV and very-forward scattering angles including 0° has been established as a tool for the study of electric-dipole strength distributions in nuclei. The present work reports a systematic investigation of the chain of stable even-mass Nd isotopes representing a transition from spherical to quadrupole-deformed nuclei. Extraction of the equivalent photo-absorption cross sections and analysis of their fine structure in the energy region of the isovector giant dipole resonance (IVGDR). Method: Proton inelastic scattering reactions of 200 MeV protons were measured at the iThemba Laboratory for Accelerator Based Sciences in Cape Town, South Africa. The scattering products were momentum-analyzed by the K600 magnetic spectrometer positioned at θ Lab = 0°. Using dispersion-matching techniques, energy resolutions of ΔE ≈ 40–50 keV (full width at half maximum) were obtained. After subtraction of background and contributions from other multipoles, the spectra were converted to photoabsorption cross sections using the equivalent virtual-photon method. Wavelet-analysis techniques are used to extract characteristic energy scales of the fine structure of the IVGDR from the experimental data. Fine structure of the IVGDR is observed even for the most deformed nuclei studied. Comparisons between the extracted experimental energy scales and those energy scales obtained from the quasiparticle-phonon model (QPM) and Skyrme separable random phase approximation (SSRPA) predictions provide insight into the role of different giant-resonance damping mechanisms. It can be seen that the scales in the spherical and most likely also in the deformed nuclei mainly result from the fragmentation of the one-particle-one-hole (1p1h) strength into several dominant transitions serving as doorway states. In cases where calculations beyond the 1p1h level are available, some impact of the spreading due to coupling of the two-particle-two-hole (2p2h) states to the 1p1h doorway states is observed. New virtual-photon absorption data for the chain of stable Nd isotopes and 152 Sm are presented, with a focus on the phenomenon of nonstatistical cross-section fluctuations, referred to as fine structure, in the energy region of the IVGDR. The wavelet-analysis techniques used allowed for the features of the fine structure to be quantified in the form of characteristic scales. Finally, comparisons between experimental results and model predictions indicate that Landau damping seems to be the main source of the fine structure in both the spherical and deformed nuclei, but calculations including 2p2h degrees of freedom would be beneficial to confirm this for the deformed cases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A comparison of several tapered wings designed to avoid tip stalling

Optimum proportions of tapered wings were investigated by a method that involved a comparison of wings designed to be aerodynamically equal. The conditions of aerodynamic equality were equality in stalling speed, in induced drag at a low speed, and in the total drag at cruising speed. After the wings were adjusted to aerodynamic equivalence, the weights of the wings were calculated as a convenient method of indicating the optimum wing. The aerodynamic characteristics were calculated from wing theory and test data for the airfoil sections. Various combinations of washout, camber increase in the airfoil sections from the center to the tips, and sharp leading edges at the center were used to bring about the desired equivalence of maximum lift and center-stalling characteristics. In the calculation of the weights of the wings, a simple type of spar structure was assumed that permitted an integration across the span to determine the web and the flange weights. The covering and the remaining weight were taken in proportion to the wing area. The total weights showed the wings with camber and washout to have the lowest weights and indicated the minimum for wings with a taper ratio between 1/2 and 1/3.

Anderson, Raymond F↗

A general approach to seismic inversion with automatic differentiation

Imaging Earth structure or seismic sources from seismic data involves minimizing a target misfit function, and is commonly solved through gradient-based optimization. The adjoint-state method has been developed to compute the gradient efficiently; however, its implementation can be time-consuming and difficult. We develop a general seismic inversion framework to calculate gradients using reverse-mode automatic differentiation. The central idea is that adjoint-state methods and reverse-mode automatic differentiation are mathematically equivalent. Here, the mapping between numerical PDE simulation and deep learning allows us to build a seismic inverse modeling library, ADSeismic, based on deep learning frameworks, which supports high performance reverse-mode automatic differentiation on CPUs and GPUs. We demonstrate the performance of ADSeismic on inverse problems related to velocity model estimation, rupture imaging, earthquake location, and source time function retrieval. ADSeismic has the potential to solve a wide variety of inverse modeling applications within a unified framework.

58 GEOSCIENCES↗

The equivalent angle-of-attack concept for engineering analysis

An account is given of methods for the estimation of a wing-body-tail missile configuration's aerodynamic performance by means of the 'component buildup' method, in which the overall aerodynamic loads for an airframe are built up from the assumed characteristics of isolated airframe components and then subjected to additional loads generated by component interference. Attention is given to the cases of missile airframes in steady flow at arbitrary angles of attack and bank; the unifying thread of the present treatment is slender body theory, together with its nonlinear extensions through the equivalent angle-of-attack concept. The estimation of the forces and moments acting on each of the fins is emphasized, so that control cross-coupling effects as well as longitudinal and lateral directional effects can be determined.

Hemsch, M. J.↗

On Formulations of Discontinuous Galerkin and Related Methods for Conservation Laws

A formulation for the discontinuous Galerkin (DG) method that leads to solutions using the differential form of the equation (as opposed to the standard integral form) is presented. The formulation includes (a) a derivative calculation that involves only data within each cell with no data interaction among cells, and (b) for each cell, corrections to this derivative that deal with the jumps in fluxes at the cell boundaries and allow data across cells to interact. The derivative with no interaction is obtained by a projection, but for nodal-type methods, evaluating this derivative by interpolation at the nodal points is more economical. The corrections are derived using the approximate (Dirac) delta functions. The formulation results in a family of schemes: different approximate delta functions give rise to different methods. It is shown that the current formulation is essentially equivalent to the flux reconstruction (FR) formulation. Due to the use of approximate delta functions, an energy stability proof simpler than that of Vincent, Castonguay, and Jameson (2011) for a family of schemes is derived. Accuracy and stability of resulting schemes are discussed via Fourier analyses. Similar to FR, the current formulation provides a unifying framework for high-order methods by recovering the DG, spectral difference (SD), and spectral volume (SV) schemes. It also yields stable, accurate, and economical methods.

Huynh, H. T.↗

Structural dynamics and vibrations of damped, aircraft-type structures

Engineering preliminary design methods for approximating and predicting the effects of viscous or equivalent viscous-type damping treatments on the free and forced vibration of lightly damped aircraft-type structures are developed. Similar developments are presented for dynamic hysteresis viscoelastic-type damping treatments. It is shown by both engineering analysis and numerical illustrations that the intermodal coupling of the undamped modes arising from the introduction of damping may be neglected in applying these preliminary design methods, except when dissimilar modes of these lightly damped, complex aircraft-type structures have identical or nearly identical natural frequencies. In such cases, it is shown that a relatively simple, additional interaction calculation between pairs of modes exhibiting this 'modal response' phenomenon suffices in the prediction of interacting modal damping fractions. The accuracy of the methods is shown to be very good to excellent, depending on the normal natural frequency separation of the system modes, thereby permitting a relatively simple preliminary design approach. This approach is shown to be a natural precursor to elaborate finite element, digital computer design computations in evaluating the type, quantity, and location of damping treatment.

Young, Maurice I.↗

Equivalent plate modeling for conceptual design of aircraft wing structures

This paper describes an analysis method that generates conceptual-level design data for aircraft wing structures. A key requirement is that this data must be produced in a timely manner so that is can be used effectively by multidisciplinary synthesis codes for performing systems studies. Such a capability is being developed by enhancing an equivalent plate structural analysis computer code to provide a more comprehensive, robust and user-friendly analysis tool. The paper focuses on recent enhancements to the Equivalent Laminated Plate Solution (ELAPS) analysis code that significantly expands the modeling capability and improves the accuracy of results. Modeling additions include use of out-of-plane plate segments for representing winglets and advanced wing concepts such as C-wings along with a new capability for modeling the internal rib and spar structure. The accuracy of calculated results is improved by including transverse shear effects in the formulation and by using multiple sets of assumed displacement functions in the analysis. Typical results are presented to demonstrate these new features. Example configurations include a C-wing transport aircraft, a representative fighter wing and a blended-wing-body transport. These applications are intended to demonstrate and quantify the benefits of using equivalent plate modeling of wing structures during conceptual design.

Giles, Gary L.↗

On High-Order Upwind Methods for Advection

In the fourth installment of the celebrated series of five papers entitled "Towards the ultimate conservative difference scheme", Van Leer (1977) introduced five schemes for advection, the first three are piecewise linear, and the last two, piecewise parabolic. Among the five, scheme I, which is the least accurate, extends with relative ease to systems of equations in multiple dimensions. As a result, it became the most popular and is widely known as the MUSCL scheme (monotone upstream-centered schemes for conservation laws). Schemes III and V have the same accuracy, are the most accurate, and are closely related to current high-order methods. Scheme III uses a piecewise linear approximation that is discontinuous across cells, and can be considered as a precursor of the discontinuous Galerkin methods. Scheme V employs a piecewise quadratic approximation that is, as opposed to the case of scheme III, continuous across cells. This method is the basis for the on-going "active flux scheme" developed by Roe and collaborators. Here, schemes III and V are shown to be equivalent in the sense that they yield identical (reconstructed) solutions, provided the initial condition for scheme III is defined from that of scheme V in a manner dependent on the CFL number. This equivalence is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The finding also shows a key connection between the approaches of discontinuous and continuous polynomial approximations. In addition to the discussed equivalence, a framework using both projection and interpolation that extends schemes III and V into a single family of high-order schemes is introduced. For these high-order extensions, it is demonstrated via Fourier analysis that schemes with the same number of degrees of freedom 𝐾 per cell, in spite of the different piecewise polynomial degrees, share the same sets of eigenvalues and thus, have the same stability and accuracy. Moreover, these schemes are accurate to order 2𝐾−1, which is higher than the expected order of 𝐾.

high-order methods↗

Self-Consistent Simulation of the Brownian Stage of Dust Growth

It is a widely accepted view that in proto-planetary accretion disks the collision and following sticking of dust particles embedded in the gas eventually leads to the formation of planetesimals (coagulation). For the smallest dust grains, Brownian motion is assumed to be the dominant source of their relative velocities leading to collisions between these dust grains. As the dust grains grow they eventually couple to the turbulent motion of the gas which then drives the coagulation much more efficiently. Many numerical coagulation simulations have been carried out to calculate the fractal dimension of the aggregates, which determines the duration of the ineffective Brownian stage of growth. Predominantly on-lattice and off-lattice methods were used. However, both methods require simplification of the astrophysical conditions. The aggregates found by those methods had a fractal dimension of approximately 2 which is equivalent to a constant, mass-independent friction time. If this value were valid for the conditions in an accretion disk, this would mean that the coagulation process would finally 'freeze out' and the growth of a planetesimal would be impossible within the lifetime of an accretion disk. In order to investigate whether this fractal dimension is model independent, we simulate self-consistently the Brownian stage of the coagulation by an N-particle code. This method has the advantage that no further assumptions about homogeneity of the dust have to be made. In our model, the dust grains are considered as aggregates built up of spheres. The equation of motion of the dust grains is based on the probability density for the diffusive transport within the gas atmosphere. Because of the very low number density of the dust grains, only 2-body-collisions have to be considered. As the Brownian stage of growth is very inefficient, the system is to be simulated over long periods of time. In order to find close particle pairs of the system which are most likely to undergo a collision, we use a particle-in-cell (PIC) method for the early stages of the simulation where the system is still very homogeneous and a tree method later when the particles are more clustered.

Kempf, S.↗

Distributive, Non-destructive Real-time System and Method for Snowpack Monitoring

A ground-based system that provides quasi real-time measurement and collection of snow-water equivalent (SWE) data in remote settings is provided. The disclosed invention is significantly less expensive and easier to deploy than current methods and less susceptible to terrain and snow bridging effects. Embodiments of the invention include remote data recovery solutions. Compared to current infrastructure using existing SWE technology, the disclosed invention allows more SWE sites to be installed for similar cost and effort, in a greater variety of terrain; thus, enabling data collection at improved spatial resolutions. The invention integrates a novel computational architecture with new sensor technologies. The invention's computational architecture is based on wireless sensor networks, comprised of programmable, low-cost, low-powered nodes capable of sophisticated sensor control and remote data communication. The invention also includes measuring attenuation of electromagnetic radiation, an approach that is immune to snow bridging and significantly reduces sensor footprints.

Frolik, Jeff↗

Generalization of Equivalent Crystal Theory to Include Angular Dependence

In the original Equivalent Crystal Theory, each atomic site in the real crystal is assigned an equivalent lattice constant, in general different from the ground state one. This parameter corresponds to a local compression or expansion of the lattice. The basic method considers these volumetric transformations and, in addition, introduces the possibility that the reference lattice is anisotropically distorted. These distortions however, were introduced ad-hoc. In this work, we generalize the original Equivalent Crystal Theory by systematically introducing site-dependent directional distortions of the lattice, whose corresponding distortions account for the dependence of the energy on anisotropic local density variations. This is done in the spirit of the original framework, but including a gradient term in the density. This approach is introduced to correct a deficiency in the original Equivalent Crystal Theory and other semiempirical methods in quantitatively obtaining the correct ratios of the surface energies of low index planes of cubic metals (100), (110), and (111). We develop here the basic framework, and apply it to the calculation of Fe (110) and Fe (111) surface energy formation. The results, compared with first principles calculations, show an improvement over previous semiempirical approaches.

Ferrante, John↗

Atom Equivalent Energies for the Rapid Estimation of the Heat of Formation of Explosive Molecules from Density Functional Tight Binding Theory

Atom equivalent energies have been derived from which the gas-phase heat of formation of explosive molecules can be estimated from fast, semiempirical density functional tight binding total energy calculations. The root-mean-square deviation and maximum deviation of the heats of formation from the experimental values for the set of 45 energetic molecules compiled by Byrd and Rice [ J. Phys. Chem. A, 2006, 110, 1005–1013] are 10.4 and 25.5 kcal/mol, respectively, using 4 atom equivalent energies and 7.4 and 15.0 kcal/mol, respectively, using 7 atom equivalent energies. These errors are around a factor of 2–3 larger than those obtained from density functional theory calculations but are smaller than those obtained from other semiempirical electronic structure methods. Heats of formation calculated with density functional tight binding theory using the 4 and 7 atom equivalent energies, the Byrd and Rice scheme, and the atom pair contribution method for a new set of 531 energetic molecules that contain only carbon, hydrogen, nitrogen, and oxygen are provided.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A 24-channel ultra-low-noise preamplifier for dN/dx measurements with drift tube detectors

Cluster counting $(dN/dx)$ is a promising method to enhance particle identification for gaseous detectors, especially in next-generation collider experiments like the FCC-ee, where good $π/K$ separation over a broad momentum range is essential. However, its implementation in large-scale systems has been limited by the challenging requirements for high-resolution signal amplification and readout. This paper presents a 24-channel ultra-low-noise preamplifier board designed for drift tube detectors to enable $dN/dx$ measurements. The three-stage amplification topology employs SiGe transistors and integrates dedicated noise-minimization techniques, achieving a charge gain of 21.11 mV/fC from 0.3 fC to 50 fC, a bandwidth of 542 MHz, and a voltage gain of 47.8 dB. The measured voltage noise density is 0.35 nV/$\sqrt{\textrm{Hz}}$ , surpassing most of the state-of-the-art preamplifiers for gaseous and silicon detectors. Validation tests conducted on the sMDT chambers at the CERN Proton Synchrotron test beam facility demonstrate that the proposed design meets the stringent preamplifier requirements for implementing the $dN/dx$ method in drift-tube detector systems, achieving an equivalent noise charge of 0.14 fC and a signal-to-noise ratio of 73 when operated with a He:iC 4 H 10 (90:10) gas mixture. The design also shows promise for broader application in other gaseous or semiconductor detectors.

Drift tube detector↗

Ferroelectric HfO 2 and the importance of strain

Ferroelectric oxides based on HfO 2 show tremendous promise for the next generation of memory and logic devices. The ferroelectric polymorph is one of several that can be derived from the high symmetry cubic fluorite structure of HfO 2 . A single grain of HfO 2 may consist of a coherent mixture of multiple orientational and translational variants of different polymorphs. Here, we use symmetry-adapted strain-order parameters to elucidate the relationship between the different HfO 2 polymorphs and their symmetrically equivalent variants. We use first-principles electronic structure methods to identify minimum energy pathways and map them in subspaces of the symmetry-adapted strain order parameters. We next investigate the atomic structure of domain boundaries that separate coexisting variants of ferroelectric HfO 2 . Further, we rely on Gibbsian excess quantities and a precise specification of mechanical boundary conditions to describe the thermodynamic properties of domain boundaries. Our first-principles calculations show that the O and Hf shuffle arrangement within a domain boundary is closely related to the intermediate shuffle patterns of the homogeneous pathways between ferroelectric variants. Furthermore, the preferred structure within a boundary is very sensitive to local strain constraints imposed by the adjacent ferroelectric variants, leading to highly anisotropic domain boundary energies.

36 MATERIALS SCIENCE↗

Effect of frequency and applied voltage of an atmospheric-pressure dielectric-barrier discharge on breakdown and hydroxyl-radical generation with a liquid electrode

This work examines the effect of the frequency and peak applied voltage on hydroxyl-radical generation in a dielectric-barrier plasma discharge between a metallic needle electrode and one electrode covered with dielectric. The authors examine a system that can expose up to 96 liquid samples in an automated fashion without human intervention beyond setting the initial software configuration. Then, hydroxylradical concentration, measured through coumarin fluorescence, was measured for 5 s plasma exposures generated under different highvoltage conditions with frequencies from 2 to 16 kHz and amplitudes from 4 to 9 kV. Their results show that an increase in frequency and/or applied voltage, within the range prescribed above and the limits of the high-voltage power supply, can yield up to a 150% increase in fluorescence with an equivalent hydroxyl-radical increase. Applications using typical previous methods, such as the Fenton Reaction, are limited in that they continuously generate hydroxyl radicals over millisecond and longer intervals. These results establish the electrical parameters that can now be applied to polymers, like proteins, which show three dimensional structures that are flexible and fluctuate on a microsecond and nanosecond time scale, with hydroxyl-radical generation on this time scale using this device. Additionally, plasma exposures may be optimized for a great variety of proteins, devices and techniques, where hydroxyl-radical generation is of utmost importance, reducing exposure time and potential subjection of samples to harmful side effects.

42 ENGINEERING↗

Asynchronous Truncated Multigrid-Reduction-in-Time

In this paper, we present the new “asynchronous truncated multigrid-reduction-in-time” (AT-MGRIT) algorithm for introducing time parallelism to the solution of discretized time-dependent problems. The new algorithm is based on the multigrid-reduction-in-time (MGRIT) approach, which, in certain settings, is equivalent to another common multilevel parallel-in-time method, Parareal. In contrast to Parareal and MGRIT that both consider a global temporal grid over the entire time interval on the coarsest level, the AT-MGRIT algorithm uses truncated local time grids on the coarsest level, each grid covering certain temporal subintervals. Further, these local grids can be solved completely in an independent way from each other, which reduces the sequential part of the algorithm and, thus, increases parallelism in the method. Here, we study the effect of using truncated local coarse grids on the convergence of the algorithm, both theoretically and numerically, and show, using challenging nonlinear problems, that the new algorithm consistently outperforms classical Parareal/MGRIT in terms of time to solution.

97 MATHEMATICS AND COMPUTING↗