SEARCH · Search NASA

Results for “explicit solution”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Structure Errors in System Identification

An approach to system identification is presented which explicitly takes structure errors into account and hence provides a systematic way for answering questions concerning the magnitude of estimated parameter errors resulting from structural errors. It is indicated that, from this point of view, it is possible to define near equivalence between process and model and to obtain meaningful theoretical results on solution error system identification. It remains to apply these results to large realistic problems such as those involving models of complex man machine systems.

Bekey, G. A.↗

NASA NTRS1988

The next generation

The next generation of computer systems are studied by examining the processes and methodologies. The present generation is ok for small projects, but not so good for large projects. They are not good for addressing the iterative nature of requirements, resolution, and implementation. They do not address complexity issues of requirements stabilization. They do not explicitly address reuse opportunities, and they do not help with people shortages. Therefore, there is a need to define and automate improved software engineering processes. Some help may be gained by reuse and prototyping, which are two sides of the same coin. Reuse library parts are used to generate good approximations to desired solutions, i.e., prototypes. And rapid prototype composition implies use of preexistent parts, i.e., reusable parts.

Yudkin, Howard↗

NASA NTRS1992

Parabolized reduced Navier-Stokes computational techniques

A review is presented of methods in which composite or reduced Navier-Stokes (RNS) equations are treated with a pressure-gradient-based flux-vector splitting. The methods are similar to large Re asymptotic formulations, and streamwise diffusion terms are ignored in favor of an explicit deferred corrector based on higher-order diffusion terms. The methods can be used for 2D and 3D supersonic flows in which the effects of real gas are incorporated. Several examples of the procedure are given, and subsonic and supersonic flows are handled well with relaxation procedures that incorporate multigrid acceleration. Effective solutions are described for problems ranging from incompressible flows and supersonic flows to sharp shocks and reverse-flow capturing.

Rubin, Stanley G.↗

NASA NTRS1993

Asymptotic theory of intermediate- and high-degree solar acoustic oscillations

A second-order asymptotic approximation is developed for adiabatic nonradial p-modes of a spherically symmetric star. The exact solutions of adiabatic oscillations are assumed in the outermost layers, where the asymptotic description becomes invalid, which results in a eigenfrequency equation with model-dependent surface phase shift. For lower degree modes, the phase shift is a function of frequency alone; for high-degree modes, its dependence on the degree is explicitly taken into account.

Brodsky, M.↗

DOE OSTI2022

High level requirements for an IFE development program

This White Paper lays out a framework to help guide the formulation of an IFE technology program that explicitly takes account of the needs of the end-user power plant design, construction, licensing, commissioning and operational phases. This viewpoint seeks to ensure that near-term decisions are consistent with the long-term goals, and seeks to minimize effort spent on developments that are incompatible with an integrated system solution.

Dunne, Mike↗

DOE OSTI2022

The Schwarz alternating method for transient solid dynamics

Abstract In our earlier work, we formulated the Schwarz alternating method as a means for concurrent multiscale coupling in finite deformation solid mechanics for quasi‐static problems. Herein, we advance this method for the study of transient dynamic multiscale solid mechanics problems where information is exchanged back and forth between small and large scales. The extension to dynamics relies on the notion of a global time stepper. Within each global time step, the subdomains are coupled by the standard Schwarz iterative process. Remarkably, each subdomain can use its own time step or even its own time integrator to advance its solution in time, provided that they synchronize at each global time step. We study the performance of the Schwarz method on several examples designed for this purpose. Our numerical experiments demonstrate that the method is capable of coupling regions with different mesh resolutions, different element types, and different time integration schemes (e.g., implicit and explicit), all without introducing any artifacts that afflict other coupling methods for transient dynamics. Finally, we apply the dynamic Schwarz alternating method to the simulation of a bolted joint subjected to dynamic loading, as a demonstration of the performance of the method in a realistic scenario.

Mota, Alejandro↗

DOE OSTI2019

Economic analysis of a nuclear hybrid energy system in a stochastic environment including wind turbines in an electricity grid

The electricity market in the U.S. is moving toward an energy mix including more variable renewable energy (VRE) sources like wind and solar. This increasing VRE penetration is altering the profile of the net electricity demand. Normally, flexible generators (e.g. gas turbines) are used to absorb fluctuations in the electricity demand. Therefore, an increase in VRE also leads to an increase in installed capacity for such flexible generators which, in turn, may increase the price of electricity. Nuclear-Renewable Hybrid Energy Systems (N-R HES) are seen as a solution to this problem. The N-R HES contains a nuclear power plant electrically and/or thermally linked to an industrial process (IP). This paper demonstrates a new methodology and presents the corresponding RAVEN/Modelica-based software framework to evaluate the financial performance of N-R HES. The novelty of the method is that it explicitly incorporates the stochastic nature of N-R HES inputs such as electricity demand and wind speeds. The method is applied to a generic (rather than region-specific) N-R HES. The analysis shows that under the right market conditions, economic benefits can be achieved by adding a suitable IP to the base-load generator. In particular, there exists an optimal relationship between the IP capacity and the demand. Furthermore, global trends indicate VRE penetration raises the cost of electricity due to the added volatility, an issue only mitigated when operating with high profit margins for the IP.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A pseudospectral implicit particle-in-cell method with exact energy and charge conservation

The standard particle-in-cell (PIC) method employs explicit finite-difference (FD) methods (e.g. the leap-frog scheme) for both spatial and temporal integrations. Here, we employ a pseudospectral method for solving the Poisson equation and a fully implicit time integration to achieve exact energy conservation. The advantage of a pseudospectral field solver is its spectral accuracy in solving field solutions. Earlier studies of implicit time integration of PIC FD equations can enforce exact energy exchange between field and particles, resulting in exact energy-conserving schemes. Here, we prove that the exact energy conservation property can be carried over to the pseudospectral scheme. Simultaneously, we provide a solution to ensure a pseudospectral charge continuity equation. We demonstrate the new scheme in a 2D electrostatic PIC code. In conclusion, theoretical results are confirmed via numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1984

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.↗

NASA NTRS1994

Structural optimization with approximate sensitivities

Computational efficiency in structural optimization can be enhanced if the intensive computations associated with the calculation of the sensitivities, that is, gradients of the behavior constraints, are reduced. Approximation to gradients of the behavior constraints that can be generated with small amount of numerical calculations is proposed. Structural optimization with these approximate sensitivities produced correct optimum solution. Approximate gradients performed well for different nonlinear programming methods, such as the sequence of unconstrained minimization technique, method of feasible directions, sequence of quadratic programming, and sequence of linear programming. Structural optimization with approximate gradients can reduce by one third the CPU time that would otherwise be required to solve the problem with explicit closed-form gradients. The proposed gradient approximation shows potential to reduce intensive computation that has been associated with traditional structural optimization.

Patnaik, S. N.↗

NASA NTRS1996

On the Retrieval of Lightning Radio Sources from Time-of-Arrival Data

We examine the problem of retrieving three-dimensional lightning locations from radio frequency Time-Of-Arrival (TOA) measurements. Arbitrary antenna locations are considered. By judiciously differencing measurements that are related to the location of the antennas and their excitation times, the problem is converted from the initial spherical nonlinear form to a system of linear equations. In the linear formalism, the source location and time-of-occurrence is viewed geometrically as an intersection of hyperplanes in the four-dimensional Minkowski space (x,y,z,t). The linear equations are solved to obtain explicit analytic expressions for the location and time variables. Retrieval errors are not interpreted with conventional Geometrical Dilution of Precision (GDOP) arguments as discussed by Holmes and Reedy (1951), but with more recent inversion analyses considered by Twomey (1977). Measurement errors are propagated analytically so that the specific effect of these errors on the solution is clarified. The sensitivity of the solution on the number of antennas used, antenna network geometry, source position, and measurement differencing schemes are discussed in terms of the eigenvalues of the linear system.

Koshak, William J.↗

NASA NTRS2004

Combined In-Plane and Through-the-Thickness Analysis for Failure Prediction of Bolted Composite Joints

Although two-dimensional methods provide accurate predictions of contact stresses and bolt load distribution in bolted composite joints with multiple bolts, they fail to capture the effect of thickness on the strength prediction. Typically, the plies close to the interface of laminates are expected to be the most highly loaded, due to bolt deformation, and they are usually the first to fail. This study presents an analysis method to account for the variation of stresses in the thickness direction by augmenting a two-dimensional analysis with a one-dimensional through the thickness analysis. The two-dimensional in-plane solution method based on the combined complex potential and variational formulation satisfies the equilibrium equations exactly, and satisfies the boundary conditions and constraints by minimizing the total potential. Under general loading conditions, this method addresses multiple bolt configurations without requiring symmetry conditions while accounting for the contact phenomenon and the interaction among the bolts explicitly. The through-the-thickness analysis is based on the model utilizing a beam on an elastic foundation. The bolt, represented as a short beam while accounting for bending and shear deformations, rests on springs, where the spring coefficients represent the resistance of the composite laminate to bolt deformation. The combined in-plane and through-the-thickness analysis produces the bolt/hole displacement in the thickness direction, as well as the stress state in each ply. The initial ply failure predicted by applying the average stress criterion is followed by a simple progressive failure. Application of the model is demonstrated by considering single- and double-lap joints of metal plates bolted to composite laminates.

Kradinov, V.↗

DOE OSTI2019

Numerical Simulation of Critical Factors Controlling Heat Extraction from Geothermal Systems Using a Closed-Loop Heat Exchange Method

Closed-loop heat exchange for geothermal energy production involves injecting working fluid down a well that extends through the geothermal resource over a significant length to absorb heat by conduction through the well pipe. The well then needs to return to the surface for energy recovery and fluid re-injection to complete the cycle. We have carried out mixed convective-conductive fluid-flow modeling using a wellbore flow model for TOUGH2 called T2Well to investigate the critical factors that control closed-loop geothermal energy recovery. T2Well solves a mixed explicit-implicit set of momentum equations for flow in the pipe with full coupling to the implicit three-dimensional integral finite difference equations for Darcy flow in the porous medium. T2Well has the option of modeling conductive heat flow from the porous medium to the pipe by means of a semi-analytical solution which makes the computation very efficient because the porous medium does not have to be discretized. When the fully three-dimensional option is chosen, the porous medium is discretized and heat flow to the pipe is by conduction and convection, depending on reservoir permeability and other factors. Simulations of the closed-loop system for a variety of parameter values have been carried out to elucidate the heat recovery process. To the extent that convection may occur to aid in heat delivery to the pipe, the permeability of the geothermal reservoir, whether natural or stimulated, is an important property in heat extraction. The injection temperature and flow rate of the working fluid strongly control the ultimate energy recovery. Pipe diameter also plays a strong role in heat extraction, but is correlated with flow rate. Similarly, the choice of working fluid plays an important role, with water showing better heat extraction than CO2 for certain flow rates, while the CO2 has higher pressure at the production wellhead which can aid in surface energy recovery. In general, we find complex interactions between the critical factors that will require advanced computational approaches to fully optimize.

Oldenburg, Curtis↗

NASA NTRS1975

Approximate solutions of radiative transfer in dusty nebulae. I - Pure hydrogen nebulae

This is the first of a series of papers dealing with transfer of radiation in dusty nebulae. The formulae for the general problem are presented. This paper discusses nebulae containing dust and hydrogen. Approximate solutions, accurate to within 50 percent, are obtained for spherical nebulae containing purely absorbing dust. Less accurate solutions are presented for the general case of scattering and absorbing dust. The results are presented in a form convenient for calculating the required number of ionizing photons, the fractional ionization of hydrogen, and the fractions of ionizing photons absorbed by dust and gas in terms of size, hydrogen density, dust albedo, and optical depth. Explicit relations between these quantities are obtained for uniform nebulae. Simple integrals are presented for calculating the above quantities for spherically symmetric nebulae with nonuniform gas and dust distribution.

Petrosian, V.↗

NASA NTRS1985

Implicit TVD schemes for hyperbolic conservation laws in curvilinear coordinates

The Harten (1983, 1984) total variation-diminishing (TVD) schemes, constituting a one-parameter explicit and implicit, second-order-accurate family, have the property of not generating spurious oscillations when applied to one-dimensional, nonlinear scalar hyperbolic conservation laws and constant coefficient hyperbolic systems. These methods are presently extended to the multidimensional hyperbolic conservation laws in curvilinear coordinates. Means by which to linearize the implicit operator and solution strategies, in order to improve the computation efficiency of the implicit algorithm, are discussed. Numerical experiments with steady state airfoil calculations indicate that the proposed linearized implicit TVD schemes are accurate and robust.

Yee, H. C.↗

NASA NTRS1988

Technique for measuring the dielectric constant of thin materials

A practical technique for measuring the dielectric constant of vegetation leaves and similarly thin materials is presented. A rectangular section of the leaf is placed in the tranverse plane in a rectangular waveguide and the magnitude and phase of the reflection coefficient are measured over the desired frequency band using a vector network analyzer. By treating the leaf as an infinitesimally thin resistive sheet, an explicit expression for its dielectric constant is obtained in terms of the reflection coefficient. Because of the thin-sheet approximation, however, this approach is valid only at frequencies below 1.5 GHz. To extend the technique to higher frequencies, higher order approximations are derived and their accuracies are compared to the exact dielectric-slab solution. For a material whose thickness is 0.5 mm or less, the proposed technique was found to provide accurate values of its dielectric constant up to frequencies of 12 GHz or higher. The technique was used to measure the 8 to 12 GHz dielectric spectrum for vegetation leaves, teflon and rock samples.

Sarabandi, K.↗

NASA NTRS1993

A numerical study of confined turbulent jets

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Zhu, J.↗

NASA NTRS2015

Estimation of Surface Temperature and Heat Flux by Inverse Heat Transfer Methods Using Internal Temperatures Measured While Radiantly Heating a Carbon/Carbon Specimen up to 1920 F

The ability to solve the heat conduction equation is needed when designing materials to be used on vehicles exposed to extremely high temperatures; e.g. vehicles used for atmospheric entry or hypersonic flight. When using test and flight data, computational methods such as finite difference schemes may be used to solve for both the direct heat conduction problem, i.e., solving between internal temperature measurements, and the inverse heat conduction problem, i.e., using the direct solution to march forward in space to the surface of the material to estimate both surface temperature and heat flux. The completed research first discusses the methods used in developing a computational code to solve both the direct and inverse heat transfer problems using one dimensional, centered, implicit finite volume schemes and one dimensional, centered, explicit space marching techniques. The developed code assumed the boundary conditions to be specified time varying temperatures and also considered temperature dependent thermal properties. The completed research then discusses the results of analyzing temperature data measured while radiantly heating a carbon/carbon specimen up to 1920 F. The temperature was measured using thermocouple (TC) plugs (small carbon/carbon material specimens) with four embedded TC plugs inserted into the larger carbon/carbon specimen. The purpose of analyzing the test data was to estimate the surface heat flux and temperature values from the internal temperature measurements using direct and inverse heat transfer methods, thus aiding in the thermal and structural design and analysis of high temperature vehicles.

Pizzo, Michelle↗