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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Deformation analysis of boron/aluminum specimens by moire interferometry

Whole-field surface deformations were measured for two slotted tension specimens from multiply laminates, one with 0 deg fiber orientation in the surface ply and the other with 45 deg orientation. Macromechanical and micromechanical details were revealed using high-sensitivity moire interferometry. Although global deformations of all plies were essentially equal, numerous random or anomalous features were observed. Local deformations of adjacent 0 deg and 45 deg plies were very different, both near the slot and remote from it, requiring large interlaminar shear strains for continuity. Shear strains were concentrated in the aluminum matrix. For 45 deg plies, a major portion of the deformation was by shear; large plastic slip of matrix occurred at random locations in 45 deg plies, wherein groups of fibers slipped relative to other groups. Shear strains in the interior, between adjacent fibers, were larger than the measured surface strains.

Post, Daniel↗

NASA NTRS1990

Turbulence in astrophysics

Some of the astrophysical scenarios in which turbulence plays an important role are discussed in view of the comparative advantages of currently available models of turbulence phenomena; attention is given to a specific model that has undergone continuous refinement since 1985. The desideratum in this inquiry is a turbulence model which incorporates the best features of an a priori deterministic model, as well as a redundant set of results from full numerical simulations for a wide variety of turbulent flows; there should also be a simplification of the former, and an enlargement of the complexities of the latter.

Canuto, V. M.↗

NASA NTRS1993

A Spectral Algorithm for Envelope Reduction of Sparse Matrices

The problem of reordering a sparse symmetric matrix to reduce its envelope size is considered. A new spectral algorithm for computing an envelope-reducing reordering is obtained by associating a Laplacian matrix with the given matrix and then sorting the components of a specified eigenvector of the Laplacian. This Laplacian eigenvector solves a continuous relaxation of a discrete problem related to envelope minimization called the minimum 2-sum problem. The permutation vector computed by the spectral algorithm is a closest permutation vector to the specified Laplacian eigenvector. Numerical results show that the new reordering algorithm usually computes smaller envelope sizes than those obtained from the current standard algorithms such as Gibbs-Poole-Stockmeyer (GPS) or SPARSPAK reverse Cuthill-McKee (RCM), in some cases reducing the envelope by more than a factor of two.

Barnard, Stephen T.↗

NASA NTRS1998

Algorithms and Libraries

This exploratory study initiated our inquiry into algorithms and applications that would benefit by latency tolerant approach to algorithm building, including the construction of new algorithms where appropriate. In a multithreaded execution, when a processor reaches a point where remote memory access is necessary, the request is sent out on the network and a context--switch occurs to a new thread of computation. This effectively masks a long and unpredictable latency due to remote loads, thereby providing tolerance to remote access latency. We began to develop standards to profile various algorithm and application parameters, such as the degree of parallelism, granularity, precision, instruction set mix, interprocessor communication, latency etc. These tools will continue to develop and evolve as the Information Power Grid environment matures. To provide a richer context for this research, the project also focused on issues of fault-tolerance and computation migration of numerical algorithms and software. During the initial phase we tried to increase our understanding of the bottlenecks in single processor performance. Our work began by developing an approach for the automatic generation and optimization of numerical software for processors with deep memory hierarchies and pipelined functional units. Based on the results we achieved in this study we are planning to study other architectures of interest, including development of cost models, and developing code generators appropriate to these architectures.

Dongarra, Jack↗

NASA NTRS2001

Modeling and Analysis of Wrinkled Membranes: An Overview

Thin-film membranes are basic elements of a variety of space inflatable/deployable structures. Wrinkling degrades the performance and reliability of these membrane structures, and hence has been a topic of continued interest. Wrinkling analysis of membranes for general geometry and arbitrary boundary conditions is quite challenging. The objective of this presentation is two-fold. Firstly, the existing models of wrinkled membranes and related numerical solution methods are reviewed. The important issues to be discussed are the capability of a membrane model to characterize taut, wrinkled and slack states of membranes in a consistent and physically reasonable manner; the ability of a wrinkling analysis method to predict the formation and growth of wrinkled regions, and to determine out-of-plane deformation and wrinkled waves; the convergence of a numerical solution method for wrinkling analysis; and the compatibility of a wrinkling analysis with general-purpose finite element codes. According to this review, several opening issues in modeling and analysis of wrinkled membranes that are to be addressed in future research are summarized, The second objective of this presentation is to discuss a newly developed membrane model of two viable parameters (2-VP model) and associated parametric finite element method (PFEM) for wrinkling analysis are introduced. The innovations and advantages of the proposed membrane model and PFEM-based wrinkling analysis are: (1) Via a unified stress-strain relation; the 2-VP model treat the taut, wrinkled, and slack states of membranes consistently; (2) The PFEM-based wrinkling analysis has guaranteed convergence; (3) The 2-VP model along with PFEM is capable of predicting membrane out-of-plane deformations; and (4) The PFEM can be integrated into any existing finite element code. Preliminary numerical examples are also included in this presentation to demonstrate the 2-VP model and PFEM-based wrinkling analysis approach.

Yang, B.↗

NASA NTRS2003

VLA Imaging of Protoplanetary Environments

Our program uses the techniques of millimeter interferometry to make high resolution observations of dust continuum emission to study the structure of protoplanetary disks around nearby stars. We continue to use the Very Large Array (VLA) of the National Radio Astronomy Observatories (NRAO) to make produce high angular resolution images, and we are improving our capabilities to interpret these observations with detailed numerical models. We have also made the first millimeter interferometer observations of protoplanetary disks in the southern hemisphere, using the recently upgraded Australia Telescope Compact Array (ATCA), taking advantage of the techniques we have developed in our work on northern sources.

Wilner, David J.↗

NASA NTRS2008

LiGa(OTf)(sub 4) as an Electrolyte Salt for Li-Ion Cells

Lithium tetrakis(trifluoromethane sulfo - nato)gallate [abbreviated "LiGa(OTf)4" (wherein "OTf" signifies trifluoro - methanesulfonate)] has been found to be promising as an electrolyte salt for incorporation into both liquid and polymer electrolytes in both rechargeable and non-rechargeable lithium-ion electrochemical cells. This and other ingredients have been investigated in continuing research oriented toward im proving the performances of rechargeable lithium-ion electrochemical cells, especially at low temperatures. This research at earlier stages, and the underlying physical and chemical principles, were reported in numerous previous NASA Tech Briefs articles. As described in more detail in those articles, lithiumion cells most commonly contain nonaqueous electrolyte solutions consisting of lithium hexafluorophosphate (LiPF6) dissolved in mixtures of cyclic and linear alkyl carbonates, including ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonate (DMC), diethyl carbonate (DEC), and ethyl methyl carbonate (EMC). Although such LiPF6-based electrolyte solutions are generally highly ionically conductive and electrochemically stable, as needed for good cell performance, there is interest in identifying alternate lithium electrolyte salts that, relative to LiPF6, are more resilient at high temperature and are less expensive. Experiments have been performed on LiGa(OTf)4 as well as on several other candidate lithium salts in pursuit of this interest. As part of these experiments, LiGa(OTf)4 was synthesized by the reaction of Ga(OTf)3 with an equimolar portion of LiOTf in a solvent consisting of anhydrous acetonitrile. Evaporation of the solvent yielded LiGa(OTf)4 as a colorless crystalline solid. The LiGa(OTf)4 and the other salts were incorporated into solutions with PC and DMC. The resulting electrolyte solutions exhibited reasonably high ionic conductivities over a relatively wide temperature range down to 40 C (see figure). In cyclic voltammetry measurements, LiGa(OTf)4 and the other salts exhibited acceptably high electrochemical stability over the relatively wide potential window of 0 to 5 V versus Li+/Li. 13C nuclear-magneticresonance measurements yielded results that suggested that in comparison with the other candidate salts, LiGa(OTf)4 exhibits less ion pairing. Planned further development will include optimization of the salt and solvent contents of such electrolyte solutions and incorporation of LiGa(OTf)4 into gel and solid-state polymer electrolytes. Of the salts, LiGa(OTf)4 is expected to be especially desirable for incorporation into lithium polymer electrolytes, wherein decreased ion pairing is advantageous and the large delocalized anions can exert a plasticizing effect.

Reddy, V. Prakash↗

NASA NTRS2009

The Capillary Flow Experiments Aboard the International Space Station: Increments 9-15

This report provides a summary of the experimental, analytical, and numerical results of the Capillary Flow Experiment (CFE) performed aboard the International Space Station (ISS). The experiments were conducted in space beginning with Increment 9 through Increment 16, beginning August 2004 and ending December 2007. Both primary and extra science experiments were conducted during 19 operations performed by 7 astronauts including: M. Fincke, W. McArthur, J. Williams, S. Williams, M. Lopez-Alegria, C. Anderson, and P. Whitson. CFE consists of 6 approximately 1 to 2 kg handheld experiment units designed to investigate a selection of capillary phenomena of fundamental and applied importance, such as large length scale contact line dynamics (CFE-Contact Line), critical wetting in discontinuous structures (CFE-Vane Gap), and capillary flows and passive phase separations in complex containers (CFE-Interior Corner Flow). Highly quantitative video from the simply performed flight experiments provide data helpful in benchmarking numerical methods, confirming theoretical models, and guiding new model development. In an extensive executive summary, a brief history of the experiment is reviewed before introducing the science investigated. A selection of experimental results and comparisons with both analytic and numerical predictions is given. The subsequent chapters provide additional details of the experimental and analytical methods developed and employed. These include current presentations of the state of the data reduction which we anticipate will continue throughout the year and culminate in several more publications. An extensive appendix is used to provide support material such as an experiment history, dissemination items to date (CFE publication, etc.), detailed design drawings, and crew procedures. Despite the simple nature of the experiments and procedures, many of the experimental results may be practically employed to enhance the design of spacecraft engineering systems involving capillary interface dynamics.

Jenson, Ryan M.↗

NASA NTRS2011

Collocation and Galerkin Time-Stepping Methods

We study the numerical solutions of ordinary differential equations by one-step methods where the solution at tn is known and that at t(sub n+1) is to be calculated. The approaches employed are collocation, continuous Galerkin (CG) and discontinuous Galerkin (DG). Relations among these three approaches are established. A quadrature formula using s evaluation points is employed for the Galerkin formulations. We show that with such a quadrature, the CG method is identical to the collocation method using quadrature points as collocation points. Furthermore, if the quadrature formula is the right Radau one (including t(sub n+1)), then the DG and CG methods also become identical, and they reduce to the Radau IIA collocation method. In addition, we present a generalization of DG that yields a method identical to CG and collocation with arbitrary collocation points. Thus, the collocation, CG, and generalized DG methods are equivalent, and the latter two methods can be formulated using the differential instead of integral equation. Finally, all schemes discussed can be cast as s-stage implicit Runge-Kutta methods.

Huynh, H. T.↗

NASA NTRS2011

Inversion of Magnetic Measurements of the CHAMP Satellite Over the Pannonian Basin

The Pannonian Basin is a deep intra-continental basin that formed as part of the Alpine orogeny. In order to study the nature of the crustal basement we used the long-wavelength magnetic anomalies acquired by the CHAMP satellite. The anomalies were distributed in a spherical shell, some 107,927 data recorded between January 1 and December 31 of 2008. They covered the Pannonian Basin and its vicinity. These anomaly data were interpolated into a spherical grid of 0.5 x 0.5, at the elevation of 324 km by the Gaussian weight function. The vertical gradient of these total magnetic anomalies was also computed and mapped to the surface of a sphere at 324 km elevation. The former spherical anomaly data at 425 km altitude were downward continued to 324 km. To interpret these data at the elevation of 324 km we used an inversion method. A polygonal prism forward model was used for the inversion. The minimum problem was solved numerically by the Simplex and Simulated annealing methods; a L2 norm in the case of Gaussian distribution parameters and a L1 norm was used in the case of Laplace distribution parameters. We INTERPRET THAT the magnetic anomaly WAS produced by several sources and the effect of the sable magnetization of the exsolution of hemo-ilmenite minerals in the upper crustal metamorphic rocks.

Kis, K. I.↗

DOE OSTI2022

A flexible linear diffusion acceleration to k-eigenvalue neutron transport with SN discontinuous finite element method

In this paper, we derive a flexible linear diffusion acceleration (LDA) for k-eigenvalue neutron transport discretized with discontinuous finite element method (DFEM) and discrete ordinates(SN). This LDA is based on our two pieces of previous works: the flexible non linear diffusion acceleration (NDA) for DFEM-SN and LDA for k-eigenvalue neutron transport using pre-conditioned Jacobian-free Newton-Krylov with self-adjoint angular flux (SAAF), continuous finite element method(CFEM), and SN. We point out the differences between LDA and NDA for DFEM-SN and the difference between DFEM-SN and SAAF-CFEM-SN for LDA. Numerical tests are presented to compare the convergence behaviour of NDA and LDA. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1983

A comparison of numerical results of Arctic Sea ice modeling with satellite images

Observations made from 1972 to 1976 with the Electrically Scanning Microwave Radiometer on board the Nimbus 5 Satellite provide sequential synoptic information of the Arctic sea ice cover. This 4 year data set was used to construct a fairly continuous series of 3 day average 19-GHz passive microwave images which has become a valuable source of polar information, yielding many anticipated and unanticipated discoveries of the sea ice canopy observed in its entirety through the clouds and during the polar night. Interpretation of the passive microwave satellite data set was performed by comparing selected sequential images with correspondin microwave profiles and images acquired by the NASA CV-990 airborne laboratory and with various in situ microwave and physical observations.

Ling, C. H.↗

DOE OSTI2022

Report of the Snowmass 2021 Topical Group on Lattice Gauge Theory

Lattice gauge theory continues to be a powerful theoretical and computational approach to simulating strongly interacting quantum field theories, whose applications permeate almost all disciplines of modern-day research in High-Energy Physics. Whether it is to enable precision quark- and lepton-flavor physics, to uncover signals of new physics in nucleons and nuclei, to elucidate hadron structure and spectrum, to serve as a numerical laboratory to reach beyond the Standard Model, or to invent and improve state-of-the-art computational paradigms, the lattice-gauge-theory program is in a prime position to impact the course of developments and enhance discovery potential of a vibrant experimental program in High-Energy Physics over the coming decade. This projection is based on abundant successful results that have emerged using lattice gauge theory over the years: on continued improvement in theoretical frameworks and algorithmic suits; on the forthcoming transition into the exascale era of high-performance computing; and on a skillful, dedicated, and organized community of lattice gauge theorists in the U.S. and worldwide. The prospects of this effort in pushing the frontiers of research in High-Energy Physics have recently been studied within the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021), and the conclusions are summarized in this Topical Report.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRSDate not supplied

Turbulence, Transition, and Numerical Method Technologies

The close collaboration in the validation experimentation effort was excellent, the complementary research efforts related tofundamental numerical model development, and the activities related to effective HPC utilization on near term architectures that are coming down the pipeline were exactly what NASA should be doing. Even though the TQR panel cannot say that the success criteria were fully met, we do feel that the research funded by TTT/RCA was successful and useful in determining the next suite of CFD validation tests and the problematic flow conditions that continue to be difficult to predict using our current suite of computational tools. Given the TQR Panel’s review criteria: 1. The deliverables meet the technical requirements; 2.The appropriate technical approaches were followed in producing the deliverables; and 3. All technical caveats and concerns have been identified and the deliverables: a.Development of more accurate physics-based methods (e.g., higher moment closure); b. Large eddy simulation (LES); c. Advanced numerical methods; d. Transition prediction and modeling; e.Validation experiments; f. Multidisciplinary analysis and design (high fidelity). TQR panel feels that the technical processes were all completed and that the review criteria and all the deliverables were met. Thus, the Milestone TACP-2016-001 has been completed.

Mujeeb R Malik↗

DOE OSTI2024

Continuous polyamorphic transition in high-entropy metallic glass

Polyamorphic transition (PT) is a compelling and pivotal physical phenomenon in the field of glass and materials science. Understanding this transition is of scientific and technological significance, as it offers an important pathway for effectively tuning the structure and property of glasses. In contrast to the PT observed in conventional metallic glasses (MGs), which typically exhibit a pronounced first-order nature, herein we report a continuous PT (CPT) without first-order characteristics in high-entropy MGs (HEMGs) upon heating. This CPT behavior is featured by the continuous structural evolution at the atomic level and an increasing chemical concentration gradient with temperature, but no abrupt reduction in volume and energy. The continuous transformation is associated with the absence of local favorable structures and chemical heterogeneity caused by the high configurational entropy, which limits the distance and frequency of atomic diffusion. As a result of the CPT, numerous glass states can be generated, which provides an opportunity to understand the nature, atomic packing, formability, and properties of MGs. Moreover, this discovery highlights the implication of configurational entropy in exploring polyamorphic glasses with an identical composition but highly tunable structures and properties.

36 MATERIALS SCIENCE↗

DOE OSTI2022

IDAC-Bio, A Software for Internal Dosimetry Based on the New ICRP Biokinetic Models and Specific Absorbed Fractions

Radiation dosimetry is central to virtually all radiation safety applications, optimization, and research. It relates to various individuals and population groups and to miscellaneous exposure situations—including planned, existing, and emergency situations. The International Commission on Radiological Protection (ICRP) has developed a new computational framework for internal dose estimations. Important components are more detailed and improved anatomical models and more realistic biokinetic models than before. The ICRP is currently producing new organ dose and effective dose coefficients for occupational intakes of radionuclides (OIR) and environmental intakes of radionuclides (EIR), which supersede the earlier dose coefficients in Publication 68 and the Publication 72 series, respectively. However, the ICRP only publishes dose coefficients for a single acute intake of a radionuclide and for an integration period of 50 years for intake by adults and to age 70 years for intakes by pre-adults. The new software, IDAC-Bio, performs committed absorbed dose and effective dose calculations for a selectable intake scenario, e.g., for a continuous intake or an intake during x hours per day and y days per week, and for any selected integration time. The software uses the primary data and models of the ICRP biokinetic models and numerically solves the biokinetic model and calculates the absorbed doses to organs and tissues in the ICRP reference human phantoms. The software calculates absorbed dose using the nuclear decay data in ICRP publication 107. IDAC-Bio is a further development and an important addition to the internal dosimetry program IDAC-Dose2.1. Here, the results generated by the software were validated against published ICRP dose coefficients. The potential of the software is illustrated by dose calculations for a nuclear power plant worker who had been exposed to varying levels of 60 Co and who had undergone repeated whole-body measurements, and for a hypothetical member of the public subject to future releases of 148 Gd from neutron spallation in tungsten at the European Spallation Source.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2025

Evolutionary and functional relationships between plant and microbial C 1 metabolism in terrestrial ecosystems

One-carbon (C 1 ) metabolism, centered on the universal methyl donor S-adenosyl methionine (SAM), plays critical roles in biosynthesis, redox regulation, and stress responses across plants and microbes. A recently proposed photosynthetic C 1 pathway links SAM methyl groups directly to RuBisCO-mediated CO 2 assimilation and integrates with nitrogen and sulfur metabolism. Light-dependent SAM synthesis may regulate the methylation of biopolymers and specialized metabolites and help mitigate photorespiratory stress under elevated temperature and drought. Phylogenetic analysis of two core enzymes suggests evolutionary continuity from methylotrophic microbes to land plants, supporting microbial origins via endosymbiotic gene transfer. Beyond intracellular roles, C 1 metabolism drives biosphere–atmosphere exchange via gases such as methane, methanol, formic acid, and formaldehyde, and numerous specialized volatiles synthesized through SAM methylation. S-methylmethionine, a mobile C 1 metabolite, may mediate phloem transport of reduced sulfur, nitrogen, and methyl groups, linking above- and belowground C 1 cycling in plants. Advances in real-time gas sensing now allow the high-frequency quantification of C 1 fluxes from leaves, stems, and soils, highlighting C 1 metabolism as a critical and underrecognized component of terrestrial carbon and nutrient cycling. Given its microbial ancestry and the production of diverse volatile biosignatures, C 1 metabolism may also offer unique insights into life's origins and biosignature detection on exoplanets.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

ABSTRACT Fault stepovers are prime examples of geometric complexity in natural fault zones that may affect seismic hazard by determining whether an earthquake rupture continues propagating or abruptly stops. However, the long-term pattern of seismicity near-fault stepovers and underlying mechanisms of rupture jumping in the context of earthquake cycles are rarely studied. Leveraging a hybrid numerical scheme combining the finite element and the spectral boundary integral methods, FEBE, we carry out fully dynamic simulations of sequences of earthquakes and aseismic slip for both compressive and tensile stepovers with off-fault plasticity. We consider a rate-and-state friction law for the fault friction and pressure-sensitive Drucker–Prager plasticity for the off-fault bulk response. We observe that the accumulation of plastic deformation, an indication of off-fault damage, is significantly different in the two cases, with more plastic deformation projected in the overlapping region for the tensile stepover. The seismic pattern for a tensile stepover is more complex than for a compressive stepover, and incorporating plasticity also increases complexity, relative to the elastic case. A tensile stepover with off-fault plasticity shows rupture segmentation, temporal clustering, and frequent rupture jumping from one fault to another. These results shed light on possible mechanisms of rupture jumping in fault stepovers as well as the long-term evolution of the fault zone.

Geochemistry & Geophysics↗