Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Partially explicit”

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
NASA NTRS1979

Generation of linear dynamic models from a digital nonlinear simulation

The results and methodology used to derive linear models from a nonlinear simulation are presented. It is shown that averaged positive and negative perturbations in the state variables can reduce numerical errors in finite difference, partial derivative approximations and, in the control inputs, can better approximate the system response in both directions about the operating point. Both explicit and implicit formulations are addressed. Linear models are derived for the F 100 engine, and comparisons of transients are made with the nonlinear simulation. The problem of startup transients in the nonlinear simulation in making these comparisons is addressed. Also, reduction of the linear models is investigated using the modal and normal techniques. Reduced-order models of the F 100 are derived and compared with the full-state models.

Daniele, C. J.↗
NASA NTRS1983

Adaptive grid methods for partial differential equations

A number of techniques for constructing adaptive mesh generators for use in solving partial differential equations are reviewed in this paper. Techniques reviewed include methods based on steady grid generation schemes and those which are explicitly designed to determine grid speeds in a time-dependent or space-marching problem. Results for candidate methods are included and suggestions for areas of future research are suggested.

Anderson, D. A.↗
NASA NTRS1992

Bleed-boundary conditions for numerically simulated mixed-compression supersonic inlet flow

A numerical study was performed to investigate boundary conditions for bleed openings that are needed in the numerical simulation of critical flow through turbojet inlets with bleed. Nine different boundary conditions (BCs) were devised and examined. Of the BCs examined, three were able to stabilize the terminal shock near the inlet's throat. This numerical study was based on the ensemble-averaged, unsteady, three-dimensional, compressible Navier-Stokes equations closed by the Baldwin-Lomax algebraic turbulence model. Solutions to the Navier-Stokes equations were obtained by using a hybrid implicit-explicit method (Beam and Warming/MacCormack). The multiblock grid system used in the solution procedure was generated by solving a set of elliptic partial differential equations. Solutions obtained were compared with available experimental data.

Chyu, W. J.↗
DOE OSTI2025

Characterization of admissible quasisymmetries

We solve ‘half’ the problem of finding three-dimensional quasisymmetric magnetic fields that do not necessarily satisfy magnetohydrostatic force balance. This involves determining which hidden symmetries are admissible as quasisymmetries, and then showing explicitly how to construct quasisymmetric magnetic fields given an admissible symmetry. The admissibility conditions take the form of a system of overdetermined nonlinear partial differential equations involving second derivatives of the symmetry's infinitesimal generator.

Burby, J. W. (ORCID:0000000317727647)↗
NASA NTRS1983

Second-order explicit finite-difference methods for transient-flow analysis

Three second-order accurate numerical methods - MacCormack's method, Lambda scheme and Gabutti scheme - are introduced to solve the quasi-linear, hyperbolic partial differential equations describing transient flows in closed conduits. The details of these methods and the treatment of boundary conditions are presented and the results computed by using these methods for a typical piping system are compared. It is shown that for the same accuracy, second-order methods require considerably lesser number of computational nodes and computer time as compared to those required by the first-order methods.

Chaudhry, M. H.↗
NASA NTRS1994

Sloshing dynamics modulated angular momentum fluctuations actuated by gravity gradient acceleration associated with spacecraft slew motion

The mathematical formulation of orbital spacecraft sloshing dynamics for a partially filled cryogenic superfluid liquid helium II in a dewar container actuated by the gravity gradient acceleration associated with slew motion is studied. The Advanced X-Ray Astrophysics Facility-Spectroscopy (AXAF-S) spacecraft is chosen as a practical example in this study. Explicit mathematical expressions which manage orbital galaxy gradient acceleration with slew motion which is acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics is based on the non-inertia frame spacecraft bound coordinates and the solution of time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. This study discloses the capillary effect of sloshing dynamics governed liquid-vapor interface fluctuations, angular momentum and moment fluctuations of fluid system, and also buble mass center fluctuations driven by the gravity gradient acceleration associated with slew motion which affects the stability of the orbital spacecraft fluid system in a microgravity environment.

Hung, R. J.↗
DOE OSTI2024

Using Machine Learning to Predict Cloud Turbulent Entrainment–Mixing Processes

Different turbulent entrainment–mixing mechanisms between clouds and environment are essential to cloud–related processes; however, accurate representation of entrainment–mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment–mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment–mixing. Specifically, number concentration after entrainment–mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment–mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment–mixing is negatively correlated with liquid water content affected by entrainment–mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment–mixing in weather/climate models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Inclusion of exact exchange in the noniterative partial-differential-equation method of electron-molecule scattering - Application to e-N2

The noniterative partial-differential-equation (PDE) approach to electron-molecule scattering of Onda and Temkin (1983) is modified to account for the effects of exchange explicitly. The exchange equation is reduced to a set of inhomogeneous equations containing no integral terms and solved noniteratively in a difference form; a method for propagating the solution to large values of r is described; the changes in the polarization potential of the original PDE method required by the inclusion of exact static exchange are indicated; and the results of computations for e-N2 scattering in the fixed-nuclei approximation are presented in tables and graphs and compared with previous calculations and experimental data. Better agreement is obtained using the modified PDE method.

Weatherford, C. A.↗
NASA NTRS2020

A 3D Model to Predict Explicit Morphologies and Volume Fraction of Lack-of-Fusion Pores Generated in Selective Laser Melting Processes

The performance of an additively manufactured (AM) component is dependent on the distribution of process-induced defects in addition to the complex microstructure of the material, surface roughness of the component and the process-induced residual stresses. For instance, it has been well demonstrated that lack-of-fusion (LoF) pores produced in the selective laser melting (SLM) AM process can significantly limit the fatigue performance of the material. Although two-dimensional (2D) models exist to predict the 2D profiles of LoF pores, the 2D pore profiles cannot be directly inserted into a three-dimensional (3D) microstructure domain that is output from several prevailing process simulation packages. A few commercial packages that simulate the SLM process can predict LoF pores in a 3D domain, but the morphologies of LoF pores are voxelated and hence do not capture sharp corners of the pores, thereby obviating their use in fatigue crack initiation studies. In order to address the aforementioned gaps, a high-fidelity model that predicts not only the volume fraction, but also the explicit 3D morphologies and spatial distributions of LoF pores has been developed using a computer aided design-based environment. The model has been partially validated for Ti-6Al-4V alloy by comparing the predictions of the volume fraction of LoF pores predicted by the model with experimental data obtained from the literature. Absolute error in predicted volume fraction of LoF pores varied between 5.16% and 1.87% for energy density values between 13 J/mm3 and 45 J/mm3 where a significant amount (over 3% volume fraction) of LoF porosity was measured. The absolute error was within 1.87% for energy density values greater than 45 J/mm3.

Saikumar R. Yeratapally↗
DOE OSTI2022

Time-explicit Darwin PIC algorithm

A new approach to Darwin particle-in-cell plasma simulation is described. Using a finite-element approach, the vector potential is assured to be exactly solenoidal. This allows writing the system of multi-species particles as an action-at-a-distance Hamiltonian system. Applying recently-developed explicit symplectic methods to this gives a time-explicit algorithm with all desired properties. Elliptic systems for the electrostatic and magneto-static portions of the problem are inverted efficiently by an algebraic multi-grid algorithm. The algorithm is implemented in a two-dimensional Cartesian-geometry code and tested by application to Weibel instability and to Alfvén wave dynamics. Results show the effectiveness of this approach. Extensions to arbitrary meshes and to a partially time-implicit scheme are briefly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ordinary and exotic mesons in the extended Linear Sigma Model

The extended Linear Sigma Model (eLSM) is a hadronic model based on the global symmetries of QCD and the corresponding explicit, anomalous, and spontaneous breaking patterns. In its basic three-flavor form, its mesonic part contains the dilaton/glueball as well as the nonets of pseudoscalar, scalar, vector, and axial–vector mesons, thus chiral symmetry is linearly realized. In the chiral limit and neglecting the chiral anomaly, only one term – within the dilaton potential – breaks dilatation invariance, and all terms are chirally symmetric. Spontaneous symmetry breaking is implemented by a generalization of the Mexican-hat potential, with explicit symmetry breaking responsible for its tilting. The overall mesonic phenomenology up to ~2 GeV is in agreement with the PDG compilation of masses and partial and total decay widths. The eLSM was enlarged in a straightforward way to include other conventional quark–antiquark nonets (pseudovector and orbitally excited vector mesons, tensor and axial-tensor mesons, radially excited (pseudo)scalar mesons, etc.), as well as two nonets of hybrid mesons, the lightest one with exotic quantum numbers J $\mathscr{PC}$ = 1 -+ not allowed for $\bar{q}q$ objects, such as the resonance π 1 (1600) and the recently discovered η 1 (1855) . In doing so, different types of chiral multiplets are introduced: heterochiral and homochiral multiplets, which differ in the way they transform under chiral transformations. Moreover, besides the scalar glueball that is present from the beginning as dilaton, other glueballs, the tensor, the pseudoscalar and the vector glueballs were coupled to the eLSM: the scalar resonance ƒ 0 (1710) turns out to be mostly gluonic, the tensor glueball couples strongly to vector mesons, and the pseudoscalar glueball couples sizably to $ππη'$ and can be assigned to X (2370) or X (2600). In all cases above, masses and decays can be analyzed allowing for a better understanding of both conventional and non-conventional mesons: whenever data are available, a comparison is performed and, when this is not the case, predictions of decay widths and decay ratios are outlined. The eLSM contains chiral partners on an equal footing and is therefore well suited for studies of chiral symmetry restoration at nonzero temperature and densities: this is done by coupling it to the Polyakov loop. Finally, the QCD phase diagram and the location of the critical endpoint were investigated within this framework.

Chiral symmetry↗
DOE OSTI2025

Pion-nucleon scattering in baryon chiral perturbation theory combined with the 1/𝑁 𝑐 expansion

This work implements the combined baryon chiral perturbation theory (BChPT) and 1/𝑁 𝑐 expansions for pion-nucleon elastic scattering. The effective theory is based on the baryon sector dynamical spin-flavor 𝑆⁢𝑈⁡(4) symmetry emergent in the large 𝑁 𝑐 limit, whose breaking is controlled by the 1/𝑁 𝑐 expansion. The noncommutativity of the chiral and 1/𝑁 𝑐 expansions in unitarity corrections (loops) requires a linking of both expansions. As it was shown in the case of baryon masses and currents, the natural linking is the 𝜉 expansion, in which 𝒪⁡(𝑝) = 𝒪⁡(1/𝑁 𝑐 ) = 𝒪⁡(𝜉). The spin-flavor symmetry requires that the ground state baryons span an 𝑆⁢𝑈⁡(4) symmetric irreducible representation which implies that in particular 𝑁 and Δ are active degrees of freedom in the effective theory. The scattering amplitude is expanded to the next-to-next-to leading order in the 𝜉 expansion, corresponding to the one-loop contributions with the leading-order Lagrangian. The results are given for generic 𝑁 𝑐 in order to demonstrate the consistency of the framework. The spin-flavor symmetry plays a central role in maintaining the consistency of the effective theory with respect to the 1/𝑁 𝑐 expansion. This consistency manifests itself in an improvement in the convergence of the low energy expansion with respect to the case of the ordinary BChPT without an explicit dynamical Δ, which is known to be inconsistent with the constraints of 𝑁 𝑐 scaling. Fits to the 𝜋⁢𝑁 → 𝜋⁢𝑁, 𝑆, 𝑃, and 𝐷 partial wave amplitudes from the SAID data base are finally used to test the framework and to determine the energy range of its applicability.

Jayakodige, D. [Hampton Univ., Hampton, VA (United↗
NASA NTRS1989

Modified Chebyshev pseudospectral method with O(N exp -1) time step restriction

The extreme eigenvalues of the Chebyshev pseudospectral differentiation operator are O(N exp 2) where N is the number of grid points. As a result of this, the allowable time step in an explicit time marching algorithm is O(N exp -2) which, in many cases, is much below the time step dictated by the physics of the partial differential equation. A new set of interpolating points is introduced such that the eigenvalues of the differentiation operator are O(N) and the allowable time step is O(N exp -1). The properties of the new algorithm are similar to those of the Fourier method. The new algorithm also provides a highly accurate solution for non-periodic boundary value problems.

Kosloff, Dan↗
NASA NTRS1993

A modified Chebyshev pseudospectral method with an O(N exp -1) time step restriction

The extreme eigenvalues of the Chebyshev pseudospectral differentiation operator are O(N exp 2) where N is the number of grid points. As a result of this, the allowable time step in an explicit time marching algorithm is O(N exp -2) which, in many cases, is much below the time step dictated by the physics of the partial differential equation. A new set of interpolating points is introduced such that the eigenvalues of the differentiation operator are O(N) and the allowable time step is O(N exp -1). The properties of the new algorithm are similar to those of the Fourier method. The new algorithm also provides a highly accurate solution for non-periodic boundary value problems.

Kosloff, Dan↗
NASA NTRS1981

The nonlinear modified equation approach to analyzing finite difference schemes

The nonlinear modified equation approach is taken in this paper to analyze the generalized Lax-Wendroff explicit scheme approximation to the unsteady one- and two-dimensional equations of gas dynamics. Three important applications of the method are demonstrated. The nonlinear modified equation analysis is used to (1) generate higher order accurate schemes, (2) obtain more accurate estimates of the discretization error for nonlinear systems of partial differential equations, and (3) generate an adaptive mesh procedure for the unsteady gas dynamic equations. Results are obtained for all three areas. For the adaptive mesh procedure, mesh point requirements for equal resolution of discontinuities were reduced by a factor of five for a 1-D shock tube problem solved by the explicit MacCormack scheme.

Klopfer, G. H.↗
DOE OSTI2024

Node-degree aware edge sampling mitigates inflated classification performance in biomedical random walk-based graph representation learning

Motivation: Graph representation learning is a family of related approaches that learn low-dimensional vector representations of nodes and other graph elements called embeddings. Embeddings approximate characteristics of the graph and can be used for a variety of machine-learning tasks such as novel edge prediction. For many biomedical applications, partial knowledge exists about positive edges that represent relationships between pairs of entities, but little to no knowledge is available about negative edges that represent the explicit lack of a relationship between two nodes. For this reason, classification procedures are forced to assume that the vast majority of unlabeled edges are negative. Existing approaches to sampling negative edges for training and evaluating classifiers do so by uniformly sampling pairs of nodes. Results: We show here that this sampling strategy typically leads to sets of positive and negative examples with imbalanced node degree distributions. Using representative heterogeneous biomedical knowledge graph and random walk-based graph machine learning, we show that this strategy substantially impacts classification performance. If users of graph machine-learning models apply the models to prioritize examples that are drawn from approximately the same distribution as the positive examples are, then performance of models as estimated in the validation phase may be artificially inflated. We present a degree-aware node sampling approach that mitigates this effect and is simple to implement. Availability and implementation: Our code and data are publicly available at https://github.com/monarch-initiative/negativeExampleSelection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Extending Symmetry-Preserving Attention Networks (SPANet) for Jet Assignment in Fully Hadronic \(t\bar{t}\) Events

Fully hadronic \(t\bar{t}\) reconstruction requires assigning reconstructed jets to the two top-quark decay branches, \(t\bar{t}\to(Wb)(Wb)\to(qqb)(qqb)\). This is a combinatorially large, symmetry-rich, and frequently underconstrained set-assignment problem due to detector effects and partial reconstructability. SPANet provides a strong symmetry-preserving baseline for this task, but its standard inference eagerly commits to a single hypothesis and does not explicitly treat branch-level reconstructability decisions. We present a two-stage extension of the SPANet pipeline for fully hadronic \(t\bar{t}\) events: (i) SPANet is modified into a proposal model that generates a shortlist of candidate jet assignments; (ii) a custom dual-head set transformer is trained to re-rank these candidates and identify reconstructible branches. Our extension gives modest improvements in candidate-selection efficiency, with SPANet achieving \(69.9\%\) and our extension achieving \(71.1\%\). Finally, an oracle study of the generated shortlist shows a substantial upper bound of \(88.5\%\) on the evaluated fully reconstructable subset, suggesting that the shortlist contains significant residual information that is not fully exploited.

Lisondodi Tada, Mateo [Puerto Rico U., Mayaguez] (↗
NASA NTRS1972

An analytical study of reduced-gravity liquid reorientation using a simplified marker and cell technique

A computer program called HOPI was developed to predict reorientation flow dynamics, wherein liquids move from one end of a closed, partially filled, rigid container to the other end under the influence of container acceleration. The program uses the simplified marker and cell numerical technique and, using explicit finite-differencing, solves the Navier-Stokes equations for an incompressible viscous fluid. The effects of turbulence are also simulated in the program. HOPI can consider curved as well as straight walled boundaries. Both free-surface and confined flows can be calculated. The program was used to simulate five liquid reorientation cases. Three of these cases simulated actual NASA LeRC drop tower test conditions while two cases simulated full-scale Centaur tank conditions. It was concluded that while HOPI can be used to analytically determine the fluid motion in a typical settling problem, there is a current need to optimize HOPI. This includes both reducing the computer usage time and also reducing the core storage required for a given size problem.

Betts, W. S., Jr.↗