Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “normalization 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

Normal and radial impact of composites with embedded penny-shaped cracks

A method is developed for the dynamic stress analysis of a layered composite containing an embedded penny-shaped crack and subjected to normal and radial impact. The material properties of the layers are chosen such that the crack lies in a layer of matrix material while the surrounding material possesses the average elastic properties of a two-phase medium consisting of a large number of fibers embedded in the matrix. Quantitatively, the time-dependent stresses near the crack border can be described by the dynamic stress intensity factors. Their magnitude depends on time, on the material properties of the composite and on the relative size of the crack compared to the composite local geometry. Results obtained show that, for the same material properties and geometry of the composite, the dynamic stress intensity factors for an embedded (penny-shaped) crack reach their peak values within a shorter period of time and with a lower magnitude than the corresponding dynamic stress intensity factors for a through-crack.

Sih, G. C.↗

Reliability of the finite element method for calculating free edge stresses in composite laminates

The interlaminar normal stress distributions along the interface between the +45 deg and -45 deg plies of a graphite/epoxy laminate, obtained by various investigators, were found to disagree in both magnitude and sign. The reliability of the displacement-formulated finite element method in analyzing the edge-stress problem of a composite laminate is investigated. The history of the edge-stress problem is reviewed, and two well-known elasticity problems, one involving a stress discontinuity and one a singularity, are analyzed. The finite element analysis in these problems yields accurate stress distributions everywhere except in two elements closest to the stress discontinuity or singularity. Stress distributions for a + or -45 deg ply laminate near the singularity were similar to those of the two elasticity problems, demonstrating the methods, accuracy for calculating interlaminar stresses in composite laminates. The disagreement between the numerical methods was attributed to the unsymmetric stress tensor at singularity.

Whitcomb, J. D.↗

Far-infrared luminosity functions of normal galaxies

A volume-limited sample is constructed from the Zwicky catalog and IRAS data base to examine the FIR luminosity functions of normal galaxies, and to investigate possible relationships between FIR emission and galaxy morphology. Quantitative and unbiased treatment is provided by 'survival analysis' statistical methods. It is found that the FIR distributions of normal galaxies are better fit by lognormal than Schechter functions. The total FIR emissivity (8 to 115 microns) of normal galaxies is approximately equal to half their emission in the B plus V optical bands. Normal galaxy FIR emission is uncorrelated with the basic S0-Sm Hubble sequence of spiral galaxy morphology, but appears to be affected by de Vaucouleurs' (1959) revised morphological classifications based on inner rings and S-shaped arms. Spirals with bars and inner rings are systematically fainter than unbarred spirals. It is suggested that bars and rings reduce the amount or spatially confine the dust in spiral disks, resulting in lower efficiency conversion of optical and UV photons into the IR.

Isobe, Takashi↗

Simulation-driven optimization of high-order meshes in ALE hydrodynamics

Here we propose tools for high-order mesh optimization and demonstrate their benefits in the context of multi-material Arbitrary Lagrangian-Eulerian (ALE) compressible shock hydrodynamic applications. The mesh optimization process is driven by information provided by the simulation which uses the optimized mesh, such as shock positions, material regions, known error estimates, etc. These simulation features are usually represented discretely, for instance, as finite element functions on the Lagrangian mesh. The discrete nature of the input is critical for the practical applicability of the algorithms we propose and distinguishes this work from approaches that strictly require analytical information. Our methods are based on node movement through a high-order extension of the Target-Matrix Optimization Paradigm (TMOP). The proposed formulation is fully algebraic and relies only on local Jacobian matrices, so it is applicable to all types of mesh elements, in 2D and 3D, and any order of the mesh. We discuss the notions of constructing adaptive target matrices and obtaining their derivatives, reconstructing discrete data in intermediate meshes, node limiting that enables improvement of global mesh quality while preserving space-dependent local mesh features, and appropriate normalization of the objective function. The adaptivity methods are combined with automatic ALE triggers that can provide robustness of the mesh evolution and avoid excessive remap procedures. The benefits of the new high-order TMOP technology are illustrated on several simulations performed in the high-order ALE application BLAST.

97 MATHEMATICS AND COMPUTING↗

Graph neural networks for mechanical property prediction of 2D fiber composites

This work investigates the ability of graph neural networks (GNNs) to homogenize 2D fiber composite microstructures. We use different inhomogeneity and anisotropy indices to motivate and show that the Volume Elements (VEs) used in ML methods should ideally be far from their Representative Volume Element (RVE) size limit and, consequently, are notably anisotropic. Hence, training only the isotropic limit properties may not be acceptable. Another aspect is the need to normalize elastic stiffness values for ML, especially when high elastic contrast ratios are encountered between composite phases or in the material set. We introduce a normalization technique based on the mean-field method (MFM) to handle such high contrast ratios and train for the entire stiffness tensor. We show that the proposed GNN approaches exhibit high accuracy and efficiency compared to traditional methods and convolutional neural networks, utilizing unstructured graphs constructed from microstructure topology. Our model successfully predicts the stiffness tensor, peak strength under bulk damage, and brittle fracture initiation strength across diverse microstructure configurations while maintaining high accuracy even for extreme material contrasts and volume fractions. We also present a method to improve prediction accuracy for small dataset sizes using Voronoi partitioning.

Brittle strength↗

Cross Calibration of TOMS, SBUV/2 and SCIAMACHY Radiances from Ground Observations

We have shown that validation of radiances is a very effective means for correcting absolute accuracy and long term drifts of backscatter type satellite measurements. This method by-passes the algorithms used for both satellite and ground based measurements which are normally used to validate and correct the satellite data. A new method for satellite validation is planned which will compliment measurements from the existing ground-based networks. This method will employ very accurate comparisons between ground based zenith sky radiances and satellite nadir radiances. These comparisons will rely heavily on the experience derived from the Shuttle SBUV (SSBUV) program which provided a reference standard of radiance measurements for SBUV/2, TOMS, and GOME. This new measurement program, called 'Skyrad', employs two well established capabilities at the Goddard Space Flight Center, 1) the SSBUV calibration facilities and 2) the radiative transfer codes used for the TOMS and SBUV/2 algorithms and their subsequent refinements. Radiative transfer calculations show that ground based zenith sky and satellite nadir backscatter ultraviolet comparisons can be made very accurately under certain viewing conditions. The Skyrad instruments (SSBUV, Brewer spectrophotometers, and possibly others) will be calibrated and maintained to a precision of a few tenths of a percent. Skyrad data will then enable long term calibration of upcoming satellite instruments such as QuickTOMS, SBUV/2s and SCIAMACHY with a high degree of precision. This technique can be further employed to monitor the performance of future instruments such as GOMEZ, OMI, and OMPS. Additional information is included in the original extended abstract.

Hilsenrath, Ernest↗

Monthly Mean DNI and GTI Derived from Monthly Mean GHI and DHI Using Two Methods: Comparisons with the BSRN Data and the Results Derived from the CERES Hourly Data

Monthly mean Global Horizontal Irradiances (GHI) and Diffuse Horizontal Irradiances (DHI) are more widely available than monthly mean Direct Normal Irradiances (DNI) and Global Tilted Irradiances (GTI). Empirical methods have been developed to derive monthly mean DNI and GTI from monthly mean GHI or from GHI and DHI. In this paper, we evaluate two such methods. The first one was the Whitlock Method developed by Charles H. Whitlock (2005) for the NASA POWER project by means of regression of the BSRN data. The method expresses the monthly mean DHI-to-GHI ratio as polynomial functions of monthly mean clearness index, sunset hour angle and noon solar elevation angle on the monthly-average-day. The monthly mean DNI is calculated by dividing the monthly mean GHI-DHI difference, or DirHI, by the cosine of the solar zenith angle at the mid-time between sunrise and solar noon on the monthly-average-day. The second method is the LJCR Method developed by Liu and Jorden (1960) and Collares-Pereira and Rabl (1979), and this method empirically splits monthly mean GHI and DHI into hourly means on the monthly-average-day, and the resulting hourly mean GHI and DHI and their difference, DirHI, can then be used to compute the monthly mean DNI, GTI and the global solar tracker irradiance (GTrI). This method is also used by RETScreen. We recently produced a set of hourly DNI and DHI by bias-correcting the CERES hourly DNI and DHI, and computed hourly GTI and GTrI as well. The data span twenty plus years from March 2000 to near present on a 1 by 1 grid system. The monthly mean CERES GHI and the corrected DHI are used as inputs to the above two methods to compute monthly mean DNI, GTI and GTrI. Through comparisons with the BSRN data, it is found that the Whitlock Method, with slight modification, and the LJCR Method can produce results that are nearly as good as the results derived from the CERES hourly data.

Taiping Zhang↗

Using a laser aureole to invert lidar return

An aureole generated by a laser beam was studied. The strength of the signal redirected towards a sensor high above the surface by a combination of one scattering event in the marine boundary layer (mbl) and one single reflection event from the ocean surface was estimated. A model of mbl aerosol size distributions was used to estimate Mie scattering for a wide range of meteorolocial conditions. The sea surface reflection was determined from a Gaussian model of the wave slopes. These laser aureoles which were estimated over the wide range of conditions and were normalized by the reflected laser light were found to be highly correlated with the optical depth of the boundary layer. By estimating optical depth from the aureole, the Bernoulli-Riccati inversion of lidar return could be constrained and the inversion accuracy improved. A Monte Carlo program was developed to study the laser aureole generated by up to 8 orders of reflection and scattering. The aureole was generated by a narrow, 10 nsec laser pulse at 1.06 microns and measured by a receiver 10 km above the ocean surface. The original theoretical computation compared well with the Monte Carlo method. When multiple scattereffects were included, the normalized aureole was still highly correlated with the mbl optical depth over the range of conditions.

Hooper, W. P.↗

The Determination of the Geometries of Multiple-Element Airfoils Optimized for Maximum Lift Coefficient

Optimum airfoils in the sense of maximum lift coefficient are obtained by a newly developed method. The maximum lift coefficient is achieved by requiring that the turbulent skin friction be zero in the pressure rise region on the upper surface. Under this constraint, the pressure distribution is optimized. The optimum pressure distribution consists of a uniform stagnation pressure on the lower surface, a uniform minimum pressure on the upper surface immediately downstream of the front stagnation point followed by a Stratford zero skin friction pressure rise. When multiple-element airfoils are under consideration, this optimum pressure distribution appears on every element. The parameters used to specify the pressure distribution on each element are the Reynolds number and the normalized trailing edge velocity. The newly developed method of design computes the velocity distribution on a given airfoil and modifies the airfoil contour in a systematic manner until the desired velocity distribution is achieved. There are no limitations on how many elements the airfoil to be designed can have.

Chen, A. W.↗

Black hole normal modes - A semianalytic approach

A new semianalytic technique for determining the complex normal mode frequencies of black holes is presented. The method is based on the WKB approximation. It yields a simple analytic formula that gives the real and imaginary parts of the frequency in terms of the parameters of the black hole and of the field whose perturbation is under study, and in terms of the quantity (n + 1/2), where n = 0, 1, 2,... and labels the fundamental mode, first overtone mode, and so on. In the case of the fundamental gravitational normal modes of the Schwarzschild black hole, the WKB estimates agree with numerical results to better than 7 percent in the real part of the frequency and 0.7 percent in the imaginary part, with the relative agreement improving with increasing angular harmonic. Carried to higher order the method may provide an accurate and systematic means to study black hole normal modes.

Schutz, B. F.↗

Estimating the Properties of Hard X-Ray Solar Flares by Constraining Model Parameters

We wish to better constrain the properties of solar flares by exploring how parameterized models of solar flares interact with uncertainty estimation methods. We compare four different methods of calculating uncertainty estimates in fitting parameterized models to Ramaty High Energy Solar Spectroscopic Imager X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method is also based on the difference between the data and the model, but instead uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the Geostationary Operational Environmental Satellite X1.3 class flare of 2005 January 19, and the other from the X4.8 flare of 2002 July 23.We find that the four methods give approximately the same uncertainty estimates for the 2005 January 19 spectral fit parameters, but lead to very different uncertainty estimates for the 2002 July 23 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that can differ greatly depending on the shape of the hypersurface. The hypersurface arising from the 2005 January 19 analysis is consistent with a normal distribution; therefore, the assumptions behind the three non- Bayesian uncertainty estimation methods are satisfied and similar estimates are found. The 2002 July 23 analysis shows that the hypersurface is not consistent with a normal distribution, indicating that the assumptions behind the three non-Bayesian uncertainty estimation methods are not satisfied, leading to differing estimates of the uncertainty. We find that the shape of the hypersurface is crucial in understanding the output from each uncertainty estimation technique, and that a crucial factor determining the shape of hypersurface is the location of the low-energy cutoff relative to energies where the thermal emission dominates. The Bayesian/MCMC approach also allows us to provide detailed information on probable values of the low-energy cutoff, Ec, a crucial parameter in defining the energy content of the flare-accelerated electrons. We show that for the 2002 July 23 flare data, there is a 95% probability that Ec lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. Further, the low-energy cutoff is more likely to be in the range 25-35 keV than in any other 10 keV wide energy range. The low-energy cutoff for the 2005 January 19 flare is more tightly constrained to 107 +/- 4 keV with 68% probability.

X-rays↗

Development of non-conventional control methods for high angle of attack flight using vortex manipulation

Methods of manipulating the forebody vortices on a generic fighter model to produce controlled yawing moments at high angles of attack were investigated. Forces and moments were measured on the complete configuration and on the isolated forebody. Two schemes for vortex control on the forebody were evaluated: individually-controlled tip strakes and individually-controlled blowing ports. The effectiveness of the forebody strakes in controlling forebody side forces was strongly dependent on the size and location of the strakes. A yawing moment can be produced by deploying the forebody strakes asymmetrically, or can be eliminated by deploying the strakes symmetrically. The most effective strake position was found to be near the primary separation point, between 105 and 120 deg from windward. Blowing on the surface of the model was also shown to have a strong effect on the yawing moment. Blowing either forward or aft tangential to the surface appears to be more effective than blowing normal to the surface. The most effective method to control the yawing moment on the forebody was to minimize the natural asymmetry with a pair of small symmetrically mounted tip strakes and to perturb the vortex system away from the symmetric condition with asymmetric blowing.

Malcolm, Gerald N.↗

Gauge R & R studies for angular measurements

Angular measurements lie on the circumference of a circle and have different characteristics than standard scalar measurements. For applications involving angular data, treating the measured values as scalars can lead to misinterpretation of results if its wrap-around nature is not taken into account. In this article, we propose a variance components wrapped normal model for angular measurements that is analogous to the standard normal model for continuous measurements. This model allows decomposition of contributions to the overall variance to be separated and compared to understand the drivers of the spread of the data. In this work, we analyze gauge R & R study data using Bayesian methods and illustrate the use of this wrapped normal model with simulated and real data. We also performed a small simulation study in considering the design of gauge R & R studies with angular measurements.

42 ENGINEERING↗

Life and death in the microcirculation: a role for angiotensin II

OBJECTIVE: Angiotensin II (ANGII) plays a critical role in the maintenance of the microcirculation and in the anatomical loss of microvessels (rarefaction) that occurs in low renin forms of hypertension and in animals fed a high-salt diet. Elevations in sodium intake can trigger a series of hemodynamic and hormonal responses culminating in a substantial rarefaction of small arterioles and capillaries in both normal and reduced renal mass hypertensive rats. METHODS: Immunohistochemistry, Northern blot, and reverse transcription-polymerase chain reaction (RT-PCR) analysis of microdissected blood vessels were used to localize ANGII receptors in the microcirculation. Chronic infusion of ANGII and other physiologic and pharmacologic manipulations of the reninangiotensin system in rats was combined with morphologic and mathematical analysis of the network architecture. RESULTS: We have shown that rarefaction of the microcirculation can cause an increase in total peripheral resistance, reduced tissue perfusion, decreased oxygen delivery, and impaired organ function. Although the mechanisms by which this occurs are not well understood, a number of key observations point to a role for the renin-angiotensin system in this effect. First, ANGII infused systemically at subpressor levels, or locally into the skeletal muscle interstitium, can induce significant microvessel growth. Second, localization of ANGII receptor proteins by immunohistochemistry and Western blotting and RNA localization by RT-PCR confirm the presence of AT1 receptors, which are growth-stimulatory, and AT2 receptors, which are growth-inhibitory in the microcirculation. Third, maintenance of ANGII at normal levels during periods of hypertension or high-salt diet completely eliminates rarefaction. CONCLUSIONS: Taken together, these results support the hypothesis that ANGII acting through AT1- and AT2-receptor mechanisms modulate vessel density during high-salt diet and hypertension.

Review↗

A Method of Calibrating Airspeed Installations on Airplanes at Transonic and Supersonic Speeds by the Use of Accelerometer and Attitude-Angle Measurements

A method is described for calibrating airspeed installation on airplanes at transonic and supersonic speeds in vertical-plane maneuvers in which use is made of measurements of normal and longitudinal accelerations and attitude angle. In this method all the required instrumentation is carried within the airplane. An analytical study of the effects of various sources of error on the accuracy of an airspeed calibration by the accelerometer method indicated that the required measurements can be made accurately enough to insure a satisfactory calibration.

Zalovick, John A↗

Wave Normal and Poynting Vector Calculations using the Cassini Radio and Plasma Wave Instrument

Wave normal and Poynting vector measurements from the Cassini radio and plasma wave instrument (RPWS) are used to examine the propagation characteristics of various plasma waves during the Earth flyby on August 18, 1999. Using the five-channel waveform receiver (WFR), the wave normal vector is determined using the Means method for a lightning-induced whistler, equatorial chorus, and a series of low-frequency emissions observed while Cassini was in the magnetosheath. The Poynting vector for these emissions is also calculated from the five components measured by the WFR. The propagation characteristics of the lightning-induced whistler were found to be consistent with the whistler wave mode of propagation, with propagation antiparallel to the magnetic field (southward) at Cassini. The sferic associated with this whistler was observed by both Cassini and the Stanford VLF group at the Palmer Station in Antarctica. Analysis of the arrival direction of the sferic at the Palmer Station suggests that the lightning stroke is in the same sector as Cassini. Chorus was observed very close (within a few degrees) to the magnetic equator during the flyby. The chorus was found to propagate primarily away from the magnetic equator and was observed to change direction as Cassini crossed the magnetic equator. This suggests that the source region of the chorus is very near the magnetic equator. The low-frequency emission in the magnetosheath has many of the characteristics of lion roars. The average value of the angle between the wave normal vector and the local magnetic field was found to be 16 degrees, and the emissions ranged in frequency from 0. 19 to 0.75 f(sub ce), where f(sub ce) is the electron cyclotron frequency. The wave normal vectors of these waves were primarily in one direction for each individual burst (either parallel or antiparallel to the local field) but varied in direction throughout the magnetosheath. This suggests that the sources of the emissions are far from the spacecraft and that there are multiple source regions.

Hospodarsky, G. B.↗

Axisymmetric elastodynamic response from normal and radial impact of layered composites with embedded penny-shaped cracks

A method is developed for the dynamic stress analysis of a layered composite containing an embedded penny-shaped crack and subjected to normal and radial impact. Quantitatively, the time-dependent stresses near the crack border can be described by the dynamic stress intensity factors. Their magnitude depends on time, on the material properties of the composite and on the relative size of the crack compared to the composite local geometry. Results obtained show that, for the same material properties and geometry of the composite, the dynamic stress intensity factors for an embedded (penny-shaped) crack reach their peak values within a shorter period of time and with a lower magnitude than the corresponding dynamic stress factors for a through-crack.

Sih, G. C.↗

High-dimensional maximum-entropy phase space tomography

Reconstructing 4D or 6D phase space distributions from 1D or 2D measurements is a challenging inverse problem encountered in particle accelerators. Entropy maximization is an established method to incorporate prior information in the reconstruction, but it is typically infeasible in high-dimensional spaces. In this paper, I review two recent approaches to high-dimensional entropy maximization. The first approach utilizes differentiable simulations and a class of generative models known as normalizing flows, whereas the second approach employs the method of Lagrange multipliers and Markov Chain Monte Carlo (MCMC) sampling. My aim is to provide a short explanation of each method using a common notation. I conclude by mentioning several unsolved problems in phase space tomography.

Hoover, Austin [ORNL] (ORCID:0000000153136962)↗