Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Elastodynamics #91

Open
wants to merge 109 commits into
base: devel
Choose a base branch
from
Open

Elastodynamics #91

wants to merge 109 commits into from

Conversation

lynnmunday
Copy link
Collaborator

No description provided.

mguddati and others added 30 commits December 7, 2023 19:25
    This becomes important when we have multiple parts of tensor
    mechanics, e.g. real and imaginary parts of displacement and
    their coupling.
two vectors of inner products computed to evaluate the gradient
for inverse optimization.
optimization with general optimzation reporter with separate
objective function and gradient computation, and (b) sampler
multiapp facility through stochastic tools.

This (1-D wave propagation) example is currently limited to
gradient check, with actual inversion to follow.
optimization with general optimzation reporter with separate
objective function and gradient computation, and (b) sampler
multiapp facility through stochastic tools.

This (1-D wave propagation) example is currently limited to
gradient check, with actual inversion to follow.
The convergence of course depends on the initial estimate.
with respect to convergence (narrower region and slower convergence)
…uencies)

Convergence region is much smaller, owing to the oscillatory nature of high
frequency misfit
first step. Need to check reasonableness first and then accuracy.
first step. Need to check reasonableness first and then accuracy.
first step. Need to check reasonableness first and then accuracy.
moosetest and others added 30 commits January 10, 2024 23:57
Using Hex8 elements for now, but need to use
Hex27 or others for real problems, perhaps
frequencies (synthetic data). Still requires some manual work.
Also, cleaned up various things to streamline better
a) typo in the kernels
b) clean up of the suffixes for complex arithmetic
c) chnge the velocity and frequencies to be more realistic
d) inject realistic viscoelasticity and reduce the demand
   on absorbing boundary conditions.
e) make the vectorpostprocessor addition construct in the
   gradient computation consistent with vectormath construct
   that is added recently
and has embedded dependency that lower resolution is successful before
higher resolution is attempted
and has embedded dependency that lower resolution is successful before
higher resolution is attempted

Also fixing some inconsitencies with problem change wrt domain size
and material properties
Fixing a problem with viscoelastic gradients
could be time sinks. Still need to fix the sampler batch issue.
abandon the !include in one location, but more importantly, after trial
and error, we got multiple qsub files for the specific geometry. Need to
figure out why they may be, but for now we need to do trial and error to
figure out these parameters. Otherwise the arcane memory/communication
error may turn up again, but don't know why yet...
Otherwise (normal) fails with the recent hack.
not hitting the arcane error in parallel communication, apparently
triggered by the large number of processors
3 plane 2 frequency synthetic generation using
- small mesh
- custom mesh material input - DONE (fixed the GrGaussian for rectangular mesh)
- large mesh - DONE
- check_grad - DONE
- inversion TBD
SWE synthetic inversion with rotating measurement planes and associated
pushes from acoustic simulations. The push is at the center and the
relaxation time is large (1ms), consistent with Rheospectric measureemnts
…o ISOPOD

This enables correlation objective, which may be done using general optimization reporter in the futre.
(facilitated inclusion of synthetic noise for testing as well)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants