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 631 records · Page 35

Capillary Channel Flow (CCF) EU2-02 on the International Space Station (ISS): An Experimental Investigation of Passive Bubble Separations in an Open Capillary Channel

It would be signicantly easier to design fluid systems for spacecraft if the fluid phases behaved similarly to those on earth. In this research an open 15:8 degree wedge-sectioned channel is employed to separate bubbles from a two-phase flow in a microgravity environment. The bubbles appear to rise in the channel and coalesce with the free surface in much the same way as would bubbles in a terrestrial environment, only the combined effects of surface tension, wetting, and conduit geometry replace the role of buoyancy. The host liquid is drawn along the channel by a pump and noncondensible gas bubbles are injected into it near the channel vertex at the channel inlet. Control parameters include bubble volume, bubble frequency, liquid volumetric flow rate, and channel length. The asymmetrically confined bubbles are driven in the cross-flow direction by capillary forces until they at least become inscribed within the section or until they come in contact with the free surface, whereupon they usually coalesce and leave the flow. The merging of bubbles enhances, but does not guarantee, the latter. The experiments are performed aboard the International Space Station as a subset of the Capillary Channel Flow experiments. The flight hardware is commanded remotely and continuously from ground stations during the tests and an extensive array of experiments is conducted identifying numerous bubble flow regimes and regime transitions depending on the ratio and magnitude of the gas and liquid volumetric flow rates. The breadth of the publicly available experiments is conveyed herein primarily by narrative and by regime maps, where transitions are approximated by simple expressions immediately useful for the purposes of design and deeper analysis.

fluid dynamics↗

CERES : a cryo-EM re-refinement system for continuous improvement of deposited models

The field of electron cryomicroscopy (cryo-EM) has advanced quickly in recent years as the result of numerous technological and methodological developments. This has led to an increase in the number of atomic structures determined using this method. Recently, several tools for the analysis of cryo-EM data and models have been developed within the Phenix software package, such as phenix.real_space_refine for the refinement of atomic models against real-space maps. Also, new validation metrics have been developed for low-resolution cryo-EM models. To understand the quality of deposited cryo-EM structures and how they might be improved, models deposited in the Protein Data Bank that have map resolutions of better than 5 Å were automatically re-refined using current versions of Phenix tools. The results are available on a publicly accessible web page (https://cci.lbl.gov/ceres). The implementation of a Cryo-EM Re-refinement System ( CERES ) for the improvement of models deposited in the wwPDB, and the results of the re-refinements, are described. Based on these results, contents are proposed for a `cryo-EM Table 1', which summarizes experimental details and validation metrics in a similar way to `Table 1' in crystallography. The consistent use of robust metrics for the evaluation of cryo-EM models and data should accompany every structure deposition and be reported in scientific publications.

59 BASIC BIOLOGICAL SCIENCES↗

Numerical study of hydrogen-air supersonic combustion by using elliptic and parabolized equations

The two-dimensional Navier-Stokes and species continuity equations are used to investigate supersonic chemically reacting flow problems which are related to scramjet-engine configurations. A global two-step finite-rate chemistry model is employed to represent the hydrogen-air combustion in the flow. An algebraic turbulent model is adopted for turbulent flow calculations. The explicit unsplit MacCormack finite-difference algorithm is used to develop a computer program suitable for a vector processing computer. The computer program developed is then used to integrate the system of the governing equations in time until convergence is attained. The chemistry source terms in the species continuity equations are evaluated implicitly to alleviate stiffness associated with fast chemical reactions. The problems solved by the elliptic code are re-investigated by using a set of two-dimensional parabolized Navier-Stokes and species equations. A linearized fully-coupled fully-implicit finite difference algorithm is used to develop a second computer code which solves the governing equations by marching in spce rather than time, resulting in a considerable saving in computer resources. Results obtained by using the parabolized formulation are compared with the results obtained by using the fully-elliptic equations. The comparisons indicate fairly good agreement of the results of the two formulations.

Chitsomboon, T.↗

Iterative and function-continuation Fourier deconvolution methods for enhancing mass spectrometer resolution

Mass spectrometer data in the form of ion current versus mass-to-charge ratio often include overlapping mass peaks, especially in low- and medium-resolution instruments. Numerical deconvolution of such data effectively enhances the resolution by decreasing the overlap of mass peaks. In this paper two approaches to deconvolution are presented: a function-domain iterative technique and a Fourier transform method which uses transform-domain function-continuation. Both techniques include data smoothing to reduce the sensitivity of the deconvolution to noise. The efficacy of these methods is demonstrated through application to representative mass spectrometer data and the deconvolved results are discussed and compared to data obtained from a spectrometer with sufficient resolution to achieve separation of the mass peaks studied. A case for which the deconvolution is seriously affected by Gibbs oscillations is analyzed.

Ioup, J. W.↗

A simple introduction to the SiMPL method for density-based topology optimization

We introduce a novel method for solving density-based topology optimization problems: Sigmoidal Mirror descent with a Projected Latent variable (SiMPL). The SiMPL method (pronounced as “the simple method”) optimizes a design using only first-order derivative information of the objective function. The bound constraints on the density field are enforced with the help of the (negative) Fermi–Dirac entropy, which is also used to define a non-symmetric distance function called a Bregman divergence on the set of admissible designs. This Bregman divergence leads to a simple update rule that is further simplified with the help of a so-called latent variable. Because the SiMPL method involves discretizing the latent variable, it produces a sequence of pointwise-feasible iterates, even when high-order finite elements are used in the discretization. Numerical experiments demonstrate that the method outperforms other popular first-order optimization algorithms. In conclusion, to outline the general applicability of the technique, we include examples with (self-load) compliance minimization and compliant mechanism optimization problems.

Calculus of Variations and Optimization↗

Approach for Propagating Radiometric Data Uncertainties Through NASA Ocean Color Algorithms

Spectroradiometric satellite observations of the ocean are commonly referred to as “ocean color” remote sensing. NASA has continuously collected, processed, and distributed ocean color datasets since the launch of the Sea-viewing Wide-field-of-view Sensor (SeaWiFS) in 1997. While numerous ocean color algorithms have been developed in the past two decades that derive geophysical data products from sensor-observed radiometry, few papers have clearly demonstrated how to estimate measurement uncertainty in derived data products. As the uptake of ocean color data products continues to grow with the launch of new and advanced sensors, it is critical that pixel-by-pixel data product uncertainties are estimated during routine data processing. Knowledge of uncertainties can be used when studying long-term climate records, or to assist in the development and performance appraisal of bio-optical algorithms. In this method paper we provide a comprehensive overview of how to formulate first-order first-moment (FOFM) calculus for propagating radiometric uncertainties through a selection of bio-optical models. We demonstrate FOFM uncertainty formulations for the following NASA ocean color data products: chlorophyll-a pigment concentration (Chl), the diffuse attenuation coefficient at 490 nm (K(sub d,490)), particulate organic carbon (POC), normalized fluorescent line height (nflh), and inherent optical properties (IOPs). Using a quality-controlled in situ hyperspectral remote sensing reflectance (R(sub rs,i)) dataset, we show how computationally inexpensive, yet algebraically complex, FOFM calculations may be evaluated for correctness using the more computationally expensive Monte Carlo approach. We compare bio-optical product uncertainties derived using our test R(sub rs) dataset assuming spectrally-flat, uncorrelated relative uncertainties of 1, 5, and 10%. We also consider spectrally dependent, uncorrelated relative uncertainties in R(sub rs). The importance of considering spectral covariances in R(sub rs), where practicable, in the FOFM methodology is highlighted with an example SeaWiFS image. We also present a brief case study of two POC algorithms to illustrate how FOFM formulations may be used to construct measurement uncertainty budgets for ecologically-relevant data products. Such knowledge, even if rudimentary, may provide useful information to end-users when selecting data products or when developing their own algorithms.

Bio-optics↗

Numerical generation of composite three dimensional grids by quasilinear elliptic systems

A three-dimensional grid generation technique designed to numerically construct a boundary-conforming grid within a three-dimensional region bounded by a closed surface is described. The boundary values are generated numerically by a modified elliptic system and are used to compute grid control parameters that are contained in the elliptic systems. The interior grid distribution is governed by the distribution of points on the boundary as well as by the boundary's geometric shape. The composite three-dimensional grid remains both continuous and smooth across the surface of juncture between any two adjoining subregions. The details of the method and its implementation are presented, displaying numerical results for three-dimensional composite grid for a wing-body combination and surface grids.

Thomas, P. D.↗

Orion Nebula and Planetary Nebulae

This report summarizes the research performed at Rice University related to NASA-Ames University consortium grant NCC2-5199 during the two year period 1996 September 1 through 1998 August 31. The research program, titled Orion Nebula and Planetary Nebulae, involved the analysis of Hubble Space Telescope (HST) imagery and spectroscopy of the Orion Nebula and of the planetary nebulae NGC 6818 and NGC 6210. In addition, we analyzed infrared spectra of the Orion Nebula taken with the Infrared Space Observatory (ISO) The primary collaborators at NASA-Ames were Drs. R. H. Rubin, A. G. C. M. Tielens, S. W. J. Colgan, and S. D. Lord (Tielens & Lord has since changed institutions). Other collaborators include Drs. P. G. Martin (CITA, Toronto), G. J. Ferland (U. KY), J. A. Baldwin (CTIO, Chile), J. J. Hester (ASU), D. K. Walter (SCSU), and P. Harrington (U. MD). In addition to the Principal Investigator, Professor Reginald J. Dufour of the Department of Space Physics & Astronomy, the research also involved two students, Mr. Matthew Browning and Mr. Brent Buckalew. Mr. Browning will be graduating from Rice in 1999 May with a B.A. degree in Physics and Mr. Buckalew continues as a graduate student in our department, having recently received a NASA GSRP research fellowship (sponsored by Ames). The collaboration was very productive, with two refereed papers already appearing in the literature, several others in preparation, numerous meeting presentations and two press releases. Some of our research accomplishments are highlighted below. Attached to the report are copies of the two major publications. Note that this research continues to date and related extensions of it recently has been awarded time with the HST for 1999-2000.

Dufour, Reginald J.↗

A well-posed numerical method to track isolated conformal map singularities in Hele-Shaw flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.

Baker, Gregory↗

Landau-Forbidden Quantum Criticality in Rydberg Quantum Simulators

The Landau-Ginzburg-Wilson theory of phase transitions precludes a continuous transition between two phases that spontaneously break distinct symmetries. However, quantum mechanical effects can intertwine the symmetries, giving rise to an exotic phenomenon called deconfined quantum criticality (DQC). Here, in this Letter, we study the ground state phase diagram of a one-dimensional array of individually trapped neutral atoms interacting strongly via Rydberg states, and demonstrate through extensive numerical simulations that it hosts a variety of symmetry-breaking phases and their transitions including DQC. We show how an enlarged, emergent continuous symmetry arises at the DQCs, which can be experimentally observed in the joint distribution of two distinct order parameters, obtained within measurement snapshots in the standard computational basis. Our findings highlight quantum simulators of Rydberg atoms not only as promising platforms to experimentally realize such exotic phenomena, but also as unique ones allowing access to physical properties not obtainable in traditional experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Composite Analysis for Low-Level Waste Disposal in the Hanford Site Central Plateau (FY 2020) (Rev. 2)

This document provides an updated Hanford Site composite analysis (CA). A Hanford Site CA was prepared and issued in 1998 (PNNL-11800) with an addendum provided in 2002 (PNNL-11800 Addendum 1). The CA was approved in 2002 (Frei, 2002) and has been maintained (DOE/RL-2000-29 and subsequent revisions) to support low-level waste disposal performance assessments (PAs) and disposal authorizations for facilities at the Hanford Site, including the following: Continued operation of the Environmental Restoration Disposal Facility (ERDF) and the 200 East and 200 West Low-Level Burial Grounds; Construction of the Integrated Disposal Facility (IDF); Forthcoming closure of tank residual waste systems such as Waste Management Area (WMA) C. The CA maintenance program resulted in a determination in 2015 (DOE/RL-2015-66) that the Hanford Site CA needed an update for the following reasons: While the initial Hanford Site CA has been maintained since 2001, the accumulation of basis changes reported in the annual summary reports over the succeeding 14 years merit evaluation in an updated analysis; The U.S. Department of Energy (DOE) Headquarters requested in a memorandum in 2015 (Gilbertson and Marcinowski, 2015) that “as soon as the relevant PAs are complete, the CA will be revised to account for all of the new information.” This updated Hanford Site CA accounts for the following new information: 1. Inclusion of a detailed Hanford Site baseline disposition that projects remedial activities through site closure. There have been significant changes through decision making in the Comprehensive Environmental Response, Compensation, and Liability Act of 1980 (CERCLA) process that were not available when the original CA was produced; 2. Inclusion of an updated inventory basis, new modeling capabilities, and new decisions reached in the associated record of decision (ROD) that was provided by issuance of DOE/EIS-0391 in fiscal year (FY) 2013. Development of a Hanford Site groundwater model from the baseline provided in a technical transfer of models for the Final Tank Closure and Waste Management Environmental Impact Statement (DOE/EIS-0391) commenced in FY 2014 and resulted in the plateau-to-river (P2R) groundwater flow model that is used in this updated Hanford Site CA; 3. Accounting for pump and treat systems, which were not evaluated in the initial CA. Pump and treat systems have had significant impact on groundwater flow system behavior, contaminant transport, and contaminant removal from Hanford Site groundwater. This process is accounted for in the use of the P2R model for this updated Hanford Site CA; 4. Inclusion of water-level data collected since the initial CA was completed. Data collected as the unconfined aquifer water levels continue to recede since the cessation of large liquid discharges in the late 1990s have led to marked improvement in understanding of the flow system for future conditions, particularly regarding northward flow potential in the critical Gable Gap area. The calibration in the latest version of the P2R model accounts for this information; 5. Inclusion of results from the revised ERDF PA completed in 2013 (WCH-5209) that account for updated inventory and expansion of the ERDF facility to about twice the size that was evaluated in the original CA; 6. Use of updated geoframeworks to provide the structural basis for numerical fate and transport models in the groundwater pathway of this updated CA. The geologic basis for groundwater models has continued to improve with additional data collection and interpretation with the creation and maintenance of the Hanford South Geoframework and the Central Plateau Vadose Zone Geoframework tools; 7. Incorporation of updated tank residual inventory estimates. Tank residual inventory estimates have improved with the incorporation of tank retrieval inventory data for those tanks that have completed retrieval. The CA inventory data package includes this updated information; 8. Incorporation of WMA C PA results. Two additional PAs for tank farm closure decisions are in preparation during the period required to prepare an updated Hanford Site CA: WMA C (FY 2016) and WMA A-AX (in preparation). The WMA C PA results are incorporated into the updated Hanford Site CA, and its grouted residuals model is used as the basis for a release model to account for the other tank farm systems modeling in this updated CA; 9. An update to the IDF PA was submitted in FY 2017 and has been reviewed and approved. The results of the IDF PA are incorporated into this updated CA; 10. Updated risk assessment scenarios. The risk assessment scenarios currently in use for Hanford Site CERCLA and Resource Conservation and Recovery Act of 1976 (RCRA) analyses differ from those evaluated in the initial Hanford Site CA. The representative person exposure scenario evaluated in this updated Hanford Site CA is consistent with recent PAs and CERCLA and RCRA analyses as the Hanford Site. This updated Hanford Site CA provides the following: A comparison of the updated Hanford Site CA all-pathway dose results with the performance measures during the compliance period, which is assumed to begin with site closure in calendar 2070 with the last scheduled disposal action and continue for 1,000 years postclosure (to calendar year 3069); A comparison of Hanford Site CA all-pathway dose results with the performance measures during the postcompliance period to address potential peaks beyond the compliance period. This is accomplished by evaluating dose in the period for 9,000 years following the compliance period (i.e., from calendar year 3070 to 12070).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

High-Performance Electron Sources: Numerical Methods and Beam Dynamics at the Precision Frontier

Electron sources have a wide range of applications and there are many stakeholders that express continuing need for improvements and performance enhancements. Whether we consider ultra-cold, high-brightness, high-charge or high-average current source needs, there are some common themes from the point of view of the beam dynamics involved. These are the following ones: ability to model accurately the emission processes, including the presence of often complicated cathode and other boundary surfaces with a wide range of spatial scales; ability to deal accurately and efficiently with a large number of particles interacting pair-wise, including the stochastic part of these interactions with a wide range of spatial scales; and ability to propagate the particle distributions in time, including collisions with a wide range of temporal scales. These tasks require high precision and accuracy since the goal is usually generation, transport and preservation of very high-quality beams. This grant addressed one of the rem

43 PARTICLE ACCELERATORS↗

Efficient numerical techniques for complex fluid flows

The central feature in any flow prediction method is the treatment of the coupling between the momentum and continuity equations. In natural-convection flows, the energy equation also becomes strongly coupled with the momentum equations. Because of the nonlinear nature of the coupling, these equations are solved iteratively. Iterative methods are often prone to slow convergence, divergence, and extreme sensitivity to underrelaxation factors. The aim of the present research is to develop more efficient and reliable solution schemes for the coupled flow equations. Such schemes will significantly reduce the expense of computing complex flows encountered in combustion chambers, gas turbines, heat exchangers, and other practical equipment. In the work completed so far, a technique employing norm reduction in conjunction with the successive-substitution and Newton-Raphson techniques was developed. Also, a block-correction procedure for the flow equations is currently being formulated and tested.

Patankar, Suhas V.↗

Numerical Challenges in Cassini Maneuver Operations

Launched in 1997 to observe Saturn and its system, Cassini successfully entered Saturn orbit in 2004 and impacted the planet on September 15, 2017 after 22 orbits each skimming over Saturn’s cloud tops. The Cassini mission represents the most complex gravity-assist trajectory ever flown. As such, the Flight Path Control team encountered many difficulties along the way, resulting in a continuously evolving maneuver process. In this paper, we focus on the challenges presented by the well-known singularities in the transfer problem and the unexpected numerical instabilities in state propagations through flybys, maneuver algorithm convergence issues, and orbital element targeting difficulties.

Hahn, Yungsun↗

A buffering heat exchanger/thermal energy storage system for desalination applications

A buffering heat exchanger/thermal energy storage (BHXTES) system was designed using a fast-response phase-change material/graphite foam (PCM/GF) medium for desalination and potentially other industrial applications. Unlike a store-now-and-use-later TES system, this study focused on developing a dynamically balanced, continuously used, system that optimized three functions. A lab-scale prototype was designed and fabricated with a PureTemp 151/GF storage medium. The experimental data showed excellent thermal performance repeatability indicating minimal effects of PCM expansion and contraction on the GF as well as negligible PCM redistribution effects. The experimental data also served as a validation of the numerical model of the BHXTES system, after which optimization of the BHXTES system was conducted through numerical simulations for both PureTemp 151/GF and solar salt/GF as the storage medium. The results show that, with a Therminol 55 heat transfer fluid, continuous operation of the system can be realized with a solar energy source through cycles of 8-h changing for desalination and energy storage with the supplied heat followed by 16-h discharging for desalination with the stored energy. As a result, the working fluid outlet temperatures are generally in the ranges of approximately 128–160 °C for a PureTemp 151/GF medium and 198–220 °C for a solar salt/GF medium.

Buffering↗

Closure of Internal Porosity in Continuous Casting Bloom During Heavy Reduction Process

To investigate the closure behavior of internal porosity (also referred to as internal void) in continuous casting bloom during heavy reduction (HR) and thus provide theoretical guidance for minimizing this kind of internal defect more effectively with HR, a three-dimensional (3D) mechanical model was developed based on the predicted temperature field by a 2D heat transfer model. With this 3D mechanical model, closure behaviors of internal porosity in continuous casting bloom during HR at and after the strand solidification end under different process conditions were numerically studied. It was found that the void axis length decreased significantly along the bloom thickness direction and increased slightly along the casting and bloom width directions after HR, and the influence of the initial void size on the void closure was not obvious. With a decrease of temperature difference between the bloom surface and center, HR efficiency for minimizing internal void decreased, while the required reduction force significantly increased. Compared with blooms with a uniform temperature distribution of 1100 °C, the void closure index after HR implemented at the strand solidification end was increased by ~ 25 pct. Compared with a conventional flat roll, the application of a convex roll during HR could contribute to minimizing the internal porosity more effectively and significantly enhance the reduction capacity of the withdrawal and straightening units. The void closure index of η{sub s} and η{sub v} (where η{sub s} and η{sub v} were defined based on the variation of the void aspect ratio and the void volume, respectively) was closely related to the equivalent strain (ε{sub eq}) and the hydrostatic integration parameter (Q), respectively, and two mathematical equations were derived to quantitatively describe the relationship of η{sub s} − ε{sub eq} and η{sub v} − Q.

36 MATERIALS SCIENCE↗

Long-Lived Triplet Excited State in a Heterogeneous Modified Carbon Nitride Photocatalyst

Heterogeneous carbon nitrides have numerous advantages as photocatalysts, including strong light absorption, tunable band edges, and scalability, but their performance and continued development are limited by fast charge recombination and an under-developed mechanistic understanding of photodriven interfacial electron transfer. These shortcomings are a result of complex photophysics, leading to rate asynchrony between oxidation and reduction, as well as redox processes driven out of electronic trap states rather than excited states. We show that a well-defined triplet excited state in cyanamide-modified carbon nitride is realized with appropriately sized particles. The utility of this long-lived excited state is demonstrated by its ability to drive a hydroamidation photoredox cycle. By the tuning of the particle size of CN x , the oxidation–reduction photochemistry of carbon nitride may be balanced to achieve a redox-neutral closed photocatalytic cycle. Here, these results uncover a triplet excited state chemistry for appropriately sized CNx particles that preludes a rich energy and electron transfer photochemistry for these materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PhILMs: Collaboratory on Mathematics and Physics-Informed Learning Machines for Multiscale and Multiphysics Problems

The landscape of computational science and engineering is continually evolving, with the challenge of high-dimensional regression problems standing as a significant hurdle in numerous scientific endeavors. Addressing this challenge, our research, funded by this award, has led to the development of an innovative computational framework known as Probabilistic Partition of Unity Networks (PPOU-Nets). This initiative represents a collaborative effort to harness the potential of mathematics and physics-informed machine learning in tackling multiscale and multiphysics problems prevalent in high-dimensional spaces. Through this work, we have proposed a novel methodology that seamlessly integrates adaptive dimensionality reduction and a mixture of experts model, thereby facilitating a more efficient and accurate approximation of complex functions. This research effort has not only advanced the state of computational science but also opened new avenues for exploration in quantum computing and beyond. This report outlines the motivation, methodology, key findings, and implications of our work, underscoring our contributions to the broader scientific community and the potential pathways for future research.

97 MATHEMATICS AND COMPUTING↗