SEARCH · Engineering Papers

Results for “convergence rate”

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,189 records · Page 66

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Pantropical tropical tree xylem hydraulic functional traits, 1991 - 2013

This data package contains xylem hydraulic properties data, both newly extracted from publications plus tropics-specific data extracted from the "Xylem Functional Traits Database" (Choat et al., 2012; Gleason et al., 2016), which is available on TRY (www.try-db.org). Choat, B., Jansen, S., et al: Global convergence in the vulnerability of forests to drought, Nature, doi: 10.1038/nature11688, 2012. 2012. Gleason, S. M., Westoby, M., Jansen, S., Choat, B., et al.: Weak tradeoff between xylem safety and xylem-specific hydraulic efficiency across the world's woody plant species, New Phytologist, 209, 123-136, 2016. Measurements include xylem hydraulic conductivity per unit sapwood (Ks) and per unit leaf (Kl) area, water potential at 12%, 50%, and 88% loss of conductivity (P12, P50, and P88, respectively), leaf mass per area, wood density, and maximum rate of photosynthesis. Data span both natural tropical forests and plantations and rainfall gradients from deciduous dry tropical forests to everwet tropical forests, and come from countries all across the pantropics. Data is included in the attached zip file as three CSV files.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1998

Observational and Modeling Studies of Radiative, Chemical, and Dynamical Interactions in the Earth''s Atmosphere

A 3-dimensional model was developed to support mechanistic studies. The model solves the global primitive equations in isentropic coordinates, which directly characterize diabatic processes forcing the Brewer-Dobson circulation of the middle atmosphere. It's numerical formulation is based on Hough harmonics, which partition horizontal motion into its rotational and divergent components. These computational features, along with others, enable 3D integrations to be performed practically on RISC computer architecture, on which they can be iterated to support mechanistic studies. The model conserves potential vorticity quite accurately under adiabatic conditions. Forced by observed tropospheric structure, in which integrations are anchored, the model generates a diabatic circulation that is consistent with satellite observations of tracer behavior and diabatic cooling rates. The model includes a basic but fairly complete treatment of gas-phase photochemistry that represents some 20 chemical species and 50 governing reactions with diurnally-varying shortwave absorption. The model thus provides a reliable framework to study transport and underlying diabatic processes, which can then be compared against chemical and dynamical structure observed and in GCM integrations. Integrations with the Langley GCM were performed to diagnose feedback between simulated convection and the tropical circulation. These were studied in relation to tropospheric properties controlling moisture convergence and environmental conditions supporting deep convection, for comparison against mechanistic integrations of wave CISK that successfully reproduce the Madden-Julian Oscillation (MJO) of the tropical circulation. These comparisons were aimed at identifying and ultimately improving aspects of the convective simulation, with the objective of recovering a successful simulation of the MJO in the Langley GCM, behavior that should be important to budgets of upper-tropospheric water vapor and chemical species.

Salby, Murry↗

NASA NTRS1991

Navier-Stokes simulation of the supersonic combustion flowfield in a ram accelerator

A computational study of the ram accelerator, a ramjet-in-tube device for accelerating projectiles to ultrahigh velocities, is presented. The analysis is performed using a fully implicit TVD scheme that efficiently solves the Reynolds-averaged Navier-Stokes equations and the species continuity equations associated with a finite rate combustion model. The present results indicate that viscous effects are of primary importance in all the cases studied, shock-induced combustion always started in the boundary layer. The effects of Mach number, mixture composition, pressure and tubulence are investigated for various configurations. Two types of combustion processes, one stable and the other unstable, were observed depending on the inflow conditions. The possibility of stabilizing the detonation wave by means of a backward facing step is also investigated. Two numerical techniques were tested: vector extrapolation, to accelerate convergence, and a diagonal formulation that eliminates the expense of inverting large block matrices which arise in chemically reacting flows.

Yungster, Shaye↗

DOE OSTI2021

The simulation of 3D hypervelocity spallation using a hydrocode PAGOSA with FLIP+MPM

n this work, a hydrocode PAGOSA with FLIP+MPM is presented and exercised to investigate the fracture in ductile material. The merit of PAGOSA with FLIP+MPM to solve the advection problem is first illustrated by a solid piston periodically moving in a sealed tube with air. Furthermore, the ability of PAGOSA with FLIP+MPM to capture the fracture in material is shown by a simple stretching fracture in ductile material. In both of two benchmark problems, the PAGOSA results and analytical solutions are also presented for comparison. Then PAGOSA with FLIP+MPM is used to model complex spall in ductile material, which is a challenging problem in engineering applications. The convergences of PAGOSA with FLIP+MPM—based on both the mesh size and the marker density—are investigated by monitoring free surface velocity. To further show the ability of PAGOSA with FLIP+MPM to predict fracture in ductile material, the numerical results are compared with the experimental results and other published numerical results. Moreover, the effect of the spall parameter in PAGOSA with FLIP+MPM on the numerical simulation is also analyzed by investigating the free surface velocity. Finally, the effects of the peak compressive stress, the tensile strain rate, and the loading rate on the spallation are further investigated. The numerical results show that PAGOSA with FLIP+MPM can improve PAGOSA's performance when applied to predicting fracture in ductile materials, and is robust enough to accurately predict spall fracture in ductile materials.

42 ENGINEERING↗

NASA NTRS1976

Emission object in Orion

The note draws attention to a structural feature in the Orion nebula, consisting of a parachute-shaped object converging on 64 Orionis, noted on H alpha photographs of the Orion nebula taken with the 76-cm Allegheny refractor. The same feature shows up clearly on a photograph taken in 1909 with the Mount Wilson 152-cm reflector, and on a photograph taken with the Palomar 5-m reflector. Measurements made on these three plates indicate a possible expansion of the object of about 3.6 seconds of arc per century, corresponding to a tangential velocity of about 90 to 100 km/sec. If the expansion continued at a linear rate, an explosive event may have occurred approximately 900 years ago.

Feibelman, W. A.↗

NASA NTRS1984

Yaw Attitude Estimation for the Tracking and Data Relay Satellite System

The Tracking and Data Relay Satellite System (TDRSS) uses a groundbased attitude determination algorithm to open loop point the satellite's high data rate antennas. The spacecraft is able to measure its pitch and roll attitude, but its yaw attitude is periodically unobservable. The ground software uses a state-space estimator, an adaptation of a Luenberger observer, to predict the spacecraft yaw angle during these unobservable periods. It contains states associated with the roll/yaw dynamics and the on-board control law. The accuracy is limited by the modeling fidelity of the disturbance torques acting on the spacecraft. After initial operating problems were cleared up, the operation of the estimator has converged to predicted performance.

Staich, S.↗

DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1975

Solution of an optimal control lifting body entry problem by an improved method of perturbation functions

This paper presents a solution to a complex lifting reentry three-degree-of-freedom problem by using the calculus of variations to minimize the integral of the sum of the aerodynamics loads and heat rate input to the vehicle. The entry problem considered does not have state and/or control constraints along the trajectory. The calculus of variations method applied to this problem gives rise to a set of necessary conditions which are used to formulate a two point boundary value (TPBV) problem. This TPBV problem is then numerically solved by an improved method of perturbation functions (IMPF) using several starting co-state vectors. These vectors were chosen so that each one had a larger norm with respect to show how the envelope of convergence is significantly increased using this method and cases are presented to point this out.

Garcia, F., Jr.↗

NASA NTRS2010

Relating Convective and Stratiform Rain to Latent Heating

The relationship among surface rainfall, its intensity, and its associated stratiform amount is established by examining observed precipitation data from the Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR). The results show that for moderate-high stratiform fractions, rain probabilities are strongly skewed toward light rain intensities. For convective-type rain, the peak probability of occurrence shifts to higher intensities but is still significantly skewed toward weaker rain rates. The main differences between the distributions for oceanic and continental rain are for heavily convective rain. The peak occurrence, as well as the tail of the distribution containing the extreme events, is shifted to higher intensities for continental rain. For rainy areas sampled at 0.58 horizontal resolution, the occurrence of conditional rain rates over 100 mm/day is significantly higher over land. Distributions of rain intensity versus stratiform fraction for simulated precipitation data obtained from cloud-resolving model (CRM) simulations are quite similar to those from the satellite, providing a basis for mapping simulated cloud quantities to the satellite observations. An improved convective-stratiform heating (CSH) algorithm is developed based on two sources of information: gridded rainfall quantities (i.e., the conditional intensity and the stratiform fraction) observed from the TRMM PR and synthetic cloud process data (i.e., latent heating, eddy heat flux convergence, and radiative heating/cooling) obtained from CRM simulations of convective cloud systems. The new CSH algorithm-derived heating has a noticeably different heating structure over both ocean and land regions compared to the previous CSH algorithm. Major differences between the new and old algorithms include a significant increase in the amount of low- and midlevel heating, a downward emphasis in the level of maximum cloud heating by about 1 km, and a larger variance between land and ocean in the new CSH algorithm.

Tao, Wei-Kuo↗

NASA NTRS2008

The Interplanetary Overlay Networking Protocol Accelerator

A document describes the Interplanetary Overlay Networking Protocol Accelerator (IONAC) an electronic apparatus, now under development, for relaying data at high rates in spacecraft and interplanetary radio-communication systems utilizing a delay-tolerant networking protocol. The protocol includes provisions for transmission and reception of data in bundles (essentially, messages), transfer of custody of a bundle to a recipient relay station at each step of a relay, and return receipts. Because of limitations on energy resources available for such relays, data rates attainable in a conventional software implementation of the protocol are lower than those needed, at any given reasonable energy-consumption rate. Therefore, a main goal in developing the IONAC is to reduce the energy consumption by an order of magnitude and the data-throughput capability by two orders of magnitude. The IONAC prototype is a field-programmable gate array that serves as a reconfigurable hybrid (hardware/ firmware) system for implementation of the protocol. The prototype can decode 108,000 bundles per second and encode 100,000 bundles per second. It includes a bundle-cache static randomaccess memory that enables maintenance of a throughput of 2.7Gb/s, and an Ethernet convergence layer that supports a duplex throughput of 1Gb/s.

Pang, Jackson↗

DOE OSTI2020

Electronic couplings for singlet fission: Orbital choice and extrapolation to the complete basis set limit

Abstract For the search for promising singlet fission candidates, the calculation of the effective electronic coupling, which is required to estimate the singlet fission rate between the initially excited state (S 0 S 1 ) and the multiexcitonic state ( 1 TT, two triplets on neighboring molecules, coupled into a singlet), should be sufficiently reliable and fast enough to explore the configuration space. We propose here to modify the calculation of the effective electronic coupling using a nonorthogonal configuration interaction approach by: (a) using only one set of orbitals, optimized for the triplet state of the molecules, to describe all molecular electronic states, and (b) only taking the leading configurations into consideration. Furthermore, we also studied the basis set convergence of the electronic coupling, and we found, by comparison to the complete basis set limit obtained using the cc‐pV n Z series of basis sets, that both the aug‐cc‐pVDZ and 6–311++G** basis sets are a good compromise between accuracy and computational feasibility. The proposed approach enables future work on larger clusters of molecules than dimers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Understanding Precipitation Bias Sensitivities in E3SM-Multi-Scale Modeling Framework From a Dilution Framework

We investigate a set of Energy Exascale Earth System Model Multi-scale modeling framework (MMF) (E3SM-MMF) simulations that vary the dimensionality and momentum transport configurations of the embedded cloud-resolving models (CRMs), including unusually ambitious 3D configurations. Issues endemic to all MMF simulations include too much Intertropical Convergence Zone rainfall and too little over the Amazon. Systematic MMF improvements include more on-equatorial rainfall across the Warm Pool. Interesting sensitivities to the CRM domain are found in the regional time-mean precipitation pattern over the tropics. The 2D E3SM-MMF produces an unrealistically rainy region over the northwestern tropical Pacific; this is reduced in computationally ambitious 3D configurations that use 1,024 embedded CRM grid columns per host cell. Trajectory analysis indicates that these regional improvements are associated with desirably fewer tropical cyclones and less extreme precipitation rates. To understand why and how the representation of precipitation improved in 3D, we propose a framework that dilution is stronger in 3D. This viewpoint is supported by multiple indirect lines of evidence, including a delayed moisture-precipitation pickup, smaller precipitation efficiency, and amplified convective mass flux profiles and more high clouds. We also demonstrate that the effects of varying embedded CRM dimensionality and momentum transport on precipitation can be identified during the first few simulated days, providing an opportunity for rapid model tuning without high computational cost. Meanwhile the results imply that other less computationally intensive ways to enhance dilution within MMF CRMs may also be strategic tuning targets.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Shedding far-ultraviolet light on the donor star and evolutionary state of the neutron-star LMXB Swift J1858.6−0814

ABSTRACT The evolution of accreting X-ray binary systems is closely coupled to the properties of their donor stars. Consequently, we can constrain the evolutionary track a system is by establishing the nature of its donor. Here, we present far-ultraviolet (far-UV) spectroscopy of the transient neutron-star low-mass X-ray binary J1858 in different accretion states (low-hard, high-hard, and soft). All of these spectra exhibit anomalous N v, C iv, Si iv, and He ii lines, suggesting that its donor star has undergone CNO processing. We also determine the donor’s effective temperature, Td ≃ 5700 K, and radius, Rd ≃ 1.7 R⊙, based on photometric observations obtained during quiescence. Lastly, we leverage the transient nature of the system to set an upper limit of $\dot{M}_{\rm acc} \lesssim 10^{-8.5}~{\rm M}_{\odot }~\mathrm{ yr}^{-1}$ on the present-day mass-transfer rate. Combining these with the orbital period of the system, Porb = 21.3 h, we search for viable evolution paths. The initial donor masses in the allowed solutions span the range 1 M⊙ ≲ Md,i ≲ 3.5 M⊙. All but the lowest masses in this range are consistent with the strong CNO-processing signature in the UV line ratios. The present-day donor mass in the permitted tracks are 0.5 M⊙ ≲ Md,obs ≲ 1.3 M⊙, higher than suggested by eclipse modelling. Since Porb is close to the so-called bifurcation period, both converging and diverging binary tracks are permitted. If the former is confirmed, J1858 will end its life as an ultracompact system with a substellar donor.

Castro Segura, N. (ORCID:0000000258700443)↗

NASA NTRS1994

Separation of crack extension modes in composite delamination problems

This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.

Beuth, Jack L.↗

DOE OSTI2025

Numerical analysis of air dehumidification through electrospray-enhanced vortical flow cyclone separator

The dehumidification of air is a highly energy-intensive process, typically requiring vapor compression systems to cool the air significantly below its dew point temperature. Here, this study numerically investigates air dehumidification achieved by electrostatic and dielectrophoresis forces in a converging duct and centrifugal principles in cyclone separators. An electrospray was employed to energize water droplets, which then captured water vapor via dielectrophoresis. A novelty of the present study is the inclusion of a cyclone separator positioned downstream of the electrospray. This combination increased the dehumidification efficiency and controlled the flow of the condensate. By integrating electrospray technology and applying an electric field, this study examines the impact of particle separation efficiency, pressure drops, and cyclone performance on air dehumidification. The control parameters for air dehumidification were investigated. Three cyclones were modelled for airflow rates of 5 cubic feet per minute [cfm] (0.0024 m 3 /s), 38 cfm (0.0179 m 3 /s), and 200 cfm (0.094 m 3 /s). The results indicated air dehumidification was up to 5.7 % under environmental conditions, similar to those used to validate the CFD model. The degree of air dehumidification was primarily influenced by the electric charge acquired by the water droplets, their flight trajectory, and the vorticity intensity within the cyclone. Air dehumidification could increase by up to 13 % by increasing the number of droplets injected into the air stream, the flight time, and the electrical charge. The breakdown electric field strength threshold near the electrospray was a limiting factor for further enhancing dehumidification performance. This study examines the optimization of an electro-assisted dehumidification system that includes cyclone-based droplet separation through numerical analysis. The results provide insights into the design of energy-efficient air treatment technologies.

42 ENGINEERING↗

NASA NTRS1934

Impact buckling of thin bars in the elastic range for any end condition

Following a qualitative discussion of the complicated process involved in a short-period, longitudinal force applied to an originally not quite straight bar, the actual process is substituted by an idealized process for the purpose of analytical treatment. The simplifications are: the assumption of an infinitely high rate of propagation of the elastic longitudinal waves in the bar, limitation to slender bars, disregard of material damping and of rotatory inertia, the assumption of consistently small elastic deformations, the assumption of cross-sectional dimensions constant along the bar axis, the assumption of a shock-load constant in time, and the assumption of eccentricities on one plane. Then follow the mathematical principles for resolving the differential equation of the simplified problem, particularly the developability of arbitrary functions with steady first and second and intermittently steady third and fourth derivatives into one convergent series, according to the natural functions of the homogeneous differential equation.

Taub, Josef↗

NASA NTRS1978

The behavior of a convergent plate boundary - Crustal deformation in the South Kanto district, Japan

The northwesternmost part of the Sagami trough, a part of the Philippine Sea-Eurasian plate boundary, was ruptured during the great South Kanto earthquake in 1923. Very extensive and frequent geodetic measurements of crustal deformation have been made in the South Kanto district since the 1890's, and these constitute the most complete data set on crustal movements in the world. These data were reanalyzed and interpreted and according to our interpretation indicate the following sequence of events. The coseismic movements were due to oblique thrust and right lateral slip of about 8 m on a fault outcropping at the base of the Sagami trough. This was followed by postseismic deformation resulting from reversed afterslip of 20-60 cm that occurred at an exponentially decaying rate in time. The interseismic deformation is produced by steady subduction at a rate of about 1.8 cm/yr. During subduction the top 10-15 km of the plate boundary is apparently locked, while deeper parts slip aseismically at an irregular rate. No significant precursory deformation was observed. The recurrence time for 1923 type earthquakes is 200-300 years. The Boso and Miura peninsulas are broken into a series of fault-bound blocks that move semi-independently of the surrounding region. The subduction zone itself, where it is exposed on land, is shown to be a wide zone encompassing several faults that are active at different times.

Scholz, C. H.↗

DOE OSTI2024

Multi-scale signaling and tumor evolution in high-grade gliomas

Although genomic anomalies in glioblastoma (GBM) have been well studied for over a decade, its 5-year survival rate remains lower than 5%. We seek to expand the molecular landscape of high-grade glioma, composed of IDH-wildtype GBM and IDH-mutant grade 4 astrocytoma, by integrating proteomic, metabolomic, lipidomic, and post-translational modifications (PTMs) with genomic and transcriptomic measurements to uncover multi-scale regulatory interactions governing tumor development and evolution. Applying 14 proteogenomic and metabolomic platforms to 228 tumors (212 GBM and 16 grade 4 IDH-mutant astrocytoma), including 28 at recurrence, plus 18 normal brain samples and 14 brain metastases as comparators, reveals heterogeneous upstream alterations converging on common downstream events at the proteomic and metabolomic levels and changes in protein-protein interactions and glycosylation site occupancy at recurrence. Recurrent genetic alterations and phosphorylation events on PTPN11 map to important regulatory domains in three dimensions, suggesting a central role for PTPN11 signaling across high-grade gliomas.

60 APPLIED LIFE SCIENCES↗