Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

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 901 records · Page 50

Mauka Energy FEVER Tool Dataset

Mauka Energy’s dataset, developed under the Forestry Electric Vehicle Energy Routing (FEVER) project and funded by the U.S. Department of Energy’s Small Business Innovation Research program, is a high-resolution geospatial resource designed to support energy modeling for electric log trucks in complex forestry environments. The dataset integrates detailed spatial and road network data to enable accurate simulation of vehicle performance across varied terrain. At its core, the dataset incorporates lidar-derived elevation models, road alignments, and surface classifications from Oregon State University’s McDonald-Dunn Research Forest. These data capture fine-scale variations in slope, curvature, and surface conditions across forest road systems, allowing for vehicle-level analysis of energy consumption and recovery. The dataset also includes data collected on the surrounding public and private road networks in Benton County, Oregon, used in real-world haul routes. These connecting segments provide critical context for modeling transitions between forest operations and regional transportation infrastructure, incorporating attributes such as grade profiles, elevation change, and speed constraints. This combined dataset underpins the development of Mauka Energy’s rolldown tool, which quantifies energy use and regenerative braking potential on downhill and variable-grade segments. By leveraging high-resolution terrain and road data, the FEVER project enables more accurate assessment of electric vehicle feasibility and performance in forestry applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Improve wildfire predictability driven by extreme water cycle with interpretable physically-guided ML/AI

Focal Area(s): Predictive modeling using AI techniques and AI-derived model components; use of AI and other tools to design a prediction system comprising of a hierarchy of models (Primary); insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics or knowledge-guided AI (Secondary). Science Challenge: Wildfires modify land surface characteristics, such as vegetation composition, soil and litter carbon stocks, and surface albedo, with significant consequences for the regional carbon cycle. For example, tropical regions (i.e., African and South America) are particularly vulnerable to wildfire, account for more than 80% of the global burned area, and emit ~1.4 PgC y -1 into the atmosphere together with other dust and aerosols that strongly affect regional climate.

54 ENVIRONMENTAL SCIENCES↗

Using CESM-RESFire to understand climate–fire–ecosystem interactions and the implications for decadal climate variability

Large wildfires exert strong disturbance on regional and global climate systems and ecosystems by perturbing radiative forcing as well as the carbon and water balance between the atmosphere and land surface, while short- and long-term variations in fire weather, terrestrial ecosystems, and human activity modulate fire intensity and reshape fire regimes. The complex climate–fire–ecosystem interactions were not fully integrated in previous climate model studies, and the resulting effects on the projections of future climate change are not well understood. Here we use the fully interactive REgion-Specific ecosystem feedback Fire model (RESFire) that was developed in the Community Earth System Model (CESM) to investigate these interactions and their impacts on climate systems and fire activity. We designed two sets of decadal simulations using CESM-RESFire for present-day (2001–2010) and future (2051–2060) scenarios, respectively, and conducted a series of sensitivity experiments to assess the effects of individual feedback pathways among climate, fire, and ecosystems. Our implementation of RESFire, which includes online land–atmosphere coupling of fire emissions and fire-induced land cover change (LCC), reproduces the observed aerosol optical depth (AOD) from space-based Moderate Resolution Imaging Spectroradiometer (MODIS) satellite products and ground-based AErosol RObotic NETwork (AERONET) data; it agrees well with carbon budget benchmarks from previous studies. We estimate the global averaged net radiative effect of both fire aerosols and fire-induced LCC at –0.59±0.52 W m –2 , which is dominated by fire aerosol–cloud interactions (–0.82±0.19 W m –2 ), in the present-day scenario under climatological conditions of the 2000s. The fire-related net cooling effect increases by ~170 % to –1.60±0.27 W m –2 in the 2050s under the conditions of the Representative Concentration Pathway 4.5 (RCP4.5) scenario. Such considerably enhanced radiative effect is attributed to the largely increased global burned area (+19 %) and fire carbon emissions (+100 %) from the 2000s to the 2050s driven by climate change. The net ecosystem exchange (NEE) of carbon between the land and atmosphere components in the simulations increases by 33 % accordingly, implying that biomass burning is an increasing carbon source at short-term timescales in the future. High-latitude regions with prevalent peatlands would be more vulnerable to increased fire threats due to climate change, and the increase in fire aerosols could counter the projected decrease in anthropogenic aerosols due to air pollution control policies in many regions. We also evaluate two distinct feedback mechanisms that are associated with fire aerosols and fire-induced LCC, respectively. On a global scale, the first mechanism imposes positive feedbacks to fire activity through enhanced droughts with suppressed precipitation by fire aerosol–cloud interactions, while the second one manifests as negative feedbacks due to reduced fuel loads by fire consumption and post-fire tree mortality and recovery processes. These two feedback pathways with opposite effects compete at regional to global scales and increase the complexity of climate–fire–ecosystem interactions and their climatic impacts.

54 ENVIRONMENTAL SCIENCES↗

Three-Dimensional Navier-Stokes Simulations with Two-Equation Turbulence Models of Intersecting Shock-Waves/Turbulent Boundary Layer at Mach 8.3

An investigation of the numerical simulation with two-equation turbulence models of a three-dimensional hypersonic intersecting (SWTBL) shock-wave/turbulent boundary layer interaction flow is presented. The flows are solved with an efficient implicit upwind flux-difference split Reynolds-averaged Navier-Stokes code. Numerical results are compared with experimental data for a flow at Mach 8.28 and Reynolds number 5.3x10(exp 6) with crossing shock-waves and expansion fans generated by two lateral 15 fins located on top of a cold-wall plate. This experiment belongs to the hypersonic database for modeling validation. Simulations show the development of two primary counter-rotating cross-flow vortices and secondary turbulent structures under the main vortices and in each corner singularity inside the turbulent boundary layer. A significant loss of total pressure is produced by the complex interaction between the main vortices and the uplifted jet stream of the boundary layer. The overall agreement between computational and experimental data is generally good. The turbulence modeling corrections show improvements in the predictions of surface heat transfer distribution and an increase in the strength of the cross-flow vortices. Accurate predictions of the outflow flowfield is found to require accurate modeling of the laminar/turbulent boundary layers on the fin walls.

Bardina, J. E.↗

Flexible Skins Containing Integrated Sensors and Circuitry

Artificial sensor skins modeled partly in imitation of biological sensor skins are undergoing development. These sensor skins comprise flexible polymer substrates that contain and/or support dense one- and two-dimensional arrays of microscopic sensors and associated microelectronic circuits. They afford multiple tactile sensing modalities for measuring physical phenomena that can include contact forces; hardnesses, temperatures, and thermal conductivities of objects with which they are in contact; and pressures, shear stresses, and flow velocities in fluids. The sensor skins are mechanically robust, and, because of their flexibility, they can be readily attached to curved and possibly moving and flexing surfaces of robots, wind-tunnel models, and other objects that one might seek to equip for tactile sensing. Because of the diversity of actual and potential sensor-skin design criteria and designs and the complexity of the fabrication processes needed to realize the designs, it is not possible to describe the sensor-skin concept in detail within this article.

Liu, Chang↗

Response Surface Modeling Tolerance and Inference Error Risk Specifications: Proposed Industry Standards

This paper reviews the derivation of an equation for scaling response surface modeling experiments. The equation represents the smallest number of data points required to fit a linear regression polynomial so as to achieve certain specified model adequacy criteria. Specific criteria are proposed which simplify an otherwise rather complex equation, generating a practical rule of thumb for the minimum volume of data required to adequately fit a polynomial with a specified number of terms in the model. This equation and the simplified rule of thumb it produces can be applied to minimize the cost of wind tunnel testing.

DeLoach, Richard↗

Deep potential molecular dynamics simulations of ion-enhanced etching of silicon by atomic chlorine

The continued development of plasma-assisted processing techniques requires a fundamental understanding of plasma-surface interactions. Molecular dynamics (MD) simulations have been employed to complement experimental studies and better understand the properties of such systems. Recently, machine learning (ML) methods have enabled the development of ab initio-based interatomic potentials, which can be generalized to complex combinations of multiple atom types. In this work, we use ML potentials developed using the Deep Potential Molecular Dynamics (DeepMD) framework to provide a model of ion-enhanced etching of Si by Cl atoms. We demonstrate the importance of proper selection of the training data set to the accuracy of the DeepMD model and compare our results to MD results using empirical potentials, as well as to experimental measurements. Exposure of undoped Si at 300 K to thermal Cl atoms yields a steady-state Cl coverage of 1.25 monolayers, which is slightly lower than the value obtained in previous experimental studies. Predictions of Si etch yields by simultaneous Cl atom and Ar + ion impacts as a function of ion energy, neutral to ion flux ratio, and angle of incidence of the ions are in reasonably good agreement with classical MD results and experimental measurements. Finally, etch yields and SiCl x mixed layer thicknesses during simultaneous bombardment of the Si(100) surface by Cl atoms and Cl + ions are in good agreement with experimental data. In conclusion, the present work is a necessary condition for the extension of the DeepMD procedure to more complex systems of interest in plasma-surface interactions.

Artificial neural networks↗

Stability of mineral-organic matter associations under varying biogeochemical conditions

It is well-recognized and documented in the literature that a considerable fraction of the soil organic matter (SOM) is associated with Fe-oxide minerals. Recent studies have also showed that the persistence of various classes of SOM may depend on soil Fe mineralogy. Other laboratory studies with model mineral-associated organic matter (MAOM) moieties, which were primarily focused on ferrihydrite, either coated or coprecipitated with natural organic matter (NOM), have suggested that the mineral phase transformation and adsorption-desorption of SOM is rather complex but appears to be largely controlled by the C/Fe ratio of the MAOM moiety and geochemical conditions. It is also clear from multiple recent studies that the nature or type of SOM (microbial vs. plant-derived) adsorbed on mineral surfaces and the rate and extent of SOM biodegradation depend on the type of mineral, and/or metal ions, that are associated with the SOM. Modeling studies that are primarily centered on achieving a better understanding of the effect of mineral type on SOM interactions with minerals and degradation are also reported. These studies have improved the understanding of the nature and stability of MAOM in soils of varying types and under different conditions, positioning researchers to better predict the properties and stabilities of organo-mineral associations in soils.

54 ENVIRONMENTAL SCIENCES↗

Catalytic resonance theory: Circumfluence of programmable catalytic loops

Chemical reactions on heterogeneous catalyst surfaces exhibit complex networks of elementary reactions with multiple pathways to fluid phase products, sometimes leading to surface reaction loops consisting of a closed cycle reaction pathway. While conventional catalysts at steady state exhibit zero net flux in either direction around a catalytic loop, the loop turnover frequency of three-species surface loops was evaluated in this work via microkinetic modeling to assess the reaction loop behavior resulting from a catalytic surface oscillating between two or more surface energy states. For dynamic heterogeneous catalysts undergoing applied oscillations of surface energy (i.e., programs), surface reaction loops of three species were shown to exhibit non-zero net flow of molecules around the loop, with the extent of loop turnover frequency varying with the applied frequency and amplitude of a sinusoidal or square wave oscillation. Alternatively, some dynamic surface reactions exhibited oscillations only between two surface species or resulted in surfaces covered by a single species at all times. Here, more complex behavior was observed for dynamic surface programs of three distinct electronic states, with the temporal order of applied states controlling the direction of flow of molecules within a three-molecule surface loop. Catalytic loops have the potential to limit overall catalytic reaction rates and use energy in programmable catalysts, while some applications could purposely impose non-zero loop turnover frequency for improved surface reaction control.

Catalysis↗

Infrared flow visualization of V/STOL aircraft

During low altitude operations of vertical/short takeoff and landing (V/STOVL) aircraft complex aeropropulsive flowfields are produced by the strong interactions between the engine exhaust jet or propeller wash and the runway and aircraft surfaces. Measurement of these unsteady 3 dimensional flowfields has not proven practical with conventional measurement techniques. It has been demonstrated that modern infrared imaging systems provide a practical means of obtaining flow visualization images of these complex aeropropulsive flowfields. Infrared images obtained with a dual wavelength scanning system are presented which provide detailed information on the flowfields about the NASA Ames XV-15 tiltrotor aircraft, the E-7 STOVL research model and the NASA Ames YAV-8B Harrier aircraft. It was found that the use of a dual wavelength infrared system provides valuable information about both the engine exhaust flowfield and the resulting heating of aircraft and runway surfaces.

Birckelbaw, Larry D.↗

Cryptate binding energies towards high throughput chelator design: metadynamics ensembles with cluster–continuum solvation

A tiered forcefield/semiempirical/meta-GGA pipeline together with a thermodynamic scheme designed with error cancellation in mind was developed to calculate binding energies of [2.2.2] cryptate complexes of mono- and divalent cations. Stable complexes of Na, K, Rb, Ca, Zn and Pb were generated, revealing consistent cation–N lengths but highly variable cation–O lengths and an amine stacking mechanism potentially augmenting the cation size selectivity. Metadynamics, used for searching the high-dimensional potential energy surface, together with a cluster–continuum model for affordable – yet accurate – solvation modeling, enabled the discovery of more stable geometries than those previously reported. Similar solvation energy curve shapes for lone vs. coordinated ions enabled rapid solvation convergence via the cancellation of errors stemming from finite cluster sizes. In conclusion, an R 2 of 0.850 vs. experimental aqueous binding energies was obtained, validating this scheme as the backbone of a high-throughput workflow for chelator design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A generalized behavioral model for rotating short period comets with spectral orbital elements and axial orientation

A generalized model for short period comets is developed which integrates in a fairly rigorous manner the isolation history of regions on rotating comets with specified axial orientation and the complex feedback processes involving heat, gas and dust transport, dust mantle development and coma opacity. Attention is focused on development, reconfiguration and partial or complete launching of dust mantles and the reciprocal effects of these three processes on ice surface temperature and gas and dust production. The dust mantle controls the H2O flux not only by its effect on the temperature at the ice interface but (dominantly) by its dynamic stability which strongly influences vapor diffusivity. The model includes the effects of latitude, rotation and spin axis orientation are included and applied to an initially homogeneous sphere of H2O ice and silicate using the orbital parameters of comet Encke. Numerous variations of the model, using combinations of grain size distribution, dust-to-ice ratio, latitude and spin axis orientation, are presented and discussed. Resulted for a similar nonrotating, constant Sun orientation models are also included.

Fanale, F. P.↗

Modeling Deformable Linear Objects for Autonomous Robotic Outfitting of Lunar Surface Systems

The construction of in-space facilities is a significant capability for the establishment of a long-term human presence in space. Autonomous robotic technologies are a critical tool to enabling the construction and maintenance of such permanent facilities. This paper evaluates the outfitting portion of the construction process, focusing on modeling required for robotic manipulation of cable-like objects, referred to as Deformable Linear Objects (DLOs). DLOs contain a high degree of variability, which makes accurate real-time autonomous operations a difficult task. Different modeling methods for DLOs such as discrete mass-spring systems and Cosserat elastic rod models range in problem complexity and accuracy, a trade-off that must be understood to fully realize autonomous cable routing technologies. This paper validates different DLO models through simulation and a hardware experiment, highlighting the size of the state space and accuracy of different approaches. The development of DLO manipulation models for outfitting enables an autonomous architecture for the construction of lunar surface systems.

Amy M Quartaro↗

Modeling Deformable Linear Objects for Autonomous Robotic Outfitting of Lunar Surface Systems

The construction of in-space facilities is a significant capability for the establishment of a long-term human presence in space. Autonomous robotic technologies are a critical tool to enabling the construction and maintenance of such permanent facilities. This paper evaluates the outfitting portion of the construction process, focusing on modeling required for robotic manipulation of cable-like objects, referred to as Deformable Linear Objects (DLOs). DLOs contain a high degree of variability, which makes accurate real-time autonomous operations a difficult task. Different modeling methods for DLOs such as discrete mass-spring systems and Cosserat elastic rod models range in problem complexity and accuracy, a trade-off that must be understood to fully realize autonomous cable routing technologies. This paper validates different DLO models through simulation and a hardware experiment, highlighting the size of the state space and accuracy of different approaches. The development of DLO manipulation models for outfitting enables an autonomous architecture for the construction of lunar surface systems.

Deformable Modeling↗

Control-structure-thermal interactions in analysis of lunar telescopes

The lunar telescope project was an excellent model for the CSTI study because a telescope is a very sensitive instrument, and thermal expansion or mechanical vibration of the mirror assemblies will rapidly degrade the resolution of the device. Consequently, the interactions are strongly coupled. The lunar surface experiences very large temperature variations that range from approximately -180 C to over 100 C. Although the optical assemblies of the telescopes will be well insulated, the temperature of the mirrors will inevitably fluctuate in a similar cycle, but of much smaller magnitude. In order to obtain images of high quality and clarity, allowable thermal deformations of any point on a mirror must be less than 1 micron. Initial estimates indicate that this corresponds to a temperature variation of much less than 1 deg through the thickness of the mirror. Therefore, a lunar telescope design will most probably include active thermal control, a means of controlling the shape of the mirrors, or a combination of both systems. Historically, the design of a complex vehicle was primarily a sequential process in which the basic structure was defined without concurrent detailed analyses or other subsystems. The basic configuration was then passed to the different teams responsible for each subsystem, and their task was to produce a workable solution without requiring major alterations to any principal components or subsystems. Consequently, the final design of the vehicle was not always the most efficient, owing to the fact that each subsystem design was partially constrained by the previous work. This procedure was necessary at the time because the analysis process was extremely time-consuming and had to be started over with each significant alteration of the vehicle. With recent advances in the power and capacity of small computers, and the parallel development of powerful software in structural, thermal, and control system analysis, it is now possible to produce very detailed analyses of intermediate designs in a much shorter period of time. The subsystems can thus be designed concurrently, and alterations in the overall design can be quickly adopted into each analysis; the design becomes an iterative process in which it is much easier to experiment with new ideas, configurations, and components. Concurrent engineering has the potential to produce efficient, highly capable designs because the effect of one subystem on another can be assessed in much more detail at a very early point in the program. The research program consisted of several tasks: scale a prototype telescope assembly to a 1 m aperture, develop a model of the telescope assembly by using finite element (FEM) codes that are available on site, determine structural deflections of the mirror surfaces due to the temperature variations, develop a prototype control system to maintain the proper shape of the optical elements, and most important of all, demonstrate the concurrent engineering approach with this example. In addition, the software used for the finite element models and thermal analysis was relatively new within the Program Development Office and had yet to be applied to systems this large or complex; understanding the software and modifying it for use with this project was also required. The I-DEAS software by Structural Dynamics Research Corporation (SDRC) was used to build the finite element models, and TMG developed by Maya Heat Transfer Technologies, Ltd. (which runs as an I-DEAS module) was used for the thermal model calculations. All control system development was accomplished with MATRIX(sub X) by Integrated Systems, Inc.

Thompson, Roger C.↗

Numerical framework for integrated additive manufacturing-compression molding (AM-CM) of thermoplastic composites

Additive manufacturing-compression molding (AM-CM) has emerged as a transformative technology in advanced composite manufacturing. Additive manufacturing (AM) offers high design flexibility and the ability to produce complex geometries with precisely aligned fibers in the preferred orientation. Compression molding (CM) enhances composite materials by providing excellent dimensional stability, reduced porosity, high production rates, and a smooth surface finish. Despite these advantages, extensive integrated analysis is required to optimize processing conditions for improved fiber orientation distribution (FOD) and porosity control. Here, this study develops a comprehensive numerical model to simulate the AM-CM manufacturing process. The model isolates the effects of both the AM and CM phases while also capturing their integration. Additionally, it accounts for heat transfer, temperature-dependent viscosity, and fiber orientation in the extruded fiber-filled polymer, accurately representing material behavior during processing. This approach enables the analysis of interactions between deposited beads of complex strand shapes and their interface regions after full compression. Moreover, the model predicts key parameters such as polymer flowability, fiber orientation, and temperature evolution in AM-CM parts. By optimizing processing conditions, it facilitates a controlled and predictable microstructure.

36 MATERIALS SCIENCE↗

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Precision Characterization of Seismicity from the 2022 Hunga Tonga-Hunga Ha'apai Volcanic Eruption

The earthquake swarm accompanying the January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption includes a large number of posteruptive moderate-magnitude seismic events and presents a unique opportunity to use remote monitoring methods to characterize and compare seismic activity with other historical caldera-forming eruptions. We compute improved epicentroid locations, magnitudes, and regional moment tensors of seismic events from this earthquake swarm using regional to teleseismic surface-wave cross correlation and waveform modeling. Precise relative locations of 91 seismic events derived from 59,047 intermediate-period Rayleigh- and Love-wave cross-correlation measurements collapse into a small area surrounding the volcano and exhibit a southeastern time-dependent migration. Regional moment tensors and observed waveforms indicate that these events have a similar mechanism and exhibit a strong positive compensated linear vector dipole component. Precise relative magnitudes agree with regional moment tensor moment magnitude ($M_w$) estimates while also showing that event sizes and frequency increase during the days after the eruption followed by a period of several weeks of less frequent seismicity of a similar size. The combined information from visual observation and early geologic models indicate that the observed seismicity may be the result of a complex series of events that occurred after the explosive eruption on 15 January, possibly involving rapid resupply of the magma chamber shortly after the eruption and additional faulting and instability in the following weeks. In addition, we identify and characterize an $M_w$ 4.5 event five days before the paroxysmal explosion on 15 January, indicating that additional seismic events preceding the main eruption could have been identified with improved local monitoring. As a result, our analysis of the HTHH eruption sequence demonstrates the value of potentially utilizing teleseismic surface-wave cross correlation and waveform modeling methods to assist in the detailed analysis of remote volcanic eruption sequences.

58 GEOSCIENCES↗