Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Geometry”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Finite Element Analysis of a Deployable Space Structure

The dynamic characteristics of the Solar Array Flight Experiment (SAFE) structure during deployment and retraction are investigated. The SAFE structure consists of a deployable mast with an attached solar blanket designed with accordion type folds to permit packaging in a small volume. The planar form of the blanket geometry during deployment is maintained by a blanket tension/guidewire system. Structurally, the mast is modeled as an Euler beam column with inplane and out of plane bending and finite torsional stiffness. For out of plane motion, the blanket is modeled as a distributed mass uniformly supported by the three guidewires. For inplane motion the blanket displacements are assumed to vary linearly from the mast base to the mast tip. The mathematical model uses a virtual work formulation, required because the axial loading on the mast is nonconservative, combined with assumed beam modes to derive the differential equations of motion. Consideration of the time dependent boundary conditions results in an infinite set of ODE with time dependent coefficients. Finally, correlation of mast tip accelerations to mast base bending moments for specified modal motions are indicated.

Weeks, G. E.↗
NASA NTRS1991

Computational attributes of the integral form of the equation of transfer

Difficulties can arise in radiative and neutron transport calculations when a highly anisotropic scattering phase function is present. In the presence of anisotropy, currently used numerical solutions are based on the integro-differential form of the linearized Boltzmann transport equation. This paper, departs from classical thought and presents an alternative numerical approach based on application of the integral form of the transport equation. Use of the integral formalism facilitates the following steps: a reduction in dimensionality of the system prior to discretization, the use of symbolic manipulation to augment the computational procedure, and the direct determination of key physical quantities which are derivable through the various Legendre moments of the intensity. The approach is developed in the context of radiative heat transfer in a plane-parallel geometry, and results are presented and compared with existing benchmark solutions. Encouraging results are presented to illustrate the potential of the integral formalism for computation. The integral formalism appears to possess several computational attributes which are well-suited to radiative and neutron transport calculations.

Frankel, J. I.↗
NASA NTRSDate not supplied

Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories

Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.

entry↗
DOE OSTI2018

Unwinding the amplituhedron in binary

We present new, fundamentally combinatorial and topological characterizations of the amplituhedron. Upon projecting external data through the amplituhedron, the resulting configuration of points has a specified (and maximal) generalized “winding number”. Equivalently, the amplituhedron can be fully described in binary: canonical projections of the geometry down to one dimension have a specified (and maximal) number of “sign flips” of the projected data. The locality and unitarity of scattering amplitudes are easily derived as elementary consequences of this binary code. Minimal winding defines a natural “dual” of the amplituhedron. This picture gives us an avatar of the amplituhedron purely in the configuration space of points in vector space (momentum-twistor space in the physics), a new interpretation of the canonical amplituhedron form, and a direct bosonic understanding of the scattering super-amplitude in planar N = 4 SYM as a differential form on the space of physical kinematical data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

High-energy Gamma Rays from the Milky Way: Three-dimensional Spatial Models for the Cosmic-Ray and Radiation Field Densities in the Interstellar Medium

High-energy γ -rays of interstellar origin are produced by the interaction of cosmic-ray (CR) particles with the diffuse gas and radiation fields in the Galaxy. The main features of this emission are well understood and are reproduced by existing CR propagation models employing 2D galactocentric cylindrically symmetrical geometry. However, the high-quality data from instruments like the Fermi Large Area Telescope reveal significant deviations from the model predictions on few to tens of degrees scales, indicating the need to include the details of the Galactic spiral structure and thus requiring 3D spatial modeling. In this paper, the high-energy interstellar emissions from the Galaxy are calculated using the new release of the GALPROP code employing 3D spatial models for the CR source and interstellar radiation field (ISRF) densities. Three models for the spatial distribution of CR sources are used that are differentiated by their relative proportion of input luminosity attributed to the smooth disk or spiral arms. Two ISRF models are developed based on stellar and dust spatial density distributions taken from the literature that reproduce local near- to far-infrared observations. The interstellar emission models that include arms and bulges for the CR source and ISRF densities provide plausible physical interpretations for features found in the residual maps from high-energy γ -ray data analysis. The 3D models for CR and ISRF densities provide a more realistic basis that can be used for the interpretation of the nonthermal interstellar emissions from the Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

Aerodynamic Shape Sensitivity Analysis and Design Optimization of Complex Configurations Using Unstructured Grids

A three-dimensional unstructured grid approach to aerodynamic shape sensitivity analysis and design optimization has been developed and is extended to model geometrically complex configurations. The advantage of unstructured grids (when compared with a structured-grid approach) is their inherent ability to discretize irregularly shaped domains with greater efficiency and less effort. Hence, this approach is ideally suited for geometrically complex configurations of practical interest. In this work the nonlinear Euler equations are solved using an upwind, cell-centered, finite-volume scheme. The discrete, linearized systems which result from this scheme are solved iteratively by a preconditioned conjugate-gradient-like algorithm known as GMRES for the two-dimensional geometry and a Gauss-Seidel algorithm for the three-dimensional; similar procedures are used to solve the accompanying linear aerodynamic sensitivity equations in incremental iterative form. As shown, this particular form of the sensitivity equation makes large-scale gradient-based aerodynamic optimization possible by taking advantage of memory efficient methods to construct exact Jacobian matrix-vector products. Simple parameterization techniques are utilized for demonstrative purposes. Once the surface has been deformed, the unstructured grid is adapted by considering the mesh as a system of interconnected springs. Grid sensitivities are obtained by differentiating the surface parameterization and the grid adaptation algorithms with ADIFOR (which is an advanced automatic-differentiation software tool). To demonstrate the ability of this procedure to analyze and design complex configurations of practical interest, the sensitivity analysis and shape optimization has been performed for a two-dimensional high-lift multielement airfoil and for a three-dimensional Boeing 747-200 aircraft.

Taylor, Arthur C., III↗
DOE OSTI2018

Formation of the Conducting Filament in TaO x -Resistive Switching Devices by Thermal-Gradient-Induced Cation Accumulation

Here, the distribution of tantalum and oxygen ions in electroformed and/or switched TaO x -based resistive switching devices has been assessed by high-angle annular dark-field microscopy, X-ray energy-dispersive spectroscopy, and electron energy-loss spectroscopy. The experiments have been performed in the plan-view geometry on the cross-bar devices producing elemental distribution maps in the direction perpendicular to the electric field. The maps revealed an accumulation of +20% Ta in the inner part of the filament with a 3.5% Ta-depleted ring around it. The diameter of the entire structure was approximately 100 nm. The distribution of oxygen was uniform with changes, if any, below the detection limit of 5%. We interpret the elemental segregation as due to diffusion driven by the temperature gradient, which in turn is induced by the spontaneous current constriction associated with the negative differential resistance-type I–V characteristics of the as-fabricated metal/oxide/metal structures. A finite-element model was used to evaluate the distribution of temperature in the devices and correlated with the elemental maps. In addition, a fine-scale (~5 nm) intensity contrast was observed within the filament and interpreted as due phase separation of the functional oxide in the two-phase composition region. Understanding the temperature-gradient-induced phenomena is central to the engineering of oxide memory cells.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Slow manifold reduction as a systematic tool for revealing the geometry of phase space

Many non-dissipative reduced plasma models can be derived from more fundamental non-dissipative models by restricting to an approximate invariant manifold. I present a general systematic procedure for finding the Hamiltonian formulation of a plasma model that can be derived in this manner. Several illustrative examples are considered in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Tip-based proximity ferroelectric switching and piezoelectric response in wurtzite multilayers

Proximity ferroelectricity is a paradigm for inducing ferroelectricity when a nonferroelectric polar material (such as Al⁢ N), which is unswitchable with an external field below the dielectric breakdown field, becomes a practically switchable ferroelectric in direct contact with a thin switchable ferroelectric layer (such as Al 1−𝑥 ⁢Sc 𝑥 ⁢N). Here, we develop a Landau-Ginzburg-Devonshire approach to study the proximity effect of local piezoelectric response and polarization reversal in wurtzite ferroelectric multilayers under a sharp electrically biased tip. Using finite-element modeling, we analyze the probe-induced nucleation of nanodomains, the features of local polarization hysteresis loops and coercive fields in the Al 1−𝑥 ⁢Sc 𝑥 ⁢N/Al⁢ N bilayers and three-layers. Similar to the wurtzite multilayers sandwiched between two parallel electrodes, the regimes of “proximity switching” (when all layers collectively switch) and the regime of “proximity suppression” (when they collectively do not switch) are the only two possible regimes in the probe-electrode geometry. However, the parameters and asymmetry of the local piezoresponse and polarization hysteresis loops depend significantly on the sequence of the layers with respect to the probe. The physical mechanism of proximity ferroelectricity in the local probe geometry is a depolarizing electric field determined by the polarization of the layers and their relative thickness. The field, whose direction is opposite to the polarization vector in the layer(s) with the larger spontaneous polarization (such as Al⁢ N), renormalizes the double-well ferroelectric potential to lower the steepness of the switching barrier in the “otherwise unswitchable” polar layers. Tip-based control of domains in otherwise nonferroelectric layers using proximity ferroelectricity can provide nanoscale control of domain reversal in memory, actuation, sensing, and optical applications. The ability of the tip-induced proximity switching to differentially switch multilayers, based on the order of the layers, provides a powerful tool for selective domain engineering.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Expansion of epicyclic gear dynamic analysis program

The multiple mesh/single stage dynamics program is a gear tooth analysis program which determines detailed geometry, dynamic loads, stresses, and surface damage factors. The program can analyze a variety of both epicyclic and single mesh systems with spur or helical gear teeth including internal, external, and buttress tooth forms. The modifications refine the options for the flexible carrier and flexible ring gear rim and adds three options: a floating Sun gear option; a natural frequency option; and a finite element compliance formulation for helical gear teeth. The option for a floating Sun incorporates two additional degrees of freedom at the Sun center. The natural frequency option evaluates the frequencies of planetary, star, or differential systems as well as the effect of additional springs at the Sun center and those due to a flexible carrier and/or ring gear rim. The helical tooth pair finite element calculated compliance is obtained from an automated element breakup of the helical teeth and then is used with the basic gear dynamic solution and stress postprocessing routines. The flexible carrier or ring gear rim option for planetary and star spur gear systems allows the output torque per carrier and ring gear rim segment to vary based on the dynamic response of the entire system, while the total output torque remains constant.

Boyd, Linda Smith↗
NASA NTRS1992

Is the Oort A-value a universal growth rate limit for accretion disk shear instabilities?

A weak-field local MHD instability that is of importance to accretion disks is examined. The maximum growth rate of the instability is found to be not only independent of the magnetic field strength but independent of field geometry as well. In particular, all Keplerian disks are unstable in the presence of any weak poloidal field, with the ratio of the maximum growth rate to disk angular velocity given by 3/4. The maximum growth rate of any weak field configuration that is not purely toroidal is given by the local Oort A-value of the disk. The behavior is studied by using a form of the dynamical Hill equations. It is conjectured that the Oort A-value is an upper bound to the growth rate of any instability feeding upon the free energy of differential rotation.

Balbus, Steven A.↗
NASA NTRS1993

Performance evaluation of the JPL TOPSAR system: An across-track interferometric SAR system for topographic mapping

An across-track interferometric (XTI) SAR has the capability of determining the three dimensional target location from the recorded radar echos thus providing high spatial resolution both horizontally and vertically. The output of an XTI SAR will typically include both a radar backscatter map and a digital elevation model (DEM). The accuracy of the derived DEM is determined by factors such as data acquisition geometry, signal-to-noise ratio, knowledge of platform position and attitude as well as the accuracy of the processing system. In the summer of 1992 the authors acquired TOPSAR data over the Ft. Irwin area in California, a desert area with significant relief (height standard deviation 150 m). Very accurate DEM's have been derived for this area by the Topographic Engineering Center (TEC) using digital correlation methods on 1:20,000 scale digitized photographs. Corner reflectors were deployed in the area, and their locations were determined to cm accuracies by the Defense Mapping Agency (DMA) using differential GPS techniques. DEM's generated from the acquired radar data have been rotated and translated to overlay the reference DEM's provided by TEC. A detailed description of the errors and their characteristics will be given. The standard deviations of the radar data sets we acquired measured over a 5.6 by 7 km area were 1.9 and 2.3 m. The corresponding numbers for flat areas were 1.1 and 2.0 m and for mountain areas 3.3 and 2.2 m.

Madsen, Soren N.↗
DOE OSTI2013

A GENERALIZED STOCHASTIC COLLOCATION APPROACH TO CONSTRAINED OPTIMIZATION FOR RANDOM DATA IDENTIFICATION PROBLEMS

We present a scalable, parallel mechanism for stochastic identification/control for problems constrained by partial differential equations with random input data. Several identification objectives will be discussed that either minimize the expectation of a tracking cost functional or minimize the difference of desired statistical quantities in the appropriate $L^p$ norm, and the distributed parameters/control can both deterministic or stochastic. Given an objective we prove the existence of an optimal solution, establish the validity of the Lagrange multiplier rule and obtain a stochastic optimality system of equations. The modeling process may describe the solution in terms of high dimensional spaces, particularly in the case when the input data (coefficients, forcing terms, boundary conditions, geometry, etc) are affected by a large amount of uncertainty. For higher accuracy, the computer simulation must increase the number of random variables (dimensions), and expend more effort approximating the quantity of interest in each individual dimension. Hence, we introduce a novel stochastic parameter identification algorithm that integrates an adjoint-based deterministic algorithm with the sparse grid stochastic collocation FEM approach. This allows for decoupled, moderately high dimensional, parameterized computations of the stochastic optimality system, where at each collocation point, deterministic analysis and techniques can be utilized. The advantage of our approach is that it allows for the optimal identification of statistical moments (mean value, variance, covariance, etc.) or even the whole probability distribution of the input random fields, given the probability distribution of some responses of the system (quantities of physical interest). Our rigorously derived error estimates, for the fully discrete problems, will be described and used to compare the efficiency of the method with several other techniques. Numerical examples illustrate the theoretical results and demonstrate the distinctions between the various stochastic identification objectives.

Webster, Clayton G↗
DOE OSTI2013

Quantum-mechanical predictions of electron-induced ionization cross sections of DNA components

Ionization of biomolecules remains still today rarely investigated on both the experimental and the theoretical sides. In this context, the present work appears as one of the first quantum mechanical approaches providing a multi-differential description of the electron-induced ionization process of the main DNA components for impact energies ranging from the target ionization threshold up to about 10 keV. The cross section calculations are here performed within the 1st Born approximation framework in which the ejected electron is described by a Coulomb wave whereas the incident and the scattered electrons are both described by a plane wave. The biological targets of interest, namely, the DNA nucleobases and the sugar-phosphate backbone, are here described by means of the GAUSSIAN 09 system using the restricted Hartree-Fock method with geometry optimization. The theoretical predictions also obtained have shown a reasonable agreement with the experimental total ionization cross sections while huge discrepancies have been pointed out with existing theoretical models, mainly developed within a semi-classical framework.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1992

Evaluation and application of the Baldwin-Lomax turbulence model in two-dimensional, unsteady, compressible boundary layers with and without separation in engine inlets

There is a practical need to model high speed flows that exist in jet engine inlets. The boundary layers that form in these inlets may be turbulent or laminar and either separated or attached. Also, unsteady supersonic inlets may be subject to frequent changes in operating conditions. Some changes in the operating conditions of the inlets may include varying the inlet geometry, bleeds and bypasses, and rotating or translating the centerbody. In addition, the inlet may be either started or unstarted. Therefore, a CFD code, used to model these inlets, may have to run for several different cases. Also, since the flow conditions through an unsteady inlet may be continually fluctuating, the CFD code which models these flows may have to be run over many time steps. Therefore, it would be beneficial that the code run quickly. Many turbulence models, however, are cumbersome to implement and require a lot of computer time to run, since they add to the number of differential equations to be solved to model a flow. The Baldwin-Lomax turbulence model is a popular model. It is an algebraic, eddy viscosity model. The Baldwin-Lomax model is used in many CFD codes because it is quick and easy to implement. In this paper, we will discuss implementing the Baldwin-Lomax turbulence model for both steady and unsteady compressible flows. In addition, these flows may be either separated or attached. In order to apply this turbulence model to flows which may be subjected to these conditions, certain modifications should be made to the original Baldwin-Lomax model. We will discuss these modifications and determine whether the Baldwin-Lomax model is a viable turbulence model that produces reasonably accurate results for high speed flows that can be found in engine inlets.

Sakowski, Barbara↗
NASA NTRS1993

Two wavelength satellite laser ranging using SPAD

When ranging to satellites with lasers, there are several principal contributions to the error budget: from the laser ranging system on the ground, from the satellite retroarray geometry, and from the atmosphere. Using a single wavelength, we have routinely achieved a ranging precision of 8 millimeters when ranging to the ERS-1 and Starlette satellites. The systematic error of the atmosphere, assuming the existing dispersion models, is expected to be of the order of 1 cm. Multiple wavelengths ranging might contribute to the refinement of the existing models. Taking into account the energy balance, the existing picosecond lasers and the existing receiver and detection technology, several pairs or multiple wavelengths may be considered. To be able to improve the atmospheric models to the subcentimeter accuracy level, the differential time interval (DTI) has to be determined within a few picoseconds depending on the selected wavelength pair. There exist several projects based on picosecond lasers as transmitters and on two types of detection techniques: one is based on photodetectors, like photomultipliers or photodiodes connected to the time interval meters. Another technique is based on the use of a streak camera as an echo signal detector, temporal analyzer, and time interval vernier. The temporal analysis at a single wavelength using the streak camera showed the complexity of the problem.

Prochazka, Ivan↗
DOE OSTI2015

Preparation and characterization of novel foamed porous glass-ceramics

Foamed glass-ceramics without using foaming agent have been synthesized in a novel glass system of SrO-CaO-Al{sub 2}O{sub 3}-TiO{sub 2}-B{sub 2}O{sub 3}-SiO{sub 2}-P{sub 2}O{sub 5}-M{sub x}O{sub y} (where M = Ba, Mg, La, Ce and Ni) by a simple process of powder sintering. The glass and glass-ceramics are characterized by dilatometry, differential scanning calorimetry (DSC), heating stage microscopy (HSM), X-ray diffraction (XRD), field emission scanning electron microscopy (FESEM), optical microscopy and Fourier transformed infrared spectroscopy (FTIR). All the glasses formed are amorphous and the glass transition temperature and dilatometric softening temperature of these glasses are found to be in the range 673–678 °C and 706–728 °C respectively. The glasses are highly stable as indicated by the DSC evaluated glass stability parameters of the range 195–240 °C. Quantitative sintering study of glass powder compacts revealed swelling in the samples with NiO and CeO{sub 2} corresponding to a geometry change of 75 and 108% around 900 °C respectively. With reference to this finding the glass powder compacts are heated to 900 °C and the foamed glass-ceramics are obtained. Characteristic crystalline silicate phases have been identified in the XRD studies and their microstructures are recorded by FESEM. Optical microscope study of the foamed samples revealed formation of bigger foamed cavity with residual pores in samples with NiO and CeO{sub 2} in comparison to samples with BaO, MgO and La{sub 2}O{sub 3}. The mean pore diameters of the samples with NiO and CeO{sub 2} are determined to be 43 and 32 μm, and their respective porosities are 2.34 and 1.82 cm{sup 3}/g respectively. Thus NiO and CeO{sub 2} are found to be very effective to obtain foamed glass-ceramics without using foaming agent by the viscous flow sintering of fine glass powder compacts along with the reduction of the respective polyvalent ions. - Highlights: • Synthesis of foamed porous glass-ceramics without foaming agent by sintering method • Only powder compact yielded foamed porous glass-ceramics but bulk glass did not. • Glasses containing NiO and CeO{sub 2} exhibited significant foaming efficiency. • Bloating of entrapped gas during viscous flow sintering is the origin of foaming. • Residual void created pores in the sintered glass-ceramics as evidenced in FESEM.

36 MATERIALS SCIENCE↗