Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled processes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

Identification of hot water end-use process of electric water heaters from energy measurements

This paper presents an algorithm for the identification of parameters for a stochastic hot water end-use process that drives a homogeneous population of thermostatically controlled electric water heaters (EWH). Usually, only metered interval consumption data (kWh) is collected and the hot water end-use process is unobservable to utility and aggregators. However, the availability of EWHs for demand response (DR) is closely coupled with the hot water end-use process. In this context, the hot water end-use process is modeled as a two-state Markov chain (Use / No use), which causes the thermostatic ON-OFF switching process to behave as a Markov renewal process (MRP). A set of first passage-time problems is developed to obtain the moments of the transition probability densities of the MRP. These problems are addressed by establishing a system of coupled partial differential equations characterizing the temperature evolution of the EWH population. A key quantity in the methodology for estimating the parameters is the total time an EWH is ON within a period of interest. It is referred to as the total busy time. Total busy time in this approach is a random variable for which analytical expressions of the moments are developed as a function of the metered window length. The latter expressions become the basis of a hot water demand model identification algorithm which is validated using agent-based simulations of EWHs.

42 ENGINEERING↗

Leptonic scalars and collider signatures in a UV-complete model

We study the non-standard interactions of neutrinos with light leptonic scalars ( ϕ ) in a global ( B – L )-conserved ultraviolet (UV)-complete model. The model utilizes Type-II seesaw motivated neutrino interactions with an SU(2) L -triplet scalar, along with an additional singlet in the scalar sector. This UV-completion leads to an enriched spectrum and consequently new observable signatures. We examine the low-energy lepton flavor violation constraints, as well as the perturbativity and unitarity constraints on the model parameters. Then we lay out a search strategy for the unique signature of the model resulting from the leptonic scalars at the hadron colliders via the processes H ±± → W ± W ± ϕ and H ± → W ± ϕ for both small and large leptonic Yukawa coupling cases. We find that via these associated production processes at the HL-LHC, the prospects of doubly-charged scalar H ±± can reach up to 800 (500) GeV and 1.1 (0.8) TeV at the 2 σ (5 σ ) significance for small and large Yukawa couplings, respectively. A future 100 TeV hadron collider will further increase the mass reaches up to 3.8 (2.6) TeV and 4 (2.7) TeV, at the 2 σ (5 σ ) significance, respectively. We also demonstrate that the mass of ϕ can be determined at about 10% accuracy at the LHC for the large Yukawa coupling case even though it escapes as missing energy from the detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A numerical model for coupling between atomization and spray dynamics in liquid rocket thrust chambers

This paper describes a novel method of coupling the atomization and spray combustion processes encountered in coaxial injection elements of liquid rocket engine thrust chambers. This method is based on the Jet-Embedding technique in which the liquid jet core equations and the gas phase equations are solved separately. The liquid and gas phase solutions, however, are coupled through the boundary conditions at the interface between the phases. The computational grid for the gas phase calculations are adapted to the shape of the liquid jet core. The axial variation of droplet sizes are calculated using a stability analysis appropriate for the atomization regime of liquid jet breakup. The predictions of this method have been validated with experimental data on low speed water jets. Using this method, calculations are performed for the SSME fuel preburner single injector flow field. The results obtained are in good agreement with the predictions of the volume-of-fluid method.

Giridharan, M. G.↗

Process-based modeling of soil nitrous oxide emissions from United States corn fields under different management and climate scenarios coupled with evaluation using regional estimates

Direct emissions of soil nitrous oxide during a growing season (N 2 O gs ) can be quantified with process-based models considering interactions between management, climate, and soil moisture when key data are available. We used an adapted “parameterized CENTURY/DAYCENT-model” ( p CENTURY) calibrated with crop growth and soil organic matter decay coefficients at the county-level for the estimation of N 2 O gs in the United States Corn Belt. Model estimated N 2 O-emissions from corn-based biofuels scenarios considering crop rotation, fertilizer inputs, tillage, and weather were compared against meta-summary of field observations from 55 studies. Both model and meta-summary ranked N 2 O gs -emissions to be corn > wheat > soybean phase while model likely underestimated cover crop N 2 O gs -emissions. The N 2 O gs -emissions and the associated emission factors (EFs) were modeled and summarized to be greater after anhydrous ammonia than urea application and from conventional tilled than non-tilled fields. Modeled and observed N 2 O gs -emissions after organic and inorganic fertilizer amendment did not differ due to high variability associated with the treatments. However, the organic fertilizer associated EFs were greater according to meta-summary data because of N input rates. Regionalized weather scenarios indicate hotspots for N 2 O gs -emissions can occur where crop N uptake is limited during dry years and in eastern states also during normal or wet seasons. The p CENTURY-derived N 2 O gs EFs (0.91 ± 0.19%) for counties investigated were only slightly lower than literature (1.07 ± 0.57%) or Tier-1 (1%) values. Our preliminary evaluation of regional soil moisture estimates showed reasonable agreement between monthly soil moisture estimates and the North American Soil Moisture Dataset during the growing season, but overestimation of soil moisture in winter-spring can influence the estimates of annual N 2 O emissions so future work is needed to calibrate soil moisture-associated model parameters. Our work provided scenario-based estimates of climate and management impacts on soil N 2 O gs -emissions together with valuable spatial insights into EFs that will be improved by more accurate information of fertilizer inputs and more temporally refined model evaluation.

54 ENVIRONMENTAL SCIENCES↗

High-Performance Piezoelectric MEMS Manufacturing & Application on Micro Power Generators (CRADA Final Report)

This project aimed to develop and characterize a micro-fabrication technology for wafer-scale heterogeneous integration of bulk piezoelectric materials on silicon substrates that could enable manufacturing of high-performance micro transducers. Typical technical challenges for piezoelectric microfabrication technologies are high-temperature processing, non-uniform film quality, and low process repeatability, which will be addressed during this project. Progress was made on preliminary research towards development of a technology platform that will enable a new micro-fabrication process to provide the highest electromechanical coupling among any other thin/thick-film deposition methods available today. This work has tested and characterized samples that we hope will further our goal of refining and enabling existing micro transducers to achieve unprecedentedly high performance. This project focused solely on the characterization of the proposed microfabrication process technology, and it did not consist of design or manufacturing of new transducers, such as sensors, actuators or energy harvesters based on this fabrication process. Future work will be aimed at the research and development of process steps specifically to enable integration of on-chip electronics with piezoelectric materials.

36 MATERIALS SCIENCE↗

Partitioning and packing mathematical simulation models for calculation on parallel computers

The development of multiprocessor simulations from a serial set of ordinary differential equations describing a physical system is described. Degrees of parallelism (i.e., coupling between the equations) and their impact on parallel processing are discussed. The problem of identifying computational parallelism within sets of closely coupled equations that require the exchange of current values of variables is described. A technique is presented for identifying this parallelism and for partitioning the equations for parallel solution on a multiprocessor. An algorithm which packs the equations into a minimum number of processors is also described. The results of the packing algorithm when applied to a turbojet engine model are presented in terms of processor utilization.

Arpasi, D. J.↗

Tailored metal matrix composites for high-temperature performance

A multi-objective tailoring methodology is presented to maximize stiffness and load carrying capacity of a metal matrix cross-ply laminated at elevated temperatures. The fabrication process and fiber volume ratio are used as the design variables. A unique feature is the concurrent effects from fabrication, residual stresses, material nonlinearity, and thermo-mechanical loading on the laminate properties at the post-fabrication phase. For a (0/90)(sub s) graphite/copper laminate, strong coupling was observed between the fabrication process, laminate characteristics, and thermo-mechanical loading. The multi-objective tailoring was found to be more effective than single objective tailoring. Results indicate the potential to increase laminate stiffness and load carrying capacity by controlling the critical parameters of the fabrication process and the laminate.

Morel, M. R.↗

Tailored metal matrix laminates for high-temperature performance

A multi-objective tailoring methodology is presented to maximize stiffness and load carrying capacity of a metal matrix cross-ply laminated at elevated temperatures. The fabrication process and fiber volume ratio are used as the design variables. A unique feature is the concurrent effects from fabrication, residual stresses, material nonlinearity, and thermo-mechanical loading on the laminate properties at the post-fabrication phase. For a (0.90) (sub s) graphite/copper laminate, strong coupling was observed between the fabrication process, laminate characteristics, and thermo-mechanical loading. The multi-objective tailoring was found to be more effective than single objective tailoring. Results indicate the potential to increase laminate stiffness and load carrying capacity by controlling the critical parameters of the fabrication process and the laminate.

Morei, Michael R.↗

Multi-channel front-end ASIC for a 3D position-sensitive detector

Arrays of 3D position-sensitive detectors (3DPSD), operating at room temperature and using cadmium zinc telluride (CZT) and thallium bromide (TIBr) sensors, are suitable for gamma-ray spectrometry in many applications. One detector configuration, the 3D position-sensitive Virtual Frisch-Grid detector (VFG), is particularly advantageous for integrating into large area arrays. The signals generated inside each detector of the array are captured with the anode, cathode and four pads that enable the reconstruction of the position and energy of the ionizing interaction by measurements of amplitude and timing of the signals. For these applications, a low-noise front-end ASIC has been developed, capable of processing bipolar signals (needed because of AC-coupling of certain electrodes). The ASIC can be coupled to an ADC in order to form a compound “waveform digitizer” capable of post-processing the analog signals and determining amplitude and timing information. This paper describes a 32-channel front-end ASIC that is suitable for reading out a 3 × 3 or 4 × 4 element matrix in the VFG configuration. Each channel is composed of a low-noise charge amplifier with an adaptive continuous reset feedback circuit suitable for both positive and negative charge, a first order shaper and a single-to-differential converter output stage. Voltage and current references are all internally generated by 10-bit DACs and the chip is fully controllable with the I 2 C communication protocol. The readout channel response has been verified using the implemented injection circuit. Linear behavior up to ~75 ke ± with the gain of ~80 mV/fC, and up to ~200 ke ± with the gain of ~30 mV/fC was demonstrated. In conclusion, the first test result waveforms using a 137 Cs radioactive source on a 5 × 5 × 12 mm 3 TIBr crystal are reported.

47 OTHER INSTRUMENTATION↗

An investigation of voids formation mechanisms and their effects on freeze and thaw processes of lithium and lithium fluoride

The mechanisms of void formation during the cooldown and freezing of lithium coolant within the primary loop of SP-100 type systems are investigated. These mechanisms are: (1) homogeneous nucleation; (2) heterogeneous nucleation; (3) normal segregation of helium gas dissolved in liquid lithium; and (4) shrinkage of lithium during freezing. To evaluate the void formation potential due to segregation, a numerical scheme that couples the freezing and mass diffusion processes in both the solid and liquid regions is developed. The results indicated that the formation of He bubbles is unlikely by either homogeneous or heterogeneous nucleation during the cooldown process. However, homogeneous nucleation of He bubbles following the segregation of dissolved He in liquid lithium ahead of the solid-liquid interface is likely to occur. Results also show that total volume of He void is insignificant when compared to that of shrinkage voids. In viewing this, the subsequent research focuses on the effects of shrinkage void forming during freezing of lithium on subsequent thaw processes are investigated using a numerical scheme that is based on a single (solid/liquid) cell approach. The cases of lithium-fluoride are also investigated to show the effect of larger volume shrinkage upon freezing on the freeze and thaw processes. Results show that a void forming at the wall appreciably reduces the solid-liquid interface velocity, during both freeze and thaw, and causes a substantial rise in the wall temperature during thaw. However, in the case of Li, the maximum wall temperature was much lower than the melting temperature of PWC-11, which is used as the structure material in the SP-100 system. Hence, it is included that a formation of hot spots is unlikely during the startup or restart of the SP-100 system.

El-Genk, Mohamed S.↗

Coupled Reactor and Engine Nuclear Thermal Propulsion Modeling Methodology

The design and development process of a Nuclear Thermal Propulsion (NTP) system requires extensive multiphysics modeling to couple the neutron physics and thermal feedback effects to determine the reactor’s power shape. Propulsion system performance codes utilize this power shape to determine NTP key performance parameters. While the power shape is heavily dependent on the temperature profile and geometry of the reactor, many analyses either assume a constant power shape, or use neutronics analysis to determine a power shape for a specific reactor configuration. The development of a coupling interface for a propulsion system performance code and a Monte Carlo neutron transport code (OpenMC) allows for the reactor power shape to be calculated in a Picard iteration.

Jacob Stonehill↗

Integrated modelling: Coupling of surface evolution and plasma-impurity transport

During the interaction of the scrape off layer (SOL) plasma with the first wall the evolution of both wall and plasma are tightly coupled: The erosion of the first wall leads to an impurity concentration in the plasma which affects the particle and power balance in the plasma. In turn the impurities, when leaving the plasma via transport, can form deposits and mixed materials, far away from their initial source location. These deposits can be eroded, allowing the impurities to stepwise migrate through the fusion device until they end up at a location where the plasma at the wall is cold enough and no further erosion occurs. To describe these processes an integrated model of surface evolution and plasma transport of impurities is needed. The WallDYN code achieves this required coupling of processes by parameterising the output of surface evolution- and plasma impurity-migration-codes by analytical models. For a given fixed background plasma it evolves the surface composition, derives impurity flux into and from the plasma and can from this derive the impurity densities in the plasma. This paper will show the importance of including this recycling of impurities at the wall in impurity migration modelling: The 13 CH 4 seeding experiment performed in the Wendelstein 7-X Stellerator is modelled using the recent extension of the WallDYN code to 3D plasma and wall geometries. A comparison with post mortem analysis of the 13 C deposition shows both qualitative and quantitative agreement with the WallDYN calculations.

36 MATERIALS SCIENCE↗

High Temperature Steam Electrolysis Process Model

The code models the operation of a high temperature steam electrolysis (HTSE) plant coupled with a nuclear power plant, and provides calculations of the process energy requirements; water feed stream flow rate requirements; process cooling requirements; process conditions such as temperature, pressure, and flow rate, throughout the process; product purity and flow rate, as well as data that can be used to determine the sizing, performance, and cost specifications of individual equipment components that comprise the system.

Wendt, DanielS (0000000344914296)↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

The Schwarz alternating method for transient solid dynamics

Abstract In our earlier work, we formulated the Schwarz alternating method as a means for concurrent multiscale coupling in finite deformation solid mechanics for quasi‐static problems. Herein, we advance this method for the study of transient dynamic multiscale solid mechanics problems where information is exchanged back and forth between small and large scales. The extension to dynamics relies on the notion of a global time stepper. Within each global time step, the subdomains are coupled by the standard Schwarz iterative process. Remarkably, each subdomain can use its own time step or even its own time integrator to advance its solution in time, provided that they synchronize at each global time step. We study the performance of the Schwarz method on several examples designed for this purpose. Our numerical experiments demonstrate that the method is capable of coupling regions with different mesh resolutions, different element types, and different time integration schemes (e.g., implicit and explicit), all without introducing any artifacts that afflict other coupling methods for transient dynamics. Finally, we apply the dynamic Schwarz alternating method to the simulation of a bolted joint subjected to dynamic loading, as a demonstration of the performance of the method in a realistic scenario.

Mota, Alejandro↗

Characterizing the Roles of Biogeochemical Cycling and Ocean Circulation in Regulating Marine Copper Distributions

Copper (Cu) is a key micronutrient for marine phytoplankton. Its oceanic biogeochemical cycle has elicited considerable attention due to dissolved Cu exhibiting a unique linear profile with depth. Several processes have been proposed for explaining this behavior. Here, in this study, we characterize the relationships between the observed Cu, PO 4 3– , and Si on a global scale. We find that the depth profiles of Cu resemble those of Si more than of PO 4 3– in the global ocean. To understand their relationships, we couple the biogeochemical and internal circulation processes in a model, fitting optimal Cu:PO 4 3– uptake ratios and remineralization length-scales to replicate the marine Cu distributions. The modeling results suggest that Cu uptake needs in the Southern Ocean are substantially higher than those in other oceanic regions. In addition, our modeling results indicate a deep Cu remineralization in the global ocean. We offer an alternative mechanism that relies on biogeochemical cycling and internal circulation to produce the linear depth profiles of dissolved Cu. Our results suggest that diatoms are likely the major phytoplankton dominating oceanic Cu cycling.

54 ENVIRONMENTAL SCIENCES↗

Interaction of a finite-length ion beam with a background plasma - Reflected ions at the quasi-parallel bow shock

The coupling of a finite-length, field-aligned, ion beam with a uniform background plasma is investigated using one-dimensional hybrid computer simulations. The finite-length beam is used to study the interaction between the incident solar wind and ions reflected from the earth's quasi-parallel bow shock, where the reflection process may vary with time. The coupling between the reflected ions and the solar wind is relevant to ion heating at the bow shock and possibly to the formation of hot, flow anomalies and re-formation of the shock itself. Consistent with linear theory, the waves which dominate the interaction are the electromagnetic right-hand polarized resonant and nonresonant modes. However, in addition to the instability growth rates, the length of time that the waves are in contact with the beam is also an important factor in determining which wave mode will dominate the interaction. It is found that interaction will result in strong coupling, where a significant fraction of the available free energy is converted into thermal energy in a short time, provided the beam is sufficiently dense or sufficiently long.

Onsager, T. G.↗

Quasiparticle Couplings in Transport of Heat, Charge, and Spin for Novel Energy Materials. Final report

Understanding the microscopic processes involved in the transport and conversion of energy from the atomic-scale to the meso-scale is critical for the development of next-generation materials for energy conversion and information processing. The goal of this program was to elucidate the couplings between microscopic degrees-of-freedom of atomic vibrations, spins, and electrons, which underpin heat, charge, and spin transport. This project systematically investigated these quasi-particle interactions, using state-of-the-art neutron scattering techniques, combined with first-principles computer simulations. This research program uncovered new information and understanding, which helps to rationalize materials properties based on the interplay of their atomistic degrees of freedom.

36 MATERIALS SCIENCE↗