Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Taylor approximation”

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 109 records · Page 6

Second-order reconstruction of the inflationary potential

To first order in the deviation from scale invariance the inflationary potential and its first two derivatives can be expressed in terms of the spectral indices of the scalar and tensor perturbations, n and n(sub T), and their contributions to the variance of the quadrupole CBR temperature anisotropy, S and T. In addition, there is a 'consistency relation' between these quantities: n(sub T) = (-1/ 7)(T/S). We derive the second-order expressions for the inflationary potential and its first two derivatives and the first-order expression for its third derivative, in terms, of n, n(sub T), S, T, and dn/d ln gamma. We also obtain the second-order consistency relation, n(sub T) = (-1/7)(T/S)(1 + 0.11(T/S) + 0.15(n-1)). As an example we consider the exponential potential, the only known case where exact analytic solutions for the perturbation spectra exist. We reconstruct the potential via Taylor expansion (with coefficients calculated at both first and second order), and introduce the Pade approximate as a greatly improved alternative.

Liddle, Andrew R.

The Mantle and Basalt-Crust Interaction Below the Mount Taylor Volcanic Field, New Mexico

The Mount Taylor Volcanic Field (MTVF) lies on the Jemez Lineament on the southeastern margin of the Colorado Plateau. The field is centered on the Mt. Taylor composite volcano and includes Mesa Chivato to the NE and Grants Ridge to the WSW. MTVF magmatism spans approximately 3.8-1.5 Ma (K-Ar). Magmas are dominantly alkaline with mafic compositions ranging from basanite to hy-basalt and felsic compositions ranging from ne-trachyte to rhyolite. We are investigating the state of the mantle and the spatial and temporal variation in basalt-crustal interaction below the MTVF by examining mantle xenoliths and basalts in the context of new mapping and future Ar-Ar dating. The earliest dated magmatism in the field is a basanite flow south of Mt. Taylor. Mantle xenolith-bearing alkali basalts and basanites occur on Mesa Chivato and in the region of Mt. Taylor, though most basalts are peripheral to the main cone. Xenolith-bearing magmatism persists at least into the early stages of conebuilding. Preliminary examination of the mantle xenolith suite suggests it is dominantly lherzolitic but contains likely examples of both melt-depleted (harzburgitic) and melt-enriched (clinopyroxenitic) mantle. There are aphyric and crystal-poor hawaiites, some of which are hy-normative, on and near Mt. Taylor, but many of the more evolved MTVF basalts show evidence of complex histories. Mt. Taylor basalts higher in the cone-building sequence contain >40% zoned plagioclase pheno- and megacrysts. Other basalts peripheral to Mt. Taylor and at Grants Ridge contain clinopyroxene and plagioclase megacrysts and cumulate-textured xenoliths, suggesting they interacted with lower crustal cumulates. Among the questions we are addressing: What was the chemical and thermal state of the mantle recorded by the basaltic suites and xenoliths and how did it change with time? Are multiple parental basalts (Si-saturated vs. undersaturated) represented and, if so, what changes in the mantle or in the tectonic regime allowed their coexistence or caused the transition?

Schrader, Christian M.

Jupiter's great red spot revisited

On the original Taylor column theory of Jupiter's Great Red Spot, the fixed latitude of the Spot is taken to imply that the Taylor column in Jupiter's atmosphere is associated with a disturbance such as a topographic feature of the surface Q underlying the atmosphere. The alternative suggestion that the Taylor column is produced by a solid raft floating at depth in the atmosphere is somewhat easier to reconcile with the approximately 10s difference between the respective rotation periods P sub S and P sub R of the Red Spot and of the radio sources, but it does not account so readily for the fixed latitude of the Spot unless it can be shown that the raft is in stable equilibrium under the north-south components of the dynamical forces, including wind effects, acting upon it. A slight wavering of the upper end of the Taylor column relative to the lower end could account at least in part for the most rapid variations in P sub S, but the slow large-amplitude variations in P sub S must reflect changes in the longitudinal motion of either the surface Q or of the raft. By generalizing the Proudman-Taylor theorem to the case of a non-homogeneous fluid it is shown that the Taylor column theory does not imply very special and therefore unlikely horizontal and vertical temperature variations in Jupiter's atmosphere, thus refuting a widely-held belief to the contrary.

Hide, R.

A semiempirical method for analysis of the reflectance spectra of binary mineral mixtures

A simple semiempirical method is presented for determination of the spectral reflectance of a powdered binary mineral mixture. This technique uses a two-stream radiative transfer model (a modified Kubelka-Munk model) on a particulate medium of isotroic scatterers The particles are assumed to be much larger than the wavelength of light under consideration. This same method can be used to determine the relative proportion of components in a mixture for which the spectral reflectance is known. Binary mixtures of olivine, two pyroxenes, and magnetite are used to test this model. The theoretical and empirical results agree approximately within experimental errors.

Johnson, P. E.

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics

An implicit and stiffly stable finite element CFD algorithm for unsteady aerodynamics

A stable and accurate finite element CFD algorithm for hyperbolic/incompletely parabolic conservation law systems is described and verified. It combines a Taylor weak statement FEM, an optimal implicit Runge-Kutta time integration algorithm, and a matrix tensor product approximate factorization linear algebra procedure. The results of computational experiments show that the developed algorithm is robust.

Baker, A. J.

On a finite element CFD algorithm for compressible, viscous and turbulent aerodynamic flows

This paper develops and analyses individual construction aspects of an efficient and accurate finite element algorithm for prediction of viscous and turbulent flow fields of impact in aerodynamics. The theoretical construction employs a Taylor weak statement (TWS) for coincident embedding of stability mechanisms within a classic Galerkin finite element formulation of semidiscrete approximation error orthogonalization. A wide variety of the stabilizing mechanisms of independently derived CFD algorithms are contained within the TWS theory. An implicit construction that meets the requirement of efficient convergence to steady state is developed. The theoretical asymptotic error estimates of the TWS finite element algorithm for supersonic and viscous boundary layer flows are verified. Application to a three-dimensional turbulent flow is cited.

Baker, A. J.

Shape optimization of 3D continuum structures via force approximation techniques

The existing need to develop methods whereby the shape design efficiency can be improved through the use of high quality approximation methods is addressed. An efficient approximation method for stress constraints in 3D shape design problems is proposed based on expanding the nodal forces in Taylor series with respect to shape variations. The significance of this new method is shown through elementary beam theory calculations and via numerical computations using 3D solid finite elements. Numerical examples including the classical cantilever beam structure and realistic automotive parts like the engine connecting rod are designed for optimum shape using the proposed method. The numerical results obtained from these methods are compared with other published results, to assess the efficiency and the convergence rate of the proposed method.

Vanderplaats, Garret N.

Spacelab experiments on convection in a rotating spherical shell with radial gravity

Data from the geophysical fluid flow cell experiment on Spacelab 3 revealed the diverse forms of fluid motion that can occur in geometrically simple models of global convection with rotation and radial gravity. The following types of convection were observed: midlatitude waves interacting with the low-latitude columnar convection rolls or 'banana cells', spiral waves' near the poles when latitudinal heating gradients are present on the bounding surfaces, and 'triangular waves' coupling midlatitude and equatorial disturbances under similar differential heating. It is believed that concepts based on the Taylor-Proudman theorem have a central role in many flows realized in the rapidly rotating experiments when H is approximately 0; however, the resulting banana cells are subject to secondary instabilities that produce intricate time dependence and eventually turbulent flows in which only hints of the simpler patterns are recognizable.

Toomre, Juri

Armalcolite/ilmenite - Mineral chemistry, paragenesis, and origin of textures

In order to understand the origin of differing ilmenite/armalcolite textural relationships, bulk rock composition, oxygen fugacity and cooling rate have been experimentally investigated. Results show that armalcolite and ilmenite compositions are dependent upon the liquid composition and temperature at the time of crystallization. Furthermore, the different armalcolite/ilmenite textures rely to a greater extent on f sub O2 than on cooling rate or bulk rock composition, and the Ti(3+) content of armalcolite can be used as an indicator of the approximate f sub O2 prevailing at the time of crystallization.

Stanin, F. T.

Venus cloud structure and water vapor abundance from Mariner 10 observations

Observations of the Venus atmosphere with the infrared radiometer on Mariner 10 have been analyzed by Taylor (1975) in terms of the vertical distribution of opacity at wavelengths near 11 microns and 45 microns in the thermal infrared. In this paper, we discuss models of the Venus atmosphere which are consistent with the inferred opacity structure. Either a two-layer cloud structure, or a single cloud deck overlaid by a layer containing approximately 40 precipitable microns of water vapor, would have the required limb-darkening characteristics at the wavelengths of observation.

Taylor, F. W.

A globally well-posed finite element algorithm for aerodynamics applications

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems.

Iannelli, G. S.

The stability of cosmic-ray-dominated shocks - A secondary instability

Previous studies in one spatial dimension showed that the precursors of cosmic-ray-dominated shocks are unstable against traveling acoustic disturbances. Here we report that a secondary, Rayleigh-Taylor type instability can exist in cosmic-ray-dominated media influenced by the acoustic instability. Using the local WKB approximation, the growth rate of the secondary instability is shown to be comparable to that of the one-dimensional acoustic instability itself in the cases we have considered. The nonlinear development has been followed numerically with a two-dimensional PPM hydrodynamics code that also incorporates the two-fluid cosmic-ray energy equation. We show that the secondary instability may cause the precursor and postshock flows to become highly turbulent.

Ryu, Dongsu

High-resolution mid-infrared spectra of Co II, Ni I, and Fe II in SN 1987A

Ground-based infrared observations of SN 1987A on day 612 after the explosion have yielded resolved line profiles of Co II, Ni I, Fe II at 10.52, 11.31, and 17.94 micron, respectively. The spectra were taken at a resolving power of about 1000 with an array grating spectrometer on the 4 m telescope of Cerro Tololo Inter-American Observatory. Based on the observed line intensities we have estimated the minimum mass of each ion: M(Co II) = (6.0 +/- 1.8) x 10 exp -5 solar mass; M(Ni I) = (1.1 +/- 0.1) x 10 exp -3 solar mass; and M(Fe II) = (8.0 +/- 1.5) x 10 exp -3 solar mass. From these we infer total masses for cobalt, nickel, and iron in the ejecta. The nickel and iron line profiles are markedly asymmetric. We interpret these as arising from two components, one centered on the stellar rest velocity with an approximately 3250 km/s full width, and the second at about +1200 km/s with an approximately 1100 km/s full width. The asymmetry may represent a large-scale fracturing of the ejecta by Rayleigh-Taylor instabilities.

Jennings, D. E.

Effects of Strain and Strain Rate on Dynamic Grain Growth and Subgrain Evolution During Plastic Deformation of an Interstitial-Free Steel at 850 ° C

Here, the effects of strain and strain rate on dynamic grain growth (DGG) and subgrain evolution are reported for an interstitial-free steel deformed at 850 ° C. Microstructures produced during tension tests at true-strain rates of 10 -4 and to 10 -3 s -1 true strains ranging from 0.02 to 0.2 were preserved following deformation. These were characterized using electron backscatter diffraction (EBSD), including the application of spherical harmonic transform indexing to produce high-angular-resolution EBSD (HR-EBSD) data. HR-EBSD data resolved the small misorientation angles of subgrain boundaries while imaging much larger data fields than possible with previously available techniques. The resulting data confirmed that steady-state flow stress is inversely proportional to the average subgrain size and that subgrain boundary misorientation angle increases with strain. The following new observations are reported. The rate of DGG increased with respect to time but decreased with respect to strain as strain rate increased. This behavior is rationalized through a simple model using separate rate parameters for the effects of time and strain. Subgrain size was not constant during steady-state deformation, but decreased slowly with increasing strain. Subgrain size distributions and subgrain boundary misorientation angle distributions were measured, and both remained approximately log-normal during steady-state deformation. Subgrain evolution demonstrated no dependence on parent grain size, crystallographic orientation, or Taylor factor. These new data suggest that steady-state flow stress is more likely controlled by the dislocation density internal to subgrains than by the spacing between subgrain boundaries.

dynamic grain growth

A new propagation method for the radial Schroedinger equation

A new method for propagating the solution of the radial Schroedinger equation is derived from a Taylor series expansion of the wavefunction and partial re-summation of the infinite series. Truncation of the series yields an approximation to the exact propagator which is applied to a model calculation and found to be highly convergent.

Devries, P. L.

Adaptive Power Flow Approximations With Second-Order Sensitivity Insights

The power flow equations are fundamental to power system planning, analysis, and control. However, the inherent non-linearity and non-convexity of these equations present formidable obstacles in problem-solving processes. To mitigate these challenges, recent research has proposed adaptive power flow linearizations that aim to achieve accuracy over wide operating ranges. The accuracy of these approximations inherently depends on the curvature of the power flow equations within these ranges, which necessitates considering second-order sensitivities. In this paper, we leverage second-order sensitivities to both analyze and improve power flow approximations. We evaluate the curvature across broad operational ranges and subsequently utilize this information to inform the computation of various sample-based power flow approximation techniques. Additionally, we leverage second-order sensitivities to guide the development of rational approximations that yield linear constraints in optimization problems. In conclusion, this approach is extended to enhance accuracy beyond the limitations of linear functions across varied operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION