Calculation of magnetic fields and currents in axisymmetric systems of inductively coupled moving conductors

An algorithm for calculating magnetic fields and currents in axisymmetric systems of inductively coupled moving and stationary conductors is developed using a hybrid method which combines the finite and boundary element methods. The finite element method is used to approximate the unsteady diffusion equation for the θ-component of the magnetic vector potential in the conductors, and the boundary-element method is employed to eliminate the space around the conductors. The proposed method takes into account the connections of the conductors with each other or with external energy sources by means of ideal electrical circuits with lumped parameters R, L, and C. An effective method is developed to take into account the external circuits by an appropriate modification of the mass matrix and the source vector of the obtained system of ordinary differential equations. Examples of using the method to calculate the fields of single- and multi-turn solenoids, magnetic flux concentrators, and induction accelerators with various methods of delivering external electromagnetic energy are considered. The high computational efficiency of the method is shown, in particular, for the case of constant electrothermal properties and sizes of the conductors. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 1, pp. 22–29, January–February, 2009.     

Natural oscillations of a liquid of finite electrical conduction in the presence of an external magnetic field

The effect of the finite electrical conduction (finiteness of the magnetic Reynolds number), which is considered a dissipative factor, on small natural oscillations of an ideal heavy liquid of finite depth whose free surface borders on vacuum is studied. A constant external horizontal magnetic field is applied to the liquid. The energy-balance equation is derived, and the theorem of wave attenuation with time is proved. Numerical calculations and the resulting asymptotic formulas give a complete pattern of the spectrum, including its continuous part. The amplitude-frequency characteristics of the wave modes are presented. Rostov State University, Rostov-on-Don 344007. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 2, pp. 3–10, March–April, 2000.     

Electromechanical Characterization of Au Thin Films using Micro-tensile Testing

This paper presents the test technique about measurement of electrical resistance changes of thin films during tensile testing. In this work, we used a real-time digital image correlation strain measurement system coupled with micro-tensile testing unit and voltage/current sourcemeter. This system has the advantage of real time displacement monitoring with a resolution of 50 nm during the micro-tensile testing, with the ability to measure the variation in electrical resistance of the specimen at the same time. We obtained the complete testing data for the stress–strain curve and associated electrical resistance-strain curve for 1 and 2 μm-thick freestanding gold films. Young’s modulus was about 61~69 GPa and 0.2% offset yield strength was about 361~402 MPa. In case of the electrical resistance, rapid change was observed under the elastic regime, while less obvious under the plastic regime. We also conducted finite element analysis, and this result implied that the electrical resistivity would not be constant during micro-tensile testing.     

The Effect of Pore Structure on the Electrical Conductivity of Ti

Open-pore Ti foam samples with porosity in the range of 10–70% and average pore size of 150–400 μm was fabricated by powder metallurgy method using polymethyl methacrylate (PMMA) as space holder initially. The resulting foam is anisotropic: the pores are spheroidal, being shorter along the pressing direction than in the pressing plane. The pore anisotropy decreases as the size of the polymethyl methacrylate (PMMA) particles used increases and hence with pore size, which leads to a higher conductivity in the plane of the pressing. As the porosity increases, the conductivity of porous Ti decreases dramatically. The porosity e{\varepsilon} dependence of the electrical conductivity σ could be well described by Maxwell approximation, while the differential effective medium approximation is only suitable to porous Ti with finite size of 400 μm in the porosity range of 40–70%, i.e., high porosity metal with randomly oriented spheroids.     

The axisymmetric lamb’s problem for a finite prestrained half-space covered with a finite prestretched layer

The piecewise-homogeneous body model and the three-dimensional linearized theory of elastic waves in prestressed bodies are used to solve the axisymmetric time-harmonic Lamb’s problem for a finite prestrained half-space covered with a finite prestretched layer. It is assumed that the half-space and layer are incompressible and their deformation is described by the Treloar potential. The normal stress at the interface is calculated Published in Prikladnaya Mekhanika, Vol. 43, No. 3, pp. 132–143, March 2007.     

Finite volume element method for analysis of unsteady reaction–diffusion problems

A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element method together. The finite volume method is used to discretize the unsteady reaction-diffusion equation, while the finite element method is applied to estimate the gradient quantities at cell faces. Robustness and efficiency of the combined method have been evaluated on uniform rectangular grids by using available numerical solutions of the two-dimensional reaction-diffusion problems. The numerical solutions demonstrate that the combined method is stable and can provide accurate solution without spurious oscillation along the high-gradient boundary layers.     

Discussion of Antiplane Singularities of Piezoelectric–Dielectric and Piezoelectric–Conductor Wedges

This technical note deals with two special topics from our previous paper (Chue and Chen in Arch Appl Mech 72 673–685, 2003) in Archive of Applied Mechanics: the effects of electrical conditions imposed on the edges and bonded interfaces of piezoelectric–dielectric and piezoelectric–conductor wedges on antiplane problems. After employing relatively realistic electrical conditions, we found that stress and electric displacement singularities are altered when boundary conditions and/or continuity conditions are changed, and we compared the results with those of previous studies.     

Averaged rotations at finite plane strain

The averaged rotations and other mechanical parameters at finite plane strains of an elastic material, which are characterized by a linear relation between the Cauchy stresses and the Almansi strains, are studied. The form of the elastic potential is determined. The displacement problem is reduced to a boundary-value problem for complex potentials, which is solved in terms of Cauchy-type integrals for the specified boundarys displacements. The results obtained are compared with the linear solution. Novosibirsk State University, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 3, pp. 187–196, May–June, 2000.     

Structure of the long-range fields of plane symmetrical subsonic flows

Plane subsonic potential flows near finite and semi-infinite bodies, symmetrical about thex axis directed along the velocity of the incident flow, are considered. The shape of the isolines of the velocity modulus and the angle of velocity vector inclination to the symmetry axis at large distances from the bodies is found. Kazan. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 132–144, May–June, 2000. The work was financially supported by the Russian Foundation for Basic Research (project No. 99-01-00169).     

Fully Coupled THMC Modeling of Wellbore Stability with Thermal and Solute Convection Considered

Wellbore stability analysis is an important topic in petroleum geomechanics. Analytical and numerical analysis of wellbore stability involves the study of interactions among pressure, temperature and chemical changes, and the mechanical response of the rock, a coupled thermal–hydraulic–mechanical–chemical (THMC) process. Thermal and solute convection have usually been overlooked in numerical models. This is appropriate for shales with extremely low permeability, but for shales with intermediate and high permeability (e.g., shale with a disseminated microfissure network), thermal and solute convection should be considered. The challenge of considering advection lies in the numerical oscillation encountered when implementing the traditional Galerkin finite element approach for transient advection–diffusion problems. In this article, we present a fully coupled THMC model to analyze the stress, pressure, temperature, and solute concentration changes around a wellbore. In order to overcome spurious spatial temperature oscillations in the convection-dominated thermal advection–diffusion problem, we place the transient problem into an advection– diffusion-reaction problem framework, which is then efficiently addressed by a stabilized finite element approach, the subgrid scale/gradient subgrid scale method (SGS/GSGS).     

Multiphase Thermohaline Convection in the Earth’s Crust: I. A New Finite Element – Finite Volume Solution Technique Combined With a New Equation of State for NaCl–H2O

We present a new finite element – finite volume (FEFV) method combined with a realistic equation of state for NaCl–H₂O to model fluid convection driven by temperature and salinity gradients. This method can deal with the nonlinear variations in fluid properties, separation of a saline fluid into a high-density, high-salinity brine phase and low-density, low-salinity vapor phase well above the critical point of pure H₂O, and geometrically complex geological structures. Similar to the well-known implicit pressure explicit saturation formulation, this approach decouples the governing equations. We formulate a fluid pressure equation that is solved using an implicit finite element method. We derive the fluid velocities from the updated pressure field and employ them in a higher-order, mass conserving finite volume formulation to solve hyperbolic parts of the conservation laws. The parabolic parts are solved by finite element methods. This FEFV method provides for geometric flexibility and numerical efficiency. The equation of state for NaCl–H₂O is valid from 0 to 750°C, 0 to 4000 bar, and 0–100 wt.% NaCl. This allows the simulation of thermohaline convection in high-temperature and high-pressure environments, such as continental or oceanic hydrothermal systems where phase separation is common.     

Elastic–plastic contact force history and response characteristics of circular plate subjected to impact by a projectile

A new elastic–plastic impact–contact model is proposed in this paper. By adopting the principle of minimum acceleration for elastic–plastic continue at finite deformation, and with the aid of finite difference method, the proposed model is applied in the problem of dynamic response of a clamped thin circular plate subjected to a projectile impact centrally. The impact force history and response characteristics of the target plate is studied in detail. The theoretical predictions of the impact force and plate deflection are in good agreements with those of LDA experimental data. Linear expressions of the maximum impact force/transverse deflection versus impact velocity are given on the basis of the theoretical results. The project supported by the National Natural Science Foundation of China (10532020).     

Asymptotic solution of the problem of the action of a stamp on an elastic layer lying on the surface of a compressible fluid of infinite depth

This paper considers a two-dimensional linear unsteady problem of rigid-stamp indentation on an elastic layer of finite thickness lying on the surface of a compressible fluid of infinite depth. The Lamé equations holds for the elastic layer, and the wave equation for the fluid velocity potential. Using the Laplace and Fourier transforms, the problem is reduced to determining the contact stresses under the stamp from the solution of an integral equation of the first kind, whose kernel has a logarithmic singularity. An asymptotic solution of the problem is constructed for large times of interaction. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 2, pp. 131–142, March–April, 2008.     

Correlation between the high-frequency elastic modulus and the interparticle interaction potential in zirconium oxide colloidal suspensions

In this work, we study the high-frequency elastic modulus of aqueous suspensions made with two kinds of zirconium oxide particles, one commercially available and the other synthesized as monodisperse spheres. The effect of volume fraction of solid, ionic strength (sodium chloride as indifferent electrolyte) and particle geometry is taken into account in the study on this viscoelastic property of the suspensions. Frequency sweeps were performed at a fixed value of the applied shear-stress in order to obtain the frequency-limiting value of the elastic modulus by rheometrical methods. On the other hand, the high-frequency modulus is theoretically calculated independently by means of the models proposed by Buscall and co-workers, Wagner and Bergenholtz and co-workers, which correlate the interaction potential between particles with this rheological parameter. The approach to the interparticle potential is the extended DLVO theory, which considers the electrical repulsion between charged colloidal particles, the van der Waals attraction and the acid–base interaction that can be attractive or repulsive depending on the thermodynamic nature of the solid–liquid interface.     

Non-smooth modelling of electrical systems using the flux approach

The non-smooth modelling of electrical systems, which allows for idealised switching components, is described using the flux approach. The formulations and assumptions used for non-smooth mechanical systems are adopted for electrical systems using the position–flux analogy. For the most important non-smooth electrical elements, like diodes and switches, set-valued branch relations are formulated and related to analogous mechanical elements. With the set-valued branch relations, the dynamics of electrical circuits are described as measure differential inclusions. For the numerical solution, the measure differential inclusions are formulated as a measure complementarity system and discretised with a difference scheme, known in mechanics as time-stepping. For every time-step a linear complementarity problem is obtained. Using the example of the DC–DC buck converter, the formulation of the measure differential inclusions, state reduction and their numerical solution using the time-stepping method is shown for the flux approach.     

Development of a Generalized Version of the Poisson– Nernst–Planck Equations Using the Hybrid Mixture Theory: Presentation of 2D Numerical Examples

A numerical scheme for the transient solution of a generalized version of the Poisson–Nernst–Planck (PNP) equations is presented. The finite element method is used to establish the coupled non-linear matrix system of equations capable of solving the present problem iteratively. The PNP equations represent a set of diffusion equations for charged species, i.e. dissolved ions, present in the pore solution of a rigid porous material in which the surface charge can be assumed neglectable. These equations are coupled to the ‘internally’ induced electrical field and to the velocity field of the fluid. The Nernst–Planck equations describing the diffusion of the ionic species and Gauss’ law in use are, however, coupled in both directions. The governing set of equations is derived from a simplified version of the so-called hybrid mixture theory (HMT). The simplifications used here mainly concerns ignoring the deformation and stresses in the porous material in which the ionic diffusion occurs. The HMT is a special version of the more ‘classical’ continuum mixture theories in the sense that it works with averaged equations at macroscale and that it includes the volume fractions of phases in its structure. The background to the PNP equations can by the HMT approach be described by using the postulates of mass conservation of constituents together with Gauss’ law used together with consistent constitutive laws. The HMT theory includes the constituent forms of the quasistatic version of Maxwell’s equations making it suitable for analyses of the kind addressed in this work. Within the framework of HTM, constitutive equations have been derived using the postulate of entropy inequality together with the technique of identifying properties by Lagrange multipliers. These results will be used in obtaining a closed set of equations for the present problem.     

A Two-Scale Model for Coupled Electro-Chemo-Mechanical Phenomena and Onsager’s Reciprocity Relations in Expansive Clays: II Computational Validation

In Part I Moyne and Murad [Transport in Porous Media 62, (2006), 333–380] a two-scale model of coupled electro-chemo-mechanical phenomena in swelling porous media was derived by a formal asymptotic homogenization analysis. The microscopic portrait of the model consists of a two-phase system composed of an electrolyte solution and colloidal clay particles. The movement of the liquid at the microscale is ruled by the modified Stokes problem; the advection, diffusion and electro-migration of monovalent ions Na⁺ and Cl⁻ are governed by the Nernst–Planck equations and the local electric potential distribution is dictated by the Poisson problem. The microscopic governing equations in the fluid domain are coupled with the elasticity problem for the clay particles through boundary conditions on the solid–fluid interface. The up-scaling procedure led to a macroscopic model based on Onsager’s reciprocity relations coupled with a modified form of Terzaghi’s effective stress principle including an additional swelling stress component. A notable consequence of the two-scale framework are the new closure problems derived for the macroscopic electro-chemo-mechanical parameters. Such local representation bridge the gap between the macroscopic Thermodynamics of Irreversible Processes and microscopic Electro-Hydrodynamics by establishing a direct correlation between the magnitude of the effective properties and the electrical double layer potential, whose local distribution is governed by a microscale Poisson–Boltzmann equation. The purpose of this paper is to validate computationally the two-scale model and to introduce new concepts inherent to the problem considering a particular form of microstructure wherein the clay fabric is composed of parallel particles of face-to-face contact. By discretizing the local Poisson–Boltzmann equation and solving numerically the closure problems, the constitutive behavior of the diffusion coefficients of cations and anions, chemico-osmotic and electro-osmotic conductivities in Darcy’s law, Onsager’s parameters, swelling pressure, electro-chemical compressibility, surface tension, primary/secondary electroviscous effects and the reflection coefficient are computed for a range particle distances and sat concentrations.     

Some Properties of a New Model for Slow Flow of Granular Materials

Some properties of a new continuum model for the bulk flow of a dense granular material in which neighbouring grains are in contact for a finite duration of time and in which the contact force is non-impulsive – the so called slow flow regime – are presented. The model generalises both the plastic potential and double-shearing models and contains an additional kinematic quantity – the intrinsic spin. The stress tensor is, in general, non-symmetric and separate yield conditions govern translational and rotational yield. We consider homogeneous, quasi-static loadings for the symmetric part of the stress and dynamic loading for the anti-symmetric part of the stress. A solution for the stress state in terms of a single parameter, namely the major principal direction of the symmetric part of the stress, is presented. This direction itself is determined by a consideration of the flow equations in the context both dilatant and isochoric simple shear flows. These simple flows are used to complete the characterisation of the relationship between the anti-symmetric part of the stress and the intrinsic spin.     

Predictions of the excess pressure drop of Boger fluids through a 2:1:2 contraction–expansion geometry using non-equilibrium molecular dynamics

Non-equilibrium molecular dynamics are used to generate the flow of polymer solutions, specifically of Boger fluids, through a planar 2:1:2 contraction–expansion geometry. The solvent molecules are represented by Lennard–Jones particles, while linear molecules are described by spring-monomers with a finite extensible non-linear elastic spring potential. The equations for Poiseuille flow are solved using a multiple time-scale algorithm extended to non-equilibrium situations. Simulations are performed at constant temperature using Nose–Hoover dynamics. At simulation conditions, changes in concentration show no significant effect on molecular conformation, velocity profiles, and stress fields, while variations in the Deborah number have a strong influence on fluid response. Increasing the magnitude of the Deborah number (De), larger deformation rates are developed in the flow region. For a Deborah number of one, the non-dimensional pressure drop presents values lower than the correspondent Newtonian case. However, for large Deborah numbers, the pressure drop increases above the Newtonian reference. An effective excess pressure drop above the Newtonian value is predicted for Boger fluids along this geometry.     

Averaging of fuzzy differential equations with delay

For fuzzy differential equations with delay, we substantiate the schemes of complete and partial averaging on a finite interval. Translated from Neliniini Kolyvannya, Vol. 11, No. 3, pp. 316–328, July–September, 2008.