Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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 739 records · Page 41

Nonlinear Reduced-Order Simulation Using An Experimentally Guided Modal Basis

A procedure is developed for using nonlinear experimental response data to guide the modal basis selection in a nonlinear reduced-order simulation. The procedure entails using nonlinear acceleration response data to first identify proper orthogonal modes. Special consideration is given to cases in which some of the desired response data is unavailable. Bases consisting of linear normal modes are then selected to best represent the experimentally determined transverse proper orthogonal modes and either experimentally determined inplane proper orthogonal modes or the special case of numerically computed in-plane companions. The bases are subsequently used in nonlinear modal reduction and dynamic response simulations. The experimental data used in this work is simulated to allow some practical considerations, such as the availability of in-plane response data and non-idealized test conditions, to be explored. Comparisons of the nonlinear reduced-order simulations are made with the surrogate experimental data to demonstrate the effectiveness of the approach.

Rizzi, Stephen A.↗

Modeling of Melt Growth During Carbothermal Processing of Lunar Regolith

The carbothermal processing of lunar regolith has been proposed as a means to produce carbon monoxide and ultimately oxygen to support human exploration of the moon. In this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. Carbon monoxide is further processed in other reactors downstream to ultimately produce oxygen. The amount of oxygen produced crucially depends on the amount of regolith that is molten. In this paper we develop a model of the heat transfer in carbothermal processing. Regolith in a suitable container is heated by a heat flux at its surface such as by continuously shining a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and its temperature rises until it attains the melting point. The energy from the heat flux is then used for the latent heat necessary to change phase from solid to liquid, after which the temperature continues to rise. Thus a small melt pool appears under the heated zone shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat and supplies the energy required to melt more of the regolith, and the size of the molten zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the melt is balanced by radiative losses from the surface. In this paper, we model the melting and the growth of the melt zone with time in a bed of regolith when a portion of its surface is subjected to a constant heat flux. The heat flux is assumed to impinge on a circular area. Our model is based on an axisymmetric three-dimensional variation of the temperature field in the domain. Heat transfer occurs only by conduction, and effects of convective heat transport are assumed negligible. Radiative heat loss from the surface of the melt and the regolith to the surroundings is permitted. We perform numerical computations to determine the shape and the mass of the melt at steady state and its time evolution. We first neglect the volume change upon melting, and subsequently perform calculations including it. Predictions from our model are compared to test data to determine the effective thermal conductivities of the regolith and the melt that are compatible with the data

Balasubramaniam, R.↗

Blocking Losses With a Photon Counter

It was not known how to assess accurately losses in a communications link due to photodetector blocking, a phenomenon wherein a detector is rendered inactive for a short time after the detection of a photon. When used to detect a communications signal, blocking leads to losses relative to an ideal detector, which may be measured as a reduction in the communications rate for a given received signal power, or an increase in the signal power required to support the same communications rate. This work involved characterizing blocking losses for single detectors and arrays of detectors. Blocking may be mitigated by spreading the signal intensity over an array of detectors, reducing the count rate on any one detector. A simple approximation was made to the blocking loss as a function of the probability that a detector is unblocked at a given time, essentially treating the blocking probability as a scaling of the detection efficiency. An exact statistical characterization was derived for a single detector, and an approximation for multiple detectors. This allowed derivation of several accurate approximations to the loss. Methods were also derived to account for a rise time in recovery, and non-uniform illumination due to diffraction and atmospheric distortion of the phase front. It was assumed that the communications signal is intensity modulated and received by an array of photon-counting photodetectors. For the purpose of this analysis, it was assumed that the detectors are ideal, in that they produce a signal that allows one to reproduce the arrival times of electrons, produced either as photoelectrons or from dark noise, exactly. For single detectors, the performance of the maximum-likelihood (ML) receiver in blocking is illustrated, as well as a maximum-count (MC) receiver, that, when receiving a pulse-position-modulated (PPM) signal, selects the symbol corresponding to the slot with the largest electron count. Whereas the MC receiver saturates at high count rates, the ML receiver may not. The loss in capacity, symbol-error-rate (SER), and count-rate were numerically computed. It was shown that the capacity and symbol-error-rate losses track, whereas the count-rate loss does not generally reflect the SER or capacity loss, as the slot-statistics at the detector output are no longer Poisson. It is also shown that the MC receiver loss may be accurately predicted for dead times on the order of a slot.

Moision, Burce E.↗

Zero-gravity Mean Free Surface Curvature of a Confined Liquid in a Radially-Vaned Container

A variety of increasingly intricate container geometries are under consideration for the passive manipulation of liquids aboard spacecraft where the impact of gravity may be neglected. In this study we examine the mean curvature of a liquid volume confined within a radial array of disconnected vanes of infinite extent. This particular geometry possesses a number of desirable characteristics relevant to waste water treatment aboard spacecraft for life support. It is observed that under certain conditions the slender shape of the free surface approaches an asymptote, which can be predicted analytically using new hybrid boundary conditions proposed herein. This contribution represents possibly the final extension of what has been referred to as the method of de Lazzer et al. (1996). The method enables the integration of the Young-Laplace equation over a domain with its boundaries, including the wetted portion of the solid boundaries, symmetry planes, and circular arcs representing free surfaces at the center plane of the liquid body. Asymptotic solutions at several limits are obtained and the analysis is confirmed with numerical computations.

Chen, Yongkang↗

Entropy Stable Wall Boundary Conditions for the Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite volume, finite difference, discontinuous Galerkin, and flux reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗

A Numerical Theory for Impedance Education in Three-Dimensional Normal Incidence Tubes

A method for educing the locally-reacting acoustic impedance of a test sample mounted in a 3-D normal incidence impedance tube is presented and validated. The unique feature of the method is that the excitation frequency (or duct geometry) may be such that high-order duct modes may exist. The method educes the impedance, iteratively, by minimizing an objective function consisting of the difference between the measured and numerically computed acoustic pressure at preselected measurement points in the duct. The method is validated on planar and high-order mode sources with data synthesized from exact mode theory. These data are then subjected to random jitter to simulate the effects of measurement uncertainties on the educed impedance spectrum. The primary conclusions of the study are 1) Without random jitter the method is in excellent agreement with that for known impedance samples, and 2) Random jitter that is compatible to that found in a typical experiment has minimal impact on the accuracy of the educed impedance.

Watson, Willie R.↗

Mechanical Prototyping and Manufacturing Internship

The internship was located at the Johnson Space Center (JSC) Innovation Design Center (IDC), which is a facility where the JSC workforce can meet and conduct hands-on innovative design, fabrication, evaluation, and testing of ideas and concepts relevant to NASA's mission. The tasks of the internship included mechanical prototyping design and manufacturing projects in service of research and development as well as assisting the users of the IDC in completing their manufacturing projects. The first project was to manufacture hatch mechanisms for a team in the Systems Engineering and Project Advancement Program (SETMAP) hexacopter competition. These mechanisms were intended to improve the performance of the servomotors and offer an access point that would also seal to prevent cross-contamination. I also assisted other teams as they were constructing and modifying their hexacopters. The success of this competition demonstrated a proof of concept for aerial reconnaissance and sample return to be potentially used in future NASA missions. I also worked with Dr. Kumar Krishen to prototype an improved thermos and a novel, portable solar array. Computer-aided design (CAD) software was used to model the parts for both of these projects. Then, 3D printing as well as conventional techniques were used to produce the parts. These prototypes were then subjected to trials to determine the success of the designs. The solar array is intended to work in a cluster that is easy to set up and take down and doesn't require powered servomechanisms. It could be used terrestrially in areas not serviced by power grids. Both projects improve planetary exploration capabilities to future astronauts. Other projects included manufacturing custom rail brackets for EG-2, assisting engineers working on underwater instrument and tool cases for the NEEMO project, and helping to create mock-up parts for Space Center Houston. The use of the IDC enabled efficient completion of these projects at significantly reduced cost. I acquired and improved manufacturing and prototyping skills during my tour including learning about a CAD (Computer-Aided Design) program called Creo (Creo Parametric; design software), gaining valuable conventional machining experience with lathes, CNC (Computer Numerical Control) milling machines and various other tools, and improving my engineering project communication and collaboration skills. The internship also allowed me to better understand operations at NASA. I plan to work in the aerospace industry or do academic research benefitting space science and exploration, and this internship experience will enable me to have insight into manufacturing processes for research and development.

Grenfell, Peter↗

On the Development of the Marshall Grazing Incidence X-ray Spectrograph (MaGIXS) Mirrors

The Marshall Grazing Incidence X-ray Spectrograph (MaGIXS) is a sounding rocket experiment that will obtain spatially resolved soft X-ray spectra of the solar corona from 0.5 - 2 keV. The optical system comprises a Wolter-I telescope mirror, a slit spectrograph, and a CCD camera. The spectrograph has a finite conjugate paraboloid pair, which re-images the slit, and a varied line-space planar reflection grating. Both the Wolter-I mirror and paraboloid pair are being fabricated at the NASA Marshall Space Flight Center (MSFC), using nickel replication. The MaGIXS mirror mandrels have been diamond turned, polished, and have yielded a set of engineering mirrors. Unlike other grazing incidence instruments, such as FOXSI, ART-XC, and IXPE, the MaGIXS prescriptions have large departure from a cone. This property exacerbates challenges with conventional lap polishing techniques and interferometric metrology. Here we discuss the progression of the optical surfaces of the mandrels through lap polishing, X-ray data from the replicated shells obtained in the MSFC Stray Light Facility (SLF), and our transition to using the ZEEKO computer numerical controlled (CNC) polisher for figure correction.

Champey, Patrick↗

Deterministic Polishing of Replicating Grazing-Incidence Mandrels

In an effort to manufacture high-angular-resolution, grazing-incidence, x-ray optics, Marshall Space Flight Center (MSFC) is taking measures to improve its electroformed replicated optics. A key development is the use of computer-numerical control (CNC) polishing to deterministically improve the surface of electroless nickel mandrels used to replicate grazing- incidence optics. Metrology, control software and polishing parameters must function together seamlessly to reach the specifications required to replicate sub-arcsecond optics. Each change in polishing parameters effects the wear pattern of the polishing head. Using Richardson-Lucy deconvolution, the controller software fits the wear pattern to metrology data to calculate the changing feedrates across the mandrel. Here we present an overview of our process, and early results showing the effectiveness of deterministic polishing for replicated optics.

Davis, Jacqueline↗

X-ray Evaluation of the Marshall Grazing Incidence X-Ray Spectrometer (MaGIXS) Nickel-Replicated Mirrors

X-ray observations of astronomical objects provides diagnostics not available in any other wavelength regime, however the capability of making these observation at a high spatial resolution has proven challenging. Recently, NASA Marshall Space Flight Center (MSFC) has made good progress in employing computer numerical control (CNC) polishing techniques on electroless nickel mandrels as part of our replicated grazing incidence optics program. CNC polishing has afforded the ability to deterministically refine mandrel figure, thereby improving mirror performance. The Marshall Grazing Incidence X-ray Spectrometer (MaGIXS) is a MSFC-led sounding rocket instrument that is designed to make the first ever soft x-ray spectral observations of the Sun spatially resolved along a narrow slit. MaGIXS incorporates some of the first mirrors produced at MSFC using this polishing technique. Here we present the predicted mirror performance obtained from metrology, after completion of CNC polishing, as well as the results of X-ray tests performed on the MaGIXS telescope mirror before and after mounting.

MaGIXS X-ray Tests↗

Machine Learning-Based Predictive Analytics for Aircraft Engine Conceptual Design

Big data and artificial intelligence/machine learning are transforming the global business environment. Data is now the most valuable asset for enterprises in every industry. Companies are using data-driven insights for competitive advantage. With that, the adoption of machine learning-based data analytics is rapidly taking hold across various industries, producing autonomous systems that support human decision-making. This work explored the application of machine learning to aircraft engine conceptual design. Supervised machine-learning algorithms for regression and classification were employed to study patterns in an existing, open-source database of production and research turbofan engines, and resulting in predictive analytics for use in predicting performance of new turbofan designs. Specifically, the author developed machine learning-based analytics to predict cruise thrust specific fuel consumption (TSFC) and core sizes of high-efficiency turbofan engines, using engine design parameters as the input. The predictive analytics were trained and deployed in Keras, an open-source neural networks application program interface (API) written in Python, with Google’s TensorFlow (an open source library for numerical computation) serving as the backend engine. The promising results of the predictive analytics show that machine-learning techniques merit further exploration for application in aircraft engine conceptual design.

deep-learning↗

Analysis of the Amplitude Response of a Two-Point and a Multi-Point Focused Laser Differential Interferometer

Experiments were performed with two different focused laser differential interferometer (FLDI) instruments to assess the relative amplitude response of each instrument to a traveling shockwave generated by a laser spark in air. The first FLDI instrument generated two slightly-separated measurement points, with each measurement point providing sensitivity to density fluctuations induced by the shockwave. By performing a cross-correlation between the signals obtained from each measurement point, the phase velocity of density fluctuations can be obtained. The second FLDI instrument generated a measurement line that was oriented parallel to the shockwave direction of travel with the resulting interference pattern sampled at multiple equally-spaced points along the line. As a result, this instrument provides density fluctuation measurement capability at multiple points simultaneously. When the measurement line is oriented parallel to the shockwave direction of travel, the phase velocity, rate of change of the phase velocity, and acceleration of density fluctuations traveling along each line can be obtained by performing a cross-correlation between points along the line. Numerical computations of the shockwave generated by the laser spark are used to simulate the response of each instrument and are compared to the experimental results. High-speed schlieren imaging has also been performed and is compared with the FLDI measurements and computational results.

Brett F. Bathel↗

Analysis of the Amplitude Response of a Two-Point and a Multi-Point Focused Laser Differential Interferometer

Experiments were performed with two different focused laser differential interferometer (FLDI) instruments to assess the relative amplitude response of each instrument to a traveling shockwave generated by a laser spark in air. The first FLDI instrument generated two slightly-separated measurement points, with each measurement point providing sensitivity to density fluctuations induced by the shockwave. By performing a cross-correlation between the signals obtained from each measurement point, the phase velocity of density fluctuations can be obtained. The second FLDI instrument generated a measurement line that was oriented parallel to the shockwave direction of travel with the resulting interference pattern sampled at multiple equally-spaced points along the line. As a result, this instrument provides density fluctuation measurement capability at multiple points simultaneously. When the measurement line is oriented parallel to the shockwave direction of travel, the phase velocity, rate of change of the phase velocity, and acceleration of density fluctuations traveling along each line can be obtained by performing a cross-correlation between points along the line. Numerical computations of the shockwave generated by the laser spark are used to simulate the response of each instrument and are compared to the experimental results. High-speed schlieren imaging has also been performed and is compared with the FLDI measurements and computational results.

Interferometry↗

Formal Methods in the Development of Highly Assured Software for Unmanned Aircraft Systems

In traditional software development methodologies, operational and functional requirements of systems are often specified in structured natural language notations. These restricted notations provide good documentation support, but only provide limited support for semantic analysis. These notations are generally not rich enough to unambiguously specify the requirements of safety-critical systems that, for example, involve complex numerical computations or that interact with the physical environment. Examples of these safety-critical systems are autonomous vehicles such as unmanned aircraft systems. This talk advocates the use of expressive formal logics, such as higher-order logic, to specify the operational and functional requirement of unmanned systems and to prove the correctness of these requirements. Semantic analysis of requirements written in higher-order logic is supported through the use of interactive theorem provers. Formal models serve as ideal reference implementations of functional requirements. Hence, formal logics enable software validation techniques where software implementations can be checked against functional requirements in a mechanical way. The Formal Methods group in the Safety-Critical Avionics Systems Branch at NASA Langley Research Center has conducted research on the development and application of formal verification techniques to safety-critical applications of interest to NASA for more than 30 years. This talk illustrates the use of formal methods in the development of highly-assured autonomous unmanned aircraft systems.

Formal Methods↗

An Expanded Role for Python in Expediting System Simulation Development

The Python programming language has traditionally been used as a “scripting” language and is not generally recognized as a language for building system simulations, where C++, Java, and other compiled languages are typically used. While very powerful and gaining favor for a large set of programming tasks, Python is generally regarded as not having the speed to directly do the extensive numeric computations required in large-scale simulation. This paper aims to present an approach, supported by case studies, to re-think this assertion. Directly using Python in this new role has the potential to expedite an agile, spiral cycle of development of design-test-modify for systems simulation.

programming language↗

X-57 “Maxwell” High-Lift Propeller Testing and Model Development

NASA’s X-57 “Maxwell” distributed electric propulsion flight demonstrator has a high-lift system that includes 12 fixed-pitch high-lift propellers located upstream of the wing leading edge for lift augmentation at low speeds. These high-lift propellers are not required at higher speeds and are folded conformally along the nacelles to reduce drag when not in operation. Aircraft performance models and flight simulations incorporate propeller performance to predict thrust, moments, and power consumption and propeller model accuracy is important in identifying safe operating regimes for this aircraft. Current high-lift propeller performance models have been verified and calibrated against numerous computational fluid dynamics analyses under a variety of flight conditions. We performed a series of full-scale wind tunnel tests at the NASA Langley Research Center Low Speed Aeroacoustic Wind Tunnel to further validate these models, to identify any adverse operating conditions for these propellers, and to assess the accuracy of the facility’s new Propeller Test Stand. The results indicate that the models accurately predicted performance and that the right-handed propeller showed lower torque and thrust for a given propeller speed compared to the left-handed version. Acoustic data support these measurements, showing higher tonal noise for the left versus right propeller. We suspect this is due to slight differences in the material properties between the early production (right) and more recent (left) blade sets leading to greater detwisting in the right-handed blades under load. Both propellers demonstrated very stable operation throughout the test including during deployment, stowing, and windmilling. Furthermore, high-speed imagery of the spinning propeller indicated no blade oscillations or instabilities for any test case. We also leveraged these images to analyze the effects of propeller loading on blade detwisting and tip position. Unfortunately, the thrust measurements varied significantly during the test due to cross-loading of the Propeller Test Stand load cell, though trends in the thrust data were still observable. The torque measurement, however, remained very precise for all test points and was used for performance comparisons. As a result of this test, we validated and improved the X-57 high-lift propeller performance models, obtained detailed propeller acoustic data and high-speed imagery, observed the propeller operating safely during a variety of dynamic events, and identified two potentially adverse, low-propeller-speed conditions which will be avoided during normal aircraft operation.

X-57↗

Extinction of the Stagnation Point Diffusion Flame: Effect of Conductive Heat Loss into Solid Interior

The planned microgravity experiment Growth and Extinction Limit (GEL) is designed to study the material flammability and extinction of thermally thick solid with an emphasis to quantify the effect of conductive heat loss into the solid interior by measuring the in-depth temperature profile. To obtain a more detailed understanding of the extinction phenomena, a previously developed one-dimensional axisymmetric stagnation-point diffusion flame model is expanded to account for the excessive heat loss into the solid beyond the steady burning state. This is important because for thick solids, this is important because the burning period can have a long and slow solid transient with a quasi-steady gas -phase. A non-dimensional excessive conductive heat loss parameter Ψ is defined in the paper and extensive numerical computations have been performed to determine the extinction boundary and the characteristics of near-limit flames. The classical U-shaped flammability boundary in the oxygen vs. stretch rate plot now becomes a series of U-shaped boundaries as a function of Ψ. For a given oxygen ambient, extinction boundary is a function of Ψ. The boundary consists of a blowoff branch and a radiative quenching branch. This is similar to the steady burning flammability boundary (Ψ =0) in the oxygen vs stretch rate plot. Extensive computed data including flame and pyrolysis temperatures, burning rate, flame standoff distance, and reaction rate are presented along the extinction boundaries. In particular, we note that the non-dimensional gas reaction rate is the ratio of diffusion reaction time to the reaction diffusion time, so it has the classical Damhkoler number character. By examining the maximum non-dimensional reaction rates, we observe that along the extinction boundaries, they are approximate constant across the entire range of stretch rates (from quenching to blow-off). These critical values are of the order unity but may vary in different oxygen ambient. The implication will be discussed in the presentation together with the near-limit flame structures.

Chengyao Li↗