Abstract: Modeling dam-break flows over non-flat beds requires
an accurate representation of the topography which is the main
source of uncertainty in the model. Therefore, developing robust
and accurate techniques for reconstructing topography in this class
of problems would reduce the uncertainty in the flow system. In
many hydraulic applications, experimental techniques have been
widely used to measure the bed topography. In practice, experimental
work in hydraulics may be very demanding in both time and cost.
Meanwhile, computational hydraulics have served as an alternative
for laboratory and field experiments. Unlike the forward problem,
the inverse problem is used to identify the bed parameters from the
given experimental data. In this case, the shallow water equations
used for modeling the hydraulics need to be rearranged in a way
that the model parameters can be evaluated from measured data.
However, this approach is not always possible and it suffers from
stability restrictions. In the present work, we propose an adaptive
optimal control technique to numerically identify the underlying bed
topography from a given set of free-surface observation data. In this
approach, a minimization function is defined to iteratively determine
the model parameters. The proposed technique can be interpreted
as a fractional-stage scheme. In the first stage, the forward problem
is solved to determine the measurable parameters from known data.
In the second stage, the adaptive control Ensemble Kalman Filter is
implemented to combine the optimality of observation data in order to
obtain the accurate estimation of the topography. The main features
of this method are on one hand, the ability to solve for different
complex geometries with no need for any rearrangements in the
original model to rewrite it in an explicit form. On the other hand, its
achievement of strong stability for simulations of flows in different
regimes containing shocks or discontinuities over any geometry.
Numerical results are presented for a dam-break flow problem over
non-flat bed using different solvers for the shallow water equations.
The robustness of the proposed method is investigated using different
numbers of loops, sensitivity parameters, initial samples and location
of observations. The obtained results demonstrate high reliability and
accuracy of the proposed techniques.
Abstract: A two-step multigrid approach is proposed to solve the inverse heat conduction problem in a 3-D object under laser irradiation. In the first step, the location of the laser center is estimated using a coarse and uniform grid system. In the second step, the front-surface temperature is recovered in good accuracy using a multiple grid system in which fine mesh is used at laser spot center to capture the drastic temperature rise in this region but coarse mesh is employed in the peripheral region to reduce the total number of sensors required. The effectiveness of the two-step approach and the multiple grid system are demonstrated by the illustrative inverse solutions. If the measurement data for the temperature and heat flux on the back surface do not contain random error, the proposed multigrid approach can yield more accurate inverse solutions. When the back-surface measurement data contain random noise, accurate inverse solutions cannot be obtained if both temperature and heat flux are measured on the back surface.
Abstract: The noise requirements for naval and research vessels
have seen an increasing demand for quieter ships in order to fulfil
current regulations and to reduce the effects on marine life. Hence,
new methods dedicated to the characterization of propeller noise,
which is the main source of noise in the far-field, are needed. The
study of cavitating propellers in closed-section is interesting for
analyzing hydrodynamic performance but could involve significant
difficulties for hydroacoustic study, especially due to reverberation
and boundary layer noise in the tunnel. The aim of this paper
is to present a numerical methodology for the identification of
hydroacoustic sources on marine propellers using hydrophone arrays
in a large hydrodynamic tunnel. The main difficulties are linked to the
reverberation of the tunnel and the boundary layer noise that strongly
reduce the signal-to-noise ratio. In this paper it is proposed to estimate
the reflection coefficients using an inverse method and some reference
transfer functions measured in the tunnel. This approach allows to
reduce the uncertainties of the propagation model used in the inverse
problem. In order to reduce the boundary layer noise, a cleaning
algorithm taking advantage of the low rank and sparse structure of the
cross-spectrum matrices of the acoustic and the boundary layer noise
is presented. This approach allows to recover the acoustic signal even
well under the boundary layer noise. The improvement brought by
this method is visible on acoustic maps resulting from beamforming
and DAMAS algorithms.
Abstract: This article presents a method of using the one
dimensional piezo-electric patch on beam model for structural
identification. A hybrid element constituted of one dimensional
beam element and a PZT sensor is used with reduced material
properties. This model is convenient and simple for identification
of beams. Accuracy of this element is first verified against a
corresponding 3D finite element model (FEM). The structural
identification is carried out as an inverse problem whereby
parameters are identified by minimizing the deviation between
the predicted and measured voltage response of the patch, when
subjected to excitation. A non-classical optimization algorithm
Particle Swarm Optimization is used to minimize this objective
function. The signals are polluted with 5% Gaussian noise to
simulate experimental noise. The proposed method is applied on
beam structure and identified parameters are stiffness and damping.
The model is also validated experimentally.
Abstract: This paper presents an effective model updating strategy for damage localization and quantification in frames by defining damage detection problem as an optimization issue. A generalized version of the Modal Residual Force (MRF) is employed for presenting a new damage-sensitive cost function. Then, Grey Wolf Optimization (GWO) algorithm is utilized for solving suggested inverse problem and the global extremums are reported as damage detection results. The applicability of the presented method is investigated by studying different damage patterns on the benchmark problem of the IASC-ASCE, as well as a planar shear frame structure. The obtained results emphasize good performance of the method not only in free-noise cases, but also when the input data are contaminated with different levels of noises.
Abstract: In this paper, we investigate certain spaces of
generalized functions for the Fourier and Fourier type integral
transforms. We discuss convolution theorems and establish certain
spaces of distributions for the considered integrals. The new Fourier
type integral is well-defined, linear, one-to-one and continuous with
respect to certain types of convergences. Many properties and an
inverse problem are also discussed in some details.
Abstract: This article presents a numerical method to find the
heat flux in an inhomogeneous inverse heat conduction problem with
linear boundary conditions and an extra specification at the terminal.
The method is based upon applying the satisfier function along with
the Ritz-Galerkin technique to reduce the approximate solution of the
inverse problem to the solution of a system of algebraic equations.
The instability of the problem is resolved by taking advantage of
the Landweber’s iterations as an admissible regularization strategy.
In computations, we find the stable and low-cost results which
demonstrate the efficiency of the technique.
Abstract: In this paper the influence of errors of function derivatives in initial time which have been obtained by experiment (uncontrollable inaccuracy) to the results of inverse problem solution was investigated. It was shown that these errors distort the inverse problem solution as a rule near the beginning of interval where the solutions are analyzed. Several methods for removing the influence of uncontrollable inaccuracy have been suggested.
Abstract: In this paper, the specific sound Transmission Loss
(TL) of the Laminated Composite Plate (LCP) with different material
properties in each layer is investigated. The numerical method to
obtain the TL of the LCP is proposed by using elastic plate theory. The
transfer matrix approach is novelty presented for computational
efficiency in solving the numerous layers of dynamic stiffness matrix
(D-matrix) of the LCP. Besides the numerical simulations for
calculating the TL of the LCP, the material properties inverse method
is presented for the design of a laminated composite plate analogous to
a metallic plate with a specified TL. As a result, it demonstrates that
the proposed computational algorithm exhibits high efficiency with a
small number of iterations for achieving the goal. This method can be
effectively employed to design and develop tailor-made materials for
various applications.
Abstract: The formulated problem of optimization of the
technological process of water treatment for thermal power plants is
considered in this article. The problem is of multiparametric nature.
To optimize the process, namely, reduce the amount of waste water, a
new technology was developed to reuse such water. A mathematical
model of the technology of wastewater reuse was developed.
Optimization parameters were determined. The model consists of a
material balance equation, an equation describing the kinetics of ion
exchange for the non-equilibrium case and an equation for the ion
exchange isotherm. The material balance equation includes a
nonlinear term that depends on the kinetics of ion exchange. A direct
problem of calculating the impurity concentration at the outlet of the
water treatment plant was numerically solved. The direct problem
was approximated by an implicit point-to-point computation
difference scheme. The inverse problem was formulated as relates to
determination of the parameters of the mathematical model of the
water treatment plant operating in non-equilibrium conditions. The
formulated inverse problem was solved. Following the results of
calculation the time of start of the filter regeneration process was
determined, as well as the period of regeneration process and the
amount of regeneration and wash water. Multi-parameter
optimization of water treatment process for thermal power plants
allowed decreasing the amount of wastewater by 15%.
Abstract: In this study, a comparative analysis of the approaches
associated with the use of neural network algorithms for effective
solution of a complex inverse problem – the problem of identifying
and determining the individual concentrations of inorganic salts in
multicomponent aqueous solutions by the spectra of Raman
scattering of light – is performed. It is shown that application of
artificial neural networks provides the average accuracy of
determination of concentration of each salt no worse than 0.025 M.
The results of comparative analysis of input data compression
methods are presented. It is demonstrated that use of uniform
aggregation of input features allows decreasing the error of
determination of individual concentrations of components by 16-18%
on the average.
Abstract: The sound pressure level (SPL) of the moving-coil
loudspeaker (MCL) is often simulated and analyzed using the lumped
parameter model. However, the SPL of a MCL cannot be simulated
precisely in the high frequency region, because the value of cone
effective area is changed due to the geometry variation in different
mode shapes, it is also related to affect the acoustic radiation mass and
resistance. Herein, the paper presents the inverse method which has a
high ability to measure the value of cone effective area in various
frequency points, also can estimate the MCL electroacoustic
parameters simultaneously. The proposed inverse method comprises
the direct problem, adjoint problem, and sensitivity problem in
collaboration with nonlinear conjugate gradient method. Estimated
values from the inverse method are validated experimentally which
compared with the measured SPL curve result. Results presented in
this paper not only improve the accuracy of lumped parameter model
but also provide the valuable information on loudspeaker cone design.
Abstract: The aim of the current work is to present a comparison among three popular optimization methods in the inverse elastostatics problem (IESP) of flaw detection within a solid. In more details, the performance of a simulated annealing, a Hooke & Jeeves and a sequential quadratic programming algorithm was studied in the test case of one circular flaw in a plate solved by both the boundary element (BEM) and the finite element method (FEM). The proposed optimization methods use a cost function that utilizes the displacements of the static response. The methods were ranked according to the required number of iterations to converge and to their ability to locate the global optimum. Hence, a clear impression regarding the performance of the aforementioned algorithms in flaw identification problems was obtained. Furthermore, the coupling of BEM or FEM with these optimization methods was investigated in order to track differences in their performance.
Abstract: A new technique of topological multi-scale analysis is
introduced. By performing a clustering recursively to build a
hierarchy, and analyzing the co-scale and intra-scale similarities, an
Iterated Function System can be extracted from any data set. The study
of fractals shows that this method is efficient to extract
self-similarities, and can find elegant solutions the inverse problem of
building fractals. The theoretical aspects and practical
implementations are discussed, together with examples of analyses of
simple fractals.
Abstract: In this paper a numerical algorithm is described for solving the boundary value problem associated with axisymmetric, inviscid, incompressible, rotational (and irrotational) flow in order to obtain duct wall shapes from prescribed wall velocity distributions. The governing equations are formulated in terms of the stream function ψ (x,y)and the function φ (x,y)as independent variables where for irrotational flow φ (x,y)can be recognized as the velocity potential function, for rotational flow φ (x,y)ceases being the velocity potential function but does remain orthogonal to the stream lines. A numerical method based on the finite difference scheme on a uniform mesh is employed. The technique described is capable of tackling the so-called inverse problem where the velocity wall distributions are prescribed from which the duct wall shape is calculated, as well as the direct problem where the velocity distribution on the duct walls are calculated from prescribed duct geometries. The two different cases as outlined in this paper are in fact boundary value problems with Neumann and Dirichlet boundary conditions respectively. Even though both approaches are discussed, only numerical results for the case of the Dirichlet boundary conditions are given. A downstream condition is prescribed such that cylindrical flow, that is flow which is independent of the axial coordinate, exists.
Abstract: Displacement measurement was conducted on compact normal and shear specimens made of acrylic homogeneous material subjected to mixed-mode loading by digital image correlation. The intelligent hybrid method proposed by Nishioka et al. was applied to the stress-strain analysis near the crack tip. The accuracy of stress-intensity factor at the free surface was discussed from the viewpoint of both the experiment and 3-D finite element analysis. The surface images before and after deformation were taken by a CMOS camera, and we developed the system which enabled the real time stress analysis based on digital image correlation and inverse problem analysis. The great portion of processing time of this system was spent on displacement analysis. Then, we tried improvement in speed of this portion. In the case of cracked body, it is also possible to evaluate fracture mechanics parameters such as the J integral, the strain energy release rate, and the stress-intensity factor of mixed-mode. The 9-points elliptic paraboloid approximation could not analyze the displacement of submicron order with high accuracy. The analysis accuracy of displacement was improved considerably by introducing the Newton-Raphson method in consideration of deformation of a subset. The stress-intensity factor was evaluated with high accuracy of less than 1% of the error.
Abstract: An inverse problem of doubly center matrices is discussed. By translating the constrained problem into unconstrained problem, two iterative methods are proposed. A numerical example illustrate our algorithms.
Abstract: The overall objective of this paper is to retrieve soil
surfaces parameters namely, roughness and soil moisture related to
the dielectric constant by inverting the radar backscattered signal
from natural soil surfaces.
Because the classical description of roughness using statistical
parameters like the correlation length doesn't lead to satisfactory
results to predict radar backscattering, we used a multi-scale
roughness description using the wavelet transform and the Mallat
algorithm. In this description, the surface is considered as a
superposition of a finite number of one-dimensional Gaussian
processes each having a spatial scale. A second step in this study
consisted in adapting a direct model simulating radar backscattering
namely the small perturbation model to this multi-scale surface
description. We investigated the impact of this description on radar
backscattering through a sensitivity analysis of backscattering
coefficient to the multi-scale roughness parameters.
To perform the inversion of the small perturbation multi-scale
scattering model (MLS SPM) we used a multi-layer neural network
architecture trained by backpropagation learning rule. The inversion
leads to satisfactory results with a relative uncertainty of 8%.
Abstract: In this paper, we establish existence and uniqueness of
solutions for a class of inverse problems of degenerate differential
equations. The main tool is the perturbation theory for linear operators.
Abstract: A thin layer on the component surface can be found
with high tensile residual stresses, due to turning operations, which
can dangerously affect the fatigue performance of the component. In
this paper an analytical approach is presented to reconstruct the
residual stress field from a limited incomplete set of measurements.
Airy stress function is used as the primary unknown to directly solve
the equilibrium equations and satisfying the boundary conditions. In
this new method there exists the flexibility to impose the physical
conditions that govern the behavior of residual stress to achieve a
meaningful complete stress field. The analysis is also coupled to a
least squares approximation and a regularization method to provide
stability of the inverse problem. The power of this new method is
then demonstrated by analyzing some experimental measurements
and achieving a good agreement between the model prediction and
the results obtained from residual stress measurement.