SEARCH · Engineering Papers
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.
ZIP3D analyses cracks in three-dimensional elastic or elastic/plastic bodies. Program calculates mixed-mode fracture-mechanics parameters in cracked solids. Bodies without cracks analyzed to obtain stresses, strains, and displacement fields. Written in FORTRAN V.
The mixed-mode bending (MMB) test uses a lever to apply simultaneously mode I and mode II loading to a split-beam specimen. An iterative analysis that accounts for the geometric nonlinearity of the MMB test was developed. The analysis accurately predicted the measured load-displacement response and the strain energy release rate, G, of an MMB test specimen made of AS4/PEEK. The errors in G when calculated using linear theory were found to be as large as 30 percent in some cases. Because it would be inconvenient to use a nonlinear analysis to analyze MMB data, the MMB apparatus was redesigned to minimize the nonlinearity. With the improved apparatus, loads are applied just above the midplane of the test specimen through a roller attached to the lever. This apparatus was demonstrated by measuring the mixed-mode delamination fracture toughhess of the test specimen. The nonlinearity errors associated with testing this tough composite material were less than +/- 3 percent. The data from the improved MMB apparatus analyzed with a linear analysis were similar to those found with the original apparatus and the nonlinear analysis.
A test rig for testing a thick split cantilever beam for scissoring delamination (mode 3) fracture toughness was developed. A 3-D finite element analysis was conducted on the test specimen to determine the strain energy release rate, G, distribution along the delamination front. The virtual crack closure technique was used to calculate the G components resulting from interlaminar tension, GI, interlaminar sliding shear, GII, and interlaminar tearing shear, GIII. The finite element analysis showed that at the delamination front no GI component existed, but a GII component was present in addition to a GIII component. Furthermore, near the free edges, the GII component was significantly higher than the GIII component. The GII/GIII ratio was found to increase with delamination length but was insensitive to the beam depth. The presence of GII at the delamination front was verified experimentally by examination of the failure surfaces. At the center of the beam, where the failure was in mode 3, there was significant fiber bridging. However, at the edges of the beam where the failure was in mode 3, there was no fiber bridging and mode 2 shear hackles were observed. Therefore, it was concluded that the split cantilever beam configuration does not represent a pure mode 3 test. The experimental work showed that the mode 2 fracture toughness, GIIc, must be less than the mode 3 fracture toughness, GIIIc. Therefore, a conservative approach to characterizing mode 3 delamination is to equate GIIIc to GIIc.
Dislocation densities are generated in GaAs single crystals due to the excessive thermal stresses induced by temperature variations during growth. A viscoplastic material model for GaAs, which takes into account the movement and multiplication of dislocations in the plastic deformation, is developed according to Haasen's theory. The dislocation density is expressed as an internal state variable in this dynamic viscoplastic model. The deformation process is a nonlinear function of stress, strain rate, dislocation density and temperature. The dislocation density in the GaAs crystal during vertical Bridgman growth is calculated using a nonlinear finite element model. The dislocation multiplication in GaAs crystals for several temperature fields obtained from thermal modeling of both the GTE GaAs experimental data and artificially designed data are investigated.
We present results of a series of numerical experiments designed to test hypothetical mechanisms that derive deformation in the New Madrid seismic zone. Experiments are constrained by subtle topography and the distribution of seismicity in the region. We use a new boundary element algorithm that permits calcuation of the three-dimensional deformation field. Surface displacement fields are calculated for the New Madrid zone under both far-field (plate tectonics scale) and locally derived driving strains. Results demonstrate that surface displacement fields cannot distinguish between either a far-field simple or pure shear strain field or one that involves a deep shear zone beneath the upper crustal faults. Thus, neither geomorphic nor geodetic studies alone are expected to reveal the ultimate driving mechanism behind the present-day deformation. We have also tested hypotheses about strain accommodation within the New Madrid contractional step-over by including linking faults, two southwest dipping and one vertical, recently inferred from microearthquake data. Only those models with step-over faults are able to predict the observed topography. Surface displacement fields for long-term, relaxed deformation predict the distribution of uplift and subsidence in the contractional step-over remarkably well. Generation of these displacement fields appear to require slip on both the two northeast trending vertical faults and the two dipping faults in the step-over region, with very minor displacements occurring during the interseismic period when the northeast trending vertical faults are locked. These models suggest that the gently dippling central step-over fault is a reverse fault and that the steeper fault, extending to the southeast of the step-over, acts as a normal fault over the long term.
This summer at Kennedy Space Center, I spent 10 weeks as an intern working at the Prototype Development Lab. During this time I learned about the design and machining done here at NASA. I became familiar with the process from where a design begins in Pro/Engineer and finishes at the hands of the machinists. As an intern I was given various small jobs to do and then one project of my own. My personal project was a job for the Applied Physics Lab; in their work they use an X-Ray Backscatter machine. Previously it was resting atop a temporary frame that limited the use of the machine. My job was to design a frame for the machine to rest upon that would allow a full range of sample sizes. The frame was required to support the machine and provide a strain relief for the cords attached to the machine as it moved in the x and y directions. Calculations also had to be done to be sure the design would be able to withstand any loads or outside sources of stress. After the calculations proved the design to be ready to withstand the requirements, the parts were ordered or fabricated, as required. This helped me understand the full process of jobs sent to the Prototype Development Lab.
Efficient, semi-analytical models of bonded joints are utilized for initial joint design and sizing due to their convenience and speed. Many of these models are based on beam and plate theories, which are most applicable when in-plane stresses/strains dominate. While this can be a sensible assumption for obtaining the deformations of the joint, out of plane (transverse) stresses are often of critical importance for failure in bonded joints composed of fiber reinforced composite materials. Since through-thickness strains are neglected in these semi-analytical theories, the present study shows a method to back-calculate the out of plane stresses based on stress balance equations, which use the in-plane stresses and displacements that have already been obtained. The method was implemented within the Joint Element Designer (JED) software, and verification and validation examples have been shown. As is typical for these types of models, corner stress singularities are not captured by the JED model, but the general trends of the stress and deformation fields agreed with a higher fidelity finite element model. As such, this method for obtaining out of plane stresses, when coupled with the JED code, is very useful for the complete stress analysis of bonded composite joints.
A new method is discussed that may be used to correct wind tunnel strain-gage balance load predictions for the influence of residual temperature effects at the location of the strain-gages. The method was designed for the iterative analysis technique that is used in the aerospace testing community to predict balance loads from strain-gage outputs during a wind tunnel test. The new method implicitly applies temperature corrections to the gage outputs during the load iteration process. Therefore, it can use uncorrected gage outputs directly as input for the load calculations. The new method is applied in several steps. First, balance calibration data is analyzed in the usual manner assuming that the balance temperature was kept constant during the calibration. Then, the temperature difference relative to the calibration temperature is introduced as a new independent variable for each strain--gage output. Therefore, sensors must exist near the strain--gages so that the required temperature differences can be measured during the wind tunnel test. In addition, the format of the regression coefficient matrix needs to be extended so that it can support the new independent variables. In the next step, the extended regression coefficient matrix of the original calibration data is modified by using the manufacturer specified temperature sensitivity of each strain--gage as the regression coefficient of the corresponding temperature difference variable. Finally, the modified regression coefficient matrix is converted to a data reduction matrix that the iterative analysis technique needs for the calculation of balance loads. Original calibration data and modified check load data of NASA's MC60D balance are used to illustrate the new method.
Device-quality tensile-strained Ge (ε-Ge) grown on a large bandgap semiconductor with superior electrical and optical carrier confinement is essential for group-IV-based optoelectronics. Properties of ε-Ge active layers synthesized on In 0.24 Ga 0.76 As buffers with two different surface terminations─arsenic-rich and indium-rich─were experimentally demonstrated, highlighting the factors not considered in theoretical calculations. High-resolution X-ray diffraction and Raman spectroscopy analyses of these ε-Ge/In 0.24 Ga 0.76 As heterostructures confirmed the fully strained (1.6%) and partially relaxed (0.82%) nature of the ε-Ge bonded with arsenic-terminated (Ge As-terminated ) and indium-terminated (Ge In-terminated ) In 0.24 Ga 0.76 As stressors, respectively. High-resolution cross-sectional transmission electron microscopy showed a coherent, sharp, and fully strained ε-Ge/In 0.24 Ga 0.76 As heterointerface in the Ge As-terminated heterostructure, whereas microtwin defects were present in the Ge In-terminated heterostructure. These heterostructures were further characterized by evaluating the minority carrier lifetimes, high for Ge As-terminated (525 ns) and low for Ge In-terminated (69 ns), using the photoconductive decay technique. Moreover, band alignment was constructed using X-ray photoelectron spectroscopy, where the Ge As-terminated heterostructure revealed that both holes and electrons were confined within the ε-Ge active layer as a type-I band alignment with ΔE V, As-terminated = 0.22 eV and ΔE C,As-terminated = 0.38 eV. On the other hand, the Ge In-terminated heterostructure exhibited a type-II band alignment with ΔE V,In-terminated = – 0.02 eV and ΔE C,In-terminated = 0.53 eV. Furthermore, the magnetotransport properties revealed high mobility (321 cm 2 /(V s)) with single-electron transport in Ge As-terminated heterostructure and low mobility (3.34 cm 2 /(V s)) with multihole transport in the Ge In-terminated heterostructure. Therefore, preferring the ε-Ge on the arsenic-rich surface of In 0.24 Ga 0.76 As stressor over the indium-rich surface during material synthesis offers device-quality materials with high carrier lifetime and superior carrier confinement, which can provide an opportunity to fabricate efficient group-IV-based optoelectronic devices.
The track structure model of heavy ion cross sections was developed by Katz and co-workers in the 1960s. In this model the action cross section is evaluated by mapping the dose-response of a detector to gamma rays (modeled from biological target theory) onto the radial dose distribution from delta rays about the path of the ion. This is taken to yield the radial distribution of probability for a "hit" (an interaction leading to an observable end-point). Radial integration of the probability yields the cross section. When different response from ions of different Z having the same stopping power is observed this model may be indicated. Since the 1960s there have been several developments in the computation of the radial dose distribution, in the measurement of these distributions, and in new radiobiological data against which to test the model. The earliest model, by Butts and Katz made use of simplified delta ray distribution functions, of simplified electron range-energy relations, and neglected angular distributions. Nevertheless it made possible the calculation of cross sections for the inactivation of enzymes and viruses, and allowed extension to tracks in nuclear emulsions and other detectors and to biological cells. It set the pattern for models of observable effects in the matter through which the ion passed. Here we outline subsequent calculations of radial dose which make use of improved knowledge of the electron emission spectrum, the electron range-energy relation, the angular distribution, and some considerations of molecular excitation, of particular interest both close to the path of the ion and the outer limits of electron penetration. These are applied to the modeling of action cross sections for the inactivation of several strains of E-coli and B. subtilis spores where extensive measurements in the "thin-down" region have been made with heavy ion beams. Such calculations serve to test the radial dose calculations at the outer limit of electron penetration. We lack data from which to test these calculations in regions close to the path of the ion aside from our earliest work on latent tracks in plastics, though it appears that the criterion then suggested for the threshold of track formation, of a minimal dose at a minimal distance (of about 20 angstroms, in plastics), remains valid.
In this work, the influence of temperature and stacking fault energy (SFE) on the strain-hardening behavior and critical resolved shear stress for twinning was investigated for three Fe–22/25/28Mn–3Al–3Si wt.% transformation- and twinning-induced plasticity (TRIP/TWIP) steels. The SFEs were calculated by two different methods, density functional theory and statistical thermodynamic modeling. The dislocation structure, observed at low levels of plastic deformation, transitions from “planar” to “wavy” dislocation glide with an increase in temperature, Mn content, and/or SFE. The change in dislocation glide mechanisms from planar to wavy reduces the strain hardening rate, in part due to fewer planar obstacles and greater cross slip activity. Additionally, the alloys exhibit a large decrease in strength and ductility with increasing temperature from 25 to 200 °C, attributed to a substantial reduction in the thermally activated component of the flow stress, predominate suppression of TRIP and TWIP, and a significant increase in the critical resolved shear stress for mechanical twinning. Interestingly, the increase in SFE with temperature had a rather minor influence on the critical resolved shear stress for mechanical twinning, and other temperature dependent factors which likely play a more dominant role are discussed.
Proper assessment of the sustainability of algal products requires an understanding of algal growth rate with regional and temporal resolution. This study developed a bulk growth model that utilizes geospatial data, reactor geometry inputs, and a maximum of six strain-specific parameters to predict algal productivity with temporal and regional resolution. The model requires local weather data to calculate the time-resolved culture temperature which is combined with a biological model based on maximum theoretical carbon fixation to predict the time-resolved algal concentration. Validation of the model against experimental data illustrates an average ac- curacy of 0.27% ± 5.32% when comparing algal concentration at harvest over 45 days of cultivation data for Galdieria sulphuraria 5587.1, and Galdieria sulphuraria Soos grown in photobioreactors (PBRs) in Mesa, AZ (33.4152° N, 111.8315° W). For open raceway ponds (ORPs), an accuracy of 0.9% ± 2.35% was achieved when comparing algal concentration at harvest over 50 days of cultivation data for Chlorella vulgaris, Desmodesmus intermedius, Nannochloropsis oceanica grown in Atlanta, GA (33.7490° N, 84.3880° W) and Mesa, AZ (33.4152° N, 111.8315° W). The validation trials spanned three seasons (winter, spring, and summer) to ensure model ac- curacy throughout the year, and concentration between harvests was modeled with hourly resolution. The validated model is used to assess the productivity at several case-study locations in the United States with Chlorella vulgaris, grown in open raceway ponds and Galdieria sulphuraria grown in vertical flat panel photo- bioreactors. In conclusion, the model has been published as an open-source tool, in an effort to increase the fidelity of future studies that rely on outdoor algae cultivation for feedstock production and allow for time resolved results for future techno-economic and life cycle assessments.
In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.
Shear assisted processing and extrusion (ShAPE) is a solid-phase processing technique that adds an additional shear force as compared with a conventional extrusion approach. Recently, ShAPE has demonstrated the capability of extruding high-performance aluminum alloy 7075 (AA7075) tubes at speeds up to 12.2 m/min without surface tearing. However, the relationship among the ShAPE processing parameters, thermomechanical conditions, contact conditions, heat generation, and microstructure evolution remains primarily empirical because an insightful understanding of the associated physics is still lacking. To help elucidate these relationships, this work proposes a thermomechanical meshfree model for the first time for ShAPE processing of AA7075 using the smoothed particle hydrodynamics (SPH) method. The meshfree model is first validated thoroughly by experimental data in terms of material flow, die face temperature, and extrusion force with various processing parameters. The validated model is then used to analyze the steady-state contact conditions and heat generation rates during ShAPE processing. Distributions of the average grain size of AA7075 being extruded are calculated using the SPH model output. The meshfree model results reveal that extrusions conducted at lower temperatures and higher strain rates yield more refined grains and possibly higher material strength, which is also consistent with the experimental observations.
Abstract Two-dimensional (2D) materials have attracted attention for quantum information science due to their ability to host single-photon emitters (SPEs). Although the properties of atomically thin materials are highly sensitive to surface modification, chemical functionalization remains unexplored in the design and control of 2D material SPEs. Here, we report a chemomechanical approach to modify SPEs in monolayer WSe 2 through the synergistic combination of localized mechanical strain and noncovalent surface functionalization with aryl diazonium chemistry. Following the deposition of an aryl oligomer adlayer, the spectrally complex defect-related emission of strained monolayer WSe 2 is simplified into spectrally isolated SPEs with high single-photon purity. Density functional theory calculations reveal energetic alignment between WSe 2 defect states and adsorbed aryl oligomer energy levels, thus providing insight into the observed chemomechanically modified quantum emission. By revealing conditions under which chemical functionalization tunes SPEs, this work broadens the parameter space for controlling quantum emission in 2D materials.
Power system simulations with long-term data tend to have large time steps varying from one second to a few minutes. However, for PV inverter semiconductors, the minimum thermal stresses cycle is with line frequency. This requires the time step of the fatigue simulation to be much smaller than the line period. This small time step results in poor simulation speed, especially for long-term simulations. This paper proposes a fast fatigue simulation for inverter semiconductors using the quasi-static time series (QSTS) simulation concept. The fatigue analysis typically focuses on the peak and valley values of a strain and neglect the transients from peaks to valleys. The proposed simulation utilizes this property of fatigue analysis and calculates the steady state of the semiconductor junction temperature only. The resulting time step of the fatigue simulation is 15 minutes, which is consistent with the solar dataset without losing accuracy.
The three-dimensional (3D) local atomic structures and crystal defects at the interfaces of heterostructures control their electronic, magnetic, optical, catalytic, and topological quantum properties but have thus far eluded any direct experimental determination. Here, we use atomic electron tomography to determine the 3D local atomic positions at the interface of a MoS 2 -WSe 2 heterojunction with picometer precision and correlate 3D atomic defects with localized vibrational properties at the epitaxial interface. We observe point defects, bond distortion, and atomic-scale ripples and measure the full 3D strain tensor at the heterointerface. By using the experimental 3D atomic coordinates as direct input to first-principles calculations, we reveal new phonon modes localized at the interface, which are corroborated by spatially resolved electron energy-loss spectroscopy. We expect that this work will pave the way for correlating structure-property relationships of a wide range of heterostructure interfaces at the single-atom level.
A comparison is made of results obtained from mathematical programming and optimality criteria procedures for the minimum-mass design of typical aircraft wing structures to satisfy prescribed flutter requirements. The mathematical programming method is based on an interior penalty function approach. A rigorous optimality criterion and an intuitive optimality criterion based on uniform strain energy density are considered. An intuitive resizing procedure is used for both optimality criterion solutions. All results are calculated using the same computer program, changing only the optimization procedure. Both high- and low-aspect-ratio wings are examined. Finite elements are used for structural modeling, and the generalized coordinates for the flutter solution are based on the natural vibration modes of the structure.