Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Order priors for Bayesian network discovery with an application to malware phylogeny

Here, Bayesian networks have been used extensively to model and discover dependency relationships among sets of random variables. We learn Bayesian network structure with a combination of human knowledge about the partial ordering of variables and statistical inference of conditional dependencies from observed data. Our approach leverages complementary information from human knowledge and inference from observed data to produce networks that reflect human beliefs about the system as well as to fit the observed data. Applying prior beliefs about partial orderings of variables is an approach distinctly different from existing methods that incorporate prior beliefs about direct dependencies (or edges) in a Bayesian network. We provide an efficient implementation of the partial-order prior in a Bayesian structure discovery learning algorithm, as well as an edge prior, showing that both priors meet the local modularity requirement necessary for an efficient Bayesian discovery algorithm. In benchmark studies, the partial-order prior improves the accuracy of Bayesian network structure learning as well as the edge prior, even though order priors are more general. Our primary motivation is in characterizing the evolution of families of malware to aid cyber security analysts. For the problem of malware phylogeny discovery, we find that our algorithm, compared to existing malware phylogeny algorithms, more accurately discovers true dependencies that are missed by other algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Error Study of FlecSim Using the Method of Manufactured Solutions [Slides]

Error was calculated using exact solution chosen with the method of manufactured solutions. Hydrodynamics match expected order of accuracy for spatial error within 1.5% for all variables tested. Radiation transport matches expected order of temporal accuracy within 0.5% for all variables tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Corrosion Behavior of AZ31 Magnesium Alloy Produced by Continuous Variable Cross Section Direct Extrusion

In order to study the corrosion behavior of magnesium alloys of high performance in the simulated seawater by continuous variable cross section direct extrusion (CVCDE), an electrochemical method and weight loss method were used to study the difference between AZ31 magnesium alloys produced by CVCDE with different numbers of interim dies to examine the electrochemical corrosion behavior and corrosion rate at 3.5 wt.% NaCl solution in this paper. Scanning electron microscopy and x-ray diffraction were used to analyze the morphology of corrosion surfaces and the composition of the corrosion products of samples soaked for different times. The results show that after CVCDE with different numbers of interim dies, the size of the microstructure is refined gradually, and the grain orientation is also different; with the increase in the number of interim dies, the thermodynamic corrosion tendencies are different, and the corrosion resistance of samples extruded by conventional extrusion is relatively good; the charge transfer and material transfer process are the two main factors affecting the corrosion rate, which determines the direction of the electrochemical reaction; the corrosion form of the sample extruded by CVCDE is the local corrosion which is the filiform corrosion and pitting corrosion, and the corrosion product is mainly water-insoluble Mg(OH){sub 2}.

36 MATERIALS SCIENCE↗
DOE OSTI2019

On mass conservation and solvability of the discretized variable-density zero-Mach Navier-Stokes equations

This paper presents a study of a variable-density zero-Mach Navier-Stokes discretization and subtle difficulties associated with the pressure. A consistent compatibility condition between the temperature and velocity is first derived. It is shown that one can accurately and efficiently solve all primitive variables with the uniform order of accuracy of the specified discretization and time-marching schemes. It is shown that higher density ratios than those commonly reported are achievable, and it is, in part, traced to the solvability of the pressure Poisson equation. Numerical properties are assessed with the method of manufactured solutions, and several classical testing cases are carried out to validate the code, ranging from variable-density mixing problems to premixed flames with forcing.

42 ENGINEERING↗
DOE OSTI2022

A high accuracy/resolution spectral element/Fourier–Galerkin method for the simulation of shoaling non-linear internal waves and turbulence in long domains with variable bathymetry

A high-order hybrid continuous-Galerkin numerical method, designed for the simulation of non-linear, non -hydrostatic internal waves and turbulence in long computational domains with complex bathymetry, is presented. The spatial discretization in the non-periodic wave-propagating directions, utilizes the nodal spectral element method. Such a high-order element-based discretization allows the highly accurate representation of complex domain geometry along with the flexibility of concentrating resolution in areas of interest. Under the assumption of the normal-to-isobath propagation of non-linear internal waves, a third periodic direction is incorporated via a Fourier-Galerkin discretization. The distinct non-hydrostatic nature of non-linear internal waves and, any instabilities and turbulence therein, necessitates the numerically challenging solution of the pressure Poisson problem. A defining feature of this work is the application of a domain decomposition approach, combined with block-Jacobi/deflation-based preconditioning to the pressure Poisson problem. Such a combined approach is particularly suitable for the long high aspect-ratio complex domains of interest and enables the efficient high-accuracy reproduction of the non-hydrostatic dynamics of non-linear internal waves. Implementation details are also described in the context of the stability of the solver and its parallelization strategy. A series of benchmarks of increasing complexity demonstrate the robustness of the flow solver. The benchmarks culminate with the three-dimensional simulation of a convectively breaking mode-one non-linear internal wave over a realistic South-China-Sea bathymetric transect and background current/stratification profiles.

Deflation↗
DOE OSTI2017

The two-dimensional kinetic ballooning theory for ion temperature gradient mode in tokamak

The two-dimensional (2D) kinetic ballooning theory is developed for the ion temperature gradient mode in an up-down symmetric equilibrium (illustrated via concentric circular magnetic surfaces). The ballooning transform converts the basic 2D linear gyro-kinetic equation into two equations: (1) the lowest order equation (ballooning equation) is an integral equation essentially the same as that reported by Dong et al., [Phys. Fluids B 4, 1867 (1992)] but has an undetermined Floquet phase variable, (2) the higher order equation for the rapid phase envelope is an ordinary differential equation in the same form as the 2D ballooning theory in a fluid model [Xie et al., Phys. Plasmas 23, 042514 (2016)]. The system is numerically solved by an iterative approach to obtain the (phase independent) eigen-value. The new results are compared here to the two earlier theories. We find a strongly modified up-down asymmetric mode structure, and non-trivial modifications to the eigen-value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Symmetry-dependent ultrafast manipulation of nanoscale magnetic domains

Symmetry is a powerful concept in physics, but its applicability to far-from-equilibrium states is still being understood. Recent attention has focused on how far-from-equilibrium states lead to spontaneous symmetry breaking. Conversely, ultrafast optical pumping can be used to drastically change the energy landscape and quench the magnetic order parameter in magnetic systems. Here, we find a distinct symmetry-dependent ultrafast behaviour by use of ultrafast x-ray scattering from magnetic patterns with varying degrees of isotropic and anisotropic symmetry. After pumping with an optical laser, the scattered intensity reveals a radial shift exclusive to the isotropic component and exhibits a faster recovery time from quenching for the anisotropic component. These features arise even when both symmetry components are concurrently measured, suggesting a correspondence between the excitation and the magnetic order symmetry. Our results underline the importance of symmetry as a critical variable to manipulate the magnetic order in the ultrafast regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Phase Field Dislocation Dynamics (PFDD)

Phase Field Dislocation Dynamics (PFDD) is a phase field approach for tracking the motion and interactions of individual dislocations in materials, mostly face-centered cubic (fcc) metals. The phase field approach relies on scalar order parameters or phase field variables to track a quantity of interest. In PFDD, these phase field variables represent the passage of perfect dislocations across each active slip system, with one order parameter defined in each slip system. Partial dislocations are represented with linear combinations of the order parameters. The system is evolved through minimization of the total system energy, which is directly dependent on the phase field variables. All of the physics is described by the total energy functional that is derived from a thermodynamic basis, and the system kinetics are evolved using a Ginzburg-Landau (GL) equation.

Hunter, Abigail↗
DOE OSTI2019

Approach to qualification using E-PBF in-situ process monitoring in Ti-6Al-4V

Traditional design and qualification methodologies for parts manufactured by traditional methods are being applied to Additive Manufacturing (AM) without understanding the nuances of the machines. While mapping process variables and tracking build data is helpful, some variables such as build geometry, support structure, and part melt order have not been researched in depth. Changing these variables can result in significant variations in material properties and defect structure such that the process appears to be unreliable compared to traditional manufacturing. Hence, this research focuses on the need to understand the effects of overlooked variables such as melt order and nested geometry on the distribution of defects and bulk material properties in Ti-6Al-4 V alloy builds manufactured using the Arcam AB ® electron beam powder bed fusion process 1 . This work collected and analyzed process log data and near infrared (NIR) images for every layer to correlate trends in porosity formation and mechanical performance. The location of pores, while naturally stochastic, is greatly influenced by the cross-sectional area as detected by NIR images and correlates with the failure sites from uniaxial testing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine-learning-based investigation of the variables affecting summertime lightning occurrence over the Southern Great Plains

Lightning is affected by many factors, many of which are not routinely measured, well understood, or accounted for in physical models. Several commonly used machine learning (ML) models have been applied to analyze the relationship between Atmospheric Radiation Measurement (ARM) data and lightning data from the Earth Networks Total Lightning Network (ENTLN) in order to identify important variables affecting lightning occurrence in the vicinity of the Southern Great Plains (SGP) ARM site during the summer months (June, July, August and September) of 2012 to 2020. Testing various ML models, we found that the random forest model is the best predictor among common classifiers. When convective clouds were detected, it predicts lightning occurrence with an accuracy of 76.9 % and an area under the curve (AUC) of 0.850. Using this model, we further ranked the variables in terms of their effectiveness in nowcasting lightning and identified geometric cloud thickness, rain rate and convective available potential energy (CAPE) as the most effective predictors. The contrast in meteorological variables between no-lightning and frequent-lightning periods was examined for hours with CAPE values conducive to thunderstorm formation. Besides the variables considered for the ML models, surface variables and mid-altitude variables (e.g., equivalent potential temperature and minimum equivalent potential temperature, respectively) have statistically significant contrasts between no-lightning and frequent-lightning hours. For example, the minimum equivalent potential temperature from 700 to 500 hPa is significantly lower during frequent-lightning hours compared with no-lightning hours. Finally, a notable positive relationship between the intracloud (IC) flash fraction and the square root of CAPE ($\sqrt{CAPE}$) was found, suggesting that stronger updrafts increase the height of the electrification zone, resulting in fewer flashes reaching the surface and consequently a greater IC flash fraction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

High-order gradients with the shifted boundary method: An embedded enriched mixed formulation for elliptic PDEs

Here, we propose an extension of the embedded boundary method known as “shifted boundary method” to elliptic diffusion equations in mixed form (e.g., Darcy flow, heat diffusion problems with rough coefficients, etc.). Our aim is to obtain an improved formulation that, for linear finite elements, is at least second-order accurate for both flux and primary variable, when either Dirichlet or Neumann boundary conditions are applied. Following previous work of Nishikawa and Mazaheri in the context of residual distribution methods, we consider the mixed form of the diffusion equation (i.e., with Darcy-type operators), and introduce an enrichment of the primary variable. This enrichment is obtained exploiting the relation between the primary variable and the flux variable, which is explicitly available at nodes in the mixed formulation. The proposed enrichment mimics a formally quadratic pressure approximation, although only nodal unknowns are stored, similar to a linear finite element approximation. We consider both continuous and discontinuous finite element approximations and present two approaches: a non-symmetric enrichment, which, as in the original references, only improves the consistency of the overall method; and a symmetric enrichment, which enables a full error analysis in the classical finite element context. Combined with the shifted boundary method, these two approaches are extended to high-order embedded computations, and enable the approximation of both primary and flux (gradient) variables with second-order accuracy, independently on the type of boundary conditions applied. We also show that the primary variable is third-order accurate, when pure Dirichlet boundary conditions are embedded.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Fluid Dynamics Effects on Microstructure Prediction in Single Laser Tracks for Additive Manufacturing

The Laser Powder Bed Fusion Additive Manufacturing (LPBFAM) is one of the most important processes for the production of lightweight, cost-effective, complex, and high-performance enduse parts. At present, the cost and time associated with LPBFAM process development is very high due to a lack of fundamental process understanding. In this project, a multi-physics model was developed on a highly parallel open-source code, Truchas, with the ultimate goal of providing experimentally validated process maps for tailoring microstructure to achieve desired performance for LPBFAM. As a critical step towards fully LPBFAM modeling, modeling of single-track laser fusion (STLF) were conducted. Multi-physics simulations were conducted using Truchas for STLF by considering heat transfer, phase-change, fluid dynamics, surface tension phenomena, and evaporation. In order to assess the effect of the fluid dynamics on the solidification and ensuing microstructure, a heat transfer-and-solidification-only (HTS) model and a fully coupled heat-transfer-solidification and fluid-dynamics (HTS-FD) model were considered. In the HTS_FD model, the fluid flow was considered to be laminar while the molten alloy surface is assumed to be flat and non-deformable. This study is one the first attempt to understand the effect of the fluid flow on microstructure in STLF and LPBFAM. The results show that the fluid flow affects the solidification and ensuing microstructure in STLF. In order to validate Truchas for STLF and LPBFAM modeling, experimental data, which was obtained at GE Global Research (GEGR) for liquid pool shape and microstructure, were compared with those from numerical simulation results. For the keyhole regime case, numerical simulation results indicate that the melt-pool shape was typical to that of the conduction case with very large deviations from the measured melt-pool depths. Two sensitivity studies were conducted by varying the evaporation flux and the surface tension coefficient at fixed laser absorptivities of 0.5 and 0.3, respectively. The minimum value for the overall combined error, between the calculated values with the HTS_FD model and measured values for both the width and height of the melt-pool, was attained for a surface tension coefficient of -0.75e-4 N/m. An analytical solidification model for the columnar-to-equiaxed transition (CET) in rapid solidification was used to assess the microstructure variation within the melt-pool. From numerical simulation results, thermal gradient, G, and solidification velocity, V, were obtained in order to predict the microstructure type (e.g., dendritic, cellular). For the fluid-flow HTS_FD model, G was found to exhibit a maximum at the extremity of the solidified pool (i.e., at the free surface). By contrast, for HTS simulations, G was found to exhibit a maximum around the entire edge of the solidified pool. For the HTS_FD simulations, the minimum values of the cooling rate, GV, were found to be approximately half than their corresponding values for HTS simulations. By contrast with the min GV values, the maximum values for GV were found to occur for the HTS_FD simulations. Thus, HTS_FD simulations were found to exhibit a wider range of cooling rates than the HTS simulations. Concerning the solidification maps, the variation of the solidification velocity, V, as function of the thermal gradient, G, was obtained. It was found that the fluid flow model results (HTS_FD) exhibited an increased spread in the V(G) variation within the melt pool with respect to the HTS model results (without the fluid flow). A preliminary correlation model for the primary dendrite arm spacing (PDAS) based on a power law dependence on thermal gradient and solidification velocity was used to estimate PDAS from HTS and HTS_FD simulations. For low laser powers and low laser speeds, the PDAS obtained 6 with the fluid dynamics model (HTS_FD) was larger by more than 30% with respect to the PDAS calculated with the simple HTS model. In the second part of the project, fifty-seven simulations were conducted in order to obtain data on microstructure variables that can be used for process map development. In order to cover the entire processing space, thirty-three cases were selected STLF process simulations at six power levels and six laser scanning speeds. The 57 STLF process simulations were conducted as follows: 33 simulations with the heat-transfer-only (HTS) model and 24 simulations with the fluid flow model (HTS_FD). The solidification map data shows that for all simulations, a columnar dendritic microstructure would be expected. It was found that the minimum, average, and maximum thermal gradient exhibit exponential variations with respect to a process variable defined as the ratio between the power and square root of scan speed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analysis of Long-Term Quality Control Data for a 137 Cs Dosimetry Calibration Source

Strict quality assurance programs are required for many radiological applications, but these seldom exist for verifying dosimetry calibration sources. After initial characterization of a dosimetry calibration facility, quality control procedures are recommended to ensure the early detection of any changes or malfunctions. These also result in refined knowledge about average dose rate and experimental variations in dose delivery. This paper describes the implementation of a phase I quality control protocol for a 137 Cs dosimetry calibration source and includes an analysis of the resulting data collected over a 24-mo period. During this time, substantial data was collected to establish trial control limits. Air kerma rate measurements were obtained using an ion chamber and were adjusted for decay, corrected for ambient temperature, pressure and humidity, and then analyzed using quality control charts. Three variations of rational subgrouping methods were used in order to find assignable causes of error, and Nelson's Rules were followed to detect any non-random statistical variations. Measurements were subgrouped according to same-day measurements in order to detect positional errors as well as atmospheric correction errors. Additionally, measurements were subgrouped according to analogous experimental setups in order to detect failure in equipment or incorrect settings. Both were analyzed using the X-bar and R chart method. Similarly, individuals and moving ranges charts were used to carefully examine each position in order to observe any situational errors that may occur which include timing, positional, or interference errors. Each method was successful in identifying unique out-of-control data points that occurred during the phase I application of forming control limits. Furthermore, over the 24-mo period, enough data points were deemed in-control to establish reliable trial limits. Future experiments will include the phase II application of gaining more reliable measurements in order to fine-tune the limits, as well as performing a designed experiment, where variables are purposefully changed in order to test the variation of the data.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2019

On the Kernel of a Two-Point Problem for a Partial Differential Equation of the Second Order in Time

We study the problem for a homogeneous partial differential equation of the second order with respect to time with given homogeneous two-point conditions in this variable and, in general, of the infinite order in the other (space) variable. It is proved that the analyzed problem possesses solely the trivial solution if the characteristic determinant is not identically equal to zero. In the case where the set of zeros of the characteristic determinant of this problem is nonempty, we propose a method for the construction of nontrivial solutions of the problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Effect of local chemical order on the irradiation-induced defect evolution in CrCoNi medium-entropy alloy

High- (and medium-) entropy alloys have emerged as potentially suitable structural materials for nuclear applications, particularly as they appear to show promising irradiation resistance. Recent studies have provided evidence of the presence of local chemical order (LCO) as a salient feature of these complex concentrated solid-solution alloys. However, the influence of such LCO on their irradiation response has remained uncertain thus far. Here, in this work, we combine ion irradiation experiments with large-scale atomistic simulations to reveal that the presence of chemical short-range order, developed as an early stage of LCO, slows down the formation and evolution of point defects in the equiatomic medium-entropy alloy CrCoNi during irradiation. In particular, the irradiation-induced vacancies and interstitials exhibit a smaller difference in their mobility, arising from a stronger effect of LCO in localizing interstitial diffusion. This effect promotes their recombination as the LCO serves to tune the migration energy barriers of these point defects, thereby delaying the initiation of damage. These findings imply that local chemical ordering may provide a variable in the design space to enhance the resistance of multi-principal element alloys to irradiation damage.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme

We present here a Lagrangian discontinuous Galerkin (DG) hydrodynamic method that is up to third-order accurate using subcell mesh stabilization (SMS) for compressible flows on quadratic meshes in two-dimensional (2D) Cartesian coordinates. Similar to the second-order accurate Lagrangian DG method with linear meshes, the physical evolution equations for the specific volume, velocity, and specific total energy are discretized using a modal DG method with Taylor series polynomials. The Riemann velocity at the vertices of a curvilinear cell, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem. Curvilinear cells (e.g., quadratic quadrilateral meshes in this work) have many deformational degrees of freedom, and with these cells, they can deform in unphysical ways. Likewise, the Riemann solution at an edge vertex differs from the one at the corner of a cell. With SMS, each quadratic quadrilateral cell is decomposed into four quadrilateral subcells, that move in a Lagrangian manner. The edge vertex is surrounded by four subcells so that it is similar to the vertex at the cell corner. SMS can detect inconsistent density fields between the cell and subcells. The difference between these two density fields is used to correct the stress (pressure) input to the Riemann solver. This SMS scheme enables stable mesh motion and accurate solutions in the context of a Lagrangian high-order DG method that is up to third-order with quadratic cells. We also present effective limiting strategies that ensure monotonicity of the primitive variables with the high-order DG method. This Lagrangian DG hydrodynamic method with SMS conserves mass, momentum, and total energy. A suite of test problems are calculated to demonstrate the designed order of accuracy (up to third-order accurate) of this method, and that the Lagrangian DG method using SMS preserves cylindrical symmetry on 1D radial flow problems with an equal-angle polar quadratic mesh.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗