SEARCH · Search NASA

Results for “analytic 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 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Synergistic data analytics for electromechanical oscillation in electric power systems

Accurate real-time estimation of the four electromechanical oscillation properties, i.e., dominant oscillation modes, mode shapes, participation factors, and coherent groups, is of great importance to assess and mitigate potential electromechanical oscillations in interconnected power systems. While eigenvalue analysis can realize such estimation, it requires precise linearized dynamic models and accurate parameters, which are highly difficult to obtain in practice. Data fusion-based modal estimation methods can extract the properties of electromechanical oscillations from measurement data without the power system model and parameters, but most of the time only one or two property assessments can be accomplished each time. To overcome this challenge, this paper presents a synergistic data analytics solution to characterize the dynamic behaviors of electromechanical oscillations from real-time measurement data. Here, the proposed method uses optimized variable projection, and it is capable of estimating all four electromechanical oscillation behavior properties from measured responses. Case studies are performed using the simulated measurement data of a 16-generator 68-bus test system and the field measurements collected by the PMUs deployed in the Yunnan Power Grid. The results demonstrate that the proposed synergistic data analytics solution can achieve satisfactory performance in capturing the properties of electromechanical oscillations from measurement data and exhibit strong robustness against measurement noise when compared with existing measurement-based methods.

42 ENGINEERING↗

DOE OSTI2020

Verification of the DIF3D Software to Support Fast Reactor Analysis

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the polynomial order used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Verification of the DIF3D Software to Support Fast Reactor Analysis

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. As some additional changes have been made to the ARC software since this point all versions between release 3012 and 3266 can be considered verified as version 3253 was used for all updates in this revision. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the methodology used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative noting that DIF3D will properly notify the user when the peaking values it is providing are potentially incorrect.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Verification of the DIF3D Software to Support Fast Reactor Analysis (Rev. 3)

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. As some additional changes have been made to the ARC software since this point all versions between release 3012 and 3266 can be considered verified as version 3253 was used for all updates in this revision. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the methodology used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative noting that DIF3D will properly notify the user when the peaking values it is providing are potentially incorrect.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1995

A 3D High-Order Unstructured Finite-Volume Algorithm for Solving Maxwell's Equations

A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Liu, Yen↗

NASA NTRS1990

Comparison of analytical and physical modeling of planetesimal accumulation

Three cases for which the analytic solutions to the planetary embryo coagulation equation are known, namely those of constant coagulation rate, of rate proportional to the sum of the masses of the two colliding bodies, and of rate proportional to the product of their masses, are presently used to test the mathematical validity of the Wetherill and Stewart (1989) physical model for the evolution of planetary embryos' mass distribution. In all cases, excellent agreement is found between numerical physical modeling results and those of the analytic solutions. The treatment of the runaway case can proceed via extension of the efforts of Trubnikov (1971).

Wetherill, George W.↗

NASA NTRS1995

Oxidation kinetics of a continuous carbon phase in a nonreactive matrix

Analytical solutions of and experimental results on the oxidation kinetics of carbon in a pore are presented. Reaction rate, reaction sequence, oxidant partial pressure, total system pressure, pore/crack dimensions, and temperature are analyzed with respect to the influence of each on an overall linear-parabolic rate relationship. Direct measurement of carbon recession is performed using two microcomposite model systems oxidized in the temperature range of 700 to 1200 C, and for times to 35 h. Experimental results are evaluated using the derived analytical solutions. Implications on the oxidation resistance of continuous-fiber-reinforced ceramic-matrix composites containing a carbon constituent are discussed.

Eckel, Andrew J.↗

DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗

DOE OSTI2023

Electrothermally Driven Reconfiguration of Microrobotic Beam Structures for the ChipSail System

Solar sailing enables efficient propellant-free attitude adjustment and orbital maneuvers of solar sail spacecraft with high area-to-mass ratios. However, the heavy supporting mass for large solar sails inevitably leads to low area-to-mass ratios. Inspired by chip-scale satellites, a chip-scale solar sail system named ChipSail, consisting of microrobotic solar sails and a chip-scale satellite, was proposed in this work. The structural design and reconfigurable mechanisms of an electrothermally driven microrobotic solar sail made of Al\Ni 50 Ti 50 bilayer beams were introduced, and the theoretical model of its electro-thermo-mechanical behaviors was established. The analytical solutions to the out-of-plane deformation of the solar sail structure appeared to be in good agreement with the finite element analysis (FEA) results. A representative prototype of such solar sail structures was fabricated on silicon wafers using surface and bulk microfabrication, followed by an in-situ experiment of its reconfigurable property under controlled electrothermal actuation. The experimental results demonstrated significant electro-thermo-mechanical deformation of such microrobotic bilayer solar sails, showing great potential in the development of the ChipSail system. Analytical solutions to the electro-thermo-mechanical model, as well as the fabrication process and characterization techniques, provided a rapid performance evaluation and optimization of such microrobotic bilayer solar sails for the ChipSail.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

Earthquake-induced impact of base-isolated buildings: theory, numerical modeling, and design solutions

Earthquake shaking more intense than that used to size the horizontal clearance between a base–isolated building and near–rigid perimeter moat wall will result in hard impact, producing high–frequency, high–amplitude acceleration response in the structure and supported equipment. Herein this paper provides a design solution for the damaging effects of hard impact by installing a compliant engineered element in the load path between the base–isolated building and the moat wall, resulting in soft impact and a much smaller acceleration response. The engineered element assumed herein is a commercial–off–the–shelf marine fender with mechanical properties determined by physical testing. The attachment of a flexible engineered element, with well–defined stiffness and damping, to a near–rigid moat wall, simplifies the numerical modeling of the building–moat wall system and eliminates the need to bound the lateral stiffness of the wall for impact calculations. The simple model of the engineered element can be implemented in commercial finite element codes. Theory is developed for two–sided impact of a single–degree–of–freedom oscillator. Analytical solutions are derived for the shifted first–mode frequency of the impacted oscillator and for its free–vibration response. The shifted first–mode frequency is a function of the composite lateral stiffness of the isolator–engineered element assembly and its earthquake–induced displacement. Local peaks in the spectral response of the impacted oscillator form at odd integer multiples of the shifted first–mode frequency. The analytical solutions can be used to verify, in part, the numerical model used for impact analysis.

42 ENGINEERING↗

DOE OSTI2017

Modeling electrokinetic flows by consistent implicit incompressible smoothed particle hydrodynamics

We present an efficient implicit incompressible smoothed particle hydrodynamics (I2SPH) discretization of Navier-Stokes, Poisson-Boltzmann, and advection-diffusion equations subject to Dirichlet or Robin boundary conditions. It is applied to model various two and three dimensional electrokinetic flows in simple or complex geometries. The I2SPH's accuracy and convergence are examined via comparison with analytical solutions, grid-based numerical solutions, or empirical models. The new method provides a framework to explore broader applications of SPH in microfluidics and complex fluids with charged objects, such as colloids and biomolecules, in arbitrary complex geometries.

Pan, Wenxiao↗

NASA NTRS1980

The effect of convection on the propagation of relativistic galactic electrons

The paper presents exact analytical solutions which describe steady-state transport of relativistic electrons subject to diffusion, convection, and radiation losses. The diffusion coefficient is spatially and energy dependent; the bulk convection velocity is considered to be constant and directed outward from the galactic plane. The electron spectral behavior with respect to galactic height and energy indicates that convection provides a 'break' in the steady-state energy spectrum without recourse to distribution of sources. Radio observations of NGC 891 and NGC 4631 galaxies can be used to estimate their diffusion coefficient and its energy dependence. Application of analytical solutions to our Galaxy indicates that the diffusion coefficient and outflow velocity deduced from relativistic electron measurements and radio observations of their synchrotron radiation agree with values based on the transport of the cosmic ray nucleonic component.

Lerche, I.↗

DOE OSTI2020

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗

DOE OSTI2018

Effective anisotropy due to the surface of magnetic nanoparticles

Analytical solution has been discovered for the second-order effective anisotropy of magnetic nanoparticles of a cubic shape due to the surface anisotropy (SA) of the Néel type. Similarly to the spherical particles, for the simple cubic lattice the grand-diagonal directions (±1,±1,±1) are favored by the effective cubic anisotropy but the effect is twice as strong. Uniaxial core anisotropy and applied magnetic field cause screening of perturbations from the surface at the distance of the domain-wall width and reduce the effect of SA near the energy minima. Yet, screening disappears near the uniaxial energy barrier, and the uniform barrier state of larger particles may become unstable. For these effects the analytical solution is obtained as well, and the limits of the additive formula with the uniaxial and effective cubic anisotropies for the particle are established. Thermally-activated magnetization-switching rates have been computed by the pulse-noise technique for the stochastic Landau-Lifshitz equation for a system of spins.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2011

The two modes extension to the Berk-Breizman equation: Delayed differential equations and asymptotic solutions

The integro-differential Berk-Breizman (BB) equation, describing the evolution of particle-driven wave mode is transformed into a simple delayed differential equation form {nu}{partial_derivative}a({tau})/{partial_derivative}{tau}=a({tau}) -a{sup 2}({tau}- 1) a({tau}- 2). This transformation is also applied to the two modes extension of the BB theory. The obtained solutions are presented together with the derived asymptotic analytical solutions and the numerical results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2002

Small Particle Response to Fluid Motion using Tethered Particles to Simulate Microgravity

This paper reports on ground based work conducted to support the Spaceflight Definition project SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus). SHIVA will advance our understanding of the movement of a particle in a fluid. Gravity usually dominates the equations of motion, but in microgravity as well as on earth other terms can become important. Before two members of our team found an analytical solution of the equations, numerical methods and/or neglecting terms were required. The general solution predicts that the usually neglected history term becomes important when the characteristic viscous time is in the same order as the vibration period and peaks when the two times are equal. In this case three force terms, the Stokes drag, the added mass, and the history drag must all be included in predicting particle movement. We also developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the "virtual apparatus". Using state-of-the-art methods in holography we will quantify the three-dimensional motion of sets of particles, allowing us to test and apply the new analytical solutions. The motion of particles up to 4 mm in diameter in a fluid that oscillates at frequencies up to 100 Hz with amplitudes up to 200 microns is being examined. Ground studies to support the flight development program have employed various schemes to simulate microgravity. One of the most reliable and meaningful methods uses spheres tethered to a fine hair suspended in the fluid. We have also investigated particles with nearly neutral buoyancy. Recordings are made at the peak amplitudes of vibration of the cell providing a measure of the ratio of fluid to particle amplitude. The experiment requires precise location of the particle at the time of recording. The hologram of the particle provides microscopic images of the particle that are used for finding the position with an accuracy of a few microns. To make the experiment more versatile, the spaceflight system will record holograms both on film and electronically. The electronic holograms can be downlinked providing real time data. Results of the ground experiments, the flight experiment design, and data analysis procedures are reported.

Trolinger, J. D.↗

NASA NTRS2002

Small Particle Response to Fluid Motion using Tethered Particles to Simulate Microgravity

Trolinger, J. D.↗

NASA NTRS2002

Finite Element Modeling of the Buckling Response of Sandwich Panels

A comparative study of different modeling approaches for predicting sandwich panel buckling response is described. The study considers sandwich panels with anisotropic face sheets and a very thick core. Results from conventional analytical solutions for sandwich panel overall buckling and face-sheet-wrinkling type modes are compared with solutions obtained using different finite element modeling approaches. Finite element solutions are obtained using layered shell element models, with and without transverse shear flexibility, layered shell/solid element models, with shell elements for the face sheets and solid elements for the core, and sandwich models using a recently developed specialty sandwich element. Convergence characteristics of the shell/solid and sandwich element modeling approaches with respect to in-plane and through-the-thickness discretization, are demonstrated. Results of the study indicate that the specialty sandwich element provides an accurate and effective modeling approach for predicting both overall and localized sandwich panel buckling response. Furthermore, results indicate that anisotropy of the face sheets, along with the ratio of principle elastic moduli, affect the buckling response and these effects may not be represented accurately by analytical solutions. Modeling recommendations are also provided.

Rose, Cheryl A.↗