Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “second moment 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1924

The Influence of the Form of a Wooden Beam on Its Stiffness and Strength I : Deflection of Beams with Special Reference to Shear Deformations

The purpose of this investigation was to determine to what extent ordinary deflection formulas, which neglect shear deformations, are in error when applied to beams of various sections, and to develop reasonably accurate yet comparatively simple formulas which take into account such deformations. A great many tests were made to determine the amount of shear deformation for beams of various sections tested over many different spans. As the span over which the beam is tested is increased the error introduced by neglecting shear deformations becomes less, and the values obtained by substituting measured deflections in the ordinary formulas approach more nearly the modulus of elasticity in tension and compression. For short spans the error is considerable and increases rapidly as the span is reduced. Two formulas were developed for estimating the magnitude of shear deformations, both of which have been verified by tests. The first assumes the parabolic distribution of shear on a cross section of a beam and, starting with a differential volume, the distortion due to shear is determined by the ordinary methods of summarizing the work. The second assumes that the deflections due to shear in any two beams of the same length, height, and moment of inertia, which are similarly loaded, are proportional to the summations of the shear stresses on their respective vertical sections. Both formulas check experimental results very closely when the calculations are made with great refinement.

Newlin, J A↗
NASA NTRS1982

Static calibration for measurement of loads in wind tunnel tests on Shuttle surface tiles

A description is given of the procedures used for the static calibrations of the Coe (1981) sensor array to determine normal force, pitching moment, and rolling moment from the ensemble of sensor outputs. The use of independent data to confirm the validity of the calibration is also described. Two methods are used in calibrating the Coe sensors. In the first, calibration loads and readings are arranged in matrix form and are assumed to be related in a linear fashion. In the second, each Coe sensor reading is regarded as a sum of contributions from each of the calibration loads. It is shown that the Coe load sensors, in conjunction with either of two data analysis methods, can be used to determine the forces and moments experienced by a Shuttle TPS tile.

Rummler, D. R.↗
NASA NTRSDate not supplied

AMSR-2 Daily Snow Depth Data Product Using a Neural Network Algorithm Trained by Collocated ICESat-2 Measurements

By using diffusion theory and Monte Carlo lidar radiative transfer simulations, Hu 1 et al. (2022b) has derived snow depth from the first-, second- and third-order moments of the lidar backscattering pathlength distribution. Lu 2 et al. (2022,2024) calculated the snow depth by applying the methods to the satellite ICESat-2 lidar measurements over the Arctic sea ice, as well as land surfaces of Northern Hemisphere. In this paper, a neural network (NN) algorithm, employing several channels from AMSR-2 and the humidity vertical profiles GMAO GEOS-IT, is trained to determine snow depth identified by time and geolocation matched 2019 ICESat-2 snow-depth data during winter months over the Arctic sea ice. The trained NN snow-depth was applied to 2014-2020 AMSR-2 clear pixel data, although the algorithms perform reasonably well in thinner clouds. This paper used AMSR-2 data, a passive microwave instrument to generate a wide range of snow depth data, covering extensive spatial areas in the cross-orbit direction.

Neural Network↗
DOE OSTI2024

Efficient Unitary Designs from Random Sums and Permutations

A unitary k-design is an ensemble of unitaries that matches the first k moments of the Haar measure. In this work, we provide two efficient constructions of k-designs on n-qubits using new random matrix theory techniques. Our first construction is based on exponentiating sums of random i.i.d. Hermitian matrices and uses O(k2n2)-many gates. In the spirit of central limit theorems, we show that this random sum approximates the Gaussian Unitary Ensemble (GUE). We then show that the product of just two exponentiated GUE matrices is already approximately Haar random. Our second construction is based on products of exponentiated sums of random permutations and uses Õ(k poly (n)) many gates. The k dependence is optimal (up to polylogarithmic factors) and is inherited from the efficiency of existing k-wise independent permutations. Furthermore, replacing random permutations with quantum-secure pseudorandom permutations (PRPs), we also obtain a pseudorandom unitary (PRU) ensemble that is secure under nonadaptive queries. A central feature of both proofs is a new connection between the polynomial method in quantum query complexity and the large-dimension (N) expansion in random matrix theory. In particular, the first construction uses the polynomial method to control high moments of certain random matrix ensembles without requiring delicate Weingarten calculations. In doing so, we define and solve a moment problem on the unit circle, asking whether a finite number of equally weighted points can reproduce a given set of moments. In our second construction, the key step is to exhibit an orthonormal basis for irreducible representations of the partition algebra that has a low-degree large-N expansion. This allows us to show that the distinguishing probability is a low-degree rational polynomial of the dimension N.

algebra↗
NASA NTRS1989

Nonprincipal-plane scattering from flat plates: Second-order and corner diffraction and pattern control of horn antennas

Several high-frequency models for nonprincipal-plane scattering from a rectangular, perfectly conducting plate are examined. Two methods, the Method of Equivalent Currents and corner diffraction coefficients, are considered. Formulations for second-order Physical Theory of Diffraction equivalent currents and for corner diffracted fields are presented. Comparisons are made among plate models. Results away from grazing are accurate using only first-order terms. Near grazing, second-order and corner diffraction terms improve the results for many cases. The pattern control of horn antennas using lossy materials to coat the inner walls of the horn is also investigated. Integral Equation and Moment Method techniques are used to formulate the problem. It is clearly demonstrated that side lobe level reduction can be achieved using impedance surfaces on the inner walls of the horn.

Balanis, Constantine A.↗
NASA NTRS2024

Development of a Data Fusion Methodology for Lineload Aerodynamic Databases for a Launch Vehicle during Liftoff and Transition

The need for databases for the distributed loading on launch vehicles during the early portion of flight necessitates the use of expensive computational flows in regimes where wake effects dominate. While also being expensive, this is a regime that computational tools tend to historically have problems simulating accurately. To help tackle this problem, a method of data fusion to combine computational results to wind tunnel derived force and moment data is developed. Using this method, significant reduction in computational costs and increases in confidence of the final product is possible and has been used to generate several databases for the Space Launch System (SLS) at NASA. While the full details of database generation are not part of this work, the crucial method at its core is developed here. Two SLS geometries are used throughout the work to demonstrate the techniques. These are two of the larger geometries and represent both planned crewed missions to the Moon as well as potential cargo missions to deep space. The method uses principal component analysis (PCA) to generate a reduced ordered model (ROM) to help fill in the full parameter space. Other similar techniques are explored, but were not found to have a significant result on the predictions of the ROM. Because the full number of components are kept to generate the model, this lack of difference is expected. This method is then extended to ensure that predicted surfaces match trusted force and moment data derived from wind tunnel testing. This extension is done by setting up a constrained optimization problem in order to minimize the deviation from the surface resolved computational data while still integrating to the desired values. When generating the constrained optimization problem, a weighting factor to balance these competing needs is introduced. The work compares previously introduced weighting terms from similar work to the proposed terms and shows that the previously used terms do not have as desirable behavior in this flow regime. This method is then expanded by developing a technique to incorporate uncertainty quantification into the developed data fusion methodology. This expansion takes a two pronged approach. One examines transferring the uncertainties in the force and moment database and characterizes how those adjustments change the predicted lineloads. The second looks at model form error and looks how rebuilding the model using slightly different data changes the predictions. These two terms are then combined in order to create an uncertainty model that takes both effects into account. The limitations of the proposed methods is then discussed as well as possible techniques to address these shortcomings.

Launch Vehicles↗
DOE OSTI2026

Propagation of partially spatially coherent laser beams in instantaneous Kerr media

The propagation of intense, partially spatially coherent laser beams in a medium with instantaneous third-order susceptibility is studied analytically and numerically. For sufficiently high power relative to that required for nonlinear self-focusing, the propagation initially proceeds in two stages. In the first stage, spatial coherence builds up, and in the second stage, the number of speckles reduces. Once the degree of coherence is sufficiently high, whole-beam self-focusing occurs. The beam power is mostly confined within the initial spot radius. Two analytical approaches for describing the evolution of the beam are presented. The method of moments leads to an analytical solution for the rms spot radius that is in excellent agreement with simulations. This method does not require any knowledge of the field statistics beyond the initial conditions and provides no information about the evolution of the individual speckles. The other approach employs a self-similar solution for the second-order coherence function of the field and assumes that the fourth-order coherence function is factorizable and obeys complex circular Gaussian random statistics. The latter method also leads to an analytical expression for the spot radius, but its predictions for the qualitative evolution of the speckles disagree with wave-optics simulations.

lasers↗
NASA NTRS2003

A Stochastic Model of Space-Time Variability of Mesoscale Rainfall: Statistics of Spatial Averages

A characteristic feature of rainfall statistics is that they depend on the space and time scales over which rain data are averaged. A previously developed spectral model of rain statistics that is designed to capture this property, predicts power law scaling behavior for the second moment statistics of area-averaged rain rate on the averaging length scale L as L right arrow 0. In the present work a more efficient method of estimating the model parameters is presented, and used to fit the model to the statistics of area-averaged rain rate derived from gridded radar precipitation data from TOGA COARE. Statistical properties of the data and the model predictions are compared over a wide range of averaging scales. An extension of the spectral model scaling relations to describe the dependence of the average fraction of grid boxes within an area containing nonzero rain (the "rainy area fraction") on the grid scale L is also explored.

Kundu, Prasun K.↗
NASA NTRS1995

Out-of-plane properties

This paper summarizes three areas of research which were performed to characterize out-of-plane properties of composite materials. In the first investigation, a series of tests was run to characterize the through-the-thickness tensile strength for a variety of composites that included 2D braids, 2D and 3D weaves, and prepreg tapes. A new test method based on a curved beam was evaluated. Failures were significantly different between the 2D materials and the 3D weaves. The 2D materials delaminated between layers due to out-of-plane tensile stresses while the 3D weaves failed due to the formation of radial cracks between the surface plies caused by high circumferential stresses along the inner radius. The strength of the 2D textile composites did not increase relative to the tapes. Final failure in the 3D weaves was caused by a circumferential crack similar to the 2D materials and occurred at a lower bending moment than in other materials. The early failures in the 3D weaves were caused by radial crack formation rather than a low through-the-thickness strength. The second investigation focused on the development of a standard impact test method to measure impact damage resistance. The only impact tests that currently exist are compression after impact (CAI) tests which incorporate elements of both damage resistance and damage tolerance. A new impact test method is under development which uses a quasi-static indentation (QSI) test to directly measure damage resistance. Damage resistance is quantified in terms of the contact force to produce a unit of damage where a metric for damage may be area in C-scan, depth of residual dent , penetration, damage growth, etc. A final draft of an impact standard that uses a QSI test method will be presented to the ASTM Impact Task Group on impact. In the third investigation, the impact damage resistance behavior of a variety of textile materials was studied using the QSI test method. In this study, the force where large damage initiates was measured and the delamination size as a function of force was determined. The force to initiate large damage was significantly lower in braids and weaves. The delamination diameter - impact forace relationship was quanitfied using a damage resistance parameter, Q(*), which related delamination diameter to imapct force over a range of delamination sizes. Using this Q(*) parameter to rate the materials, the stitched uniweaves, toughened epoxy tapes, and through-the-thickness orthogonal interlock weave were the most damage resistant.

Jackson, Wade C.↗
NASA NTRS1991

Prediction of forces and moments for hypersonic flight vehicle control effectors

Developing methods of predicting flight control forces and moments for hypersonic vehicles, included a preliminary assessment of subsonic/supersonic panel methods and hypersonic local flow inclination methods for such predictions. While these findings clearly indicate the usefulness of such methods for conceptual design activities, deficiencies exist in some areas. Thus, a second phase of research was proposed in which a better understanding is sought for the reasons of the successes and failures of the methods considered, particularly for the cases at hypersonic Mach numbers. To obtain this additional understanding, a more careful study of the results obtained relative to the methods used was undertaken. In addition, where appropriate and necessary, a more complete modeling of the flow was performed using well proven methods of computational fluid dynamics. As a result, assessments will be made which are more quantitative than those of phase 1 regarding the uncertainty involved in the prediction of the aerodynamic derivatives. In addition, with improved understanding, it is anticipated that improvements resulting in better accuracy will be made to the simple force and moment prediction.

Maughmer, Mark D.↗
NASA NTRS1986

A second-order accurate kinetic-theory-based method for inviscid compressible flows

An upwind method for the numerical solution of the Euler equations is presented. This method, called the kinetic numerical method (KNM), is based on the fact that the Euler equations are moments of the Boltzmann equation of the kinetic theory of gases when the distribution function is Maxwellian. The KNM consists of two phases, the convection phase and the collision phase. The method is unconditionally stable and explicit. It is highly vectorizable and can be easily made total variation diminishing for the distribution function by a suitable choice of the interpolation strategy. The method is applied to a one-dimensional shock-propagation problem and to a two-dimensional shock-reflection problem.

Deshpande, Suresh M.↗
NASA NTRS1973

Statistical turbulence theory and turbulence phenomenology

The application of deductive turbulence theory for validity determination of turbulence phenomenology at the level of second-order, single-point moments is considered. Particular emphasis is placed on the phenomenological formula relating the dissipation to the turbulence energy and the Rotta-type formula for the return to isotropy. Methods which deal directly with most or all the scales of motion explicitly are reviewed briefly. The statistical theory of turbulence is presented as an expansion about randomness. Two concepts are involved: (1) a modeling of the turbulence as nearly multipoint Gaussian, and (2) a simultaneous introduction of a generalized eddy viscosity operator.

Herring, J. R.↗
NASA NTRS1989

Study of the integration of wind tunnel and computational methods for aerodynamic configurations

A study was conducted to determine the effectiveness of using a low-order panel code to estimate wind tunnel wall corrections. The corrections were found by two computations. The first computation included the test model and the surrounding wind tunnel walls, while in the second computation the wind tunnel walls were removed. The difference between the force and moment coefficients obtained by comparing these two cases allowed the determination of the wall corrections. The technique was verified by matching the test-section, wall-pressure signature from a wind tunnel test with the signature predicted by the panel code. To prove the viability of the technique, two cases were considered. The first was a two-dimensional high-lift wing with a flap that was tested in the 7- by 10-foot wind tunnel at NASA Ames Research Center. The second was a 1/32-scale model of the F/A-18 aircraft which was tested in the low-speed wind tunnel at San Diego State University. The panel code used was PMARC (Panel Method Ames Research Center). Results of this study indicate that the proposed wind tunnel wall correction method is comparable to other methods and that it also inherently includes the corrections due to model blockage and wing lift.

Browne, Lindsey E.↗
NASA NTRS1977

An optimal Space Shuttle ascent trajectory for the first orbital flight test

An optimal solution of the ascent trajectory of the Space Shuttle for the first orbital flight test is presented; the optimization is a minimum propellant, four-control problem in yaw angle, roll angle, pitch angle and vacuum thrust of each Space Shuttle main engine. Piecewise linear segments with juncture points treated as parameters are employed to model the controls. Equations of motion for a three-dimensional flight with pitch plane moment balance about an oblate are integrated numerically with a fourth-order Runge-Kutta method; two- and one-dimensional cubic spline function curve fits of aerodynamic coefficients are used during the first and second stages, respectively. The constraint minimization problem is solved with the Davidon-Fletcher-Powell function method.

Johnson, I. L., Jr.↗
NASA NTRS1993

A theoretical and numerical study of polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum

In this paper, theoretical and numerical results of the polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum are presented for the remote sensing of ocean and soil surfaces. The polarimetric scattered field for rough dielectric surfaces is derived to the second order by the small perturbation method (SPM). It is found that the second-order scattered field is coherent in nature, and its coefficients for different polarizations present the lowest-order corrections to the Fresnel reflection coefficients of the surfaces. In addition, the cross-polarized (HV and VH) components of the coherent fields are reciprocal and not zero for surfaces with anisotropic directional spectrum when the azimuth angle of the incident direction is not aligned with the symmetry directions of surfaces. In order to verify the energy conservation condition of the theoretical results, which is important if the theory is to be applied to the passive polarimetry of rough surfaces, a Monte Carlo simulation is performed to numerically calculate the polarimetric reflectivities of one-dimensional random rough surfaces which are generated with a prescribed power-law spectrum in the spectral domain and transformed to the spatial domain by the FFT. The surfaces simulated by this approach are periodic with the period corresponding to the low-wavenumber cutoff. To calculate the scattering from periodic dielectric surfaces, the authors present a new numerical technique which applies the Floquet theorem to reduce the problem to one period and does not require the evaluation of one-dimensional periodic Green's function used in the conventional method of moment formulation. Once the scattering coefficients are obtained, the polarimetric Stokes vectors for the emission from the random surfaces are then calculated according to the Kirchhoff's law and are illustrated as functions of relative azimuth observation and row directions. The second-order SPM is also reduced to this one-dimensional case, and theoretical results for all four Stokes parameters are in excellent agreement with the numerical results obtained from the Monte Carlo simulation. Moreover, the second order coherent fields are found to be necessary for the theoretical evaluation of the third and fourth Stokes parameters. Without that, the reflectivities of random rough surfaces are significantly over-estimated, and the signs of the third and the fourth Stokes parameters are found to be incorrect, if calculated by using the first-order incoherent fields alone.

Yueh, S. H.↗
DOE OSTI2025

Dark Energy Survey Year 3 results: $w$CDM cosmology from simulation-based inference with persistent homology on the sphere

We present cosmological constraints from Dark Energy Survey Year 3 (DES Y3) weak lensing data using persistent homology, a topological data analysis technique that tracks how features like clusters and voids evolve across density thresholds. For the first time, we apply spherical persistent homology to galaxy survey data through the algorithm TopoS2, which is optimized for curved-sky analyses and HEALPix compatibility. Employing a simulation-based inference framework with the Gower Street simulation suite, specifically designed to mimic DES Y3 data properties, we extract topological summary statistics from convergence maps across multiple smoothing scales and redshift bins. After neural network compression of these statistics, we estimate the likelihood function and validate our analysis against baryonic feedback effects, finding minimal biases (under $0.3σ$) in the $Ω_\mathrm{m}-S_8$ plane. Assuming the $w$CDM model, our combined Betti numbers and second moments analysis yields $S_8 = 0.821 \pm 0.018$ and $Ω_\mathrm{m} = 0.304\pm0.037$-constraints 70% tighter than those from cosmic shear two-point statistics in the same parameter plane. Our results demonstrate that topological methods provide a powerful and robust framework for extracting cosmological information, with our spherical methodology readily applicable to upcoming Stage IV wide-field galaxy surveys.

Prat, J. [Nordita; Royal Inst. Tech., Sodertalje; ↗
NASA NTRS1996

LES on unstructured deforming meshes: Towards reciprocating IC engines

A variable explicit/implicit characteristics-based advection scheme that is second-order accurate in space and time has been developed recently for unstructured deforming meshes (O'Rourke & Sahota 1996a). To explore the suitability of this methodology for Large-Eddy Simulation (LES), three subgrid-scale turbulence models have been implemented in the CHAD CFD code (O'Rourke & Sahota 1996b): a constant-coefficient Smagorinsky model, a dynamic Smagorinsky model for flows having one or more directions of statistical homogeneity, and a Lagrangian dynamic Smagorinsky model for flows having no spatial or temporal homogeneity (Meneveau et al. 1996). Computations have been made for three canonical flows, progressing towards the intended application of in-cylinder flow in a reciprocating engine. Grid sizes were selected to be comparable to the coarsest meshes used in earlier spectral LES studies. Quantitative results are reported for decaying homogeneous isotropic turbulence, and for a planar channel flow. Computations are compared to experimental measurements, to Direct-Numerical Simulation (DNS) data, and to Rapid-Distortion Theory (RDT) where appropriate. Generally satisfactory evolution of first and second moments is found on these coarse meshes; deviations are attributed to insufficient mesh resolution. Issues include mesh resolution and computational requirements for a specified level of accuracy, analytic characterization of the filtering implied by the numerical method, wall treatment, and inflow boundary conditions. To resolve these issues, finer-mesh simulations and computations of a simplified axisymmetric reciprocating piston-cylinder assembly are in progress.

Haworth, D. C.↗
DOE OSTI2025

Robust error calibration for serial crystallography

Serial crystallography is an important technique with unique abilities to resolve enzymatic transition states, minimize radiation damage to sensitive metalloenzymes and perform de novo structure determination from micrometre-sized crystals. This technique requires the merging of data from thousands of crystals, making manual identification of errant crystals unfeasible. cctbx.xfel.merge uses filtering to remove problematic data. However, this process is imperfect, and data reduction must be robust to outliers. We add robustness to cctbx.xfel.merge at the step of uncertainty determination for reflection intensities. This step is a critical point for robustness because it is the first step where the data sets are considered as a whole, as opposed to individual lattices. Robustness is conferred by reformulating the error-calibration procedure to have fewer and less stringent statistical assumptions and incorporating the ability to down-weight low-quality lattices. We then apply this method to five macromolecular XFEL data sets and observe the improvements to each. The appropriateness of the intensity uncertainties is demonstrated through internal consistency. This is performed through theoretical CC 1/2 and I /σ relationships and by weighted second moments, which use Wilson's prior to connect intensity uncertainties with their expected distribution. This work presents new mathematical tools to analyze intensity statistics and demonstrates their effectiveness through the often underappreciated process of uncertainty analysis.

Mittan-Moreau, David W.↗