Search results

The purpose of this paper is to present the development and application of a numerical formulation for the structural dynamics and aeroelastic analysis of new generation helicopter and tiltrotor rotor blades. These are characterized by a curvilinear elastic axis, typically with the presence of tip sweep and anhedral angles.

The structural dynamics model implemented is based on nonlinear, flap‐lag‐torsion, rotating beam equations that are valid for slender, homogeneous, isotropic, non‐uniform, twisted blades undergoing moderate displacements. A second‐order approximation scheme for strain‐displacement is adopted. Aerodynamic contributions for aeroelastic applications are derived from sectional theories, with inclusion of wake inflow models to take into account three‐dimensional effects. The numerical integration is obtained through implementation within the COMSOL Multiphysics Finite‐Element‐Method (FEM) software code, considering the elastic axis of arbitrary curvilinear shape.

The computational tool developed is validated by comparisons with results available in the literature. These demonstrate the capability of the tool to accurately predict structural dynamics and aeroelastic behavior of curved‐axis rotor blades. In particular, the influence of sweep and anhedral angles at the blade tip is successfully captured.

The numerical tool developed is limited to the analysis of isotropic blades, with a simple sectional aerodynamic modeling for aeroelastic applications. However, the flexibility of the process through which the proposed tool has been developed is such that a moderate effort is required for its extension to composite blades and more accurate aerodynamic loads predictions.

The proposed computational solver is a reliable tool for preliminary design and optimal design processes of helicopter and tiltrotor rotor blades.

Computational tools for rotors with advanced‐geometry blades are not commonly available. Therefore, the presentation of a successful way to implement structural dynamics/aeroelastic mathematical formulations for rotor blades with curvilinear elastic axis in highly flexible, multiphysics, FEM‐based, commercial software may be of interest for designers and researchers.

The purpose of this paper is to provide a framework for the implementation of an adjoint sensitivity formulation for least‐squares partial differential equations constrained optimization problems exploiting a multiphysics finite elements package. The estimation of the diffusion coefficient in a Poisson‐type diffusion equation is used as an example.

The authors derive the adjoint formulation in a continuous setting allowing to attribute to the direct and adjoint states the role of different fields to be solved for. They are one‐way coupled through the mismatch between measured and direct states acting as a source term in the adjoint equation. Having solved for the direct and adjoint state, the sensitivity of the cost function with respect to the design variables can then be obtained by a suitable post‐processing procedure. This sensitivity can then be used to efficiently solve the least‐squares problem.

The authors derived the adjoint formulation in a continuous setting allowing the direct and adjoint states to be attributed the role of different fields to be solved. They are one‐way coupled through the mismatch between measured and direct states acting as a source term in the adjoint equation. It is found that, having solved for the direct and adjoint state, the sensitivity of the cost function with respect to the design variables can then be obtained by a suitable post‐processing procedure.

This paper implies that modern multiphysics finite elements packages provide a flexible and extendable software environment for the experimentation with different adjoint formulations. Such tools are therefore expected to become increasingly important in solving notoriously difficult partial differential equation (PDE)‐constrained least‐squares problems. The framework also provides the possibility of experimentation with different regularization techniques (total variation and multiscale techniques for instance) to handle the ill‐posedness of the problem.

In this paper the adjoint sensitivity computation is casted as a multiphysics problem allowing for a flexible and extendable implementation.

The purpose of this paper is to employ the Generalized Integral Transform Technique in the analysis of conjugated heat transfer in micro-heat exchangers, by combining this hybrid numerical-analytical approach with a reformulation strategy into a single domain that envelopes all of the physical and geometric sub-regions in the original problem. The solution methodology advanced is carefully validated against experimental results from non-intrusive techniques, namely, infrared thermography measurements of the substrate external surface temperatures, and fluid temperature measurements obtained through micro Laser Induced Fluorescence.

The methodology is applied in the hybrid numerical-analytical treatment of a multi-stream micro-heat exchanger application, involving a three-dimensional configuration with triangular cross-section micro-channels. Space variable coefficients and source terms with abrupt transitions among the various sub-regions interfaces are then defined and incorporated into this single domain representation for the governing convection-diffusion equations. The application here considered for analysis is a multi-stream micro-heat exchanger designed for waste heat recovery and built on a PMMA substrate to allow for flow visualization.

The methodology here advanced is carefully validated against experimental results from non-intrusive techniques, namely, infrared thermography measurements of the substrate external surface temperatures and fluid temperature measurements obtained through Laser Induced Fluorescence. A very good agreement among the proposed hybrid methodology predictions, a finite elements solution from the COMSOL code, and the experimental findings has been achieved. The proposed methodology has been demonstrated to be quite flexible, robust, and accurate.

The hybrid nature of the approach, providing analytical expressions in all but one independent variable, and requiring numerical treatment at most in one single independent variable, makes it particularly well suited for computationally intensive tasks such as in optimization, inverse problem analysis, and simulation under uncertainty.

The purpose of this paper is to propose and analyze a combined heat treatment of metal materials, consisting in classic induction pre-heating and/or post-heating and full heating by laser beam. This technology is prospective for some kinds of surface hardening and welding because its application leads to lowering of temperature gradients at the heated spots, which substantially reduces local residual mechanical strains and stresses.

The task was solved like the 3D hard-coupled problem for electromagnetic field, temperature field and field of displacements. It was solved numerically using the techniques based on the FEM. For solution was used commercial software COMSOL Multiphysics, some parts were solved using own scripts in the software Agros.

In the paper are shown results of the numerical solution and experimental measured data. Due the work the authors found that the influence of the pre-heating and post-heating really leads to limit the temperature gradients and from other measurements is clear that also to improving of the welding.

The paper presents fully 3D nonlinear and nonstationary mathematical model of hybrid laser welding, its numerical solution experimental verification.

The purpose of this paper was to promote the use of residual moringa seed powder (RMSP) for the enhancement of cereal-based products. RMSP is usually discarded after seed-oil extraction. This work also promotes zero-waste and rheological approaches.

In search of novel and sustainable food products with high nutritional value, cold-pressed Moringa oleifera Lam. seeds residue (RMSP) was used for incorporation in muffin formulations. Wheat flour was partially substituted (0%, 1%, 3%, 5%, 7% and 9%) by RMSP. Sodium (Na), calcium (Ca) and iron (Fe) contents were quantified through atomic absorption spectrometry; protein, through the Kjeldahl method followed by AACC Method 46–13.01; and, fat content, by a modified version of AACC Method 30–25.01a. Analysis of variance (ANOVA) and Tukey tests were performed to determine significant differences between formulations at 95% reliability using Minitab® software. Furthermore, simple viscosity studies of the dough mixture were carried in a Brookfield DV-III Ultra Rheometer; Matlab® curve-fitting tool was used for obtaining the best non-Newtonian equation that modeled experimental data. Subsequently, computational fluid dynamics (CFD) simulations of non-Newtonian fluids along a segmented pipe were carried out in Comsol Multiphysics® software to depict the importance of modeling non-Newtoning fluids for downstream processes.

RMSP significantly (p < 0.0001) increased protein and Fe content for the 7% and 9% formulations; it dramatically changed Na and Ca content in all formulations, while fat remained constant.

The development of this type of product is an opportunity for communities that grow and harvest moringa as well as for food industries which can take advantage of moringa by-products for several subsequent processing.

For the first time, it was found that dough formulations with RMSP presented a pseudo-plastic and thixotropic behavior. In addition, the use of lignocellulosic by-products such as RMSP incorporates an added value to food products. In this case, it was demonstrated that moringa seed residue enhanced nutritional value to muffins and provided coagulant/flocculant action, which is essential during dough preparation.

The purpose of this study is to use a weak light source with spatial distribution to realize light-driven fluid by adding high-absorbing nanoparticles to the droplets, thereby replacing a highly focused strong linear light source acting on pure droplets.

First, Fe₃O₄ nanoparticles with high light response characteristics were added to the droplets to prepare nanofluid droplets, and through the Gaussian light-driven flow experiment, the Marangoni effect inside a nanofluid droplet was studied, which can produce the surface tension gradient on the air/liquid interface and induce the vortex motion inside a droplet. Then, the numerical simulation method of multiphysics field coupling was used to study the effects of droplet height and Gaussian light distribution on the flow characteristics inside a droplet.

Nanoparticles can significantly enhance the light absorption, so that the Gaussian light is enough to drive the flow, and the formation of vortex can be regulated by light distribution. The multiphysics field coupling model can accurately describe this problem.

This study is helpful to understand the flow behavior and heat transfer phenomenon in optical microfluidic systems, and provides a feasible way to construct the rapid flow inside a tiny droplet by light.

Among the various applications involving the use of microwave energy, its growing utility within the mining industry is particularly noteworthy. Conventional grinding processes are often overburdened by energy inefficiencies that are directly related to machine wear, pollution and rising project costs. In this work, we numerically investigate the effects of microwave pretreatment through a series of compression tests as a means to help mitigate these energy inefficiencies.

We investigate the effects of microwave pretreatment on various rock samples, as quantified by uniaxial compression tests. In particular, we assign sample heterogeneity based on a Gaussian statistical distribution and invoke a damage model for elemental tensile and compressive stresses based on the maximum tensile stress and the Mohr–Coulomb theories, respectively. We further couple the electromagnetic, thermal and solid displacement relations using finite element modeling.

(1) Increased power intensity during microwave pretreatment results in decreased axial compressive stress. (2) Leveraging statistics to induce variable compressive and tensile strength can greatly facilitate sample heterogeneity and prove necessary for damage modeling. (3) There exists a nonlinear trend to the reduction in smax with increasing power levels, implying an optimum energy output efficiency to create the maximum degradation-power cost relationship.

Previous research in this area has been largely limited to two-dimensional thermo-electric models. The onset of high-performance computing has allowed for the development of high-fidelity, three-dimensional models with coupled equations for electromagnetics, heat transfer and solid mechanics.

The purpose of this paper is to extend the penalty concept to treat partial slip, free surface, contact and related boundary conditions in viscous flow simulation.

The penalty partial‐slip formulation is analysed and related to the classical Navier slip condition. The same penalty scheme also allows partial penetration through a boundary, hence the implementation of porous wall boundaries. The finite element method is used for investigating and interpreting penalty approaches to boundary conditions.

The generalised penalty approach is verified by means of a novel variant of the circular‐Couette flow problem, having partial slip on one of the cylindrical boundaries, for which an analytic solution is derived. Further verificationis provided by consideration of viscous flow over a sphere with partial slip on the surface, and comparison of numerical and classical solutions. Numerical studies illustrate the versatility of the approach.

The penalty approach is applied to some different boundaries: partial slip and partial penetration with no/full slip/penetration as limiting cases; free surface; space‐ and time‐varying boundary conditions which allow progressive contact over time. Application is made to curved and inclined boundaries. Sensitivity of flow to penalty parameters is an avenue for continued research, as is application of the penalty approach for non‐Newtonian flows.

This is the first work to show the relation between penalty formulation of boundary conditions and physical boundary conditions. It provides a method that overcomes past difficulties in implementing partial slip on boundaries of general shape, and which handles progressive contact. It also provides useful benchmark problems for future studies.

For the past decade, plasma actuators have been identified as a subset in the realm of active flow control devices. As research into plasma actuators continues to mature, computational modelling is needed to complement the investigation of the actuators. This paper seeks to address these issues.

In this study, the Suzen‐Huang model is chosen because of its ability to simulate both the charge density and Lorentz body force. Its advantages and limitations have been identified with a parametric study of two constants used in the modelling: the Debye length (λD) and the maximum charge density value (ρc* ). By varying the two scalars, the effects of charge density, body force and induced velocity are examined.

The results show that the non‐dimensionalised body force (Fb*) is nonlinearly dependent on Debye length. However, a linear variation of Fb* is observed with increasing values of maximum charge density. The optimized form of the Suzen‐Huang model shows better agreement in the horizontal velocity profile but still points to inaccuracy when compared to vertical velocity profile.

The results indicate that the body force still has to be modelled more extensively above the encapsulated electrode, so that the horizontal and vertical components of induced velocities are accurately obtained.

The purpose of the paper is to clarify the influence of introducing magnetic concentrators on the performance of Hall effect based current sensors and to obtain dependencies of the sensor characteristics on the conductor position.

The finite element method and Comsol software are used for analysis of the three-dimensional magnetic field of the constructions of Hall effect based current sensor with different types of magnetic concentrators – closed-core (of rectangular and toroidal type) and open-core of toroidal type – with additional larger air gap. The Hall plate is also included in the model with its real dimensions and the magnetic flux density is obtained by integrating over its volume.

It has been found that there is dependence of the output signal (proportional to the magnetic flux density) of Hall effect based current sensor with both closed- and open-core magnetic concentrators on the position of the current carrying conductor. Distribution of the magnetic flux density and dependencies of its value in the Hall plate on the conductor position and on the additional air gap have been obtained. Optimization is carried out with respect to the additional air gap and cross-section dimensions of the concentrator.

Estimation of the influence of the introducing magnetic concentrators is made with respect to relationships between the output signal and conductor position for different constructions of the magnetic core of the concentrators.