Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HOMOGENIZATION METHODS”

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 523 records · Page 29

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗

Topics in spectral methods

After detailing the construction of spectral approximations to time-dependent mixed initial boundary value problems, a study is conducted of differential equations of the form 'partial derivative of u/partial derivative of t = Lu + f', where for each t, u(t) belongs to a Hilbert space such that u satisfies homogeneous boundary conditions. For the sake of simplicity, it is assumed that L is an unbounded, time-independent linear operator. Attention is given to Fourier methods of both Galerkin and pseudospectral method types, the Galerkin method, the pseudospectral Chebyshev and Legendre methods, the error equation, hyperbolic partial differentiation equations, and time discretization and iterative methods.

Gottlieb, D.↗

Quadrature imposition of compatibility conditions in Chebyshev methods

This paper handles a problem encountered in the solution by way of Chebyshev spectral collocation discretization of the Poisson equation with homogeneous Neumann boundary conditions. When pseudospectral Chebyshev methodology is utilized to discretize the partial differential equation these integrals have to be approximated by an appropriate quadrature formula. The correct quadrature formula is developed for these problems, taking into account the degree of the polynomials involved. It is shown that this formula leads to a well-conditioned Chebyshev approximation to the differential equations and that the compatibility condition is automatically satisfied.

Gottlieb, D.↗

Remote sensing of optical properties in continuously stratified waters

The radiative transfer equation is solved by Monte Carlo methods for natural waters in which the optical properties are distributed with depth. It is demonstrated that interpreting the reflectance of a continuously stratified ocean in terms of an equivalent homogeneous ocean yields the average of a particular combination of the water's optical properties over the dimensionless penetration depth. Although in general the dimensionless penetration depth cannot be remotely measured, a method is presented for estimating the actual penetration depth from the remote observations if the medium's absorption coefficient is known, independent of depth, and sufficiently large. The application of this to the remote measurement of the vertical distribution of suspended sediments is discussed in detail.

Gordon, H. R.↗

Advanced single crystal for SSME turbopumps

The objective of this program was to evaluate the influence of high thermal gradient casting, hot isostatic pressing (HIP) and alternate heat treatments on the microstructure and mechanical properties of a single crystal nickel base superalloy. The alloy chosen for the study was PWA 1480, a well characterized, commercial alloy which had previously been chosen as a candidate for the Space Shuttle Main Engine high pressure turbopump turbine blades. Microstructural characterization evaluated the influence of casting thermal gradient on dendrite arm spacing, casting porosity distribution and alloy homogeneity. Hot isostatic pressing was evaluated as a means of eliminating porosity as a preferred fatigue crack initiation site. The alternate heat treatment was chosen to improve hydrogen environment embrittlement resistance and for potential fatigue life improvement. Mechanical property evaluation was aimed primarily at determining improvements in low cycle and high cycle fatigue life due to the advanced processing methods. Statistically significant numbers of tests were conducted to quantitatively demonstrate life differences. High thermal gradient casting improves as-cast homogeneity, which facilitates solution heat treatment of PWA 1480 and provides a decrease in internal pore size, leading to increases in low cycle and high cycle fatigue lives.

Fritzemeier, L. G.↗

Direct numerical simulations of activation and deactivation in turbulent atmospheric clouds

Significant knowledge gaps remain in our understanding of turbulence–cloud–aerosol interactions in the Earth's atmosphere, and direct numerical simulation (DNS) has increasingly become an indispensable tool to fill such critical knowledge gaps. Here, this study is an extension of our previous DNS model [Gao et al., J. Geophys. Res.: Atmos., 123(4), 2194–2214 (2018)], with a focus on the activation of aerosol particles into cloud droplets and deactivation of cloud droplets into aerosol particles in a microscale cloud environment. The effects of turbulence intensity, particle curvature, and solute, as well as the initial distributions of the aerosol particles (monodisperse or polydisperse) are investigated. The governing equations for the flow of air, temperature, and water vapor mixing ratio are solved numerically in the Eulerian fashion, assuming homogeneous and isotropic turbulence. The dynamics of the aerosol and cloud particles are calculated with the Lagrangian particle tracking method. The results show that the deviations of the thermodynamic variables from their respective means are significantly reduced, the activation process appears to be delayed, and the deactivation process occurs more rapidly, as the turbulence intensity is increased. The inclusion of particle curvature and solute effects, as well as polydispersity, tends to retard the activation of aerosols into cloud droplets. It is also observed that fluctuations in supersaturation broaden the spread of particle radii, and the broadening is followed by a narrowing as turbulent homogenization reduces thermodynamic fluctuations over time.

54 ENVIRONMENTAL SCIENCES↗

Assessment of mechanical, thermal properties and crystal shapes of monoclinic tricalcium silicate from atomistic simulations

The two most common polymorphs in industrial alite, M1 and M3, were characterized at the molecular scale. Different methods were employed and discussed to assess mechanical properties and specific heat of both polymorphs. The calculated homogenized elastic moduli and specific heat were found in good agreement with experimental measurements. A comparative analysis of spacial Youngs modulus reveal isotropic and anisotropic spacial distribution for M{sub 1} and M{sub 3} respectively. A more isotropic compressive strength is also reported for M{sub 1} when compared to M{sub 3} polymorph. Cleavage energies computation allowed to proposed equilibrium shapes for both polymorph, with significant differences. While the lowest cleavage energies were found along (100) and (001) for both polymorphs, the constructed M1 crystal possesses 3 independent facets, against seven for the M3 polymorph.

36 MATERIALS SCIENCE↗

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗

Solid Phase Gradient Alloying Method via ShAPE

Alloying has been used to confer desirable properties to various metallic systems since the bronze age. For example, Ni, Mn, and Cr are alloyed into an Fe matrix to make steel, which has improved strength and corrosion resistance (as displayed in Fig. 1). In addition to the type of alloying elements, the amount of alloying elements added into a matrix is also very critical for a material’s mechanical property, physical property, machinability, and cost. Conventionally, new alloys are discovered by combining alloy precursors and elemental constituents into a mixture and melting them together. The solidified product or ingot represents a usually non-homogeneous combination of the chemistry with a microstructure dictated by the physics of solidification. It is time and energy intensive to optimize the alloying amounts via the casting process and hence, economically not conducive to rapid alloy discovery. In addition, the cast microstructure seen in the as-cast ingot is often not the ideal, nor even desired structure for optimized performance. Cast materials often need further processing such as homogenization, annealing, and deformation work put in through rolling, forging, extruding, etc. to have favorable microstructures and properties. A faster method to discover new alloy combinations and evaluate them in their worked or processed form is needed. Recently, Laser Engineered Net Shaping was applied to fabricate gradient compositional material for alloy designing. However, the mechanical properties of this gradient alloy were poor due to the existence of impurities, oxidation and cracks that are inherent in the melt-solidification process. Because of these limitations, we propose to use a solid-phase process (ShAPE) to create bulk materials (extrudates) that vary in composition from one end of the extruded solid to the other. Various manufacturing methods for alloy design including conventional casting method, laser based combinatorial method and current solid phase gradient alloying method are displayed and compared in Fig. 2. ShAPE machine and a schematic of ShAPE process are displayed in Fig. 3. Starting material is processed by a rotating and plunging die to form an extrudate. We are proposing through this methodology to invent a new, bulk-scale combinatorial technique. With this novel technique, alloying element can be dissolved into the matrix with a continued gradient without bulk melting. Formation of inter-metallics and defects originating from melt processing can be avoided. In addition, the ShAPE process creates the mixed alloy chemistry, and meanwhile subjects the material to severe plastic strain, which induces the favorable “worked” microstructure. Products from the ShAPE process can be tested in hardness directly providing a path to rapid evaluation of mechanical properties. The ability to create graded structures using friction stir processing (FSP) has been demonstrated, however using the ShAPE process we believe will lead to much faster and higher fidelity results. When combined with high throughput screening methods of physical, mechanical and microstructural property characterization, efficiency and accuracy of alloy design can be significantly improved.

36 MATERIALS SCIENCE↗

Probing basal planes and edge sites in polygranular nuclear graphite by gas adsorption: Estimation of active surface area

Polygranular nuclear graphite, manufactured at 2800–3000 °C from a carbonized filler and binder mix, has high graphitization degree, crystalline structure close to perfect graphite, and about 20% porosity. The pore surfaces expose large regions of rough, defective prismatic edges of graphite crystallites which are the locus of graphite materials surface sites active for oxidation, chemisorption, and electron transfer. However, we show that high-resolution N 2 and Kr first monolayer adsorption on polygranular graphite (P/P 0 < 0.015) occurs in many ways like adsorption on graphitized carbon blacks. This proves the presence of energetically homogeneous basal planes domains in graphite porosity, which was not fully acknowledged before. Using classical analysis methods (Langmuir, Hill-de Boer, adsorption potential distribution) we quantity the basal plane area (BPA) of several polygranular graphite types and correlate it with their microstructure. Finally, we propose that gas adsorption is uniquely positioned to reliably estimate the active surface area of graphite by subtracting BPA from the BET total surface area. Direct estimation based on adsorption is preferable to indirect calculations based on microstructural information.

36 MATERIALS SCIENCE↗

Optimization-Based Azeotropic Distillation System Synthesis Using Geometric Insights

The synthesis of azeotropic distillation systems is challenging due to the existence of compartments in the residue curve map coupled with the combinatorial aspects from numerous possible system configurations. In this work, an optimization-based approach is introduced to synthesize homogeneous azeotropic distillation systems. The approach employs a network-based representation generated via a matrix method. To design the distillation columns, the modified Underwood equations are adopted, in which pseudocomponent-based compositions are used. For cases where separatrices are significantly curved, corrections via piece-wise linear functions and collinearity properties are implemented. Two examples are presented to illustrate the proposed approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Preparation, analysis, and application of coated glass targets for the Wendelstein 7-X laser blow-off system

Coated glass targets are a key component of the Wendelstein 7-X laser blow-off system that is used for impurity transport studies. The preparation and analysis of these glass targets as well as their performance is examined in this paper. The glass targets have a high laser damage threshold and are coated via physical vapor deposition with μm thick films. In addition, nm-thin layers of Ti are used as an interface layer for improved ablation efficiency and reduced coating stress. Hence, the metallic or ceramic coating has a lateral homogeneity within 2% and contaminants less than 5%, being optimal for laser ablation processing. With this method, a short (few ms) and well defined pulse of impurities with about 1017 particles can be injected close to the last closed flux surface of Wendelstein 7-X. In particular, a significant amount of atoms with a velocity of about 1 km/s enters the plasma within 1 ms. The atoms are followed by a negligible concentration of slower clusters and macro-particles. This qualifies the use of the targets and applied laser settings for impurity transport studies with the laser blow-off system in Wendelstein 7-X.

47 OTHER INSTRUMENTATION↗

The role of breakup and coalescence in fine-scale bubble-induced turbulence. II. Kinematics

This second part of our research explores the kinematic aspect of fine-scale bubble-induced turbulence (BIT) to (i) present the effect of bubble breakup and coalescence and (ii) compare it against the universal kinematic fine-scale turbulence characteristics reported in the literature. To this end, we simulate a dilute bubbly system of 0.5% void fraction using two distinct numerical simulations. In the volume-of-fluid (VoF) simulation, bubbles undergo breakup and coalescence. In the immersed boundary method (IBM) simulation, however, they act as rigid spheres. We also perform a simulation of classical homogeneous isotropic turbulence (HIT). The first important outcome of this study is that BIT is radically different from HIT in terms of its kinematic fine-scale characteristics. In the vorticity-dominating regions, BIT exhibits a weak vortex stretching. This weak vortex stretching is due to (a) the intermediate strain-rate eigenvalues skewed weakly to positive and (b) the extensive strain-rate eigenvector aligning perpendicular to the vorticity vector. The BIT has, on average, not only a weak enstrophy production but also a weak strain production in strain-dominating regions. The weak strain production is due to (a) the presence of vortex stretching in highly strained fluid elements and (b) the absolute magnitude of compressive strain-rate eigenvalue being as close to the extensive strain-rate eigenvalue. Thus, none of the kinematic fine-scale HIT characteristics is noted for BIT. The second important conclusion is that bubble breakup and coalescence play little to no influence on the kinematics of fine-scale BIT as VoF and IBM simulations produce similar results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Radiation and scattering by thin-wire structures in the complex frequency domain

Piecewise-sinusoidal expansion functions and Galerkin's method are employed to formulate a solution for an arbitrary thin-wire configuration in a homogeneous conducting medium. The analysis is performed in the real or complex frequency domain. In antenna problems, the solution determines the current distribution, impedance, radiation efficiency, gain and far-field patterns. In scattering problems, the solution determines the absorption cross section, scattering cross section and the polarization scattering matrix. The electromagnetic theory is presented for thin wires and the forward-scattering theorem is developed for an arbitrary target in a homogeneous conducting medium.

Richmond, J. H.↗

Generation of ULF waves by electric or magnetic dipoles

The generation of ULF waves by ground-based magnetic and electric dipoles is studied with a simplified model consisting of three adjoining homogeneous regions representing the groud, the vacuum (free space) region, and the ionosphere. The system is assumed to be immersed in a homogeneous magnetic field with an arbitrary tilt angle. By the use of Fourier techniques and the method of stationary phase, analytic expressions are obtained for the field strength of the compressional Alfven waves in the ionosphere. Expressions are also obtained for the strength of the torsional Alfven wave in the ionosphere and the ULF magnetic field at ground level. Numerical results are obtained for the compressional Alfven-wave field strength in the ionosphere with a nonvertical geomagnetic field and for the ULF magnetic field at ground level for a vertical geomagnetic field.

Harker, K. J.↗

Sulfur cycling and metabolism of phototrophic and filamentous sulfur bacteria

Phototrophic sulfur bacteria taken from different habitate (Alum Rock State Park, Palo Alto salt marsh, and Big Soda Lake) were grown on selective media, characterized by morphological and pigment analysis, and compared with bacteria maintained in pure culture. A study was made of the anaerobic reduction of intracellular sulfur globules by a phototrophic sulfur bacterium (Chromatium vinosum) and a filamentous aerobic sulfur bacterium (Beggiatoa alba). Buoyant densities of different bacteria were measured in Percoll gradients. This method was also used to separate different chlorobia in mixed cultures and to assess the relative homogeneity of cultures taken directly or enriched from natural samples (including the purple bacterial layer found at a depth of 20 meters at Big Soda Lake.) Interactions between sulfide oxidizing bacteria were studied.

Guerrero, R.↗

Significance of foundation-soil separation in dynamic soil-structure interaction

THe dynamic response of flexible surface strip-foundations allowed to uplift is numerically obtained for externally applied forces of a transient time variation. The soil medium is represented by an isotropic, homogeneous and linear half-space. The soil is treated by a time domain boundary element method, while the flexible foundation is treated by the finite element method. It was concluded that intermediate relative stiffness leads to moderate deformations when uplift is permitted. Very flexible footings produce higher deformations in unilateral contact compared to bilateral contact, and thus should be considered in their design. Unilateral contact does not significantly increase deformations for stiff footings subjected to concentrated central loading. However, relatively large deformation differences occur when the loading is eccentric, necessitating consideration of uplift in their design.

Spyrakos, Constantine C.↗

AIS-2 radiometry and a comparison of methods for the recovery of ground reflectance

A field experiment and its results involving Airborne Imaging Spectrometer-2 data are described. The radiometry and spectral calibration of the instrument are critically examined in light of laboratory and field measurements. Three methods of compensating for the atmosphere in the search for ground reflectance are compared. It was found that laboratory determined responsitivities are 30 to 50 percent less than expected for conditions of the flight for both short and long wavelength observations. The combined system atmosphere surface signal to noise ratio, as indexed by the mean response divided by the standard deviation for selected areas, lies between 40 and 110, depending upon how scene averages are taken, and is 30 percent less for flight conditions than for laboratory. Atmospheric and surface variations may contribute to this difference. It is not possible to isolate instrument performance from the present data. As for methods of data reduction, the so-called scene average or log-residual method fails to recover any feature present in the surface reflectance, probably because of the extreme homogeneity of the scene.

Conel, James E.↗