Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor sampling”

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.

33 records · Page 2

Instrumentation of the variable-angle magneto-optic ellipsometer and its application to M-O media and other non-magnetic films

A new and comprehensive dielectric tensor characterization instrument is presented for characterization of magneto-optical recording media and non-magnetic thin films. Random and systematic errors of the system are studied. A series of TbFe, TbFeCo, and Co/Pt samples with different composition and thicknesses are characterized for their optical and magneto-optical properties. The optical properties of several non-magnetic films are also measured.

Zhou, Andy F.↗

Correlation function study for sea ice

For active and passive microwave remote sensing of sea ice, a correlation function of exponential form is extracted from the photograph of a horizontal thin section taken from a sample of artificially grown saline ice that closely resembled Arctic congelation sea ice. It is found that the extracted correlation lengths are consistent with the published average size of brine pockets. With the application of strong fluctuation theory and the bilocal approximation, the effective permittivity tensor is derived in the low-frequency limit for an unbounded uniaxial random medium with two-phase mixtures. Using the extracted correlation lengths, the effective permittivity tensor is computed as a function of fractional volume of brine inclusions and compared with in situ measurements at 4.8 and 9.5 GHz.

Lin, F. C.↗

Preliminary Modeling of CH3D From 4000 to 4550 cm (exp-1)

A new study of 12CH3D line positions and intensities was performed for the upper portion of the Enneadecad polyad between 4000 and 4550/cm. For this, FTIR spectra were recorded with D-enriched methane samples (at 80 K with a Bruker 125 IFS at 0.005/cm resolution and at 291 K with the McMath-Pierce FTS at 0.011/cm resolution, respectively). Line positions and intensities were retrieved by least square curve-fitting procedures and analyzed using the effective Hamiltonian and the effective dipole moment expressed in terms of irreducible tensor operators adapted to symmetric top molecules. Initially, only the cold spectrum was used to identify quantum assignments and predict 12CH3D relative intensities in this region. To assign higher quanta up to J equal 14, additional line positions and intensities were obtained from two room temperature spectra. In the final stage, measured intensities from both the cold and the room temperature data were normalized to corresponding values at 296 K and then averaged. Combining the two temperature datasets confirmed the assumed quantum assignments and also demonstrated the relative accuracies to be better than 70.0002/cm for line positions and at least +/-6% for intensities so that approx.1160 features were selected. Including additional assignments from the room temperature spectra alone permitted 1362 line intensities of 11 bands (involving 23 vibrational symmetry components) to be reproduced with an RMS of 9%. Over 4085 selected positions for 12 bands were modeled to 0.008/cm. Nevertheless a number of known assignments could not be modeled to within our experimental precisions. More work is needed to obtain a complete characterization of this complex polyad.

Nikitin, A. V.↗

Two Dimensional Symmetric Correlation Functions of the S Operator and Two Dimensional Fourier Transforms: Considering the Line Coupling for P and R Lines of Linear Molecules

The refinement of the Robert-Bonamy (RB) formalism by considering the line coupling for isotropic Raman Q lines of linear molecules developed in our previous study [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)] has been extended to infrared P and R lines. In these calculations, the main task is to derive diagonal and off-diagonal matrix elements of the Liouville operator iS1 − S2 introduced in the formalism. When one considers the line coupling for isotropic Raman Q lines where their initial and final rotational quantum numbers are identical, the derivations of off-diagonal elements do not require extra correlation functions of the ^S operator and their Fourier transforms except for those used in deriving diagonal elements. In contrast, the derivations for infrared P and R lines become more difficult because they require a lot of new correlation functions and their Fourier transforms. By introducing two dimensional correlation functions labeled by two tensor ranks and making variable changes to become even functions, the derivations only require the latters' two dimensional Fourier transforms evaluated at two modulation frequencies characterizing the averaged energy gap and the frequency detuning between the two coupled transitions. With the coordinate representation, it is easy to accurately derive these two dimensional correlation functions. Meanwhile, by using the sampling theory one is able to effectively evaluate their two dimensional Fourier transforms. Thus, the obstacles in considering the line coupling for P and R lines have been overcome. Numerical calculations have been carried out for the half-widths of both the isotropic Raman Q lines and the infrared P and R lines of C2H2 broadened by N2. In comparison with values derived from the RB formalism, new calculated values are significantly reduced and become closer to measurements.

infrared radiation↗

Stress Concentration and Fracture at Inter-variant Boundaries in an Al-Li Alloy

Delamination fracture has limited the use of lightweight Al-Li alloys. Studies of secondary, delamination cracks in alloy 2090, L-T fracture toughness samples showed grain boundary failure between variants of the brass texture component. Although the adjacent texture variants, designated B(sub s1) and B(sub s2), behave similarly during rolling, their plastic responses to mechanical tests can be quite different. EBSD data from through-thickness scans were used to generate Taylor factor maps. When a combined boundary normal and shear tensor was used in the calculation, the delaminating grains showed the greatest Taylor Factor differences of any grain pairs. Kernel Average Misorientation (KAM) maps also showed damage accumulation on one side of the interface. Both of these are consistent with poor slip accommodation from a crystallographically softer grain to a harder one. Transmission electron microscopy was used to confirm the EBSD observations and to show the role of slip bands in the development of large, interfacial stress concentrations. A viewgraph presentation accompanies the provided abstract.

Crooks, Roy↗

An orthogonal coordinate grid following the three-dimensional viscous flow over a concave surface

Swept wings designed for laminar flow control exhibit both centrifugal and crossflow instabilities which produce streamwise vortices that can lead to early transition from laminar to turbulent flow in the presence of Tollmien-Schlichting waves. This paper outlines an iterative algorithm for generation of an orthogonal, curvilinear, coordinate grid following the streamlines of the three-dimensional viscous flow over a swept, concave surface. The governing equations for the metric tensor are derived from the Riemann-Christoffel tensor for an Euclidian geometry. Unit vectors along streamline, normal and binormal directions are determined. The governing equations are not solved directly, but are employed only as compatibility equations. The scale factor for the streamline coordinate is obtained by an iterative integration scheme on a 200 x 100 x 5 grid, while the other two scale factors are determined from definitions. Sample results are obtained which indicate that the compatibility equation error decreases linearly with grid step size. Grids smaller than 200 x 100 x 5 are found to be inadequate to resolve the grid curvature.

Dagenhart, J. R↗

Statistical analysis of near-wall structures in turbulent channel flow

The two-point correlation tensor is used to analyze near-wall structures in turbulent channel flow. Simulation results have been utilized to obtain the correlation tensor with sufficient spatial resolution to permit accurate differentiation for computing the vorticity correlation, and the resolution of sharp features such as the shear layer and its corresponding region of spanwise vorticity. Proper orthogonal decomposition and stochastic estimation were observed to yield similar results when appropriate conditions were employed for the estimation. Resulting structures were such that they would be detected by standard conditional sampling methods. The results of decomposing the vorticity field and the velocity field were significantly different, indicating that if a coherent structure is found that dominates the velocity fields, the curl of that velocity structure will not dominate the vorticity field.

Moser, Robert D.↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Multi-Scale Computational Modeling of Two-Phased Metal Using GMC Method

A multi-scale computational model for determining plastic behavior in two-phased CMSX-4 Ni-based superalloys is developed on a finite element analysis (FEA) framework employing crystal plasticity constitutive model that can capture the microstructural scale stress field. The generalized method of cells (GMC) micromechanics model is used for homogenizing the local field quantities. At first, GMC as stand-alone is validated by analyzing a repeating unit cell (RUC) as a two-phased sample with 72.9% volume fraction of gamma'-precipitate in the gamma-matrix phase and comparing the results with those predicted by finite element analysis (FEA) models incorporating the same crystal plasticity constitutive model. The global stress-strain behavior and the local field quantity distributions predicted by GMC demonstrated good agreement with FEA. High computational saving, at the expense of some accuracy in the components of local tensor field quantities, was obtained with GMC. Finally, the capability of the developed multi-scale model linking FEA and GMC to solve real life sized structures is demonstrated by analyzing an engine disc component and determining the microstructural scale details of the field quantities.

nickel alloys↗

Enhanced Quantification of 3D Woven Composites Via Fourier Analysis and Structure Tensors Applied to CT Scans

NASA's development of advanced 3D woven carbon composites, such as HEEET and 3MDCP, represents a significant leap in materials designed to endure extreme atmospheric entry conditions. However, the internal complexity of these materials introduces a variety of structures and novel defect types. Traditional characterization methods rely on labor-intensive, manual examination of X-ray CT scans, demanding considerable time and expertise. This process, while thorough, is not sustainable for large-scale analysis and defect tracking. Addressing this critical challenge, this presentation introduces an innovative suite of automated tools and analytical techniques, pioneered by the Mars Sample Return (MSR) Earth Entry System (EES) Thermal Protection System (TPS) project. These advancements offer a more efficient and accurate approach for quantifying CT data and characterizing the defects of 3D woven composites.

Magnus Haw↗

Quantitative Characterization of 3D Woven Composites Via Fourier Analysis and Structure Tensors Applied to CT Scans

NASA's development of advanced 3D woven carbon composites, such as HEEET and 3MDCP, represents a significant leap in materials designed to endure extreme atmospheric entry conditions. However, the internal complexity of these materials introduces a variety of structures and novel defect types. Traditional characterization methods rely on labor-intensive, manual examination of X-ray CT scans, demanding considerable time and expertise. This process, while thorough, is not sustainable for large-scale analysis and defect tracking. Addressing this critical challenge, this presentation introduces an innovative suite of automated tools and analytical techniques, pioneered by the Mars Sample Return (MSR) Earth Entry System (EES) Thermal Protection System (TPS) project. These advancements offer a more efficient and accurate approach for quantifying CT data and characterizing the defects of 3D woven composites.

3D Woven↗

An Experimental Investigation of a Wing-Fuselage Junction Model in the NASA Langley 14- by 22-Foot Subsonic Tunnel

Current turbulence models, such as those employed in Reynolds-averaged Navier-Stokes CFD, are unable to reliably predict the onset and extent of the three-dimensional separated flow that typically occurs in wing-fuselage junctions. To critically assess, as well as to improve upon, existing turbulence models, experimental validation-quality flow-field data in the junction region is needed. In this report, we present an overview of experimental measurements on a wing-fuselage junction model that addresses this need. The experimental measurements were performed in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The model was a full-span wing-fuselage body that was configured with truncated DLR-F6 wings, both with and without leading-edge extensions at the wing root. The model was tested at a fixed chord Reynolds number of 2.4 million, and angles-of-attack ranging from -10 degrees to +10 degrees were considered. Flow-field measurements were performed with a pair of miniature laser Doppler velocimetry (LDV) probes that were housed inside the model and attached to three-axis traverse systems. One LDV probe was used to measure the separated flow field in the trailing-edge junction region. The other LDV probe was alternately used to measure the flow field in the leading-edge region of the wing and to measure the incoming fuselage boundary layer well upstream of the leading edge. Both LDV probes provided measurements from which all three mean velocity components, all six independent components of the Reynolds-stress tensor, and all ten independent components of the velocity triple products were calculated. In addition to the flow-field measurements, static and dynamic pressures were measured at selected locations on the wings and fuselage of the model, infrared imaging was used to characterize boundary-layer transition, oil-flow visualization was used to visualize the separated flow in the leading- and trailing-edge regions of the wing, and unsteady shear stress was measured at limited locations using capacitive shear-stress sensors. Sample results from the measurement techniques employed during the test are presented and discussed.

Kegerise, Michael A.↗

Prediction of Turbulence-Generated Noise in Unheated Jets: JeNo Technical Manual (Version 1.0) - Part 1

The model-based approach, used by the JeNo code to predict jet noise spectral directivity, is described. A linearized form of Lilley's equation governs the non-causal Green s function of interest, with the non-linear terms on the right hand side identified as the source. A Reynolds-averaged Navier-Stokes (RANS) solution yields the required mean flow for the solution of the propagation Green s function in a locally parallel flow. The RANS solution also produces time- and length-scales needed to model the non-compact source, the turbulent velocity correlation tensor, with exponential temporal and spatial functions. It is shown that while an exact non-causal Green s function accurately predicts the observed shift in the location of the spectrum peak with angle as well as the angularity of sound at low to moderate Mach numbers, the polar directivity of radiated sound is not entirely captured by this Green s function at high subsonic and supersonic acoustic Mach numbers. Results presented for unheated jets in the Mach number range of 0.51 to 1.8 suggest that near the peak radiation angle of high-speed jets, a different source/Green s function convolution integral may be required in order to capture the peak observed directivity of jet noise. A sample Mach 0.90 heated jet is also discussed that highlights the requirements for a comprehensive jet noise prediction model.

Khavaran, Abbas↗

Design of a Low-Cost, Submersible, Digital Holographic Microscope for in Situ Microbial Imaging

The methodologies for studying marine microbiology typically consist of utilizing instrumentation within a laboratory. This typically requires extracting a sample from its place of origin prior to examination, which may be days after collection. Oftentimes, the solution is to bring the lab to the ocean, which may be costly and provide further limitations for a sterile and stable laboratory environment. Here we present a low-cost, submersible, digital holographic microscope (DHM) designed to image marine microorganisms (such as bacteria and plankton) in their natural underwater environment. Our instrument eliminates the need to transport samples and allows for instantaneous data collection of microbes in-situ. The DHM achieves sub-micron spatial resolution and is paired with artificial intelligence for the detection and tracking of specimens to reduce the overall collected data. This instrument also aims to reduce the cost of manufacturing and field expenses relative to marine microbiological research. “Off the shelf” components were selected in the design process of this instrument which allows us to achieve precise results without sacrificing data quality. The DHM itself costs under one thousand dollars and features a low-cost high-resolution camera, the Arducam MT9J001. Included in the design were five main subsystems: optical, mechanical, electrical, power, and machine learning. Our on board computer and artificial intelligence consist of a Raspberry Pi 4 (8 Gb) and Google Coral USB Tensorflow accelerator. Instrument testing has successfully proven our abilities of data acquisition for at least two hours in depths of at least forty meters below sea level. Additionally, our artificial intelligence system is currently capable of tracking up to ten areas of interest in a fraction of a second with over ninety percent confidence via the neural net driven by the tensor cores on the Google Coral. Furthermore, we have demonstrated the versatility of our instrument by mounting it on an ocean-going ROV, the BlueRov2 by BlueRobotics. Our tests on the BlueRov2 exemplified the cost-effective nature of a submersible and reusable instrument that can be implemented in moderate environments and on most vessels.

Wallace, James Kent↗