Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geometric modelling”

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 1,063 records · Page 59

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗

Methodology to Generate the Decay Gamma Source for a Second Target Station Target Wedge

This report details the specifics of applying the Position-Averaged Method methodology for the calculation of the rotating lasagna-style target wedge activation and the resulting decay gamma source terms. The methodology uses MCNP® Code Version 6.2.0 with the RNUCS patch coupled with CINDER2008 from the AARE V1.0 package. The process begins with model preparations that include ensuring impurities are included in the materials in the regions of interest along with segmentation of the geometry in order to capture the decay gamma source gradients. This report also details the physics models and other MCNP®, CINDER2008, and ADVANTG 3.2.0 related options used to calculate the decay gamma source terms. ADVANTG 3.2.0 is used to generate a weight window set to aid in convergence of the neutron fluxes in the initial MCNP® protons-on-target particle transport calculation to calculate the neutron fluxes and spallation products in the regions of interest. The neutron fluxes and spallation products are then used in the CINDER2008 transmutation calculation, and the decay gamma sources are generated by the AARE_GAMMA_SOURCE_SCRIPT as a part of the AARE V1.0 package. These decay gamma sources include the energy distributions as well as the necessary distributions to sample the location of the source using the cell rejection method. The gamma sources that are produced using the methodology may be used for a wide variety of analysis involving the target wedge. This report includes the material specifications, physics options, proton source term description, geometrical configuration, and post-processing details.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Physics basis for design of 3D coils in tokamaks

Recent progress in 3D tokamak modeling is now leveraged to create a conceptual design of new external 3D field coils for the DIII-D tokamak. In this work generalized perturbed equilibrium code is used to determine optimally efficient spectrum for driving total, core, and edge neoclassical toroidal viscosity torque. These fundamental modes of 3D control are shown to have consistent outboard structures across a wide variety of plasma scenarios and machines. Given these target spectra, the currents and 3D geometry of multiple coils can be optimized to increase efficient drive for the physics of interest without undesired secondary effects. Here, this nonlinear optimization is demonstrated using the flexible optimized coils using space-curves code. The optimized coils are individually distorted in space, creating toroidal 'arrays' containing a variety of shapes that often wrap around a significant poloidal extent of the machine. Importantly, efficient coupling can be maintained even when enforcing large distances between coils and the plasma during the geometric optimization of coil designs. Finally, the physics-driven optimization presented here thus provides a practical path to utilizing coils built on the exterior of the vacuum chamber in future reactors to obtain the powerful 3D field benefits demonstrated on current machines with close, internal coils.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optical information transfer through random unknown diffusers using electronic encoding and diffractive decoding

Free-space optical information transfer through diffusive media is critical in many applications, such as biomedical devices and optical communication, but remains challenging due to random, unknown perturbations in the optical path. We demonstrate an optical diffractive decoder with electronic encoding to accurately transfer the optical information of interest, corresponding to, e.g., any arbitrary input object or message, through unknown random phase diffusers along the optical path. This hybrid electronic-optical model, trained using supervised learning, comprises a convolutional neural network-based electronic encoder and successive passive diffractive layers that are jointly optimized. After their joint training using deep learning, our hybrid model can transfer optical information through unknown phase diffusers, demonstrating generalization to new random diffusers never seen before. The resulting electronic-encoder and optical-decoder model was experimentally validated using a 3D-printed diffractive network that axially spans <70λ, where λ = 0.75 mm is the illumination wavelength in the terahertz spectrum, carrying the desired optical information through random unknown diffusers. The presented framework can be physically scaled to operate at different parts of the electromagnetic spectrum, without retraining its components, and would offer low-power and compact solutions for optical information transfer in free space through unknown random diffusive media.

36 MATERIALS SCIENCE↗

Computational Modeling of Water Infiltration From Ditches Into Roadway Embankments

Ditches are built along roadways to convey the stormwater runoff. They collect the runoff from the pavement and adjacent drainage areas. Water can move freely in a ditch away from the road surface, thus making it safer for the vehicular traffic. South Carolina DOT expressed interest in studying infiltration of water from ditches into roadway embankments. As part of the current design criteria, it is recommended that the road subgrades be 1.0 foot above the design event water surface elevation in the ditch. The distance between the water surface and bottom of the subbase is referred to as “freeboard”. The purpose of this requirement is to lower the risk of pavement damage by moving the flow away from the pavement, containing the captured flow, and conveying flow received from adjacent properties to outfalls. With changes in roadway design standards, larger shoulders and flatter cross-slopes have moved ditches further away from the pavement and the subbase. These changes could allow for a lower freeboard standard without reducing safety margins due to the increased lateral distance from the subbase to the ditch. The objective of this study is to evaluate whether the current SCDOT criterion requiring the roadway subgrade to be 1 ft above the design-event ditch water surface can be relaxed under selected geometric and soil conditions, particularly where larger shoulder widths and flatter side slopes increase horizontal separation between ditch water and pavement structure.

42 ENGINEERING↗

A microchannel shell-and-tube absorber for ammonia-water absorption

A shell-and-tube absorber with microscale geometric features that yields significant volumetric heat duties is investigated. A vertically-oriented baffled shell-and-tube heat exchanger with a tightly packed tube bundle is designed to serve as the absorber for a 10.5 kW nominal cooling capacity absorption chiller. The heat and mass transfer in this absorber is investigated experimentally at operating conditions spanning a wide range of working fluid and coupling fluid flow rates and ambient temperatures. An experimental facility is fabricated to closely simulate the operating conditions in an absorption refrigeration system. Absorption heat duties in the range of 5-15 kW are measured. In addition to operating the absorber in the baseline co-current flow configuration with premixed inlet, fluid ports are provided to enable operation in counter-current flow, where the dilute solution enters from the top and vapor enters from the bottom. The performance of these two flow configurations is compared. A preliminary heat and mass transfer model is developed to quantify the local temperatures, concentrations, and mass flow rates of the solution and vapor streams along the length of the absorber. The model predictions for total heat duty agree with the experimental results within 10% for almost all the data points, with an AAD of ~3.4%. Finally, the results of this study, along with the modeling framework, can guide further development of highly compact and efficient absorbers absorption systems.

42 ENGINEERING↗

An adaptive Hessian approximated stochastic gradient MCMC method

Bayesian approaches have been successfully integrated into training deep neural networks. One popular family is stochastic gradient Markov chain Monte Carlo methods (SG-MCMC), which have gained increasing interest due to their ability to handle large datasets and the potential to avoid overfitting. Although standard SG-MCMC methods have shown great performance in a variety of problems, they may be inefficient when the random variables in the target posterior densities have scale differences or are highly correlated. Here, we present an adaptive Hessian approximated stochastic gradient MCMC method to incorporate local geometric information while sampling from the posterior. The idea is to apply stochastic approximation (SA) to sequentially update a preconditioning matrix at each iteration. The preconditioner possesses second-order information and can guide the random walk of a sampler efficiently. Instead of computing and saving the full Hessian of the log posterior, we use limited memory of the samples and their stochastic gradients to approximate the inverse Hessian-vector multiplication in the updating formula. Moreover, by smoothly optimizing the preconditioning matrix via SA, our proposed algorithm can asymptotically converge to the target distribution with a controllable bias under mild conditions. To reduce the training and testing computational burden, we adopt a magnitude-based weight pruning method to enforce the sparsity of the network. Our method is user-friendly and demonstrates better learning results compared to standard SG-MCMC updating rules. The approximation of inverse Hessian alleviates storage and computational complexities for large dimensional models. Numerical experiments are performed on several problems, including sampling from 2D correlated distribution, synthetic regression problems, and learning the numerical solutions of heterogeneous elliptic PDE. The numerical results demonstrate great improvement in both the convergence rate and accuracy.

97 MATHEMATICS AND COMPUTING↗

Impact of particle arrays on phase separation composition patterns

We examine the symmetry-breaking effect of fixed constellations of particles on the surface-directed spinodal decomposition of binary blends in the presence of particles whose surfaces have a preferential affinity for one of the components. Our phase-field simulations indicate that the phase separation morphology in the presence of particle arrays can be tuned to have a continuous, droplet, lamellar, or hybrid morphology depending on the interparticle spacing, blend composition, and time. In particular, when the interparticle spacing is large compared to the spinodal wavelength, a transient target pattern composed of alternate rings of preferred and non-preferred phases emerges at early times, tending to adopt the symmetry of the particle configuration. We reveal that such target patterns stabilize for certain characteristic length, time, and composition scales characteristic of the pure phase-separating mixture. To illustrate the general range of phenomena exhibited by mixture-particle systems, we simulate the effects of single-particle, multi-particle, and cluster–particle systems having multiple geometrical configurations of the particle characteristic of pattern substrates on phase separation. Our simulations show that tailoring the particle configuration, or substrate pattern configuration, a relative fluid-particle composition should allow the desirable control of the phase separation morphology as in block copolymer materials, but where the scales accessible to this approach of organizing phase-separated fluids usually are significantly larger. Limited experiments confirm the trends observed in our simulations, which should provide some guidance in engineering patterned blend and other mixtures of technological interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-Richtmyer–Meshkov instability ejecta production based on shallow bubble collapse

The study of shock-driven ejecta production has focused on Richtmyer–Meshkov instability (RMI) growth from geometric features of the material surface. Extensive study of this mechanism under both single- and multiple-shock conditions has found that the ejected mass tends to be closely associated with the shocked surface phase, and its temperature is not dramatically greater than the hydrodynamic shock temperature of the bulk. In this work, we propose and demonstrate a new ejecta production mechanism that can occur under multiple-shock conditions based on the collapse of bubbles near the free surface of the material. This mechanism produces ejected mass that is much greater in quantity than observed in the RMI case. The particles are much hotter than predicted by the shock Hugoniot state, and the ejected mass does not appear to be strongly dependent upon initial surface finish. The ejecta source extends into the material with no clear remaining free surface. We name this mechanism Shallow Bubble Collapse (SBC) and discuss the conditions under which it activates. We demonstrate resolved modeling methods that enable the calculation, design, and study of SBC as a mechanism and perform a series of experiments to compare with the models. Under some multiple-shock conditions, SBC ejection produces ten times more ejected mass than RMI growth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Revealing melt flow instabilities in laser powder bed fusion additive manufacturing of aluminum alloy via in-situ high-speed X-ray imaging

Laser metal additive manufacturing technologies enable the fabrication of geometrically and compositionally complex parts unachievable by conventional manufacturing methods. However, the certification and qualification of additively manufactured parts are greatly hindered by the stochastic melt flow instabilities intrinsic to the process, which has not been explicitly revealed by direct observation. Here, we report the mechanisms of the melt flow instabilities in laser powder bed fusion additive manufacturing process revealed by in-situ high-speed high-resolution synchrotron X-ray imaging. Here we identified powder/droplet impact, significant keyhole oscillation, and melting-mode switching as three major mechanisms for causing melt flow instabilities. We demonstrated the detrimental consequences of these instabilities brought to the process, and presented new understanding on the melt flow evolution and keyhole oscillation. This work provides critical insights into process instabilities during laser metal additive manufacturing, which may guide the development of instability mitigation approaches. The results reported here are also important for the development and validation of high-fidelity computational models.

36 MATERIALS SCIENCE↗

Crossing Symmetric Dispersion Relations in Quantum Field Theories

For 2-2 scattering in quantum field theories, the usual fixed t dispersion relation exhibits only two-channel symmetry. This Letter considers a crossing symmetric dispersion relation, reviving certain old ideas from the 1970s. Rather than the fixed t dispersion relation, this needs a dispersion relation in a different variable z, which is related to the Mandelstam invariants s, t, u via a parametric cubic relation making the crossing symmetry in the complex z plane a geometric rotation. The resulting dispersion is manifestly three-channel crossing symmetric. We give simple derivations of certain known positivity conditions for effective field theories, including the null constraints, which lead to two sided bounds and derive a general set of new nonperturbative inequalities. We show how these inequalities enable us to locate the first massive string state from a low energy expansion of the four dilaton amplitude in type II string theory. We also show how a generalized (numerical) Froissart bound, valid for all energies, is obtained from this approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Inferring Stochastic Rates from Heterogeneous Snapshots of Particle Positions

Many imaging techniques for biological systems—like fixation of cells coupled with fluorescence microscopy—provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. In this study, these snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

59 BASIC BIOLOGICAL SCIENCES↗

Machine Learning for Prediction of Thermodynamic Descriptors

Our objective is to apply machine learning (ML) algorithms for the prediction of molecular catalysis descriptors from geometric properties derived from experimental crystallographic databases. Catalysis is often considered a “low-data” discipline that is poorly suited for ML methods. An exception is the extensive structural information that is available for molecular catalysts through the Cambridge Structural Database (CSD), which contains atomically precise molecular structures from X-ray diffraction analysis for >600K metal complexes. As a proof-of-principle, we targeted the prediction of hydricity, a thermodynamic property that provides understanding and control of catalytic hydride transfer. We built a training set composed of ~100 molecular complexes with a known hydricity and structural information from the CSD. This data set was converted into a machine-readable format using the smooth overlap of atomic positions (SOAP) representation and further labeled with simple electronic descriptors for the metal centers. Multiple different neural networks were trained on this data set, and the accuracy of the hydricity predictions ranged from < 2 kcal/mol to 20 kcal/mol. The accuracy of each model was highly sensitive to which compounds were in the train versus test set, underscoring the challenges associated with small and chemically diverse data sets. Finally, to further augment the data set, we attempted to experimentally measure several new hydricity values, however these experiments were unsuccessful due to undesired chemical reactivity of the selected complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Generalising G 2 geometry: involutivity, moment maps and moduli

We analyse the geometry of generic Minkowski N = 1, D = 4 flux compactifications in string theory, the default backgrounds for string model building. In M-theory they are the natural string theoretic extensions of G2 holonomy manifolds. In type II theories, they extend the notion of Calabi-Yau geometry and include the class of flux backgrounds based on generalised complex structures first considered by Graña et al. (GMPT). Using E 7(7) × $\mathbb{R}$ + generalised geometry we show that these compactifications are characterised by an SU(7) ⊂ E 7(7) structure defining an involutive subbundle of the generalised tangent space, and with a vanishing moment map, corresponding to the action of the diffeomorphism and gauge symmetries of the theory. The Kähler potential on the space of structures defines a natural extension of Hitchin’s G 2 functional. Using this framework we are able to count, for the first time, the massless scalar moduli of GMPT solutions in terms of generalised geometry cohomology groups. It also provides an intriguing new perspective on the existence of G 2 manifolds, suggesting possible connections to Geometrical Invariant Theory and stability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Convergence study of wakefield simulations with GdfidL and ECHO3D

The interaction of charged particle beams with vacuum chamber components gives rise to electromagnetic wakefields, whose frequency-domain representation is known as beam coupling impedance. Geometric impedance arising from discontinuities and transitions in the vacuum chamber is the focus of this study. Minimizing this impedance is essential to mitigate adverse collective effects in modern storage rings operating with high-intensity particle beams. Accurate and reliable impedance simulations is a key factor of the vacuum chamber design. This paper presents the results of a convergence study of two widely used electromagnetic solvers, GdfidL and ECHO3D, applied to key vacuum-chamber components of the National Synchrotron Light Source II (NSLS-II) storage ring. Detailed comparisons are performed for several geometries, including flange absorbers, RF bellows, button-type beam position monitors, and an in-vacuum undulator (IVU). The results show notable differences in convergence and computational efficiency between the two codes. While GdfidL provides highly resolved results and serves as a common benchmark tool, ECHO3D yields consistent results with coarser meshes, significantly reducing simulation time and memory demands. Simulations with a full-geometry IVU model demonstrate that simplified taper-transition models can miss important impedance contributions. In conclusion, these findings provide practical guidelines for efficient and accurate impedance modeling to optimize design of vacuum chamber components for accelerators.

36 MATERIALS SCIENCE↗

Properties of interfaces between copper and copper sulphide/oxide films

We compare different means of predicting and rationalizing properties of corrosion films aiming at gaining deeper understanding of the behaviour of copper in aqueous and sulphide-containing environments. Purely geometrical tools ranging from mean bulk information to anisotropic and facet-specific approaches are discussed, and their performances compared against DFT-based stability and property evaluations of a range of interfaces arising from combining low-index crystallographic planes of the compounds Cu, Cu 2 O, and Cu 2 S. We demonstrate the necessity to consider facet-specific interactions to understand nanolevel differences between Cu 2 S and Cu 2 O film behaviour, and that thin films cannot be directly described by bulk properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Embedded Sensing in Additive Manufacturing Metal and Polymer Parts: A Comparative Study of Integration Techniques and Structural Health Monitoring Performance

This study presents a comparative evaluation of post-process sensor integration in additively manufactured (AM) metal and the in-situ process for polymer structures for structural health monitoring (SHM), with an emphasis on embedded sensors. Geometrically identical specimens were fabricated using copper via metal fused filament fabrication (FFF) and PLA via polymer FFF, with piezoelectric transducers (PZTs) inserted into internal cavities to assess the influence of material and placement on sensing fidelity. Mechanical testing under compressive and point loads generated signals that were transformed into time–frequency spectrograms using a Short-Time Fourier Transform (STFT) framework. An engineered RGB representation was developed, combining global amplitude scaling with an amplitude-envelope encoding to enhance contrast and highlight subtle wave features. These spectrograms served as inputs to convolutional neural networks (CNNs) for classification of load conditions and detection of damage-related features. Results showed reliable recognition in both copper and PLA specimens, with CNN classification accuracies exceeding 95%. Embedded PZTs were especially effective in PLA, where signal damping and environmental sensitivity often hinder surface-mounted sensors. This work demonstrates the advantages of embedded sensing in AM structures, particularly when paired with spectrogram-based feature engineering and CNN modeling, advancing real-time SHM for aerospace, energy, and defense applications.

additive manufacturing↗

Virtual Flow Solver - Geophysics: A 3D Incompressible Navier-Stokes Solver

Virtual Flow Solver - Geophysics (VFS-Geophysics) is a three-dimensional (3D) incompressible Navier-Stokes solver based on the Curvilinear Immersed Boundary (CURVIB) method. The CURVIB is a sharp interface type of immersed boundary (IB) method that enables the simulation of fluid flows in the presence of geometrically complex moving bodies. The CURVIB method can be applied to wind/MHK turbine simulations and energy applications. VFS-Geophysics is the result of many years of research work by several graduate students, post-docs, and research associates that have been involved in the Computational Hydrodynamics and Biofluids Laboratory directed by Professor Fotis Sotiropoulos. The preparation of the present manual has been supported by the U.S. Department of Energy (DE-EE 0005482).

16 TIDAL AND WAVE POWER↗