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 1,603 records · Page 89

Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks

Genetic and evolutionary algorithms have been applied to solve numerous problems in engineering design where they have been used primarily as optimization procedures. These methods have an advantage over conventional gradient-based search procedures became they are capable of finding global optima of multi-modal functions and searching design spaces with disjoint feasible regions. They are also robust in the presence of noisy data. Another desirable feature of these methods is that they can efficiently use distributed and parallel computing resources since multiple function evaluations (flow simulations in aerodynamics design) can be performed simultaneously and independently on ultiple processors. For these reasons genetic and evolutionary algorithms are being used more frequently in design optimization. Examples include airfoil and wing design and compressor and turbine airfoil design. They are also finding increasing use in multiple-objective and multidisciplinary optimization. This lecture will focus on an evolutionary method that is a relatively new member to the general class of evolutionary methods called differential evolution (DE). This method is easy to use and program and it requires relatively few user-specified constants. These constants are easily determined for a wide class of problems. Fine-tuning the constants will off course yield the solution to the optimization problem at hand more rapidly. DE can be efficiently implemented on parallel computers and can be used for continuous, discrete and mixed discrete/continuous optimization problems. It does not require the objective function to be continuous and is noise tolerant. DE and applications to single and multiple-objective optimization will be included in the presentation and lecture notes. A method for aerodynamic design optimization that is based on neural networks will also be included as a part of this lecture. The method offers advantages over traditional optimization methods. It is more flexible than other methods in dealing with design in the context of both steady and unsteady flows, partial and complete data sets, combined experimental and numerical data, inclusion of various constraints and rules of thumb, and other issues that characterize the aerodynamic design process. Neural networks provide a natural framework within which a succession of numerical solutions of increasing fidelity, incorporating more realistic flow physics, can be represented and utilized for optimization. Neural networks also offer an excellent framework for multiple-objective and multi-disciplinary design optimization. Simulation tools from various disciplines can be integrated within this framework and rapid trade-off studies involving one or many disciplines can be performed. The prospect of combining neural network based optimization methods and evolutionary algorithms to obtain a hybrid method with the best properties of both methods will be included in this presentation. Achieving solution diversity and accurate convergence to the exact Pareto front in multiple objective optimization usually requires a significant computational effort with evolutionary algorithms. In this lecture we will also explore the possibility of using neural networks to obtain estimates of the Pareto optimal front using non-dominated solutions generated by DE as training data. Neural network estimators have the potential advantage of reducing the number of function evaluations required to obtain solution accuracy and diversity, thus reducing cost to design.

Rai, Man Mohan↗

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗

What are grana in chloroplasts of vascular plants good for?

All plants and green algae contain stacked grana thylakoid membranes in their chloroplasts, underscoring an evolutionary pressure to maintain this unique structural feature. In addition, numerous studies have demonstrated that particular lateral and vertical dimensions of grana facilitate the function, regulation and repair of the photosynthetic machinery responsible for energy conversion. In this review, we present an updated overview of our understanding concerning the structure of grana thylakoids, the forces that contribute to their formation and their architectural dynamics. After establishing the structural foundation, we continue by reviewing the implications of grana formation on light harvesting, electron transport and protein maintenance in the thylakoid membranes of vascular plants. Regarding light harvesting, we discuss the implications of grana formation on energy spillover, macromolecular crowding, lateral versus vertical excitation energy transfer, and state transition. With respect to electron transport, we summarize recent findings regarding the functional connectivity of spatially separated photosystems facilitated by grana formation through small mobile electron carriers. We also explore how the dynamic responses of grana architecture regulate electron transport. Finally, we address how grana formation contributes to the organization of protein repair and biogenesis within thylakoid membranes.

59 BASIC BIOLOGICAL SCIENCES↗

Finite-Time Teleportation Phase Transition in Random Quantum Circuits

How long does it take to entangle two distant qubits in a quantum circuit evolved by generic unitary dynamics? Here, we show that if the time evolution is followed by measurements of all but two infinitely separated test qubits, then the entanglement between them can undergo a phase transition and become nonzero at a finite critical time t c . The fidelity of teleporting a quantum state from an input qubit to an infinitely distant output qubit shows the same critical onset. Specifically, these finite-time transitions occur in short-range interacting two-dimensional random unitary circuits and in sufficiently long-range interacting one-dimensional circuits. The phase transition is understood by mapping the random continuous-time evolution to a finite-temperature thermal state of an effective spin Hamiltonian, where the inverse temperature equals the evolution time in the circuit. In this framework, the entanglement between two distant qubits at times t > t c corresponds to the emergence of long-range ferromagnetic spin correlations below the critical temperature. We verify these predictions using numerical simulation of Clifford circuits and propose potential realizations in existing platforms for quantum simulation.

36 MATERIALS SCIENCE↗

The COMBLE campaign: a study of marine boundary-layer clouds in Arctic cold-air outbreaks

One of the most intense air mass transformations on Earth happens when cold air flows from frozen surfaces to much warmer open water in cold-air outbreaks (CAOs), a process captured beautifully in satellite imagery. Despite the ubiquity of the CAO cloud regime over high-latitude oceans, we have a rather poor understanding of its properties, its role in energy and water cycles, and its treatment in weather and climate models. The Cold-air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) was conducted to better understand this regime and its representation in models. COMBLE aimed to examine the relations between surface fluxes, boundary-layer structure, aerosol, cloud and precipitation properties, and mesoscale circulations in marine CAOs. Processes affecting these properties largely fall in a range of scales where boundary-layer processes, convection, and precipitation are tightly coupled, which makes accurate representation of the CAO cloud regime in numerical weather prediction and global climate models most challenging. COMBLE deployed an Atmospheric Radiation Measurement Mobile Facility at a coastal site in northern Scandinavia (69°N), with additional instruments on Bear Island (75°N), from December 2019 to May 2020. CAO conditions were experienced 19% (21%) of the time at the main site (on Bear Island). A comprehensive suite of continuous in situ and remote sensing observations of atmospheric conditions, clouds, precipitation, and aerosol were collected. Because of the clouds’ well-defined origin, their shallow depth, and the broad range of observed temperature and aerosol concentrations, the COMBLE dataset provides a powerful modeling test bed for improving the representation of mixed-phase cloud processes in large-eddy simulations and large-scale models.

54 ENVIRONMENTAL SCIENCES↗

Temperature Effects on Droplet Oscillation Decay with Application to Fuel Property Measurement

Observation of oscillation decay in droplets has been shown to be an effective approach in determining physical properties such as viscosity and surface tension for emerging biofuels using μl quantities. Herein this work extends the approach to higher temperatures relevant to fuel injection conditions for internal combustion engines. Experiments are conducted that use high-resolution strobed imaging of moving, heated, μm-sized, fuel droplets to capture shape oscillation decay through image processing and analysis. Two fuels are investigated, iso-butanol and a primary reference fuel (PRF 84), which is a mixture of iso-octane and n-heptane. A piezoelectric droplet generator is used to generate a continuous train of single droplets, which are given an initial perturbation and subsequently undergo damped oscillations, captured using strobed imaging. Surface tension and viscosity calculated at four different temperatures ranging from 30°C-55°C using the frequency and decay time associated with the fundamental mode are found to be within ~ 10% of reference values obtained from the literature. Complementary numerical simulations are performed that utilize a volume-of-fluid approach to track transient droplet oscillation phenomena along with heat and mass transfer in a 2D axisymmetric domain. Simulations, where droplet size and temperature can be independently varied, capture an expected decrease in oscillation frequency and increase in decay time, with increase in fuel temperature. The simulations are further used to study the relative contribution to deviations in surface tension and viscosity predictions due to heat and mass transfer from the droplet, as well as viscous effects violating the inviscid flow assumption in the droplet oscillation theory. Mass loss effects are found to be negligible. Temperature change due to heat transfer is found to have the next highest sensitivity, particularly for the more volatile PRF 84, which has a lower viscosity and associated Ohnesorge number. For isobutanol, which has a higher viscosity and Ohnesorge number, viscous effects contribute the most to deviation in fuel property predictions.

09 BIOMASS FUELS↗

Buckling analysis for structural sections and stiffened plates reinforced with laminated composites.

A classical buckling analysis is developed for stiffened, flat plates composed of a series of linked flat plate and beam elements. Plates are idealized as multilayered orthotropic elements; structural beads and lips are idealized as beams. The loaded edges of the stiffened plate are simply supported and the conditions at the unloaded edges can be prescribed arbitrarily. The plate and beam elements are matched along their common junctions for displacement continuity and force equilibrium in an exact manner. Offsets between elements are considered in the analysis. Buckling under uniaxial compressive load for plates, sections and stiffened plates is investigated. Buckling loads are found as the lowest of all possible general and local failure modes and the mode shape is used to determine whether buckling is a local or general instability. Numerical correlations with existing analysis and test data for plates, sections and stiffened plates including boron-reinforced structures are discussed. In general, correlations are reasonably good.

Viswanathan, A. V.↗

On signal design by the R sub 0 criterion for non-white Gaussian noise channels

The use of the R sub 0 criterion for modulation system design is investigated for channels with non-white Gaussian noise. A signal space representation of the waveform channel is developed, and the cut-off rate R sub 0 for vector channels with additive nonwhite Gaussian noise and unquantized demodulation is derived. When the signal unput to the channel is a continuous random vector, maximization of R sub 0 with constrained average signal energy leads to a water-filling interpretation of optimal energy distribution in signal space. The necessary condition for a finite signal set to maximize R sub 0 with constrained energy and an equally likely probability assignment of signal vectors is presented, and an algorithm is outlined for numerically computing the optimum signal set. A necessary condition on a constrained energy, finite signal set is found which maximizes a Taylor series approximation of R sub 0. This signal set is compared with the finite signal set which has the water-filling average energy distribution.

Bordelon, D. L.↗

An explicit solution to the exoatmospheric powered flight guidance and trajectory optimization problem for rocket propelled vehicles

A derivation of an explicit solution to the two point boundary-value problem of exoatmospheric guidance and trajectory optimization is presented. Fixed initial conditions and continuous burn, multistage thrusting are assumed. Any number of end conditions from one to six (throttling is required in the case of six) can be satisfied in an explicit and practically optimal manner. The explicit equations converge for off nominal conditions such as engine failure, abort, target switch, etc. The self starting, predictor/corrector solution involves no Newton-Rhapson iterations, numerical integration, or first guess values, and converges rapidly if physically possible. A form of this algorithm has been chosen for onboard guidance, as well as real time and preflight ground targeting and trajectory shaping for the NASA Space Shuttle Program.

Jaggers, R. F.↗

Simulation of Heavy Lift Airship dynamics over large ranges of incidence and speed

A nonlinear, multibody, six-degrees-of-freedom digital simulation has been developed to study generic Heavy Lift Airship (HLA) dynamics and control. The basic aerodynamic functions developed to model the hull, tail, and rotor loads continuously over all incidence ranges are reviewed and applied to a Quadrotor HLA with a low fineness ratio hull and a small vee-tail. Trim calculations for a test vehicle suggest control power deficiencies in crosswind stationkeeping for the unloaded vehicle. Gust responses show the importance of correctly calculating loads due to accelerated relative motion of air and hull. Numerically linearized dynamics for the test vehicle show the existence of a divergent yaw mode, and an oscillatory pitch mode whose stability characteristics are sensitive to flight speed. A considerable improvement in the vehicle's stability and response results from a simple multi-axis closed-loop control system operating on the rotors and propeller blades.

Tischler, M. B.↗

Thin-plate spline quadrature of geodetic integrals

Thin-plate spline functions (known for their flexibility and fidelity in representing experimental data) are especially well-suited for the numerical integration of geodetic integrals in the area where the integration is most sensitive to the data, i.e., in the immediate vicinity of the evaluation point. Spline quadrature rules are derived for the contribution of a circular innermost zone to Stoke's formula, to the formulae of Vening Meinesz, and to the recursively evaluated operator L(n) in the analytical continuation solution of Molodensky's problem. These rules are exact for interpolating thin-plate splines. In cases where the integration data are distributed irregularly, a system of linear equations needs to be solved for the quadrature coefficients. Formulae are given for the terms appearing in these equations. In case the data are regularly distributed, the coefficients may be determined once-and-for-all. Examples are given of some fixed-point rules. With such rules successive evaluation, within a circular disk, of the terms in Molodensky's series becomes relatively easy. The spline quadrature technique presented complements other techniques such as ring integration for intermediate integration zones.

Vangysen, Herman↗

Prediction of vortical flows on wings using incompressible Navier-Stokes equations

Numerical solutions to the incompressible Navier-Stokes equations are in greater demand than ever before as the field of computational fluid dynamics (CFD) increases its impact as an engineering tool. An implicit upwind-relation finite-difference scheme is developed for solving unsteady incompressible Navier-Stokes equations incorporating an artificial time derivative of the pressure to the continuity equation. The essential features observed in experiments for complicated vortical flows past three round-edged low-aspect-ratio wings at high angles of attack are reproduced both qualitatively and quantitatively.

Hsu, C.-H.↗

Re-forming supercritical quasi-parallel shocks. II - Mechanism for wave generation and front re-formation

This paper continues the study of Thomas et al. (1990) in which hybrid simulations of quasi-parallel shocks were performed in one and two spatial dimensions. To identify the wave generation processes, the electromagnetic structure of the shock is examined by performing a number of one-dimensional hybrid simulations of quasi-parallel shocks for various upstream conditions. In addition, numerical experiments were carried out in which the backstreaming ions were removed from calculations to show their fundamental importance in reformation process. The calculations show that the waves are excited before ions can propagate far enough upstream to generate resonant modes. At some later times, the waves are regenerated at the leading edge of the interface, with properties like those of their initial interactions.

Winske, D.↗

Current and future uses of accelerators in particle astrophysics

Beams of artificially accelerated heavy ions, protons, antiprotons, electrons, and positrons currently available at (and planned for) numerous facilities around the world are a valuable resource to the Cosmic Ray community. Such beams have been used to test detector concepts, calibrate balloon-borne and space flight experiments and to measure fundamental nuclear physics parameter necessary for the interpretation of Cosmic Ray data. As new experiments are flown the quality and extent of Cosmic Ray measurements will continue to improve. It will be necessary to increase activity at ground based accelerators in order to test/calibrate these new instruments and to maintain (or possibly improve) the ability to interpret these data. In this area, the newly formed Transport Collaboration, supported by NASA, will be providing new nuclear interaction cross section measurements for beams with Z less than or = 58 and supporting new instrument calibrations at the Lawrence Berkeley Laboratory Bevalac accelerator.

Guzik, T. G.↗

Identification of Gravity-Related Effects on Crystal Growth From Melts With an Immiscibility Gap

This work involves an experimental-numerical approach to study the effects of natural and Marangoni convections on solidification of single crystals from a silicate melt with a liquid-liquid immiscibility gap. Industrial use of crystals grown from silicate melts is becoming increasingly important in electronic, optical, and high temperature structural applications. Even the simplest silicate systems like Al203-SiO2 have had, and will continue to have, a significant role in the development of traditional and advanced ceramics. A unique feature of crystals grown from the silicate systems is their outstanding linear electro-optic properties. They also exhibit exceptionally high optical rotativity. As a result, these crystals are attractive materials for dielectric, optical, and microwave applications. Experimental work in our laboratory has indicated that directional solidification of a single crystal mullite appears to be preceded by liquid-liquid phase separation in the melt. Disruption of the immiscible state results in crystallization of a two phase structure. There is also evidence that mixing in the melt caused by density-driven convection can significantly affect the stability of the immiscible liquid layers and result in poly-crystalline growth. On earth, the immiscible state has only been observed for small diameter crystals grown in float zone systems where natural convection is almost negligible. Therefore, it is anticipated that growth of large single crystals from silicate melts would benefit from microgravity conditions because of the reduction of the natural convective mixing. The main objective of this research is to determine the effects of transport processes on the phase separation in the melt during growth of a single crystal while addressing the following issues: (1) When do the immiscible layers form and are they real?; (2) What are the main physical characteristics of the immiscible liquids?; and (3) How mixing by natural or Marangoni convection affects the stability of the phase separated melt.

Kassemi, M.↗

Monitoring a Complex Physical System using a Hybrid Dynamic Bayes Net

The Reverse Water Gas Shift system (RWGS) is a complex physical system designed to produce oxygen from the carbon dioxide atmosphere on Mars. If sent to Mars, it would operate without human supervision, thus requiring a reliable automated system for monitoring and control. The RWGS presents many challenges typical of real-world systems, including: noisy and biased sensors, nonlinear behavior, effects that are manifested over different time granularities, and unobservability of many important quantities. In this paper we model the RWGS using a hybrid (discrete/continuous) Dynamic Bayesian Network (DBN), where the state at each time slice contains 33 discrete and 184 continuous variables. We show how the system state can be tracked using probabilistic inference over the model. We discuss how to deal with the various challenges presented by the RWGS, providing a suite of techniques that are likely to be useful in a wide range of applications. In particular, we describe a general framework for dealing with nonlinear behavior using numerical integration techniques, extending the successful Unscented Filter. We also show how to use a fixed-point computation to deal with effects that develop at different time scales, specifically rapid changes occuring during slowly changing processes. We test our model using real data collected from the RWGS, demonstrating the feasibility of hybrid DBNs for monitoring complex real-world physical systems.

Lerner, Uri↗

Website of the Systems and Analysis Branch Supported Projects and Graph Analysis

Throughout the past few weeks I have learned a great amount of information about many interesting aspects that go on here at NASA Glenn Research Center Branch 7820. Branch 7820 is the Systems and Analysis Branch. The people involved in this Branch deal with in a nutshell the analysis of propulsion systems for Earth to orbit and space transportation systems. The first project that I had worked on was helping my mentor learn more about lunar geography and the most recommended way to maintain communication for our future lunar missions. During this time I studied the craters of the moon, especially the South Pole, to provide her with information so that she can make decisions. I also researched to provide her with contact information on those people who are specialized in lunar geography so that she may talk to them to find out more in depth information. Most of my time spent here has been helping to develop a comprehensive explanation and background of the different projects our Branch has supported. When I first came to NASA Glenn and started working with the 7820 Branch there website had many holes that needed to be filled in. I have spent numerous weeks researching information about topics such as Project Prometheus itself and one of its components Jupiter s Icy Moons Orbiter (JIMO). I have also done a large amount of research on propulsion systems and how different kinds work. I have learned many facts about Nuclear Electric Propulsion (NEP) all the way to Nuclear Thermal Propulsion (NTP) systems. I will continue to do this until all the holes are filled and find out about Global Integrated Design Environment (GLIDE) and Next Generation Launch technology (NGLT). Since most of my job was providing information to go onto a website I has to learn how to put my information into a HTML format. I had no previous knowledge on how to do that kind of task and had to study how to do it and am now able to create a document in HTML format. There has been reorganizing done here at NASA Glenn and our Branch was moved to another building. Therefore, our library had to be moved with us. I spent time helping to put together the boxes, pack the library, and label them accordingly. This was not an easy task but was an experience in itself. I was able to see old posters that NASA had produced about different space missions and look at Russian map of the US and books on space missions. I was also about to see what was in the library in terms of reference material helped because now I can make use of the information for my research on the website. Throughout my internship my mentor will provide me with graphs to analyze and recreate so that she may use them to her advantage. I will learn from every piece of data that comes my way. Later, I will study and analysis gravity-loss for Earth departure trajectories. Since I haven't done that yet I cannot really describe what 1 will learn or what exactly the project entails. The whole experience has been great and I have no doubt that it will exceed every expectation previously thought.

Kellerman, Corinne↗