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 1,279 records · Page 71

G2Aero (Separable shape tensors for aerodynamic design) [SWR-22-44]

G2Aero is a flexible and practical tool for design and deformation of 2D airfoils and 3D blades using data-driven approaches. G2Aero utilizes the geometry of matrix manifolds—specifically the Grassmannian—to build a novel framework for representing physics-based separable deformations of shapes. G2Aero offers the flexibility to generate perturbations in a customizable way over any portion of the blade. The G2Aero framework utilizes data-driven methods based on a curated database of physically relevant airfoils. Specific tools include: - principal geodesic analysis over normal coordinate neighborhoods of matrix manifolds; - a variety of data-regularized deformations to nominal 2D airfoil shapes; - Riemannian interpolation connecting a sequence of airfoil cross-sections to build 3D blades from 2D data; - consistent perturbations over the span of interpolated 3D blades based on dominant modes from the data-driven analysis.

Doronina, Olga↗

Effects of Chlorine Capture and a Proposed Density Law on the Reactivity of Plutonium Solution Systems

During fissionable material processing, all normal and credible abnormal conditions must remain safely subcritical. Nuclear Criticality Safety (NCS) uses a number of methods to determine subcriticality, one of which is the use of neutron transport codes such as MCNP6. In order to create models for use with MCNP6, both the geometry and materials in fissionable material processes must be known, or assumptions must be made and quantified for the impact to bias. One of the systems with a significant amount of bias due to material modeling assumptions is in the area of aqueous plutonium processing. These solutions are typically plutonium nitrate solutions or plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution and there is no current predictive capability approved for use at Los Alamos National Laboratory (LANL) for modeling them. This research is currently underway to fill the gap and develop an algorithm for use with MCNP6 to model the density of plutonium chloride solutions. The method is to be validated with experimental data for density, and also validated with critical experiments using MCNP6. Note that the Chlorine Worth Study (CWS) was performed in December 2021 to help bridge the gap in chlorine data for critical experiments, and is currently awaiting International Criticality Safety Benchmark Evaluation Project (ICSBEP) review. This study was performed by LANL at the National Criticality Experiments Research Center (NCERC) at the Nevada National Security Site (NNSS). Additional information regarding this experiment may currently be found in LA-UR- 22-29180. Additionally, the Chemistry-Actinide Analytical Chemistry (C-AAC) at LANL has performed a number of solution density measurements for PuCl 3 -HC 1 -H 2 O, allowing for such data be used to create a semi-empirical density via the Pitzer method. The published dataset for the measurements is documented in LA-UR-22-25454. This method has already been tested successfully for aqueous plutonium nitrate solutions in SCALE. Current solution density measurements exist of plutonium concentrations of 0-~142g/L, all at 2M HC1, at temperatures 20-40°C. Additional data was taken for HC1-corrected density values, which essentially mimics the data for a pure PuCl x -water solution. The calculations in this report aim to support the current research by demonstrating the difference in system reactivity for the current modeling method when compared to the new proposed modeling with a density law implementation, which is being written as a Python tool to be used with MCNP6.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Interpreting the Strain Tensor During Injection into a Reservoir

A suite of well tests was conducted to evaluate the feasibility of measuring and interpreting the strain tensor during injection. Water was injected at 1 liter/s for a few days to a few weeks into a well completed at a depth of 530 m in the Bartlesville sandstone in the North Avant Field in Oklahoma. The strain tensor was measured using strainmeters at a depth of 30m and a distance of 220m from the well. We used a Gladwin strainmeter, a state-of-the-art geodetic strainmeter, as well as novel borehole strainmeters developed for the project. After processing, the horizontal strains were tensile and the vertical strain was compressive during injection. The radial strain was less than the circumferential strain. Strain magnitudes increased during injection and recovered afterward. The strain rate was as fast as 100 nanostrain/d at the start of injection, and decreased to 10 nanostrain/d after a few days. The normal strain signals were approximately proportional to the log of time during the first few days of injection, and the semi-log slope was 50 to 75 nanostrain. Data from the Gladwin strainmeter were essentially identical to data from a strainmeter we developed using optical fiber sensors. Tilt was measured using two tiltmeters, but the tilts were erratic compared to the normal strain and were of limited value. Strain data were analyzed using a strain type-curve method, an analytical solution to a poroelastic inhomogeneity, and a stochastic inversion method. Poroelastic simulations explain the data well. Results from the interpretation methods are consistent with each other, and with independent information characterizing the area. The hydraulic diffusivity is approximately 0.5 m2/s and the permeability is between 100 and 500 mD, according to analysis of the strain signals. This is consistent with measurements of permeability from cores, and with analyses of pressure transients. The analyses indicate that Young’s modulus of the reservoir is between 2 and 6 GPa. Strain interpretations indicate the well intersects a permeable lens with a maximum dimension of 500 to 1,000 m. An isopach map of coarse-grained sand confirms a feature roughly 1,000 m in lateral extent. The results demonstrate that the strain signal measured at shallow (30m) depth can be used to infer properties and heterogeneities in an underlying reservoir or aquifer.

Murdoch, Larry↗

Shock-wave reflexion in a relaxing gas

Plane normal shock wave reflection in relaxing gas for shock tube endwall upstream and downstream dynamic pressures, using method of characteristics

Hanson, R. K.↗

The application of the phase space time evolution method to electron shielding

A computer technique for treating the motion of charged and neutral particles and called the phase space time evolution method was developed. This technique employs the computer's bookkeeping capacity to keep track of the time development of a phase space distribution of particles. This method was applied to a study of the penetration of electrons. A 1 MeV beam of electrons normally incident on a semi-infinite slab of aluminum was used. Results of the calculation were compared with Monte Carlo calculations and experimental results. Time-dependent PSTE electron penetration results for the same problem are presented.

Cordaro, M. C.↗

A viscous/potential flow interaction analysis method for multi-element infinite swept wings, volume 1

An analysis method and computer program have been developed for the calculation of the viscosity dependent aerodynamic characteristics of multi-element infinite swept wings in incompressible flow. The wing configuration consisting at the most of a slat, a main element and double slotted flap is represented in the method by a large number of panels. The inviscid pressure distribution about a given configuration in the normal chord direction is determined using a two dimensional potential flow program employing a vortex lattice technique. The boundary layer development over each individual element of the high lift configuration is determined using either integral or finite difference boundary layer techniques. A source distribution is then determined as a function of the calculated boundary layer displacement thickness and pressure distributions. This source distribution is included in the second calculation of the potential flow about the configuration. Once the solution has converged (usually after 2-5 iterations between the potential flow and boundary layer calculations) lift, drag, and pitching moments can be determined as functions of Reynolds number.

Dvorak, F. A.↗

The effect of bending on the stresses in adhesive joints

The problem of stress distribution in adhesive joints where two orthotropic plates are bonded through a flexible adhesive layer is analyzed. It is shown that the effect of bending of the adherends on the stresses in the adhesive layer is very significant. The transverse shear deformations of the adherends appear to have little influence on the adhesive layer stresses. The maximum transverse normal stress in the adhesive is shown to be larger than the maximum longitudinal shear stress. The method of solution is applied to several examples of specific joint geometries and material combinations, and is proven to be applicable to other related problems.

Yuceoglu, U.↗

Reduction of wind-tunnel wall interference by controlled wall flow

Corrections for wind tunnel wall interferences are applied successfully to high lift models when the model to tunnel size ratio is small. The accuracy of the corrections becomes poorer when larger models are tested. An alternate method of testing was developed in which flow through the porous walls of the tunnel was actively controlled so as to approximate free air conditions in the neighborhood of the model during the test. The amount and distribution of the controlled flow through the walls is computed using a potential flow representation of the model based on the measured lift. Theoretical analysis is presented to prove the convergence of the method to free air conditions and to substantiate the general three dimensional theory of operation when the normal flow distribution is continuous. A two dimensional tunnel was constructed to evaluate the concept. Results show that substantial reduction of wall interference may be achieved with relatively low values of porosity of actively controlled walls.

Bernstein, S.↗

Scintillation statistics caused by atmospheric turbulence and speckle in satellite laser ranging

We study the statistics of scintillation at the ground-based receiver for the earth-space-earth retroreflector configuration of satellite laser ranging. These statistics are governed by the joint effects of atmospheric turbulence and speckle produced by the retroreflector array. An expression for the probability density function of scintillation is obtained and evaluated numerically. Comparison of the normalized variance of scintillation calculated by using this function shows good agreement with results obtained by other methods.

Bufton, J. L.↗

Drag reduction effects in turbulent boundary layers over wavy walls

Two dimensional incompressible flow over wavy surfaces are analyzed numerically by spectral methods. Algorithms for periodic flows (Fourier modes in the periodic flow direction and Chebycheff modes in the normal direction), and inflow-outflow boundary conditions (Chebycheff modes used in both directions) are described. Results obtained using both codes are reported for laminar flows. Comparisons with known theoretical and experimental results are made.

Balasubramanian, R.↗

Spacecraft outgas ambient flow interaction

An investigation has been made of the extent to which outgassed or emitted molecules return to a spacecraft as a result of intermolecular collisions with ambient freestream molecules. The governing parameters are discussed and heuristic arguments are used to predict the combination of parameters that is most likely to describe the return flux ratio. This is then tested by an extensive set of numerical calculations using the direct simulation Monte Carlo method. Computations have been made for the flow past spheres and past circular cylinders with their axes normal in the stream. The numerical results lead to empirical expressions that can be used for engineering estimates in spacecraft contamination studies.

Bird, G. A.↗

Testing solar models with global solar oscillations in the 5-minute band

Frequencies of solar oscillation for normal modes described by spherical harmonics with l-values between 0 and 4 are computed. The method of computation is discussed and some of the theoretical uncertainties are examined. It is shown that the standard solar model has eigenfrequencies which do not agree with the frequencies observed for the low l-modes to within the estimated accuracy of either the observed or theoretical frequencies. Four non-standard models are considered: (1) the interior Z abundance is lower than the surface abundance; (2) the interior Z abundance is higher than the surface abundance; (3) the interior Z abundance is altered by mixing; and (4) a large primordial magnetic field remains in the solar core. The effect of all these models on the solar neutrino flux is considered, with the result that the high-Z model is rejected. The conclusions of Bahcall and Ulrich (1971) that a primordial magnetic field increases the neutrino flux are disputed.

Ulrich, R. K.↗

Parametric study of turbine blade platform friction damping using the lumped parameter analysis

The hardware configuration used in the present study of turbine blade planform friction damping, by means of the lumped parameter analysis, is the first turbine stage of the Space Shuttle Main Engine's High Pressure Fuel Turbopump. The analysis procedure solves the nonlinear equations of motion for a turbine blade that is acted on by a platform friction damper, using an iterative matrix method. Attention is given to the effects on blade deflection response of variations in friction coefficient, the normal force on the friction surface interface, blade hysteretic damping, the blade-to-blade phase angle of the harmonic forcing function, and the amplitude of the forcing function.

Dominic, R. J.↗

Calculation of intensity of high energy muon groups observed deep underground

The intensity of narrow muon groups observed in Kolar Gold Field (KGF) at the depth of 3375 m.w.e. was calculated in terms of quark-gluon strings model for high energy hadron - air nuclei interactions by the method of direct modeling of nuclear cascade in the air and muon propagation in the ground for normal primary cosmic ray composition. The calculated intensity has been found to be approx. 10 to the 4 times less than one observed experimentally.

Vavilov, Y. N.↗

Engineering analysis of slender-body aerodynamics using Sychev similarity parameters

The similarity parameters deduced by Sychev for inviscid hypersonic flow over slender bodies are reviewed and used to correlate flowfield, surface-pressure, normal-force and center-of-pressure data for supersonic flow over thin slender wings and smooth slender bodies at low-to-high angles of attack. Although Sychev expected similarity to hold only for hypersonic freestream and cross flows, it is demonstrated empirically that similarity holds for any value of cross-flow Mach number if the axial flow component is supersonic. It is also shown for thin wings that similarity holds for much larger values of aspect ratio than Sychev supposed. A one-term power-law expression is found to fit all of the normal-force correlations which suggests that it may be possible to develop a simple semi-empirical method for estimating the aerodynamic characteristics of arbitrary slender airframes.

Hemsch, Michael J.↗

Solution of the surface Euler equations for accurate three-dimensional boundary-layer analysis of aerodynamic configurations

The three-dimensional boundary-layer equations in the limit as the normal coordinate tends to infinity are called the surface Euler equations. The present paper describes an accurate method for generating edge conditions for three-dimensional boundary-layer codes using these equations. The inviscid pressure distribution is first interpolated to the boundary-layer grid. The surface Euler equations are then solved with this pressure field and a prescribed set of initial and boundary conditions to yield the velocities along the two surface coordinate directions. Results for typical wing and fuselage geometries are presented. The smoothness and accuracy of the edge conditions obtained are found to be superior to the conventional interpolation procedures.

Iyer, V.↗

GRID3D-v2: An updated version of the GRID2D/3D computer program for generating grid systems in complex-shaped three-dimensional spatial domains

In order to generate good quality systems for complicated three-dimensional spatial domains, the grid-generation method used must be able to exert rather precise controls over grid-point distributions. Several techniques are presented that enhance control of grid-point distribution for a class of algebraic grid-generation methods known as the two-, four-, and six-boundary methods. These techniques include variable stretching functions from bilinear interpolation, interpolating functions based on tension splines, and normalized K-factors. The techniques developed in this study were incorporated into a new version of GRID3D called GRID3D-v2. The usefulness of GRID3D-v2 was demonstrated by using it to generate a three-dimensional grid system in the coolent passage of a radial turbine blade with serpentine channels and pin fins.

Steinthorsson, E.↗

Discrete ordinate solution of the radiative transfer equation in the 'polarization normal wave representation'

The transfer equations for normal waves in finite, inhomogeneous and plane-parallel magnetoactive media are solved using the discrete ordinate method. The physical process of absorption, emission, and multiple scattering are accounted for, and the medium may be forced both at the top and bottom boundary by anisotropic radiation as well as by internal anisotropic sources. The computational procedure is numerically stable for arbitrarily large optical depths, and the computer time is independent of optical thickness.

Kylling, A.↗