Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Broad-band Spectroscopy of Hercules X-1 with Suzaku

Hercules X-1 was observed with Suzaku in the main-on state from 2005 to 2010. The 0.4–100 keV wide-band spectra obtained in four observations showed a broad hump around 4–9 keV in addition to narrow Fe lines at 6.4 and 6.7 keV. The hump was seen in all the four observations regardless of the selection of the continuum models. Thus it is considered a stable and intrinsic spectral feature in Her X-1. The broad hump lacked a sharp structure like an absorption edge. Thus it was represented by two different spectral models: an ionized partial covering or an additional broad line at 6.5 keV. The former required a persistently existing ionized absorber, whose origin was unclear. In the latter case, the Gaussian fitting of the 6.5-keV line needs a large width of σ = 1.0–1.5 keV and a large equivalent width of 400–900 eV. If the broad line originates from Fe fluorescence of accreting matter, its large width may be explained by the Doppler broadening in the accretion flow. However, the large equivalent width may be inconsistent with a simple accretion geometry.

Asami, Fumi↗
NASA NTRS2023

Radiation Interaction with Particulates in Planetary Entry Flows

A hypersonic particle phase solver based on a continuum model is implemented in the NASA HyperSolve CFD code and used to predict particle phase conditions in a dusty Mars entry scenario. Black-body radiative emission from high temperature dust particles is predicted for a Mars relevant entry problem, and the emission integrated to an entry vehicle using a ray tracing approach to estimate surface heat flux. The surface heating due to this mechanism is found to be small relative to convective and radiation from the gas phase emission. Additionally, the effect of atmospheric haze particulates vaporizing in the Titan atmosphere is studied. The increase in carbon availability due to the composition of the haze particles is found to increase the gas radiative heating to the vehicle by up to 3% on the forebody of a representative Titan entry vehicle with approximately 0.15% haze in the freestream by mass.

Andrew Hinkle↗
NASA NTRS2025

Assessment of Radiation Scattering by Dust in Mars Entry

The impact of gas phase radiation scattering by dust particle presence in Mars entry scenarios is assessed using a suite of numerical tools. A hypersonic gas and dilute particle multiphase flow is simulated using a mixed continuum model and an efficient size distribution quadrature scheme. Radiative transport is simulated with a hybrid ray-tracing and P1 spherical harmonics model. The NASA HARA radiation code and a Lorentz-Mie scattering model are used to predict the gas-particle mixture radiative properties. Application of the numerical suite to a representative case based on the Mars 2020 entry shows that, even at conditions in excess of the strongest recorded dust storm, the augmentation in surface radiative heat flux due to scattering alone is lower than 1%, and therefore negligible. Predictions with even higher dust loading conditions confirm that even highly uncertain dust mass profiles are unlikely to be a concern for the scattering augmentation mechanism.

Andrew Hinkle↗
DOE OSTI2021

Improving scalability of electronic structure code for molecular simulations in the presence of environment

A scalable density functional electronic code with Gaussian basis set, called UTEP-NRLMOL, is developed to perform simulations of molecular systems in the presence of the environment with particular attention to the memory requirements. In the electronic structure calculations, the memory and computation time are proportional to the number of atoms. Memory requirements for density functional calculations scale as N*N, where N is the number of atoms. While the recent advances in HPC offer platforms with large numbers of cores, the limited amount of memory available on a given node and poor scalability of the electronic structure codes hinder their efficient usage of these platforms. We have introduced new scaling and parallelization paradigms using MPI-3 shared-memory functionality combined with usage of sparse algebra and storage of matrices in sparse format. This extends the range of applicability of the UTEP-NRLMOL code to large systems over 10,000 atoms, or using up to 67,000 basis functions, and making use of HPC architectures using over 6,000 processors utilizing all available cores. We have also interfaced code with effective fragment potential and polarizable continuum model libraries. The code was used in simulations of several applications which are published in reputed scientific journals.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2011

Continuum Damage Mechanics Models for the Analysis of Progressive Failure in Open-Hole Tension Laminates

The performance of a state-of-the-art continuum damage mechanics model for interlaminar damage, coupled with a cohesive zone model for delamination is examined for failure prediction of quasi-isotropic open-hole tension laminates. Limitations of continuum representations of intra-ply damage and the effect of mesh orientation on the analysis predictions are discussed. It is shown that accurate prediction of matrix crack paths and stress redistribution after cracking requires a mesh aligned with the fiber orientation. Based on these results, an aligned mesh is proposed for analysis of the open-hole tension specimens consisting of different meshes within the individual plies, such that the element edges are aligned with the ply fiber direction. The modeling approach is assessed by comparison of analysis predictions to experimental data for specimen configurations in which failure is dominated by complex interactions between matrix cracks and delaminations. It is shown that the different failure mechanisms observed in the tests are well predicted. In addition, the modeling approach is demonstrated to predict proper trends in the effect of scaling on strength and failure mechanisms of quasi-isotropic open-hole tension laminates.

Song, Kyonchan↗
NASA NTRS2002

Applicability of the Continuum-Shell Theories to the Mechanics of Carbon Nanotubes

Validity of the assumptions relating the applicability of continuum shell theories to the global mechanical behavior of carbon nanotubes is examined. The present study focuses on providing a basis that can be used to qualitatively assess the appropriateness of continuum-shell models for nanotubes. To address the effect of nanotube structure on their deformation, all nanotube geometries are divided into four major classes that require distinct models. Criteria for the applicability of continuum models are presented. The key parameters that control the buckling strains and deformation modes of these classes of nanotubes are determined. In an analogy with continuum mechanics, mechanical laws of geometric similitude are presented. A parametric map is constructed for a variety of nanotube geometries as a guide for the applicability of different models. The continuum assumptions made in representing a nanotube as a homogeneous thin shell are analyzed to identify possible limitations of applying shell theories and using their bifurcation-buckling equations at the nano-scale.

Harik, V. M.↗
NASA NTRS2017

Multiscale Fiber Kinking: Computational Micromechanics and a Mesoscale Continuum Damage Mechanics Models

In this work, the fiber kinking phenomenon, which is known as the failure mechanism that takes place when a fiber reinforced polymer is loaded under longitudinal compression, is studied. A computational micromechanics model is employed to interrogate the assumptions of a recently developed mesoscale continuum damage mechanics (CDM) model for fiber kinking based on the deformation gradient decomposition (DGD) and the LaRC04 failure criteria.

Herraez, Miguel↗
DOE OSTI2022

Tree hydrodynamic modelling of the soil–plant–atmosphere continuum using FETCH3

Modelling the water transport along the soil–plant–atmosphere continuum is fundamental to estimating and predicting transpiration fluxes. A Finite-difference Ecosystem-scale Tree Crown Hydrodynamics model (FETCH3) for the water fluxes across the soil–plant–atmosphere continuum is presented here. The model combines the water transport pathways into one vertical dimension, and assumes that the water flow through the soil, roots, and above-ground xylem can be approximated as flow in porous media. This results in a system of three partial differential equations, resembling the Richardson–Richards equation, describing the transport of water through the plant system and with additional terms representing sinks and sources for the transfer of water from the soil to the roots and from the leaves to the atmosphere. The numerical scheme, developed in Python 3, was tested against exact analytical solutions for steady state and transient conditions using simplified but realistic model parameterizations. The model was also used to simulate a previously published case study, where observed transpiration rates were available, to evaluate model performance. With the same model setup as the published case study, FETCH3 results were in agreement with observations. Through a rigorous coupling of soil, root xylem, and stem xylem, FETCH3 can account for variable water capacitance, while conserving mass and the continuity of the water potential between these three layers. FETCH3 provides a ready-to-use open access numerical model for the simulation of water fluxes across the soil–plant–atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Analysis of Open Hole Compression Specimens Using the CompDam Continuum Damage Mechanics Model

This report documents a validation study on the standard Open Hole Compression (OHC) laminate test specimen conducted as part of the NASA Advanced Composites Project (ACP). Tests were conducted on OHC specimens with hard, quasi-isotropic, and soft layups using digital image correlation and X-Ray computed tomography to capture the structural response and damage evolution. Progressive damage models were constructed for use with the CompDam continuum damage mechanics code following the best practices established during the ACP. Detailed interrogation of the analysis results and comparison with experimental measurements provide a basis for assessing the capability of the modeling approach for OHC. The structural response is found to be captured well, with strengths predicted within 3\% of the experimental values for hard and quasi-isotropic laminates. In the soft laminate, the model predicts failure to be more brittle than the nonlinear, ductile response that was measured. Damage states extracted from the models at the same load level as test measurements are overlaid to show directly the similarities and differences between test and analysis results. Studying the damage evolution predicted by the analysis reveals that the failure process is a competition between fiber damage and delamination/sub-laminate buckling, with fiber damage dominating the collapse in the hard laminate and sub-laminate buckling governing in the soft laminate. Finally, a series of parametric studies varying numerical solution parameters (mesh size, mass scaling) and physical properties (fiber direction compressive strength and toughness) reveal sensitivities and deficiencies of the model. To the authors' knowledge, this study is the first for OHC specimens to include detailed evaluation of damage mode interactions, direct overlay of predicted and measured damage states, and sensitivity of the predicted results to difficult-to-measure fiber direction material properties.

Continuum damage mechanics↗
NASA NTRS2016

A Continuum Damage Mechanics Model to Predict Kink-Band Propagation Using Deformation Gradient Tensor Decomposition

A new model is proposed that represents the kinematics of kink-band formation and propagation within the framework of a mesoscale continuum damage mechanics (CDM) model. The model uses the recently proposed deformation gradient decomposition approach to represent a kink band as a displacement jump via a cohesive interface that is embedded in an elastic bulk material. The model is capable of representing the combination of matrix failure in the frame of a misaligned fiber and instability due to shear nonlinearity. In contrast to conventional linear or bilinear strain softening laws used in most mesoscale CDM models for longitudinal compression, the constitutive response of the proposed model includes features predicted by detailed micromechanical models. These features include: 1) the rotational kinematics of the kink band, 2) an instability when the peak load is reached, and 3) a nonzero plateau stress under large strains.

Bergan, Andrew C.↗
NASA NTRS2002

Molecular Line and Continuum Opacities for Modeling of Extrasolar Giant Planet and Cool Stellar Atmospheres

The molecular line and continuum opacities are investigated in the atmospheres of cool stars and Extrasolar Giant Planets (EGPs). Using a combination of ab inito and experimentally derived potential curves and dipole transition moments, accurate data have been calculated for rovibrationally-resolved oscillator strengths and photodissociation cross sections in the B' (sup 2)Sigma+ (left arrow) X (sup 2)Sigma+ and A (sup 2)Pi (left arrow) X (sup 2)Sigma+ band systems in MgH. We also report our progress on the study of the electronic structure of LiCl and FeH.

Weck, P. F.↗
DOE OSTI2026

Continuum Correlations from CFD-DEM Modeling of Conduction Heat Transfer in Granular Flows

Heat transfer between a surface and flowing particles is analyzed to improve the accuracy of continuum models for wall-to-bed heat transfer in a fluidized bed. Discrete element modeling (DEM) is used to model a fluidized bed heat exchanger where heat enters the system through a heated wall. The DEM heat transfer predictions are validated against published experimental work (Brewster et al., 2024) with less than 15% error. In previous work by Morris et al. (2015), a continuum model was developed using data from high-fidelity DEM simulations of chute flows. In the current study, the continuum model is extended and validated for fluidized beds. The sensitivity of the continuum heat transfer model parameters, which was not quantified in previous studies, is also investigated. It is observed that for a given particle with specific properties, e.g. the particle size, roughness, and conduction lens radius, the continuum correlation developed for heat transfer from a heated boundary to the particle bed depends mainly on the solid fraction or porosity of the particle bed for a given fluid. The new continuum heat transfer model is then validated over a wide range of superficial velocities via comparisons to both discrete element and experimental data. It is shown that this correlation is valid for a large range of particle flow conditions from chute flows to fluidized beds with less than 10% error as compared to DEM predictions.

14 SOLAR ENERGY↗
NASA NTRS1996

Simulation of Two-Fluid Flows by the Least-Squares Finite Element Method Using a Continuum Surface Tension Model

In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols and the continuum surface force (CSF) model developed by Brackbill, et al. In the VOF method fluids of different properties are identified through the use of a continuous field variable (color function). The color function assigns a unique constant (color) to each fluid. The interfaces between different fluids are distinct due to sharp gradients of the color function. The evolution of the interfaces is captured by solving the convective equation of the color function. The CSF model is used as a means to treat surface tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin, is used to calculate the tension force. In the modified version, the force term is obtained by calculating the divergence of a stress tensor defined by the gradient of the color function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the interfaces. The least-squares finite element method (LSFEM) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a unified manner without any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both two dimensional planar and axisymmetric flows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer.

Wu, Jie↗
DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗
DOE OSTI2022

Kinetic modeling of neutral transport for a continuum gyrokinetic code

In this work we present the first-of-its-kind coupling of a continuum full- f gyrokinetic turbulence model with a 6D continuum model for kinetic neutrals, carried out using the Gkeyll code. Our objective is to improve the first-principle understanding of the role of neutrals in plasma fueling, detachment, and their interaction with edge plasma profiles and turbulence statistics. Our model includes only atomic hydrogen and incorporates electron-impact ionization, charge exchange, and wall recycling. These features have been successfully verified with analytical predictions and benchmarked with the DEGAS2 Monte Carlo neutral code. We carry out simulations for a scrape-off layer (SOL) with simplified geometry and National Spherical Torus Experiment parameters. We compare these results to a baseline simulation without neutrals and find that neutral interactions reduce the normalized density fluctuation levels and associated skewness and kurtosis, while increasing auto-correlation times. A flatter density profile is also observed, similar to the SOL density shoulder formation in experimental scenarios with high fueling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Correcting implicit solvation at metal/water interfaces through the incorporation of competitive water adsorption

Conventional continuum solvation models are ubiquitous in computational catalysis, including for describing metal/water interfaces, which are relevant to both solution-phase heterogeneous catalysis and electrocatalysis. Nonetheless, we find that such continuum models qualitatively fail to describe both the adsorption free energy and conformational preference for many organic molecules at such interfaces, largely due to the failure of continuum models to incorporate the role of competitive water adsorption. We develop a simple phenomenological model that accounts for competitive water adsorption and show that the model, when used in conjunction with continuum solvation, provides a dramatic improvement in the description of both adsorption and conformational preference. The model is also extended to additionally incorporate the influence of applied potential at the electrode surface, thus facilitating computationally efficient applications to scenarios including electrocatalysis.

Chemistry↗
DOE OSTI2021

Open momentum space method for the Hofstadter butterfly and the quantized Lorentz susceptibility

Here we develop a generic k ∙ p open momentum space method for calculating the Hofstadter butterfly of both continuum (moiré) models and tight-binding models, where the quasimomentum is directly substituted by the Landau level (LL) operators. By taking a LL cutoff (and a reciprocal lattice cutoff for continuum models), one obtains the Hofstadter butterfly with in-gap spectral flows. For continuum models such as the moiré model for twisted bilayer graphene, our method gives a sparse Hamiltonian, making it much more efficient than existing methods. The spectral flows in the Hofstadter gaps can be understood as edge states on a momentum space boundary, from which one can determine the two integers (t ν , s ν ) of a gap ν satisfying the Diophantine equation. The spectral flows can also be removed to obtain a clear Hofstadter butterfly. While t ν is known as the Chern number, our theory identifies s ν as a dual Chern number for the momentum space, which corresponds to a quantized Lorentz susceptibility γ xy = eBs ν .

2-dimensional systems↗
DOE OSTI2020

Neutron-rich calcium isotopes within realistic Gamow shell model calculations with continuum coupling

Based on the realistic nuclear force of the high-precision CD-Bonn potential, we have performed comprehensive calculations for neutron-rich calcium isotopes using the Gamow shell model (GSM) which includes resonance and continuum. The realistic GSM calculations produce well binding energies, one- and two-neutron separation energies, predicting that 57 Ca is the heaviest bound odd isotope and 70 Ca is the dripline nucleus. Resonant states are predicted, which provides useful information for future experiments on particle emissions in neutron-rich calcium isotopes. Shell evolutions in the calcium chain around neutron numbers $\textit{N}$ = 32, 34, and 40 are understood by calculating effective single-particle energies, the excitation energies of the first 2 + states and two-neutron separation energies. The calculations support shell closures at 52 Ca ($\textit{N}$ = 32) and 54 Ca ($\textit{N}$ = 34) but show a weakening of shell closure at 60 Ca ($\textit{N}$ = 40). Finally, the possible shell closure at 70 Ca ($\textit{N}$ = 50) is predicted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗