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,207 records · Page 67

Geotail MCA Plasma Wave Investigation Data Analysis

The primary goals of the International Solar Terrestrial Physics/Global Geospace Science (ISTP/GGS) program are identifying, studying, and understanding the source, movement, and dissipation of plasma mass, momentum, and energy between the Sun and the Earth. The GEOTAIL spacecraft was built by the Japanese Institute of Space and Astronautical Science and has provided extensive measurements of entry, storage, acceleration, and transport in the geomagnetic tail and throughout the Earth's outer magnetosphere. GEOTAIL was launched on July 24, 1992, and began its scientific mission with eighteen extensions into the deep-tail region with apogees ranging from around 60 R(sub e) to more than 208 R(sub e) in the period up to late 1994. Due to the nature of the GEOTAIL trajectory which kept the spacecraft passing into the deep tail, GEOTAIL also made 'magnetopause skimming passes' which allowed measurements in the outer magnetosphere, magnetopause, magnetosheath, bow shock, and upstream solar wind regions as well as in the lobe, magnetosheath, boundary layers, and central plasma sheet regions of the tail. In late 1994, after spending nearly 30 months primarily traversing the deep tail region, GEOTAIL began its near-Earth phase. Perigee was reduced to 10 R(sub e) and apogee first to 50 R(sub e) and finally to 30 R(sub e) in early 1995. This orbit provides many more opportunities for GEOTAIL to explore the upstream solar wind, bow shock, magnetosheath, magnetopause, and outer magnetosphere as well as the near-Earth tail regions. The WIND spacecraft was launched on November 1, 1994 and the POLAR spacecraft was launched on February 24, 1996. These successful launches have dramatically increased the opportunities for GEOTAIL and the GGS spacecraft to be used to conduct the global research for which the ISTP program was designed. The measurement and study of plasma waves have made and will continue to make important contributions to reaching the ISTP/GGS goals and solving the significant problems of sun-earth connections. Plasma waves are involved in the energization and de-energization of plasma and energetic particles via numerous wave-particle interaction processes. Plasma waves in many instances are the source for the heating or cooling of the particles. They can cause particle precipitation by scattering particles into the loss cone. They move particles across boundaries in mass and energy dependent ways. Identifying the waves and the instabilities which produce them are thus crucial for understanding the plasma processes. Wave-particle interaction processes are especially important at various boundaries between the different regions of geospace including the bow shock, magnetopause, and interfaces in the geomagnetic tail between the magnetosheath, lobe, plasmasheet, boundary layers, and neutral sheet. In addition to identifying the characteristics of the instabilities and generation mechanisms encountered, plasma wave measurement are used in conjunction with other fields and particle measurements to identify the region of space the spacecraft is in or the boundary that is being crosed.

Anderson, Roger R.↗

Multidimensional MHD Model Studies of the Ionospheres of Venus and Mars

Continuing efforts have been made towards an increased understanding of the solar wind interaction and ionospheric processes at Venus and Mars. This work centered on a systematic development of a new generation of three dimensional magnetohydrodynamic (MHD) numerical code, which models the interaction processes of the solar wind with non-magnetic planets, such as Venus and Mars. We have also worked on a number of different, more specific and discrete studies, as various opportunities arose. We have developed new numerical codes for magnetospheric and cometary studies. As a first step in this process we built an axisymmetric model in which the solar wind interacts with a hard, perfectly conducting sphere. Even that model provided, in certain respects, significant improvements over previous ones.

Nagy, Andrew↗

A Systems Model for Power Technology Assessment

A computer model is under continuing development at NASA Glenn Research Center that enables first-order assessments of space power technology. The model, an evolution of NASA Glenn's Array Design Assessment Model (ADAM), is an Excel workbook that consists of numerous spreadsheets containing power technology performance data and sizing algorithms. Underlying the model is a number of databases that contain default values for various power generation, energy storage and power management and distribution component parameters. These databases are actively maintained by a team of systems analysts so that they contain state-of-art data as well as the most recent technology performance projections. Sizing of the power subsystems can be accomplished either by using an assumed mass specific power (W/kg) or energy (Wh/kg) or by a bottoms-up calculation that accounts for individual component performance and masses. The power generation, energy storage and power management and distribution subsystems are sized for given mission requirements for a baseline case and up to three alternatives. This allows four different power systems to be sized and compared using consistent assumptions and sizing algorithms. The component sizing models contained in the workbook are modular so that they can be easily maintained and updated. All significant input values have default values loaded from the databases that can be over-written by the user. The default data and sizing algorithms for each of the power subsystems are described in some detail. The user interface and workbook navigational features are also discussed. Finally, an example study case that illustrates the model's capability is presented.

Hoffman, David J.↗

An Improved Design for Air Removal from Aerospace Fluid Loop Coolant Systems

Aerospace applications with requirements for large capacity heat removal (launch vehicles, platforms, payloads, etc.) typically utilize a liquid coolant fluid as a transport media to increase efficiency and flexibility in the vehicle design. An issue with these systems however, is susceptibility to the presence of noncondensable gas (NCG) or air. The presence of air in a coolant loop can have numerous negative consequences, including loss of centrifugal pump prime, interference with sensor readings, inhibition of heat transfer, and coolant blockage to remote systems. Hardware ground processing to remove this air is also cumbersome and time consuming which continuously drives recurring costs. Current systems for maintaining the system free of air are tailored and have demonstrated only moderate success. An obvious solution to these problems is the development and advancement of a passive gas removal device, or gas trap, that would be installed in the flight cooling system simplifying the initial coolant fill procedure and also maintaining the system during operations. The proposed device would utilize commercially available membranes thus increasing reliability and reducing cost while also addressing both current and anticipated applications. In addition, it maintains current pressure drop, water loss, and size restrictions while increasing tolerance for pressure increases due to gas build-up in the trap.

Ritchie, Stephen M. C.↗

The Distribution of Loads on Rivets Connecting a Plate to a Beam under Transverse Loads

This report gives theoretical discussion of the distribution of leads on rivets connecting a plate to a beam under transverse leads. Two methods of solution are given which are applicable to loads up to the limit of proportionality; in the first the rivets are treated as discrete members, and in the second they are replaced by a continuous system of jointing. A method of solution is also given which is applicable to the case when nonlinear deformations occur in the rivets and the plate, but not in the beam. The methods are illustrated by numerical examples, and these show that the loads carried by the rivets and the plate are less than the values given by classical theory, which does not take into account the slip of the rivets, even below the limit of proportionality. The difference is considerably accentuated when nonlinear deformations occur in the restructure and the beam then carries the greater portion of the bending moment. If the material of the beam has a higher proportional limit and a higher ultimate strength than the material of the plate, there is thus a transfer of load from weaker to stronger material, and this is to the advantage of the structure. The methods given are of simple application and are recommended for use in the design of light-alloy structures when the design lead is likely to be above the proportional limit.

Vogt, F.↗

Efficient computation of optimal low thrust gravity perturbed orbit transfers

We have developed a new method for solving low-thrust fuel-optimal orbit transfer problems in the vicinity of a large body (planet or asteroid), considering a high-fidelity spherical harmonic gravity model. The algorithm is formulated via the indirect optimization method, leading to a two-point boundary value problem (TPBVP). We make use of a hyperbolic tangent smoothing law for performing continuation on the thrust magnitude to reduce the sharpness of the control switches in early iterations and thus promote convergence. The TPBVP is solved using the method of particular solutions (MPS) shooting method and Picard-Chebyshev numerical integration. Application of Picard-Chebyshev integration affords an avenue for increased efficiency that is not available with step-by-step integrators. We demonstrate that computing the particular solutions with only a low-fidelity force model greatly increases the efficiency of the algorithm while ultimately achieving near machine precision accuracy. A salient feature of the MPS is that it is parallelizable, and thus further speedups are available. It is also shown that, for near Earth orbits and over a small number of en-route revolutions around the Earth, only the zonal perturbation terms are required in the costate equations to obtain a solution that is accurate to machine precision and optimal to engineering precision. The proposed framework can be used for trajectory design around small asteroids and also for orbit debris rendezvous and removal tasks.

Junkins, John L.↗

Growth of zinc selenide single crystals by physical vapor transport in microgravity

The goals of this research were the optimization of growth parameters for large (20 mm diameter and length) zinc selenide single crystals with low structural defect density, and the development of a 3-D numerical model for the transport rates to be expected in physical vapor transport under a given set of thermal and geometrical boundary conditions, in order to provide guidance for an advantageous conduct of the growth experiments. In the crystal growth studies, it was decided to exclusively apply the Effusive Ampoule PVT technique (EAPVT) to the growth of ZnSe. In this technique, the accumulation of transport-limiting gaseous components at the growing crystal is suppressed by continuous effusion to vacuum of part of the vapor contents. This is achieved through calibrated leaks in one of the ground joints of the ampoule. Regarding the PVT transport rates, a 3-D spectral code was modified. After introduction of the proper boundary conditions and subroutines for the composition-dependent transport properties, the code reproduced the experimentally determined transport rates for the two cases with strongest convective flux contributions to within the experimental and numerical error.

Rosenberger, Franz↗

C deg continuity elements by Hybrid Stress method

An intensive study of the assumed variable distribution necessary for the Assumed Displacement Formulation, the Hellinger-Reissner Formulation, and the Hu-Washizu Formulation is made in a unified manner. With emphasis on physical explanation, a systematic method for the Hybrid Stress element construction is outlined. The numerical examples use four and eight node plane stress elements and eight and twenty node solid elements. Computation cost study indicates that the hybrid stress element derived using recently developed Uncoupled Stress Formulation is comparable in CPU time to the Assumed Displacement element. Overall, main emphasis is placed on providing a broader understanding of the Hybrid Stress Formulation.

Kang, David Sung-Soo↗

A simple numerical method for the calculation of the laminar boundary layer

An iteration difference method for the calculation of the incompressible laminar boundary layer is described. The method uses Prandtl's boundary layer equation and the boundary conditions directly and permits the attainment of an arbitrary accuracy. The method has been tested successfully in the continuation of the Blasius profile on the flat plate, on the circular cylinder investigated by Heimenz and on an elliptical cylinder of fineness ratio 1:4. The method makes possible the testing of previously developed methods, all of which contain important assumptions.

BOUNDARY LAYER, LAMINAR - CALCULATIONS - SCHRODER ↗

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗

Risk-averse optimal control of semilinear elliptic PDEs

In this paper, we consider the optimal control of semilinear elliptic PDEs with random inputs. These problems are often nonconvex, infinite-dimensional stochastic optimization problems for which we employ risk measures to quantify the implicit uncertainty in the objective function. In contrast to previous works in uncertainty quantification and stochastic optimization, we provide a rigorous mathematical analysis demonstrating higher solution regularity (in stochastic state space), continuity and differentiability of the control-to-state map, and existence, regularity and continuity properties of the control-to-adjoint map. Our proofs make use of existing techniques from PDE-constrained optimization as well as concepts from the theory of measurable multifunctions. We illustrate our theoretical results with two numerical examples motivated by the optimal doping of semiconductor devices.

97 MATHEMATICS AND COMPUTING↗

Code for the manuscript titled "Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces"

We would like to disclose two scripts, written in Jupyter notebook, in which we implement the "Blackout Diffusion Process" described in the manuscript "Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces" (LA-UR-23-20509), to be submitted to the International Conference of Machine Learning (ICML). The abstract of the manuscript is append below. == Typical generative diffusion models rely on a Gaussian diffusion process for training the backward transformations, which can then be used generate samples from Gaussian noise. However, real world data often takes place in discrete-state spaces, which occur in many scientific applications. Here, we develop a theoretical formulation for arbitrary discrete-state Markov processes in the forward diffusion process. We relate the theory to the existing continuous-state Gaussian diffusion and identify the corresponding reverse-time stochastic process and score function in the continuous-time setting, and the reverse-time mapping in the discrete-time setting. As an example of this framework, we introduce "Blackout Diffusion", which learns to produce samples from an empty image instead of from noise. Numerical experiments on the CIFAR-10 dataset confirm the feasibility of generative diffusion modeling in a discrete space. Generalizing from specific (Gaussian) forward processes to a more general framework also sheds light on how to interpret generative diffusion models and their mathematical structure, which we comment on.

Lin, Yen Ting↗

New ASME Section III, Division 5 Creep-Fatigue Design Rules Based on EPP and SMT Approaches

The integrated elastic-perfectly plastic (EPP) and simplified model test (SMT) creep-fatigue (CF) design methodology, referred to as the EPP-SMT method, is being developed as an alternative for CF evaluation in the design of pressure boundary components for high-temperature reactors. This report reviews the conceptual basis of the EPP-SMT methodology, summarizes the SMT experimental development efforts and results, and provides the technical basis for finalizing the EPP-SMT CF design curves for Alloy 617, based on a combined experimental and numerical approach conducted in FY 2024. This report presents the effect of hold time on the CF design curves for Alloy 617 at elevated temperatures. It includes proposed EPP-SMT CF design curves and tabulated values for continuous cycling, along with the effects of maximum hold time, for the use of this EPP-SMT CF evaluation method.

36 MATERIALS SCIENCE↗

Development of a computer code for calculating the steady super/hypersonic inviscid flow around real configurations. Volume 1: Computational technique

A numerical procedure has been developed to compute the inviscid super/hypersonic flow field about complex vehicle geometries accurately and efficiently. A second order accurate finite difference scheme is used to integrate the three dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects of blunt nose entropy layers are computed in detail. Real gas effects for equilibrium air are included using curve fits of Mollier charts. Typical calculated results for shuttle orbiter, hypersonic transport, and supersonic aircraft configurations are included to demonstrate the usefulness of this tool.

Marconi, F.↗

Derivation and evaluation of an approximate analysis for three-dimensional viscous subsonic flow with large secondary velocities

An approximate analysis is presented for calculating three-dimensional, low Mach number, laminar viscous flows in curved passages with large secondary flows and corner boundary layers. The analysis is based on the decomposition of the overall velocity field into inviscid and viscous components with the overall velocity being determined from superposition. An incompressible vorticity transport equation is used to estimate inviscid secondary flow velocities to be used as corrections to the potential flow velocity field. A parabolized streamwise momentum equation coupled to an adiabatic energy equation and global continuity equation is used to obtain an approximate viscous correction to the pressure and longitudinal velocity fields. A collateral flow assumption is invoked to estimate the viscous correction to the transverse velocity fields. The approximate analysis is solved numerically using an implicit ADI solution for the viscous pressure and velocity fields. An iterative ADI procedure is used to solve for the inviscid secondary vorticity and velocity fields. This method was applied to computing the flow within a turbine vane passage with inlet flow conditions of M = 0.1 and M = 0.25, Re = 1000 and adiabatic walls, and for a constant radius curved rectangular duct with R/D = 12 and 14 and with inlet flow conditions of M = 0.1, Re = 1000, and adiabatic walls.

Anderson, O. L.↗

Experiences performing conceptual design optimization of transport aircraft

Optimum Preliminary Design of Transports (OPDOT) is a computer program developed at NASA Langley Research Center for evaluating the impact of new technologies upon transport aircraft. For example, it provides the capability to look at configurations which have been resized to take advantage of active controls and provide and indication of economic sensitivity to its use. Although this tool returns a conceptual design configuration as its output, it does not have the accuracy, in absolute terms, to yield satisfactory point designs for immediate use by aircraft manufacturers. However, the relative accuracy of comparing OPDOT-generated configurations while varying technological assumptions has been demonstrated to be highly reliable. Hence, OPDOT is a useful tool for ascertaining the synergistic benefits of active controls, composite structures, improved engine efficiencies and other advanced technological developments. The approach used by OPDOT is a direct numerical optimization of an economic performance index. A set of independent design variables is iterated, given a set of design constants and data. The design variables include wing geometry, tail geometry, fuselage size, and engine size. This iteration continues until the optimum performance index is found which satisfies all the constraint functions. The analyst interacts with OPDOT by varying the input parameters to either the constraint functions or the design constants. Note that the optimization of aircraft geometry parameters is equivalent to finding the ideal aircraft size, but with more degrees of freedom than classical design procedures will allow.

Arbuckle, P. D.↗

Computational aerodynamics - The next generation

Continued advances in the various elements that comprise the field of computational fluid dynamics (CFD) are promoting a radically different approach to the aerodynamic design and analysis of aerospace vehicles and systems. The elements of CFD generally include numerical algorithm development, transition and turbulence modeling, surface modeling, and grid generation, scientific visualization and validation methodologies. This paper discusses the research progress and prospects for the future in each of these elements within NASA's CFD and Experimental Validation Program. The applicability of computational methods for the purposes of understanding complex flow phenomena, exploring aerodynamic concepts, and providing vehicle-design input is also addressed.

Hessenius, Kristin A.↗

Navier-Stokes relaxation to sinh-Poisson states at finite Reynolds numbers

A mathematical framework is proposed in which it seems possible to justify the computationally-observed relaxation of a two-dimensional Navier-Stokes fluid to a 'most probable', or maximum entropy, state. The relaxation occurs at large but finite Reynolds numbers, and involves substantial decay of higher-order ideal invariants such as enstrophy. A two-fluid formulation, involving interpenetrating positive and negative vorticity fluxes (continuous and square integrable) is developed, and is shown to be intimately related to the passive scalar decay problem. Increasing interpenetration of the two fluids corresponds to the decay of vorticity flux due to viscosity. It is demonstrated numerically that, in two dimensions, passive scalars decay rapidly, relative to mean-square vorticity (enstrophy). This observation provides a basis for assigning initial data to the two-fluid field variables.

Montgomery, David↗