Development of a three-dimensional nonlinear viscoelastic constitutive model of solid propellant

A three dimensional nonlinear viscoelastic constitutive model for the solid propellant is developed. In their earlier work, the authors have developed an isotropic constitutive model and verified it for one dimensional case. In the present work, the validity of the model is extended to three-dimensional cases. Large deformation, dewetting and cyclic loading effects are treated as the main sources of nonlinear behavior of the solid propellant. Viscoelastic dewetting criteria is used and the softening of the solid propellant due to dewetting is treated by the modulus decrease. The nonlinearities during cyclic loading are accounted for by the functions of the octahedral shear strain measure. The constitutive equation is implemented into a finite element code for the analysis of propellant grains. A commercial finite element package ‘ABAQUS’ is used for the analysis and the model is introduced into the code through a user subroutine. The model is evaluated with different loading conditions and the predicted values are in good agreement with the measured ones. The resulting model applied to analyze a solid propellant grain for the thermal cycling load.

Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.

Spanning avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics: Field dependence and geometrical properties

Spanning avalanches in the 3D Gaussian Random Field Ising Model (3D-GRFIM) with metastable dynamics at T=0 have been studied. Statistical analysis of the field values for which avalanches occur has enabled a Finite-Size Scaling (FSS) study of the avalanche density to be performed. Furthermore, a direct measurement of the geometrical properties of the avalanches has confirmed an earlier hypothesis that several types of spanning avalanches with two different fractal dimensions coexist at the critical point. We finally compare the phase diagram of the 3D-GRFIM with metastable dynamics with the same model in equilibrium at T=0.

Diluted three-dimensional random field Ising model at zero temperature with metastable dynamics.

The influence of vacancy concentration on the behavior of the three-dimensional random field Ising model with metastable dynamics is studied. We have focused our analysis on the number of spanning avalanches which allows us a clear determination of the critical line where the hysteresis loops change from continuous to discontinuous. By a detailed finite-size scaling analysis we determine the phase diagram and numerically estimate the critical exponents along the whole critical line. Finally, we discuss the origin of the curvature of the critical line at high vacancy concentration.

Three dimensional picture of the local structure of 1,4-Polybutadiene from a complete atomistic model and neutron scattering data

An efficient method of combining neutron diffraction data over an extended Q range with detailed atomistic models is presented. A quantitative and qualitative mapping of the organization of the chain conformation in both glass and liquid phase has been performed. The proposed structural refinement method is based on the exploitation of the intrachain features of the diffraction pattern by the use of internal coordinates for bond lengths, valence angles and torsion rotations. Models are built stochastically by assignment of these internal coordinates from probability distributions with limited variable parameters. Variation of these parameters is used in the construction of models that minimize the differences between the observed and calculated structure factors. A series of neutron scattering data of 1,4-polybutadiene at the region 20320 K is presented. Analysis of the experimental data yield bond lengths for C-C and C=C of 1.54 and 1.35 Å respectively. Valence angles of the backbone were found to be at 112 and 122.8 for the CCC and CC=C respectively. Three torsion angles corresponding to the double bond and the adjacent R and β bonds were found to occupy cis and trans, s(, trans and g( and trans states, respectively. We compare our results with theoretical predictions, computer simulations, RIS models, and previously reported experimental results.

Simulation of a domestic ground source heat pump system using a three-dimensional numerical borehole heat exchanger model

Common approaches to the simulation of borehole heat exchangers (BHEs) assume heat transfer in circulating fluid and grout to be in a quasi-steady state and ignore fluctuations in fluid temperature due to transport of the fluid around the loop. However, in domestic ground source heat pump (GSHP) systems, the heat pump and circulating pumps switch on and off during a given hour; therefore, the effect of the thermal mass of the circulating fluid and the dynamics of fluid transport through the loop has important implications for system design. This may also be important in commercial systems that are used intermittently. This article presents transient simulation of a domestic GSHP system with a single BHE using a dynamic three-dimensional (3D) numerical BHE model. The results show that delayed response associated with the transit of fluid along the pipe loop is of some significance in moderating swings in temperature during heat pump operation. In addition, when 3D effects are considered, a lower heat transfer rate is predicted during steady operations. These effects could be important when considering heat exchanger design and system control. The results will be used to develop refined two-dimensional models.

A model of the three-dimensional evolution of Arctic melt ponds on first-year and multiyear sea ice

During winter the ocean surface in polar regions freezes over to form sea ice. In the summer the upper layers of sea ice and snow melts producing meltwater that accumulates in Arctic melt ponds on the surface of sea ice. An accurate estimate of the fraction of the sea ice surface covered in melt ponds is essential for a realistic estimate of the albedo for global climate models. We present a melt-pond–sea-ice model that simulates the three-dimensional evolution of melt ponds on an Arctic sea ice surface. The advancements of this model compared to previous models are the inclusion of snow topography; meltwater transport rates are calculated from hydraulic gradients and ice permeability; and the incorporation of a detailed one-dimensional, thermodynamic radiative balance. Results of model runs simulating first-year and multiyear sea ice are presented. Model results show good agreement with observations, with duration of pond coverage, pond area, and ice ablation comparing well for both the first-year ice and multiyear ice cases. We investigate the sensitivity of the melt pond cover to changes in ice topography, snow topography, and vertical ice permeability. Snow was found to have an important impact mainly at the start of the melt season, whereas initial ice topography strongly controlled pond size and pond fraction throughout the melt season. A reduction in ice permeability allowed surface flooding of relatively flat, first-year ice but had little impact on the pond coverage of rougher, multiyear ice. We discuss our results, including model shortcomings and areas of experimental uncertainty.

Numerical Solution of the Upper-Convected Maxwell Model for Three-Dimensional Free Surface Flows

This work presents a finite difference technique for simulating three-dimensional free surface flows governed by the Upper-Convected Maxwell (UCM) constitutive equation. A Marker-and-Cell approach is employed to represent the fluid free surface and formulations for calculating the non-Newtonian stress tensor on solid boundaries are developed. The complete free surface stress conditions are employed. The momentum equation is solved by an implicit technique while the UCM constitutive equation is integrated by the explicit Euler method. The resulting equations are solved by the finite difference method on a 3D-staggered grid. By using an exact solution for fully developed flow inside a pipe, validation and convergence results are provided. Numerical results include the simulation of the transient extrudate swell and the comparison between jet buckling of UCM and Newtonian fluids.

A mixture theory based method for three-dimensional modeling of reinforce concrete members with embedded crack finite elements

The paper presents a methodology to model three-dimensional reinforced concrete members by means of embedded discontinuity elements based on the Continuum Strong Discontinuous Approach (CSDA). Mixture theory concepts are used to model reinforced concrete as a 31) composite material constituted of concrete with long fibers (rebars) bundles oriented in different directions embedded in it. The effects of the rebars are modeled by phenomenological constitutive models devised to reproduce the axial non-linear behavior, as well as the bond-slip and dowel action. The paper presents the constitutive models assumed for the components and the compatibility conditions chosen to constitute the composite. Numerical analyses of existing experimental reinforced concrete members are presented, illustrating the applicability of the proposed methodology.

Effect of Metal-Ceramic or All-Ceramic Superstructure Materials on Stress Distribution in a Single Implant-Supported Prosthesis: Three-Dimensional Finite Element Analysis

Purpose: This three-dimensional finite element analysis study evaluated the effect of different material combinations on stress distribution within metal-ceramic and all-ceramic single implant-supported prostheses. Materials and Methods: Three-dimensional finite element models reproducing a segment of the maxilla with a missing left first premolar were created. Five groups were established to represent different superstructure materials: GP, porcelain fused to gold alloy; GR, modified composite resin fused to gold alloy; TP, porcelain fused to titanium; TR, modified composite resin fused to titanium; and ZP, porcelain fused to zirconia. A 100-N vertical force was applied to the contact points of the crowns. All models were fixed in the superior region of bone tissue and in the mesial and distal faces of the maxilla section. Stress maps were generated by processing with finite element software. Results: Stress distribution and stress values of supporting bone were similar for the GP, GR, TP, and ZP models (1,574.3 MPa, 1,574.3 MPa, 1,574.3 MPa, and 1,574.2 MPa, respectively) and different for the TR model (1,838.3 MPa). The ZP model transferred less stress to the retention screw (785 MPa) than the other groups (939 MPa for GP, 961 MPa for GR, 1,010 MPa for TP, and 1,037 MPa for TR). Conclusion: The use of different materials to fabricate a superstructure for a single implant-supported prosthesis did not affect the stress distribution in the supporting bone. The retention screw received less stress when a combination of porcelain and zirconia was used. Int J Oral Maxillofac Implants 2011;26:1202-1209

Mean-field model for the interference of matter-waves from a three-dimensional optical trap

Using the mean-field time-dependent Gross-Pitaevskii equation we study the formation of a repulsive Bose-Einstein condensate on a combined optical and harmonic traps in two and three dimensions and subsequent generation of the interference pattern upon the removal of the combined traps as in the experiment by, Greiner et al. [Nature (London 415 (2002) 39]. For optical traps of moderate strength, interference pattern of 27 (9) prominent bright spots is found to be formed in three. (two) dimensions on a cubic (square) lattice in agreement with experiment. Similar interference pattern can also be formed upon removal of the optical lattice trap only. The pattern so formed can oscillate for a long time in the harmonic trap which can be observed experimentally. (C) 2003 Elsevier B.V. B.V. All rights reserved.

An improved three-dimensional model for growth of oxide films induced by laser heating

In this paper we consider a three-dimensional heat diffusion model to explain the growth of oxide films which takes place when a laser beam is shined on and heats a metallic layer deposited on a glass substrate in a normal atmospheric environment. In particular, we apply this model to the experimental results obtained for the dependence of the oxide layer thickness on the laser density power for growth of TiO2 films grown on Ti-covered glass slides. We show that there is a very good agreement between the experimental results and the theoretical predictions from our proposed three-dimensional model, improving the results obtained with the one-dimensional heat diffusion model previously reported. Our theoretical results also show the occurrence of surface cooling between consecutive laser pulses, and that the oxide track surface profile closely follows the spatial laser profile indicating that heat diffusive effects can be neglected in the growth of oxide films by laser heating. © 2001 Elsevier Science B.V.

Three-dimensional analysis of reinforced concrete members via embedded discontinuity finite elements

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modelling the surface circulation and thermal structure of a tropical reservoir using three-dimensional hydrodynamic lake model and remote-sensing data

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)