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.
The most recent (2015) results for 0νββ nuclear matrix elements in the interacting boson model (IBM-2) with light and heavy neutrino exchange, including R0νECEC, are given for all nuclei of interest from {sup 48}Ca to {sup 238}U. Predictions for half-lives and limits for average neutrino mass are also made. Possible additional scenarios, such as Majoron emission, is also discussed.
A model for predicting the direction of crack growth in resin matrix composite materials is presented. The model is based upon the normal stress and the anisotropic tensile strength on arbitrary planes about the tip of a crack. The model is incorporated into a finite element solution which predicts energy release rates in addition to the direction of crack growth. Results are obtained for cracks in unidirectional off-axis slotted tensile coupons. Comparisons are made with other theories and experimental results.
A method for adaptive identification of reduced-order models for continuous stable SISO and MIMO plants is presented. The method recursively finds a model whose transfer function (matrix) matches that of the plant on a set of frequencies chosen by the designer. The algorithm utilizes the Moving Discrete Fourier Transform (MDFT) to continuously monitor the frequency-domain profile of the system input and output signals. The MDFT is an efficient method of monitoring discrete points in the frequency domain of an evolving function of time. The model parameters are estimated from MDFT data using standard recursive parameter estimation techniques. The algorithm has been shown in simulations to be quite robust to additive noise in the inputs and outputs. A significant advantage of the method is that it enables a type of on-line model validation. This is accomplished by simultaneously identifying a number of models and comparing each with the plant in the frequency domain. Simulations of the method applied to an 8th-order SISO plant and a 10-state 2-input 2-output plant are presented. An example of on-line model validation applied to the SISO plant is also presented.
An approach based on the Constant Gain Extended Kalman Filter (CGEKF) technique is investigated for the in-flight estimation of non-measurable performance parameters of aircraft engines. Performance parameters, such as thrust and stall margins, provide crucial information for operating an aircraft engine in a safe and efficient manner, but they cannot be directly measured during flight. A technique to accurately estimate these parameters is, therefore, essential for further enhancement of engine operation. In this paper, a CGEKF is developed by combining an on-board engine model and a single Kalman gain matrix. In order to make the on-board engine model adaptive to the real engine s performance variations due to degradation or anomalies, the CGEKF is designed with the ability to adjust its performance through the adjustment of artificial parameters called tuning parameters. With this design approach, the CGEKF can maintain accurate estimation performance when it is applied to aircraft engines at offnominal conditions. The performance of the CGEKF is evaluated in a simulation environment using numerous component degradation and fault scenarios at multiple operating conditions.
It is often advantageous to account for the microstructure of the material directly using multiscale modeling. For computational tractability, an idealized repeating unit cell (RUC) is used to capture all of the pertinent features of the microstructure. Typically, the RUC is dimensionless and depends only on the relative volume fractions of the different phases in the material. This works well for non-linear and inelastic behavior exhibiting a positive-definite constitutive response. Although, once the material exhibits strain softening, or localization, a mesh objective failure theories, such as smeared fracture theories, nodal and element enrichment theories (XFEM), cohesive elements or virtual crack closure technique (VCCT), can be utilized at the microscale, but the dimensions of the RUC must then be defined. One major challenge in multiscale progressive damage modeling is relating the characteristic lengths across the scales in order to preserve the energy that is dissipated via localization at the microscale. If there is no effort to relate the size of the macroscale element to the microscale RUC, then the energy that is dissipated will remain mesh dependent at the macroscale, even if it is regularized at the microscale. Here, a technique for mapping characteristic lengths across the scales is proposed. The RUC will be modeled using the generalized method of cells (GMC) micromechanics theory, and local failure in the matrix constituent subcells will be modeled using the crack band theory. The subcell characteristic lengths used in the crack band calculations will be mapped to the macroscale finite element in order to regularize the local energy in a manner consistent with the global length scale. Examples will be provided with and without the regularization, and they will be compared to a baseline case where the size and shape of the element and RUC are coincident (ensuring energy is preserved across the scales).
Missions to Mars will differ from previous human spaceflight missions in that the crew of astronauts will be required to operate in an Earth-independent manner due to the long communication delays. Without a systematic, repeatable process to determine the number and composition of crew necessary to successfully accomplish these missions, NASA increases the risk that crew sizes may be too small to meet primary mission objectives under nominal conditions and, more consequentially, that crewmembers may not have the expertise needed to successfully respond to unforeseen failures without the real-time expertise of the Mission Control Central (MCC) team NASA currently relies upon. The NASA Engineering and Safety Center (NESC) developed a methodology for assessing the trade space of factors that affect the number of crew for Mars missions based on human performance and expertise modeling. This methodology includes the consideration of results from three human performance models developed using the Improved Performance Research and Integration Tool (IMPRINT) modeling platform as well as a custom-built model on expertise trained within the crew based on NASA’s crew qualification and responsibilities matrix (CQRM). We present example model results and discuss implications on trade space analyses for Mars mission crew size.
The physical basis for how macromolecules regulate the onset of mineral formation in calcifying tissues is not well established. A popular conceptual model assumes the organic matrix provides a stereochemical match during cooperative organization of solute ions. In contrast, another uses simple binding assays to identify good promoters of nucleation. Here, we reconcile these two views and provide a mechanistic explanation for template-directed nucleation by correlating heterogeneous nucleation barriers with crystal–substrate-binding free energies. We first measure the kinetics of calcite nucleation onto model substrates that present different functional group chemistries (carboxyl, thiol, phosphate, and hydroxyl) and conformations (C11 and C16 chain lengths). We find rates are substrate-specific and obey predictions of classical nucleation theory at supersaturations that extend above the solubility of amorphous calcium carbonate. Analysis of the kinetic data shows the thermodynamic barrier to nucleation is reduced by minimizing the interfacial free energy of the system, ?. We then use dynamic force spectroscopy to independently measure calcite–substrate-binding free energies, ?Gb. Moreover, we show that within the classical theory of nucleation, ? and ?Gb should be linearly related. The results bear out this prediction and demonstrate that low-energy barriers to nucleation correlate with strong crystal–substrate binding. This relationship is general to all functional group chemistries and conformations. These findings provide a physical model that reconciles the long-standing concept of templated nucleation through stereochemical matching with the conventional wisdom that good binders are good nucleators. The alternative perspectives become internally consistent when viewed through the lens of crystal–substrate binding.
Super-critical carbon dioxide (S-CO2) is a promising thermodynamic cycle for advanced nuclear reactors and solar energy conversion applications. Dynamic control of the proposed recompression S-CO2 cycle is accomplished with input from resistance temperature detector (RTD) measurements of the process fluid. One of the challenges in practical implementation of S-CO2 cycle is high corrosion rate of component and sensor materials. In this paper, we develop a mathematical model of RTD sensing using eigendecomposition model of radial heat transfer in a layered long cylinder. We show that the value of RTD time constant primarily depends on the rate of heat transfer from the fluid to the outer wall of RTD. We also show that for typical material properties, RTD time constant can be calculated as the sum of reciprocal eigen-values of the heat transfer matrix. Using the computational model and a set of RTD and CO2 fluid thermophysical parameter values, we calculate the value of time constant of thermowell-mounted RTD sensor at the hot side of the precooler in the S-CO2 cycle. The eigendecomposition model of RTD will be used in future studies to model sensor degradation and its impact on control of S-CO2. (C) 2016 Elsevier B.V. All rights reserved.
We study the operation of an 8.5 μm quantum cascade laser based on GaInAs/AlInAs lattice matched to InP using three different simulation models based on density matrix (DM) and non-equilibrium Green's function (NEGF) formulations. The latter advanced scheme serves as a validation for the simpler DM schemes and, at the same time, provides additional insight, such as the temperatures of the sub-band carrier distributions. We find that for the particular quantum cascade laser studied here, the behavior is well described by simple quantum mechanical estimates based on Fermi's golden rule. As a consequence, the DM model, which includes second order currents, agrees well with the NEGF results. Both these simulations are in accordance with previously reported data and a second regrown device.
Stable isotope based metabolic flux analysis is currently the unique methodology that allows the experimental study of the integrated responses of metabolic networks. This method primarily relies on isotope labeling and modeling, which could be a challenge in both experimental and computational biology. Specifically, the algorithm implementation for isotope simulation is a critical step, limiting extensive usage of this powerful approach. In this work, we introduce EMUlator a Python-based isotope simulator which is developed on Elementary Metabolite Unit (EMU) algorithm, an efficient and powerful algorithm for isotope modeling. We propose a novel adjacency matrix method to implement EMU modeling and exemplify it stepwise. This method is intuitively straightforward and can be conveniently mastered for various customized purposes. We apply this arithmetic pipeline to understand the phosphoketolase flux in the metabolic network of an industrial microbe Clostridium acetobutylicum . The resulting design enables a high-throughput and non-invasive approach for estimating phosphoketolase flux in vivo . Our computational insights allow the systematic design and prediction of isotope-based metabolic models and yield a comprehensive understanding of their limitations and potentials.
A method for constructing the generalized system-level admittance matrix for use with a Newton-Raphson power flow is presented. The network modeling technique presented does not use the standard pi-equivalent models, which assume a lossless return path, for the transmission line and transformer. If the return path cannot be assumed lossless, then the standard algorithms for constructing the system admittance matrix cannot be used. The method presented here uses concepts from linear graph theory to combine network modules to form the system-level admittance matrix. The modeling technique is presented, and the resulting matrix is used with a standard Newton-Raphson power flow to calculate all system voltages and current (power) flows.
Poor processability of fiber reinforced high performance polyimide thermoplastic resin composites is a well recognized issue which, in many cases, prohibits the fabrication of composite parts with satisfactorily consolidated quality. Without modifying the resin matrix chemistry, improved compression modeling procedures were proposed and investigated with the AS-4/LaRC-TPI 1500 High Flow Grade (HFG) prepreg system. Composite panels with excellent C-scans can be consistently molded by this method under 700 F and a consolidation pressure as low as 100 psi. A mechanism for the consolidation of the composite under this improved molding technique is discussed. This mechanism reveals that a certain degree of matrix shear and tow filament slippage and nesting between plies occur during consolidation, which leads to a reduction of the consolidating pressure necessary to offset the otherwise intimate inter fiber-fiber contact and consequently achieves a better consolidation quality. Outstanding short beam shear strength and flexural strength were obtained from the molded panels. A prolonged consolidation step under low pressure, i.e., 100 psi at 700 F for 75 minutes, was found to significantly enhance the composite mechanical properties.
We show evidence of an antioxidant mechanism for vitamin E which correlates strongly with its physical location in a model lipid bilayer. This data addresses the overlooked problem of the physical distance between the vitamin's reducing hydrogen and lipid acyl chain radicals. Our combined data from neutron diffraction, NMR and UV spectroscopy experiments, all suggest that reduction of reactive oxygen species and lipid radicals occurs specifically at the membrane's hydrophobic-hydrophilic interface. The latter is possible when the acyl chain snorkels to the interface from the hydrocarbon matrix. Moreover, not all model lipids are equal in this regard, as indicated by the small differences in vitamin's location. The present result is a clear example of the importance of lipid diversity in controlling the dynamic structural properties of biological membranes. Importantly, our results suggest that measurements of aToc oxidation kinetics, and its products, should be revisited by taking into consideration the physical properties of the membrane in which the vitamin resides.
We investigate the effect of surface disorder on the chiral surface states of a three-dimensional quantum Hall system. Utilizing a transfer-matrix method, we find that the localization length of the surface state along the magnetic field decreases with the surface disorder strength in the weak disorder regime but increases anomalously in the strong disorder regime. In the strong disorder regime, the surface states mainly locate at the first inward layer to avoid the strong disorder in the outmost layer. The anomalous increase of the localization length can be explained by an effective model, which maps the strong disorder on the surface layer to the weak disorder on the first inward layer. Our work demonstrates that surface disorder can be an effective way to control the transport behavior of the surface states along the magnetic field. We also investigate the effect of surface disorder on the full distribution of conductances P(g) of the surface states in the quasi-one-dimensional (1D) regime for various surface disorder strengths. In particular, we find that P(g) is Gaussian in the quasi-1D metal regime and log-normal in the quasi-1D insulator regime. In the crossover regime, P(g) exhibits highly nontrivial forms, whose shapes coincide with the results obtained from the Dorokhov-Mello-Pereyra-Kumar equation of a weakly disordered quasi-1D wire in the absence of time-reversal symmetry. Our results suggest that P(g) is fully determined by the average conductance, independent of the details of the system, in agreement with the single-parameter scaling hypothesis.
We present a new linearization of T-Matrix and Mie computations for light scattering by non-spherical and spherical particles, respectively. In addition to the usual extinction and scattering cross-sections and the scattering matrix outputs, the linearized models will generate analytical derivatives of these optical properties with respect to the real and imaginary parts of the particle refractive index, and (for non-spherical scatterers) with respect to the ''shape'' parameter (the spheroid aspect ratio, cylinder diameter/height ratio, Chebyshev particle deformation factor). These derivatives are based on the essential linearity of Maxwell's theory. Analytical derivatives are also available for polydisperse particle size distribution parameters such as the mode radius. The T-matrix formulation is based on the NASA Goddard Institute for Space Studies FORTRAN 77 code developed in the 1990s. The linearized scattering codes presented here are in FORTRAN 90 and will be made publicly available.
A comprehensive, point-information-to-continuum-level analysis framework is presented in this paper to accurately characterize the behavior of carbon nanotube (CNT)-enhanced composite materials. Molecular dynamics (MD) simulations are performed to study sub-nanoscale interactions of the CNT with the polymeric phase of the nanocomposite. The effect of cross-linking between the epoxy resin and the hardener on the mechanical properties of the polymer is investigated; furthermore, the effect of CNT weight fraction on the probability distribution of polymer cross-linking degree is also studied through stochastic models. The stochastic distributions obtained from MD simulations provide a basis to simulate local variations in the matrix properties in the continuum model at the microscale. The inclusion of an atomistically informed elastic-plastic model at the microscale reveals a significant deviation of the mechanical properties from those obtained based on classical homogenization approaches. Microstructural variability arising from heterogeneous cross-linking degree in the polymer phase and variations in fiber geometry and spacing is also found to cause deviations in the mechanical response when compared to the assumption of a perfectly ordered fiber-matrix microstructure.
The complexity of photocathode designs and detector materials, and the need to model their performance for short pulse durations, response to high frequency photons, presence of coatings and/or thinness of the absorptive layer, necessitates modifications to three-step and moments models of photoemission that are used in simulation codes. Here, methods to include input from computationally intensive approaches, such as density functional theory to model optical properties and transfer matrix approaches to treat emission from the surface or transport past coatings, by means of parametric models are demonstrated. First, a technique to accurately represent optical behavior so as to model reflectivity and penetration depth is given. Second, modifications to bulk models arising from the usage of thin film architectures, and a means to rapidly calculate them, is provided. Third, a parameterization to model the impact of wells associated with coatings and surface layers on the transmission probably is given. In all cases, the methods are computationally efficient and designed to allow for including input from numerically intensive approaches that would otherwise be unavailable for simulations.
A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The material stiffness and strength values required for the constitutive model are determined based on coupon level tests on the braided composite. Simulations of quasi-static coupon tests of a representative braided composite are conducted. Varying the strength values that are input to the material model is found to have a significant influence on the effective material response predicted by the finite element analysis, sometimes in ways that at first glance appear non-intuitive. A parametric study involving the input strength parameters provides guidance on how the analysis model can be improved.