Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 955 records · Page 53

MMIC Power Amplifier Puts Out 40 mW From 75 to 110 GHz

A three-stage monolithic microwave integrated circuit (MMIC) W-band amplifier has been constructed and tested in a continuing effort to develop amplifiers as well as oscillators, frequency multipliers, and mixers capable of operating over wide frequency bands that extend above 100 GHz. There are numerous potential uses for MMICs like these in scientific instruments, radar systems, communication systems, and test equipment operating in this frequency range.

Samoska, Lorene↗

Applying Motion Constraints Based on Test Data

MSC ADAMS is a simulation software that is used to analyze multibody dynamics. Using user subroutines, it is possible to apply motion constraints to the rigid bodies so that they match the motion profile collected from test data. This presentation describes the process of taking test data and passing it to ADAMS using user subroutines, and uses the Morpheus free-flight 4 test as an example of motion data used for this purpose. Morpheus is the name of a prototype lander vehicle built by NASA that serves as a test bed for various experimental technologies (see backup slides for details) MSC.ADAMS"TM" is used to play back telemetry data (vehicle orientation and position) from each test as the inputs to a 6-DoF general motion constraint (details in backup slides) The MSC.ADAMS"TM" playback simulations allow engineers to examine and analyze flight trajectory as well as observe vehicle motion from any angle and at any playback speed. This facilitates the development of robust and stable control algorithms, increasing reliability and reducing development costs of this developmental engine The simulation also incorporates a 3D model of the artificial hazard field, allowing engineers to visualize and measure performance of the developmental autonomous landing and hazard avoidance technology ADAMS is a multibody dynamics solver. It uses forces, constraints, and mass properties to numerically integrate equations of motion. The ADAMS solver will ask the motion subroutine for position, velocity, and acceleration values at various time steps. Those values must be continuous over the whole time domain. Each degree of freedom in the telemetry data can be examined separately; however, linear interpolation of the telemetry data is invalid, since there will be discontinuities in velocity and acceleration.

Burlone, Michael↗

Students Planting Research Of Use To Space (SPROUTS) Networking and Engagement of Interns at Kennedy Space Center

While Kennedy Space Center (KSC) has had summer interns in the life sciences for many years, beginning in 2016, the internship program was re-envisioned. Rather than having many interns in the summer with numerous innovative small projects created primarily for this purpose, the current model is to offer internship opportunities every semester (Spring, Summer, and Fall), so that research is continuous and sustained. Since 2016 we have had 64 internship opportunities for undergraduate and graduate students in the areas of space crop production and microgravity simulation support. Some students have participated in more than one opportunity, so 52 separate individuals have participated in this program. While interns primarily worked with individual mentors on a main project, most also worked on group projects with other interns or developed smaller secondary projects on their own initiative. In addition, all interns participated in a weekly Lunch, Learn, and Discuss seminar series where they were encouraged to present their research, learn from scientists at KSC, and engage with members of the space life sciences community. When the face-to-face spring internship of 2020 ended prematurely, we resumed the weekly seminar series as a virtual event, and this led to the concept of engaging former interns. Current and former interns were invited to participate in the on-going virtual series, and 46 asked to be included on the mailing list. Our current seminar series involves remote seminars or hosted discussions from contacts within and outside of the agency, and engages a large network of current scientists and former interns with an weekly participation averaging 40-50 people. Speakers include experts from industry, academia, and other space agencies working on topics of interest to KSC space life science researchers. In addition, this series has served as a brain trust and brainstorming mechanism for KSC scientists to invite and discuss new ideas and concepts. The advent of remote work and virtual platforms has allowed us to engage with interns in a new way that will continue as future interns are integrated into this network of individuals interested in space research. This work is funded by NASA’s Space Biology Program.

Internships↗

Tracking Blade Tip Vortices for Numerical Flow Simulations of Hovering Rotorcraft

Blade tip vortices generated by a helicopter rotor blade are a major source of rotor noise and airframe vibration. This occurs when a vortex passes closely by, and interacts with, a rotor blade. The accurate prediction of Blade Vortex Interaction (BVI) continues to be a challenge for Computational Fluid Dynamics (CFD). Though considerable research has been devoted to BVI noise reduction and experimental techniques for measuring the blade tip vortices in a wind tunnel, there are only a handful of post-processing tools available for extracting vortex core lines from CFD simulation data. In order to calculate the vortex core radius, most of these tools require the user to manually select a vortex core to perform the calculation. Furthermore, none of them provide the capability to track the growth of a vortex core, which is a measure of how quickly the vortex diffuses over time. This paper introduces an automated approach for tracking the core growth of a blade tip vortex from CFD simulations of rotorcraft in hover. The proposed approach offers an effective method for the quantification and visualization of blade tip vortices in helicopter rotor wakes. Keywords: vortex core, feature extraction, CFD, numerical flow visualization

Numerical↗

Machine Learning-Assisted High-Temperature Reservoir Thermal Energy Storage Optimization: Numerical Modeling and Machine Learning Input and Output Files

This data set includes the numerical modeling input files and output files used to synthesize data, and the reduced-order machine learning models trained from the synthesized data for reservoir thermal energy storage site identification. In this study, a machine-learning-assisted computational framework is presented to identify High-Temperature Reservoir Thermal Energy Storage (HT-RTES) site with optimal performance metrics by combining physics-based simulation with stochastic hydrogeologic formation and thermal energy storage operation parameters, artificial neural network regression of the simulation data, and genetic algorithm-enabled multi-objective optimization. A doublet well configuration with a layered (aquitard-aquifer-aquitard) generic reservoir is simulated for cases of continuous operation and seasonal-cycle operation scenarios. Neural network-based surrogate models are developed for the two scenarios and applied to generate the Pareto fronts of the HT-RTES performance for four potential HT-RTES sites. The developed Pareto optimal solutions indicate the performance of HT-RTES is operation-scenario (i.e., fluid cycle) and reservoir-site dependent, and the performance metrics have competing effects for a given site and a given fluid cycle. The developed neural network models can be applied to identify suitable sites for HT-RTES, and the proposed framework sheds light on the design of resilient HT-RTES systems. All the simulations and the neural network model were done by Idaho National Laboratory. A detailed description of the work was reported in publication linked below.

15 GEOTHERMAL ENERGY↗

Study of a Plug Nozzle for Supersonic Aircraft Concepts

Nozzles with external plugs are candidates for the exhaust system of supersonic aircraft concepts in an ongoing research effort at NASA. As part of that effort a fundamental study is being conducted at NASA Glenn Research Center (GRC). Experimental and limited numerical simulation results from this study, for a plug nozzle with a convergent nozzle (cowl),have been reported in recent publications[1,2]. Conventional method-of-characteristics (MoC) design for the plug requires that the flow is choked at the nozzle exit. This is why a sharp convergence of the cowl is invoked in such a design in most previous studies [3-5], leading to the earlier choice of the convergent cowl. Such a nozzle, however, is not suitable for supersonic flights because of large boat tail drag suffered by its outer surface. The boat tail drag consideration dictates that the cowl be more or less cylindrical in shape. While various other shapes have been studied in the past [6], essentially a cylindrical shape is adopted in the larger research program at NASA. The cylindrical shape has certain advantages. With the plug crown located somewhat inside the nozzle, the flow experiences convergence and divergence as with a C-D nozzle. By translating the plug, different throat-to-exit area ratio could be achieved to obtain fully expanded condition at different values of NPR. This offers a relatively easy active control strategy for the nozzle geometry throughout the flight regime from landing and takeoff (LTO) to cruise conditions. A plug nozzle model mimicking the geometry being considered in the NASA program, i.e., having a cylindrical exit, have been fabricated recently for continued fundamental studies. Some of the plugs used with the earlier convergent cowl [1] are readily adopted to this model. In addition, a few other plugs have been obtained following a numerical optimization study [7]. The various configurations are to be investigated for thrust, noise and flow fields, with complementing numerical simulations. Thrust data for a variety of configurations as well as limited noise data are in hand. Further explorations including schlieren visualization are to be conducted in the near future. The proposed paper is to summarize all these results. In the following, a few key results obtained so far are described.

jets↗

Steering of Solar Sails Using Optical Lift Force

Optical wing structures were theoretically and numerically analyzed, and prototype arrays of wings called optical flying carpets were fabricated with solar sail material clear polyimide (CP1). This material was developed at NASA Langley to better withstand damaging ultraviolet radiation found in outer space. Various optical wing sizes and shapes were analyzed to develop design strategies for thrust and torque applications. The developed ray-tracing model has undergone continual advancement, and stands as an effective tool for modeling most types of solar sails. To our understanding, such a model does not exist else where. The distributed forces and torques have been reduced to a simple theoretical whereby the fundamental mechanics may be understood in terms of the numerically determined center of pressure off set from the center of mass. This description applies to any type of solar sail, affording our ray-tracing model a general utility. This research has established a foundation for understanding the force and torque afforded by optical wings. The study began by considering transparent wings and ended by considering wings having a reflecting face. The latter was found to afford the advantages of high thrust and both intrinsic and extrinsic torque. Our discovery of the intrinsic torque on optical wings (meaning that a moment arm is not required) has no analogy for a flat reflective solar sail, and therefore provides an extra degree of control that may be useful for sail craft attitude and navigation purposes.

Swartzlander, Gover↗

Aerodynamic influence coefficient method using singularity splines

A numerical lifting surface formulation, including computed results for planar wing cases is presented. This formulation, referred to as the vortex spline scheme, combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of loading function methods. The formulation employes a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfied all the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise. Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral. The current formulation uses the elementary horseshoe vortex as the basic singularity and is therefore restricted to linearized potential flow. As part of the study, a non planar development was considered, but the numerical evaluation of the lifting surface concept was restricted to planar configurations. Also, a second order sideslip analysis based on an asymptotic expansion was investigated using the singularity spline formulation.

Mercer, J. E.↗

Tomographic X-ray particle tracking velocimetry

We investigate the feasibility of in-laboratory tomographic X-ray particle tracking velocimetry (TXPTV) and consider creeping flows with nearly density matched flow tracers. Specifically, in these proof-of-concept experiments we examined a Poiseuille flow, flow through porous media and a multiphase flow with a Taylor bubble. For a full 360° computed tomography (CT) scan we show that the specially selected 60 micron tracer particles could be imaged in less than 3 seconds with a signal-to-noise ratio between the tracers and the fluid of 2.5, sufficient to achieve proper volumetric segmentation at each time step. In the pipe flow, continuous Lagrangian particle trajectories were obtained, after which all the standard techniques used for PTV or PIV (taken at visible wave lengths) could also be employed for TXPTV data. And, with TXPTV we can examine flows inaccessible with visible wave lengths due to opaque media or numerous refractive interfaces. In the case of opaque porous media we were able to observe material accumulation and pore clogging, and for flow with Taylor bubble we can trace the particles and hence obtain velocities in the liquid film between the wall and bubble, with thickness of liquid film itself also simultaneously obtained from the volumetric reconstruction after segmentation. While improvements in scan speed are anticipated due to continuing improvements in CT system components, we show that for the flows examined even the presently available CT systems could yield quantitative flow data with the primary limitation being the quality of available flow tracers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the Equilibrium Electrostatic Potential and Light-Induced Charge Redistribution in Halide Perovskite Structures

Lead halide perovskites are semiconductor materials which are employed as nonintentionally doped absorbers inserted between two selective carrier transport layers (SCTL), realizing a p-i-n or n-i-p heterojunction. In our study, we have developed and investigated a lateral device, based on methylammonium lead iodide (MAPbI3) in which the p-i-n heterojunction develops in the horizontal direction. Our research suggests that the effective doping level in the MAPbI3 film should be very low, below 1012 cm-3. Along the vertical direction, this doping level is not enough to screen the electric field of the buried heterojunction with the SCTL. The perovskite work function is therefore affected by the work function of the SCTL underneath. From drift-diffusion simulations, we show that intrinsic perovskite-SCTL structures develop mV range surface photovoltages (SPVs) under continuous illumination. However, perovskite-SCTL structures can develop SPVs of hundreds of mV, as confirmed by our measurements. We therefore analyzed the compatibility between low doping and low defect densities in the perovskite layer and such high SPV values using numerical modeling. It is shown that these high SPV values could originate from electronic processes due to large band offsets in the buried perovskite-SCTL heterojunctions, or at the SCTL-transparent conductive oxide (TCO) buried heterojunction. However, such electronic processes can hardly explain the long SPV persistence after switching off the illumination.

defects↗

Large-Scale and Global Hydrology

Powered by the sun, water moves continuously between and through Earths oceanic, atmospheric, and terrestrial reservoirs. It enables life, shapes Earths surface, and responds to and influences climate change. Scientists measure various features of the water cycle using a combination of ground, airborne, and space-based observations, and seek to characterize it at multiple scales with the aid of numerical models. Over time our understanding of the water cycle and ability to quantify it have improved, owing to advances in observational capabilities, the extension of the data record, and increases in computing power and storage. Here we present some of the most recent estimates of global and continental ocean basin scale water cycle stocks and fluxes and provide examples of modern numerical modeling systems and reanalyses.Further, we discuss prospects for predicting water cycle variability at seasonal and longer scales, which is complicated by a changing climate and direct human impacts related to water management and agriculture. Changes to the water cycle will be among the most obvious and important facets of climate change, thus it is crucial that we continue to invest in our ability to monitor it.

land hydrology↗

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗

Magnon-Magnon Interaction Induced by Dynamic Coupling in a Hybrid Magnonic Crystal

We report a combined experimental and numerical investigation of spin-wave dynamics in a hybrid magnonic crystal consisting of a CoFeB artificial spin ice (ASI) of stadium-shaped nanoelements patterned atop a continuous NiFe film separated by a 5 nm Al 2 O 3 spacer. Using Brillouin light scattering spectroscopy, we probe the frequency dependence of thermal spin waves as functions of applied magnetic field and wavevector, revealing the decisive role of interlayer dipolar coupling in the magnetization dynamics. Micromagnetic simulations complement the experiments, showing a strong interplay between ASI edge modes and backward volume modes in the NiFe film. The contrast in saturation magnetization between CoFeB and NiFe enhances this coupling, leading to a pronounced hybridization manifested as a triplet of peaks in the BLS spectra predicted by simulations and observed experimentally. This magnon−magnon coupling persists over a wide magnetic field range, shaping both the spin-wave dispersion at fixed fields and the full frequency-field response throughout the magnetic hysteresis loop. Our findings establish how ASI geometry can selectively enhance specific spin-wave wavelengths in the underlying film, thereby boosting their amplitude and identifying them as preferential channels for spin wave transmission and manipulation.

Brillouin light scattering spectroscopy↗

Benchmark for two-dimensional large scale coherent structures in partially magnetized E × B plasmas—community collaboration & lessons learned

Low-temperature plasmas (LTPs) are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of LTP codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E × B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development.

benchmarking↗

Modeling and Market Design Considerations for Conventional and Decarbonized Resources

The following paper reviews various market design issues that will need to be reconsidered due to anticipated changes in the resources that supply energy in wholesale electricity markets. As the energy supply continues towards decarbonization, the change in resource technologies will affect the fundamentals of production scheduling and the policies to address supply variability and uncertainty. We show through simple numerical examples that the existing market design approach based on fuel costs will result in $0/MWh prices with intermittent price spikes during reserve shortages, but that incorporating more granular reserve pricing, energy storage participation, and price-responsive demand participation can restore efficient market clearing with reasonable pricing outcomes. The market design approach fundamentally shifts from fuel based to opportunity cost based. We briefly review alternative market design frameworks that can support this shift, including intraday markets, decentralized markets, flexibility options, and swing contracts. Market designs may also be required to accommodate or support various out-of-market policies and agreements; ideally, these external factors can be integrated into the market design to facilitate efficient exchanges across longer time scales, between other markets, and in support of public policy goals. The paper concludes by discussing how decarbonization trends may affect the design and use of production cost models for short term operations and long term planning studies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Electronic Interactions Between the Receptor-Binding Domain of Omicron Variants and Angiotensin-Converting Enzyme 2: A Novel Amino Acid–Amino Acid Bond Pair Concept

SARS-CoV-2 remains a severe threat to worldwide public health, particularly as the virus continues to evolve and diversify into variants of concern (VOCs). Among these VOCs, Omicron variants exhibit unique phenotypic traits, such as immune evasion, transmissibility, and severity, due to numerous spike protein mutations and the rapid subvariant evolution. These Omicron subvariants have more than 15 mutations in the receptor-binding domain (RBD), a region of the SARS-CoV-2 spike protein that is important for recognition and binding with the angiotensin-converting enzyme 2 (ACE2) human receptor. To address the impact of these high numbers of Omicron mutations on the binding process, we have developed a novel method to precisely quantify amino acid interactions via the amino acid–amino acid bond pair (AABP). We applied this concept to investigate the interface interactions of the RBD–ACE2 complex in four Omicron Variants (BA.1, BA.2, BA.5, and XBB.1.16) with its Wild Type counterpart. Based on the AABP analysis, we have identified all the sites that are affected by mutation and have provided evidence that unmutated sites are also impacted by mutation. We have calculated that the binding between RBD and ACE2 is strongest in OV BA.1, followed by OV BA.2, WT, OV BA.5, and OV XBB.1.16. We also present the partial charge values for all 311 residues across these five models. Our analysis provides a detailed understanding of changes caused by mutation in each Omicron interface complex.

Biochemistry & Molecular Biology↗

Periodic orbits in the Planar General Three-Body Problem

The article contains a numerical study of periodic solutions of the Planar General Three-Body Problem. Several new periodic solutions have been discovered and are described. In particular, there is a continuous family with variable masses, extending all the way from the elliptic restricted problem to the general problem with three equal masses. All our examples have special symmetry properties which are described in detail. Finally we also suggest some important applications to the natural satellites of the solar system.

Broucke, R.↗

Compliant wall surface motion and its effect on the structure of a turbulent boundary layer

The status of the continuing compliant wall drag reduction research at NASA-Langley Research Center is discussed. Preliminary surface motion calculations are reported along with compliant surface design concepts and their numerical models. A compliant drag reduction theory based on stabilizing the turbulent substructure is proposed and previous experiments have been examined relative to that theory. Results of recent low speed compliant surface experiments have been reported which include initial attempts to measure local compliant surface motion.

Ash, R. L.↗