Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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,045 records · Page 58

The effects of magnetic fields on the growth of thermal instabilities in cooling flows

The effects of heat conduction and magnetic fields on the growth of thermal instabilities in cooling flows are examined using a time-dependent hydrodynamics code. It is found that, for magnetic field strengths of roughly 1 micro-Gauss, magnetic pressure forces can completely suppress shocks from forming in thermally unstable entropy perturbations with initial length scales as large as 20 kpc, even for initial amplitudes as great as 60 percent. Perturbations with initial amplitudes of 50 percent and initial magnetic field strengths of 1 micro-Gauss cool to 10,000 K on a time scale which is only 22 percent of the initial instantaneous cooling time. Nonlinear perturbations can thus condense out of cooling flows on a time scale substantially less than the time required for linear perturbations and produce significant mass deposition of cold gas while the accreting intracluster gas is still at large radii.

David, Laurence P.↗

Conceptual design for a Mercury relativity satellite

It was shown earlier that 1 x 10 to the -14th Doppler data and 3 cm accuracy range measurements to a small Mercury Relativity Satellite in a polar orbit with four-hour period can give high-accuracy tests of gravitational theory. A particular conceptual design has been developed for such a satellite, which would take less than 10 percent of the approach mass for a possible future Mercury Orbiter Mission. The spacecraft is similar to the Pioneer Venus Orbiter, but scaled down by about a factor four in linear dimensions. A despun antenna 30 cm in diameter is used for tracking. The transmitted power is roughly 0.2 watts at K-band and 0.5 watts at X-band. The orbit parameters for individual eight-hour arcs and the gravity field of Mercury through degree and order 10 are determined mainly from the Doppler data. A 50 MHz K-band sidetone system provides the basic ranging accuracy. The spacecraft mass is 50 kg or less.

Bender, P. L.↗

Rugged Direct-Current Transducer

Direct-current transducer withstands thermal and mechanical shocks. Operates from single, nominally 15-V power supply and nearly insensitive to both large variations in temperature and variations of as much as plus or minus 5 V in supply voltage. Its output voltage highly linear function of sensed current, with full-scale value of about 3 Vdc and offset of about 0.1 Vdc at 0 sensed current. Ruggedness of transducer due in large part to novel magnetic core, machined from solid block of nickel/iron high-permeability, low-hysteresis alloy.

Mclyman, W. T.↗

An experimental analysis of critical factors involved in the breakdown process of leading edge vortex flows

Experimental crosswire measurements of the flowfield above a 70 and 75 degree flat plate delta wing were performed at a Reynolds number of 250,000. Survey grids were taken normal to the platform at a series of chordwise locations for angles of attack of 20 and 30 degrees. Axial and azimuthal vorticity distributions were derived from the velocity fields. The dependence of circulation on distance from the vortex core as well as on chordwise location was examined. The effects of nondimensionalization in comparison with other experimental data was made. The circulation distribution scales with the local semispan and grows approximately linearly in the chordwise direction. For regions of the flow outside of the vortex subcore, the circulation at any chordwise station was observed to vary logarithmically with distance from the vortex axis. The circulation was also found to increase linearly with angle of incidence at a given chordwise station. A reduction in the local circulation about the vortex axis occurred at breakdown. The spanwise distribution of axial vorticity was severely altered through the breakdown region and the spanwise distribution of axial vorticity present appeared to reach a maximum immediately preceding breakdown. The local concentration of axial vorticity about the vortex axis was reduced while the magnitude of the azimuthal vorticity decreased throughout the breakdown zone. The axial vorticity components with a negative sense, found in the secondary vortex, remained unaffected by changes in wing sweep or angle of attack, in direct contrast to the positive components. The inclusion of the local wing geometry into a previously derived correlation parameter indicated that the circulation of growing leading edge vortex flows were similar at corresponding radii from the vortex axis. It was concluded that the flow over a delta wing, upstream of the breakdown regions and away from the apex and trailing edge regions, is conical. In addition, the dominating factors leading to the onset of breakdown are felt to be the local circulation of the vortex and the accompanying pressure field.

Visser, Kenneth D.↗

Efficient multitasking of Choleski matrix factorization on CRAY supercomputers

A Choleski method is described and used to solve linear systems of equations that arise in large scale structural analysis. The method uses a novel variable-band storage scheme and is structured to exploit fast local memory caches while minimizing data access delays between main memory and vector registers. Several parallel implementations of this method are described for the CRAY-2 and CRAY Y-MP computers demonstrating the use of microtasking and autotasking directives. A portable parallel language, FORCE, is used for comparison with the microtasked and autotasked implementations. Results are presented comparing the matrix factorization times for three representative structural analysis problems from runs made in both dedicated and multi-user modes on both computers. CPU and wall clock timings are given for the parallel implementations and are compared to single processor timings of the same algorithm.

Overman, Andrea L.↗

Parallel variable-band Choleski solvers for computational structural analysis applications on vector multiprocessor supercomputers

A Choleski method used to solve linear systems of equations that arise in large scale structural analyses is described. The method uses a novel variable-band storage scheme and is structured to exploit fast local memory caches while minimizing data access delays between main memory and vector registers. Several parallel implementations of this method are described for the CRAY-2 and CRAY Y-MP computers demonstrating the use of microtasking and autotasking directives. A portable parallel language, FORCE, is also used for two different parallel implementations, demonstrating the use of CRAY macrotasking. Results are presented comparing the matrix factorization times for three representative structural analysis problems from runs made in both dedicated and multi-user modes on both the CRAY-2 and CRAY Y-MP computers. CPU and wall clock timings are given for the various parallel methods and are compared to single processor timings of the same algorithm. Computation rates over 1 GIGAFLOP (1 billion floating point operations per second) on a four processor CRAY-2 and over 2 GIGAFLOPS on an eight processor CRAY Y-MP are demonstrated as measured by wall clock time in a dedicated environment. Reduced wall clock times for the parallel methods relative to the single processor implementation of the same Choleski algorithm are also demonstrated for runs made in multi-user mode.

Poole, E. L.↗

Low-frequency mobility response functions for the central plasma sheet with application to tearing modes

Consideration is given to the effect of constant cross-tail magnetic field By on the collisionless conductivity produced by chaotic scattering and stochastic diffusion of particles in the current sheet for a parabolic geometry. It is shown that the correlation time scales as (By/Bz)-squared, and from this strong By scaling a strong tendency toward stabilization of the linear tearing modes with increasing values of By is inferred. This effect of increased dawn-dusk mobility is particularly dramatic when electrons are introduced in the calculation, and is in agreement with the results of kinetic particle simulations. The collisionless conductivity is expressed in terms of the ensemble-averaged power spectrum of the single particle trajectories, which makes it possible to calculate directly the linear conductivity instead of deriving it from the calculation of the irreversible heating rates.

Hernandez, J.↗

The embedded young stars in the Taurus-Auriga molecular cloud. II - Models for scattered light images

We describe NIR imaging observations of embedded young stars in the Taurus-Auriga molecular cloud. We find a large range in J-K and H-K colors for these class I sources. The bluest objects have colors similar to the reddest T Tauri stars in the cloud; redder objects lie slightly above the reddening line for standard ISM dust and have apparent K extinctions of up to 5 mag. Most of these sources also show extended NIR emission on scales of 10-20 arcsec which corresponds to linear sizes of 1500-3000 AU. The NIR colors and nebular morphologies for this sample and the magnitude of linear polarization in several sources suggest scattered light produces most of the NIR emission in these objects. We present modeling results that suggest mass infall rates that agree with predictions for cold clouds and are generally consistent with rates estimated from radiative equilibrium models. For reasonable dust grain parameters, the range of colors and extinctions require flattened density distributions with polar cavities evacuated by bipolar outflows. These results support the idea that infall and outflow occur simultaneously in deeply embedded bipolar outflow sources. The data also indicate fairly large centrifugal radii and large inclinations to the rotational axis for a typical source.

Kenyon, Scott J.↗

Genetic and molecular dosimetry of HZE radiation (US-1 RADIAT)

In order to estimate radiation exposure in space, experiments were conducted during the 1st International Microgravity Laboratory (IML-1) mission in order to isolate genetic changes in animal cells caused by cosmic rays. The space measurements were evaluated against results from synthetic cosmic rays produced by particle accelerators on the ground. The biological material used was the tiny soil nematode, Caenorhabditis elegans. The measurements were made by thermoluminescent detectors and plastic nuclear track detectors. The development and the chromosome mechanics in microgravity were studied, and the mutagenesis induced by radiation exposure was analyzed. The results showed that there are no obvious differences in the development, behavior and chromosome mechanics, as a function of gravity unloading (reproduction, self-fertilization and mating of males with hermaphrodites, gross anatomy, symmetry and gametogenesis, pairing, disjoining and recombination of chromosomes). A variety of mutants were isolated, and it was noted that mutants isolated from regions of identified high particles were more severely affected than those isolated by random screening. Linear energy transfer particles seem to favor large scale genetic lesions.

Nelson, Gregory A.↗

Internal velocity and mass distributions in simulated clusters of galaxies for a variety of cosmogonic models

The mass and velocity distributions in the outskirts (0.5-3.0/h Mpc) of simulated clusters of galaxies are examined for a suite of cosmogonic models (two Omega(sub 0) = 1 and two Omega(sub 0) = 0.2 models) utilizing large-scale particle-mesh (PM) simulations. Through a series of model computations, designed to isolate the different effects, we find that both Omega(sub 0) and P(sub k) (lambda less than or = 16/h Mpc) are important to the mass distributions in clusters of galaxies. There is a correlation between power, P(sub k), and density profiles of massive clusters; more power tends to point to the direction of a stronger correlation between alpha and M(r less than 1.5/h Mpc); i.e., massive clusters being relatively extended and small mass clusters being relatively concentrated. A lower Omega(sub 0) universe tends to produce relatively concentrated massive clusters and relatively extended small mass clusters compared to their counterparts in a higher Omega(sub 0) model with the same power. Models with little (initial) small-scale power, such as the hot dark matter (HDM) model, produce more extended mass distributions than the isothermal distribution for most of the mass clusters. But the cold dark matter (CDM) models show mass distributions of most of the clusters more concentrated than the isothermal distribution. X-ray and gravitational lensing observations are beginning providing useful information on the mass distribution in and around clusters; some interesting constraints on Omega(sub 0) and/or the (initial) power of the density fluctuations on scales lambda less than or = 16/h Mpc (where linear extrapolation is invalid) can be obtained when larger observational data sets, such as the Sloan Digital Sky Survey, become available.

Cen, Renyue↗

The Sun's Impact on Climate

We provide an overview of the impact of the Sun on the Earth atmosphere and climate system, focused on heating of Earth's atmosphere and oceans. We emphasize the importance of the spectral measurements of SIM and SOLSTICE- that we must know how solar variations are distributed over ultraviolet, visible, and infrared wavelengths, since these have separate characteristic influences on Earth's ozone layer, clouds, and upper layers of the oceans. Emphasis is also given to understanding both direct and indirect influences of the Sun on the Earth, which involve feedbacks between Earth's stratosphere, troposphere, and oceans, each with unique time scales, dynamics, chemistry, and biology, interacting non-linearly. Especially crucial is the role of all three phases of water on Earth, water vapor being the primary greenhouse gas in the atmosphere, the importance of trace gases such as CO2 arising from their absorption in the "water vapor window" at 800 - 1250/cm (12.5 to 8 microns). Melting of polar ice is one major response to the post-industrial global warming, enhanced due to "ice-albedo" feedback. Finally, water in liquid form has a major influence due to cloud albedo feedback, and also due to the oceans' absorption of solar radiation, particularly at visible wavelengths, through the visible "liquid water window" that allows penetration of visible light deep into the mixed layer, while nearby ultraviolet and infrared wavelengths do not penetrate past the upper centimeter ocean surface skin layer. A large fraction of solar energy absorbed by the oceans goes into the latent heat of evaporation. Thus the solar heating of the atmosphere-ocean system is strongly coupled through the water cycle of evaporation, cloud formation, precipitation, surface runoff and ice formation, to Earth's energy budget and climate, each different climate component responding to variations in different solar spectral bands, at ultraviolet, visible and infrared wavelengths.

Cahalan, Robert↗

Modeling of Interfacial Modification Effects on Thermal Conductivity of Carbon Nanotube Composites

The effect of functionalization of carbon nanotubes on the thermal conductivity of nanocomposites has been studied using a multi-scale modeling approach. These results predict that grafting linear hydrocarbon chains to the surface of a single wall carbon nanotube with covalent chemical bonds should result in a significant increase in the thermal conductivity of these nanocomposites. This is due to the decrease in the interfacial thermal (Kapitza) resistance between the single wall carbon nanotube and the surrounding polymer matrix upon chemical functionalization. The nanocomposites studied here consist of single wall carbon nanotubes in a bulk poly(ethylene vinyl acetate) matrix. The nanotubes are functionalized by end-grafting linear hydrocarbon chains of varying length to the surface of the nanotube. The effect which this functionalization has on the interfacial thermal resistance is studied by molecular dynamics simulation. Interfacial thermal resistance values are calculated for a range of chemical grafting densities and with several chain lengths. These results are subsequently used in an analytical model to predict the resulting effect on the bulk thermal conductivity of the nanocomposite.

Clancy, Thomas C.↗

CoLoRe-2LPT: Lyman-$α$ mock catalogues for the validation of DESI cosmological analyses

The Lyman-$α$ (Ly$α$) forest has become a crucial probe for studying the large-scale structure of the universe at high redshift ($z > 2$), providing powerful constraints on Baryon Acoustic Oscillations (BAO) and the full-shape (FS) clustering of matter. As a key ingredient for upcoming BAO and FS analyses, we present a new generation of fast cosmological Ly$α$ mocks based on second-order Lagrangian perturbation theory (2LPT). These new mocks significantly improve upon previous log-normal approaches, both at accurately capturing small scale clustering and at recovering the non-linear broadening of the BAO peak. They are able to reproduce Ly$α$ statistics within $10\%$ of the latest DESI measurement; including the Ly$α$ bias and the redshift-space distortion $β$ parameter, mean transmitted flux, and 1D power spectrum. The corresponding quasar (QSO) clustering is also improved with respect to previous approaches, calibrated against high-resolution Abacus simulations, recovering the observational QSO linear bias to less than $5\%$ and improving redshift-space distortions via 2LPT velocities and the addition of Fingers-of-God effects. Furthermore, these mocks incorporate high column density systems and metal lines, allowing us to explore the effects and systematics induced by these astrophysical contaminants. This new set of mocks has been key for enhancing the modeling and validation of the DESI DR2 Ly$α$ full shape cosmological analysis. This work provides a physically motivated and computationally efficient tool for simulating current and next-generation Ly$α$ surveys and validating FS and BAO analysis.

Bernal, M.F. Ruiz-Herrera [Madrid, CIEMAT] (ORCID:↗

Stochastic Model for High Temperature Oxidation of Cr–Ni Austenitic Steels Assisted by Spallation

Abstract Cr–Ni austenitic steels offer significant high temperature corrosion protection by forming a surface oxide layer. However, above critical service conditions (temperature, atmosphere, thermal cycling), oxidized surface can experience intensive degradation because of scale spallation, which could be detrimental to the in-service life. To predict the effect of scale spallation on oxidation kinetics, a simulation was implemented using a stochastic model. The model considers topological parameters and intensity of spallation which can occur, while delivering a true oxidation constant. The experimental procedure identified the amount of formed spalled scale and topology of spallation based on the use of element mapping of the surface. This information was used to determine a true kinetic constant for a corresponding spallation intensity in oxidized Cr–Ni austenitic steel. To illustrate the capability of the stochastic model, a parametric analysis was performed. The model verified how the spallation parameters could change the oxidation processes from parabolic growth of an adhered oxide layer without spallation to a mixed linear-parabolic, or with a constant thickness of residual scale at high spallation intensity. The spallation model will be used in a separate article to characterize high temperature surface degradation of several Cr–Ni austenitic steels during harsh oxidation environments.

36 MATERIALS SCIENCE↗

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗

Retrieving the Balanced Winds on the Globe as a Generalized Inverse Problem

A generalized inverse technique is applied to retrieve two types of balanced winds that characterize the large-scale dynamics of the atmosphere: rotational winds based upon the linear balance equation, and divergent winds based upon the vorticity budget equation. Both balance equations are singular at or near the equator. The balance equations are transformed in spherical harmonic function space to an under-determined system, for which the scale-weighed least-squares solution consists of a sum of principal and singular components. The principal components represent the response to the source function for the regular eigenmodes, while the singular components are determined by the projection of an independent measurement on the singular eigenmodes. The method was tested with the NCEP/NCAR reanalysis data in which a quasi-balance condition exists. A realistic balanced wind field is retrievable when the singular components are computed based upon the reanalyzed wind data.

Lu, Huei-Iin↗

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗