Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “formulation systems”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Numerical simulation using vorticity-vector potential formulation

An accurate and efficient computational method is needed for three-dimensional incompressible viscous flows in engineering applications. On solving the turbulent shear flows directly or using the subgrid scale model, it is indispensable to resolve the small scale fluid motions as well as the large scale motions. From this point of view, the pseudo-spectral method is used so far as the computational method. However, the finite difference or the finite element methods are widely applied for computing the flow with practical importance since these methods are easily applied to the flows with complex geometric configurations. However, there exist several problems in applying the finite difference method to direct and large eddy simulations. Accuracy is one of most important problems. This point was already addressed by the present author on the direct simulations on the instability of the plane Poiseuille flow and also on the transition to turbulence. In order to obtain high efficiency, the multi-grid Poisson solver is combined with the higher-order, accurate finite difference method. The formulation method is also one of the most important problems in applying the finite difference method to the incompressible turbulent flows. The three-dimensional Navier-Stokes equations have been solved so far in the primitive variables formulation. One of the major difficulties of this method is the rigorous satisfaction of the equation of continuity. In general, the staggered grid is used for the satisfaction of the solenoidal condition for the velocity field at the wall boundary. However, the velocity field satisfies the equation of continuity automatically in the vorticity-vector potential formulation. From this point of view, the vorticity-vector potential method was extended to the generalized coordinate system. In the present article, we adopt the vorticity-vector potential formulation, the generalized coordinate system, and the 4th-order accurate difference method as the computational method. We present the computational method and apply the present method to computations of flows in a square cavity at large Reynolds number in order to investigate its effectiveness.

Tokunaga, Hiroshi↗
NASA NTRS1990

Compact scheme for systems of equations applied to fundamental problems of mechanics of continua

Compact scheme formulation was used in the treatment of boundary conditions for a system of coupled diffusion and Poisson equations. Models and practical solutions of specific engineering problems arising in solid mechanics, chemical engineering, heat transfer and fuid mechanics are described and analyzed for efficiency and accuracy. Only 2-D cases are discussed and a new method of numerical treatment of boundary conditions common in the fundamental problems of mechanics of continua is presented.

Klimkowski, Jerzy Z.↗
NASA NTRS2001

Institutionalizing Lessons Learned

The NASA Integrated Action Team (NIAT) was formed by the NASA Administrator in March 2000. The purpose of this team was to identify the actions that NASA must take to address systemic findings reported in 4 different anomaly investigations. Team membership represented senior managers from all the field centers and NASA Headquarters. NIAT report addressed 165 findings and developed 17 action plans that are described in five themes: people and teams, technology, risk, formulation rigor, and communications. The NIAT actions present a systems solution for strengthening formulation and implementation of programs and improving the environment for their support. NIAT results included: enhancing success by avoiding failures that could have been prevented through good planning and sound practice; ensuring that prudent risks do not compromise safety; and ensuring that mission risks are objectively assessed, appropriately mitigated and consciously accepted by the program team and customers. Definitions of Faster, Better, Cheaper and Success Criteria were also developed and included as part of the NIAT report. As a result of the NIAT report, program and project management process changes were incorporated into NASA's quality system documentation, including NPG 7120.513, "NASA Program and Project Management Processes and Requirements. This paper describes the NIAT results and the resulting updates to NPG 7120.5 that keep this program and project management description a living process.

McBrayer, Robert O.↗
DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗
DOE OSTI2016

A difference-equation formalism for the nodal domains of separable billiards

Recently, the nodal domain counts of planar, integrable billiards with Dirichlet boundary conditions were shown to satisfy certain difference equations in Samajdar and Jain (2014). The exact solutions of these equations give the number of domains explicitly. For complete generality, we demonstrate this novel formulation for three additional separable systems and thus extend the statement to all integrable billiards.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Variational, stable, and self-consistent coupling of 3D electromagnetics to 1D transmission lines in the time domain

This work presents a new multiscale method for coupling the 3D Maxwell's equations to the 1D telegrapher's equations. While Maxwell's equations are appropriate for modeling complex electromagnetics in arbitrary-geometry domains, simulation cost for many applications (e.g. pulsed power) can be dramatically reduced by representing less complex transmission line regions of the domain with a 1D model. By assuming a transverse electromagnetic (TEM) ansatz for the solution in a transmission line region, we reduce the Maxwell's equations to the telegrapher's equations. Here, we propose a self-consistent finite element formulation of the fully coupled system that uses boundary integrals to couple between the 3D and 1D domains and supports arbitrary unstructured 3D meshes. Additionally, by using a Lagrange multiplier to enforce continuity at the coupling interface, we allow for an absorbing boundary condition to also be applied to non-TEM modes on this boundary. We demonstrate that this feature reduces non-physical reflection and ringing of non-TEM modes off of the coupling boundary. By employing implicit time integration, we ensure a stable coupling, and we introduce an efficient method for solving the resulting linear systems. We demonstrate the accuracy of the new method on two verification problems, a transient O-wave in a rectilinear prism and a steady-state problem in a coaxial geometry, and show the efficiency and weak scalability of our implementation on a cold test of the Z-machine MITL and post-hole convolute.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Adsorption of Mixed Micelles of Polysorbate 80 and Oleic Acid to the Air–Water Interface

The solubilization of long chain amphiphiles with limited solubility into mixed micelles composed of highly soluble surfactants plays a crucial role in modulating the stability and functionality of formulations in pharmaceutical and food systems. Subsequently, the mixed micelles adsorb from solution to the air−water surface, defining the transport mechanism of the insoluble amphiphiles to the interface. We use X-ray reflectivity to measure the composition and provide insight into the structure of these mixed monolayers, and, alongside interfacial tension measurements, we provide an understanding of how the bulk composition determines the surface composition and the dynamics of tension reduction. We use a model system consisting of an insoluble fatty acid, oleic acid (OA), and a soluble micelle-forming surfactant, polysorbate 80 (PS80). PS80 forms spherical micelles, and, above the critical micelle concentration (cmc), OA is readily solubilized inside the micelles. We show that the adsorption of PS80/OA mixed micelles to the air−water interface rapidly reduces the tension and lowers the equilibrium tension in proportion to the OA concentration. X-ray reflectivity data, fit using Parratt-slab models, quantitively demonstrates that the monolayers become enriched with OA, and we show how the OA intercalates into the PS80 monolayers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A two-stage service restoration method for electric power distribution systems

Improving the reliability of power distribution systems is critically important for both utilities and customers. This calls for an efficient service restoration module within a distribution management system to support the implementation of self-healing smart grid networks. Although the emerging smart grid technologies, including distributed generators (DGs) and remote-controlled switches, enhance the self-healing capability and allow faster recovery, they still pose additional complexity to the service restoration problem, especially under cold load pickup (CLPU) conditions. Herein, a novel two-stage restoration framework is proposed to generate a restoration solutions with a sequence of control actions. The first stage generates a restoration plan that supports both the traditional service restoration using feeder reconfiguration and the grid-forming DG-assisted intentional islanding methods. The second stage generates an optimal sequence of switching operations to bring the outaged system quickly to the final restored configuration. The problem is formulated as a mixed-integer linear program that incorporates system connectivity, operating constraints, and the CLPU models. It is demonstrated that on using a multi-feeder test case, the proposed framework is effective in utilizing all available resources to quickly restore the service and generate an optimal sequence of switching actions to be used by the operator to reach the desired optimal configuration.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

Implicit versus explicit momentum relaxation time solution for semiconductor nanowires

We discuss the necessity of the exact implicit Momentum Relaxation Time (MRT) solution of the Boltzmann transport equation in order to achieve reliable carrier mobility results in semiconductor nanowires. Firstly, the implicit solution for a 1D electron gas with a isotropic bandstructure is presented resulting in the formulation of a simple matrix system. Using this solution as a reference, the explicit approach is demonstrated to be inaccurate for the calculation of inelastic anisotropic mechanisms such as polar optical phonons, characteristic of III-V materials. Its validity for elastic and isotropic mechanisms is also evaluated. Finally, the implications of the MRT explicit approach inaccuracies on the total mobility of Si and III-V NWs are studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Action principle for Coulomb collisions in plasmas

In this study, an action principle for Coulomb collisions in plasmas is proposed. Although no natural Lagrangian exists for the Landau-Fokker-Planck equation, an Eulerian variational formulation is found considering the system of partial differential equations that couple the distribution function and the Rosenbluth-MacDonald-Judd potentials. Conservation laws are derived after generalizing the energy-momentum stress tensor for second order Lagrangians and, in the case of a test-particle population in a given plasma background, the action principle is shown to correspond to the Langevin equation for individual particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Precise laser gyroscope for autonomous inertial navigation

Requirements to gyroscopes of strapdown inertial navigation systems for aircraft application are formulated. The construction of a ring helium – neon laser designed for autonomous navigation is described. The processes that determine the laser service life and the relation between the random error of the angular velocity measurement and the surface relief features of the cavity mirrors are analysed. The results of modelling one of the promising approaches to processing the laser gyroscope signals are presented. (laser gyroscopes)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Magnetic dipole interactions in crystals

The influence of magnetic dipole interactions (MDIs) on the magnetic properties of local-moment Heisenberg spin systems is investigated. A general formulation is presented for calculating the eigenvalues λ and eigenvectors μ ˆ of the MDI tensor of the magnetic dipoles in a line (one dimension, 1D), within a circle (2D) or a sphere (3D) of radius r surrounding a given moment μ → i for given magnetic propagation vectors k for collinear and coplanar noncollinear magnetic structures on both Bravais and non-Bravais spin lattices. Results are calculated for collinear ordering on 1D chains, 2D square and simple-hexagonal (triangular) Bravais lattices, 2D honeycomb and kagomé non-Bravais lattices, and 3D cubic Bravais lattices. The λ and μ ˆ values are compared with previously reported results. Calculations for collinear ordering on 3D simple tetragonal, body-centered tetragonal, and stacked triangular and honeycomb lattices are presented for c/a ratios from 0.5 to 3 in both graphical and tabular form to facilitate comparison of experimentally determined easy axes of ordering on these Bravais lattices with the predictions for MDIs. Comparisons with the easy axes measured for several illustrative collinear antiferromagnets (AFMs) are given. The calculations are extended to the cycloidal noncollinear 120 ° AFM ordering on the triangular lattice where λ is found to be the same as for collinear AFM ordering with the same k. The angular orientation of the ordered moments in the noncollinear coplanar AFM structure of GdB 4 with a distorted stacked 3D Shastry-Sutherland spin-lattice geometry is calculated and found to be in disagreement with experimental observations, indicating the presence of another source of anisotropy. Similar calculations for the undistorted 2D and stacked 3D Shastry-Sutherland lattices are reported. The thermodynamics of dipolar magnets are calculated using the Weiss molecular field theory for quantum spins, including the magnetic transition temperature T m and the ordered moment, magnetic heat capacity, and anisotropic magnetic susceptibility χ versus temperature T . The anisotropic Weiss temperature θ p in the Curie-Weiss law for T>T m is calculated. A quantitative study of the competition between FM and AFM ordering on cubic Bravais lattices versus the demagnetization factor in the absence of FM domain effects is presented. The contributions of Heisenberg exchange interactions and of the MDIs to T m and to θ p are found to be additive, which simplifies analysis of experimental data. Some properties in the magnetically-ordered state versus T are presented, including the ordered moment and magnetic heat capacity and, for AFMs, the dipolar anisotropy of the free energy and the perpendicular critical field. The anisotropic χ for dipolar AFMs is calculated both above and below the Néel temperature T N and the results are illustrated for a simple tetragonal lattice with c/a>1, c/a=1 (cubic), and c/a<1 , where a change in sign of the χ anisotropy is found at c/a=1 . Finally, following the early work of Keffer [Phys. Rev. 87, 608 (1952)], the dipolar anisotropy of χ above T N =69 K of the prototype collinear Heisenberg-exchange-coupled tetragonal compound MnF 2 is calculated and found to be in excellent agreement with experimental single-crystal literature data above 130 K, where the smoothly increasing deviation of the experimental data from the theory on cooling from 130 K to T N is deduced to arise from dynamic short-range collinear c -axis AFM ordering in this temperature range driven by the exchange interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effects of composition and growth parameters on phase formation in multicomponent aluminum garnet crystals

The effects of com­position on the phase formation of multicom­ponent garnet crystals grown via directional solidification by the micro-pulling-down method are studied. A relatively wide range of rare-earth (RE) average ionic radii (AIR) is explored by formulating ten com­positions from the system (Lu,Y,Ho,Dy,Tb,Gd) 3 Al 5 O 12 . Crystals were grown at either 0.05 or 0.20 mm min -1 . The hypothesis is that multicom­ponent com­pounds with large AIR will form secondary phases as the single-RE aluminum garnets formed by larger Tb 3+ or Gd 3+ ; this will result in crystals of poor optical quality. Crystals with large AIR have a central opaque region in optical microscopy images, which is responsible for their reduced transparency com­pared to crystals with small AIR. Slow pulling rates suppress the formation of the opaque region in crystals with inter­mediate AIR. Powder and single-crystal X-ray diffraction and electron probe microanalysis results indicate that the opaque region is a perovskite phase. Scanning electron microscopy and energy dispersive spectroscopy measurements reveal eutectic inclusions at the outer surface of the crystals. Finally, the concentration of the eutectic inclusions increases with increasing AIR.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A Plug-and-Play Priors Approach for Solving Nonlinear Imaging Inverse Problems

In the past two decades, nonlinear image reconstruction methods have led to substantial improvements in the capabilities of numerous imaging systems. Such methods are traditionally formulated as optimization problems that are solved iteratively by simultaneously enforcing data consistency and incorporating prior models. Recently, the Plug-and-Play Priors (PPP) framework suggested that by using more sophisticated denoisers, not necessarily corresponding to an optimization objective, it is possible to improve the quality of reconstructed images. Here in this letter, we show that the PPP approach is applicable beyond linear inverse problems. In particular, we develop the fast iterative shrinkage/thresholding algorithm variant of PPP for model-based nonlinear inverse scattering. The key advantage of the proposed formulation over the original ADMM-based one is that it does not need to perform an inversion on the forward model. We show that the proposed method produces high quality images using both simulated and experimentally measured data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Optimal Scheduling of an Isolated Microgrid With Battery Storage Considering Load and Renewable Generation Uncertainties

By modeling the uncertainty of spinning reserves provided by energy storage with probabilistic constraints, a new optimal scheduling mode is proposed in this paper for minimizing the operating costs of an isolated microgrid (MG) by using chance-constrained programming. The model is transformed into a readily solvable mixed integer linear programming formulation in general algebraic modeling system (GAMS) via a proposed discretized step transformation approach and finally solved by applying the CPLEX solver. By properly setting the confidence levels of the spinning reserve probability constraints, the MG operation can achieve a tradeoff between reliability and economy. The test results on the modified Oak Ridge National Laboratory Distributed Energy Control and Communication lab MG test system reveal that the proposal significantly exceeds the commonly used hybrid intelligent algorithm with much better and more stable optimization results and significantly reduced calculation times.

Li, Yang↗
DOE OSTI2018

On Lagrangian Construction of Massive Higher Spins Supermultiplets in 3D AdS

We give the component Lagrangian construction of the on-shell N=1 massive higher spin supermultiplets in three dimensional (3D) anti de Sitter (AdS ) space. The approach is based on the frame-like gauge invariant formulation, where massive higher spin fields are realized through a system of massless ones. We construct massive supermultiplets using a formulation of the massive fields in terms of the set of gauge invariant objects (curvatures) in the process of their consistent supersymmetric deformation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗