Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global shape 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Comparison of Two IRI Plasmasphere Extensions with GPS-TEC Observations

Two plasmasphere extensions of the International Reference Ionosphere are made available for the users. It is aimed to estimate the effect of charged particles on technical devices in the Earth's environment and to define the ionosphere-plasmasphere operational conditions compatible with existing and future systems of radio communication, radio navigation and other relevant radio technologies in the ranges of medium and higher frequencies. The Global Core Plasma Model (GCPM-2000) of Gallagher et al. (2000) is an empirical description of thermal plasma densities in the plasmasphere, plasmapause, magnetospheric trough, and polar cap. GCPM-2000 uses the Kp index and is coupled to IRI in the transition region 500-600 km. The IZMIRAN plasmasphere model (Chasovitin et al., 1998; Gulyaeva et al., 2002) is an empirical model based on whistler and satellite observations. It presents global vertical analytical profiles of electron density smoothly fitted to IRI electron density profile at 1000 km altitude and extended towards the plasmapause (up to 36,000 km). For the smooth fitting of the two models, the shape of the IRI topside electron density profile is improved using ISIS 1, ISIS 2, and IK19 satellite inputs (Gulyaeva, 2003). The plasmasphere model depends on solar activity and magnetic activity (kp-index). The two IRI plasmasphere extensions are compared in the present study with the total electron content derived from records of Global Positioning Satellites (GPS-TEC) observations for different latitudinal, solar activity, magnetic activity, diurnal and seasonal conditions. The differences of model TEC with observed TEC in the topside ionosphere and plasmasphere are discussed.

Gulyacva, Tamara↗
DOE OSTI2024

Impact of ITCZ width on global climate: ITCZ-MIP

Abstract. The width of the Inter-Tropical Convergence Zone (ITCZ) affects tropical rainfall, Earth's albedo, large-scale circulation, and climate sensitivity. To better understand the ITCZ width and its effects on global climate, we present a protocol to force changes in ITCZ width in climate models. Starting from an aquaplanet configuration with a slab ocean, adding surface heat fluxes in the deep tropics forces the ITCZ to narrow, and subtracting them causes it to widen. The protocol successfully generates changes in ITCZ width in all four climate models used in this study. ITCZ width in each model responds linearly to forcing magnitude and sign. Comparing across the four climate models, a response to varying ITCZ width that is remarkably consistent among models is the ITCZ strength, which is greater the narrower the ITCZ. On the other hand, the effect of varying ITCZ width on climate sensitivity is divergent among our four models, varying even in sign. Results from this pilot study highlight the connections between surface fluxes, ITCZ width, and the wider climate. A comprehensive model intercomparison project (MIP) has the potential to advance understanding of both the physical processes shaping ITCZ width and its influence on remote atmospheric circulations and global climate.

Geology↗
DOE OSTI2021

Gravitational Lensing Formalism in a Curved Arc Basis: A Continuous Description of Observables and Degeneracies from the Weak to the Strong Lensing Regime

Gravitationally lensed curved arcs provide a wealth of information about the underlying lensing distortions. Extracting precise lensing information from extended sources is a key component in many studies aiming to answer fundamental questions about the universe. To maintain accuracy with increased precision, it is of vital importance to characterize and understand the impact of degeneracies inherent in lensing observables. In this work, we present a formalism to describe the gravitational lensing distortion effects resulting in curved extended arcs based on the eigenvectors and eigenvalues of the local lensing Jacobian and their directional differentials. We identify a nonlocal and nonlinear extended deflector basis that inherits these local properties. Our parameterization is tightly linked to observable features in extended sources and allows one to accurately extract the lensing information of extended images without imposing an explicit global deflector model. We quantify what degeneracies can be broken based on specific assumptions about the local lensing nature and assumed intrinsic source shape. Our formalism is applicable from the weak linear regime to the semi-linear regime and all the way up to the highly nonlinear regime of highly magnified arcs of multiple images. Furthermore, the methodology and implementation presented in this work provides a framework to assessing systematics, to guide inference efforts in the right choices in complexity based on the data at hand, and to quantify the lensing information extracted in a model-independent way

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Conservation laws and effective hadronization models

Hadronization models based on local string-breaking dynamics are typically Markovian by construction, yet the physical ensemble of final states is shaped by global constraints that couple the entire fragmentation trajectory. Recasting hadronization as a conditioned stochastic diffusion process provides a precise mathematical resolution to this tension. In particular, this language reveals explicitly that constraints stemming from conservation laws induce non-Markovian correlations between otherwise independent fragmentation steps, and that these correlations can be absorbed exactly into a renormalization of the local dynamics through a Doob $h$-transform. We develop this formalism for a $q\bar{q}$ string in the chiral limit, where the longitudinal-transverse factorization of the Lund kernel becomes exact, enabling systematic power counting and clean ultraviolet (UV)/infrared (IR) separation. The dynamics organize naturally into a tower of effective theories distinguished by the remaining string mass, spanning a UV fixed point with scale-invariant transport coefficients, an intermediate regime where transverse phase space induces controlled running, and an IR boundary layer where non-local effects enter at leading order. The tower exhibits genuine Wilsonian structure, including $β$-functions, anomalous dimensions, and systematic matching conditions. The resulting framework achieves a clean factorization of universal microscopic fragmentation dynamics from infrared constraint effects, and opens new directions for both the theoretical analysis and practical simulation of hadronization.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗
NASA NTRS1996

A Near-Global Survey of Cirrus Particle Size Using ISCCP

Cirrus is the most frequently occurring and widely distributed cloud type. The average annual frequency of occurrence for cirrus is 34% and its global coverage is about 20-30% (Warren et al. 1985). It strongly influences weather and climate processes through its effects on the radiation budget of the earth and the atmosphere (Liou 1986). Microphysics of cirrus is a critical component in understanding cloud-climate radiative interactions. For example, ice water content feedback is positive from a 1-D model study. But the feedback is substantially reduced upon the inclusion of small ice crystals (Sinha and Shine 1994). Due to the complexity caused by the non-spherical shape of ice crystals in cirrus, retrievals of cirrus properties are difficult. In recent years, advances have been made both in models and in case studies (e.g., Takano and Liou 1989, Young et al. 1994), but no global scale survey has been conducted. Similar to our previous near-global survey of droplet sizes of liquid water clouds (Han et al. 1994), a survey of cirrus ice crystal sizes is conducted over both continental and oceanic areas. We describe a method for retrieving cirrus particle size information on a near-global scale 50 deg S to 50 deg N using currently available satellite data from ISCCP. To retrieve cirrus particle size, we use a radiative transfer model that includes all major absorbing gases and cloud scattering/absorption to compute synthetic radiances as a function of satellite viewing geometry. Ice crystal shapes are assumed to be hexagonal columns and plates. The model results have been validated against clear sky observations and are consistent with the observed radiance range under cloudy conditions.

Han, Qingyuan↗
NASA NTRS1993

System identification for Space Station Freedom using observer/Kalman filter Markov parameters

The Modal Identification Experiment (MIE) is a proposed experiment to define the dynamic characteristics of Space Station Freedom. Previous studies emphasized free-decay modal identification. The feasibility of using a forced response method (Observer/Kalman Filter Identification (OKID)) is addressed. The interest in using OKID is to determine the input mode shape matrix which can be used for controller design or control-structure interaction analysis, and investigate if forced response methods may aid in separating closely spaced modes. A model of the SC-7 configuration of Space Station Freedom was excited using simulated control system thrusters to obtain acceleration output. It is shown that an 'optimum' number of outputs exists for OKID. To recover global mode shapes, a modified method called Global-Local OKID was developed. This study shows that using data from a long forced response followed by free-decay leads to the 'best' modal identification. Twelve out of the thirteen target modes were identified for such an output.

Papadopoulos, Michael↗
DOE OSTI2025

Implications of a weakening N = 126 shell closure away from stability for r -process astrophysical conditions

The formation of the third r-process abundance peak near A ∼ 195 is highly sensitive to both nuclear structure far from stability and the astrophysical conditions that produce the heaviest elements. In particular, the N = 126 shell closure plays a crucial role in shaping this peak. Experimental data hints that the shell weakens as proton number departs from Z = 82, a trend largely missed by global mass models. To investigate its impact on r-process nucleosynthesis, we employ both standard global models with strong closures and modified Duflo-Zuker (DZ) models that reproduce the weakening, combined with three sets of β − -decay rates. Strong shell closures generate sharply peaked abundances, whereas weakened closures consistent with the experimental trend produce broader, flatter patterns. Accurately reproducing the solar third peak under weakened shell strength requires sufficiently neutron-rich conditions that significant fission occurs, and slower decay rates. These results demonstrate that a weakening N = 126 shell closure away from stability imposes significant constraints on the astrophysical environments of the r-process and underscores the need for precise mass measurements and improved characterization of β − -decay properties in this region.

N = 126 closed shells↗
NASA NTRS2023

Salt Distribution from Freezing Intrusions in Ice Shells on Ocean Worlds: Application to Europa

Several icy moons and dwarf planets appear to have hosted subsurface liquid water. Liquid water intruding upwards into the icy outer shells of these worlds freezes, forming ice and (from ocean solutes) non-ice solids. Here, we model concentrated aqueous solutions below 273 K to simulate the compositional evolution of freezing spherical intrusions. Starting solutions are based on five previously reported compositional end members for Europa's ocean. For moderate-pH end members dominated by chloride, sulfate, and/or carbonate, the solids formed include Ca-, Mg-, and Na-sulfates and -carbonates, as well as Na- and K-chlorides. For silica-rich, high-pH end members, abundant amorphous silica forms with, potentially, similarly abundant NaOH and KOH. We further develop a new numerical model to compute the spatial distribution of the formed solids and residual brine as freezing progresses. If non-ice solids settle to the bottom, their deposits tend to have stacked hourglass shapes, widening each time the crystallization temperature of a new solid is reached. We discuss the applicability of this model to vertical fractures and global freezing of a subsurface ocean. These results inform (i) how compositional heterogeneities may affect the thermophysical properties of ice shells, which in turn influence convective and cryovolcanic transport, (ii) the compatibility of brine pockets with physicochemical conditions suitable for microbial life, and (iii) possible measurements of compositional heterogeneities within ice shells by spacecraft such as NASA's Europa Clipper and ESA's JUICE missions. The methodology developed here is applicable to other ice-covered ocean worlds.

Planetary structure↗
NASA NTRS2016

Contributions of the ARM Program to Radiative Transfer Modeling for Climate and Weather Applications

Accurate climate and weather simulations must account for all relevant physical processes and their complex interactions. Each of these atmospheric, ocean, and land processes must be considered on an appropriate spatial and temporal scale, which leads these simulations to require a substantial computational burden. One especially critical physical process is the flow of solar and thermal radiant energy through the atmosphere, which controls planetary heating and cooling and drives the large-scale dynamics that moves energy from the tropics toward the poles. Radiation calculations are therefore essential for climate and weather simulations, but are themselves quite complex even without considering the effects of variable and inhomogeneous clouds. Clear-sky radiative transfer calculations have to account for thousands of absorption lines due to water vapor, carbon dioxide, and other gases, which are irregularly distributed across the spectrum and have shapes dependent on pressure and temperature. The line-by-line (LBL) codes that treat these details have a far greater computational cost than can be afforded by global models. Therefore, the crucial requirement for accurate radiation calculations in climate and weather prediction models must be satisfied by fast solar and thermal radiation parameterizations with a high level of accuracy that has been demonstrated through extensive comparisons with LBL codes. See attachment for continuation.

climatology↗
NASA NTRS2018

Magnetized Disk Winds in NGC 3783

We analyze a 900 kilosecond stacked Chandra/HETG (High-Energy Transmission Grating) spectrum of NGC 3783 in the context of magnetically driven accretion-disk wind models in an effort to provide tight constraints on the global conditions of the underlying absorbers. Motivated by the earlier measurements of its absorption measure distribution (AMD) indicating X-ray-absorbing ionic columns that decrease slowly with decreasing ionization parameter, we employ 2-dimension (2-D) magnetohydrodynamic (MHD) disk wind models to describe the global outflow. We compute its photoionization structure along with the wind kinematic properties, allowing us to further calculate in a self-consistent fashion the shapes of the major X-ray absorption lines. With the wind radial density profile determined by the AMD, the profiles of the ensemble of the observed absorption features are determined by the two global parameters of the MHD wind; i.e., disk inclination theta (sub obs) and wind density normalization n (sub o). Considering the most significant absorption features in the approximately 1.8-20 angstrom range, we show that the MHD wind is best described by n(r) approximately equal to 6.9 times 10 (sup 11) (r/r (sub o)) (sup - 1.15) cubic centimeters and theta (sub obs). We argue that winds launched by X-ray heating or radiation pressure, or even MHD winds but with steeper radial density profiles, are strongly disfavored by data. Considering the properties of Fe K-band absorption features (i.e., Fe XXV and Fe XXVI), while typically prominent in the active galactic nucleus X-ray spectra, they appear to be weak in NGC 3783. For the specific parameters of our model obtained by fitting the AMD and the rest of the absorption features, these features are found to be weak, in agreement with observations.

Fukumura, Keigo↗
NASA NTRS1990

Application of mixing-controlled combustion models to gas turbine combustors

Gas emissions were studied from a staged Rich Burn/Quick-Quench Mix/Lean Burn combustor were studied under test conditions encountered in High Speed Research engines. The combustor was modeled at conditions corresponding to different engine power settings, and the effect of primary dilution airflow split on emissions, flow field, flame size and shape, and combustion intensity, as well as mixing, was investigated. A mathematical model was developed from a two-equation model of turbulence, a quasi-global kinetics mechanism for the oxidation of propane, and the Zeldovich mechanism for nitric oxide formation. A mixing-controlled combustion model was used to account for turbulent mixing effects on the chemical reaction rate. This model assumes that the chemical reaction rate is much faster than the turbulent mixing rate.

Nguyen, Hung Lee↗
NASA NTRS2021

Limits on Runoff Episode Duration for Early Mars: Integrating Lake Hydrology and Climate Models

Fluvio-lacustrine features on the martian surface attest to a climate that was radically different in the past. Since climate models have difficulty sustaining a liquid hydrosphere at the surface, multiple cycles of runoff episodes may have characterized the ancient Mars climate. A fundamental question thus remains: what was the duration of these runoff-producing episodes? Here we use morphometric measurements from newly identified coupled lake systems, containing both an open- and a closed-basin lake (n=7). We combined hydrological balances with precipitation outputs from climate models, and found that breaching runoff episodes likely lasted 102−105yr; other episodes may have been shorter but could not be longer. Runoff episode durations are model-dependent and spatiallyvariable, and no27climate model scenario can satisfy a unique duration for all coupled systems. In the near future, these quantitative constraints on early Mars lake persistence may be tested through in situ observations from Perseverance rover. Plain Language Summary Ancient river and lake beds on the martian surface point to a climate that was radically different in the past, unlike the dry, cold desert that Mars is today. It experienced significant liquid water flowing across the surface (i.e., ‘runoff’) and hosted long-lived bodies of water (lakes and seas). However, scientists that generate models to simulate the ancient Mars climate have difficulty sustaining this proposed liquid water cycle at the surface. One possible hypothesis is that multiple cycles of individual liquid water episodes may have characterized the ancient martian climate (i.e. runoff episodes). A fundamental question thus remains: what was the duration of these individual runoff-producing episodes? Here we use an approach that takes advantage of unique martian lake system, ‘coupled’ lake systems, and their shapes. We combined these measurements with precipitation results from previous climate models, and found that runoff episodes likely lasted 100-10,000 years, requiring a global climate mechanism that generated repeated cycles of this duration. Ultimately, our results can be tested through observations made by NASA’s Perseverance rover, shedding light on the uninterrupted duration of surface water on Mars and its potential to host life in the past.

Gaia Stucky de Quay↗
DOE OSTI2022

A Typology for Characterizing Human Action in MultiSector Dynamics Models

The role of individual and collective human agency is increasingly recognized as a prominent and arguably paramount determinant in shaping the behavior, trajectory, and vulnerability of multisector systems. This human influence operates at multiple scales: from short-term (hourly to daily) to long-term (annually to centennial) timescales, and from the local to the global, pushing systems towards either desirable or undesirable outcomes. However, the effort to represent human systems in multisector models has been fragmented across philosophical, methodological, and disciplinary lines. To cohere insights across diverse modeling approaches, we present a new typology for classifying how human actors are represented in the broad suite of coupled human-natural system models that are applied in MultiSector Dynamics (MSD) research. The typology conceptualizes a “sector” as a system-of-systems that includes a diverse group of human actors, defined across individual to collective social levels, involved in governing, provisioning, and utilizing products, goods, or services towards some human end. We trace the salient features of modeled representations of human systems by organizing the typology around three key questions: 1) Who are the modeled actors in MSD systems? 2) What are their modeled actions? 3) How and for what purpose are these actors and actions operationalized in a computational model? The typology is used to critically examine existing models and chart the frontier of future human systems modeling for MSD research.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Introduction to and comparison of deep learning and optimization approaches to analytical wake modeling of a tilted wind turbine

This paper introduces innovative optimization and deep learning techniques to enhance the prediction of complex wake dynamics in the downstream wind velocity of tilted wind turbines. Traditional methods for calibrating the Bastankhah wake model often lead to increased errors in wind velocity distribution due to overfitting of the local wake characteristics. To address this issue, we propose an additional global optimization step to reduce errors in wind velocity predictions with respect to various wake parameters. Despite this improvement, the Bastankhah model's axisymmetric Gaussian wake shape limits its accuracy for complex wake structures. Therefore, we also propose a deep learning approach, which demonstrates promising results by accurately modeling complex wake shapes across a broader range of tilt angles with minimal computational cost. The deep learning approach achieves near-identical predictions to high-fidelity large-eddy simulations, representing a promising advancement in wake modeling.

17 WIND ENERGY↗
NASA NTRS1985

Development of a global 3-D magnetohydrodynamic computational model for solar wind/cometary and planetary studies

Progress is reported in the development of a global, 3-D magnetohydrodynamic computational model to quantitatively describe the detailed continuum field and plasma interaction process of the solar wind with cometary and planetary magneto/ionopause shapes. Activities described include enhancements to the present Level 1 model, development of the 3-D gas dynamic flow field solvers for the Level 2 model, and continued development of the 3-D frozen magnetic field solver for the Level 2 model.

Stahara, S. S.↗
NASA NTRS1984

The effects of dipole tilt on the structure of the magnetosphere

The structure of the magnetosphere as a function of dipole tilt was studied by using a global magnetohydrodynamics (MHD) model. We show here that the cusp current sheet in the MHD model is related to the tail current sheet and that the calculated shape and positions of the tail current sheet are consistent with observations. In addition, we discuss the difference between the MHD model and the Chapman-Ferraro type model which manifests as a function of tilt angle.

Wu, C. C.↗
DOE OSTI2025

Knowledge Oriented Graph Unified Transformer (KOGUT) v0.1

KOGUT — Knowledge Oriented Graph Unified Transformer KOGUT implements the Relational Graph Transformer (RelGT) architecture for knowledge graph link prediction in biological domains, with a primary focus on microbial growth media prediction. While the original RelGT (arXiv:2505.10960) targets relational tables, time series, and multi-table databases, KOGUT adapts this architecture for heterogeneous biological knowledge graphs, providing first-in-class AI predictive models for microbial cultivation. Key Adaptations Beyond Original RelGT: - Knowledge Graph Focus: Applied to biological KGs with semantic node types (taxa, chemicals, media, phenotypes, environments) versus generic relational database tables, trained on the KG-Microbe knowledge graph (1.3M entities, 2.9M edges, 24 relation types). - Multimodal Node Encoding: Integrates node labels, categories, descriptions, and synonyms from KG metadata through learned embedding layers—adapting relational column features to graph node attributes with textual semantics. - Extended K-Hop Subgraph Strategy: Optimized neighborhood sampling (3-hop default, configurable up to 200 nodes) tuned for sparse biological networks, building on the original local-global attention framework with biological relation preservation. - Biolink Predicate Preservation: Type-specific transformations for 24 biological edge semantics (occurs_in, consumes, produces, has_phenotype, subclass_of) beyond standard relational foreign keys, enabling multi-relation link prediction. - Inductive Learning Support: Enables zero-shot predictions for novel taxa through feature-based embeddings (temperature, oxygen requirements, gram stain, cell shape), extending the original transductive relational benchmark scope to uncultured microorganisms. CheapSOTA Performance Optimizations (This Distribution): - VQ-EMA Centroid Attention: Vector quantization with exponential moving average for improved global context modeling (+5-10% MRR improvement). - HDF5 Precomputed Data Loading: One-time preprocessing of k-hop subgraphs to eliminate redundant graph traversals (2-5× training speedup). - Distributed Data Parallel Training: Multi-GPU support for scaling to larger knowledge graphs (tested on 4× NVIDIA A100 GPUs at NERSC Perlmutter). - Mixed Precision Training: Automatic mixed precision (AMP) for memory efficiency and faster training. Advantages Over Standard Knowledge Graph Embedding Models: Combines RelGT's proven multi-element tokenization (features, type, hop, structure) with graph-native biological representations, enabling interpretable link prediction across heterogeneous entities that standard embedding models (TransE, RotatE, ComplEx) and table-based transformers cannot directly model. Achieves near-perfect performance on microbial growth media prediction (MRR: 0.9966, Precision@1: 0.9932, Hit@10: 1.0000) while maintaining explainability through attention-based reasoning over biological pathways. Training Data: - KG-Microbe merged knowledge graph: 1,379,337 nodes, 2,960,472 edges - 24 biological relation types including taxonomic hierarchies, metabolic interactions, phenotype associations, and environmental relationships - Primary prediction task: Growth media suitability for microbial taxa (biolink:occurs_in, 50K edges) - Multi-relation capability: Predicts links for any of the 24 relation types, including chemical consumption/production, phenotype associations, and taxonomic classification Citation: Original RelGT Architecture: Dwivedi et al., "Relational Graph Transformer", arXiv:2505.10960, 2025 KOGUT Implementation: Knowledge Oriented Graph Unified Transformer for Microbial Growth Media Prediction Developed at Lawrence Berkeley National Laboratory (LBNL) Trained on NERSC Perlmutter supercomputer

Joachimiak, Marcin [Lawrence Berkeley National Lab↗
DOE OSTI2024

Shape restricted additive hazards models: Monotone, unimodal, and U‐shape hazard functions

We consider estimation of the semiparametric additive hazards model with an unspecified baseline hazard function where the effect of a continuous covariate has a specific shape but otherwise unspecified. Such estimation is particularly useful for a unimodal hazard function, where the hazard is monotone increasing and monotone decreasing with an unknown mode. A popular approach of the proportional hazards model is limited in such setting due to the complicated structure of the partial likelihood. Our model defines a quadratic loss function, and its simple structure allows a global Hessian matrix that does not involve parameters. Thus, once the global Hessian matrix is computed, a standard quadratic programming method can be applicable by profiling all possible locations of the mode. However, the quadratic programming method may be inefficient to handle a large global Hessian matrix in the profiling algorithm due to a large dimensionality, where the dimension of the global Hessian matrix and number of hypothetical modes are the same order as the sample size. We propose the quadratic pool adjacent violators algorithm to reduce computational costs. The proposed algorithm is extended to the model with a time‐dependent covariate with monotone or U‐shape hazard function. In simulation studies, our proposed method improves computational speed compared to the quadratic programming method, with bias and mean square error reductions. We analyze data from a recent cardiovascular study.

Mathematical & Computational Biology↗