SEARCH · Engineering Papers

Results for “numerical calculations: variational”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Configuration and impurity quantification of AmLi sources using radiography and gamma spectroscopy

This paper presents the non-destructive assay of two AmLi sources used for detector characterization and active interrogation at PNNL, identified as MRC-101 and MRC-103. First, a detailed description of the internal configuration of 2724-BT encapsulated MRC sources is established. X-ray radiographs of the AmLi sources are reported for the first time to verify the dimensions and orientation of the encapsulation. Notably, a density variation was observed as a contrast change across the height of the source cavities. Contrary to commonly made assumptions in literature, it is indirectly observed that the AmLi source/target material only occupies a portion of the source cavity. The intensity of the prompt gamma-rays resulting from the inelastic scatter of alpha particles on 7Li was used to calculate the mass ratios of LiH:AmO2. The gamma-ray spectra revealed the presence of numerous impurities in both source and target material (243Am, 237Np, 154Eu, 23Na, and 9Be). Neutron-gamma production ratios were applied to the intensity of the Doppler-broadened peaks from alpha-induced reactions with 7Li, 9Be, and 23Na to calculate neutron rates for both sources. A neutron/alpha production yield (derived in this work) was then related to the atomic concentration of these isotopes to find the approximate masses within each source. Ultimately, the AmLi source density was estimated with prompt gamma analysis based on partial source volumes as 0.45±0.08 g/cm^3.

DTRA, AmLi, neutron↗

NASA NTRS1987

Vortex simulation of three-dimensional mixing layers

Spanwise structures of 'nominally' two-dimensional (2D), spatially growing, turbulent mixing layers are simulated numerically by a three-dimensional (3D) vortex method. Small-amplitude 3D disturbances are introduced into an otherwise 2D flow field. Results show that a large-scale spanwise variation of the flow field is produced because of amplification of initial disturbances. Pairs of counterrotating streamwise vortices are formed as a result of stretching of primary spanwise vortices. The streamwise vortices are formed at a fixed spanwise location which depends on initial disturbances. The calculated magnitude of streamwise vorticity is close to that of the spanwise vorticity. The results also suggest that the presence of background disturbances may be essential for the formation of streamwise vortices. The results are in good qualitative agreement with experiments.

Inoue, Osamu↗

NASA NTRS2004

Shock-Wave/Boundary-Layer Interactions in Hypersonic Low Density Flows

Results of numerical simulations of Mach 10 air flow over a hollow cylinder-flare and a double-cone are presented where viscous effects are significant. The flow phenomena include shock-shock and shock- boundary-layer interactions with accompanying flow separation, recirculation, and reattachment. The purpose of this study is to promote an understanding of the fundamental gas dynamics resulting from such complex interactions and to clarify the requirements for meaningful simulations of such flows when using the direct simulation Monte Carlo (DSMC) method. Particular emphasis is placed on the sensitivity of computed results to grid resolution. Comparisons of the DSMC results for the hollow cylinder-flare (30 deg.) configuration are made with the results of experimental measurements conducted in the ONERA RSCh wind tunnel for heating, pressure, and the extent of separation. Agreement between computations and measurements for various quantities is good except that for pressure. For the same flow conditions, the double- cone geometry (25 deg.- 65 deg.) produces much stronger interactions, and these interactions are investigated numerically using both DSMC and Navier-Stokes codes. For the double-cone computations, a two orders of magnitude variation in free-stream density (with Reynolds numbers from 247 to 24,7 19) is investigated using both computational methods. For this range of flow conditions, the computational results are in qualitative agreement for the extent of separation with the DSMC method always predicting a smaller separation region. Results from the Navier-Stokes calculations suggest that the flow for the highest density double-cone case may be unsteady; however, the DSMC solution does not show evidence of unsteadiness.

Moss, James N.↗

DOE OSTI2021

Multiscale Electrothermal Design of a Modular Multilevel Converter for Grid-Tied Applications

As a key feature of modular multilevel converters (MMCs), a large number of semiconductor devices are employed in the converter and distributed over a stack of submodules. Each submodule has a strict temperature limit that imposes constraints on the operating range of the converter. Different loading conditions/losses in the submodules lead to unavoidable temperature variations inside the MMC, which consequently affect the system-level performance and reliability. This paper is focused on the electrothermal analysis and design of a medium-voltage silicon carbide (SiC)-based MMC system, from submodule power semiconductors to the overall MMC system integration of multiple submodules. Loss calculations are performed to estimate the cooling requirements and aid thermal design at different levels. The performance of forced-air cooling approach is analyzed within a numerical modeling framework. Maximum temperature of the SiC power modules is predicted using numerical tools. Thermal design of the MMC cabinet with various arrangement of air inlet(s) and outlet(s) is investigated and compared from a cooling performance perspective.

CFD↗

NASA NTRS1973

Effects of inclination and eccentricity on optimal trajectories between earth and Venus

The true optimal transfers, including the effects of the inclination and eccentricity of the planets' orbits, between earth and Venus are presented as functions of the corresponding idealized Hohmann transfers. The method of determining the optimal transfers using the calculus of variations is presented. For every possible Hohmann window, specified as a continuous function of the longitude of perihelion of the Hohmann trajectory, the corresponding numerically exact optimal two-impulse transfers are given in graphical form. The cases for which the optimal two-impulse transfer is the absolute optimal, and those for which a three-impulse transfer provides the absolute optimal transfer are indicated. This information furnishes everything necessary for quick and accurate orbit calculations for preliminary Venus mission analysis. This makes it possible to use the actual optimal transfers for advanced planning in place of the standard Hohmann transfers.

Gravier, J.-P.↗

NASA NTRS1982

Multiple scattering of electromagnetic waves by rain

As the operating frequencies of communications systems move higher into the millimeter wave region, the effects of multiple scattering in precipitation media become more significant. In this paper, general formulations are presented for single, first-order multiple, and complete multiple scattering. Included specifically are distributions of particle size, shape, and orientation angle, as well as variation in the medium density along the direction of wave propagation. Calculations are performed for rain. It is shown that the effects of higher-order scattering are not noticeable in either attenuation or channel isolation on a dual-polarized system until frequencies of about 30 GHz are reached. The complete multiple-scattering formulation presented gives accurate results at high millimeter wave frequencies as well as including realistic medium parameter distributions. Furthermore, it is numerically efficient.

Tsolakis, A.↗

NASA NTRS1990

Stress Dependence of Ultrasonic Velocity in Unidirectional Graphite/Epoxy Composites for Longitudinal Waves Propagating Along the Direction of Stress

The first measurements of the stress induced velocity changes for propagation directions along the direction of applied stress in gr/ep composites have been presented. For propagation and stress direction perpendicular to the fiber direction, the data demonstrated a linear relation between normalized velocity shift and stress. After corrections for the delay line were made, the slope or SAC was determined and compared favorably with the expected value calculated from the previously determined nonlinear coefficients of this material. The ratio of the SAC to the elastic compliance for this direction of loading was evaluated and found to have a value similar to numerous other materials which have very different linear elastic properties. Measurements with stress and propagation along the fibers yielded unusual behavior. The curves were very nonlinear and even shifted direction at higher loads. The large scatter in the data due to bond variations made separation of material effects from bond induced artifacts impossible. Thus the SAC, R, and the remaining two unknown TOEC's could not be determined for this direction of propagation. These measurements further expand the basis of determining nonlinear elastic properties of composite materials. These properties may be useful in developing much needed NDE techniques to determine such important parameters as residual stress after cure and residual strength after impact damage. Additional study is needed to measure the nonlinear behavior in other composite materials including angle ply laminates. Also, other techniques to measure elastic nonlinearity such as harmonic generation should be applied to composites to improve the understanding of these properties and their importance.

Prosser, William H.↗

NASA NTRS1979

On the Angular Variation of Solar Reflectance of Snow

Spectral and integrated solar reflectance of nonhomogeneous snowpacks were derived assuming surface reflection of direct radiation and subsurface multiple scattering. For surface reflection, a bidirectional reflectance distribution function derived for an isotropic Gaussian faceted surface was considered and for subsurface multiple scattering, an approximate solution of the radiative transfer equation was studied. Solar radiation incident on the snowpack was decomposed into direct and atmospherically scattered radiation. Spectral attenuation coefficients of ozone, carbon dioxide, water vapor, aerosol and molecular scattering were included in the calculation of incident solar radiation. Illustrative numerical results were given for a case of North American winter atmospheric conditions. The calculated dependence of spectrally integrated directional reflectance (or albedo) on solar elevation was in qualitative agreement with available observations.

Chang, A. T. C.↗

NASA NTRS1981

On the angular variation of solar reflectance of snow

Spectral and integrated solar reflectance of nonhomogeneous snowpacks are derived assuming surface reflection of direct radiation and subsurface multiple scattering. For surface reflection, a bidirectional reflectance distribution function derived for an isotropic Gaussian faceted surface is considered, and for subsurface multiple scattering an approximate solution of the radiative transfer equation is studied. Solar radiation incident on the snowpack is decomposed into direct and atmospherically scattered radiation. Spectral attenuation coefficients of ozone, carbon dioxide, water vapor, aerosol, and molecular scattering are included in the calculation of incident solar radiation. Illustrative numerical results are given for a case of North American winter atmospheric condition. The calculated dependence of spectrally integrated directional reflectance (or albedo) on solar elevation is in qualitative agreement with available observations.

Choudhury, B. J.↗

NASA NTRS1972

Optimization of simple structures with higher mode frequency constraints.

Results of a study of least weight optimization of simple structures with a single natural frequency constraint. Variational techniques were used to derive the necessary equations, and numerical methods were used, where required, to find solutions to the resulting nonlinear, two-point boundary value problems. The numerical results suggest a similarity between the fixed fundamental frequency solution and solutions in which a single frequency other than fundamental is held fixed. Such a similarity is shown to occur under special circumstances, and the knowledge of the fundamental solution in these cases makes it possible to calculate the solution to problems with constraints on frequencies other than the fundamental.

Weisshaar, T. A.↗

NASA NTRS1973

A numerical procedure for the parametric optimization of three dimensional scramjet nozzles

A numerical procedure permitting the rapid determination of the internal performance of a class of scramjet nozzle configurations is presented. The approach developed is based on the construction of quasi two dimensional simple wave networks, where lateral expansion effects are incorporated by one dimensional approximations. A numerical procedure following this approach has has been developed and results obtained are highly comparable to those obtained employing a characteristic procedure. The numerical program developed permits the parametric variation of cowl length, turning angles on the cowl and vehicle undersurface and lateral expansion and is subject to fixed constraints such as the vehicle length and nozzle exit height. The program requires uniform initial conditions at the burner exit station and yields the location of all predominant wave zones, accounting for lateral expansion effects. In addition, the program yields the detailed pressure distribution on the cowl and vehicle undersurface and calculates the nozzle thrust, lift and pitching moment.

Dash, S.↗

NASA NTRS1978

On the supposed anticorrelation of solar polar and equatorial rotation rates

The anticorrelation between two of the three parameters used to calculate the angular velocity of the sun's differential rotation is here said to be due to numerical coupling. A computer simulation technique shows that the relationship between the two parameters is caused by the effect of noise on the least-squares analysis used to obtain the three parameters used to determine the angular velocity in terms of the heliographic latitude. The computer simulation technique is described. The supposed anticorrelation had been used to infer that variations of the sun's polar and equatorial rotation rates are anticorrelated.

Duvall, T. L., Jr.↗

NASA NTRS2000

Asymptotic Solutions for Optical Properties of Large Particles with Strong Absorption

For scattering calculations involving nonspherical particles such as ice crystals, we show that the transverse wave condition is not applicable to the refracted electromagnetic wave in the context of geometric optics when absorption is involved. Either the TM wave condition (i.e., where the magnetic field of the refracted wave is transverse with respect to the wave direction) or the TE wave condition (i.e., where the electric field is transverse with respect to the propagating direction of the wave) may be assumed for the refracted wave in an absorbing medium to locally satisfy the electromagnetic boundary condition in the ray tracing calculation. The wave mode assumed for the refracted wave affects both the reflection and refraction coefficients. As a result, a nonunique solution for these coefficients is derived from the electromagnetic boundary condition. In this study we have identified the appropriate solution for the Fresnel reflection/refraction coefficients in light scattering calculation based on the ray tracing technique. We present the 3 x 2 refraction or transmission matrix that completely accounts for the inhomogeneity of the refracted wave in an absorbing medium. Using the Fresnel coefficients for an absorbing medium, we derive an asymptotic solution in an analytical format for the scattering properties of a general polyhedral particle. Numerical results are presented for hexagonal plates and columns with both preferred and random orientations. The asymptotic theory can produce reasonable accuracy in the phase function calculations in the infrared window region (wavelengths near 10 micron) if the particle size (in diameter) is on the order of 40 micron or larger. However, since strong absorption is assumed in the computation of the single-scattering albedo in the asymptotic theory, the single scattering albedo does not change with variation of the particle size. As a result, the asymptotic theory can lead to substantial errors in the computation of single-scattering albedo for small and moderate particle sizes. However, from comparison of the asymptotic results with the FDTD solution, it is expected that a convergence between the FDTD results and the asymptotic theory results can be reached when the particle size approaches 200 micron. We show that the phase function at side-scattering and backscattering angles is insensitive to particle shape if the random orientation condition is assumed. However, if preferred orientations are assumed for particles, the phase function has a strong dependence on scattering azimuthal angle. The single-scattering albedo also shows very strong dependence on the inclination angle of incident radiation with respect to the rotating axis for the preferred particle orientations.

Yang, Ping↗

NASA NTRS1992

Numerical study of low pressure nuclear thermal rockets

The flowfields and performance of low pressure nuclear thermal rockets, which use hydrogen as a propellant, are studied by solving the Navier-Stokes equations and the species equations. A finite-rate chemistry model is used in the species equations, and the turbulence is simulated by the Baldwin-Lomax turbulence model with a modified van Driest's damping constant. The calculated results for the chamber temperatures of 3200 K and 4000 K with a chamber pressure range of 0.1 atm to 6 atm are presented as contours, centerline variations, and exit profiles. The performance values from the present calculations, such as the vacuum specific impulse and thrust, are compared with those from the 1D, inviscid equilibrium and frozen flow code.

Kim, Suk C.↗

NASA NTRS1977

A catalog of 0.2 A resolution far-ultraviolet stellar spectra measured with Copernicus

Spectra between 1000 and 1450 A for 60 O- and B-type stars observed by Copernicus at 0.2-A resolution are presented in three forms: tables containing the numerical data, plots showing renormalized spectra, and synthetic photographic spectra. The data have been corrected for all instrument effects of importance for the photometric accuracy except fluctuations in continuum level caused by small variations in spacecraft guidance. Spectrometer sensitivity curves are provided for use in converting to absolute fluxes. It is expected that this catalog will be of use for research on many aspects of stellar UV spectra, including spectral classification, line identification, abundance determinations, spectrum synthesis, model atmosphere calculations, flux distributions, bolometric corrections, stellar winds, and mass loss.

Snow, T. P., Jr.↗

NASA NTRS1981

A sensitivity model for energy consumption in buildings. Part 1: Effect of exterior environment

A simple analytical model is developed for the simulation of seasonal heating and cooling loads of any class of buildings to complement available computerized techniques which make hourly, daily, and monthly calculations. An expression for the annual energy utilization index, which is a common measure of rating buildings having the same functional utilization, is derived to include about 30 parameters for both building interior and exterior environments. The sensitivity of a general class building to either controlled or uncontrolled weather parameters is examined. A hypothetical office type building, located at the Goldstone Space Communication Complex, Goldstone, California, is selected as an example for the numerical sensitivity evaluations. Several expressions of variations in local outside air temperature, pressure, solar radiation, and wind velocity are presented.

Lansing, F. L.↗

NASA NTRS1981

Observations and models of inertial waves in the deep ocean

A study of the structure of the inertial peak in deep ocean kinetic energy is presented, based on records taken from Polymode arrays deployed in the western North Atlantic Ocean. Results are interpreted in terms of both local sources and turning point effects on internal waves generated at lower latitudes, and it is found that three classes of environment and their corresponding spectra emerge from peak height variations: (1) the 1500-m level near the Mid-Atlantic Ridge, with the greatest peak height of 18 dB; (2) the upper and deep ocean over rough topography and the deep ocean underneath the Gulf Stream, with the intermediate peak height of 11.5 dB; and (3) the deep ocean over smooth topography, with the lowest peak height of 7.5 dB. Using the globally valid wave functions obtained by Munk and Phillips (1968), frequency spectra near f are calculated numerically. The model is latitudinally dependent, with the frequency shift and bandwidth of the inertial peak decreasing with latitude.

Fu, L.-L.↗

NASA NTRS1994

An Interactive Educational Tool for Compressible Aerodynamics

A workstation-based interactive educational tool was developed to aid in the teaching of undergraduate compressible aerodynamics. The tool solves for the supersonic flow past a wedge using the equations found in NACA 1135. The student varies the geometry or flow conditions through a graphical user interface and the new conditions are calculated immediately. Various graphical formats present the variation of flow results to the student. One such format leads the student to the generation of some of the graphs found in NACA-1135. The tool includes interactive questions and answers to aid in both the use of the tool and to develop an understanding of some of the complexities of compressible aerodynamics. A series of help screens make the simulator easy to learn and use. This paper will detail the numerical methods used in the tool and describe how it can be used and modified.

Benson, Thomas J.↗