Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “First Order Second Moment”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

On the Accuracy of Probabilistic Bucking Load Prediction

The buckling strength of thin-walled stiffened or unstiffened, metallic or composite shells is of major concern in aeronautical and space applications. The difficulty to predict the behavior of axially compressed thin-walled cylindrical shells continues to worry design engineers as we enter the third millennium. Thanks to extensive research programs in the late sixties and early seventies and the contributions of many eminent scientists, it is known that buckling strength calculations are affected by the uncertainties in the definition of the parameters of the problem such as definition of loads, material properties, geometric variables, edge support conditions, and the accuracy of the engineering models and analysis tools used in the design phase. The NASA design criteria monographs from the late sixties account for these design uncertainties by the use of a lump sum safety factor. This so-called 'empirical knockdown factor gamma' usually results in overly conservative design. Recently new reliability based probabilistic design procedure for buckling critical imperfect shells have been proposed. It essentially consists of a stochastic approach which introduces an improved 'scientific knockdown factor lambda(sub a)', that is not as conservative as the traditional empirical one. In order to incorporate probabilistic methods into a High Fidelity Analysis Approach one must be able to assess the accuracy of the various steps that must be executed to complete a reliability calculation. In the present paper the effect of size of the experimental input sample on the predicted value of the scientific knockdown factor lambda(sub a) calculated by the First-Order, Second-Moment Method is investigated.

Arbocz, Johann↗
DOE OSTI2020

Dynamics of Subsiding Shells in Actively Growing Clouds with Vertical Updrafts

The dynamics of a subsiding shell at the edges of actively growing shallow cumulus clouds with updrafts is analyzed using direct numerical simulation. The actively growing clouds have a fixed in-cloud buoyancy and velocity. Turbulent mixing and evaporative cooling at the cloud edges generate a subsiding shell that grows with time. A self-similar regime is observed for first- and second-order moments when normalized with respective maximum values. Internal scales derived from integral properties of the flow problem are identified. A self-similarity analysis using these scales reveals that contrary to classical self-similar flows, the turbulent kinetic energy budget terms and velocity moments scale according to the buoyancy and not with the mean velocity. The shell thickness is observed to increase linearly with time. The buoyancy scale remains time invariant and is set by the initial cloud–environment thermodynamics. The shell accelerates ballistically with a magnitude set by the saturation value of the buoyancy of the cloud–environment mixture. In this regime, the shell is buoyancy driven and independent of the in-cloud velocity. Relations are obtained for predicting the shell thickness and minimum velocities by linking the internal scales with external flow parameters. The values thus calculated are consistent with the thickness and velocities observed in typical shallow cumulus clouds. The entrainment coefficient is a function of the initial state of the cloud and the environment, and is shown to be on the same order of magnitude as fractional entrainment rates calculated for large-scale models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1972

Nonlinear estimation theory applied to orbit determination

The development of an approximate nonlinear filter using the Martingale theory and appropriate smoothing properties is considered. Both the first order and the second order moments were estimated. The filter developed can be classified as a modified Gaussian second order filter. Its performance was evaluated in a simulated study of the problem of estimating the state of an interplanetary space vehicle during both a simulated Jupiter flyby and a simulated Jupiter orbiter mission. In addition to the modified Gaussian second order filter, the modified truncated second order filter was also evaluated in the simulated study. Results obtained with each of these filters were compared with numerical results obtained with the extended Kalman filter and the performance of each filter is determined by comparison with the actual estimation errors. The simulations were designed to determine the effects of the second order terms in the dynamic state relations, the observation state relations, and the Kalman gain compensation term. It is shown that the Kalman gain-compensated filter which includes only the Kalman gain compensation term is superior to all of the other filters.

Choe, C. Y.↗
NASA NTRS1990

Estimation of multidimensional precipitation parameters by areal estimates of oceanic rainfall

The parameters of the multidimensional precipitation model proposed by Waymire et al. (1984) are estimated using the areal-averaged radar measurements of precipitation of the Global Atlantic Tropical Experiment (GATE) data set. The procedure followed was the fitting of the first- and second-order moments at different aggregation scales by nonlinear regression techniques. The numerical estimates of the parameters using different subsets of GATE information were reasonably stable, i.e., they were not affected by changes of the area-averaging size, temporal length of the records, and percentage of areal coverage of rainfall. This suggests that the estimation procedure is relatively robust and suitable to estimate the parameters of the multidimensional model in areas of sparse density of rain gages. The use of the space-time spectrum of rainfall to help in the determination of sampling errors due to intermittent visits of future space-borne low-altitude sensors of precipitation is also discussed.

Valdes, J. B.↗
DOE OSTI2021

Lagrangian analysis for turbulent transport in variable-density turbulence

In this report Lagrangian analysis of materially conserved scalars is applied to the problem of turbulent transport in variable-density flows. The consequences of an additional material conserved quantity, the density, is generally not acknowledged and leads to significant and meaningfully different expressions for turbulent transport in the moment equations. The formal Lagrangian analysis produces gradient transport expressions substantially different from those obtained by the physically intuitive “argument by analogy” method used in computational models. Various intuitive arguments, in Favre and Reynolds averaged settings, are contrasted to the formal Lagrangian results. Using expressions from the formal analysis, we derive consistent gradient transport closures for the turbulent transport terms that appear in the first- and second-order Favre moment equations. Results for coupled multispecies turbulent transport are given. The analysis is limited to variable-density turbulence in which the dilatation of the fluctuating velocity is small. The results are applicable to turbulent combustion and to stellar convection problems in which the density fluctuations are on the order of the mean density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
NASA NTRSDate not supplied

Snow Depth from AMSR-2 Using Multispectral Satellite Data in an Artificial Neural Network

By using diffusion theory and Monte Carlo lidar radiative transfer simulations, Hu et al. (2022b) has derived snow depth from the first-, second- and third-order moments of the lidar backscattering pathlength distribution. Lu et al. (2022) 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, an artificial neural network (ANN) algorithm, employing several channels from Advanced Microwave Scanning Radiometer 2 (AMSR-2) and the humidity vertical profiles from Global Modeling and Assimilation Office (GMAO) Goddard Earth Observing System for Instrument Teams (GEOS-IT) product, 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 ANN snow-depth was applied to 2018 AMSR-2 clear pixel data, although the algorithms perform reasonably well in thinner clouds. The validation data (different from the training set) of ANN snow depth from AMSR-2 showed a good agreement with time matched and co-located snow-depth values from ICESat-2. The bias was near zero, with mean absolute error (MAE) 0.05 cm and a root-mean-square-error (RMSE) 0.08 cm. Prior applying the trained ANN snow depth to AMSR-2 data, a cloud screening algorithm was developed with a similar approach. A separate ANN cloud mask was trained to determine an AMSR-2 pixel is clear or cloudy with time and geolocation matched 2015 CALIOP Vertical Feature Mask (VFM) over Arctic sea ice. The ANN cloud mask from AMSR-2 under-estimated cloud fraction by 3-6% compared to CALIOP . The additional research is needed to conclusively evaluate the ANN cloud mask accuracy. Finally, this paper will lay the foundation for a sustained long-term snowfall and snow-storm monitoring system. The future Cloud Aerosol LIdar for Global scale Observations of the ocean-Land Atmosphere system (CALIGOLA) mission will provide a means to calculate snow depth from the lidar backscattering pathlength distribution, benefiting from the UV, visible and infrared pulses. With the calculated snow depth as the truth one could develop a machine learning algorithm, as it was done in this paper, using a passive microwave instrument available at that time to generate a wide range of snow depth data, covering extensive spatial areas in the cross-orbit direction.

Neural Network↗
DOE OSTI2026

Second Order Closures for the Radiative Transfer Equation: Some Are Unstable

The largest existing simulations of cosmic reionization model radiative transfer with moment methods that require a closure relation. The two most commonly used closure relations are M1 and OTVET; both close the moment hierarchy at the first moment. We explore the properties of a higher, second-order closure. We show that direct generalizations of M1 and OTVET to one higher order are physically unstable - i.e., the closure equations themselves result in unstable solutions, not just their numerical implementation. In fact, a generalization of OTVET to any order higher than the first one is unstable. We are also able to show that any local (i.e., depending only on the local moments of the radiation field, like M1) second-order closure that depends only on the radiation intensity and radiation flux, but does not explicitly depend on the radiation pressure, is physically unstable. This result restricts the choice of possible second-order closure relations.

Gnedin, Nickolay Y. [Fermilab; Chicago U., Astron.↗
NASA NTRS1996

The Hot Bands of Methane Between 5 and 10(micro)m

Experimental line intensities of 1727 transitions arising from nine hot bands in the pentad-dyad system of methane are fitted to first and second order using the effective dipole moment expansion in the polyad scheme. The observed bands are V3-V2, V3-V4, V1-V2, V1-V4, 2V4-V4, V2+V4-V2, V2+V4-V4, 2V2-V4, and the intensities are obtained from long-path spectra recorded with the Fourier transform spectrometer located at Kitt Peak National Observatory. For the second order model, some of the 27 intensity parameters are not linearly independant, and so two methods (extrapolation and effective parameters) are propopsed to model the intensities of the hot bands.

methane hot bands dipole moment expansion Kitt Pea↗
NASA NTRS1998

The Real-Time Wall Interference Correction System of the NASA Ames 12-Foot Pressure Wind Tunnel

An improved version of the Wall Signature Method was developed to compute wall interference effects in three-dimensional subsonic wind tunnel testing of aircraft models in real-time. The method may be applied to a full-span or a semispan model. A simplified singularity representation of the aircraft model is used. Fuselage, support system, propulsion simulator, and separation wake volume blockage effects are represented by point sources and sinks. Lifting effects are represented by semi-infinite line doublets. The singularity representation of the test article is combined with the measurement of wind tunnel test reference conditions, wall pressure, lift force, thrust force, pitching moment, rolling moment, and pre-computed solutions of the subsonic potential equation to determine first order wall interference corrections. Second order wall interference corrections for pitching and rolling moment coefficient are also determined. A new procedure is presented that estimates a rolling moment coefficient correction for wings with non-symmetric lift distribution. Experimental data obtained during the calibration of the Ames Bipod model support system and during tests of two semispan models mounted on an image plane in the NASA Ames 12 ft. Pressure Wind Tunnel are used to demonstrate the application of the wall interference correction method.

Ulbrich, Norbert↗
NASA NTRS1992

Polarimetric signatures of a canopy of dielectric cylinders based on first and second order vector radiative transfer theory

Complete polarimetric signatures of a canopy of dielectric cylinders overlying a homogeneous half space are studied with the first and second order solutions of the vector radiative transfer theory. The vector radiative transfer equations contain a general nondiagonal extinction matrix and a phase matrix. The energy conservation issue is addressed by calculating the elements of the extinction matrix and the elements of the phase matrix in a manner that is consistent with energy conservation. Two methods are used. In the first method, the surface fields and the internal fields of the dielectric cylinder are calculated by using the fields of an infinite cylinder. The phase matrix is calculated and the extinction matrix is calculated by summing the absorption and scattering to ensure energy conservation. In the second method, the method of moments is used to calculate the elements of the extinction and phase matrices. The Mueller matrix based on the first order and second order multiple scattering solutions of the vector radiative transfer equation are calculated. Results from the two methods are compared. The vector radiative transfer equations, combined with the solution based on method of moments, obey both energy conservation and reciprocity. The polarimetric signatures, copolarized and depolarized return, degree of polarization, and phase differences are studied as a function of the orientation, sizes, and dielectric properties of the cylinders. It is shown that second order scattering is generally important for vegetation canopy at C band and can be important at L band for some cases.

Tsang, Leung↗
DOE OSTI2024

On the statistical theory of self-gravitating collisionless dark matter flow: High order kinematic and dynamic relations

Dark matter, if it exists, accounts for five times as much as ordinary baryonic matter. To better understand the self-gravitating collisionless dark matter flow on different scales, a statistical theory involving kinematic and dynamic relations must be developed for different types of flow, e.g., incompressible, constant divergence, and irrotational flow. This is mathematically challenging because of the intrinsic complexity of dark matter flow and the lack of a self-closed description of flow velocity. Here, this paper extends our previous work on second-order statistics Xu to kinematic relations of any order for any type of flow. Dynamic relations were also developed to relate statistical measures of different orders. The results were validated by N-body simulations. On large scales, we found that (i) third-order velocity correlations can be related to density correlation or pairwise velocity; (ii) the pth-order velocity correlations follow ∝ a (p+2)/2 for odd p and ∝ a p/2 for even p, where a is the scale factor; (iii) the overdensity δ is proportional to density correlation on the same scale, $\langle$δ$\rangle$∝$\langle$δδ'$\rangle$; (iv) velocity dispersion on a given scale r is proportional to the overdensity on the same scale. On small scales, (i) a self-closed velocity evolution is developed by decomposing the velocity into motion in haloes and motion of haloes; (ii) the evolution of vorticity and enstrophy are derived from the evolution of velocity; (iii) dynamic relations are derived to relate second- and third-order correlations; (iv) while the first moment of pairwise velocity follows $\langle$Δu L $\rangle$=-Har (H is the Hubble parameter), the third moment follows $\langle$(Δu L ) 3 $\rangle$ ∝ ε u ar that can be directly compared with simulations and observations, where ε u ≈ 10 -7 m 2 /s 3 is the constant rate for energy cascade; (v) the pth order velocity correlations follow ∝ a (3p-5)/4 for odd p and ∝ a 3p/4 for even p. Finally, the combined kinematic and dynamic relations lead to exponential and one-fourth power-law velocity correlations on large and small scales, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Second-order closure PBL model with new third-order moments: Comparison with LES data

This paper contains two parts. In the first part, a new set of diagnostic equations is derived for the third-order moments for a buoyancy-driven flow, by exact inversion of the prognostic equations for the third-order moment equations in the stationary case. The third-order moments exhibit a universal structure: they all are a linear combination of the derivatives of all the second-order moments, bar-w(exp 2), bar-w theta, bar-theta(exp 2), and bar-q(exp 2). Each term of the sum contains a turbulent diffusivity D(sub t), which also exhibits a universal structure of the form D(sub t) = a nu(sub t) + b bar-w theta. Since the sign of the convective flux changes depending on stable or unstable stratification, D(sub t) varies according to the type of stratification. Here nu(sub t) approximately equal to wl (l is a mixing length and w is an rms velocity) represents the 'mechanical' part, while the 'buoyancy' part is represented by the convective flux bar-w theta. The quantities a and b are functions of the variable N(sub tau)(exp 2), where N(exp 2) = g alpha derivative of Theta with respect to z and tau is the turbulence time scale. The new expressions for the third-order moments generalize those of Zeman and Lumley, which were subsequently adopted by Sun and Ogura, Chen and Cotton, and Finger and Schmidt in their treatments of the convective boundary layer. In the second part, the new expressions for the third-order moments are used to solve the ensemble average equations describing a purely convective boundary laye r heated from below at a constant rate. The computed second- and third-order moments are then compared with the corresponding Large Eddy Simulation (LES) results, most of which are obtained by running a new LES code, and part of which are taken from published results. The ensemble average results compare favorably with the LES data.

Canuto, V. M.↗
DOE OSTI2021

Field moment expansion method for interacting bosonic systems

In this work, we introduce a numerical method and python package, https://github.com/andillio/CHiMES, that simulates quantum systems initially well approximated by mean field theory using a second order extension of the classical field approach. We call this the field moment expansion method. In this way, we can accurately approximate the evolution of first and second field moments beyond where the mean field theory breaks down. This allows us to estimate the quantum breaktime of a classical approximation without any calculations external to the theory. We investigate the accuracy of the field moment expansion using a number of well studied quantum test problems. Interacting Bosonic systems similar to scalar field dark matter are chosen as test problems. We find that successful application of this method depends on two conditions: the quantum system must initially be well described by the classical theory, and that the growth of the higher order moments be hierarchical.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
DOE OSTI2021

Harmonic Infrared and Raman Spectra in Molecular Environments Using the Polarizable Embedding Model

We present a fully analytic approach to calculate infrared (IR) and Raman spectra of molecules embedded in complex molecular environments modeled using the fragment-based polarizable embedding (PE) model. We provide the theory for the calculation of analytic second-order geometric derivatives of molecular energies and first-order geometric derivatives of electric dipole moments and dipole–dipole polarizabilities within the PE model. The derivatives are implemented using a general open-ended response theory framework, thus allowing for an extension to higher-order derivatives. The embedding-potential parameters used to describe the environment in the PE model are derived through first-principles calculations, thus allowing a wide variety of systems to be modeled, including solvents, proteins, and other large and complex molecular environments. Here, we present proof-of-principle calculations of IR and Raman spectra of acetone in different solvents. This work is an important step toward calculating accurate vibrational spectra of molecules embedded in realistic environments.

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

Employing Sensitivity Derivatives to Estimate Uncertainty Propagation in CFD

Two methods that exploit the availability of sensitivity derivatives are successfully employed to predict uncertainty propagation through Computational Fluid Dynamics (CFD) code for an inviscid airfoil problem. An approximate statistical second-moment method and a Sensitivity Derivative Enhanced Monte Carlo (SDEMC) method are successfully demonstrated on a two-dimensional problem. First- and second-order sensitivity derivatives of code output with respect to code input are obtained through an efficient incremental iterative approach. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables); these sensitivity derivatives enable one to formulate first- and second-order Taylor Series approximations for the mean and variance of CFD output quantities. Additionally, incorporation of the first-order sensitivity derivatives into the data reduction phase of a conventional Monte Carlo (MC) simulation allows for improved accuracy in determining the first moment of the CFD output. Both methods are compared to results generated using a conventional MC method. The methods that exploit the availability of sensitivity derivatives are found to be valid when considering small deviations from input mean values.

Putko, Michele M.↗