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 703 records · Page 39

An Integrated Platform for Mission Performance: Goddard Space Flight Center's Meta Information System

Successful management of NASA Goddard Space Flight Centers portfolio of mission projects (i.e. spacecraft, spacecraft instruments / subsystems, and ground systems) depends upon numerous processes, data and information resources coming together throughout the life-cycle of the projects to enable sound decision-making, risk reduction, conformance with requirements and continual improvement. In the summer of 2012, GSFC SMA management recognized that the set of legacy data systems used to support the SMA mission and responsibilities at GSFC were not effective, integrated nor on par with best practices. Accordingly, GSFC SMA management launched a project to transform the GSFC electronic management system to create and deploy an integrated, state-of-the-market information system for process performance, data management and analytics that is currently operating in support of GSFC overall mission performance. The information system, known as Meta, provides controlled, process-based applications, integrated data management, and analytics for a robust set of integrated applications. Since the inception of the project, a lean, award-winning SMA project team has thus far designed, configured and managed the system to serve as an integrated platform for 18 applications that support key processes for GSFC and the Science Mission Directorate at NASA headquarters, such as internal Management System assessments, Supplier Assessments, Supplier Research & Analysis, Problem Reporting for in-house mission projects, On-orbit Anomaly Reporting for GSFC-operated spacecraft, GSFC Project Life-Cycle Reviews, Risk Management and Reporting, and SMA Project Team Reporting. The system has experienced steady growth in utilization with 1,917 users (NASA employees and support contractor personnel) and 44 gigabytes of data and related files managed by the system as of the end of April 2018. Overall, Meta has achieved increasing levels of process management and real-time data / information integration to enable improved process performance, richer analytics and informed decision-making. This paper will present and examine the need, goals, approach and design driving the management and operation of the Meta information system and its ongoing growth as an integrated platform for advancing NASA mission performance.

Integration↗

Advanced Technology Infusion into a Spacesuit System

Advancement in technology drives our future. Moreover, the successful implementation of a technology drives its possibilities. The National Aeronautics and Space Administration (NASA) has invested in numerous technologies over the decades. A quote from NASA’s Space Technology Mission Directorates indicates the importance of technology: “Historically, technology has driven humanity’s progress and will continue to define our future. Our nation chooses to invest in new technology not only to maintain our edge in the global economy but also because technology helps us: Redefine the possible; Create a technologically advanced future; and Drive economic growth.” It is difficult for a technology to satisfy these goals unless it can be successfully infused into a system. For a technology to continue to evolve, become a reality, and infuse into NASA’s missions, there must exist a compilation of success-oriented factors for the technology to reach fruition. Understanding these factors could help decrease the complexity of technology infusion and bridge the gap between technology developers and system integrators. Ultimately, the knowledge gained could facilitate the design, development, and infusion of a technology to be more effective and efficient. Successful technology infusion is complex. Even complexity at the interfaces can be daunting when advanced technology’s infuse into complex systems such as a spacesuit. Additionally, unforeseen emergent behavior and challenges can surface upon integration. Overall, there is a need to understand and measure the success of infusing an advanced technology into a complex system. Industry and academia desire to quantify performance benefits of technologies integrated into complex systems. This paper will focus advanced technology infusion into a spacesuit system. This research could help NASA and industry’s project managers and system mangers integrate advanced technologies more effectively and efficiently.

Advanced Technology↗

Clustering Approach To Identifying Low Lunar Frozen Orbits In A High-Fidelity Model

Low lunar frozen orbits are of continued interest in the astrodynamics community for trajectory design and space domain awareness. This paper presents a data-driven approach to analyzing a wide variety of numerically-generated lunar trajectories in a 100 × 100 lunar gravity model with the point mass gravity of the Earth and Sun. First, clustering is used to extract a summary of these trajectories: within each cluster, trajectories possess a geometrically similar evolution of perilune but varying drift and lifetimes. Within some clusters, trajectories with a bounded perilune evolution are also identified to produce candidates for lunar frozen orbits of distinct geometries.

lunar frozen orbit gravity model data clustering↗

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗

A numerical model of magnetosphere-ionosphere coupling Preliminary results

A three-dimensional simulation model was developed to study magnetosphere-ionosphere coupling in the auroral region. One-fluid MHD equations are adopted to model the magnetosphere, and current density continuity equations were solved consistently to model the ionosphere. In the preliminary simulation runs described here, the electrodynamics of region 1 field-aligned currents were modeled. Initially, the electric field is taken to be in the magnetospheric equatorial plane. Alfven waves then propagate down to the ionosphere, accompanying a field-aligned current and exciting the electrostatic potential in the ionosphere by electrodynamic coupling. The results of the preliminary runs directly correspond with the fundamental characteristics of global magnetosphere-ionosphere coupling. These characteristics include an ionospheric electrostatic potential which varies in its development in time, depending upon the ratio of ionospheric resistance to magnetospheric impedance, the flowing of field-aligned currents into the ionosphere on the dawnside and out on the duskside, and the distribution along the geomagnetic field of the amplitude of the field-aligned current density in proportion to the field intensity.

Watanabe, K.↗

Numerical simulation of turbomachinery flows with advanced turbulence models

A three dimensional full Navier-Stokes (FNS) code is used to simulate complex turbomachinery flows. The code incorporates an explicit multistep scheme and solves a conservative form of the density averaged continuity, momentum, and energy equations. A compressible low Reynolds number form of the k-epsilon turbulence model, and a q-omega model and an algebraic Reynolds stress model have been incorporated in a fully coupled manner to approximate Reynolds stresses. The code is used to predict the viscous flow field in a backswept transonic centrifugal compressor for which laser two focus data is available. The code is also used to simulate the tip clearance flow in a cascade. The code has been extended to include unsteady Euler solutions for predicting the unsteady flow through a cascade due to incoming wakes, simulating rotor-stator interactions.

Lakshminarayana, B.↗

The NASA Low-Pressure Turbine Flow Physics Program: A Review

An overview of the NASA Glenn Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The flow in the LPT is unique for the gas turbine. It is characterized by low Reynolds number and high freestream turbulence intensity and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The flow of most interest is on the suction surface, where large losses are generated due to separation. The LPT is a large, multistage, heavy, jet engine component that suffers efficiency degradation between takeoff and cruise conditions due to decrease in Reynolds number with altitude. The performance penalty is around 2 points for large commercial bypass engines and as much as 7 points for small, high cruise altitude, military engines. The gas-turbine industry is very interested in improving the performance of the LPT and in reducing its weight, part count and cost. Many improvements can be accomplished by improved airfoil design, mainly by increasing the airfoil loading that can yield reduction of airfoils and improved performance. In addition, there is a strong interest in reducing the design cycle time and cost. Key enablers of the needed improvements are computational tools that can accurately predict LPT flows. Current CFD tools in use cannot yet satisfactorily predict the unsteady, transitional and separated flow in the LPT. The main reasons are inadequate transition & turbulence models and incomplete understanding of the LPT flow physics. NASA Glenn has established its LPT program to answer these needs. The main goal of the program is to develop and assess models for unsteady CFD of LPT flows. An approach that consists of complementing and augmenting experimental and computational work elements has been adopted. The work is performed in-house and by several academic institutions, in cooperation and interaction with industry. The program was reviewed at the Minnowbrook II meeting in 1997. This review will summarize the progress that was made since and will introduce newly started projects. The LPT program is focused on three areas: acquisition of experimental and numerical databases and on modeling and computation. Priority was initially given to experiments. There are three classes of experiments: simulated LPT passages, linear cascade, both with and without wakes, and low-speed rotating rig. They are being conducted as follows: At NASA GRC on a flat surface with blade pressure distribution, at the US Naval Academy on a curved surface. The addition of wakes is studied at the University of Minnesota in a curved passage with a retractable wake generator, and at Texas A&M University in a linear cascade with continuously running wake generator. The pressure distribution of the Pratt & Whitney blade 'Pak B' is used in all these experiments. Experiments have been performed also in the GEAE Low-Speed Rotating Turbine (LSRT) rig with GE-designed airfoils. Work on numerically generated database is in progress at the University of Kentucky, using the DNS/LES code LESTool developed there. Turbulence/transition model assessment and development is performed also at the University of Kentucky, where a new intermittency transport model was developed and many experimental test cases have been numerically computed. Assessments of models using simulations of multistage LPT experiments were performed at Virginia Commonwealth University using the Corsair code. Work on suction surface separation delay, using passive and active flow-control, has also been initiated. Following the overview, Principal Investigators attending the workshop will present in detail several of the projects supported by NASA.

Ashpis, David E.↗

Research Progress and Accomplishments on ISS

The first research payloads reached the International Space Station (ISS) more than two years ago, with research operating continuously since March 2001. Seven research racks are currently on-orbit, with three more arriving soon to expand science capabilities. Through the first five expeditions, 60 unique NASA-managed investigations from 11 nations have been supported, many continuing into later missions. More than 90,000 experiment hours have been completed, and more than 1,000 hours of crew time have been dedicated to research, numbers that grow daily. The multidisciplinary program includes research in life sciences, physical sciences, biotechnology, Earth sciences, technology demonstrations as well as commercial endeavors and educational activities. The Payload Operations and Integration Center monitors the onboard activities around the clock, working with numerous Principal Investigators and Payload Developers at their remote sites. Future years will see expansion of the station with research modules provided by the European Space Agency and Japan, which will be outfitted with additional research racks. The first research payloads arrived at ISS more than two years ago, and continuous science has been ongoing for more than one and a half years. During this time, the research capabilities have been tremendously increased, even as assembly of the overall platform continues. Despite significant challenges along the way, ISS continues to successfully support a large number of investigations in a variety of research disciplines. The results of some of the early investigations are reaching the publication stage. The near future looms with new challenges, but experience to date and dedicated efforts give reason to be optimistic that the challenges will be overcome and that new and greater successes will be added to past ones.

Roe, Lesa B.↗

Three-dimensional stability of an elliptical vortex in a straining field

The three-dimensional linear stability of a rectilinear vortex of elliptical cross-section existing as a steady state in an irrotational straining field is studied numerically in the case of finite strain. It is shown that the instability predicted analytically for weak strain persists for finite strain and that the weak-strain results continue to be quantitatively valid for finite strain. The dependence of the growth rates of the unstable modes on the strain and the axial-disturbance wavelength is discussed. It is also shown that a three-dimensional instability is always more unstable than a two-dimensional instability in the range of parameters of most interest.

Robinson, A. C.↗

A time-dependent theoretical model of the polar wind Preliminary results

The coupled, time-dependent continuity, momentum, and energy equations of a two-ion O(+) and H(+) quasi-neutral plasma were solved in order to extend understanding of polar wind behavior. This numerical code allows studies of the time dependent behavior of polar wind-type flows into and out of the ionosphere. Initial studies indicate that the typical time constants for electron and ion temperature changes are of the order of minutes and tens of minutes, respectively. The response time of the minor high altitude ion O(+) is less than an hour, whereas that of the major ion, H(+), is many hours. The initial test runs also demonstrate the fact that temporary supersonic flows of both O(+) and H(+) are possible, especially in the presence of significant ion heating.

Gombosi, T. I.↗

A spline-based parameter estimation technique for static models of elastic structures

The problem of identifying the spatially varying coefficient of elasticity using an observed solution to the forward problem is considered. Under appropriate conditions this problem can be treated as a first order hyperbolic equation in the unknown coefficient. Some continuous dependence results are developed for this problem and a spline-based technique is proposed for approximating the unknown coefficient, based on these results. The convergence of the numerical scheme is established and error estimates obtained.

Dutt, P.↗

A spline-based parameter estimation technique for static models of elastic structures

The problem of identifying the spatially varying coefficient of elasticity using an observed solution to the forward problem is considered. Under appropriate conditions this problem can be treated as a first order hyperbolic equation in the unknown coefficient. Some continuous dependence results are developed for this problem and a spline-based technique is proposed for approximating the unknown coefficient, based on these results. The convergence of the numerical scheme is established and error estimates obtained.

Dutt, P.↗

Compensating for pneumatic distortion in pressure sensing devices

A general numerical technique for obtaining unsteady pressure measurements using conventional pressure sensing technology has been developed. A pneumatic distortion model, based on the Navier-Stokes equations of momentum and continuity, was reduced to a low-order, state-variable model retaining most of the dynamic characteristics of the full model. The reduced-order model is coupled with results from minimum variance estimation theory to develop an algorithm to compensate the effects of pneumatic distortion. Both postflight and real-time algorithms were developed and evaluated using simulated and flight data.

Whitmore, Stephen A.↗

A time-dependent dusty gas dynamic model of axisymmetric cometary jets

The present time-dependent, axisymmetric dusty gas dynamical model of inner cometary atmospheres solves the coupled and time-dependent equations of continuity, momentum, and energy for a gas-dust mixture between the surface of the nucleus and 100 km, using an axisymmetric 40 x 40 grid structure. A novel numerical method employing a second-order accurate Godunov-type scheme with dimensional splitting is used to solve the time-dependent pde system. It is established that a subsolar dust spike not predicted by previous calculations is generated by narrow axisymmetric jets, together with a jet cone whose opening angle depends on the jet length.

Korosmezey, A.↗

Nonlinearly stable compact schemes for shock calculations

The applications of high-order, compact finite difference methods in shock calculations are discussed. The main concern is to define a local mean which will serve as a reference for introducing a local nonlinear limiting to control spurious numerical oscillations while maintaining the formal accuracy of the scheme. For scalar conservation laws, the resulting schemes can be proven total-variation stable in one space dimension and maximum-norm stable in multiple space dimensions. Numerical examples are shown to verify accuracy and stability of such schemes for problems containing shocks. These ideas can also be applied to other implicit schemes such as the continuous Galerkin finite element methods.

Cockburn, Bernardo↗

Transport Phenomena During Equiaxed Solidification of Alloys

Recent progress in modeling of transport phenomena during dendritic alloy solidification is reviewed. Starting from the basic theorems of volume averaging, a general multiphase modeling framework is outlined. This framework allows for the incorporation of a variety of microscale phenomena in the macroscopic transport equations. For the case of diffusion dominated solidification, a simplified set of model equations is examined in detail and validated through comparisons with numerous experimental data for both columnar and equiaxed dendritic growth. This provides a critical assessment of the various model assumptions. Models that include melt flow and solid phase transport are also discussed, although their validation is still at an early stage. Several numerical results are presented that illustrate some of the profound effects of convective transport on the final compositional and structural characteristics of a solidified part. Important issues that deserve continuing attention are identified.

Beckermann, C.↗

A Preconditioning Method for Shape Optimization Governed by the Euler Equations

We consider a classical aerodynamic shape optimization problem subject to the compressible Euler flow equations. The gradient of the cost functional with respect to the shape variables is derived with the adjoint method at the continuous level. The Hessian (second order derivative of the cost functional with respect to the shape variables) is approximated also at the continuous level, as first introduced by Arian and Ta'asan (1996). The approximation of the Hessian is used to approximate the Newton step which is essential to accelerate the numerical solution of the optimization problem. The design space is discretized in the maximum dimension, i.e., the location of each point on the intersection of the computational mesh with the airfoil is taken to be an independent design variable. We give numerical examples for 86 design variables in two different flow speeds and achieve an order of magnitude reduction in the cost functional at a computational effort of a full solution of the analysis partial differential equation (PDE).

Arian, Eyal↗