Mixed-dimensional coupling in finite element models

In many finite element analysis models it would be desirable to combine reduced- or lower-dimensional element types with higher-dimensional element types in a single model. In order to achieve compatibility of displacements and stress equilibrium at the junction or interface between the differing element types, it is important in such cases to integrate into the analysis some scheme for coupling the element types. A novel and effective scheme for establishing compatibility and equilibrium at the dimensional interface is described and its merits and capabilities are demonstrated. Copyright (C) 2000 John Wiley & Sons, Ltd.

The Development of a Two-Dimensional Transient Catalyst Model for Direct Injection Two-Stroke Applications.

This paper describes the development of a two-dimensional transient catalyst model. Although designed primarily for two-stroke direct injection engines, the model is also applicable to four-stroke lean burn and diesel applications. The first section describes the geometries, properties and chemical processes simulated by the model and discusses the limitations and assumptions applied. A review of the modeling techniques adopted by other researchers is also included. The mathematical relationships which are used to represent the system are then described, together with the finite volume method used in the computer program. The need for a two-dimensional approach is explained and the methods used to model effects such as flow and temperature distribution are presented. The problems associated with developing surface reaction rates are discussed in detail and compared with published research. Validation and calibration of the model is achieved by comparing predictions with measurements from a flow reactor. While an extensive validation process, involving detailed measurements of gas composition and thermal gradients, has been completed, the analysis is too detailed for publication here and is the subject of a separate technical paper.

Isometric embedding of finite ultrametric spaces in Banach spaces

We prove that for any finite ultrametric space M and any infinite-dimensional Banach space B there exists an isometric embedding of M into B.

A Virtual Inspection Technique for Assessing the Dimensional Accuracy of Forged Compressor Blades Using FE Modeling and CMM Inspection

This paper presents research for developing a virtual inspection system that evaluates the dimensional tolerance of forged aerofoil blades formed using the finite element (FE) method. Conventional algorithms adopted by modern coordinate measurement processes have been incorporated with the latest free-form surface evaluation techniques to provide a robust framework for the dimensional inspection of FE aerofoil models. The accuracy of the approach had been verified with a strong correlation obtained between the virtual inspection data and coordinate measurement data from corresponding aerofoil components.

Three-Dimensional Analysis of Seepage below and around Hydraulic Structures

This paper investigates the problem of seepage under the floor of hydraulic structures considering the compartment of flow that seeps through the surrounding banks of the canal. A computer program, utilizing a finite-element method and capable of handling three-dimensional (3D) saturated–unsaturated flow problems, was used. Different ratios of canal width/differential head applied on the structure were studied. The results produced from the two-dimensional (2D) analysis were observed to deviate largely from that obtained from 3D analysis of the same problem, despite the fact that the porous medium was isotropic and homogeneous. For example, the exit gradient obtained from 3D analysis was as high as 2.5 times its value obtained from 2D analysis. Uplift force acting upwards on the structure has also increased by about 46% compared with its value obtained from the 2D solution. When the canal width/ differential head ratio was 10 or higher, the 3D results were comparable to the 2D results. It is recommended to construct a core of low permeability soil in the banks of canal to reduce the seepage losses, uplift force, and exit gradient.

Two-dimensional hygrothermal transfer in porous building materials

A two-dimensional mathematical model for evaluating the simultaneous heat and moisture migration in porous building materials was proposed. Vapor content and temperature were chosen as the principal driving potentials. The numerical solution was based on the control volume finite difference technique with fully implicit scheme in time. Two validation experiments were developed in this study. The evolution of transient moisture distributions in both one-dimensional and two-dimensional cases was measured. A comparison between experimental results and those obtained by the numerical model proves that they are fully consistent with each other. The model can be easily integrated into a whole building heat, air and moisture transfer model. Another main advantage of the present numerical method lies in the fact that the required moisture transport properties are comparatively simple and easy to determine.

Modelling the mechanical response of the Reed-mouthpiece-lip system of a clarinet. Part I. A one-dimensional distributed model

The motion of a clarinet reed that is clamped to a mouthpiece and supported by a lip is simulated in the time-domain using finite difference methods. The reed is modelled as a bar with non-uniform cross section, and is described using a one-dimensional, fourth-order partial differential equation. The interactions with the mouthpiece Jay and the player's lip are taken into account by incorporating conditional contact forces in the bar equation. The model is completed by clamped-free boundary conditions for the reed. An implicit finite difference method is used for discretising the system, and values for the physical parameters are chosen both from laboratory measurements and by accurate tuning of the numerical simulations. The accuracy of the numerical system is assessed through analysis of frequency warping effects and of resonance estimation. Finally, the mechanical properties of the system are studied by analysing its response to external driving forces. In particular, the effects of reed curling are investigated.

Automated mixed dimensional modelling from 2D and 3D CAD models

The motivation for this paper is to present procedures for automatically creating idealised finite element models from the 3D CAD solid geometry of a component. The procedures produce an accurate and efficient analysis model with little effort on the part of the user. The technique is applicable to thin walled components with local complex features and automatically creates analysis models where 3D elements representing the complex regions in the component are embedded in an efficient shell mesh representing the mid-faces of the thin sheet regions. As the resulting models contain elements of more than one dimension, they are referred to as mixed dimensional models. Although these models are computationally more expensive than some of the idealisation techniques currently employed in industry, they do allow the structural behaviour of the model to be analysed more accurately, which is essential if appropriate design decisions are to be made. Also, using these procedures, analysis models can be created automatically whereas the current idealisation techniques are mostly manual, have long preparation times, and are based on engineering judgement. In the paper the idealisation approach is first applied to 2D models that are used to approximate axisymmetric components for analysis. For these models 2D elements representing the complex regions are embedded in a 1D mesh representing the midline of the cross section of the thin sheet regions. Also discussed is the coupling, which is necessary to link the elements of different dimensionality together. Analysis results from a 3D mixed dimensional model created using the techniques in this paper are compared to those from a stiffened shell model and a 3D solid model to demonstrate the improved accuracy of the new approach. At the end of the paper a quantitative analysis of the reduction in computational cost due to shell meshing thin sheet regions demonstrates that the reduction in degrees of freedom is proportional to the square of the aspect ratio of the region, and for long slender solids, the reduction can be proportional to the aspect ratio of the region if appropriate meshing algorithms are used.

Three-dimensional pattern formation, multiple homogeneous soft modes, and nonlinear dielectric electroconvection

Patterns forming spontaneously in extended, three-dimensional, dissipative systems are likely to excite several homogeneous soft modes (approximate to hydrodynamic modes) of the underlying physical system, much more than quasi-one- (1D) and two-dimensional (2D) patterns are. The reason is the lack of damping boundaries. This paper compares two analytic techniques to derive the pattern dynamics from hydrodynamics, which are usually equivalent but lead to different results when applied to multiple homogeneous soft modes. Dielectric electroconvection in nematic liquid crystals is introduced as a model for 3D pattern formation. The 3D pattern dynamics including soft modes are derived. For slabs of large but finite thickness the description is reduced further to a 2D one. It is argued that the range of validity of 2D descriptions is limited to a very small region above threshold. The transition from 2D to 3D pattern dynamics is discussed. Experimentally testable predictions for the stable range of ideal patterns and the electric Nusselt numbers are made. For most results analytic approximations in terms of material parameters are given. [S1063-651X(00)09512-X].

Detailed analytical investigation of magnetic field line random walk in turbulent plasmas: I. Two-component slab/two-dimensional turbulence

The random displacement of magnetic field lines in the presence of magnetic turbulence in plasmas is investigated from first principles. A two-component (slab/two-dimensional composite) model for the turbulence spectrum is employes. An analytical investigation of the asymptotic behavior of the field-line mean square displacement (FL-MSD) is carried out. It is shown that the magnetic field lines behave superdifusively for every large values of the position variable z, since the FL-MSD sigma varies as sigma similar to z(4/3). An intermediate diffusive regime may also possible exist for finite values of z under conditions which are explicitly determined in terms of the intrinsic turbulent plasma parameters. The superdiffusie asymptotic result is confirmed numerically via an iterative algorithm. The relevance to previous resuslts is discussed.

Boson pairs in a one-dimensional split trap

We describe the properties of a pair of ultracold bosonic atoms in a one-dimensional harmonic trapping potential with a tunable zero-ranged barrier at the trap center. The full characterization of the ground state is done by calculating the reduced single-particle density, the momentum distribution, and the two-particle entanglement. We derive several analytical expressions in the limit of infinite repulsion (Tonks-Girardeau limit) and extend the treatment to finite interparticle interactions by numerical solution. As pair interactions in double wells form a fundamental building block for many-body systems in periodic potentials, our results have implications for a wide range of problems.

Global quantum correlations in finite-size spin chains

We perform an extensive study of the properties of global quantum correlations in finite-size one-dimensional quantum spin models at finite temperature. By adopting a recently proposed measure for global quantum correlations (Rulli and Sarandy 2011 Phys. Rev. A 84 042109), called global discord, we show that critical points can be neatly detected even for many-body systems that are not in their ground state. We consider the transverse Ising model, the cluster-Ising model where three-body couplings compete with an Ising-like interaction, and the nearest-neighbor XX Hamiltonian in transverse magnetic field. These models embody our canonical examples showing the sensitivity of global quantum discord close to criticality. For the Ising model, we find a universal scaling of global discord with the critical exponents pertaining to the Ising universality class.

Effect of the crown design and interface lute parameters on the stress-state of a machined crown-tooth system: A finite element analysis

The effect of preparation design and the physical properties of the interface lute on the restored machined ceramic crown-tooth complex are poorly understood. The aim of this work was to determine, by means of three-dimensional finite element analysis (3D FEA) the effect of the tooth preparation design and the elastic modulus of the cement on the stress state of the cemented machined ceramic crown-tooth complex. The three-dimensional structure of human premolar teeth, restored with adhesively cemented machined ceramic crowns, was digitized with a micro-CT scanner. An accurate, high resolution, digital replica model of a restored tooth was created. Two preparation designs, with different occlusal morphologies, were modeled with cements of 3 different elastic moduli. Interactive medical image processing software (mimics and professional CAD modeling software) was used to create sophisticated digital models that included the supporting structures; periodontal ligament and alveolar bone. The generated models were imported into an FEA software program (hypermesh version 10.0, Altair Engineering Inc.) with all degrees of freedom constrained at the outer surface of the supporting cortical bone of the crown-tooth complex. Five different elastic moduli values were given to the adhesive cement interface 1.8 GPa, 4 GPa, 8 GPa, 18.3 GPa and 40 GPa; the four lower values are representative of currently used cementing lutes and 40 GPa is set as an extreme high value. The stress distribution under simulated applied loads was determined. The preparation design demonstrated an effect on the stress state of the restored tooth system. The cement elastic modulus affected the stress state in the cement and dentin structures but not in the crown, the pulp, the periodontal ligament or the cancellous and cortical bone. The results of this study suggest that both the choice of the preparation design and the cement elastic modulus can affect the stress state within the restored crown-tooth complex.