Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic equations”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A New Method for Deriving Composition of S-type Asteroids from Noisy and Incomplete Near-infrared Spectra

The surface composition of S-type asteroids can be determined using band parameters extracted from their near-infrared (NIR) spectra (0.7–2.50 μm) along with spectral calibrations derived from laboratory samples. In the past, these empirical equations have been obtained by combining NIR spectra of meteorite samples with information about their composition and mineral abundance. For these equations to give accurate results, the characteristics of the laboratory spectra they are derived from should be similar to those of asteroid spectral data (i.e., similar signal-to-noise ratio (S/N) and wavelength range). Here we present new spectral calibrations that can be used to determine the mineral composition of ordinary chondrite-like S-type asteroids. Contrary to previous work, the S/N of the ordinary chondrite spectra used in this study has been decreased to recreate the S/N typically observed among asteroid spectra, allowing us to obtain more realistic results. In addition, the new equations have been derived for five wavelength ranges encompassed between 0.7 and 2.50 μm, making it possible to determine the composition of asteroids with incomplete data. The new spectral calibrations were tested using band parameters measured from the NIR spectrum of asteroid (25143) Itokawa, and comparing the results with laboratory measurements of the returned samples. We found that the spectrally derived olivine and pyroxene chemistry, which are given by the molar contents of fayalite (Fa) and ferrosilite (Fs), are in excellent agreement with the mean values measured from the samples (Fa{sub 28.6±1.1} and Fs{sub 23.1±2.2}), with a maximum difference of 0.6 mol% for Fa and 1.4 mol% for Fs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A High-Order Eulerian–Lagrangian Runge–Kutta Finite Volume (EL–RK–FV) Method for Scalar Nonlinear Conservation Laws

Abstract We present a class of high-order Eulerian–Lagrangian Runge–Kutta finite volume methods that can numerically solve Burgers’ equation with shock formations, which could be extended to general scalar conservation laws. Eulerian–Lagrangian (EL) and semi-Lagrangian (SL) methods have recently seen increased development and have become a staple for allowing large time-stepping sizes. Yet, maintaining relatively large time-stepping sizes post shock formation remains quite challenging. Our proposed scheme integrates the partial differential equation on a space-time region partitioned by linear approximations to the characteristics determined by the Rankine–Hugoniot jump condition. We trace the characteristics forward in time and present a merging procedure for the mesh cells to handle intersecting characteristics due to shocks. Following this partitioning, we write the equation in a time-differential form and evolve with Runge–Kutta methods in a method-of-lines fashion. High-resolution methods such as ENO and WENO-AO schemes are used for spatial reconstruction. Extension to higher dimensions is done via dimensional splitting. Numerical experiments demonstrate our scheme’s high-order accuracy and ability to sharply capture post-shock solutions with large time-stepping sizes.

Chen, Jiajie↗
NASA NTRS1977

Acceleration of protons at 32 Jovian radii in the outer magnetosphere of jupiter

During the inbound pass of Pioneer 10, a rapid ten-fold increase of the 0.2 to MeV proton flux was observed at 32 Jovian radii (R sub J). The total event lasted for 30 minutes and was made up of a number of superimposed individual events. At the time, the spacecraft was in the outer magnetosphere about 7 R sub J below the magnetic equator. Before and after the event, the proton flux was characteristic of the low flux level normally encountered between crossings of the magnetic equator. Flux changes at different energies were coherent within 1 minute; a time comparable to the time resolution of the data. The angular distributions were highly anisotropic with protons streaming towards Jupiter. A field-aligned dumbbell distribution was observed initially, and a pancake distribution just before the flux decayed to its pre-event value. The alpha particle flux changed as rapidly as the proton flux but peaked at different times. The energetic electron flux behaved differently; it increased gradually throughout the period.

Schardt, A. W.↗
NASA NTRS1985

Computational analysis of scramjet dual mode operation

One critical element in the design of a Scramjet is the detailed understanding of the complex flow field in the engine during various phases of operation. One area of interest is the computation of chemically reacting flows in the vicinity of flame holders. The characteristics of a method for solving the Navier-Stokes equations with chemical reactions are proposed. Also of interest are the flame holding characteristics of simple ramp and rearward facing steps. Both of these configurations are considered candidates for Scramjet flame holders.

Source record↗
NASA NTRS1978

Calculation of the flow field in supersonic mixed-compression inlets at angle of attack using the three-dimensional method of characteristics with discrete shock wave fitting

The influence of molecular transport is included in the computation by treating viscous and thermal diffusion terms in the governing partial differential equations as correction terms in the method of characteristics scheme. The development of a production type computer program is reported which is capable of calculating the flow field in a variety of axisymmetric mixed-compression aircraft inlets. The results agreed well with those produced by the two-dimensional method characteristics when axisymmetric flow fields are computed. For three-dimensional flow fields, the results agree well with experimental data except in regions of high viscous interaction and boundary layer removal.

Vadyak, J.↗
NASA NTRS1978

Synthesis and dynamic characteristics of large structures with rotating substructures

The variational equations of motion of large structures with rotating substructures are derived by the substructure synthesis approach, whereby a discretization procedure akin to the Rayleigh-Ritz method is used to represent the elastic motion of every substructure by a suitable set of admissible functions. Using an inclusion principle for gyroscopic systems, the effects on the predicted system dynamic characteristics of truncating the number of admissible functions used for each substructure can be assessed. Criteria for rational selection of admissible functions are presented.

Meirovitch, L.↗
NASA NTRS1948

Equations for the design of two-dimensional supersonic nozzles

Equations are presented for obtaining the wall coordinates of two-dimensional supersonic nozzles. The equations are based on the application of the method of characteristics to irrotational flow of perfect gases in channels. Curves and tables are included for obtaining the parameters required by the equations for the wall coordinates. A brief discussion of characteristics as applied to nozzle design is given to assist in understanding and using the nozzle-design method of this report. A sample design is shown.

Pinkel, I Irving↗
NASA NTRS1974

Sensitivity analysis of torsional vibration characteristics of helicopter rotor blades. Part 1: Structural dynamics analysis

A theoretical investigation of structural vibration characteristics of rotor blades was carried out. Coupled equations of motion for flapwise bending and torsion were formulated for rotor blades with noncollinear elastic and mass axes. The finite element method was applied for a detailed representation of blade structural properties. Coupled structural mass and stiffness coefficients were evaluated. The range of validity of a set of coupled equations of motion linearized with respect to eccentricity between elastic and mass axes was investigated. The sensitivity of blade vibration characteristics to torsion were evaluated by varying blade geometric properties, boundary conditions, and eccentricities between mass and elastic axes.

Bratanow, T.↗
NASA NTRS1993

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗
NASA NTRS1980

Vibration characteristics of a steadily rotating slender ring

Partial differential equations are derived to describe the structural vibrations of a uniform homogeneous ring which is very flexible because the radius is very large compared with the cross sectional dimensions. Elementary beam theory is used and small deflections are assumed in the derivation. Four sets of structural modes are examined: bending and compression modes in the plane of the ring; bending modes perpendicular to the plane of the ring; and twisting modes about the centroid of the ring cross section. Spatial and temporal characteristics of these modes, presented in terms of vibration frequencies and ratios between vibration amplitudes, are demonstrated in several figures. Given a sufficiently high rotational rate, the dynamics of the ring approach those of a vibrating string. In this case, the velocity of traveling wave in the material of the ring approaches in velocity of the material relative to inertial space, resulting in structural modes which are almost stationary in space.

Lallman, F. J.↗
NASA NTRS1974

Sensitivity analysis of torsional vibration characteristics of helicopter rotor blades. Part 2: Aerodynamics and sensitivity analysis

A theoretical investigation of dynamic response characteristics of rotor blades was carried out with special emphasis on torsional degrees-of-freedom. Coupled equations of motion for flapwise bending and torsion were formulated at varying azimuth positions for rotor blades with noncollinear aerodynamic, elastic and mass axes. Both structural and aerodynamic mass, damping and stiffness coefficients were included. The variations of a sample blade at different flight conditions were investigated from these equations. The obtained numerical results were illustrated. The sensitivity of overall blade vibration characteristics to torsional oscillations was also investigated from the equations of motion for the sample blade. The illustrated results show the importance of torsional degrees-of-freedom in rotor blade analysis. Various possibilities of improving the overall response by tuning blade geometric, structural and aerodynamic characteristics are discussed.

Bratanow, T.↗
NASA NTRS1987

The method of projected characteristics for the evolution of magnetic arches

A numerical method of solving fully nonlinear MHD equation is described. In particular, the formulation based on the newly developed method of projected characteristics (Nakagawa, 1981) suitable to study the evolution of magnetic arches due to motions of their foot-points is presented. The final formulation is given in the form of difference equations; therefore, the analysis of numerical stability is also presented. Further, the most important derivation of physically self-consistent, time-dependent boundary conditions (i.e. the evolving boundary equations) is given in detail, and some results obtained with such boundary equations are reported.

Nakagawa, Y.↗
NASA NTRS1950

The solution of the laminar-boundary-layer equation for the flat plate for velocity and temperature fields for variable physical properties and for the diffusion field at high concentration

In connection with Pohlhausen's solution for the temperature field on the flat plate, a series of formulas were indicated by means of which the velocity and temperature field for variable physical characteristics can be computed by an integral equation and an iteration method based on it. With it, the following cases were solved: On the assumption that the viscosity simply varies with the temperature while the other fluid properties remain constant, the velocity and temperature field on the heated and cooled plate, respectively, was computed at the Prandtl numbers 12.5 and 100 (viscous fluids). A closer study of these two cases resulted in general relations: The calculations for a gas of Pr number 0.7 (air) were conducted on the assumption that all fluid properties vary with the temperature, and the velocities are low enough for the heat of friction to be discounted. The result was a thickening of the boundary layers, but no appreciable modification in shearing stress or heat-transfer coefficient.

BOUNDARY LAYER, LAMINAR - PLATES, FLAT↗
DOE OSTI2022

The high-density equation of state in heavy-ion collisions: constraints from proton flow

Abstract A set of different equations of state is implemented in the molecular dynamics part of a non-equilibrium transport simulation (UrQMD) of heavy-ion collisions. It is shown how different flow observables are affected by the density dependence of the equation of state. In particular, the effects of a phase transition at high density are explored, including an expected reduction in mean $$m_T$$ m T . We also show that an increase in $$v_2$$ v 2 is characteristic for a strong softening of the equation of state. The phase transitions with a low coexistence density, $$n_{\text {CE}}<4 n_0$$ n CE < 4 n 0 , show a distinct minimum in the slope of the directed flow as a function of the beam energy, which would be a clear experimental signal. By comparing our results with experimental data, we can exclude any strong phase transition at densities below $$4n_0$$ 4 n 0 .

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1949

Theory of Characteristics

The theory of characteristics will be presented generally for quasilinear differential equations of the second order in two variables. This is necessary because of the manifold requirements to be demanded from the theory of characteristics.

Tollmien, W.↗
NASA NTRS1985

Modeling taper charge with a non-linear equation

Work aimed at modeling the charge voltage and current characteristics of nickel-cadmium cells subject to taper charge is presented. Work reported at previous NASA Battery Workshops has shown that the voltage of cells subject to constant current charge and discharge can be modeled very accurately with the equation: voltage = A + (B/(C-X)) + De to the -Ex where A, B, D, and E are fit parameters and x is amp-hr of charge removed during discharge or returned during charge. In a constant current regime, x is also equivalent to time on charge or discharge.

Mcdermott, P. P.↗
NASA NTRS2001

An Implicit Characteristic Based Method for Electromagnetics

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Beggs, John H.↗
DOE OSTI2022

Recent progress on the APOLLO3 calculation scheme for the FLUOLE2 experiment

Studies of the neutron fluence on reactor pressure vessel and the correlated reactor lifetime issues require very good knowledge of the neutron source distribution in the reactor core and especially in the peripheral region. The FLUOLE2 experiment offers an appropriate starting configuration to investigate the possibility of using the APOLLO3 deterministic code with a two-step calculation scheme. It enables us to calculate a pin-by-pin neutron source distribution in the core, then to carry out the neutron propagation from the core to neutron detectors. In this paper, we focus on a calculation scheme dedicated to the determination of core neutron source distribution using APOLLO3 with the SHEM281 energy mesh, a Method of Characteristics (MOC) to solve the transport equation for the lattice computation step and a 20-groups 3D computation using the MINARET SN solver of APOLLO3. Reactivity and pin-by-pin reaction rate results have been analyzed. Results obtained with both lattice and core computations show a good agreement with those produced through a reference calculation using the TRIPOLI-4 Monte Carlo transport code. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗