Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing continuum”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

White light and radio studies of the coronal transient of 14-15 September 1973. I - Material motions and magnetic field

Observations of a coronal transient event were obtained in white light by the Skylab coronagraph and at metric wavelengths by the radioheliograph and spectrograph at Culgoora, Australia, and the spectrograph-interferometer at Boulder, Colo. The continuum radio burst was found to originate above the outward-moving white-light loop, a region of compressed material headed by a bow wave. The computed density in the region of radio emission, based on either gyrosynchrotron or harmonic plasma radiation mechanisms, was approximately 10 times the ambient coronal density; this is compatible with the density deduced from the white-light observations. The magnetic-energy density derived from the radio observations was greater than 10 times the thermal energy density, marginally larger than the kinetic energy density in the fastest-moving portion of the transient, and considerably larger in most other regions. The ambient medium, the white-light front, the compression region, the loop, and the slower massive flow of material behind are each examined. It is found that the plasma was magnetically controlled throughout and that magnetic forces provided the principal mechanism for acceleration of the transient material from the sun.

Dulk, G. A.↗
NASA NTRS1986

The mathematical modeling of rapid solidification processing

The detailed formulation of and the results obtained from a continuum mechanics-based mathematical model of the planar flow melt spinning (PFMS) rapid solidification system are presented and discussed. The numerical algorithm proposed is capable of computing the cooling and freezing rates as well as the fluid flow and capillary phenomena which take place inside the molten puddle formed in the PFMS process. The FORTRAN listings of some of the most useful computer programs and a collection of appendices describing the basic equations used for the modeling are included.

Gutierrez-Miravete, E.↗
NASA NTRS1986

Design sensitivity analysis using EAL. Part 2: Shape design parameters

A numerical implementation of shape design sensitivity analysis of built-up structures is presented, using the versatility and convenience of an existing finite element structural analysis code and its data base management system. This report is a continuation of a previous report on conventional design parameters. The finite element code used in the implementation presented is the Engineering Analysis Language (EAL), which is based on a hybrid analysis method. It has been shown that shape design sensitivity computations can be carried out using the database management system of EAL, without writing a separate program and a separate data base. The material derivative concept of continuum mechanics and an adjoint variable method of design sensitivity analysis are used to derive shape design sensitivity information of structural performances. A domain method of shape design sensitivity analysis and a design component method are used. Displacement and stress functionals are considered as performance criteria.

Dopker, B.↗
NASA NTRS1997

Aerodynamics of the Mars Microprobe Entry Vehicles

The selection of the unique aeroshell shape for the Mars Microprobes is discussed. A description of its aerodynamics in hypersonic rarefied, hypersonic continuum, supersonic and transonic flow regimes is then presented. This description is based on Direct Simulation Monte Carlo analyses in the rarefied-flow regime, thermochemical nonequilibrium Computational Fluid Dynamics in the hypersonic regime, existing wind tunnel data in the supersonic and transonic regime, additional computational work in the transonic regime, and finally, ballistic range data. The aeroshell is shown to possess the correct combination of aerodynamic stability and drag to convert the probe's initial tumbling attitude and high velocity at atmospheric-interface into the desired surface-impact orientation and velocity.

Mitcheltree, R. A.↗
NASA NTRS2018

Parameterizations of Chromospheric Condensations in dG and dMe Model Flare Atmospheres

The origin of the near-ultraviolet and optical continuum radiation in flares is critical for understanding particle acceleration and impulsive heating in stellar atmospheres. Radiative-hydrodynamic (RHD) simulations in 1D have shown that high energy deposition rates from electron beams produce two flaring layers at T approximately 10 (exp 4) K that develop in the chromosphere: a cooling condensation (downflowing compression) and heated non-moving (stationary) flare layers just below the condensation. These atmospheres reproduce several observed phenomena in flare spectra, such as the red-wing asymmetry of the emission lines in solar flares and a small Balmer jump ratio in M dwarf flares. The high beam flux simulations are computationally expensive in 1D, and the (human) timescales for completing NLTE models with adaptive grids in 3D will likely be unwieldy for some time to come. We have developed a prescription for predicting the approximate evolved states, continuum optical depth, and emergent continuum flux spectra of RHD model flare atmospheres. These approximate prescriptions are based on an important atmospheric parameter: the column mass (m(sub ref)) at which hydrogen becomes nearly completely ionized at the depths that are approximately in steady state with the electron beam heating. Using this new modeling approach, we find that high energy flux density (>F11) electron beams are needed to reproduce the brightest observed continuum intensity in IRIS data of the 2014 March 29 X1 solar flare, and that variation in m(sub ref) from 0.001 to 0.02 g cm (exp -2) reproduces most of the observed range of the optical continuum flux ratios at the peak of M dwarf flares.

Kowalski, Adam F.↗
DOE OSTI2022

Calibrating uncertain parameters in melt pool simulations of additive manufacturing

Melt pool scale numerical modeling of additive manufacturing (AM) processes can provide predictive capabilities and theoretical insight into the process-property-structure-performance relationships for AM parts. Despite capabilities of numerical models to solve complex multi-physics problems, it is often important to consider a tradeoff between detailed physics and computational cost. Therefore, sources of uncertainty in both experimental conditions and the parameters needed for modeling require models to be validated against empirical evidence. Here, a method is proposed to calibrate uncertain parameters used in continuum-scale melt pool models for powder bed fusion (PBF) AM. Both a simplified heat transfer model and a heat transfer and fluid flow model were investigated. A surrogate model and Markov chain-based optimization algorithm calibrated melt pool geometry for models within experimental variation of the target melt pool width and depth from the NIST AM-Bench 2018-02 dataset. The melt pool temperature distributions, solidification parameters, and simulated multi-layer solidification microstructures were compared between the two models. Similar results from both models indicate that calibrated, lower fidelity numerical models may be used in place of higher fidelity models to generate melt pool solidification data. Finally, these calibrated models therefore enable lower computational cost melt pool simulations without a noticeable decrease in simulation accuracy for grain-scale microstructure simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiscale Concurrent Atomistic-Continuum (CAC) modeling of multicomponent alloys

We report strengthening in complex multicomponent systems such as solid solution alloys is controlled primarily by the dynamic interactions between dislocation lines and heterogeneously distributed solute species. Modeling of extended defect length scales in such multicomponent systems becomes prohibitively expensive, motivating the development of reduced order approaches. This work explores the application of the Concurrent Atomistic-Continuum (CAC) method to model dislocation mobility in random alloys at extended length scales. By employing recently developed average-atom interatomic potentials, the average “bulk” material response in coarse-grained regions interacts with true random solute species in the atomistic-scale domain. We demonstrate that spurious stresses in domain resolution transition regions are eliminated entirely due to the CAC formulation. Simultaneously, the key details of local stress fluctuation due to randomness in the dislocation core region are captured, and fluctuating stress smoothly decays to the long-range dislocation stress field response. Dislocation mobility calculations, for line lengths over 400 nm, are computed as a function of alloy composition in the model FeNiCr system and compared to full molecular dynamics (MD). The results capture the composition-dependent trends, while reducing degrees of freedom by nearly 40%. This approach can be readily extended to any system described by an EAM potential and facilitates the study of large-scale defect dynamics in complex solute environments to support computational alloy design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗
DOE OSTI2021

Stress-Dependent Chemo-Mechanical Performance of Amorphous Si Anodes for Li-Ion Batteries upon Lithiation

Alloying-type anodes are significantly governed by their chemo-mechanical performance during the electrochemical cycling. The reaction-induced huge volumetric change of these anodes may cause material degradation and failure under mechanical constraints. Here, we investigate the stress-dependent lithiation behavior of amorphous Si (a-Si) anodes using molecular dynamics simulations. It is indicated that a-Si anodes can sustain higher hydrostatic stress than biaxial/uniaxial ones without the occurrence of mechanical failure. Thermodynamic and electrochemical calculations demonstrate that although the lithiation procedure also affects the thermodynamic stability of a-Si anodes, it is mainly dominated by the external mean stresses. Compressive stress is confirmed to destabilize a-Si anodes and further trigger their capacity fading. Compared with our atomistic simulations, previous continuum models underestimate the open-cell potentials of a-Si anodes, due to their ignored large volumetric deformation at higher stresses and Li concentrations. Finally, this computational study provides the intensive atomic-level understanding of the stress-dependent lithiation behavior of a-Si anodes.

25 ENERGY STORAGE↗
NASA NTRS1976

Fundamental analysis of the failure of polymer-based fiber reinforced composites

A mathematical model is described which will permit predictions of the strength of fiber reinforced composites containing known flaws to be made from the basic properties of their constituents. The approach was to embed a local heterogeneous region (LHR) surrounding the crack tip into an anisotropic elastic continuum. The model should (1) permit an explicit analysis of the micromechanical processes involved in the fracture process, and (2) remain simple enough to be useful in practical computations. Computations for arbitrary flaw size and orientation under arbitrary applied load combinations were performed from unidirectional composites with linear elastic-brittle constituent behavior. The mechanical properties were nominally those of graphite epoxy. With the rupture properties arbitrarily varied to test the capability of the model to reflect real fracture modes in fiber composites, it was shown that fiber breakage, matrix crazing, crack bridging, matrix-fiber debonding, and axial splitting can all occur during a period of (gradually) increasing load prior to catastrophic fracture. The computations reveal qualitatively the sequential nature of the stable crack process that precedes fracture.

Kanninen, M. F.↗
NASA NTRS1976

Analysis of simultaneous Skylab and ground based flare observations

HeI and HeII resonance line data from Skylab were reduced, analyzed and compared with HeI D3 line intensities taken simultaneously from the Lockheed Rye Canyon Solar Observatory. Computer codes were developed for the calculation of total He line intensities and line profiles from model flare regions. These codes incorporate simultaneous solution of the line and continuum transport equations as needed together with the statistical equilibrium equations for a 30 level HeI, HeII, HeIII system. The energy level model consists of all terms through principal quantum number four. Interpretation of the observed data in terms of these parametric solutions and with simultaneous solution of the transport equations are discussed.

Kulander, J. L.↗
NASA NTRS1991

M.S.L.A.P. Modular Spectral Line Analysis Program documentation

MSLAP is a software for analyzing spectra, providing the basic structure to identify spectral features, to make quantitative measurements of this features, and to store the measurements for convenient access. MSLAP can be used to measure not only the zeroth moment (equivalent width) of a profile, but also the first and second moments. Optical depths and the corresponding column densities across the profile can be measured as well for sufficiently high resolution data. The software was developed for an interactive, graphical analysis where the computer carries most of the computational and data organizational burden and the investigator is responsible only for all judgement decisions. It employs sophisticated statistical techniques for determining the best polynomial fit to the continuum and for calculating the uncertainties.

Joseph, Charles L.↗
NASA NTRS2010

Analysis of Plume Impingement Effects from Orion Crew Service Module Dual Reaction Control System Engine Firings

Plume impingement effects on the Orion Crew Service Module (CSM) were analyzed for various dual Reaction Control System (RCS) engine firings and various configurations of the solar arrays. The study was performed using a decoupled computational fluid dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) approach. This approach included a single jet plume solution for the R1E RCS engine computed with the General Aerodynamic Simulation Program (GASP) CFD code. The CFD solution was used to create an inflow surface for the DSMC solution based on the Bird continuum breakdown parameter. The DSMC solution was then used to model the dual RCS plume impingement effects on the entire CSM geometry with deployed solar arrays. However, because the continuum breakdown parameter of 0.5 could not be achieved due to geometrical constraints and because high resolution in the plume shock interaction region is desired, a focused DSMC simulation modeling only the plumes and the shock interaction region was performed. This high resolution intermediate solution was then used as the inflow to the larger DSMC solution to obtain plume impingement heating, forces, and moments on the CSM and the solar arrays for a total of 21 cases that were analyzed. The results of these simulations were used to populate the Orion CSM Aerothermal Database.

Prisbell, Andrew↗
NASA NTRSDate not supplied

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗
NASA NTRSDate not supplied

Finite Element Modeling of Plastic Deformation During Spin Forming of Aluminum 6061-O

Spin and flow forming are metal deformation techniques in which a disk or tube of material is radially thinned and axially lengthened over a rotating mandrel. A finite-element continuum model (FEM) based on commercial software (DEFORM®) was developed to determine the stress/strain distribution and damage accumulation in a flow-formed, near-net-shape part. Computational validation and experimental verification of the model was leveraged to assess the effectiveness of computer simulations. Quantitative comparisons between experimental and computational results show promise for applying a FEM to the flow forming process. Experimental failure locations correlated well with computed damage gradients, and there was good agreement between the measured and predicted roller forces during forming.

Elizabeth Urig↗
DOE OSTI2025

Surrogate-constructed scalable-circuits adaptive variational quantum eigensolver in the Schwinger model

Inspired by recent advancements in simulating periodic systems on quantum computers, we develop an approach to further advance the simulation of these systems, named (SC) 2 -ADAPT-VQE. Our approach extends the scalable-circuits ADAPT-VQE framework, which builds an ansatz from a pool of coordinate-invariant operators defined for arbitrarily large, though not arbitrarily small, volumes. Our method uses a classically tractable “surrogate constructed” method to remove irrelevant operators from the pool, reducing the minimum size for which the scalable circuits are defined. Bringing together the scalable circuits and the surrogate constructed approaches forms the core of the (SC) 2 methodology. Our approach allows for a wider set of classical computations on small volumes, which can be used for a more robust extrapolation protocol. While developed in the context of lattice models, the surrogate construction portion is applicable to a wide variety of problems where information about the relative importance of operators in the pool is available. As an example, we use it to compute the properties of the Schwinger model—quantum electrodynamics for a single, massive fermion in 1 +1 dimensions—and show that our method can be used to accurately extrapolate to the continuum limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Jan-Dec 2021 IC Project: w20_earthquake (Yearly Report)

The simulations conducted under this IC project used an in-house code developed at Los Alamos National Laboratory called HOSS (Hybrid Optimization Software Suite) and its extended version called HOSS-FSIS (fluid-structure interaction solver). HOSS has been under development in LANL’s Geophysics Group since 2004 and was used to simulate the signal propagation in pieces of material with pre-existing defects and also to simulate the near source effects of above ground nuclear explosions. One of the key feature of HOSS is that it is superior to pure finite element method and discrete element method implementations because it allows for a natural transition from continuum to discontinuum when the material breaks or fractures.

58 GEOSCIENCES↗