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.
In this contribution we consider the recent computation of the gauge coupling \beta β -function at four loops and the Yukawa matrix \beta β -function at three loops in the most general, renormalizable and four-dimensional quantum field theory. Furthermore, we discuss ambiguities and divergences arising in Yukawa matrix \beta β -functions.
The understanding and prediction of silicon carbide (SiC) material evolution exposed to SOL plasma conditions is of prime interest because SiC represents a promising main chamber wall plasma-facing material for next-step fusion devices (low hydrogenic diffusion, good mechanical and thermal properties under neutron irradiation). Gross and net Si erosion rates from SiC surfaces in contact with a well-diagnosed L-mode plasma in the DIII-D tokamak have been simulated and the surface concentrations of impurities have been tracked as a function of time. Coupled simulation of surface model and impurity transport demonstrates amorphization of crystalline SiC exposed to L-mode plasma due to the accumulation of displacement damages under ion irradiation. This affects the lifetime of SiC plasma facing components. Surface evolution is tightly coupled to impurity transport in the plasma and therefore needs to be integrated with impurity transport simulations to effectively predict Si erosion rates and sub-surface concentrations as a function of time. The simulation workflow couples a semi-analytical surface model to the impurity transport code GITR. The surface model is a homogeneous mixed-material model that tracks physical & chemical sputtering and reflection of impurities. Gross erosion is primarily influenced by the background plasma parameters and redeposition patterns are mainly influenced by the prompt redeposition due to the gyro-orbits of impurity ions. Although crystalline form of SiC is preferable for fusion wall applications because of resistance to neutron irradiation, this work indicates that crystalline SiC will undergo amorphization under D plasma contact with implications of higher sputtering and fuel retention. These results direct us to explore the effects of amorphization on crystalline SiC and further the physics basis of SiC usage as first wall material for fusion environments.
The design and implementation of a Concurrent Image Processing Executive (CIPE), which is intended to become the support system software for a prototype high performance science analysis workstation are described. The target machine for this software is a JPL/Caltech Mark 3fp Hypercube hosted by either a MASSCOMP 5600 or a Sun-3, Sun-4 workstation; however, the design will accommodate other concurrent machines of similar architecture, i.e., local memory, multiple-instruction-multiple-data (MIMD) machines. The CIPE system provides both a multimode user interface and an applications programmer interface, and has been designed around four loosely coupled modules: user interface, host-resident executive, hypercube-resident executive, and application functions. The loose coupling between modules allows modification of a particular module without significantly affecting the other modules in the system. In order to enhance hypercube memory utilization and to allow expansion of image processing capabilities, a specialized program management method, incremental loading, was devised. To minimize data transfer between host and hypercube, a data management method which distributes, redistributes, and tracks data set information was implemented. The data management also allows data sharing among application programs. The CIPE software architecture provides a flexible environment for scientific analysis of complex remote sensing image data, such as planetary data and imaging spectrometry, utilizing state-of-the-art concurrent computation capabilities.
The Cohen & Strand model for AP composite propellants is applied as boundary conditions, one for AP and one for binder, in solving the heat conduction equation in each to compute linear and non-linear combustion response properties for each and for the aggregate propellant. Iterations couple AP and binder through the quasi-steady flame processes. Illustrative results for linear response functions (pressure-coupled and velocity-coupled) are presented for a monomodal AP propellant varying AP size, pressure and crossflow speed, and examples of non-linear responses to arbitrary waveforms are shown. A quantitative comparison with response function data is limited to one well-characterized research formulation. Mechanistic implications are discussed, including recommendations for future modeling work.
Atomistic modeling based on Desnity Functional Theory (DFT)- 1/2 method coupled with surface Green’s function has been employed to investigate the energy band alignment results in cadmium selenium telluride (CdSeTe) surfaces. The structural and electronic properties of the bulk ternary alloy were established before cleaving the CdSeTe low index surface facets. The dependency of atomic-scale electronic features on different plane orientations was explored by studying the energy band alignment diagrams in the unreconstructed and reconstructed surfaces. While the low index CdSeTe surface geometry reconstructions replicate the CdTe surface geometries, the energy band alignment features of the unreconstructed and reconstructed CdSeTe surfaces differ from those observed in CdTe low index surfaces. The structurally relaxed reconstructed CdSeTe surface results in purely flat bands as opposed to an internal cusp feature observed in CdTe surfaces. The presence of an internal energy cusp in CdTe(111) surfaces is expected to play a key role in enhancing the hole charge transport towards the back of CdTe-only solar cell device and the absence of such feature in CdSeTe surface may explain one of the reason for the lower performance of CdSeTe-only solar cell seen experimentally.
The existence in osteoblasts of the G-protein-coupled extracellular calcium (Ca(o)(2+))-sensing receptor (CaR) that was originally cloned from parathyroid and kidney remains controversial. In our recent studies, we utilized multiple detection methods to demonstrate the expression of CaR transcripts and protein in several osteoblastic cell lines, including murine MC3T3-E1 cells. Although we and others have shown that high Ca(o)(2+) and other polycationic CaR agonists modulate the function of MC3T3-E1 cells, none of these actions has been unequivocally shown to be mediated by the CaR. Previous investigations using neurons and lens epithelial cells have shown that activation of the CaR stimulates Ca(2+)-activated K(+) channels. Because osteoblastic cells express a similar type of channel, we have examined the effects of specific "calcimimetic" CaR activators on the activity of a Ca(2+)-activated K(+) channel in MC3T3-E1 cells as a way of showing that the CaR is not only expressed in those cells but is functionally active. Patch-clamp analysis in the cell-attached mode showed that raising Ca(o)(2+) from 0.75 to 2.75 mmol/L elicited about a fourfold increase in the open state probability (P(o)) of an outward K(+) channel with a conductance of approximately 92 pS. The selective calcimimetic CaR activator, NPS R-467 (0.5 micromol/L), evoked a similar activation of the channel, while its less active stereoisomer, NPSS-467 (0.5 micromol/L), did not. Thus, the CaR is not only expressed in MC3T3-E1 cells, but is also functionally coupled to the activity of a Ca(2+)-activated K(+) channel. This receptor, therefore, could transduce local or systemic changes in Ca(o)(2+) into changes in the activity of this ion channel and related physiological processes in these and perhaps other osteoblastic cells.
The cross sections for photoionization from the 2 1S, 2 3S, 2 1P and 2 3P excited states of helium are calculated for photoelectron energies below the n = 2 threshold of He(+) using Hylleraas bound state wave functions and 1s-2s-2p close coupling final state wave functions. The resonant structures associated with the lowest-lying 1S, 1P, 3P, and 1D autoionizing states of helium are found to be characterized by large values of the line profile parameter q. The cross sections and the photoelectron angular distribution asymmetry parameters for the P-states are calculated for various polarization states of the target atom and the incident photon. Experiments which would lead to the separate determinations of the S- and D- wave partial photoionization cross sections are discussed.
The design and implementation of a Concurrent Image Processing Executive (CIPE), which is intended to become the support system software for a prototype high performance science analysis workstation are discussed. The target machine for this software is a JPL/Caltech Mark IIIfp Hypercube hosted by either a MASSCOMP 5600 or a Sun-3, Sun-4 workstation; however, the design will accommodate other concurrent machines of similar architecture, i.e., local memory, multiple-instruction-multiple-data (MIMD) machines. The CIPE system provides both a multimode user interface and an applications programmer interface, and has been designed around four loosely coupled modules; (1) user interface, (2) host-resident executive, (3) hypercube-resident executive, and (4) application functions. The loose coupling between modules allows modification of a particular module without significantly affecting the other modules in the system. In order to enhance hypercube memory utilization and to allow expansion of image processing capabilities, a specialized program management method, incremental loading, was devised. To minimize data transfer between host and hypercube a data management method which distributes, redistributes, and tracks data set information was implemented.
Understanding structure-property relationships is foundational to numerous modern chemistries, such as proton-coupled electron transfer (PCET). However, an experimentally measured property is the result of the behavior from an ensemble of molecules. Neglecting ensemble effects, especially under complex chemical environments, may obfuscate these relationships and lead to discrepancies between theory and experiment. In this work, we demonstrate the impact of configurational entropy and local chemical environments on hydroxide bond dissociation free energies [BDFE- (O−H)] for a set of polyoxovanadate nanoclusters, at ambient conditions. The O−H bond strengths are investigated via density functional theory (DFT) coupled with statistical thermodynamic analysis and bilinear modeling, and compared with previous experimental results on the same systems, namely electrochemical solutions of: [V 6 O 13−x (OH) x (TRIOL R ) 2 ] −2 (x = 2, 4, 6; R = NO 2 , Me) and [V 6 O 11−x (OMe) 2 (OH) x (TRIOL NO 2 ) 2 ] −2 (x = 2, 4). Interestingly, we find that ensemble effects, even at room temperature, can account for a significant portion of the BDFE(O−H) trend with the degree of reduction via H atom binding, which cannot be fully captured by single-structure, static DFT calculations. Moreover, we find that the ensemble effects may be replicated statistically, requiring only enumeration of energetically accessible H-binding sites. With the ensemble effects resolved, we present a simple bilinear model to reconcile remaining biases between experiment and ensemble-informed theory, which corelate with clusterspecific electronic environment differences. The bilinear model achieves outstanding accuracy vs experiments with a root-mean squared error of 0.4 kcal/mol. Finally, based on the physicochemical characteristics of hydrogen interaction with polyoxometalates, we present a simple methodology that captures the BDFE(O−H) trend while dramatically reducing required DFT calculations by 98% and achieving accuracy within 1 kcal/mol. Overall, this work elucidates the roles and structural origins of configurational entropy and chemical effects on polyoxometalate hydroxide bond energies, with potential applicability to various atomically precise metal oxide systems. Importantly, it introduces models for rapid and highly accurate property calculations in connection with experiments.
Novel heterostructures created by coupling one-dimensional semiconductor nanowires with a superconducting thin film show great potential toward next-generation quantum computing. Here, by growing high-crystalline SiGe nanowires on a NbTiN thin film, the resulting heterostructure exhibits Ohmic characteristics as well as a shift of the superconducting transition temperature (T c ). The structure was characterized at atomic resolution showing a sharp SiGe/NbTiN interface without atomic interdiffusion. Lattice spacing, as calculated from large-area x-ray diffraction experiments, suggests a potential preferred d-spacing matching between (200) NbTiN and (110) SiGe grains. The observed out-of-plane compressive strain within the NbTiN films coupled with SiGe nanowires explains the downward shift of the superconductivity behavior. The presented results post scientific insights toward functional heterostructures by coupling multi-dimensional materials, which could enable tunable superconductivity that benefits the quantum science applications.
In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.
The paper presents theoretical models developed to account for the heterogeneity of composite propellants in expressing the pressure-coupled combustion response function. It is noted that the model of Lengelle and Williams (1968) furnishes a viable basis to explain the effects of heterogeneity.
In chloride containing environments, two metals in physical contact can undergo galvanic corrosion limiting the lifetime of components. Being able to accurately predict galvanic corrosion damage distributions over time has not been widely presented in literature. Therefore, Finite Element Method (FEM) corrosion models were experimentally validated for two galvanic couples as a function of governing equations, environment, anode material, anode:cathode ratio, and time. Carbon steel/stainless steel (CS/SS) and zinc/stainless steel (Zn/SS) galvanic couples were exposed to NaCl solutions at room temperatures for up to 14 days. For the galvanic couples and environments, the Laplace equation with variable conductivity and reactions was sufficient to model the experimental corrosion damage. Increasing the chloride concentration for the CS/SS galvanic couple, regardless of the anode:cathode ratio, decreased the observed and modeled corrosion damage. Decreasing the anode:cathode ratio (i.e., increasing the cathode length), increased the experimentally observed and modeled corrosion damage. For small anode:cathode ratios, the governing equations deviate over long time periods with the Nernst-Plank equation being conservative. For the Zn/SS galvanic couple, the cathode length controls the dominant cathodic reduction reaction. For small cathode lengths, the hydrogen evolution reaction is dominant. For larger cathode lengths, the oxygen reduction reaction is dominant. The results are discussed with regard to the influence of solution chemistry, ohmic drop, and governing reactions. Overall, validated FEM models were presented, and the resultant models, physics, and mechanisms can be applied to other corrosion scenarios with confidence.
The NEQAIR line-by-line radiation code has been incorporated into the DPLR Navier-Stokes flow solver such that the NEQAIR subroutines are now callable functions of DPLR. The coupled DPLR-NEQAIR code was applied to compute the convective and radiative heating rates over high-mass Mars entry vehicles. Two vehicle geometries were considered - a 15 m diameter 70-degree sphere cone configuration and a slender, mid-L/D vehicle with a diameter of 5 m called an Ellipsled. The entry masses ranged from 100 to 165 metric tons. Solutions were generated for entry velocities ranging from 6.5 to 9.1 km/s. The coupled fluids-radiation solutions were performed at the peak heating location along trajectories generated by the Traj trajectory analysis code. The impact of fluids-radiation coupling is a function of the level of radiative heating and the freestream density and velocity. For the high-mass Mars vehicles examined in this study, coupling effects were greatest for entry velocities above 8.5 km/s where the surface radiative heating was reduced by up 17%. Generally speaking, the Ellipsled geometry experiences a lower peak radiative heating rate but a higher peak turbulent convective heating rate than the MSL-based vehicle.
The NEQAIR line-by-line radiation code has been incorporated into the DPLR Navier-Stokes flow solver such that the NEQAIR subroutines are now callable functions of DPLR. The coupled DPLR-NEQAIR code was applied to compute the convective and radiative heating rates over high-mass Mars entry vehicles. Two vehicle geometries were considered - a 15 m diameter 70-degree sphere cone configuration and a slender, mid-L/D vehicle with a diameter of 5 m called an Ellipsled. The entry masses ranged from 100 to 165 metric tons. Solutions were generated for entry velocities ranging from 6.5 to 9.1 km/s. The coupled fluids-radiation solutions were performed at the peak heating location along trajectories generated by the Traj trajectory analysis code. The impact of fluids-radiation coupling is a function of the level of radiative heating and the freestream density and velocity. For the high-mass Mars vehicles examined in this study, coupling effects were greatest for entry velocities above 8.5 km/s where the surface radiative heating was reduced by up 17%. Generally speaking, the Ellipsled geometry experiences a lower peak radiative heating rate but a higher peak turbulent convective heating rate than the MSL-based vehicle.
Bayesian optimization (BO) is an efficient and flexible global optimization framework that is applicable to a very wide range of engineering applications. To leverage the capability of the classical BO, many extensions, including multi-objective, multi-fidelity, parallelization, and latent-variable modeling, have been proposed to address the limitations of the classical BO framework. In this work, we propose a novel multi-objective BO formalism, called srMO-BO-3GP, to solve multi-objective optimization problems in a sequential setting. Three different Gaussian processes (GPs) are stacked together, where each of the GPs is assigned with a different task. The first GP is used to approximate a single-objective computed from the multi-objective definition, the second GP is used to learn the unknown constraints, and the third one is used to learn the uncertain Pareto frontier. At each iteration, a multi-objective augmented Tchebycheff function is adopted to convert multi-objective to single-objective, where the regularization with a regularized ridge term is also introduced to smooth the single-objective function. Finally, we couple the third GP along with the classical BO framework to explore the convergence and diversity of the Pareto frontier by the acquisition function for exploitation and exploration. The proposed framework is demonstrated using several numerical benchmark functions, as well as a thermomechanical finite element model for flip-chip package design optimization.
The definition of some modern high demanding space systems requires a different approach to system definition and design from that adopted for traditional missions. System functionality is strongly coupled to the operational analysis, aimed at characterizing the dynamic interactions of the flight element with its surrounding environment and its ground control segment. Unambiguous functional, operational and performance requirements are to be defined for the system, thus improving also the successive development stages. This paper proposes a Petri Nets based methodology and two related prototype applications (to ARISTOTELES orbit control and to Hermes telemetry generation) for the operational analysis of space systems through the dynamic modeling of their functions and a related computer aided environment (ISIDE) able to make the dynamic model work, thus enabling an early validation of the system functional representation, and to provide a structured system requirements data base, which is the shared knowledge base interconnecting static and dynamic applications, fully traceable with the models and interfaceable with the external world.
BiVO 4 is a promising photocatalyst for efficient water oxidation, with surface reactivity determined by the structure of active catalytic sites. Surface oxidation in the presence of oxygen vacancies induces electron localization, suggesting an atomistic route to improve the charge transfer efficiency within the catalytic cycle. In this paper, we study the effect of oxygen vacancies on the electronic and optical properties at BiVO 4 surfaces upon water oxidation. We use density functional theory and many-body perturbation theory to explore the change in the electronic and quasiparticle energy levels and to evaluate the electron-hole coupling as a function of the underlying structure. Further we show that while the presence of defects alters the atomic structure and largely modifies the wave-function nature, leading to defect-localized states at the quasiparticle gap region, the optical excitations remain largely unchanged due to the substantial hybridization of defect and nondefect electron-hole transitions. Our findings suggest that defect-induced surface oxidation supports improved electron transport, both through bound and tunable electronic states and via a mixed nature of the optical transitions, expected to reduce electron-hole defect trapping.