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,333 records · Page 74

Adaptive explicit and implicit finite element methods for transient thermal analysis

The application of adaptive finite element methods to the solution of transient heat conduction problems in two dimensions is investigated. The computational domain is represented by an unstructured assembly of linear triangular elements and the mesh adaptation is achieved by local regeneration of the grid, using an error estimation procedure coupled to an automatic triangular mesh generator. Two alternative solution procedures are considered. In the first procedure, the solution is advanced by explicit timestepping, with domain decomposition being used to improve the computational efficiency of the method. In the second procedure, an algorithm for constructing continuous lines which pass only once through each node of the mesh is employed. The lines are used as the basis of a fully implicit method, in which the equation system is solved by line relaxation using a block tridiagonal equation solver. The numerical performance of the two procedures is compared for the analysis of a problem involving a moving heat source applied to a convectively cooled cylindrical leading edge.

Probert, E. J.↗

An efficient formulation of the coupled finite element-integral equation technique for solving large 3D scattering problems

It is often desirable to calculate the electromagnetic fields inside and about a complicated system of scattering bodies, as well as in their far-field region. The finite element method (FE) is well suited to solving the interior problem, but the domain has to be limited to a manageable size. At the truncation of the FE mesh one can either impose approximate (absorbing) boundary conditions or set up an integral equation (IE) for the fields scattered from the bodies. The latter approach is preferable since it results in higher accuracy. Hence, the two techniques can be successfully combined by introducing a surface that encloses the scatterers, applying a FE model to the inner volume and setting up an IE for the tangential fields components on the surface. Here the continuity of the tangential fields is used bo obtain a consistent solution. A few coupled FE-IE methods have recently appeared in the literature. The approach presented here has the advantage of using edge-based finite elements, a type of finite elements with degrees of freedom associated with edges of the mesh. Because of their properties, they are better suited than the conventional node based elements to represent electromagnetic fields, particularly when inhomogeneous regions are modeled, since the node based elements impose an unnatural continuity of all field components across boundaries of mesh elements. Additionally, our approach is well suited to handle large size problems and lends itself to code parallelization. We will discuss the salient features that make our approach very efficient from the standpoint of numerical computation, and the fields and RCS of a few objects are illustrated as examples.

Cwik, T.↗

Comparative study of energetic heliospheric electrons at wind and geotail

The current fleet of both heliospheric and magnetospheric spacecraft provides a unique opportunity of studying spatio-temporal plasma phenomena. Among the numerous topics that can be addressed by such a fleet, figures the 3D study of energetic (greater than 20 keV) solar electron) solar electron events. WIND and GEOTAIL are particulary interesting for such analysis: more than two years of data; up to 3/4 continuous days spent by GEOTAIL in the solar wind; and comparable experiments. During solar electron events, rapid electron flux changes can be associated with these structures. Sometimes they are seen at both spacecraft with a time difference corresponding to the convection time, but sometimes they are observed at one spacecraft but not at the other. Several preliminary conclusions can be drawn from such events: (1) there is direct evidence that, within distances lower than 1.5 x 10(exp 6) km at 1 A (approximately 7000 km at the sun), there are significant spatial variations of (a) the magnetic field line connections to the source region and probably of (b) the particle propagation features; (2) magnetic field lines inside and outside particle propagation structures can keep distinct access to the source region for more than one hour; (3) flux waves, of probable local origin, were observed; and (4) within one hour, spatial flux discrepancies, separated by magnetic structures, may remain or disappear.

Chaizy, P. A.↗

Video Based Sensor for Tracking 3-Dimensional Targets

Video-Based Sensor for Tracking 3-Dimensional Targets The National Aeronautics and Space Administration's (NASAs) Marshall Space Flight Center (MSFC) has been developing and testing video-based sensors for automated spacecraft guidance for several years, and the next generation of video sensor will have tracking rates up to 100 Hz and will be able to track multiple reflectors and targets. The Video Guidance Sensor (VGS) developed over the past several years has performed well in testing and met the objective of being used as the terminal guidance sensor for an automated rendezvous and capture system. The first VGS was successfully tested in closed-loop 3-degree-of-freedom (3- DOF) tests in 1989 and then in 6-DOF open-loop tests in 1992 and closed-loop tests in 1993-4. Development and testing continued, and in 1995 approval was given to test the VGS in an experiment on the Space Shuttle. The VGS flew in 1997 and in 1998, performing well for both flights. During the development and testing before, during, and after the flight experiments, numerous areas for improvement were found. The VGS was developed with a sensor head and an electronics box, connected by cables. The VGS was used in conjunction with a target that had wavelength-filtered retro-reflectors in a specific pattern, The sensor head contained the laser diodes, video camera, and heaters and coolers. The electronics box contained a frame grabber, image processor, the electronics to control the components in the sensor head, the communications electronics, and the power supply. The system works by sequentially firing two different wavelengths of laser diodes at the target and processing the two images. Since the target only reflects one wavelength, it shows up well in one image and not at all in the other. Because the target's dimensions are known, the relative positions and attitudes of the target and the sensor can be computed from the spots reflected from the target. The system was designed to work from I meter out to I 10 meters. The VGS was mounted on the Space Shuttle while its target was mounted on the Spartan free-flyer carried on the same Shuttle flight. The VGS tracked the Spartan at ranges up to 170m, and the VGS range data very closely matched the range data from the Hand-Held Laser- Rangefinder used by the astronauts on board the Shuttle. While the VGS was designed primarily as a terminal guidance sensor for an automated spacecraft, it could be applied to other uses. It could be used as an alignment aid for an operator of a remote system (giving position and attitude feedback data, as well as a camera view of the target), as a feedback system for a robotic arm, or for automated vehicle guidance. The next generation VGS, with its higher tracking rates, smaller size, and lower power could be used in more places than the original VGS, and by using LED's instead of laser diodes, the system would be eye-safe at any range. Other possible uses include tracking 3-dimensional objects with retro-reflectors mounted at various locations or motion analysis by placing several retro-reflectors on the moving object and tracking them at high speeds. There are few sensors capable of performing tasks similar to those the VGS can perform, and the next generation VGS will be even more capable than the original. Some of this work is previously presented in the papers.

Howard, R. T.↗

Solidification Using a Baffle in Sealed Ampoules (SUBSA)

Solidification Using a Baffle in Sealed Ampoules (SUBSA) will be the first materials science experiment conducted in the Microgravity Science Glovebox (MSG) Facility at the International Space Station (ISS) Alpha. The launch is schedule for May 31, 2002. Using the specially developed furnace, 10 Te and Zn-doped single crystals of InSb will be directionally solidified in microgravity. A key goal of the SUBSA investigation is to (i) clarify the origin of the melt motion in space laboratories and (ii) to reduce the magnitude of the melt motion to the point that it does not interfere with the transport phenomena. These goals will be accomplished through a special ampoule and furnace design. A disk-shaped baffle, positioned close to the freezing front, is used to reduce melt motion. Furthermore, the solidification will be visualized by using a transparent furnace, with a video camera, continuously sending images to the earth. This allows detection of bubbles and melt de-wetting that could cause surface tension driven convection. In preparation for the space experiments, 30 ground-based experiments were conducted. The results of ground based tests and numerical modeling will be presented. Based on numerical modeling, 12 mm 1D silica ampoules were selected. The small diameter ampoule favors closer placement of the baffle to the interface, without excessive radial segregation caused by forced convection while providing more damping of natural convection. The parts in the silica ampoule include 2 carbon springs made by Energy Science Laboratories, Inc., a pyrocarbon-coated graphite cylinder, pyrocarbon-coated graphite a baffle with the shaft and the InSb charge with the seed crystal grown by W.A. Bonner of Crystallod Inc.

Marin, C.↗

A Global, Multi-Waveband Model for the Zodiacal Cloud

This recently completed three-year project was undertaken by the PI at the University of Florida, NASA Goddard and JPL, and by the Co-I and Collaborator at the University of Florida. The funding was used to support a continuation of research conducted at the University of Florida over the last decade which focuses on the dynamics of dust particles in the interplanetary environment. The main objectives of this proposal were: To produce improved dynamical models of the zodiacal cloud by performing numerical simulations of the orbital evolution of asteroidal and cometary dust particles. To provide visualizations of the results using our visualization software package, SIMUL, simulating the viewing geometries of IRAS and COBE and comparing the model results with archived data. To use the results to provide a more accurate model of the brightness distribution of the zodiacal cloud than existing empirical models. In addition, our dynamical approach can provide insight into fundamental properties of the cloud, including but not limited to the total mass and surface area of dust, the size-frequency distribution of dust, and the relative contributions of asteroidal and cometary material. The model can also be used to provide constraints on trace signals from other sources, such as dust associated with the "Plutinos" , objects captured in the 2:3 resonance with Neptune.

Grogan, Keith↗

Data Assimilation of Photosynthetic Light-use Efficiency using Multi-angular Satellite Data: II Model Implementation and Validation

Spatially explicit and temporally continuous estimates of photosynthesis will be of great importance for increasing our understanding of and ultimately closing the terrestrial carbon cycle. Current capabilities to model photosynthesis, however, are limited by accurate enough representations of the complexity of the underlying biochemical processes and the numerous environmental constraints imposed upon plant primary production. A potentially powerful alternative to model photosynthesis through these indirect observations is the use of multi-angular satellite data to infer light-use efficiency (ε) directly from spectral reflectance properties in connection with canopy shadow fractions. Hall et al. (this issue) introduced a new approach for predicting gross ecosystem production that would allow the use of such observations in a data assimilation mode to obtain spatially explicit variations in ε from infrequent polar-orbiting satellite observations, while meteorological data are used to account for the more dynamic responses of ε to variations in environmental conditions caused by changes in weather and illumination. In this second part of the study we implement and validate the approach of Hall et al. (this issue) across an ecologically diverse array of eight flux-tower sites in North America using data acquired from the Compact High Resolution Imaging Spectroradiometer (CHRIS) and eddy-flux observations. Our results show significantly enhanced estimates of ε and therefore cumulative gross ecosystem production (GEP) over the course of one year at all examined sites. We also demonstrate that ε is greatly heterogeneous even across small study areas. Data assimilation and direct inference of GEP from space using a new, proposed sensor could therefore be a significant step towards closing the terrestrial carbon cycle.

Hilker, Thomas↗

Towards Quantifying and Correlating RTV Intumescence to Entry-Relevant Conditions

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also adopted to bond instrumentation plugs of temperature and pressure sensors into heatshield’s tiles. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Experimental data has also shown heating rate-dependent swelling (or intumescence) and shrinking of RTV, these volume changes may cause roughness-induced boundary layer transition. Outgassing and surface oxide formation of RTV upon decomposition also occur, which can be sources of contamination of heat shield sensors. This combination of detrimental effects motivates a critical need to develop a high-fidelity model for RTV ablation. As a first critical step in RTV modeling, a comprehensive material properties database for RTV was collected that account for the lack of data in properties such as pyrolysis char yield, shape change, virgin and char porosity. The next step is to quantify the volume change of RTV as a function of temperature. While this was done in previous campaigns, limited data was collected which showed a possible heating rate dependent intumescent behavior. To further understand the intumescence phenomenon of RTV, numerous dedicated experiments were performed at Beamline 8.3.2 of the Advanced Light Source, where RTV samples were heated while X-ray scans were taken in situ. Micro-Computed Tomography scans were taken in situ for low heating rates (≤ 60 °C/min), whereas continuous radiography scans were taken for higher heating rates (≤ 1500 ºC/min). Unconstrained samples and samples constrained in different holders (quartz, graphite, FiberForm and PICA) were tested to determine the change in intumescence due to confinement. Using deep learning techniques, the tomographies and radiographies are segmented to obtain volume change, porosity, pore size distributions and connectivity. From the results obtained, modifications to the previous RTV intumescence model are proposed, along with preliminary comparisons to testing RTV at high-enthalpy facilities.

RTV↗

Machine Learning Based Crater Detection for Terrain Relative Navigation

As Lunar exploration continues to become more commonplace, reliable methods of precise Terrain Relative Navigation (TRN) are needed. While there are many TRN techniques available, one that has received increased interest in the past few years is that of crater based navigation. Crater based navigation has numerous benefits, including being a human recognizable feature (important for crewed missions), as well as the fact that craters are often possible hazards that need to be detected and avoided. The use of crater based navigation has been limited however. This has been due to the difficulty of running such algorithms on board a spacecraft, as well as the difficulty in procuring large amounts of the required training data. This paper presents a new rendering tool for generating large amounts of high quality training data. It then looks at two recently developed machine learning techniques for crater detection and crater identification in real-time on near-future space hardware.

computer vision↗

A Scenario-Based Approach to Assess Continuity Gaps in Earth Observations

Decision analysis processes outlined in systems engineering references, such as the Systems Engineering Body of Knowledge [1] or the NASA Systems Engineering Handbook [2], recommend following a series of steps to support decision-making for engineering applications. These steps are consistent across the literature and typically involve defining objectives, defining relevant criteria against which candidate alternatives can be assessed, selecting an evaluation method, assessing alternatives, and making a recommendation based on this assessment. These approaches are well-suited to assess systems for which sets of common objectives and constraints can be identified. A vast body of literature describes their application to numerous engineering problems, and many methods have been developed since the middle of the 20th century to support the assessment of candidate alternatives for such systems [3, 4, 5, 6]. These approaches, however, show some limitations when decisions pertain to systems of systems. These collections of individual systems typically do not share a common set of objectives and constraints, tend to be highly complex, and the definition of assessable, high-value candidate alternatives poses a challenge due to interdependencies. The formulation, development, operation, and funding of Earth observing missions are such that gaps may occur between missions, therefore impacting the continuity of measurements. These Earth observing mission architectures constitute systems of systems for which a common set of objectives and constraints can be challenging to define. When tasked with assessing gaps that may occur for a series of spaceborne missions that address similar science parameters, the authors therefore proposed to depart from the traditional objective- and criteria-based approach and instead adopted a scenario-based approach. This approach does not require a prioritized set of common objectives; rather, it assesses the impact of possible decisions on the system-of-systems and characterizes its possible future states. It also offers the flexibility to adjust assumptions, update inputs, and refine supporting models over time, while enabling the rapid, early identification of challenges and key decision points. By simulating the impact of potential decisions on the entirety of the system-of-systems, the approach enables the synchronization of multiple decisions that typically occur at the level of the individual system. This paper provides an overview of the approach that was developed to assess continuity gaps for a series of operational and planned spaceborne Earth observing missions, a description of its application to the problem at hand, and a discussion of currently known limitations.

scenario earth observations continuity gaps↗

A Scenario-Based Approach to Assess Continuity Gaps in Earth Observations

Decision analysis processes outlined in systems engineering references, such as the Systems Engineering Body of Knowledge [1] or the NASA Systems Engineering Handbook [2], recommend following a series of steps to support decision-making for engineering applications. These steps are consistent across the literature and typically involve defining objectives, defining relevant criteria against which candidate alternatives can be assessed, selecting an evaluation method, assessing alternatives, and making a recommendation based on this assessment. These approaches are well-suited to assess systems for which sets of common objectives and constraints can be identified. A vast body of literature describes their application to numerous engineering problems, and many methods have been developed since the middle of the 20th century to support the assessment of candidate alternatives for such systems [3, 4, 5, 6]. These approaches, however, show some limitations when decisions pertain to systems of systems. These collections of individual systems typically do not share a common set of objectives and constraints, tend to be highly complex, and the definition of assessable, high-value candidate alternatives poses a challenge due to interdependencies. The formulation, development, operation, and funding of Earth observing missions are such that gaps may occur between missions, therefore impacting the continuity of measurements. These Earth observing mission architectures constitute systems of systems for which a common set of objectives and constraints can be challenging to define. When tasked with assessing gaps that may occur for a series of spaceborne missions that address similar science parameters, the authors therefore proposed to depart from the traditional objective- and criteria-based approach and instead adopted a scenario-based approach. This approach does not require a prioritized set of common objectives; rather, it assesses the impact of possible decisions on the system-of-systems and characterizes its possible future states. It also offers the flexibility to adjust assumptions, update inputs, and refine supporting models over time, while enabling the rapid, early identification of challenges and key decision points. By simulating the impact of potential decisions on the entirety of the system-of-systems, the approach enables the synchronization of multiple decisions that typically occur at the level of the individual system. This paper provides an overview of the approach that was developed to assess continuity gaps for a series of operational and planned spaceborne Earth observing missions, a description of its application to the problem at hand, and a discussion of currently known limitations.

scenario earth observations continuity gaps↗

Understanding and Developing Dual-Wavelength Olefin Metathesis Polymerizations for the Rapid, Continuous Additive Manufacturing of High-Performance Thermosets

This report details the goals, approach, results, and future needs of a three-year Laboratory Directed Research and Development project funded by Sandia’s Materials Science Research Foundation. In this project, we established basic principles and mechanistic understanding of orthogonal or sequential chemistries for fabricating multi-material or robust materials. We developed numerous additive manufacturing, or 3D-printing, techniques including volumetric solid-state printing, dual-wavelength printing, and lithographic regulation of polymer crystallinity. We leveraged these new printing motifs to expand upon what is possible with manufacturing to give bespoke optical and mechanical performance. Moreover, we discovered new techniques for depolymerization enabling recovery of pristine carbon-fiber and high value electronics under moderate conditions.

36 MATERIALS SCIENCE↗

Numerical Predictions of Sonic Boom Signatures for a Straight Line Segmented Leading Edge Model

A sonic boom wind tunnel test was conducted on a straight-line segmented leading edge (SLSLE) model in the NASA Langley 4- by 4- Foot Unitary Plan Wind Tunnel (UPWT). The purpose of the test was to determine whether accurate sonic boom measurements could be obtained while continuously moving the SLSLE model past a conical pressure probe. Sonic boom signatures were also obtained using the conventional move-pause data acquisition method for comparison. The continuous data acquisition approach allows for accurate signatures approximately 15 times faster than a move-pause technique. These successful results provide an incentive for future testing with greatly increased efficiency using the continuous model translation technique with the single probe to measure sonic boom signatures. Two widely used NASA codes, USM3D (Navier-Stokes) and CART3D-AERO (Euler, adjoint-based adaptive mesh), were used to compute off-body sonic boom pressure signatures of the SLSLE model at several different altitudes below the model at Mach 2.0. The computed pressure signatures compared well with wind tunnel data. The effect of the different altitude for signature extraction was evaluated by extrapolating the near field signatures to the ground and comparing pressure signatures and sonic boom loudness levels.

Elmiligui, Alaa A.↗

On the S-matrix of Ising field theory in two dimensions

We explore the analytic structure of the non-perturbative S-matrix in arguably the simplest family of massive non-integrable quantum field theories: the Ising field theory (IFT) in two dimensions, which may be viewed as the Ising CFT deformed by its two relevant operators, or equivalently, the scaling limit of the Ising model in a magnetic field. Our strategy is that of collider physics: we employ Hamiltonian truncation method (TFFSA) to extract the scattering phase of the lightest particles in the elastic regime, and combine it with S-matrix bootstrap methods based on unitarity and analyticity assumptions to determine the analytic continuation of the 2 → 2 S-matrix element to the complex s-plane. Focusing primarily on the “high temperature” regime in which the IFT interpolates between that of a weakly coupled massive fermion and the E 8 affine Toda theory, we will numerically determine 3-particle amplitudes, follow the evolution of poles and certain resonances of the S-matrix, and exclude the possibility of unknown wide resonances up to reasonably high energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Continuous conditional generative adversarial networks for data-driven solutions of poroelasticity with heterogeneous material properties

Machine learning-based data-driven modeling can allow computationally efficient time-dependent solutions of PDEs, such as those that describe subsurface multiphysical problems. In this work, our previous approach (Kadeethum et al., 2021d) of conditional generative adversarial networks (cGAN) developed for the solution of steady-state problems involving highly heterogeneous material properties is extended to time-dependent problems by adopting the concept of continuous cGAN (CcGAN). The CcGAN that can condition continuous variables is developed to incorporate the time domain through either element-wise addition or conditional batch normalization. Moreover, this framework can handle training data that contain different timestamps and then predict timestamps that do not exist in the training data. As a numerical example, the transient response of the coupled poroelastic process is studied in two different permeability fields: Zinn & Harvey transformation and a bimodal transformation. The proposed CcGAN uses heterogeneous permeability fields as input parameters while pressure and displacement fields over time are model output. Our results show that the model provides sufficient accuracy with computational speed-up. This robust framework will enable us to perform real-time reservoir management and robust uncertainty quantification in poroelastic problems.

97 MATHEMATICS AND COMPUTING↗

Effects of grain size and porosity on cladding failure in high-burnup UO 2 : A sensitivity and uncertainty study

Isotopic taggants are being studied to aid in the provenance assessment of nuclear materials. However, these taggants must be selected such that they do not adversely affect fuel performance during normal operation or accident scenarios. Taggants are known to affect the fuel’s grain size and porosity. In the work described in this paper, the BISON fuel performance code was used to assess the potential effects of taggants (i.e., grain size and porosity) on fuel rod behavior and cladding failure during a high-burnup, large-break loss-of-coolant accident. Here, 281 individual fuel rods from the same reactor core were modeled for a sensitivity study, a parametric study, and uncertainty quantification. The cladding failure predictions often exhibited stochastic behavior. After additional study, it was found that the cladding failure model is highly sensitive to residual error inherent to numerical approximation solvers. Some strategies to mitigate this sensitivity are discussed. The study found no relationship between known taggant effects and cladding failure status. However, taggants were found to affect the time and location of failure in certain rods. In conclusion, future work to continue investigating and validating these findings is briefly discussed.

Doped UO 2↗

Signature of hadron-quark crossover in binary-neutron-star mergers

Here, we study observational signatures of the hadron-quark crossover in binary-neutron-star mergers by numerical-relativity simulations with various mass configurations. We employ two equations of state (EoSs) for matter consistent with inference from the observational data. In the crossover scenario the EoS is softened in a density realized in binary-neutron-star mergers and is smoothly continued to quark matter. In the phase transition scenario without crossover, the EoS remains stiff and a first-order phase transition takes place in a density out of reach of mergers. A GW170817-like system forms a remnant massive neutron star in both scenarios, and it collapses into a black hole only in the crossover scenario due to the softening while gravitational-wave emission is strong. This difference is clearly reflected in the sudden shutdown of gravitational waves. For a given EoS, the lifetime of the merger remnant is determined primarily by the total mass of the system. Identifying these features in a variety of future events with the next generation of ground-based gravitational-wave detectors will enable us to clarify details of hadron-quark transition. The mass of the accretion disk surrounding the remnant black hole is affected not only by the lifetime of the remnant but also by the mass ratio of the system. Electromagnetic emission associated with the disk outflow will also be useful for detailed investigation of the hadron-quark transition.

QCD phase transitions↗

Landau Modes are Eigenmodes of Stellar Systems in the Limit of Zero Collisions

Abstract We consider the spectrum of eigenmodes in a stellar system dominated by gravitational forces in the limit of zero collisions. We show analytically and numerically using the Lenard–Bernstein collision operator that the Landau modes, which are not true eigenmodes in a strictly collisionless system (except for the Jeans unstable mode), become part of the true eigenmode spectrum in the limit of zero collisions. Under these conditions, the continuous spectrum of true eigenmodes in a collisionless system, also known as the Case–van Kampen modes, is eliminated. Furthermore, because the background distribution function in a weakly collisional system can exhibit significant deviations from a Maxwellian distribution function over long times, we show that the spectrum of Landau modes can change drastically even in the presence of slight deviations from a Maxwellian, primarily through the appearance of weakly damped modes that may be otherwise heavily damped for a Maxwellian distribution. Our results provide important insights for developing statistical theories to describe thermal fluctuations in a stellar system, which are currently a subject of great interest for N -body simulations as well as observations of gravitational systems.

79 ASTRONOMY AND ASTROPHYSICS↗