Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Soft X-ray measurements of transport and MHD activity in the core and edge NSTX plasma (Final Technical Report)

A comprehensive suite of innovative diagnostic instrumentation and techniques has been developed by the Johns Hopkins collaboration for measurements of transport and MHD activity in the core and edge of the National Spherical Torus Experiment (NSTX) at the Princeton Plasma Physics Laboratory (PPPL). These include ultrafast multi-energy poloidal soft X-ray arrays, high spatial resolution multi-energy tangential soft X-ray arrays, high spectral resolution transmission grating imaging VUV spectrometers and ultrafast bolometric arrays, which measure over a broad range of spatial, spectral, and temporal scales the emission from ohmic, RF, and beam driven NSTX plasmas. The instrumentation was designed to fulfill the specific measurement needs of the Spherical Torus (ST) magnetic fusion concept and is accompanied by a comprehensive set of numerical tools for data analysis and modeling, aimed at determining as many as possible plasma parameters and at extracting as much as possible physics information from NSTX experiments. The Johns Hopkins diagnostics have been and continue to be instrumental for essentially all the NSTX experiments. In addition, our own research and experiments using these diagnostics has led over the years to fundamental conclusions for the ST fusion physics, such as the establishment of near neoclassical impurity transport in the NSTX core and edge plasma, and the critical role of the Global Alfven Eigenmode (GAE) MHD activity for the ST core electron thermal transport. The Johns Hopkins Diagnostics continue to be operational on the upgraded NSTX-U and will provide key plasma information in its upcoming high field, current and power experimental campaigns. The present report summarizes the main technical characteristics and measurement capabilities of the diagnostic instrumentation and methods developed within our research and describes a key physics result achieved during this research. Success in our research was made possible also by the excellent support and cooperation of the NSTX research team and of PPPL researchers, who are warmly acknowledged here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Continuous integration data-driven platform of industrial-scale subsurface storage for real-time analytics

This project helped address the growing need for efficient and scalable models to support geological carbon and energy storage, which are crucial for achieving net-zero emissions. Traditionally accurate high-fidelity numerical models have been used to simulate relevant storage processes under a handful of processes, however such models are computationally demanding, making uncertainty quantification impractical. Consequently, we first developed a machine learning framework, based on Graph Neural Operators (GNOs), to improving the accuracy of model predictions for a fixed computational budget. We then developed an Ensemble of Improved Neural Operators (ENO), which uses bagging and Monte Carlo dropout techniques, to further improve prediction accuracy. Lastly, we developed the way to explain progressive transfer learning methods to reduce the amount of training data and computational cost of training (i.e., reduce trainable parameters) when using our models for multiple storage sites. Our numerical investigation, which used real-world case studies, demonstrated that our framework can significantly improve the safety and efficiency of geological storage operations, with potential applications in other domains such as geothermal reservoirs and climate modeling.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Vertical-structure effects on planetary microwave brightness temperature measurements - Applications to the lunar regolith

An attempt to constrain the effects of vertical variations in dielectric properties on lunar microwave observations is presented. A numerical approach for deriving the reflectivity and microwave weighting function of a vertically varying half-space is used, assuming variance in the dielectric properties with depth only, and negligible magnetic effects. The cases of continuous and stratified models of vertical structures are discussed, and a concentration of emitted energy in upper layers is found. The total emitted energy oscillates, varying with the thickness of the upper soil layer, but averaging out interference effects due to random variations in the substrate depth. Consideration is also given to the vertical structure effects on the lunation-mean disk-center brightness temperature, its variations, and the regolith electrical loss, and predicted reflectivity effects by feasible models of the lunar regolith dielectric profile.

Keihm, S. J.↗
NASA NTRS1983

A spectral numerical method for the Navier-Stokes equations with applications to Taylor-Couette flow

A new spectral method for solving the incompressible Navier-Stokes equations in a plane channel and between concentric cylinders is presented. The method uses spectral expansions which inherently satisfy the boundary conditions and the continuity equation and yield banded matrices which are efficiently solved at each time step. In addition, the number of dependent variables is reduced, resulting in a reduction in computer memory requirements. Several test problems have been computed for the channel flow and for flow between concentric cylinders, including Taylor-Couette flow with axisymmetric Taylor vortices and wavy vortices. In all cases, agreement with available experimental and theoretical results is very good.

Moser, R. D.↗
NASA NTRS1986

Modelling short channel mosfets for use in VLSI

In an investigation of metal oxide semiconductor field effect transistor (MOFSET) devices, a one-dimensional mathematical model of device dynamics was prepared, from which an accurate and computationally efficient drain current expression could be derived for subsequent parameter extraction. While a critical review revealed weaknesses in existing 1-D models (Pao-Sah, Pierret-Shields, Brews, and Van de Wiele), this new model in contrast was found to allow all the charge distributions to be continuous, to retain the inversion layer structure, and to include the contribution of current from the pinched-off part of the device. The model allows the source and drain to operate in different regimes. Numerical algorithms used for the evaluation of surface potentials in the various models are presented.

Klafter, Alex↗
NASA NTRS1988

Unsteady solution of incompressible Navier-Stokes equations

The numerical scheme used by the present time-accurate FEM numerical method for incompressible Navier-Stokes equations, using primitive variables as the unknowns, is a Crank-Nicholson implicit treatment of all equation terms with central differencing for space derivatives. The introduction of a continuous auxilliary system in pseudo-time, with artificial compressibility, yields the incompressible solution at the advanced time level; time-accurate solutions are thereby obtained for two-dimensional fluid flows in a square cavity, in the cases of both an impulsively starting lid and an oscillating lid.

Soh, W. Y.↗
NASA NTRS1992

Diffraction by a symmetric material junction simulated with generalized sheet transistion conditions

When diffraction by a material discontinuity in a thick dielectric/ferrite slab is treated for the dual integral equation approach (a variation of the Weiner-Hopf method), the resulting solution displays unknown constants, characteristic of the incompleteness associated with the generalized transition conditions employed in modeling a thick slab. The constants can be shown to depend on the geometry and properties of the discontinuity. Their explicit determination by enforcing field continuity across the slab junction is treated. The field internal to the slab is determined by first expressing the field as an expansion of discrete and continuous eigenmodes valid in the exterior and interior slab region. The expansion or mode coefficients are then identified by recasting the Weiner-Hopf solution in a suitable form, enabling the solution constants to be numerically computed by matching the field expansions to the left and right of the slab discontinuity at a discrete set of points.

Ricoy, Mark A.↗
NASA NTRS1994

Numerical aperture limits on efficient ball lens coupling of laser diodes to single-mode fibers with defocus to balance spherical aberration

The potential capabilities and limitations of single ball lenses for coupling laser diode radiation to single-mode optical fibers have been analyzed; parameters important to optical communications were specifically considered. These parameters included coupling efficiency, effective numerical apertures, lens radius, lens refractive index, wavelength, magnification in imaging the laser diode on the fiber, and defocus to counterbalance spherical aberration of the lens. Limiting numerical apertures in object and image space were determined under the constraint that the lens perform to the Rayleigh criterion of 0.25-wavelength (Strehl ratio = 0.80). The spherical aberration-defocus balance to provide an optical path difference of 0.25 wavelength units was shown to define a constant coupling efficiency (i.e., 0.56). The relative numerical aperture capabilities of the ball lens were determined for a set of wavelengths and associated fiber-core diameters of particular interest for single-mode fiber-optic communication. The results support general continuing efforts in the optical fiber communications industry to improve coupling links within such systems with emphasis on manufacturing simplicity, system packaging flexibility, relaxation of assembly alignment tolerances, cost reduction of opto-electronic components and long term reliability and stability.

Wilson, R. Gale↗
NASA NTRS1995

SEDAC information gateway plan V(1)

This annual update of the Information Gateway Plan incorporates changes recommended by the Socioeconomic Data and Applications Center (SEDAC) User Working Group (UWG) and reflects comments and suggestions from users, collaborators, and the Contracting Officer Technical Representative (COTR). The Information Gateway Plan is a concise and specific plan that outlines SEDAC activities and services in support of the earth and social sciences and other user communities. The SEDAC Information Gateway effort is a primary means by which the Earth Observing System Data and Information System (EOSDIS) can link meaningfully with a broad range of social science data sources and users in ways that lead to tangible benefits to the American people. The SEDAC Information Gateway provides interdisciplinary access to socioeconomic and physical science data and information resources held by SEDAC and numerous other institutions and networks around the world. The Plan describes the areas of research of earth scientists and socioeconomic scientists where interchange of data and information is most needed. It sets guidelines for the continued development of SEDAC's directory of social science datasets and establishes priorities for efforts to make data held by SEDAC or accessible through SEDAC available to the user community. The Plan also describes the means by which the SEDAC user community can access information products specified by the SEDAC Data and Applications Development Plan (DADP). Among other major activities, SEDAC will continue to enhance and operate a directory capability, interoperable with the Global Change Master Directory, that provides the socioeconomic community with information about earth science products and the earth science research community with information about socioeconomic data. The Information Gateway also serves as a unique and powerful access pathway for a wide range of users and potential users of socioeconomic and earth science data, including especially remote sensing data.

Chen, Robert S.↗
NASA NTRS1996

Valuable Data Provided by PASP Plus Flight Experiment After 1 Year in Orbit

Successfully launched on August 3, 1994, the Photovoltaic Array Space Power Plus Diagnostics (PASP Plus) flight experiment is a joint Air Force Phillips Laboratory/NASA Lewis Research Center program designed to test a variety of new and existing photovoltaic (PV) cell and array technologies within the space environment. The experiment consists of 12 different experimental photovoltaic modules, along with numerous diagnostic instruments that measure the space environment and interactions between the experimental modules and that environment. Power Technology Division personnel at Lewis, who had the primary responsibility for integrating the individual photovoltaic module experiments, continue to analyze and interpret the current-voltage characteristics and environmental effects data received back from the spacecraft. The major goals of the experiment included determining the long-term performance of the photovoltaic modules within a high-radiation environment (PASP Plus was launched into an elliptical, high-radiation orbit) and measuring the interactions between the experimental modules and the space plasma when the modules were biased at high positive and negative voltages. Understanding array interactions within the space plasma is critical to spacecraft power system operations and important to providing higher bus voltage designs in the future.

Source record↗
NASA NTRS2011

CFD Study of the Hydrocarbon Boost Low-Pressure Inducer and Kicker in the Presence of a Circumferential Groove

Results are presented of a computational fluid dynamics (CFD) study done in support of Marshall Space Flight Center's (MSFC) sub-scale water flow experiments of the Hydrocarbon Boost (HCB) Oxidizer Turbopump (OTP) being developed by the Air Force Research Laboratory (AFRL) and Aerojet. A circumferential groove may be added to the pump to reduce synchronous cavitation and subsequent bearing loads at a minimal performance cost. However, the energy may reappear as high order cavitation (HOC) that spans a relatively large frequency range. Thus, HOC may have implications for the full-scale OTP inducer in terms of reduced structural margin at higher mode frequencies. Simulations using the LOCI/Stream CFD program were conducted in order to explore the fluid dynamical impact of the groove on the low-pressure inducer and kicker. It was found that the circumferential groove has minimal head performance impact, but causes back-flowing high-swirl fluid to interact with the nearly-axial incoming fluid just above the inducer blades. The high-shear interface between the fluids is Kelvin-Helmholtz unstable, resulting in trains of low pressure regions or 'pearls' forming near the upstream edge of the groove. When the static pressure in these regions becomes low enough and they get cut by the blade leading edge, HOC is thought to occur. Although further work is required, the numerical models indicate that HOC will occur in the runbox of the AFRL/Aerojet HCB OTP. Comparisons to the ongoing water flow experiments will be discussed, as well as possible designs that may mitigate HOC while continuing to reduce synchronous cavitation. December 2011 MSS/LPS/SPS Joint Subcommittee Meeting ABSTRACT SUBMITTAL FORM

Coker, Robert↗
NASA NTRS1998

Acoustics and Aeroperformance of Nozzles With Screwdriver Shaped and Axisymmetric Plugs

The recent experimental and numerical tests of corrugated nozzles have shown some acoustic and thrust benefits relative to traditional round nozzles. For example, a Bluebell nozzle which was obtained by 3D nozzle design incorporating a corrugated cross section nozzle shape with a sinusoidal lip line nozzle edge, can provide an acoustic benefit up to 4dB with about a 1% thrust augmentation. In references, this effect was explained as being the result of the corrugated design producing more efficient mixing of the exhausted jet with ambient air. Based on this argument, the authors have proposed the application of this concept for a centerbody (plug) which can form several vortices downstream from the centerbody. Several different corrugated designs are proposed and described in detail in this paper. The main design is a Screwdriver shaped centerbody or plug (SCR) which was tested experimentally and numerically. The acoustic tests were conducted in the anechoic chamber of the Central AeroHydrodynamics Institute (TsAGI, Moscow) under Civilian Research and Development Foundation (CRDF) grant. These experiments have shown an essential acoustic benefit of about 10-13% with the application of the co-annular nozzles by comparison with the reference round nozzle with the same mass flow rate. However, the expected acoustic benefits with the application of the 4-petal Screwdriver shaped centerbody were not obtained by comparison with the reference axisymmetric centerbody (CON) having the same length and the same cross section areas at the same distance from the nozzle throat. Moreover, for some angles (Theta = 60 deg and 90 deg) noise increase was observed (about 1-3%). These tests will be continued with the goal of obtaining better acoustic results. In particular, acoustic characteristics are hoped to be improved by moving t lie centerbody into the nozzle and using penetrable walls for the SCR and/or for the main nozzle. Preliminary results for such approach are very promising. Aeroperformance effects were analyzed numerically. The nozzle thrust, calculations were based on a full Navier--Stokes equations solver (NSE), and both full and marching Euler codes: CFL3D, CRAFT, and Krayko-Godunov . Grid preparation and its optimization were conducted using GRIDGEN and our own codes. The general conclusion of this numerical analysis is some thrust loss with the application of SCR design (about 1 - 1.5%). But, again, some constructive features of SCR design give some promising perspectives for its application in aviation and domestic industries.

Gilinsky, M.↗
NASA NTRS1994

The luminosity functions of embedded stellar clusters. 2: Numerical results

We use the theory developed previously to calculate the evolving luminosity functions of very young stellar clusters. The luminosity function of the protostars alone quickly attains a characteristic, sharply peaked shape, which it maintains throughout the time of cluster formation. For the pre-main-sequence stars, the function continually changes but displays a pronounced 'step' near a luminosity of 10 Solar luminosity for a prolonged period. At most times, the vast majority of cluster members are pre-main-sequence stars. However, the total cluster luminosity is dominated first by the protostellar and later by the main-sequence component. At these late times, it is the brightest stars which reach the main sequence first. Finaly, a preliminary application of our model to the rho Ophiuchi embedded cluster indicates that star formation has proceeded in that region for 1 x 10(exp 6).

Fletcher, Andre B.↗
NASA NTRS2011

The 2010 Desert Rats Science Operations Test: Outcomes and Lessons Learned

The Desert RATS 2010 Team tested a variety of science operations management techniques, applying experience gained during the manned Apollo missions and the robotic Mars missions. This test assessed integrated science operations management of human planetary exploration using real-time, tactical science operations to oversee daily crew science activities, and a night shift strategic science operations team to conduct strategic level assessment of science data and daily traverse results. In addition, an attempt was made to collect numerical metric data on the outcome of the science operations to assist test evaluation. The two most important outcomes were 1) the production of significant (almost overwhelming) volume of data produced during daily traverse operations with two rovers, advanced imaging systems and well trained, scientifically proficient crew-members, and 2) the degree to which the tactical team s interaction with the surface crew enhanced science return. This interaction depended on continuous real-time voice and data communications, and the quality of science return from any human planetary exploration mission will be based strongly on the aggregate interaction between a well trained surface crew and a dedicated science operations support team using voice and imaging data from a planet s surface. In addition, the scientific insight developed by both the science operations team and the crews could not be measurable by simple numerical quantities, and its value will be missed by a purely metric-based evaluation of test outcome. In particular, failure to recognize the critical importance of this qualitative type interaction may result in mission architecture choices that will reduce science return.

Eppler, D. B.↗
DOE OSTI2024

Exact enforcement of temporal continuity in sequential physics-informed neural networks

The use of deep learning methods in scientific computing represents a potential paradigm shift in engineering problem solving. One of the most prominent developments is Physics-Informed Neural Networks (PINNs), in which neural networks are trained to satisfy partial differential equations (PDEs). While this method shows promise, the standard version has been shown to struggle in accurately predicting the dynamic behavior of time-dependent problems. To address this challenge, methods have been proposed that decompose the time domain into multiple segments, employing a distinct neural network in each segment and directly incorporating continuity between them in the loss function of the minimization problem. In this work we introduce a method to exactly enforce continuity between successive time segments via a solution ansatz. This hard constrained sequential PINN (HCS-PINN) method is simple to implement and eliminates the need for any loss terms associated with temporal continuity. The method is tested for a number of benchmark problems involving both linear and non-linear PDEs. Examples include various first order time dependent problems in which traditional PINNs struggle, namely advection, Allen–Cahn, and Korteweg–de Vries equations. Furthermore, second and third order time-dependent problems are demonstrated via wave and Jerky dynamics examples, respectively. Notably, the Jerky dynamics problem is chaotic, making the problem especially sensitive to temporal accuracy. Finally, the numerical experiments conducted with the proposed method demonstrated superior convergence and accuracy over both traditional PINNs and the soft-constrained counterparts.

42 ENGINEERING↗
DOE OSTI2021

Transport Upscaling under Flow Heterogeneity and Matrix-Diffusion in Three-Dimensional Discrete Fracture Networks

For this work, we investigate the combined effects of network scale flow variability and retention due to matrix-diffusion on the scaling behavior of transport through fractured media. Two of the principal mechanisms controlling the transport of solutes through fractured low-permeability media are broad distributions of flow velocities and retention times in the solid matrix. We study the relative impact of these two processes under different initial conditions using a set of three-dimensional discrete fracture network simulations. We use these simulations to develop and calibrate an upscaled continuous time random walk (CTRW) approach for advective transport based on an Ornstein-Uhlenbeck model for the particle velocities that accounts for the fracture-matrix coupling using a compound Poisson process. This CTRW model can be conditioned on the initial solute distribution and allows to observe late-time scaling behavior at distances beyond what is feasible to observe using high-fidelity direct numerical simulations. We determine that the initial distribution of particles leads to marked differences in the persistent long-term scale behavior in the solute travel time distributions, even those undergoing retention due to matrix diffusion through implementation and analysis of the model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1983

Calculation of afterbody flows with a composite velocity formulation

A recently developed technique for numerical solution of the Navier-Stokes equations for subsonic, laminar flows is investigated. It is extended here to allow for the computation of transonic and turbulent flows. The basic approach involves a multiplicative composite of the appropriate velocity representations for the inviscid and viscous flow regions. The resulting equations are structured so that far from the surface of the body the momentum equations lead to the Bernoulli equation for the pressure, while the continuity equation reduces to the familiar potential equation. Close to the body surface, the governing equations and solution techniques are characteristic of those describing interacting boundary layers. The velocity components are computed with a coupled strongly implicity procedure. For transonic flows the artificial compressibility method is used to treat supersonic regions. Calculations are made for both laminar and turbulent flows over axisymmetric afterbody configurations. Present results compare favorably with other numerical solutions and/or experimental data.

Swanson, R. C.↗