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

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↗

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↗

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↗

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↗

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↗

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.↗

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.↗

Carrying the mass flux terms exactly in the first and second moment equations of compressible turbulence

In compressible turbulence models, it is assumed that the Favre-mean velocities are suitable approximations to the Reynolds-mean velocities in order to close unknown terms. This neglects, in the mean momentum and energy equations, the contribution to the stress and work terms by the mean of the fluctuating Favre velocity, a quantity proportional to the turbulent mass flux. As the stress and work terms do not introduce any new unknown correlations requiring closure in either k-epsilon or Reynolds stress closures and because the exact form of the terms can, with little additional work, be carried there is no need to make any modeling assumptions. In the Reynolds stress equations the viscous terms appear naturally in Reynolds variables while the problem is posed in Favre variables. In the process of splitting the viscous terms into the viscous transport terms, carried in Favre variables, and the dissipation terms, carried in Reynolds variables, important contributions from the mass flux appear. The accurate accounting of these terms is important for any consistent near wall modeling and the retention of the mass flux terms is important in complex compressible turbulent flows.

Ristorcelli, J. R.↗

Optimal compensator structure for linear time-invariant plant with inaccessible states

The problem is considered of designing an optimal linear time-invariant dynamic compensator for the regulation of an n-th order linear time-invariant plant with m independent outputs. The initial plant state is characterized by its first and second moments, and the cost is usual quadratic infinite-time penalty on the state and control, averaged over the initial plant and compensator states. The compensator is based on a minimal-order Luenberger observer and consequently has fixed dimension n-m. Necessary and sufficient conditions are derived for optimality of the compensator gains. The optimal compensator is shown to be unique if the plant has a particular canonical form, and, in general, for any arbitrary plant, the class of all optimal compensators is precisely determined.

Blanvillain, P. J. P.↗

Coupling Between Waveguide-Fed Slot Arrays

Coupling between two waveguide-fed planar slot arrays has been investigated using full-wave analysis. The analysis employs the method-of-moments solution to the pertinent coupled integral equations for the aperture electric field of all slots. In order to compute coupling between two arrays, the input port of the first array is excited with a TE(sub 10) mode wave while the second one is match-terminated. After solving the moment method matrix equations, the aperture fields of all slots are obtained and thereby the TE(sub 10) mode wave received at the input port of the second array is determined. Coupling between two arrays is the ratio of the wave amplitude arriving in the second array port to the incident wave amplitude at the first array port. The coupling mechanism has been studied as a function of spacing between arrays in different directions, e.g. the electric field plane, the magnetic field plane, and the diagonal plane. Computed coupling values are presented for different array geometries. This work is novel since it provides a good understanding of coupling between waveguide-fed slot arrays as a function of spacing and orientation for different aperture distributions and array architectures. This serves as a useful tool for antenna design engineers and system engineers.

Rengarajan, Sembiam↗

Modelling of Collision Induced Absorption Spectra Of H2-H2 Pairs for the Planetary Atmospheres Structure: The Second Overtone Band

The main objective of the proposal was to model the collision induced, second overtone band of gaseous hydrogen at low temperatures. The aim of this work is to assist planetary scientists in their investigation of planetary atmospheres, mainly those of Uranus and Neptune. The recently completed extended database of collision induced dipole moments of hydrogen pairs allowed us, for the first time, to obtain dipole moment matrix elements responsible for the roto-vibrational collision induced absorption spectra of H2-H2 in the second overtone band. Despite our numerous attempts to publish those data, the enormous volume of the database did not allow us to do this. Instead, we deposited the data on a www site. The final part of this work has been partially supported by NASA, Division for Planetary Atmospheres. In order to use our new data for modelling purpose, we first needed to test how well we can reproduce the existing experimental data from theory, when using our new input data. Two papers resulted from this work. The obtained agreement between theoretical results and the measurements appeared to be within 10-30%. The obviously poorer agreement than observed for the first H2 overtone, the fundamental, and the rototranslational bands can be attributed to the fact that dipole moments responsible for the second overtone are much weaker, therefore susceptible to larger numerical uncertainties. At the same time, the intensity of the second overtone band is much weaker and therefore it is much harder to be measured accurately in the laboratory. We need to point out that until now, no dependable model of the 2nd overtone band was available for modelling of the planetary atmospheres. The only one, often referred to in previous works on Uranian and Neptune's atmospheres, uses only one lineshape, with one (or two) parameter(s) deduced at the effective temperature of Uranus (by fitting the planetary observation). After that, the parameter(s) was(were) made temperature dependent according to some very simple relation. Summarizing, no reliable temperature-dependent model has been available yet. Our approach was a bit different from similar attempts done earlier, on account of the poorer agreement of theory with experiment. We needed to resort to some semi-empirical procedure. While we were in a favourable position to be able to rely on the physical input data, these, apparently, did not supply the most dependable predictions (simply because the results did not agree well enough with experimental data). On the other hand, the relative deviations between the theory and experiment were comparable at 77 and at 298 K. That fact indicated that theory is capable of predicting the temperature dependence of the absorption spectra well. We have thus chosen the "middle way". We have fitted the existing measurements with many 3- parameter lineshapes, in order to achieve the closest fit.

Borysow, Aleksandra↗

Sensitivity derivatives for advanced CFD algorithm and viscous modelling parameters via automatic differentiation

The computational technique of automatic differentiation (AD) is applied to a three-dimensional thin-layer Navier-Stokes multigrid flow solver to assess the feasibility and computational impact of obtaining exact sensitivity derivatives typical of those needed for sensitivity analyses. Calculations are performed for an ONERA M6 wing in transonic flow with both the Baldwin-Lomax and Johnson-King turbulence models. The wing lift, drag, and pitching moment coefficients are differentiated with respect to two different groups of input parameters. The first group consists of the second- and fourth-order damping coefficients of the computational algorithm, whereas the second group consists of two parameters in the viscous turbulent flow physics modelling. Results obtained via AD are compared, for both accuracy and computational efficiency with the results obtained with divided differences (DD). The AD results are accurate, extremely simple to obtain, and show significant computational advantage over those obtained by DD for some cases.

Green, Lawrence L.↗

Spectral Moments of Collision-Induced Absorption of CO2 Pairs: The Role of the Intermolecular Potential

In this paper we examine the role of the anisotropy of the intermolecular potential in the rototranslational collision-induced absorption of the CO2 pairs. Using newly developed formulas that include the effects of anisotropy of the potential to all orders, we calculate the two lowest spectral moments gamma(prime), and alpha(prime), for four different classes of C02 pair potentials and compare the results with the experimental values. We assumed only multipolar induction in the process of forming the induced dipole, with the second-order contributions included. Using a site-site LJ and a site-site semi-ab initio intermolecular potentials we were able to reproduce the experimental values of gamma(prime), and alpha(prime) moments over entire temperature range from 230 to 330 K. Also, the role of an electrostatic interaction between two C02 molecules and its impact on the spectral moments is thoroughly investigated. An isotropic core with a point quadrupole centered at each molecule is shown to be an inadequate representation of the C02-CO2 potential. Additionally, we show the results obtained with the first- and second-order perturbation theory to be more than twice too small.

Gruszka, Marcin↗

Numerical solution of the Navier-Stokes equations with topography

A finite difference scheme for solving the equations of fluid motion in a generalized coordinate system has been constructed. The scheme conserves mass and all the first integral moments of the motion. The scheme also advectively 'almost conserves' second moments, in that the magnitude of implicit numerical smoothing is typically about an order smaller than explicit viscosity and diffusion. Calculations with the model support the theoretical conjecture that the difference scheme is stable whenever the analogous Cartesian scheme is stable. The scheme has been used to calculate dry atmospheric convection due to differential heating between top and bottom of mountainous terrain. The general small-scale characteristics of mountain up-slope winds have been simulated. In addition, the results have demonstrated the crucial role played by the eddy diffusivities and the environmental stability, in determining both the quantitative and the qualitative features of the circulation.

Gal-Chen, T.↗

Analyses of convective flows and spatial gradients in energetic ion observations

Based on the measured anisotropy of three-dimensional energetic ion distributions obtained from the Medium Energy Particle Instrument experiment of the ISEE 1 spacecraft, bulk flow, epsilon x B convection, and pressure gradients in magnetospheric plasmas are measured with two techniques. In the first technique, using a Compton-Getting transform into the appropriate rest frame in an iterative fashion and converging to the unique velocity that removes the first-order anisotropy from the flux angular distribution, the actual measured angular distribution, and spectrum at each angle, are employed. The second technique uses a first moment of the distribution function in which the shift in velocity space of the distribution function is obtained at two phase space density values. The convection contribution to the velocity shift is assumed to be perpendicular to the pressure gradient contribution if an MHD generator condition is assumed, rendering the two shifts separable, and a gradient scale and a lower limit to the convection electric field are derived in one calculation.

Mitchell, D. G.↗

High-Order Space-Time Methods for Conservation Laws

Current high-order methods such as discontinuous Galerkin and/or flux reconstruction can provide effective discretization for the spatial derivatives. Together with a time discretization, such methods result in either too small a time step size in the case of an explicit scheme or a very large system in the case of an implicit one. To tackle these problems, two new high-order space-time schemes for conservation laws are introduced: the first is explicit and the second, implicit. The explicit method here, also called the moment scheme, achieves a Courant-Friedrichs-Lewy (CFL) condition of 1 for the case of one-spatial dimension regardless of the degree of the polynomial approximation. (For standard explicit methods, if the spatial approximation is of degree p, then the time step sizes are typically proportional to 1/p(exp 2)). Fourier analyses for the one and two-dimensional cases are carried out. The property of super accuracy (or super convergence) is discussed. The implicit method is a simplified but optimal version of the discontinuous Galerkin scheme applied to time. It reduces to a collocation implicit Runge-Kutta (RK) method for ordinary differential equations (ODE) called Radau IIA. The explicit and implicit schemes are closely related since they employ the same intermediate time levels, and the former can serve as a key building block in an iterative procedure for the latter. A limiting technique for the piecewise linear scheme is also discussed. The technique can suppress oscillations near a discontinuity while preserving accuracy near extrema. Preliminary numerical results are shown

Huynh, H. T.↗

Near surface analysis

This is a study to assist in the understanding of earth near surface structure. Higher order moments are used to detect the density distribution as well as to seek patterns found in geological structures. It is shown how higher order moments at points outside a mass structure are determined as well as how to recover the mass distribution from the higher order moments. It is interesting to note that the first moment at a point P outside the mass structure, V(sub O)(P), is the entire mass and the second moment, V(sub 1)(P), is the potential at P due to the mass structure. Usually only the mass and the potential function are used to determine the density distribution in a body. An infinite function sequence (V(sub n)(P))(sub n=0)(sup infinity) is required to uniquely determine the density distribution.

Johnson, Gordon G.↗

Progress in wind tunnel wall interference assessment/correction procedures at the NAE

Wall corrections as a function of wall porosity in the transonic wall interference problem was assessed. Effective porosities primarily for the two dimensional case were established as follows: (1) comparison of experimental data for two geometrically similar models of different chord/height ratio, an overall value of wall porosity could be deduced; (2) theoretical development which allows for unequal porosity for the floor and ceiling and wall boundary pressure measurements, porosities for floor and ceiling could be deduced; (3) a scheme was developed which allowed unequal porosity of floor and ceiling and streamwise varying porosity. The boundary layer development along the perforated floor and ceiling under the influence of the model pressure field, variations in boundary layer thickness underlining the difficulties in deducing meaningful values of wall porosity were determined. Wall boundary pressure measurement, in combination with singularity modelling of the airfoil, was sufficient to yield required information on the wall interference flow without having to establish some value for wall porosity. The singularity modelling of the airfoil initially covered only lift and volume but was extended to include drag and pitching moment, and second order volume term. It is shown by asymptotic transonic small disturbance analysis, that the derived corrections to angle of attack and free stream Mach number are correct to the first order.

Ohman, L. H.↗